Hydrogen‐Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach

The cooperative effects of hydrogen bonding in small water clusters (H₂O)_n (n=3–6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ( r )) compliant with the quantum theory of atoms‐in‐molecules (QTAIM) and the calculation of electrostatic interactions by using one‐ and two‐electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (E_int) of these contorted molecules in (H₂O)_n. Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of E_int for a pair of water molecules in the cluster (H₂O)_n?1 become more attractive when a new H₂O unit is incorporated to generate the system (H₂O)_n with the last‐mentioned contribution being consistently the most important part of E_int throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω′), the QTAIM atomic charges, the topology of ρ( r ), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra‐ and intermolecular interactions present in hydrogen‐bonded systems.     

Study of water clusters in the n = 2–34 size regime, based on the ABEEM/MM model

Various properties of water clusters in the n = 2–34 size regime with the change of cluster size have been systemically explored based on the newly developed flexible-body and charge-fluctuating ABEEM/MM water potential model. The ABEEM/MM water model is to take ABEEM charges of all atoms, bonds, and lone-pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The computed correlating properties characterizing water clusters (H₂O) _n (n = 2–34) include optimal structures, structural parameters, ABEEM charge distributions, binding energies, hydrogen bonds, dipole moments, and so on. The study of optimal structures shows that the ABEEM/MM model can correctly predict the following important structural features, such as the transition from two-dimensional (from dimer to pentamer) to three-dimensional (for clusters larger than the hexamer) structures at hexamer region, the transition from cubes to cages at dodecamer (H₂O)₁₂, the transition from all-surface (all water molecules on the surface of the cluster) to one water-centered (one water molecule at the center of the cluster, fully solvated) structures at (H₂O)₁₇, the transition from one to two internal molecules in the cage at (H₂O)₃₃, and so on. The first three structural transitions are in good agreement with those obtained from previous work, while the fourth transition is different from that identified by Hartke. Subsequently, a systematic investigation of structural parameters, ABEEM charges, energetic properties, and dipole moments of water clusters with increasing cluster size can provide important reference for describing the objective trait of hydrogen bonds in water cluster system, and also provide a strong impetus toward understanding how the water clusters approach the bulk limit.     

Improved intermolecular force field for crystalline oxohydrocarbons including O?HO hydrogen bonding

An improved intermolecular force field of the (exp‐6‐1) type was obtained by fitting a training set of 124 observed oxohydrocarbon crystal structures and seven observed heats of sublimation. All of these structures, when energy minimized, showed cell edge length shifts of 3% or less. Especially good results were obtained for modeling carbohydrate crystal structures. The fitted crystal structures are a subset of a database of 180 structures systematically selected from the CSD library according to functional group and threshold accuracy criteria. Energy minimization results are also presented for 56 structures not in the training set, which showed cell edge length shifts larger than 3%. The previously published W99 force field, using C(4), C(3), and H(1) potentials, was slightly modified and adopted for hydrocarbon portions of the molecules. Oxygen atoms with one bond, O(1), and those with two bonds, O(2), were assigned separate parameters. Hydrogen atoms bonded to oxygen were assigned exponential repulsion functions and divided into two types: in hydroxy groups, H(2); and in carboxyl groups, H(3). Wavefunctions of HF 6‐31g** quality were calculated for every molecule and the molecular electric potential (MEP) was modeled with net atomic charges. Methylene bisector charges were used for all CH₂ and CH₃ groups, and ring center site charges were added if necessary to fit the MEP. For some structures, upward scaling of the MEP to simulate intermolecular polarization gave better results. MEP interaction generally gave a satisfactory representation of weak C HO hydrogen bonding. Medium strength C HO hydrogen bonding was modeled by reducing the repulsion of the involved hydrogen atoms. Crystal structures of several biologically interesting molecules were modeled with the force field, yielding good results. These molecules include aspirin, sucrose, β‐cellobiose, β‐lactose, progesterone, testosterone, prostaglandin E1, cholesterol acetate, and stearic acid. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 1–20, 2001     

Nonempirical Description of the Atmospherically Important Anionic Species. III. Hydrated Nitrogen Dioxide Anions

Neutral and negatively charged NO₂(H₂O) _n clusters are simulated at the second order of the Møller–Plesset perturbation theory with 6-31G basis set extended with diffuse and polarization functions on all nuclei. For better reliability, configuration interaction and multiconfiguration self-consistent field calculations with the active spaces, formed by all single and double excitations to the basic determinant, are carried out. The weak binding of a neutral NO₂ molecule to water clusters is provided by its coordination to two water molecules, either directly H bonded to each other or joined in an H-bond network via the third molecule. The presence of an excess electron strongly decreases the summary energy of the NO₂(H₂O) _n system, so that its adiabatic affinity exceeds the summary affinity of NO₂ and water system, although the excess electron is localized predominantly by NO₂ fragment.     

Atomic transferability within the exchange‐correlation density

Starting from either the exchange or the exchange‐correlation density together with Bader's definition of an atom in a molecule, an atomic hole density function can be defined. Contour maps of atomic hole density functions are able to show how the electron density of each atom in a molecule is partially delocalized into the rest of atoms in the molecule. The degree of delocalization of the atomic density ultimately depends on the nature of the atom studied and its environment. Atomic hole density functions are also used to define an atomic similarity measure, which allows for the quantitative assessment of the degree of atomic transferability in different molecular environments. In this article, contour maps for the N atom in the (N₂, CN⁻, NO⁺) series and the O atom in the (CO, H₂CO, and HCOOH) series are presented at the Hartree–Fock and CISD levels of theory. Moreover, the transferability of N and O within the two series is studied by means of atomic similarity measures. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1361–1374, 2000     

A multipole approximation of the electrostatic potential of molecules

The problem of approximating three-dimensional spatial distributions of quantum-mechanical electrostatic potentials of molecules by analytic potentials on the basis of atomic charges, real dipoles, and atomic multipoles up to quadrupoles inclusive was considered. Real dipole potentials are created by pairs of point charges of opposite signs, and the search for their arrangement in the volume of a molecule is part of the approximation problem. A FitMEP program was developed for the optimization of the parameters of models of the types specified taking into account molecular symmetry. It was shown for the example of several molecules (HF, CO, H₂O, NH₃, CH₄, formaldehyde, methanol, formamide, ethane, cyclopropane, cyclobutane, cyclohexane, tetrahedrane, cubane, adamantane, ethylene, and benzene) that the real dipole and atomic multipole models gave errors in approximated quantum-mechanical electrostatic potential values smaller by one or two orders of magnitude compared with the atomic charge model. The atomic charge model was shown to be virtually inoperative as applied to saturated hydrocarbons. Real dipole models were slightly inferior to atomic multipole models in quality but had all the advantages of the potential of point charges as concerned simplicity and compactness, and their use in potential energy calculations did not require changes in the existing program codes.     

The problem of small negatively ionized water clusters

Structures of the negatively charged water clusters (H₂O)⁻_n with n up to 6 were optimized at the UHF/4–31 + + G** level without any geometrical restriction. There was found a chainlike structure for the trimer anion and both chainlike and cyclic geometries for the larger clusters. According to the MP2 energy estimates, all the clusters with n ≥ 3 are more stable than the combinations of separated water molecules and a free electron. The energy of the cycles is lower than that of the chains of the same size. However, the latter species are energetically closer to the neutral oligomers of the same geometry. The energy of vertical detachment of an electron from the chainlike hexamer anion is already about zero. Addition of a diffuse s function centered equidistantly between the oxygen atoms proved metastability, in particular, of the cyclic (H₂O)⁻₄ structure. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 737–748, 1997     

Molecular orbital theory of the hydrogen bond.: XVI. A comparative study of pyridine and the diazines as proton acceptors in ground and excited n → π* states

Ab initio LCAO MO SCF calculations with a minimal STO-3G basis set have been performed to determine the structures and energies of dimers having pyridazine, pyrimidine, and pyrazine as proton acceptor molecules, with HF and H₂O as proton donors. The structures of these dimers are consistent with structures anticipated from the General Hybridization Model. Differences in the relative stabilities of dimers in the two series which have HF and H₂O as proton donors and pyridine and the diazines as proton acceptors are attributed to different weightings of secondary effects which influence dimer stabilities. These azabenzeme molecules form stronger hydrogen bonds than HCN and weaker hydrogen bonds than NH₃ whether HF or H₂O is the proton donor. Configuration interaction calculations indicate that vertical excitation to n → π* states of these proton aceptor molecules results in various degrees of destabilization of hydrogen bonded dimers and trimers, depending upon the excited state electron densities at the nitrogen atoms and the excited state dipole moments. With respect to the proton acceptor molecule, computed blue shifts of the n → π* bands increase in the order pyrazine < pyradizine < pyrimidine < pyridine.     

Nonempirical calculations of dipole moments of molecules in semifloating gaussian basis sets

For the example of the calculation of the dipole moments of the HF, HCl, H₂O, NH₃, CO, H₂CO, CH₃F molecules in two-exponent and three-exponent Gaussian basis sets, we have studied the effect of including floating functions in the basis, directly giving the effect of polarization of the electron shell of the atom in the molecule. We have established a weak dependence of the calculated dipole moment on the dimensionality of the basis, the number of floating functions, and also the orbital exponents of the hydrogen atoms. The correction introduced by the floating functions in molecules with polar bonds is considerably greater than the correlation correction. The proposed approach allows us to decrease the dimensionality of the orbital basis by a factor of 1.5–2 without making the agreement with experiment worse.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 26, No. 4, pp. 481–485, July–August, 1990.     

2,3,5,6‐Tetrakis[3,5‐bis(trifluoromethyl)phenoxy]‐2,5‐bis(dimethylamino)2,3,5,6‐tetrabora‐1,4‐dioxane diethyl ether 0.667‐solvate

The title compound, C₃₆H₂₆B₄F₂₄N₂O₆·0.667C₄H₁₀O, has centrosymmetric tetraboradioxane molecules, half each of three of these comprising the asymmetric unit together with a molecule of diethyl ether. Disorder affects most of the CF₃ groups and one ethyl group of the solvent molecule. The B₄O₂ rings are approximately planar and contain two B atoms with trigonal geometry and two with distorted tetrahedral geometry, the B—O bonds for the four‐coordinate B atoms being longer than those for the three‐coordinate B atoms. N—H...O hydrogen bonds link two of the crystallographically independent molecules together in chains, while the third molecule forms discrete trimolecular clusters with two solvent molecules via N—H...O hydrogen bonds. This is the first crystallographically characterized example of a tetrabora‐dioxane molecule containing both four‐ and three‐coordinate B atoms.     

Reactions in clusters of acetone and fluorinated acetones triggered by low energy electrons

Electron attachment to clusters of acetone (A), trifluoroacetone (TFA) and hexafluoroacetone (HFA) is studied in a crossed beam experiment with mass spectrometric detection of the anionic products. We find that the electron attachment properties in A change dramatically on going from isolated molecules to clusters. While single acetone is a very weak electron scavenger (via a dissociative electron attachment (DEA) resonance near 8.5 eV), clusters of A capture electrons at very low energy (close to 0 eV). The final ionic products consist of an ensemble of molecules (M) subjected to the loss of two neutral H₂ molecules ((M_n−2H₂)⁻, n ≥ 2). Their formation at low energies can only be explained by invoking new cyclic structures and polymers. In clusters of TFA, anionic complexes containing non-decomposed molecules (M_n⁻) including the monomer (M⁻) and ionic products formed by the loss of one and two HF molecules are observed. Loss of HF units is also interpreted by the formation of new cyclic structures in the anionic system. HFA is a comparatively stronger electron scavenger forming a non-decomposed anion via a narrow resonant feature near 0 eV in the gas phase. In HFA clusters, the non-decomposed parent anion is additionally observed at higher electron energies in the range 3–9 eV. The M⁻ signal carries signatures of self-scavenging processes, i.e., inelastic scattering by one molecule and capture of the completely slowed down electron by a second molecule within the same cluster. The scavenging spectrum is hence an image of the electronically excited states of the neutral molecule.     

A DFT study of H2O dissociation on metal‐precovered Fe (100) surface

The H₂O adsorption and dissociation on the Fe (100) surface with different precovered metals are studied by density functional theory. On both kinds of metal‐precovered surface, H₂O molecules prefer adsorb on hollow sites than bridge and top sites. The impurity energy difference is proportional to the adsorption energy, but the adsorbates are not sensitive to the adsorption orientation and height relative to the surface. The Hirshfeld charge analysis shows that water molecules act as an electron donor while the surface Fe atoms act as an electron acceptor. The rotation and dissociation of H₂O molecule occur on the Co‐ and Mn‐precovered surfaces. Some H₂O molecules are dissociated into OH and H groups. The energy barriers are about 0.5 to 1.0 eV, whose are consistence with the experimental data. H₂O molecules can be dissociated more easily at the top site on Co‐precovered surface 1 than that at bridge site on Mn‐precovered surface 2 because of the lower reaction barrier. The dispersion correction effects on the energies and adsorption configurations on Co‐precovered surface 1 were calculated by OBS + PW91. The dispersion contributions can improve a bit of the bond energy of adsorbates and weaken the hydrogen bond effect between adsorption molecules a little.     

Chemical bonding effect in electron scattering by gaseous molecules

The experimental intensity of 30 keV electron small angle scattering by a gaseous molecule is much different from the calculation using usual independent atom model. This is due to the rearrangement of electron distribution in a molecule by the formation of chemical bonds, and is called chemical bonding effect (CBE). The molecules studied are mainly hydrocarbons such as methane, acetylene, ethane, etc. and some non-hydrocarbons. The measurement was carried out on both elastic and total scattering and the effect was found for not only elastic but also inelastic scattering. The effect is relatively large for hydrogen rich molecules as H₂O, NH₃ and hydrocarbons, but is essentially related to the number of atoms contained in molecules. The origin of CBE will attribute mainly to the concentration of inner atomic electrons resulting from chemical bonding.     

Potential curves for proton motion in H-bonded systems: HF and H2O polymers

The presence of long range coupling between hydrogen atoms is shown for the HF and H₂O hydrogen bonded systems. The coupling of H atoms critically depends on the spatial orientation of the H atoms being considered. Explicit calculations of the potential curves of the protons are performed using as a model a ring of six HF, or H₂O, molecules. The method of calculation is the CNDO/2. The strong similarities of the results for H₂O and HF polymers supports the conclusion that the coupling is essentially due to factors such as the asymmetric equilibrium position of the H atoms, the high electronic polarizability of the system, etc.     

Off‐atomic site charges in some anions,metal cations,and complexes: Significance for electrostatic properties and binding

Electronic structures and properties of several anions, metal cations, and their complexes with neutral molecules were investigated at the HF/6‐31G** and B3LYP/6‐31G** levels of theory. Charges shifted from atomic sites due to atomic orbital hybridization called hybridization displacement charges (HDC) were investigated in detail. It has been found that many components of HDC are associated with each atom of ion that are shifted from the atomic sites, those associated with metal cations being shifted by large distances as found previously in electrically neutral systems. It is shown that atomic orbitals are appreciably rehybridized in going from neutral molecules to anions and cations. Molecular dipole moments and surface molecular electrostatic potentials (MEP) are obtained satisfactorily using HDC for the various types of species mentioned above. In the OH^?? H₂O complex, reversal of direction of shift of an HDC component associated with the hydrogen atom of H₂O involved in hydrogen bonding, indicates that the hydrogen bond between OH^? and H₂O would have some covalent character. Other atomic site‐based point charge models cannot provide such information about the nature of bonding. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem 2007     

Looking at orbitals in the laboratory: The experimental investigation of molecular wavefunctions and binding energies by electron momentum spectroscopy

The experimental technique of electron momentum spectroscopy (EMS ) (i.e., binary (e, 2e) spectroscopy) is discussed together with typical examples of its applications over the past decade in the area of experimental quantum chemistry. Results interpreted within the framework of the plane wave impulse and the target Hartree—Fock approximations provide direct measurements of, spherically averaged, orbital electron momentum distributions. Results for a variety of atoms and small molecules are compared with calculations using a range of Fourier transformed SCF position space wavefunctions of varying sophistication. Measured momentum distributions (MD ) provide a “direct” view of orbitals. In addition to offering a sensitive experimental diagnostic for semiempirical molecular wavefunctions, the MD's provide a chemically significant, additional experimental constraint to the usual variational optimization of wavefunctions. The measured MD's clearly reflect well known characteristics of various chemical and physical properties. It appears that EMS and momentum space chemistry offer the promise of supplementary perspectives and new vistas in quantum chemistry, as suggested by Coulson more than 40 years ago. Binding energy spectra in the inner valence region reveal, in many cases, a major breakdown of the simple MO model for ionization in accord with the predictions of many-body calculations. Results are considered for atomic targets, including H and the noble gases. The measured momentum distribution for H₂ is also compared with results from Compton scattering. Results for H₂ and H are combined to provide a direct experimental assessment of the bond density in H₂, which is compared with calculations. The behavior of the outer valence MD ''s for small row two and row three hydride molecules such as H₂O and H₂S, NH₃, HF, and HCl are consistent with well known differences in chemical and physical behavior such as ligand-donor activity and hydrogen bonding. MD measurements for the outermost valence orbitals of HF, H₂O and NH₃ show significant differences from those calculated using even very high-quality wavefunctions. Measurements of MD's for outer σ_g orbitals of small polyatomic molecules such as CO₂, COS, CS₂, and CF₄ show clear evidence of mixed s and p character. It is apparent that EMS is a sensitive probe of details of electronic structure and electron motion in atoms and molecules.     

Improved intermolecular force field for molecules containing H,C, N,and O atoms,with application to nucleoside and peptide crystals

A new intermolecular force field for nitrogen atoms in organic molecules was derived from a training dataset of 76 observed azahydrocarbon crystal structures and 11 observed heats of sublimation. The previously published W99 force field for hydrogen, carbon, and oxygen was thus extended to include nitrogen atoms. Nitrogen atoms were divided into four classes: N(1) for triply bonded nitrogen, N(2) for nitrogen with no bonded hydrogen (except the triple bonded case), N(3) for nitrogen with one bonded hydrogen, and N(4) for nitrogen with two or more bonded hydrogens. H(4) designated hydrogen bonded to nitrogen. Wavefunctions of 6‐31g** quality were calculated for each molecule and the molecular electric potential (MEP) was modeled with net atomic and supplementary site charges. Lone pair electron charge sites were included for nitrogen atoms where appropriate, and methylene bisector charges were used for CH₂ and CH₃ groups when fitting the MEP. X? H bond distances were set to standard values for the wave function calculation and then foreshortened by 0.1 Å for the MEP and force field fitting. Using the force field optimized to the training dataset, each azahydrocarbon crystal structure was relaxed by intermolecular energy minimization. Predicted maximum changes in unit cell edge lengths for each crystal were 3% or less. The complete force field for H, C, N, and O atoms was tested by intermolecular energy relaxation of nucleoside and peptide molecular crystals. Even though these molecules were not included in any of the training datasets for the force field, agreement with their observed crystal structures was very good, with predicted unit cell edge shifts usually less than 2%. These tests included crystal structures of representatives of all eight common nucleosides found in DNA and RNA, 15 dipeptides, four tripeptides, two tetrapeptides, and a pentapeptide with two molecules in the asymmetric unit. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1154–1166, 2001     

Structure of heteroassociates formed in the HF-(C<Subscript>2</Subscript>H<Subscript>5</Subscript>)<Subscript>2</Subscript>O liquid system

The aim of this work was to determine the structure of stable heteroassociates (HAs) with the stoichiometric ratios 1:2, 2:1, and 4:1 of molecules formed in the HF-(C₂H₅)₂O binary liquid system. The stretching frequencies of HF molecules found for each HA using a special procedure for processing IR spectra were compared with the calculated frequencies V HF of the stable molecular complexes (HF)_m ((C₂H₅)₂O)_n (m = 1, 2, 4, 8; n = 1, 2) with different topologies by the density functional method (B3LYP/6-31++G(d,p)). As a result, it was shown that the most stable (among H-bonded complexes with the same stoichiometric ratio of molecules) HAs HF((C₂H₅)₂O)₂, (HF)₄ ((C₂H₅)₂O)₂, and (HF)₈-((C₂H₅)₂O)₂ formed in HF solutions in diethyl ether. All of them had a cyclic structure and a common peculiarity of structure: only one lone electron pair of the oxygen atom of the (C₂H₅)₂O molecules is involved in hydrogen bonding.     

The structure of {[Co(pht)(bpy)(H2O)]·3H2O}n (pht is phthalate and bpy is 4,4′‐bipyridine) and the role of solvent water clusters in structure stability

The Co^II cation in poly[[aqua(μ‐benzene‐1,2‐dicarboxylato‐κ³O¹,O²:O¹)(μ‐4,4′‐bipyridine‐κ²N:N′)cobalt(II)] trihydrate], {[Co(C₈H₄O₄)(C₁₀H₈N₂)(H₂O)]·3H₂O}_n, is octahedrally coordinated by two N atoms of two 4,4′‐bipyridine ligands, three O atoms from phthalate anions and a fourth O atom from a coordinated water molecule. The packing consists of planes of coordination polymers linked by hydrogen bonds mediated by three solvent water molecules; the linkage is achieved by the water molecules forming intricate oligomeric clusters which also involve the O atoms of the phthalate ligands.     

Computer-aided study of oxygen absorption by water clusters. IR spectra of heteroclusters

Spectral characteristics of (H₂O)_n, (O₂)_m(H₂O)_n, and (O)_i(H₂O)_n cluster systems, where m≤2, i≤4, and 10 ≤ n ≤ 50, are studied with the molecular dynamics method using a flexible molecule model. The IR absorption spectra are changed substantially as a result of O₂ molecule dissociation, and in the presence of atomic oxygen in the clusters, the spectra are characterized by a deep minimum at 520 cm^?1. The absorption of oxygen causes a marked reduction in reflection coefficient R of monochromatic IR radiation. The number of peaks in the R(ω) spectra decreases to two in the case of molecular oxygen absorption and is no larger than four in the case of atomic oxygen absorption. The absorption of atomic oxygen by the clusters is also accompanied by a significant increase in the dissipation of energy accumulated by the clusters. This effect weakens when molecular oxygen is absorbed. An increase in atomic oxygen concentration in the clusters renders their radiation harder.