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     

D∞h B2(BS)2−/2− and Td B(BS)4−: Boron sulfide clusters containing BB multiple bonds and B− tetrahedral centers

A density functional theory investigation on the geometrical and electronic properties of B₄S (B₂(BS)) and B₅S (B(BS)) clusters has been performed in this work. Both the doublet B₂(BS) ([S?B? BB? B?S]^?) (D_∞h, ²Π_u) and the singlet B₂(BS) ([S?B? B?B? B?S]^2?) (D_∞h, ¹Σ) proved to have perfect linear ground‐state structures containing a multiply bonded BB core (BB or B?B) terminated with two BS groups, while T_d B(BS) turned out to possess a perfect B^? tetrahedral center directly corrected to four BS groups, similar to the corresponding boron hydride molecules of D_∞h B₂H, D_∞h B₂H, and T_d BH, respectively. B₄S₂ and B₅S₄ neutrals, however, appeared to be much different: they favor a planar fan‐shaped C_2v B₄S₂ (a di‐S‐bridged B₄ rhombus) and a planar kite‐like C_2v B₅S₄ (a di‐S‐bridged B₃ triangle bonded to two BS groups), respectively. One‐electron detachment energies and symmetrical stretching vibrational frequencies are calculated for D_∞h B₂(BS) and T_d B(BS) monoanions to facilitate their future characterizations. Neutral salts of B₂(BS)₂Li₂ with an elusive B?B triple bond and B(BS)₄Li containing a tetrahedral B^? center are predicted possible to be targeted in experiments. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Structures and basicity of LiNOH (N = 1 − 5) species

LiOH is one of the strong bases among neutral molecules. What about hydroxides of small Li_n (n = 2 ? 5) clusters? The addition of a single atom to a cluster sometimes has dramatic effects on its reactivity. This fact motivated us to perform an ab initio MP2/6‐311++G(d, p) investigation on Li_nOH species (n = 1 ? 5). These Li_nOH species are stabilized by both ionic as well as covalent interactions, and are found to be stable against elimination of LiOH and OH. We have determined their gas and aqueous phase basicity by considering hypothetical protonation reactions. The calculated proton affinities of Li_nOH (n ≥ 2) suggest their reduced basicity as compared to LiOH by 50–100 kJ/mol. The NBO charges and the highest occupied molecular orbitals also reveal the electride and alkalide characteristics of Li₂OH and Li₃OH, respectively. © 2016 Wiley Periodicals, Inc.     

Ab initio molecular dynamics studies of ionic dissolution and precipitation of sodium chloride and silver chloride in water clusters, NaCl(H2O)n and AgCl(H2O)n, n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.     

Structures,vibrational spectra,and relative energies of HXSiS (X = H,F, and Cl) isomers

The potential energy surface for the decomposition of HXSiS (X = H, F, and Cl) on the singlet state has been explored by B3LYP and CCSD(T) calculations. Five different types of reaction are proposed: (A) 1,1‐HX elimination, (B) 1,2‐H shift, (C) 1,2‐X shift, (D) H · and XSiS · radical formation, and (E) X · and HSiS · radical formation. These results show interesting trends for the HXSiS isomers. Our theoretical investigations suggest that the doubly bonded species HXSiS should be the lowest energy structure among the isomers from both kinetic and thermodynamic viewpoints. We also report theoretical predictions of molecular parameters and vibrational infrared (IR) spectra of the monohalogen‐substituted silanethione, which should be useful for future experimental observations. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 14–25, 2001     

Ab initio study of the reaction mechanism of water dissociation into the ionic species OH− and H3O+

A reaction mechanism of water dissociation is proposed where solvent effects are accounted for via a minimum stable model that considers the interaction of five water molecules. It is based on the fully self-consistent field (SCF) optimized structures of the reactant, product, and transition state, the calculations being at the Hartree–Fock and configuration interaction level [Møller–Plesset second-order perturbation (MP2) and coupled-cluster single and double excitations (CCSD)]. They were performed with four different basis sets that included polarized and diffuse orbitals. The dissociative mechanism leads to the ionic species OH⁻+H₃O⁺ as stable products and upon analysis of the energy hypersurface, a transition state is found which yields an activation barrier of 21.2 kcal/mol. This value is in good agreement with the experimentally determined enthalpy for the reaction. The contribution of the aggregation energy is emphasized. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 68: 253–259, 1998     

Electronic structure and stability of BP clusters: theoretical calculations for (BP)n (n = 2–4)

The geometry, electronic configurations, harmonic vibrational frequencies, and stability of the structural isomers of boron phosphide clusters have been investigated using density functional theory (DFT). CCSD(T) calculations show that the lowest‐energy structures are cyclic (IIt, IVs) with D_nh symmetry for dimers and trimers. The caged structure for B₄P₄ lie higher in energy than the monocyclic structure with D_2d symmetry (VIs). The B–P bond dominates the structures for many isomers, so that one preferred dissociation channel is loss of the BP monomer. The hybridization and chemical bonding in the different structures are also discussed. Comparisons with boron nitride clusters, the ground state structures of B_nP_n (n = 2, 3) clusters are analogous to those of their corresponding B_nN_n (n = 2, 3) counterparts. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

CH⋅⋅⋅N Hydrogen‐Bonding Interaction in 7‐Azaindole:CHX3 (X=F,Cl) Complexes

The C?H???Y (Y=hydrogen‐bond acceptor) interactions are somewhat unconventional in the context of hydrogen‐bonding interactions. Typical C?H stretching frequency shifts in the hydrogen‐bond donor C?H group are not only small, that is, of the order of a few tens of cm^?1, but also bidirectional, that is, they can be red or blue shifted depending on the hydrogen‐bond acceptor. In this work we examine the C?H???N interaction in complexes of 7‐azaindole with CHCl₃ and CHF₃ that are prepared in the gas phase through supersonic jet expansion using the fluorescence depletion by infra‐red (FDIR) method. Although the hydrogen‐bond acceptor, 7‐azaindole, has multiple sites of interaction, it is found that the C?H???N hydrogen‐bonding interaction prevails over the others. The electronic excitation spectra suggest that both complexes are more stabilized in the S₁ state than in the S₀ state. The C?H stretching frequency is found to be red shifted by 82 cm^?1 in the CHCl₃ complex, which is the largest redshift reported so far in gas‐phase investigations of 1:1 haloform complexes with various substrates. In the CHF₃ complex the observed C?H frequency is blue shifted by 4 cm^?1. This is at variance with the frequency shifts that are predicted using several computational methods; these predict at best a redshift of 8.5 cm^?1. This discrepancy is analogous to that reported for the pyridine‐CHF₃ complex [W. A. Herrebout, S. M. Melikova, S. N. Delanoye, K. S. Rutkowski, D. N. Shchepkin, B. J. van der Veken, J. Phys. Chem. A­ 2005 , 109, 3038], in which the blueshift is termed a pseudo blueshift and is shown to be due to the shifting of levels caused by Fermi resonance between the overtones of the C?H bending and stretching modes. The dissociation energies, (D₀), of the CHCl₃ and CHF₃ complexes are computed (MP2/aug‐cc‐pVDZ level) as 6.46 and 5.06 kcal mol^?1, respectively.     

Influence of the Li⋅⋅⋅π Interaction on the H/X⋅⋅⋅π Interactions in HOLi⋅⋅⋅C6H6⋅⋅⋅HOX/XOH (X=F,Cl, Br,I) Complexes

The influences of the Li???π interaction of C₆H₆???LiOH on the H???π interaction of C₆H₆???HOX (X=F, Cl, Br, I) and the X???π interaction of C₆H₆???XOH (X=Cl, Br, I) are investigated by means of full electronic second‐order Møller–Plesset perturbation theory calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The binding energies, binding distances, infrared vibrational frequencies, and electron densities at the bond critical points (BCPs) of the hydrogen bonds and halogen bonds prove that the addition of the Li???π interaction to benzene weakens the H???π and X???π interactions. The influences of the Li???π interaction on H???π interactions are greater than those on X???π interactions; the influences of the H???π interactions on the Li???π interaction are greater than X???π interactions on Li???π interaction. The greater the influence of Li???π interaction on H/X???π interactions, the greater the influences of H/X???π interactions on Li???π interaction. QTAIM studies show that the intermolecular interactions of C₆H₆???HOX and C₆H₆???XOH are mainly of the π type. The electron densities at the BCPs of hydrogen bonds and halogen bonds decrease on going from bimolecular complexes to termolecular complexes, and the π‐electron densities at the BCPs show the same pattern. Natural bond orbital analyses show that the Li???π interaction reduces electron transfer from C₆H₆ to HOX and XOH.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

Theoretical investigation of the structures and spectroscopic properties of (H2O4)n (n = 1–4) clusters

Density functional theory and ab initio calculations were performed to elucidate the hydrogen interactions in (H₂O₄)_n (n = 1–4) clusters. The optimized geometries, binding energies, and harmonic vibrational frequencies were predicted at various levels of theory. The trans conformer of the H₂O₄ monomer was predicted to be the most stable structure at the CCSD(T)/aug‐cc‐pVTZ level of theory. The binding energies per H₂O₄ monomer increased in absolute value by 9.0, 10.1, and 11.8 kcal/mol from n = 2 to n = 4 at the MP2/cc‐pVTZ level of theory (after the zero‐point vibrational energy and basis set superposition error corrections). This result implies that the intermolecular hydrogen bonds were stronger in the long‐chain clusters, that is, the formation of the longer chain in the (H₂O₄)_n clusters was more energetically favorable.     

Theoretical investigation of gas phase ethanol–(water)n (n = 1–5) clusters and comparison with gas phase pure water clusters (water)n (n = 2–6)

Various properties (such as optimal structures, structural parameters, hydrogen bonds, natural bond orbital charge distributions, binding energies, electron densities at hydrogen bond critical points, cooperative effects, and so on) of gas phase ethanol–(water)_n (n = 1–5) clusters with the change in the number of water molecules have been systematically explored at the MP2/aug‐cc‐pVTZ//MP2/6‐311++G(d,p) computational level. The study of optimal structures shows that the most stable ethanol‐water heterodimer is the one where exists one primary hydrogen bond (O? H…O) and one secondary hydrogen bond (C? H …O) simultaneously. The cyclic geometric pattern formed by the primary hydrogen bonds, where all the molecules are proton acceptor and proton donor simultaneously, is the most stable configuration for ethanol–(water)_n (n = 2–4) clusters, and a transition from two‐dimensional cyclic to three‐dimensional structures occurs at n = 5. At the same time, the cluster stability seems to correlate with the number of primary hydrogen bonds, because the secondary hydrogen bond was extremely weaker than the primary hydrogen bond. Furthermore, the comparison of cooperative effects between ethanol–water clusters and gas phase pure water clusters has been analyzed from two aspects. First of all, for the cyclic structure, the cooperative effect in the former is slightly stronger than that of the latter with the increasing of water molecules. Second, for the ethanol–(water)₅ and (water)₆ structure, the cooperative effect in the former is also correspondingly stronger than that of the latter except for the ethanol–(water)₅ book structure. © 2012 Wiley Periodicals, Inc.     

The arsenic clusters Asn (n = 1-5) and their anions: structures, thermochemistry, and electron affinities

The molecular structures, electron affinities, and dissociation energies of the As(n)/As(-) (n) (n = 1-5) species have been examined using six density functional theory (DFT) methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. These methods have been carefully calibrated (Chem Rev 2002, 102, 231) for the prediction of electron affinities. The geometries are fully optimized with each DFT method independently. Three different types of the neutral-anion energy separations reported in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The first dissociation energies D(e)(As(n-1)-As) for the neutral As(n) species, as well as those D(e)(As(-) (n-1)-As) and D(e) (As(n-1)-As(-)) for the anionic As(-) (n) species, have also been reported. The most reliable adiabatic electron affinities, obtained at the DZP++ BLYP level of theory, are 0.90 (As), 0.74 (As(2)), 1.30 (As(3)), 0.49 (As(4)), and 3.03 eV (As(5)), respectively. These EA(ad) values for As, As(2), and As(4) are in good agreement with experiment (average absolute error 0.09 eV), but that for As(3) is a bit smaller than the experimental value (1.45 +/- 0.03 eV). The first dissociation energies for the neutral arsenic clusters predicted by the B3LYP method are 3.93 eV (As(2)), 2.04 eV (As(3)), 3.88 eV (As(4)), and 1.49 eV (As(5)). Compared with the available experimental dissociation energies for the neutral clusters, the theoretical predictions are excellent. Two dissociation limits are possible for the arsenic cluster anions. The atomic arsenic results are 3.91 eV (As(-) (2) --> As(-) + As), 2.46 eV (As(-) (3) --> As(-) (2) + As), 3.14 eV (As(-) (4) --> As(-) (3) + As), and 4.01 eV (As(-) (5) --> As(-) (4) + As). For dissociation to neutral arsenic clusters, the predicted dissociation energies are 2.43 eV (As(-) (3) --> As(2) + As(-)), 3.53 eV (As(-) (4) --> As(3) + As(-)), and 3.67 eV (As(-) (5) --> As(4) + As(-)). For the vibrational frequencies of the As(n) series, the BP86 and B3LYP methods produce good results compared with the limited experiments, so the other predictions with these methods should be reliable.     

DFT and ab initio direct dynamics studies on the hydrogen abstraction reactions of chlorine atoms with CH(4-n)F(n) (n = 1-3)

A direct dynamics study is carried out for the hydrogen abstraction reactions Cl + CH(4-n)F(n) (n = 1-3) in the temperature range of 200-1,000 K. The minimum energy paths (MEPs) of these reactions are calculated at the BH&H-LYP/6-311G(d,p) level, and the energies along the MEPs are further refined at the QCISD(T)/6-311+G(2df,2p) and QCISD(T)/6-311+G(d,p) (single-point) level. The rate constants obtained by using the improved canonical variational transition state theory incorporating small-curvature tunneling correction (ICVT/SCT) are in good agreement with the available experimental results. It is shown that the vibrational adiabatic potential energy curves for these reactions have two barriers, a situation similar to the analogous reactions CH(3)X+Cl (X=Cl, Br). The theoretical results show that for the title reactions the variational effect should not be neglected over the whole considered temperature range, while the small-curvature tunneling effect is only important in the lower temperature range. The effects of fluorine substitution on the rate of this kind of reactions are also examined.     

Cooperativity between the Dihydrogen Bond and the N⋅⋅⋅HC Hydrogen Bond in LiH–(HCN)n Complexes

The cooperativity between the dihydrogen bond and the N???HC hydrogen bond in LiH–(HCN)_n (n=2 and 3) complexes is investigated at the MP2 level of theory. The bond lengths, dipole moments, and energies are analyzed. It is demonstrated that synergetic effects are present in the complexes. The cooperativity contribution of the dihydrogen bond is smaller than that of the N???HC hydrogen bond. The three‐body energy in systems involving different types of hydrogen bonds is larger than that in the same hydrogen‐bonded systems. NBO analyses indicate that orbital interaction, charge transfer, and bond polarization are mainly responsible for the cooperativity between the two types of hydrogen bonds.     

Carbon Dioxide Activation at Metal Centers: Evolution of Charge Transfer from Mg .+ to CO2 in [MgCO2(H2O)n].+, n=0–8

We investigate activation of carbon dioxide by singly charged hydrated magnesium cations Mg ^.+(H₂O)_n, through infrared multiple photon dissociation (IRMPD) spectroscopy combined with quantum chemical calculations. The spectra of [MgCO₂(H₂O)_n]^.+ in the 1250–4000 cm^?1 region show a sharp transition from n=2 to n=3 for the position of the CO₂ antisymmetric stretching mode. This is evidence for the activation of CO₂ via charge transfer from Mg ^.+ to CO₂ for n≥3, while smaller clusters feature linear CO₂ coordinated end‐on to the metal center. Starting with n=5, we see a further conformational change, with CO₂^.? coordination to Mg²⁺ gradually shifting from bidentate to monodentate, consistent with preferential hexa‐coordination of Mg²⁺. Our results reveal in detail how hydration promotes CO₂ activation by charge transfer at metal centers.     

Ab initio molecular orbital study on structures and energetics of CH3O−(H2O)n and CH3S−(H2O)n in gas phase

The purpose of this article was to calculate the structures and energetics of CH₃O⁻(H₂O)_n and CH₃S⁻(H₂O)_n in the gas phase; the maximum number of water molecules that can directly interact with the O of CH₃O⁻; and when n is larger, we asked how the CH₃O⁻ and CH₃S⁻ moiety of CH₃O⁻(H₂O)_n and CH₃S⁻(H₂O)_n changes and how we can reproduce experimental ΔH ⁰_{n−1, n}. Using the ab initio closed-shell self-consistent field method with the energy gradient technique, we carried out full geometry optimizations with the MP2/aug-cc-pVDZ for CH₃O⁻(H₂O)_n (n=0, 1, 2, 3) and the MP2/6–31+G(d,p) (for n=5, 6). The structures of CH₃S⁻(H₂O)_n (n=0, 1, 2, 3) were fully optimized using MP2/6–31++G(2d,2p). It is predicted that the CH₃O⁻(H₂O)₆ does not exist. We also performed vibrational analysis for all clusters [except CH₃O⁻(H₂O)₆] at the optimized structures to confirm that all vibrational frequencies are real. Those clusters have all real vibrational frequencies and correspond to equilibrium structures. The results show that the above maximum number of water molecules for CH₃O⁻ is five in the gas phase. For CH₃O⁻(H₂O)_n, when n becomes larger, the C—O bond length becomes longer, the C—H bond lengths become smaller, the HCO bond angles become smaller, the charge on the hydrogen of CH₃ becomes more positive, and these values of CH₃O⁻(H₂O)_n approach the corresponding values of CH₃OH with the n increment. The C—O bond length of CH₃O⁻(H₂O)₃ is longer than the C—O bond length of CH₃O⁻ in the gas phase by 0.044 Å at the MP2/aug-cc-pVDZ level of theory. The structure of the CH₃S⁻ moiety in CH₃S⁻(H₂O)_n does not change with the n increment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1138–1144, 1999     

六氢吡啶和氨形成的氢键团簇C5H10NH(NH3)n(n=1～3)结构的从头算研究

在RHF/6-31G(d)水平下,对C5H10NH(NH3)n(n=1～3)氢键团簇的平衡构型进行了从头算研究,优化得到各种可能的平衡构型.C5H10NH(NH3)为线型氢键结构,而C5H10NH(NH3)2为三元环结构,C5H10NH(NH3)3为四元环结构.在MP2/6-31G(d)//B3LYP/6-31G(d)水平下,对最稳定构型C5H10NH(NH3)n(Ⅰ)(n=1～3)的分子轨道进行布居分析,并且对相应的占据轨道进行指认.C5H10NH(NH3)n(Ⅰ)(n=1～3)垂直电离势的计算结果表明,形成氢键团簇后,分子的垂直电离势降低.