The effects of higher orders of perturbation theory on the correlation energy of atoms and bonds in molecules

We examine, for the first time, the effects of higher orders of Møller–Plesset perturbation theory on the individual atoms within a molecule and the bonds between them, via the topological energy partitioning method of interacting quantum atoms. In real terms (i.e., not by absolute value) MP3 decreases the correlation energy of a bond, and MP4SDQ also decreases the energy of the atoms at either end of the bond. In addition, we investigated long‐range through‐space dispersive effects on a H₂ oligomer. Overall, MP3 is the largest correction to the correlation energy, and most of that energy is allocated to chemical bonds, reducing their values in actual terms. The MP4SDQ bond correlation correction, despite being relatively small, tends to have two effects: (i) for small or negative correlation energies MP4SDQ tends to decrease the bond correlation values even more, and (ii) for large (positive) bond correlation energies MP4SDQ tends to restore the bond correlation energies from the MP3 back toward the MP2 values. Furthermore, each individual part of a molecule or complex (atom or bond) has a specific convergence pattern for the MPn series: through‐space interactions converge at MP2 but bonds converge at MP3 level. The atomic correlation energy appears to head toward convergence at the MP4 level.     

Interaction energy of a water molecule with a single-layer graphitic surface modeled by hydrogen- and fluorine-terminated clusters

In ab initio calculations a finite graphitic cluster model is often used to approximate the interaction energy of a water molecule with an infinite single-layer graphitic surface (graphene). In previous studies, the graphitic cluster model is a collection of fused benzene rings terminated by hydrogen atoms. In this study, the effect of using fluorine instead of hydrogen atoms for terminating the cluster model is examined to clarify the role of the boundary. The interaction energy of a water molecule with the graphitic cluster was computed using ab initio methods at the MP2 level of theory and with the 6-31G(d = 0.25) basis set. The interaction energy of a water molecule with graphene is estimated by extrapolation of two series of increasing size graphitic cluster models (C(6n2)H(6n) and C(6n2)F(6n), n = 1-3). Two fixed orientations of water molecule are considered: (a) both hydrogen atoms of water pointing toward the cluster (mode A) and (b) both hydrogen atoms of water pointing away from the cluster (mode B). The interaction energies for water mode A are found to be -2.39 and -2.49 kcal/mol for C(6n2)H(6n) and C(6n2)F(6n) cluster models, respectively. For water mode B, the interaction energies are -2.32 and -2.44 kcal/mol for C(6n2)H(6n) and C(6n2)F(6n) cluster models, respectively.     

A Comparison of the Interacting Quantum Atoms (IQA) Analysis of the Two-Particle Density-Matrices of MP4SDQ and CCSD

Within the quantum topological energy partitioning method called Interacting Quantum Atoms (IQA) we transition from Møller-Plesset (MP4SDQ) to CCSD in calculating intra- and interatomic electron correlation energies for a set of hydrides, diatomics, a few simple molecules and non-covalently bonded complexes, using the uncontracted basis set 6-31++G(2d,2p). CCSD-IQA allows a more rigorous analysis of atomic electron correlation than that offered by Møller-Plesset, which returns IQA contributions that are identical to Hartree–Fock counterparts except for two-electron terms. The CCSD-IQA analysis returns bond and other interatomic correlation energies that are typically much larger in magnitude than the MP4SDQ values. Crisp patterns of energy transferability are detected in water clusters, both for intra-atomic and interatomic correlation energies. CCSD determines that the intra-atomic correlation energy of an oxygen drops by 15 kJ · mol^–1 for donating a hydrogen and by 25 kJ · mol^–1 for accepting a hydrogen.     

Cooperative enhancement of water binding to crownophane by multiple hydrogen bonds: analysis by high level ab initio calculations.

The intermolecular interaction energy of the model system of the water-crownophane complex was analyzed. The water molecule has four hydrogen bonds, with the two hydrogen-donating phenolic hydroxy groups and two hydrogen-accepting oxygen atoms of the poly-oxyethylene chain of the crownophane in the complex. The MP2/6-311G(2d,2p) level calculations of the model system of the complex (hydrogen donating unit + hydrogen accepting unit + water) indicate that the binding energy of the water is 21.85 kcal/mol and that the hydrogen bond cooperativity increases the binding energy as much as 3.67 kcal/mol. The calculated interaction energies depend on the basis set, while the basis set dependence of the cooperative increment is negligible. Most of the cooperative increment is covered by the HF level calculation, which suggests that the major source of the hydrogen bond cooperativity in this system has its origin in induction. The BLYP/6-311G** and PW91/6-311G** level interaction energies of the model system are close to the MP2/6-311G** interaction energies, which suggests that the DFT calculations with these functionals are useful methods to evaluated the interactions of hydrogen bonded systems.     

N-甲替甲酰胺-水的氢键团簇的从头计算研究

The N-methylformamide（NMF）-water clusters were studied by ab inito calculations at MP2/6-31+G** and MP2 / 6-311 ++ G（d,p）levels. The equilibrium geometries and the dissociation channels and dissociation energies of both neutral and ionic NMF-H2O clusters are presented. For N-methylformamide,cis-form has lower energy than trans-form. In NMFH+,the proton prefers to link with the O atom of N-methylformamide. The results show that both cis- and trans- form of NMF can form a linear hydrogen bond with water. Although the energy of trans-NMF is higher than cis-NMF,trans-form exits more stably because it can form a double hydrogen bond with water. After the ionization of the NMF-H2O cluster,both the cis- and the trans-form will produce protonated products.     

Vacuum ultraviolet (VUV) photoionization of small water clusters

Tunable vacuum ultraviolet (VUV) photoionization studies of water clusters are performed using 10-14 eV synchrotron radiation and analyzed by reflectron time-of-flight (TOF) mass spectrometry. Photoionization efficiency (PIE) curves for protonated water clusters (H2O)(n)H+ are measured with 50 meV energy resolution. The appearance energies of a series of protonated water clusters are determined from the photoionization threshold for clusters composed of up to 79 molecules. These appearance energies represent an upper limit of the adiabatic ionization energy of the corresponding parent neutral water cluster in the supersonic molecular beam. The experimental results show a sharp drop in the appearance energy for the small neutral water clusters (from 12.62 +/- 0.05 to 10.94 +/- 0.06 eV, for H2O and (H2O)4, respectively), followed by a gradual decrease for clusters up to (H2O)23 converging to a value of 10.6 eV (+/-0.2 eV). The dissociation energy to remove a water molecule from the corresponding neutral water cluster is derived through thermodynamic cycles utilizing the dissociation energies of protonated water clusters reported previously in the literature. The experimental results show a gradual decrease of the dissociation energy for removal of one water molecule for small neutral water clusters (3 相似文献   

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. I. Estimates of MP2/CBS binding energies and comparison with empirical potentials

We report estimates of complete basis set (CBS) limits at the second-order M?ller-Plesset perturbation level of theory (MP2) for the binding energies of the lowest-lying isomers within each of the four major families of minima of (H(2)O)(20). These were obtained by performing MP2 calculations with the family of correlation-consistent basis sets up to quadruple zeta quality, augmented with additional diffuse functions (aug-cc-pVnZ, n=D, T, Q). The MP2/CPS estimates are -200.1 (dodecahedron, 30 hydrogen bonds), -212.6 (fused cubes, 36 hydrogen bonds), -215.0 (face-sharing pentagonal prisms, 35 hydrogen bonds), and -217.9 kcal/mol (edge-sharing pentagonal prisms, 34 hydrogen bonds). The energetic ordering of the various (H(2)O)(20) isomers does not follow monotonically the number of hydrogen bonds as in the case of smaller clusters such as the different isomers of the water hexamer. The dodecahedron lies ca. 18 kcal/mol higher in energy than the most stable edge-sharing pentagonal prism isomer. The TIP4P, ASP-W4, TTM2-R, AMOEBA, and TTM2-F empirical potentials also predict the energetic stabilization of the edge-sharing pentagonal prisms with respect to the dodecahedron, albeit they universally underestimate the cluster binding energies with respect to the MP2/CBS result. Among them, the TTM2-F potential was found to predict the absolute cluster binding energies to within <1% from the corresponding MP2/CBS values, whereas the error for the rest of the potentials considered in this study ranges from 3% to 5%.     

Proton transfer and autoionization in HNO3·HCl·(H2O)n particles

The structure and spectroscopic properties of clusters of HNO(3)·HCl·(H(2)O)(n), with n = 1 to 6, have been calculated at the MP2/aug-cc-pVDZ level of theory. Altogether 22 different clusters have been found as stable structures, with minima in their potential energy surfaces. The clusters can be grouped in families with the same number of water molecules, and with close aggregation energies within each family. The addition of each new water molecule increments the aggregation energy of the clusters by a nearly constant value of 76.2 ± 0.1 Hartree. The proton transfer parameter and the coordination number of HNO(3) and HCl in each cluster have been evaluated, and the wavenumber shifts for the X(-)-H(+) vibration from the corresponding mode in the isolated molecules have also been predicted. These values allow classification of the acidic species in the clusters into three types, characterized by the strength of the hydrogen bond and the degree of ionization. A correspondence is found between the coordination number of HNO(3) and the magnitude of the X(-)-H(+) vibrational shift.     

Guided ion-beam studies of the kinetic-energy-dependent reactions of Co(n)+ (n = 2-16) with D2: cobalt cluster-deuteride bond energies

The kinetic-energy-dependent cross sections for the reactions of Co(n)+ (n = 2-16) with D2 are measured as a function of kinetic energy over a range of 0-8 eV in a guided ion-beam tandem mass spectrometer. The observed products are Co(n) D+ for all clusters and Co(n)D2+ for n = 4,5,9-16. Reactions for the formation of Co(n)D+ (n = 2-16) and Co9D2+ are observed to exhibit thresholds, whereas cross sections for the formation of Co(n)D2+ (n = 4,5,10-16) exhibit exothermic reaction behavior. The Co(n)+-D bond energies as a function of cluster size are derived from the threshold analysis of the kinetic-energy dependence of the endothermic reactions and are compared to previously determined metal-metal bond energies, D0(Co(n)+-Co). The bond energies of Co(n)+-D generally increase as the cluster size increases, and roughly parallel those for Co(n)+-Co for clusters n > or = 4. These trends are explained in terms of electronic and geometric structures for the Co(n)+ clusters. The bond energies of Co(n)+-D for larger clusters (n > or = 10) are found to be very close to the value for chemisorption of atomic hydrogen on bulk-phase cobalt. The rate constants for D2 chemisorption on the cationic clusters are compared with the results from previous work on cationic and neutral cobalt clusters.     

Atomic contributions to bond dissociation energies in aliphatic hydrocarbons

This paper explores the atomic contributions to the electronic vibrationless bond dissociation enthalpy (BDE) at 0 K of the central C-C bond in straight-chain alkanes (C(n)H(2n+2)) and trans-alkenes (C(n)H(2n)) with an even number of carbon atoms, where n=2, 4, 6, 8. This is achieved using the partitioning of the total molecular energy according to the quantum theory of atoms in molecules by comparing the atomic energies in the intact molecule and its dissociation products. The study is conducted at the MP2(full)6-311++G(d,p) level of theory. It is found that the bulk of the electronic energy necessary to sever a single C-C bond is not supplied by these two carbon atoms (the alpha-carbons) but instead by the atoms directly bonded to them. Thus, the burden of the electronic part of the BDE is primarily carried by the two hydrogens attached to each of the alpha-carbons and by the beta-carbons. The effect drops off rapidly with distance along the hydrocarbon chain. The situation is more complex in the case of the double bond in alkenes, since here the burden is shared between the alpha-carbons as well as the atoms directly bonded to them, namely, again the alpha-hydrogens and the beta-carbons. These observations may lead to a better understanding of the bond dissociation process and should be taken into account when locally dense basis sets are introduced to improve the accuracy of BDE calculations.     

Atomic partitioning of the dissociation energy of the P-O(H) bond in hydrogen phosphate anion (HPO4(2-)): disentangling the effect of Mg2+

This paper has three goals: (1) to provide a first step in understanding the atomic basis of the role of magnesium in facilitating the dissociation of the P-O bond in phosphorylated biochemical fuel molecules (such as ATP or GTP), (2) to compare second-order M?ller-Plesset perturbation theory (MP2) results with those obtained at the more economical density functional theory (DFT) level for a future study of larger more realistic models of ATP/GTP, and (3) to examine the calculation of atomic total energies from atomic kinetic energies within a Kohn-Sham implemention of DFT, as compared to ab initio methods. A newly described method based on the quantum theory of atoms in molecules (QTAIM), which is termed the "atomic partitioning of the bond dissociation energy" (APBDE), is applied to a simple model of phosphorylated biological molecules (HPO42-). The APBDE approach is applied in the presence and in the absence of magnesium. It is found that the P-O(H) bond in the magnesium complex is shorter, exhibits a higher stretching frequency, and has a higher electron density at the bond critical point than in the magnesium-free hydrogen phosphate anion. Though these data would seem to suggest a stronger P-O(H) bond in the magnesium complex compared to the magnesium-free case, the homolytic breaking of the P-O(H) bond in the complex is found to be easier, i.e., has a lower BDE. This effect is the result of the balance of several atomic contributions to the BDE induced by the magnesium cation, which stabilizes the dissociation product more than it stabilizes the intact model molecule.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.     

Cooperative effects and strengths of hydrogen bonds in open-chain cis-triaziridine clusters (n = 2-8): a DFT investigation

We employ DFT/B3LYP method to investigate linear open-chain clusters (n = 2-8) of the cis-triaziridine molecule that is a candidate molecule for high energy density materials (HEDM). Our calculations indicate that the pervasive phenomena of cooperative effects are observed in the clusters of n = 3-8, which are reflected in changes in lengths of N...H hydrogen bonds, stretching frequencies, and intensities of N-H bonds, dipole moments, and charge transfers as cluster size increases. The n(N) --> sigma*(N-H) interactions, i.e., the charge transfers from lone pairs (n(N)) of the N atoms into antibonds (sigma*) of the N-H bonds acting as H-donors, can be used to explain the observed cooperative phenomena. The approaches based upon natural bond orbital (NBO) method and theory of atoms in molecule (AIM) to evaluating N...H strengths are found to be equivalent. In the process of N...H bonding, cooperative nature of n(N) --> sigma*(N-H) interactions promotes formation of stronger N...H bonds as reflected in increases in the capacities of cis-triaziridine clusters to concentrate electrons at the bond critical points of N...H bonds. The calculated nonadditive energies also show that the cooperative effects due to n(N) --> sigma*(N-H) interactions indeed provide additional stabilities for the clusters.     

Hydrogen bonding in phenol, water, and phenol-water clusters

Structure, stability, and hydrogen-bonding interaction in phenol, water, and phenol-water clusters have been investigated using ab initio and density functional theoretical (DFT) methods and using various topological features of electron density. Calculated interaction energies at MP2/6-31G level for clusters with similar hydrogen-bonding pattern reveal that intermolecular interaction in phenol clusters is slightly stronger than in water clusters. However, fusion of phenol and water clusters leads to stability that is akin to that of H(2)O clusters. The presence of hydrogen bond critical points (HBCP) and the values of rho(r(c)) and nabla(2)rho(r(c)) at the HBCPs provide an insight into the nature of closed shell interaction in hydrogen-bonded clusters. It is shown that the calculated values of total rho(r(c)) and nabla(2)rho(r(c)) of all the clusters vary linearly with the interaction energy.     

可快速计算水团簇三体作用强度的新方法

将水分子视为由2个O—H键偶极构成, 再将水分子间的三体作用视为长程诱导作用和短程校正之和, 使用Thole模型计算长程诱导作用, 通过同时考虑不同水分子间的置换和同一个水分子中2个键偶极间的置换计算短程校正, 从而提出了一个可快速计算水团簇三体作用强度的新方法. 根据已报道的12347个水三聚体的结构和CCSD(T)三体作用能, 确定了该方法所需参数. 将该方法和所确定的参数应用于67个水团簇体系, 计算这些体系的三体作用能, 并与CCSD(T), MP2, M06-2X方法的计算结果进行比较. 结果表明, 相对于CCSD(T)方法的总三体作用能, 本文方法的均方根偏差(RMSD)仅为3.32 kJ/mol, 平均相对偏差(MRD)仅为2.43%; 对较大水团簇体系, 该方法计算精度稍优于MP2方法, 明显优于M06-2X方法, 并且更快捷高效.     

A study of the reactions of molecular hydrogen with small gold clusters

This work presents a study of reactions between neutral and negatively charged Au(n) clusters (n=2,3) and molecular hydrogen. The binding energies of the first and second hydrogen molecule to the gold clusters were determined using density functional theory (DFT), second order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. It is found that molecular hydrogen easily binds to neutral Au(2) and Au(3) clusters with binding energies of 0.55 eV and 0.71 eV, respectively. The barriers to H(2) dissociation on these clusters with respect to Au(n)H(2) complexes are 1.10 eV and 0.59 eV for n=2 and 3. Although negatively charged Au(n) (-) clusters do not bind molecular hydrogen, H(2) dissociation can occur with energy barriers of 0.93 eV for Au(2) (-) and 1.39 eV for Au(3) (-). The energies of the Au(2)H(2) (-) and Au(3)H(2) (-) complexes with dissociated hydrogen molecules are lower than the energies of Au(2) (-)+H(2) and Au(3) (-)+H(2) by 0.49 eV and 0.96 eV, respectively. There is satisfactory agreement between the DFT and CCSD(T) results for binding energies, but the agreement is not as good for barrier heights.     

羟基二酰亚胺二聚体分子间相互作用的量子化学研究

采用从头算MP2方法在6-311++G~**基组水平上讨论了CP梯度校正对两种羟基 二酰亚胺异构体所的相互作用能和几何结构的影响，并利用分子中的原子理论（ Atoms in molecules，AIM）计算了五个拓扑参数：键临荷密度、电荷密度的 Iaplacian值、氢键中氢原子的体积、氢原子集居数、氢原子能量来表征氢键的形 成．种构型氢键体系中还讨论了二聚体的相互作用能与氢键临界点的电荷密度、质 子供体X-H键长的线性相关性问题．表明这种线性相关性的存在有范围限制，复合 物和其中单体的构型能够影响这种关系的存在.     

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen‐bonded amide–thymine and amide–uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen‐bonded amide–base systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded amide–base dimers quickly and accurately. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Microhydration of Alanine in Gas Phase Studied by Quantum Chemical Method and ABEEMσπ/MM Fluctuating Charge Model

A fluctuating charge interaction potential function for alanine-water was constructed in the spirit of newly developed ABEEMσπ/MM(atom-bond electronegativity equalization method at the σπ level fused into molecular mechanics). The properties of gaseous neutral alanine-(H₂O)_n(n=1-7) clusters were systematically investigated by quantum mechanics(QM) and the constructed ABEEMσπ/MM potential, such as conformations, hydrogen bonds (H-bonds), interaction energies, charge distributions, and so on. The results of ABEEMσπ/MM model are in fair agreement with those of QM and available experimental data. For isolated alanine, compared with those of experimental structure, the average absolute deviations(AAD) of bond length and bond angle are 0.002 nm and 1.4°, respectively. For alanine-water clusters, the AAD of interaction energies and H-bond lengths are only 3.77 kJ/mol and 0.012 nm, respectively, compared to the results of MP2/aug-cc-pVDZ//MP2/6-311+G^** method. The ABEEMσπ charges fluctuate with the changing conformation of the system, and can accurately and reasonably reflect the interpolarization between water and alanine. The presented alanine-water potential function may provide a basis for further simulations on related aqueous solutions of biomolecules.     

Origin of the size-dependence of the polarizability per atom in heterogeneous clusters: The case of AlP clusters

An analysis of the atomic polarizabilities α in stoichiometric aluminum phosphide clusters, computed at the MP2 and density functional theory (DFT) levels, the latter using the B3LYP functional, and partitioned using the classic and iterative versions of the Hirshfeld method, is presented. Two sets of clusters are examined: the ground-state Al(n)P(n) clusters (n=2-9) and the prolate clusters (Al(2)P(2))(N) and (Al(3)P(3))(N) (N≤6). In the ground-state clusters, the mean polarizability per atom, i.e., α/2n, decreases with the cluster size but shows peaks at n=5 and at n=7. We demonstrate that these peaks can be explained by a large polarizability of the Al atoms and by a low polarizability of the P atoms in Al(5)P(5) and Al(7)P(7) due to the presence of homopolar bonds in these clusters. We show indeed that the polarizability of an atom within an Al(n)P(n) cluster depends on the cluster size and the heteropolarity of the bonds it forms within the cluster, i.e., on the charges of the atoms. The polarizabilities of the fragments Al(2)P(2) and Al(3)P(3) in the prolate clusters were found to depend mainly on their location within the cluster. Finally, we show that the iterative Hirshfeld method is more suitable than the classic Hirshfeld method for describing the atomic polarizabilities and the atomic charges in clusters with heteropolar bonds, although both versions of the Hirshfeld method lead to similar conclusions.