Assessment of Density Functional Theory in Predicting Interaction Energies Between Water and Polycyclic Aromatic Hydrocarbons: From Water on Benzene to Water on Graphene

The interactions of water with polycyclic aromatic hydrocarbons, from benzene to graphene, are investigated using various exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. The accuracy of the different functionals is assessed through comparisons with random phase approximation (RPA) and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] calculations. Diffusion Monte Carlo (DMC) data reported in the literature are also used for comparison. Relatively large variations are found in interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules. The meta-GGA B97M-rV and range-separated hybrid, meta-GGA ωB97M-V functionals provide nearly quantitative agreement with CCSD(T) values for the water–benzene, water–coronene, and water–circumcoronene dimers, while RPA and DMC predict interaction energies that differ by up to ∼1 kcal/mol and ∼0.4 kcal/mol from the corresponding CCSD(T) values, respectively. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for larger polycyclic aromatic hydrocarbons. By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA), it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet

Assessment of Density Functional Theory in Predicting Interaction Energies Between Water and Polycyclic Aromatic Hydrocarbons: From Water on Benzene to Water on Graphene

The interactions of water with polycyclic aromatic hydrocarbons, from benzene to graphene, are investigated using various exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. The accuracy of the different functionals is assessed through comparisons with random phase approximation (RPA) and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] calculations. Diffusion Monte Carlo (DMC) data reported in the literature are also used for comparison. Relatively large variations are found in interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules. The meta-GGA B97M-rV and range-separated hybrid, meta-GGA ωB97M-V functionals provide nearly quantitative agreement with CCSD(T) values for the water–benzene, water–coronene, and water–circumcoronene dimers, while RPA and DMC predict interaction energies that differ by up to ∼1 kcal/mol and ∼0.4 kcal/mol from the corresponding CCSD(T) values, respectively. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for larger polycyclic aromatic hydrocarbons. By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA), it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet

Synthesis and reductive chemistry of bimetallic and trimetallic rare-earth metallocene hydrides with (C5H4SiMe3)1− ligands

The reductive chemistry of [Cp\u272Ln(μ–H)(THF)x]y [Ln = Y, Dy, Tb; Cp\u27 = (C5H4SiMe3)1−; x = 2, 0 and y = 2, 3] was examined to determine if these hydrides would be viable precursors for 4fn5d1 Ln2+ ions that could form 5d1-5d1 metal–metal bonded complexes. The hydrides were prepared by reaction of the chlorides, [Cp\u272Ln(μ–Cl)]2, 1-Ln, with allylmagnesium chloride to form the allyl complexes, [Cp\u272Y(η3–C3H5)(THF)], 2-Ln, which were hydrogenolyzed. The solvent-free reaction of solid 2-Ln with 60 psi of H2 gas in a Fischer-Porter apparatus produced, in the Y case, the trimetallic species, [Cp\u272Y(μ–H)]3, 3-Y, and in the Dy and Tb cases, the bimetallic complexes [Cp\u272Ln(μ–H)(THF)]2, 4-Ln (Ln = Dy, Tb). The latter complexes could be converted to 3-Dy and 3-Tb by heating under vacuum. Isopiestic data indicate that 3-Y solvates to 4-Y in THF. Reductions of 4-Y, 4-Dy, and 4-Tb with KC8 in the presence of a chelate such as 2.2.2-cryptand or 18-crown-6 all gave reaction products with intense dark colors characteristic of Ln2+ ions. In the yttrium case, with either chelating agent, the dark green product gives a rhombic EPR spectrum (g1 = 2.01, g2 = 1.99, g3 = 1.98, A = 24.1 G) at 77 K. However, the only crystallographically-characterizable products obtainable from these solutions were Ln3+polyhydride anion complexes of composition, [K(chelate)]{[Cp\u272Ln(μ–H)]3(μ–H)}

EMSL Science Theme Advisory Panel Workshop - Atmospheric Aerosol Chemistry, Climate Change, and Air Quality

This report contains the workshop scope and recommendations from the workshop attendees in identifying scientific gaps in new particle formation, growth and properties of particles and reactions in and on particles as well as the laboratory-focused capabilities, field-deployable capabilities and modeling/theory tools along with linking of models to fundamental data

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light–matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE’s functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree–Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties

A Parameter-Free Density Functional that Works for Noncovalent Interactions

The failure of semilocal density functional theory for medium- and long-range noncovalent molecular interactions is a long-standing challenge for computational chemistry. Here, we assess the performance of the random phase approximation (RPA), a parameter-free fifth-rung functional, for reaction energies governed by changes in medium- and long-range noncovalent interactions. Our benchmark data include relative energies of alkane isomers, two sets of isomerization reactions testing intramolecular dispersion, a set of dimers of biological importance, a hierarchy of n-homodesmotic reactions, and the predissociation of a ruthenium-based Grubbs catalyst with bulky ligands. The RPA results are an order of magnitude more accurate than those of popular semilocal functionals such as PBE or B3LYP and more systematic than those of semiempirical functionals parametrized for weak interactions, such as B2PLYP-D or M06-2X. In conclusion, RPA is highly promising for thermochemical applications, particularly if noncovalent interactions are important

Basis Set Convergence of Molecular Correlation Energy Differences Within the Random Phase Approximation

The basis set convergence of energy differences obtained from the random phase approximation (RPA) to the correlation energy is investigated for a wide range of molecular interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark [P. Jurecka, Phys. Chem. Chem. Phys. 8, 1985 (2006)10.1039/b600027d]. The use of very large basis sets ( quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable estimate of the binding energy for these systems. Counterpoise corrected results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramolecular dispersion, isomerization energies, and reaction energies of small organic molecules. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency