skip to main content

Explore Research Products in the NSF-PAR

Title: Progress toward a one-electron model for the non-valence correlation-bound anions of polycyclic aromatic hydrocarbons
Abstract A one-electron model Hamiltonian is used to characterize the non-valence correlation-bound (NVCB) anions of hexagonal polycyclic aromatic hydrocarbons (PAHs) C 6 n 2 H 6 n ( n = 3–7). The model incorporates atomic electrostatic moments up to the quadrupole, coupled inducible charges and dipoles, and atom-centered repulsive Gaussians to describe the interaction between the excess electron and PAH. These model components are parameterized on and validated against all-electron calculations. Good agreement is found between the static dipole polarizabilities obtained from the model and those from PBE0 density functional theory and second-order Møller–Plesset perturbation theory calculations. In the model, charge flow dominates the in-plane polarizability of PAHs larger than C 54 H 18 , yielding an approximately quadratic scaling of the mean polarizabilty with the number of carbon atoms. Inclusion of electrostatic interactions decreases the electron binding energies for the largest PAHs considered by about 20% and shift charge distribution from above and below the plane of the ring system toward the periphery. Analysis of the electrostatic and polarization interactions provides insight into qualitative trends in the electron binding energy and the charge distribution of the lowest energy NVCB anion.  more » « less
Award ID(s):
2028826
NSF-PAR ID:
10353228
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Electronic Structure
Volume:
4
Issue:
1
ISSN:
2516-1075
Page Range / eLocation ID:
014010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ( , Structural chemistry)
    null (Ed.)
    Density functional theory (DFT) calculations were used to examine the binding strength of one and two methane molecule(s) with graphene (62 and 186 carbon atoms) and model systems of aromatic hydrocarbons (benzene, pyrene, and coronene). We explored different possibilities of binding modes of methane such as one, two, and three C-H interacting with small π-systems. Two methane molecules were considered to bind from the same as well as opposite sides of the plane of benzene and other πsystems including graphene models. Our results show that methane molecule prefers to bind with three C-H…π interactions with all the π-systems except benzene. The preference of tripod configuration of methane on the surface of graphene systems strongly agrees with the neutron diffraction experiment of methane on graphitized carbon black. The binding strength is almost doubled by increasing the number of methane molecules from one to two. Importantly, two methane molecules prefer to bind on the same side rather than opposite sides of the plane of graphene due to stabilizing CH…HC interactions between them in addition to six CH…π interactions. Interestingly, binding strength contributions from CH…HC interactions (approx. 0.4–0.5 kcal/mol) of two methane molecules on the surface are analogous to methane dimer complex free from the surface of graphene. C-H stretching frequency shifts, bond lengths, and binding distances support the presence of CH…HC interactions between two methane molecules. Structures of complexes, binding energies, and C-H stretching frequency shifts agree with available experimental data 
    more » « less
  2. ; ; ; ( , Journal of Computational Chemistry)

    The model reactions CH3X + (NH—CH=O)M ➔ CH3—NH—NH═O or NH═CH—O—CH3 + MX (M = none, Li, Na, K, Ag, Cu; X = F, Cl, Br) are investigated to demonstrate the feasibility of Marcus theory and the hard and soft acids and bases (HSAB) principle in predicting the reactivity of ambident nucleophiles. The delocalization indices (DI) are defined in the framework of the quantum theory of atoms in molecules (QT‐AIM), and are used as the scale of softness in the HSAB principle. To react with the ambident nucleophile NH═CH—O−, the carbocation H3C+from CH3X (F, Cl, Br) is actually a borderline acid according to the DI values of the forming C…N and C…O bonds in the transition states (between 0.25 and 0.49), while the counter ions are divided into three groups according to the DI values of weak interactions involving M (M…X, M…N, and M…O): group I (M = none, and Me4N) basically show zero DI values; group II species (M = Li, Na, and K) have noticeable DI values but the magnitudes are usually less than 0.15; and group III species (M = Ag and Cu(I)) have significant DI values (0.30–0.61). On a relative basis, H3C+is a soft acid with respect to group I and group II counter ions, and a hard acid with respect to group III counter ions. Therefore, N‐regioselectivity is found in the presence of group I and group II counter ions (M = Me4N, Li, Na, K), while O‐regioselectivity is observed in the presence of the group III counter ions (M = Ag, and Cu(I)). The hardness of atoms, groups, and molecules is also calculated with new functions that depend on ionization potential (I) and electron affinity (A) and use the atomic charges obtained from localization indices (LI), so that the regioselectivity is explained by the atomic hardness of reactive nitrogen atoms in the transition states according to the maximum hardness principle (MHP). The exact Marcus equation is derived from the simple harmonic potential energy parabola, so that the concepts of activation free energy, intrinsic activation barrier, and reaction energy are completely connected. The required intrinsic activation barriers can be either estimated fromab initiocalculations on reactant, transition state, and product of the model reactions, or calculated from identity reactions. The counter ions stabilize the reactant through bridging N‐ and O‐site of reactant of identity reactions, so that the intrinsic barriers for the salts are higher than those for free ambident anions, which is explained by the increased reorganization parameter Δr. The proper application of Marcus theory should quantitatively consider all three terms of Marcus equation, and reliably represent the results with potential energy parabolas for reactants and all products. For the model reactions, both Marcus theory and HSAB principle/MHP principle predict the N‐regioselectivity when M = none, Me4N, Li, Na, K, and the O‐regioselectivity when M = Ag and Cu(I). © 2019 Wiley Periodicals, Inc.

     
    more » « less
  3. ; ; ( , Journal of Computational Chemistry)
    Abstract

    We report the application of our fragment‐based quantum chemistry model MIM (Molecules‐In‐Molecules) with electrostatic embedding. The method is termed “EE‐MIM (ElectrostaticallyEmbeddedMolecules‐In‐Molecules)” and accounts for the missing electrostatic interactions in the subsystems resulting from fragmentation. Point charges placed at the atomic positions are used to represent the interaction of each subsystem with the rest of the molecule with minimal increase in the computational cost. We have carefully calibrated this model on a range of different sizes of clusters containing up to 57 water molecules. The fragmentation methods have been applied with the goal of reproducing the unfragmented total energy at the MP2/6‐311G(d,p) level. Comparative analysis has been carried out between MIM and EE‐MIM to gauge the impact of electrostatic embedding. Performance of several different parameters such as the type of charge and levels of fragmentation are analyzed for the prediction of absolute energies. The use of background charges in subsystem calculations improves the performance of both one‐ and two‐layer MIM while it is noticeably important in the case of one‐layer MIM. Embedded charges for two‐layer MIM are obtained from a full system calculation at the low‐level. For one‐layer MIM, in the absence of a full system calculation, two different types of embedded charges, namely, Geometry dependent (GD) and geometry independent (GI) charges, are used. A self‐consistent procedure is employed to obtain GD charges. We have further tested our method on challenging charged systems with stronger intermolecular interactions, namely, protonated ammonia clusters (containing up to 30 ammonia molecules). The observations are similar to water clusters with improved performance using embedded charges. Overall, the performance of NPA charges as embedded charges is found to be the best.

     
    more » « less
  4. ; ; ; ; ( , Physical Chemistry Chemical Physics)
    Previous microwave studies of naphthol monomers were supplemented by measuring spectra of all 13 C mono-substituted isotopologues of the cis - and trans -conformers of 1-naphthol and 2-naphthol in their natural abundances. The resulting data were utilized to determine substitution structures and so-called semi-experimental effective structures. Results from electronic structure calculations show that the OH group of cis -1-naphthol points ≈6° out of plane, which is consistent with the inertial defect data of cis - and trans -1-naphthol. The non-planarity of cis -1-naphthol is a result of a close-contact H-atom–H-atom interaction. This type of H–H interaction has been the subject of much controversy in the past and we provide here an in-depth theoretical analysis of it. The naphthol system is particularly well-suited for such analysis as it provides internal standards with its four different isomers. The methods used include quantum theory of atoms in molecules, non-covalent interactions, independent gradient model, local vibrational mode, charge model 5, and natural bond orbital analyses. We demonstrate that the close-contact H–H interaction is neither a purely attractive nor repulsive interaction, but rather a mixture of the two. 
    more » « less
  5. ; ; ; ; ; ( , Dalton Transactions)
    null (Ed.)
    A series of complexes with low-energy Fe II to Ti IV metal-to-metal charge-transfer (MMCT) transitions, Cp 2 Ti(C 2 Fc) 2 , Cp* 2 Ti(C 2 Fc) 2 , and MeOOC Cp 2 Ti(C 2 Fc) 2 , was investigated using solvatochromism and resonance Raman spectroscopy (RRS) augmented with time-dependent density functional theory (TDDFT) calculations in order to interrogate the nature of the CT transitions. Computational models were benchmarked against the experimental UV-Vis spectra and B3LYP/6-31G(d) was found to most faithfully represent the spectra. The energy of the MMCT transition was measured in 15 different solvents and a multivariate fit to the Catalán solvent parameters – solvent polarizability (SP), solvent dipolarity (SdP), solvent basicity (SB), and solvent acidity (SA) – was performed. The effect of SP indicates a greater degree of electron delocalization in the excited state (ES) than the ground state (GS). The small negative solvatochromism with respect to SdP indicates a smaller dipole moment in the ES than the GS. The effect of SB is consistent with charge-transfer to Ti. Upon excitation into the MMCT absorption band, the RRS data show enhancement of the alkyne stretching modes and of the out-of-plane bending modes of the cyclopentadienyl ring connected to Fe and the alkyne bridge. This is consistent with changes in the oxidation states of Ti and Fe, respectively. The higher-energy transitions (350–450 nm) show enhancement of vibrational modes consistent with ethnylcyclopentadienyl to Ti ligand-to-metal charge transfer (LMCT). The RRS data is consistent with the TDDFT predicted character of these transitions. TDDFT suggests that the lowest-energy transition in Cp 2 Ti(C 2 Fc) 2 CuI, where CuI is coordinated between the alkynes, retains its Fe II to Ti IV MMCT character, in agreement with the RRS data, but that the lowest-energy transitions have significant CuI to Ti character. For Cp 2 Ti(C 2 Fc) 2 CuI, excitation into the low-energy MMCT absorption band results in selective enhancement of the symmetric alkynyl stretching mode. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email