Ab initio investigation of water clusters (H2O)n (n = 2–34)

Various properties of typical structures of water clusters in the n = 2–34 size regime with the change of cluster size have been systematically explored. Full optimizations are carried out for the structures presented in this article at the Hartree–Fock (HF) level using the 6‐31G(d) basis set by taking into account the positions of all atoms within the cluster. The influence of the HF level on the results has been reflected by the comparison between the binding energies of (H₂O)_n (n = 2–6, 8, 11, 13, 20) calculated at the HF level and those obtained from high‐level ab initio calculations at the second‐order Møller–Plesset (MP2) perturbation theory and the coupled cluster method including singles and doubles with perturbative triples (CCSD(T)) levels. HF is inaccurate when compared with MP2 and CCSD(T), but it is more practical and allows us to study larger systems. The computed properties characterizing water clusters (H₂O)_n (n = 2–34) include optimal structures, structural parameters, binding energies, hydrogen bonds, charge distributions, dipole moments, and so on. When the cluster size increases, trends of the above various properties have been presented to provide important reference for understanding and describing the nature of the hydrogen bond. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

SCF,MP2, and CEPA-1 calculations on the OH ‥ O hydrogen bonded complexes (H2O)2 and (H2O-H2CO)

Counterpoise corrected ab initio calculations are reported for (H₂O)₂ and H₂O-H₂CO. Geometry searches were done in the moment-optimized basis DZP' at the SCF, MP2, and CEPA-1 levels of theory, followed by more accurate single-point calculations in basis ESPB, which includes bondfunctions to saturate the dispersion energy. The final equilibrium binding energies obtained are ?4.7 ±0.3 kcal/mol for a near-linear (H₂O)₂ structure and ?4.6 ±0.3 kcal/mol for a strongly bent HOH ‥ OCH₂ structure. The energy difference between these systems is much smaller than in all previous ab initio work. Cyclic (C_2h) and bifurcated (C_2v) transition structures for (H₂O)₂ are located at 1.0 ±0.1 kcal/mol and 1.9 ±0.3 kcal/mol above the global minimum, respectively. A new partitioning scheme is presented that rigorously partitions the MP2 correlation interaction energy in intra and intermolecular (dispersion) contributions. These terms are large (up to 2 kcal/mol) but of opposite sign for most geometries studied and hence their overall effect upon the final structures is relatively small. The relative merits of the MP2 and CEPA-1 approaches are discussed are discussed and it is concluded that for economical reasons MP2 is to be preferred, especially for larger systems.     

QMX: A versatile environment for hybrid calculations applied to the grafting of Al2Cl3Me3 on a silica surface

We present a new software to easily perform QM:MM and QM:QM' calculations called QMX. It follows the subtraction scheme and it is implemented in the Atomic Simulation Environment (ASE). Special attention is paid to couple molecular calculations with periodic boundaries approaches. QMX inherits the flexibility and versatility of the ASE package: any combination of methods namely force field, semiempirical, first principle, and ab initio, can be used as hybrid potential energy surface (PES). Its ease of use is demonstrated by considering the adsorption of Al₂Cl₃Me₃ on silica surface and by combining different levels of theory (from standard DFT to MP2 calculations) for the so‐called High Level cluster with standard PW91 density functional theory calculations for the Low Level environment. It is shown that the High Level cluster must contain the silanol group close to the aluminum atoms. The bridging adsorption is favored by 58 kJ mol^?1 at the MP2:PW91 level with respect to the terminal position. Using large clusters at the MP2:PW91 level, it is shown that PW91 calculations are sufficient for structure optimization but that embedded methods are required for accurate energy profiles. © 2013 Wiley Periodicals, Inc.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Dispersion interactions between neighboring Bi atoms in (BiH3)2 and Te(BiR2)2

Triggered by the observation of a short Bi⋯Bi distance and a Bi Te Bi bond angle of only 86.6° in the crystal structure of bis(diethylbismuthanyl)tellurane quantum chemical computations on interactions between neighboring Bi atoms in Te(BiR₂)₂ molecules (R = H, Me, Et) and in (BiH₃)₂ were undertaken. Bi⋯Bi distances atoms were found to significantly shorten upon inclusion of the d shells of the heavy metal atoms into the electron correlation treatment, and it was confirmed that interaction energies from spin component‐scaled second‐order Møller–Plesset theory (SCS‐MP2) agree well with coupled‐cluster singles and doubles theory including perturbative triples (CCSD(T)). Density functional theory‐based symmetry‐adapted perturbation theory (DFT‐SAPT) was used to study the anisotropy of the interplay of dispersion attraction and steric repulsion between the Bi atoms. Finally, geometries and relative stabilities of syn–syn and syn–anti conformers of Te(BiR₂)₂ (R = H, Me, Et) and interconversion barriers between them were computed. © 2018 Wiley Periodicals, Inc.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Modeling of the hydration sphere of gadolinium(III) ion using density functional theory

The hydrated gadolinium(III) ion cluster Gd³⁺ (H₂O)_x, with x = 8,9, was studied using density functional theory. The different conformations of the first hydration shell were calculated. For x = 8, the results for the cubic conformation correspond to previously published Hartree-Fock and MP2 results, whereas much lower energies were found for the square antiprismatic and dodecahedral conformations. For x = 9, the energy of the tricapped trigonal prism is nearly identical to the one of the antiprism plus an extra free water molecule. © 1997 John Wiley & Sons, Inc.     

Kinetic parameters for the reaction of hydroxyl radical with CH3OCH2F (HFE‐161) in the temperature range of 200–400 K: Transition state theory and Ab initio calculations

Rate coefficients for the reaction of the hydroxyl radical with CH₃OCH₂F (HFE‐161) were computed using transition state theory coupled with ab initio methods, viz., MP2, G3MP2, and G3B3 theories in the temperature range of 200–400 K. Structures of the reactants and transition states (TSs) were optimized at MP2(FULL) and B3LYP level of theories with 6‐31G* and 6‐311++G** basis sets. The potential energy surface was scanned at both the level of theories. Five different TSs were identified for each rotamer. Calculations of Intrinsic reaction coordinates were performed to confirm the existence of all the TSs. The kinetic parameters due to all different TSs are reported in this article. The rate coefficients for the title reaction were computed to be k = (9 ± 1.08) × 10^?13 exp [?(1,713 ± 33)/T] cm³ molecule^?1 s^?1 at MP2, k = (7.36 ± 0.42) × 10^?13 exp [?(198 ± 16)/T] cm³ molecule^?1 s^?1 at G3MP2 and k = (5.36 ± 1.57) × 10^?13 exp [?(412 ± 81)/T] cm³ molecule^?1 s^?1 at G3B3 theories. The atmospheric lifetimes of CH₃OCH₂F at MP2, G3MP2, and G3B3 level of theories were estimated to be 20, 0.1, and 0.3 years, respectively. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Theoretical study of the interaction d10‐d8 between Pt(0) and M(I) on the [Pt(PH3)3?MPH3] complexes (M = Cu,Ag, Au)

Ab initio calculations suggest that a series of complexes of type [Pt(PH₃)₃? MPH₃]⁺ (M = Au, Ag, Cu) are stable. We found that changes around the equilibrium distance Pt? M and in the interaction energies are sensitive to the electron correlation potential. This effect was evaluated using several levels of theory, including HF, MP2, and B3LYP. Both the magnitude of the interaction energies and distances Pt? M indicate a formal chemical bond, the latter being ratified by orbital diagram. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Reaction mechanism of NO with hydrolysates of NAMI-A: an MD simulation by combining the QM/MM(ABEEM) with the MD-FEP method

Nitrosylation reaction mechanisms of the hydrolysates of NAMI-A and hydrolysis reactions of ruthenium nitrosyl complexes were investigated in the triplet state and the singlet state. Activation free energies were calculated by combining the QM/MM(ABEEM) method with free energy perturbation theory, and the explicit solvent environment was simulated by an ABEEMσπ polarizable force field. Our results demonstrate that nitrosylation reactions of the hydrolysates of NAMI-A occur in both the triplet and the singlet states. The Ru-N-O angle of the triplet ruthenium nitrosyl complexes is in the range of 132.0°–138.2°. However, all the ruthenium nitrosyl complexes at the singlet state show an almost linear Ru-N-O angle. The nitrosylation reaction happens prior to the hydrolysis reaction for the first-step hydrolysates. The activation free energies of the nitrosylation reactions show that the H₂O-NO exchange reaction of [RuCl₄(Im)(H₂O)] in the singlet spin sate is the most likely one. Comparing with the activation free energies of the hydrolysis reactions of the ruthenium nitrosyl complexes, the results indicate that the rate of the DMSO–H₂O exchange reaction of [RuCl₃(NO)(Im)(DMSO)] is faster than that of [RuCl₃(H₂O)(Im)(DMSO)] in both the triplet spin state and the singlet spin state. © 2018 Wiley Periodicals, Inc.     

A computational study of the non-covalent bindings in complexes pairing sulfur tetroxide (SO4(C2V)) with the nitrous oxide (NNO)

The computational investigations are carried out on heterodimers containing sulfur tetroxide (SO₄(C_2V)) with the nitrous oxide (NNO) through MP2/cc–pVDZ and MP2/aug–cc–pVTZ//MP2/cc–pVDZ levels. Eight heterodimers are located on the potential energy surface of SO₄(C_2V)–NNO system. Binding energies of heterodimers in the SO₄(C_2V)–NNO system corrected with BSSE and ZPE are in the range of 1.17–7.90 kJ/mol. The calculated results reveal that the individual interaction of NNO terminal nitrogen atom with one of oxygen atoms of OSO ring in the SO₄(C_2V) monomer leads to the formation of the more stable heterodimer of SO₄(C_2V)–NNO system. The atoms in molecules theory were applied to analyze the nature of intermolecular interactions.     

The gas phase hydrogen‐bonded dimers of HOCl: A high‐level quantum chemical study

The O···H? O and Cl···H? O hydrogen bonding interactions were analyzed for HOCl dimers by using B3LYP, MP2, CCSD, and MP4(SDTQ) methods in conjunction with the various basis sets. Five isomers were found for the HOCl dimer. The ZPE and BSSE corrected binding energies were computed at the different levels of theory. At the optimized geometries obtained at CCSD/AUG‐cc‐pVDZ level, energies were re‐evaluated at MP4(SDTQ)/AUG‐cc‐pVTZ and CCSD(T)/cc‐pVTZ levels of theory. We found an average of ?20.9 and ?9.6 kJ/mol for the strength of the O···H and Cl···H hydrogen bonding interactions, respectively. Excitation and vertical ionization energies as well as rotational constants were computed at different levels of theory. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis were used to elucidate the nature of the interactions of HOCl dimers. The interaction energies were decomposed by Morokuma methodology. We have computed Δ_fH°(HOCl) and Δ_fH°(HOCl⁺) using the atomization reactions. The Δ_fH°₂₉₈(HOCl) values are ?17.85 and ?18.05 kcal/mol by using CBS‐Q and CBS‐QB3 extrapolation models, respectively, in good agreement with the results given in JANAF tables. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

Pt原子在γ-Al2O3(001)表面的吸附及迁移

采用密度泛函理论(DFT)中广义梯度近似(GGA)方法, 对Pt原子与γ-Al₂O₃(001)面的相互作用及迁移性能进行了研究. 分析了各种可能吸附位及吸附构型的松弛和变形现象, 吸附能和迁移能垒的计算结果表明: Pt团簇能够稳定吸附在该表面. Pt原子在表面O位的吸附能明显较高, 这主要是由Pt向基底O原子转移了电子所致. 电荷布居分析表明, Pt原子显电正性, Pt和Al原子之间存在排斥作用, 导致与Al原子产生较弱相互作用. 计算的平均吸附能大小依赖于Pt团簇的大小和形状, 总体趋势是随着Pt原子数增多, 吸附能降低. Pt原子在γ-Al₂O₃(001)表面迁移过程所需克服的迁移能垒最高值为0.51 eV. 随着吸附的Pt原子数增多,更倾向于形成Pt团簇. 因此, Pt原子在γ-Al₂O₃(001)表面的吸附演变不可能形成光滑、均匀平铺的吸附构型, 而在一定条件下容易出现团聚.     

Non-quantum-confinement optics in (CdS)n clusters

Geometry optimizations and the calculation of properties of (CdS)_n clusters for n = 2–12 are performed using density functional theory (DFT) with generalized gradient approximation. The all-electron basis is used for sulfur atoms while cadmium atoms are treated using relativistic effective core potentials. It is confirmed that the lowest energy state of (CdS)_n has a ring structure for n = 3–5 and a cage geometry for n ? 6. Optical absorption spectra are simulated using the energies and oscillator strengths of the first 20 singlet–singlet transitions computed for each cluster by time-dependent DFT. It is found that the energies of the lowest optical transitions do oscillate around the bulk gap value, which contradicts the concept of quantum confinement. The congender (CdSe)₁₂ and (CdTe)₁₂ clusters are found to show similar optical properties as the (CdS)₁₂ cluster.     

Photodissociation of O2+(H2O)

The photodissociation cross section of the weakly bound positive ion cluster O₂⁺(H₂O) has been measured at 15 discrete energies between 1.833 and 2.727 eV. Measurements indicate the cross section increases smoothly from 0.6 to 6 × 10^-18 cm² over this energy range. These cross section values are the largest reported for a positive ion cluster of atmospheric importance.     

Ab initio MO calculations on the stable conformations and their binding energies of the ion-molecule complexes: Ion = H+, Li+, Na+, K+ Be2+, and molecule = CO and N2

The most stable conformation of ion-molecule complexes involving a CO molecule were surveyed by the use of Hartree-Fock (HF) MO and third-order Moller-Plesset perturbation (MP3) methods with a 6–31G* basis set ion = H⁺, Li⁺, Na⁺, K⁺, Bc²⁺, Mg²⁺, and Ca²⁺. The MP3 level of theory reveals the ion-CO conformation in which the ion bonds to a carbon atom of CO to be the most stable; these MP3 results are contrary to the HF ones. Binding energies of ion-molecule complexes involving CO and N₂ were computed; MP3 energies are in good agreement with the experimental ones. The computed binding energies of cation-N₂ are about one-third of cation-NH₃ due to the absence of dipole moment and the smaller polarizability of N₂. The decrease in binding energy in cation-CO and -N₂ complexes, with increasing cation size, is mainly caused by the decrease of the electrostatic and polarization stabilizations.     

Theoretical study of MgO(001) surfaces: Pure,doped with Fe,Ca, and Al,and with and without adsorbed water

Ab initio calculations of large cluster models have been performed in order to study water adsorption at the five‐fold coordinated adsorption site on pure Mg(001) and MgO(001) surfaces doped with Fe, Ca, and Al. The geometric parameters of the adsorbed water molecule have been optimized preparatory to analysis of binding energies, charge transfer, preferential sites of interaction, and bonding distances. We have used Mulliken population analysis methods in order to analyze charge distributions and the direction of charge transfer. We have also investigated energy gaps, HOMO energies, and SCF orbital energies as well as the acid‐base properties of our cluster model. Numerical results are compared, where possible, with experiment and interpreted in the framework of various analytical models. © 2001 John Wiley & Sons, Inc. Int J Quant Chem, 2001