A post-Hartree-Fock model of intermolecular interactions

Revisiting the Volmer-Heyrovský mechanism of hydrogen evolution on a nitrogen doped carbon nanotube: constrained molecular dynamics versus the nudged elastic band method

Density functional theory (DFT) based computational electrochemistry has the potential to serve as a tool with predictive power in the rational development and screening of electrocatalysts for renewable energy technologies. It is, however, of paramount importance that simulations are conducted rigorously at a level of theory that is sufficiently accurate in order to obtain physicochemically sensible results. Herein, we present a comparative study of the performance of the static climbing image nudged elastic band method (CI-NEB) vs. DFT based constrained molecular dynamics simulations with thermodynamic integration in estimating activation and reaction (free) energies of the Volmer-Heyrovský mechanism on a nitrogen doped carbon nanotube. Due to cancellation of errors within the CI-NEB calculations, static and dynamic activation barriers are observed to be surprisingly similar, while a substantial decrease in reaction energies is seen upon incorporation of solvent dynamics. This finding is attributed to two competing effects; (1) solvent reorganization that stabilizes the transition and, in particular, the product states with respect to the reactant state and (2) destabilizing entropic contributions due to solvent fluctuations. Our results highlight the importance of explicitly sampling the interfacial solvent dynamics when studying hydrogen evolution at solid-liquid interfaces.

Synthesis of 2,1,3-Benzoxadiazole Derivatives as New Fluorophores-Combined Experimental, Optical, Electro, and Theoretical Study

Herein, we report the synthesis and characterization of fluorophores containing a 2,1,3-benzoxadiazole unit associated with a π-conjugated system (D-π-A-π-D). These new fluorophores in solution exhibited an absorption maximum at around ~419 nm (visible region), as expected for electronic transitions of the π-π* type (ε ~2.7 × 107 L mol-1 cm-1), and strong solvent-dependent fluorescence emission (ΦFL ~0.5) located in the bluish-green region. The Stokes' shift of these compounds is ca. 3,779 cm-1, which was attributed to an intramolecular charge transfer (ICT) state. In CHCl3 solution, the compounds exhibited longer and shorter lifetimes, which was attributed to the emission of monomeric and aggregated molecules, respectively. Density functional theory was used to model the electronic structure of the compounds 9a-d in their excited and ground electronic states. The simulated emission spectra are consistent with the experimental results, with different solvents leading to a shift in the emission peak and the attribution of a π-π* state with the characteristics of a charge transfer excitation. The thermal properties were analyzed by thermogravimetric analysis, and a high maximum degradation rate occurred at around 300°C. Electrochemical studies were also performed in order to determine the band gaps of the molecules. The electrochemical band gaps (2.48-2.70 eV) showed strong correlations with the optical band gaps (2.64-2.67 eV).

Theoretical study of the mechanism behind the site- and enantio-selectivity of C-H functionalization catalysed by chiral dirhodium catalyst

The C-H functionalization is very important for the synthesis of pharmaceuticals and complex natural products. Rhodium carbenoids, obtained when a dirhodium(ii) catalyst containing a crown formed by chiral ligands reacts with diazo compounds with both an electron donating group and an electron withdrawing group, play an important part in controlling site- and enantio-selectivity for functionalization of non-activated C-H bonds. It has earlier been demonstrated that the tertiary C-H bond is more favored to be functionalized inside the crown of the dirhodium catalyst with S-configuration ligands compared with the secondary and primary C-H bonds although the latter possess weaker steric effects. We argue that the higher site- and enantio-selectivity for some types of C-H bond functionalization can be related to intermolecular hydrogen bonding, steric hindrance, and weak interactions when the dirhodium catalyst is interacting with the chiral ligands.

A consistent model for the key complex in chronic beryllium disease

A hypothesis on the structure of the key complex in chronic beryllium disease (CBD) is discussed with respect to the current knowledge on CBD, and with respect to the constraints implied by the coordination chemistry of beryllium and experimental data on the engaged protein complexes. The structure hypothesis is based on the [Be4O]6+ moiety as a coordination center, which is also found in the so called “basic beryllium carboxylates”. The structure of a small molecular model, optimized at the DFT level of theory, is used to compare the structural demands of this coordination center with a structure of the in vitro model of a beryllium immunoprotein complex determined previously by protein crystallography (Clayton & al., Cell 2014, 158, 132). 9Be NMR chemical shielding values, quadrupole coupling constants and asymmetry parameters (η) have been calculated.

Semiclassical Dispersion Corrections Efficiently Improve Multiconfigurational Theory with Short-Range Density-Functional Dynamic Correlation

Multiconfigurational wave functions are known to describe the electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamic correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually limit the ultimate accuracy achievable (most prominently, multireference perturbation theory). At the same time, the computational cost increases dramatically due to the necessity to obtain higher-order reduced density matrices. It is this disproportion that leads us here to propose an MC-srDFT-D hybrid approach of semiclassical dispersion (D) corrections to cover long-range dynamic correlation in a multiconfigurational (MC) wave function theory, which includes short-range (sr) dynamic correlation by density functional theory (DFT) without double counting. We demonstrate that the reliability of this approach is very good (at negligible cost), especially when considering that standard second-order multireference perturbation theory usually overestimates dispersion interactions.

Elastic properties of the molecular crystals of hydrocarbons from first principles calculations

I studied the elastic properties of oligoacenes and polycyclic aromatic hydrocarbons crystals within the framework of density functional theory with van der Waals interactions. The full sets of elastic constants were computed. The computed parameters have good agreement with experimental data. The study of four forms (α, β, γ and dimeric) demonstrate that the stacking of molecules leads to different character of interaction between molecules. Crystals with dimeric form have lower anisotropy. If symmetry of crystals is the same, the increase in number of aromatic rings results in anisotropy increase.

Active site manipulation in MoS2 cluster electrocatalysts by transition metal doping

The development of non-platinum group metal catalysts for the hydrogen evolution reaction (HER) in water electrolyser devices is essential for their widespread and sustainable deployment. In recent years, molybdenum disulfide (MoS2) catalysts have received significant attention as they not only exhibit good electrocatalytic HER activity but also, crucially, acid-stability. However, further performance enhancement is required for these materials to be competitive with Pt and to that end transition metal doping of MoS2 has been explored as a route to further increasing its catalytic activity. In this work, cluster beam deposition was employed to produce controlled cobalt-doped MoS2 clusters (MoS2-Co). We demonstrate that, in contrast to previous observations of performance enhancement in MoS2 resulting from nickel doping (MoS2-Ni), the introduction of Co has a detrimental effect on HER activity. The contrasting behaviours of Ni and Co doping are rationalized by density functional theory (DFT) calculations, which suggest that HER-active surface vacancies are deactivated by combination with Co dopant atoms, whilst their activity is retained, or even partially enhanced, by combination with Ni dopant atoms. Furthermore, the adatom dopant-vacancy combination kinetics appear to be more than three orders of magnitude faster in MoS2-Co than for MoS2-Ni. These findings highlight a fundamental difference in the influence of transition metal dopants on the HER performance of MoS2 electrocatalysts and stress the importance of considering surface atomic defects when predicting their behaviour.

The unusual symmetry of hexafluoro-o-xylene-A microwave spectroscopy and computational study

The rotational constants and quartic centrifugal distortion coefficients of hexafluoro-o-xylene were precisely derived from the 8 GHz to 18 GHz gas phase microwave spectrum. In addition, the rotational constants of all singly substituted ¹³C isotopologues were determined. Instead of the intuitively expected symmetry of C_2v, as in o-xylene, calculations with a variety of methods (B3LYP, CAM-B3LYP, ωB97XD, MP2, and coupled-cluster singles, doubles, and perturbative triples) predict a C₂ symmetry structure in which the two CF₃ groups rotate in opposite directions by about 16°. While the experimental results in this study are not capable of proving this unusual symmetry, they can support future microwave, gas phase electron diffraction or nuclear magnetic resonance studies. From the presented data, a preliminary r₀ structure was determined, reproducing the experimental rotational constants with deviations of no more than 15 kHz. Analysis of the interactions between the two CF₃ groups using an effective fragment potential approach identifies two major contributions to their interaction, due to exchange-repulsion and electrostatic repulsion, with electrostatic repulsion responsible for the barrier at the C_2v geometry.

Physicochemical properties of L- and DL-valine: first-principles calculations

At present, physicochemical properties of amino acid molecular crystals are of the utmost interest. The compounds where molecules have different chirality are the focus of particular interest. This paper, presents a study on the structural and electronic properties of crystalline L- and DL-valine within the framework of density functional theory including van der Waals interactions. The results of this study showed that electronic properties of the two forms of valine are similar at zero pressure. Pressure leads to different responses in these crystals which is manifested as various deformations of molecules. The pressure effect on the infrared spectra and distribution of electron density of L- and DL-valine has been studied.

Theoretical studies of infrared signatures of proton-bound amino acid dimers with homochiral and heterochiral moieties

Proton-bound homochiral and heterochiral dimers, X-H⁺ -X, of five amino acids (X = Ser, Ala, Thr, Phe, and Arg) are investigated theoretically using quantum chemical density functional theory (DFT) calculations and molecular dynamics simulations with the aim to unveil diastereomer-specific mid-infrared (mid-IR) absorption bands in the spectral range of 1000 to 1800 cm^-1 . The theoretical calculations performed in this work imply that all systems, except Ala₂ H⁺ , have distinct mid-IR absorption bands in homochiral and heterochiral configurations, which make them appropriate systems to be studied experimentally with mid-IR spectroscopy. We show that intermolecular interaction with the side chain, in the form of hydrogen bonding or cation-π interaction, is necessary for chiral effects to be present in the mid-IR spectra of proton-bound dimers of amino acids. We also report new conformers for Ala₂ H⁺ , Thr₂ H⁺ , Phe₂ H⁺ , and Arg₂ H⁺ , which were not found in earlier studies of these dimers.

Understanding chemical reactivity using the activation strain model

Understanding chemical reactivity through the use of state-of-the-art computational techniques enables chemists to both predict reactivity and rationally design novel reactions. This protocol aims to provide chemists with the tools to implement a powerful and robust method for analyzing and understanding any chemical reaction using PyFrag 2019. The approach is based on the so-called activation strain model (ASM) of reactivity, which relates the relative energy of a molecular system to the sum of the energies required to distort the reactants into the geometries required to react plus the strength of their mutual interactions. Other available methods analyze only a stationary point on the potential energy surface, but our methodology analyzes the change in energy along a reaction coordinate. The use of this methodology has been proven to be critical to the understanding of reactions, spanning the realms of the inorganic and organic, as well as the supramolecular and biochemical, fields. This protocol provides step-by-step instructions-starting from the optimization of the stationary points and extending through calculation of the potential energy surface and analysis of the trend-decisive energy terms-that can serve as a guide for carrying out the analysis of any given reaction of interest within hours to days, depending on the size of the molecular system.

Three-dimensional docking of alcohols to ketones: an experimental benchmark based on acetophenone solvation energy balances

The two hydrogen bond solvation sites exhibited by the carbonyl group in acetophenone are influenced by alkylation of the methyl group in both the acetophenone and in the prototype solvent methanol, largely due to London dispersion forces. Phenyl docking and alkyl docking preferences can be realized at will by appropriate substitution. In particular, cyclopropylation helps to stabilize the opposite phenyl docking site. In all cases, the energy gap is small enough to allow for a simultaneous detection even under low temperature conditions. This density functional prediction is checked experimentally by jet FTIR spectroscopy and largely confirmed. A spurious out-of-plane solvation preference predicted for cyclopropylphenylketone with tert-butyl alcohol by B3LYP-D3 calculations is not confirmed experimentally. It is unlikely that this discrepancy is due to zero-point energy effects. Instead, the second most stable alkyl-side solvation motif predicted with a more in-plane coordination is found in the jet expansion. Overall, the ability of carbonyl solvation balances to benchmark subtle electronic structure effects for non-covalent interactions without major nuclear motion corrections is supported.

Effect of van der Waals interactions on the adhesion strength at the interface of the hydroxyapatite–titanium biocomposite: a first-principles study

Hydroxyapatite (HAP) is frequently used as biocompatible coating on Ti-based implants. In this context, the HAP-Ti adhesion is of crucial importance. Here, we report ab initio calculations to investigate the influence of Si incorporation into the amorphous calcium-phosphate (a-HAP) structure on the interfacial bonding mechanism between the a-HAP coating and an amorphous titanium dioxide (a-TiO₂) substrate, contrasting two different density functionals: PBE-GGA, and DFT-D3, which are capable of describing the influence of the van der Waals (vdW) interactions. In particular, we discuss the effect of dispersion on the work of adhesion (W_ad), equilibrium geometries, and charge density difference (CDD). We find that replacement of P by Si in a-HAP (a-Si-HAP) with the creation of OH vacancies as charge compensation results in a significant increase in the bond strength between the coating and substrate in the case of using the PBE-GGA functional. However, including the vdW interactions shows that these forces considerably contribute to the W_ad. We show that the difference (W_ad − W_ad(vdW)) is on average more than 1.1 J m⁻² and 0.5 J m⁻² for a-HAP/a-TiO₂ and a-Si-HAP/a-TiO₂, respectively. These results reveal that including vdW interactions is essential for accurately describing the chemical bonding at the a-HAP/a-TiO₂ interface.

Dispersion interactions play a significant role in altering the bonding strength of the hydroxyapatite coating to the substrate, and such interactions is critical in molecular simulations.

The simultaneous recognition mechanism of cations and anions using macrocyclic–iodine structures: insights from dispersion-corrected DFT calculations

The design of structures to recognize ions is growing in recent years. Here, the simultaneous recognition of cations and anions by a macrocycle comprising a simple crown ether and an iodine-triazole unit has been investigated using DFT calculations.

Microhydration of verbenone: how the chain of water molecules adapts its structure to the host molecule

The microsolvation of verbenone (C₁₀H₁₄O)·(H₂O)_n (n = 1, 2, 3) was investigated in a supersonic expansion using a cavity-based Fourier transform microwave spectrometer, in the 2.8–14 GHz frequency range and by computational chemistry.

Intra-octahedral distortion on lamellar potassium niobate K4Nb6O17: a periodic DFT study of structural, electronic and vibrational properties

DFT calculation applied to K₄Nb₆O₁₇ allowed to identify and correlate its electronic and vibrational properties with [NbO₆] intraoctahedral distortion.

Quantum chemical topology and natural bond orbital analysis of M–O covalency in M(OC6H5)4 (M = Ti, Zr, Hf, Ce, Th, Pa, U, Np)

V_XC(M,O): the exchange–correlation metric quantifies covalency between M and O atomic basins in M(OC₆H₅)₄ (M = Ti, Zr, Hf, Ce, Th, Pa, U, Np).

Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm

Through a comprehensive survey, reliable double-hybrid models have been validated and proposed for spin-state energetics in transition metal complexes.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state

Parametrization of PBE+U under the D3 and D3-BJ dispersion corrections to study Fe^II and Fe^III-based Spin Crossover complexes.

Xe⋯OCS: relatively straightforward?

Spectroscopy meets theory in a study of Xe⋯OCS complex: accurate near-equilibrium structures, experimental interaction energies, and CCSD(T)/CBS results presented.

Modeling the Hydrolysis of Iron-Sulfur Clusters

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.

Non-covalent Interactions and Charge Transfer between Propene and Neutral Yttrium-Doped and Pure Gold Clusters

The dopant and size-dependent propene adsorption on neutral gold (Aun ) and yttrium-doped gold (Aun-1 Y) clusters in the n=5-15 size range are investigated, combining mass spectrometry and gas phase reactions in a low-pressure collision cell and density functional theory calculations. The adsorption energies, extracted from the experimental data using an RRKM analysis, show a similar size dependence as the quantum chemical results and are in the range of ≈0.6-1.2 eV. Yttrium doping significantly alters the propene adsorption energies for n=5, 12 and 13. Chemical bonding and energy decomposition analysis showed that there is no covalent bond between the cluster and propene, and that charge transfer and other non-covalent interactions are dominant. The natural charges, Wiberg bond indices, and the importance of charge transfer all support an electron donation/back-donation mechanism for the adsorption. Yttrium plays a significant role not only in the propene binding energy, but also in the chemical bonding in the cluster-propene adduct. Propene preferentially binds to yttrium in small clusters (n<10), and to a gold atom at larger sizes. Besides charge transfer, relaxation also plays an important role, illustrating the non-local effect of the yttrium dopant. It is shown that the frontier molecular orbitals of the clusters determine the chemical bonding, in line with the molecular-like electronic structure of metal clusters.

Old Concepts, New Application - Additive-Free Hydrogenation of Nitriles Catalyzed by an Air Stable Alkyl Mn(I) Complex

An efficient additive-free manganese-catalyzed hydrogenation of nitriles to primary amines with molecular hydrogen is described. The pre-catalyst, a well-defined bench-stable alkyl bisphosphine Mn(I) complex fac-[Mn(dpre)(CO)3(CH3)] (dpre=1,2-bis(di-n-propylphosphino)ethane), undergoes CO migratory insertion into the manganese-alkyl bond to form acyl complexes which upon hydrogenolysis yields the active coordinatively unsaturated Mn(I) hydride catalyst [Mn(dpre)(CO)2(H)]. A range of aromatic and aliphatic nitriles were efficiently and selectively converted into primary amines in good to excellent yields. The hydrogenation of nitriles proceeds at 100 °C with a catalyst loading of 2 mol % and a hydrogen pressure of 50 bar. Mechanistic insights are provided by means of DFT calculations.

Gaussian Process-Based Refinement of Dispersion Corrections

We employ Gaussian process (GP) regression to adjust for systematic errors in D3-type dispersion corrections. We refer to the associated, statistically improved model as D3-GP. It is trained on differences between interaction energies obtained from PBE-D3(BJ)/ma-def2-QZVPP and DLPNO-CCSD(T)/CBS calculations. We generated a data set containing interaction energies for 1248 molecular dimers, which resemble the dispersion-dominated systems contained in the S66 data set. Our systems represent not only equilibrium structures but also dimers with various relative orientations and conformations at both shorter and longer distances. A reparametrization of the D3(BJ) model based on 66 of these dimers suggests that two of its three empirical parameters, a₁ and s₈, are zero, whereas a₂ = 5.6841 bohr. For the remaining 1182 dimers, we find that this new set of parameters is superior to all previously published D3(BJ) parameter sets. To train our D3-GP model, we engineered two different vectorial representations of (supra-)molecular systems, both derived from the matrix of atom-pairwise D3(BJ) interaction terms: (a) a distance-resolved interaction energy histogram, histD3(BJ), and (b) eigenvalues of the interaction matrix ordered according to their decreasing absolute value, eigD3(BJ). Hence, the GP learns a mapping from D3(BJ) information only, which renders D3-GP-type dispersion corrections comparable to those obtained with the original D3 approach. They improve systematically if the underlying training set is selected carefully. Here, we harness the prediction variance obtained from GP regression to select optimal training sets in an automated fashion. The larger the variance, the more information the corresponding data point may add to the training set. For a given set of molecular systems, variance-based sampling can approximately determine the smallest subset being subjected to reference calculations such that all dispersion corrections for the remaining systems fall below a predefined accuracy threshold. To render the entire D3-GP workflow as efficient as possible, we present an improvement over our variance-based, sequential active-learning scheme [ J. Chem. Theory Comput. 2018 , 14 , 5238 ]. Our refined learning algorithm selects multiple (instead of single) systems that can be subjected to reference calculations simultaneously. We refer to the underlying selection strategy as batchwise variance-based sampling (BVS). BVS-guided active learning is an essential component of our D3-GP workflow, which is implemented in a black-box fashion. Once provided with reference data for new molecular systems, the underlying GP model automatically learns to adapt to these and similar systems. This approach leads overall to a self-improving model (D3-GP) that predicts system-focused and GP-refined D3-type dispersion corrections for any given system of reference data.

Cost-effective density functional theory (DFT) calculations of equilibrium isotopic fractionation in large organic molecules

The application of stable isotopes to address a wide range of biochemical, microbiological and environmental problems is hindered by the experimental difficulty and the computational cost of determining equilibrium isotopic fractionations (EIF) of large organic molecules. Here, we evaluate the factors that impact the accuracy of computed EIFs and develop a framework for cost-effective and accurate computation of EIFs by density functional theory (DFT). We generated two benchmark databases of experimentally determined EIFs, one for H isotopes and another for the isotopes of the heavy atoms C, N and O. The accuracy of several DFT exchange-correlation functionals in calculating EIFs was then evaluated by comparing the computational results to these experimental datasets. We find that with the def2-TZVP basis set, O3LYP had the lowest mean absolute deviation (21‰ and 3.9‰ for the isotopic fractionation of H and the heavier atoms, respectively), but the GGA/meta-GGA functionals τHCTH_D3BJ, τHCTH and HCTH have similar performances (22‰ and 4.1‰, respectively, for τHCTH_D3BJ). Leveraging the good performance of computationally efficient functionals, we provide a robust, practical, experimentally validated framework for using DFT to accurately predict EIFs of large organic molecules, including uncertainty estimates.

2,3-(Dibenzimidazol-2-yl)quinoxalines: Unexpected Dynamical Effect on Steady-State Electronic Absorption Spectra

We report on the electronic absorption spectra, conformational behavior, and intra- and intermolecular hydrogen bonds of 2,3-(dibenzimidazol-2-yl)-quinoxaline (DBIQ). The experimentally found strong solvent dependence of the absorption spectra of DBIQ solutions cannot be assigned to electronic excitations of the equilibrium ground-state DBIQ structure. Extended consideration including the nonequilibrium structures within the framework of ab initio molecular dynamics (MD) revealed the importance of torsion molecular motions not covered by the static case. The strong impact of solute-solvent hydrogen bonding on stabilization of these nonequilibrium structures and on conformational composition of DBIQ was demonstrated. A presence of twisted nonplanar geometries along the whole MD trajectory was shown to drastically influence not only energies but also characters of electronic excitations, resulting in a change of local π-π* character in a solution of 1,2-dichloroethane to charge-transfer character in polar dimethylsulfoxide.

New scaling relations to compute atom-in-material polarizabilities and dispersion coefficients: part 1. Theory and accuracy

Polarizabilities and London dispersion forces are important to many chemical processes. Force fields for classical atomistic simulations can be constructed using atom-in-material polarizabilities and C n (n = 6, 8, 9, 10…) dispersion coefficients. This article addresses the key question of how to efficiently assign these parameters to constituent atoms in a material so that properties of the whole material are better reproduced. We develop a new set of scaling laws and computational algorithms (called MCLF) to do this in an accurate and computationally efficient manner across diverse material types. We introduce a conduction limit upper bound and m-scaling to describe the different behaviors of surface and buried atoms. We validate MCLF by comparing results to high-level benchmarks for isolated neutral and charged atoms, diverse diatomic molecules, various polyatomic molecules (e.g., polyacenes, fullerenes, and small organic and inorganic molecules), and dense solids (including metallic, covalent, and ionic). We also present results for the HIV reverse transcriptase enzyme complexed with an inhibitor molecule. MCLF provides the non-directionally screened polarizabilities required to construct force fields, the directionally-screened static polarizability tensor components and eigenvalues, and environmentally screened C6 coefficients. Overall, MCLF has improved accuracy compared to the TS-SCS method. For TS-SCS, we compared charge partitioning methods and show DDEC6 partitioning yields more accurate results than Hirshfeld partitioning. MCLF also gives approximations for C8, C9, and C10 dispersion coefficients and quantum Drude oscillator parameters. This method should find widespread applications to parameterize classical force fields and density functional theory (DFT) + dispersion methods.

Dispersion Energy from Local Polarizability Density

A simple nonlocal functional for calculation of dispersion energies is proposed. Compared to a similar formula used earlier, we introduced a regularization to remove its singularities and used a dynamic polarizability density similar to those in the so-called van der Waals density functionals. The performance of the new functional is tested on dispersion energies for a set of representative dimers, and it is found that it is significantly more accurate than published nonlocal functionals.

Evaluating Transition Metal Barrier Heights with the Latest Density Functional Theory Exchange-Correlation Functionals: The MOBH35 Benchmark Database

A new database of transition metal reaction barrier heights (MOBH35) is presented. Benchmark energies (forward and reverse barriers and reaction energy) are calculated using DLPNO-CCSD(T) extrapolated to the complete basis set limit using a Weizmann-1-like scheme. Using these benchmark energies, the performance of a wide selection of density functional theory (DFT) exchange-correlation functionals, including the latest from the Martin, Truhlar, and Head-Gordon groups, is evaluated. It was found, using the def2-TZVPP basis set, that the ωB97M-V (MAD 1.7 kcal/mol), ωB97M-D3BJ (MAD 1.9 kcal/mol), ωB97X-V (MAD 2.0 kcal/mol), and revTPSS0-D4 (MAD 2.2 kcal/mol) hybrid functionals are recommended. The double-hybrid functionals B2K-PLYP (MAD 1.7 kcal/mol) and revDOD-PBEP86-D4 (MAD 1.8 kcal/mol) also performed well, but this has to be balanced by their increased computational cost.