Theoretical studies of the water-cluster anions containing the OH{e}HO structure: energies and harmonic frequencies

In addition to isomers having a dipole-bound electron, the internally bound electron isomers of trimer, tetramer and hexamer water anions are found using an ab initio molecular orbital method. The latter isomers have a characteristic OH{e}HO structure. The interaction between the excess electron {e} and the surrounding OH bonds holds the structure stable. The calculated vibrational infrared spectrum for a hexamer anion with two double proton-acceptor water molecules shows a qualitatively similar vibrational spectrum with the one observed. A strong correlation between the vertical detachment energy and the distribution of the excess electron is also found.     

Intermolecular interaction in water hexamer

The origin of the intermolecular interaction, especially the many-body interaction, in eight low-lying water hexamer structures (prism, cage, book-1, book-2, cyclic-chair, bag, cyclic-boat-1, and cyclic-boat-2) is unraveled using the localized molecular orbital energy decomposition analysis (LMO-EDA) method at the second-order M?ller-Plesset perturbation (MP2) level of theory with a large basis set. It is found that the relative stabilities of these hexamer structures are determined by delicate balances between different types of interaction. According to LMO-EDA, electrostatic and exchange interactions are strictly pairwise additive. Dispersion interaction in these water hexamer structures is almost pairwise additive, with many-body effects varying from -0.13 to +0.05 kcal/mol. Repulsion interaction is roughly pairwise additive, with many-body effects varying from -0.84 to -0.62 kcal/mol. Polarization interaction is not pairwise additive, with many-body effects varying from -13.10 to -8.85 kcal/mol.     

Ab initio determination of the ionization potentials of water clusters (H2O)n (n = 2-6)

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic ionization potentials of several water clusters: dimer, trimer, tetramer, pentamer, hexamer book, hexamer ring, hexamer cage, and hexamer prism. The present results establish reference values at a level not reported before for these systems, calibrating different computational strategies and helping to discard less reliable theoretical and experimental data. The systematic study with the increasing size of the water cluster allows obtaining some clues on the structure and reductive properties of liquid water.     

Synthesis, characterization, and application of cationic water-soluble oligofluorenes in DNA-hybridization detection

A simple and efficient approach was developed for the synthesis of a series of cationic water-soluble oligofluorenes up to a chain length of a heptamer. Bromoalkyl-substituted fluorenyl boronic esters as the key intermediates were synthesized by using a modified Miyaura reaction. With an increasing number of repeat units (trimer to hexamer), the size-specific oligomers have shown redshifts in both the absorption and emission maxima. The emission maximum reaches the limit for the hexamer in both water and buffer solution. The quantum yields of the oligomers decreased with increased oligomer size in water. Both fluorescence quenching of the oligomers by 9,10-anthraquinone-2,6-disulfonate and the fluorescence resonance energy transfer experiments with the oligomers as the donor and fluorescein (Fl)-labeled double-stranded DNA (dsDNA-Fl) as the acceptor revealed the chain-length-dependent behavior. The Stern-Volmer quenching constant increased with the molecular size, whereas the highest donor-sensitized Fl emission was observed for the hexamer. These size-specific oligomers also served as a model to study the structure-property relationships for cationic polyfluorenes.     

Theoretical infrared spectrum of the ethanol hexamer

The hydrogen bond network of ethanol clusters is among the most complex hydrogen bond networks of molecular clusters. One of the reasons of its complexity arises from the number of possible ethanol monomers (there are three isoenergetic isomers of the ethanol monomer). This leads to difficulties in the exploration of potential energy surfaces (PESs) of ethanol clusters. In this work, we have explored the PES of the ethanol hexamer at the MP2/aug-cc-pVDZ level of theory. We have provided structures and their relative stability at 0 K and for temperatures ranging from 20 to 400 K in the gas phase. These structures are used to compute the theoretical infrared (IR) spectrum of the ethanol hexamer at the MP2/aug-cc-pVDZ level of theory. As a result, 98 different structures have been investigated, and six isomers are reported to be the most isoenergetically stable structures of the ethanol hexamer. These isomers are folded cyclic structures in which the stability is enhanced by the implication of CH⋯O interactions. Our investigations show that the PES of the ethanol hexamer is very flat, yielding several isoenergetic structures. Furthermore, we have noted that several isomers contribute to the population of the ethanol hexamer at high temperatures. As far as the IR spectroscopic study is concerned, we have found that the IR spectra of the most stable structures are in good agreement with the experiment. Considering this agreement, these structures are used to assign the experimental peaks in the CH-stretching region. We concluded that the stability of the structures of the ethanol hexamer is related both to OH⋯O hydrogen bonds and CH⋯O interactions. Overall, we have found that the IR spectrum of the ethanol hexamer, calculated from the contribution of all the possible stable structures weighted by their probability, excellently reproduce the experimental spectrum of the ethanol hexamer.     

Melting of the water hexamer

Journal of Structural Chemistry - Results are presented for a molecular dynamics simulation of melting of the water hexamer from three-dimensional configurations (the book, cage, and prism...     

Low-lying isomers and finite temperature behavior of (H2O)6 -

(H2O)(6) (-) appears as a "magic" number water cluster in (H2O)(n) (-) mass spectra. The structure of the (H2O)(6) (-) isomer dominating the experimental population has been established only recently [N. I. Hammer et al., J. Phys. Chem. A 109, 7896 (2005)], and the most noteworthy characteristic of this isomer is the localization of the excess electron in the vicinity of a double-acceptor monomer. In the present work, we use a quantum Drude model to characterize the low-energy isomers and the finite temperature properties of (H2O)(6) (-). Comparison with ab initio calculations shows that the use of a water model employing distributed polarizabilities and distributed repulsive sites is necessary to correctly reproduce the energy ordering of the low-lying isomers. Both the simulations and the ab initio calculations predict that there are several isomers of (H2O)(6) (-) significantly lower in energy than the experimentally observed species, suggesting that the experimental distribution is far from equilibrium.     

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. II. Spectroscopic signatures of the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms hydrogen bonding networks

We report the first harmonic vibrational spectra for each of the lowest lying isomers within the four major families of minima of (H2O)20, namely, the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms. These were obtained at the second-order Moller-Plesset perturbation level of theory (MP2) with the augmented correlation consistent basis set of double zeta quality (aug-cc-pVDZ) at the corresponding minimum energy geometries. The computed infrared (IR) spectra are the first ones obtained from first principles for these clusters. They were found to contain spectral features, which can be directly mapped onto the distinctive spectroscopic signatures of their constituent tetramer, pentamer, and octamer fragments. The dodecahedron spectra show the richest structure in the OH stretching region and are associated with the most redshifted OH vibrations with respect to the monomer. The lowest lying edge-sharing pentagonal prism isomer displays intense IR active vibrations that are redshifted by approximately 600 cm(-1) with respect to the water monomer. Furthermore the most redshifted, IR-active OH stretching vibrations for all four networks correspond to hydrogen bonded OH groups, which exhibit the following two common characteristics: (i) they belong to fragments which have a "free" OH stretch and (ii) they act as donors to a neighboring water molecule along a "dimerlike" (strong) hydrogen bond. The zero-point energy corrected MP2/CBS (complete basis set) limit binding energies D(0) for the four isomers are -163.1 kcal/mol (edge-sharing pentagonal prism), -160.1 kcal/mol (face-sharing pentagonal prism), -157.5 kcal/mol (fused cubes), and -148.1 kcal/mol (dodecahedron).     

Ring-Stacking Water Clusters: Morphology and Stabilities

The structures and interaction energies of water clusters with ring stacking motifs are studied by using ab initio calculations. The structures of the water clusters are constructed by stacking either single rings or multi-rings of tetramer, pentamer, and hexamer. We found that, in the single-ring-stacking motif, the most stable isomers exhibit an alternative clockwise-anticlockwise stacking pattern. We also show that four-layer single-ring-stacking isomers are not energetically favorable in comparison with those of two-layer multi-ring-stacking isomers. The relative stability of the isomers is also analyzed in terms of H-bond strength and elastic distortions of the water molecules.     

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%.     

Thermodynamics of forming water clusters at various temperatures and pressures by Gaussian-2, Gaussian-3, complete basis set-QB3, and complete basis set-APNO model chemistries; implications for atmospheric chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate Delta H degrees and Delta G degrees values for neutral clusters of water, (H(2)O)(n), where n = 2-6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 x 10(14) molecules/cm(3), water trimer concentrations of 2.6 x 10(12) molecules/cm(3), tetramer concentrations of approximately 5.8 x 10(11) molecules/cm(3), and pentamer concentrations of approximately 3.5 x 10(10) molecules/cm(3) in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.     

Structure and dynamics of binary and ternary lanthanide(III) and actinide(III) tris[4,4,4-trifluoro-1-(2-thienyl)-1,3-butanedione] (TTA)-tributylphosphate (TBP) complexes. Part 3, the structure, thermodynamics and reaction mechanisms of 8- and 9-coordinated binary and ternary Y-TTA-TBP complexes studied by quantum chemical methods

Possible mechanisms for intermolecular exchange between coordinated and solvent water in the complexes Y(TTA)(3)(OH(2))(2) and Y(TTA)(3)(TBP)(OH(2)) and intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(OH(2))(HTTA) and Y(TTA)(3)(TBP)(HTTA) have been investigated using ab initio quantum chemical methods. The calculations comprise both structures and energies of isomers, intermediates and transition states. Based on these data and experimental NMR data (Part 2) we have suggested intimate reaction mechanisms for water exchange, intramolecular exchange between structure isomers and intermolecular exchange between free HTTA and coordinated TTA. A large number of isomers are possible for the complexes investigated, but only some of them have been investigated, in all of them the most stable geometry is a more or less distorted square anti-prism or bicapped trigonal prism; the energy differences between the various isomers are in general small, less than 10 kJ mol(-1). 9-coordinated intermediates play an important role in all reactions. Y(TTA)(3)(OH(2))(3) has three non-equivalent water ligands that can participate in ligand exchange reactions. The fastest of these exchanging sites has a QM activation energy of 18.1 kJ mol(-1), in good agreement with the experimental activation enthalpy of 19.6 kJ mol(-1). The mechanism for the intramolecular exchange between structure isomers in Y(TTA)(3)(OH(2))(2) involves the opening of a TTA-ring as the rate determining step as suggested by the good agreement between the QM activation energy and the experimental activation enthalpy 47.8 and 58.3 J mol(-1), respectively. The mechanism for the intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(HTTA) and Y(TTA)(3)(TBP)(HTTA) involves the opening of the intramolecular hydrogen bond in coordinated HTTA followed by proton transfer to coordinated TTA. This mechanism is supported by the good agreement between experimental activation enthalpies (within parenthesis) and calculated activation energies 68.7 (71.8) and 35.3 (38.8) kJ mol(-1). The main reason for the difference between the two systems is the much lower energy required to open the intramolecular hydrogen bond in the latter. The accuracy of the QM methods and chemical models used is discussed.     

Accurate predictions of water cluster formation, (H₂O)(n=2-10)

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2-10 water molecules. This approach reveals new low energy conformers for (H(2)O)(n=7,9,10). The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.     

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

Exploration of the Ca2+ interaction modes of the nifedipine calcium channel antagonist.

A comprehensive study is carried out using quantum chemical computation and molecular dynamics (MD) simulations to gain insight into the interaction between Ca(2+) ions and the most important class of calcium channel antagonists--nifedipine. First, the chelating structures and energetic characters of nifedipine-Ca(2+) in the gas phase are explored, and 25 isomers are found. The most favorable chelating mode is a tridentate one, that is, Ca(2+) binds to two carbonyl O atoms and one nitryl O atom, where Ca(2+) is above the plane of the three O atoms to form a pyramidal structure. Accurate geometric structures, relative stabilities, vertical and adiabatic binding energies, and charge distributions are discussed. The differences in the geometries and energies among these isomers are analyzed from the contributions of chelating sites, electrostatics and polarizations, steric repulsions, and charge distributions. The interconversions among isomers with similar geometries and energies are also investigated because of the importance of the geometric transformation in the biological system. Furthermore, certain numbers of water molecules are added to the nifedipine-Ca(2+) system to probe the effect of water. A detailed study is performed on the hydrated geometries on the basis of the most stable isomer 1. Stepwise hydration can weaken the nifedipine-Ca(2+) interaction, and the chelating sites of nifedipine are gradually replaced by the added water molecules. Hexacoordination is found to be the most favorable geometry no matter how many water molecules were added, which can be verified by the MD simulations. The transfer of water molecules from the inner shell to the outer shell is also supported by MD simulations of the hexahydrated complexes.     

povel isomers of hexasulfur: Prediction of a stable prism isomer and implications for the thermal reactivity of elemental sulfur

High-level ab initio molecular orbital calculations were employed to explore the potential energy hypersurface of hexasulfur, S(6). Twelve isomeric structures of S(6) have been identified: two unbranched rings (chair and boat), one trigonal prism of D(3h) symmetry, two singly branched rings (S(5)double bondS), three triplet chains, one singlet chain, and three doubly branched rings (Sdouble bondS(4)double bondS). The prism structure is essentially a cluster of three S(2) molecules connected via a six-center pi(*)-pi(*)-pi(*) interaction. It is by 51 kJ mol(-1) less stable than the lowest-energy chair form. The reactions to generate the boat, the prism, and the singly branched isomers from the chair form are predicted to have lower barriers than the ring opening reaction of cyclo-S(6), which requires an activation energy of 149 kJ mol(-1). The prism and singly branched isomers are found to be more reactive species than the chair form and they are potential sources of S(2) in chemical reactions involving elemental sulfur.     

Quantum effects in ice Ih

Quantum and classical simulations are carried out on ice Ih over a range of temperatures utilizing the TIP4P water model. The rigid-body centroid molecular dynamics method employed allows for the investigation of equilibrium and dynamical properties of the quantum system. The impact of quantization on the local structure, as measured by the radial and spatial distribution functions, as well as the energy is presented. The effects of quantization on the lattice vibrations, associated with the molecular translations and librations, are also reported. Comparison of quantum and classical simulation results indicates that shifts in the average potential energy are equivalent to rising the temperature about 80 K and are therefore non-negligible. The energy shifts due to quantization and the quantum mechanical uncertainties observed in ice are smaller than the values previously reported for liquid water. Additionally, we carry out a comparative study of melting in our classical and quantum simulations and show that there are significant differences between classical and quantum ice.     

Structural isomers and reactivity for Rh6 and Rh6+

The structure, energetics, and interconversion of isomers of Rh(6) and Rh(6)(+) are studied by using density functional theory with Gaussian basis sets, using guess structures derived from basin-hopping simulations, and obtained by using the Sutton-Chen potential. A large range of spin multiplicities is considered for each isomer. Our calculations suggest two low-lying structures as possible structural isomers: a square bipyramid and a trigonal prism. The reactivity of these two candidate structural isomers with respect to adsorption of nitric oxide is studied via location of reaction transition states and calculation of reaction barriers. Similarities and differences with surface reaction studies are highlighted. These data provide powerful evidence that structural isomerism, and not different spin states, is responsible for the observed biexponential reaction kinetics.     

Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III

Several thermodynamic properties of ice Ih, II, and III are studied by a quasi-harmonic approximation and compared to results of quantum path integral and classical simulations. This approximation allows to obtain thermodynamic information at a fraction of the computational cost of standard simulation methods, and at the same time permits studying quantum effects related to zero-point vibrations of the atoms. Specifically, we have studied the crystal volume, bulk modulus, kinetic energy, enthalpy, and heat capacity of the three ice phases as a function of temperature and pressure. The flexible q-TIP4P/F model of water was employed for this study, although the results concerning the capability of the quasi-harmonic approximation are expected to be valid independently of the employed water model. The quasi-harmonic approximation reproduces with reasonable accuracy the results of quantum and classical simulations showing an improved agreement at low temperatures (T< 100 K). This agreement does not deteriorate as a function of pressure as long as it is not too close to the limit of mechanical stability of the ice phases.