Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

Kohn-Sham Density Matrix and the Kernel Energy MethodSCI北大核心CSCD

The kernel energy method(KEM) has been shown to provide fast and accurate molecular energy calculations for molecules at their equilibrium geometries.KEM breaks a molecule into smaller subsets,called kernels,for the purposes of calculation.The results from the kernels are summed according to an expression characteristic of KEM to obtain the full molecule energy.A generalization of the kernel expansion to density matrices provides the full molecule density matrix and orbitals.In this study,the kernel expansion for the density matrix is examined in the context of density functional theory(DFT) Kohn-Sham(KS) calculations.A kernel expansion for the one-body density matrix analogous to the kernel expansion for energy is defined,and is then converted into a normalizedprojector by using the Clinton algorithm.Such normalized projectors are factorizable into linear combination of atomic orbitals(LCAO) matrices that deliver full-molecule Kohn-Sham molecular orbitals in the atomic orbital basis.Both straightforward KEM energies and energies from a normalized,idempotent density matrix obtained from a density matrix kernel expansion to which the Clinton algorithm has been applied are compared to reference energies obtained from calculations on the full system without any kernel expansion.Calculations were performed both for a simple proof-of-concept system consisting of three atoms in a linear configuration and for a water cluster consisting of twelve water molecules.In the case of the proof-of-concept system,calculations were performed using the STO-3 G and6-31 G(d,p) bases over a range of atomic separations,some very far from equilibrium.The water cluster was calculated in the 6-31 G(d,p) basis at an equilibrium geometry.The normalized projector density energies are more accurate than the straightforward KEM energy results in nearly all cases.In the case of the water cluster,the energy of the normalized projector is approximately four times more accurate than the straightforward KEM energy result.The KS density matrices of this study are applicable to quantum crystallography.     

An improved theoretical approach to the empirical corrections of density functional theory

An empirical correction to density functional theory (DFT) has been developed in this study. The approach, called correlation corrected atomization–dispersion (CCAZD), involves short- and long-range terms. Short-range correction consists of bond (1,2-) and angle (1,3-) interactions, which remedies the deficiency of DFT in describing the proto-branching stabilization effects. Long-range correction includes a Buckingham potential function aiming to account for the dispersion interactions. The empirical corrections of DFT were parameterized to reproduce reported ΔH _f values of the training set containing alkane, alcohol and ether molecules. The ΔH _f of the training set molecules predicted by the CCAZD method combined with two different DFT methods, B3LYP and MPWB1K, with a 6-31G* basis set agreed well with the experimental data. For 106 alkane, alcohol and ether compounds, the average absolute deviations (AADs) in ΔH _f were 0.45 and 0.51 kcal/mol for B3LYP- and MPWB1K-CCAZD, respectively. Calculations of isomerization energies, rotational barriers and conformational energies further validated the CCAZD approach. The isomerization energies improved significantly with the CCAZD treatment. The AADs for 22 energies of isomerization reactions were decreased from 3.55 and 2.44 to 0.55 and 0.82 kcal/mol for B3LYP and MPWB1K, respectively. This study also provided predictions of MM4, G3, CBS-QB3 and B2PLYP-D for comparison. The final test of the CCAZD approach on the calculation of the cellobiose analog potential surface also showed promising results. This study demonstrated that DFT calculations with CCAZD empirical corrections achieved very good agreement with reported values for various chemical reactions with a small basis set as 6-31G*.     

Ag-based bimetallic clusters as catalysts for p-nitrophenol reduction by glycerol: A DFT investigation

Reducing p-nitrophenol (PNP) to p-aminophenol is an industrially relevant synthesis. Nevertheless, only a few heterogeneous catalysts have been evaluated for the reduction of PNP by glycerol. Appropriate quantum computational studies can screen potential catalysts for this crucial green reaction. The present research investigates the catalytic activities of Pd@Ag and Ni@Ag core-shell nanogeometries toward PNP reduction by glycerol through density functional theory (DFT) calculations. The central atom of a geometry-optimized 13-atom Ag cluster was replaced by Pd and Ni atoms to create the core-shell morphologies. The interaction energies of PNP and glycerol with each of the (metal/bimetallic) clusters were evaluated by DFT calculations to find the best PNP and glycerol molecule orientation with the respective bimetallic cluster. Electrostatic potential surface and natural bond orbital analyses were performed to study the charge distribution and transfer between atomic orbitals. The frequencies of vibrational modes in isolated PNP/glycerol structures were compared to those when these molecules were in the presence of the different metal clusters to infer the effect of the interactions. All performed analyses indicated improved catalytic activity toward PNP reduction by glycerol upon Ni-doping of the Ag₁₃ cluster.     

DFT and MO calculations of atomic and molecular chemisorption energies on surface cluster models

Summary Density functional theory (DFT) (including gradient corrections) and MCPF calculations have been performed for atomic (H, C, N, O) and molecular CH _x (x = 1–3) chemisorption on cluster models of different sites of the Cu(100) surface. The DFT and MCPF results are in good agreement once the important effects of core-valence correlation have been accounted for in the MCPF calculations by including contributions from a core polarization potential (CPP); in the DFT approach the core-valence correlation is obtained directly from the total density using the functional. Very large effects on the four-fold hollow site binding energy from core-valence correlation are found for C, N and CH. Several different DFT functionals were employed and compared in the calculations.     

Binding interactions and Raman spectral properties of pyridine interacting with bimetallic silver-gold clusters.

The binding interactions between pyridine and bimetallic silver-gold clusters are investigated using density functional theory (DFT). The binding energies of pyridine-bimetallic cluster complexes indicate that the bonding depends strongly on the binding site (Au or Ag atom) and bonding molecular orbitals in a given configuration. The donation of the lone-pair electrons of the nitrogen of pyridine to an appropriate unoccupied orbital of each metal cluster plays an important role. The low-lying excited states and charge-transfer states of four stable complexes of interest are calculated on the basis of a time-dependent DFT method. In nonresonance Raman scattering processes, the influence of binding interactions on the relative Raman intensity of totally symmetric pyridine vibrational modes is discussed. These calculated relative Raman intensities are compared with observed surface-enhanced Raman spectra of pyridine adsorbed on silver-gold alloy surfaces.     

Grid-based energy density analysis: implementation and assessment

Grid-based energy density analysis (grid-EDA) that decomposes the total energy into atomic energies by a space-partitioning function is proposed. The kinetic energy, nuclear attraction, and exchange-correlation functional are evaluated on grid points and are split into atomic contributions. To reduce numerical errors in the conventional scheme of numerical integration, the electronic Coulomb and HF exchange interactions are evaluated by the pseudospectral method, which was first applied to an ab initio method by Friesner [Chem. Phys. Lett. 116, 39 (1985)], and are decomposed into atomic contributions. Grid-EDA using the pseudospectral method succeeds in ensuring less than 1 kcalmol error in total energies for small molecules and providing reliable atomic energy contributions for the problematic lithium cluster, which exhibits a strong basis-set dependence for Mulliken-type EDA. Also, site-dependent atomization energies are estimated by grid-EDA for cluster models such as Li(48), C(41)H(60), and Mg(32)O(32). Grid-EDA reveals that these models imitate crystal environments reasonably because atomization energies estimated from the inner atoms of the models are close to the experimental cohesive energies.     

Direct calculation of electron transfer parameters through constrained density functional theory

It is shown that constrained density functional theory (DFT) can be used to access diabatic potential energy surfaces in the Marcus theory of electron transfer, thus providing a means to directly calculate the driving force and the inner-sphere reorganization energy. We present in this report an analytic expression for the forces in constrained DFT and their implementation in geometry optimization, a prerequisite for the calculation of electron transfer parameters. The method is then applied to study the symmetric mixed-valence complex tetrathiafulvalene-diquinone radical anion, which is observed experimentally to be a Robin-Day class II compound but found by DFT to be in class III. Constrained DFT avoids this pitfall of over-delocalization and provides a way to find the charge-localized structure. In another application, driving forces and inner-sphere reorganization energies are calculated for the charge recombination (CR) reactions in formanilide-anthraquinone (FA-AQ) and ferrocene-formanilide-anthraquinone (Fc-FA-AQ). While the two compounds have similar reorganization energies, the driving force in FA-AQ is 1 eV larger than in Fc-FA-AQ, in agreement with experimental observations and supporting the experimental conclusion that the anomalously long-lived FA-AQ charge-separated state arises because the electron transfer is in the Marcus inverted region.     

Comprehensive analysis of the hydrogen bond network morphology and OH stretching vibrations in protonated methanol-water mixed clusters, H(+)(MeOH)1(H2O)n (n = 1-8)

Density functional theory (DFT) calculations of protonated methanol-water mixed clusters, H (+)(MeOH) 1(H 2O) n ( n = 1-8), were extensively carried out to analyze the hydrogen bond structures of the clusters. Various structural isomers were energy optimized, and their relative energies with zero point energy corrections and temperature dependence of the free energies were examined. Coexistence of different morphological isomers was suggested. Infrared spectra were simulated on the basis of the optimized structures. The infrared spectra were also experimentally measured for n = 3-9 in the OH stretching vibrational region. The observed broad bands in the hydrogen-bonded OH stretch region were assigned in comparison with the simulations. From the DFT calculations, the preferential proton location was also investigated. Clear correlations between the excess proton location and the cluster morphology were found.     

Interactions in large, polyaromatic hydrocarbon dimers: application of density functional theory with dispersion corrections

The interactions within two models for graphene, coronene and hexabenzocoronene (HBC), and (H 3C(CH 2) 5) 6-HBC, a synthesizable model for asphaltenes, were studied using density functional theory (DFT) with dispersion corrections. The corrections were implemented using carbon atom-centered effective core-type potentials that were designed to correct the erroneous long-range behavior of several DFT methods. The potentials can be used with any computational chemistry program package that can handle standard effective core potential input, without the need for software modification. Testing on a set of common noncovalently bonded dimers shows that the potentials improve calculated binding energies by factors of 2-3 over those obtained without the potentials. Binding energies are predicted to within ca. 15%, and monomer separations to within ca. 0.1 A, of high-level wave function data. The application of the present approach predicts binding energies and structures of the coronene dimer that are in excellent agreement with the results of other DFT methods in which dispersion is taken into account. Dimers of HBC show extensive binding in pi-stacking arrangements, with the largest binding energy, 44.8 kcal/mol, obtained for a parallel-displaced structure. This structure is inline with the published crystal structure. Conformations in which the monomers are perpendicular to one another are much more weakly bound and have binding energies less than 10 kcal/mol. For dimers of (H 3C(CH 2) 5) 6-HBC, which contain 336 atoms, we find that a slipped-parallel structure with C s symmetry has a binding energy of 52.4 kcal/mol, 8.9 kcal/mol lower than that of a bowl-like, C 6 v -symmetric structure.     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

Treating dispersion effects in extended systems by hybrid MP2:DFT calculations--protonation of isobutene in zeolite ferrierite

We propose use of a hybrid method to study problems that involve both bond rearrangements and van-der-Waals interactions. The method combines second-order M?ller-Plesset perturbation theory (MP2) calculations for the reaction site with density functional theory (DFT) calculations for a large system under periodic boundary conditions. Hybrid MP2:DFT structure optimisation for a cluster embedded in the periodic model is the first of three steps in a multi-level approach. The second step is extrapolation of the MP2 energy to the complete basis set limit. The third step is extrapolating the high-level (MP2) correction to the limiting case of the full periodic structure. This is done by calculating the MP2 correction for a series of cluster models of increasing size, fitting an analytic expression to these energy corrections, and applying the fitted expression to the full periodic structure. We assume that, up to a constant, the high-level correction is described by a damped dispersion expression. Combining the results of all three steps yields an estimate of the MP2 reaction energy for the full periodic system at the complete basis set level. The method is designed for a reaction between a small or medium sized substrate molecule and a very large chemical system. For adsorption of isobutene in zeolite H-ferrierite, the energies obtained for the formation of different structures, the pi-complex, the isobutoxide, the tert-butoxide, and the tert-butyl carbenium ion, are -78, -73, -48, and -21 kJ mol(-1), respectively. This corresponds to corrections of the pure DFT (PBE functional) results by -62, -70, -67, and -29 kJ mol(-1), respectively. Hence, the MP2 corrections are substantial and, perhaps more importantly, not the same for the different hydrocarbon species in the zeolite. Coupled-cluster (CCSD(T)) calculations change the MP2 energies by -4 kJ mol(-1) (tert-butyl cation) or less (below +/-1 kJ mol(-1) for the other species).     

A New Many-Body Expansion Scheme for Atomic Clusters: Application to Nitrogen Clusters

Although the many-body expansion (MBE) approach is widely applied to estimate the energy of large systems containing weak interactions, it is inapplicable to calculating the energies of covalent or metal clusters. In this work, we propose an interaction many-body expansion (IMBE) to calculate the energy of atomic clusters containing covalent bonds. In this approach, the energy of a system is expressed as the sum of the energy of atoms and the interaction energy between the atom and its surrounding atoms. The IMBE method is first applied to calculate the energies of nitrogen clusters, in which the interatomic interactions are truncated to four-body terms. The results show that the IMBE approach could significantly reduce the energy error for nitrogen clusters compared with the traditional MBE method. The weak size and structure dependence of the IMBE error with respect to DFT calculations indicates the IMBE method has good potential application in estimating energy of large covalent systems.     

CeO2 catalysed conversion of CO, NO2 and NO from first principles energetics

First principles calculations using density functional theory with corrections for on-site Coulomb interactions (DFT + U) are presented in which we compute the energy for the conversion of CO to CO(2), NO(2) to NO and NO to N(2) over ceria surfaces. The surface sensitivity is discussed on the basis of the vacancy formation energies.     

Nano-jewellery: C5Au12--a gold-plated diamond at molecular level

A mixed carbon-metal cluster is designed by combining the tetrahedral C(5) radical (with a central atom-the skeleton of the C(5)H(12) molecule) and the spherical Au(12) layer (the external atomic shell of the Au(13) cluster). The C(5)Au(12) cluster and its negative and positive ionic derivatives, C(5)Au(12)(+/-), are investigated ab initio (DFT) in terms of optimized structures and relative energies of a few spin-states, for the icosahedral-like and octahedral-like isomers. The cluster is predicted to be generally more stable in its octahedral shape (similar to C(5)H(12)) which prevails for the negative ion and may compete with the icosahedral shape for the neutral system and positive ion. Adiabatic ionization energies (AIE) and electron affinities (AEA) of C(5)Au(12), vertical electron-detachment (VDE) energies of C(5)Au(12)(-), and vertical ionization and electron-attachment energies (VIE, VEA) of C(5)Au(12) are calculated as well, and compared with those for the corresponding isomers of the Au(13) cluster. The AIE and VIE values are found to be close for the two systems, while the AEA and VDE values are significantly reduced for the radical-based species. A simple fragment-based model is proposed for the decomposition of the total interaction into carbon-gold and gold-gold components.     

羰基镍簇Ni(CO)n(n=1-4)的结构和Ni-CO键解离性质的密度泛函理论研究

在密度泛函理论框架下, 应用不同泛函计算了配合物Ni(CO)n(n=1~4)的平衡几何构型和振动频率. 考察了泛函和基组重叠误差对预测Ni—CO键解离能的影响. 计算结果表明, 用杂化泛函能得到与实验一致的优化几何构型和较合理的振动频率. 对Ni(CO)n(n=2~4)体系, 用“纯”泛函, 如BP86和BPW91, 可得到与CCSD(T)更符合、 并与实验值接近的解离能. 当解离产物出现单个金属原子或离子(如金属羰基配合物的完全解离)时, BSSE校正项的计算中应保持金属部分的电子结构一致. 只有考虑配体基组和不考虑配体基组两种情况下金属的电子构型与配合物中金属的构型一致时, 才能得到合理的BSSE校正, 从而预测合理的解离能.     

Bonding and excitation in COCu(001) from a cluster model and density functional treatments

The bonding properties and charge distributions of the COCu(001) system have been studied within density functional theory (DFT) with several density functionals. A Cu18(9,4,5)CO three layer cluster was found to give bond distances and energies in agreement with previous experimental and theoretical results for low coverage systems, provided the atomic basis set includes diffuse orbitals and d-orbitals at the Cu atoms. Charge distributions give insight on the nature of the localized adsorbate bonding. Time-dependent DFT results on excitation energies and on transition and average electric dipoles, relevant to photodesorption, are also presented.     

Binding energies of five molecular pincers calculated by explicit and implicit solvent models

Molecular pincers or tweezers are designed to hold and release the target molecule. Potential applications involve drug distribution in medicine, environment technologies, or microindustrial techniques. Typically, the binding is dominated by van der Waals forces. Modeling of such complexes can significantly enhance their design; yet obtaining accurate complexation energies by theory is difficult. In this study, density functional theory (DFT) computations combined with dielectric continuum solvent model are compared with the potential of mean force approach using umbrella sampling and the weighted histogram analysis method (WHAM) with molecular dynamics (MD) simulations. For DFT, functional and basis set effects are discussed. The computed results are compared to experimental data based on NMR spectroscopic measurements of five synthesized tweezers based on the Tröger's basis. Whereas the DFT computations correctly provided the observed trends in complex stability, they failed to produce realistic magnitudes of complexation energies. Typically, the binding was overestimated by DFT if compared to experiment. The simpler semiempirical PM6‐DH2X scheme proposed lately yielded better magnitudes of the binding energies than DFT but not the right order. The MD‐WHAM simulations provided the most realistic Gibbs binding energies, although the approximate MD force fields were not able to reproduce completely the ordering of relative stabilities of model complexes found by NMR. Yet the modeling provides interesting insight into the complex geometry and flexibility and appears as a useful tool in the tweezers' design. © 2012 Wiley Periodicals, Inc.