Prism Isomers Research Articles

Evolution of Solute-Water Interactions in the Benzaldehyde-(H2O)1-6 Clusters by Rotational Spectroscopy.

The investigation on the preferred arrangement and intermolecular interactions of gas phase solute-water clusters gives insights into the intermolecular potentials that govern the structure and dynamics of the aqueous solutions. Here, we report the investigation of hydrated coordination networks of benzaldehyde-(water)n (n = 1-6) clusters in a pulsed supersonic expansion using broadband rotational spectroscopy. Benzaldehyde (PhCHO) is the simplest aromatic aldehyde that involves both hydrophilic (CHO) and hydrophobic (phenyl ring) functional groups, which can mimic molecules of biological significance. For the n = 1-3 clusters, the water molecules are connected around the hydrophilic CHO moiety of benzaldehyde through a strong CO···HO hydrogen bond and weak CH···OH hydrogen bond(s). For the larger clusters, the spectra are consistent with the structures in which the water clusters are coordinated on the surface of PhCHO with both the hydrophilic CHO and hydrophobic phenyl ring groups being involved in the bonding interactions. The presence of benzaldehyde does not strongly interfere with the cyclic water tetramer and pentamer, which retain the same structure as in the pure water cluster. The book isomer instead of cage or prism isomers of the water hexamer is incorporated into the microsolvated cluster. The PhCHO molecule deviates from the planar structure upon sequential addition of water molecules. The PhCHO-(H2O)1-6 clusters may serve as a simple model system in understanding the solute-water interactions of biologically relevant molecules in an aqueous environment.

The structures and properties of Mon (n = 2[formula omitted]15) cluster

The stable isomers of Mon (n=2∼15) are determined by using genetic algorithm combined with DFT calculations. Some new ground-state structures are found. The 3D structure firstly appears for Mo4 and the distorted pyramid or prism isomers are energy favored for n = 4∼9, while icosahedron-like or cage-like structures are dominant for n = 10∼15. Mo3, Mo5, Mo7, Mo8, Mo11 and Mo13 are magic number clusters. There are small magnetic moments and vertical electron affinity but large vertical ionization potential for Mon. The electronic structure analysis indicates that there form strong Mo-Mo covalent bonds, which are also delocalization to a certain extent due to the 4d and 5s orbital hybridization. The IR and UV–vis spectra are acquired and the particularly large hyper polarizability of Mo9 also suggests that it can be used as a nonlinear optical material.

A joint experimental and theoretical study on structural, electronic, and magnetic properties of MnGen - (n = 3-14) clusters.

A systematic structure and property investigation of MnGen - (n = 3-14) was conducted by means of density functional theory coupled with mass-selected anion photoelectron spectroscopy. This combined theoretical and experimental study allows global minimum and coexistence structures to be identified. It is found that the pentagonal bipyramid shape is the basic framework for the nascent growth process of MnGen - (n = 3-10), and from n = 10, the endohedral structures can be found. For n = 12, the anion MnGe12 - cluster probably includes two isomers: a major isomer with a puckered hexagonal prism geometry and a minor isomer with a distorted icosahedron geometry. Specifically, the puckered hexagonal prism isomer follows the Wade-Mingos rules and can be suggested as a new kind of superatom with the magnetic property. Furthermore, the results of adaptive natural density partitioning and deformation density analyses suggest a polar covalent interaction between Ge and Mn for endohedral clusters of MnGe12 -. The spin density and natural population analysis indicate that MnGen - clusters have high magnetic moments localized on Mn. The density of states diagram visually shows the significant spin polarization for endohedral structures and reveals the weak interaction between the Ge 4p orbital and the 4s, 3d orbitals of Mn.

Assessing the Performance of the Diffusion Monte Carlo Method As Applied to the Water Monomer, Dimer, and Hexamer.

The diffusion Monte Carlo (DMC) method is applied to the water monomer, dimer, and hexamer using q-TIP4P/F, one of the most simple empirical water models with flexible monomers. The bias in the time step (Δτ) and population size (Nw) is investigated. For the binding energies, the bias in Δτ cancels nearly completely, whereas a noticeable bias in Nw remains. However, for the isotope shift (e.g, in the dimer binding energies between (H2O)2 and (D2O)2), the systematic errors in Nw do cancel. Consequently, very accurate results for the latter (within ∼0.01 kcal/mol) are obtained with moderate numerical effort (Nw ∼ 10(3)). For the water hexamer and its (D2O)6 isotopomer, the DMC results as a function of Nw are examined for the cage and prism isomers. For a given isomer, the issue of the walker population leaking out of the corresponding basin of attraction is addressed by using appropriate geometric constraints. The population size bias for the hexamer is more severe, and to maintain accuracy similar to that of the dimer, Nw must be increased by ∼2 orders of magnitude. Fortunately, when the energy difference between the cage and prism is taken, the biases cancel, thereby reducing the systematic errors to within ∼0.01 kcal/mol when using a population of Nw = 4.8 × 10(5) walkers. Consequently, a very accurate result for the isotope shift is also obtained. Notably, both the quantum and isotope effects for the prism-cage energy difference are small.

Benchmark Structures and Harmonic Vibrational Frequencies Near the CCSD(T) Complete Basis Set Limit for Small Water Clusters: (H2O)n=2, 3, 4, 5, 6.

A series of (H2O)n clusters ranging from the dimer to the hexamer have been characterized with the CCSD(T) and the 2-body:Many-body CCSD(T):MP2 methods near the complete basis set (CBS) limit to generate benchmark-quality optimized structures and harmonic vibrational frequencies for these important systems. Quadruple-ζ correlation-consistent basis sets that augment the O atoms with diffuse functions have been employed in the analytic computation of harmonic vibrational frequencies for the global minima of the dimer, trimer, tetramer, and pentamer as well as the ring, book, cage, and prism isomers of the hexamer. Prior calibration [J. Chem. Phys. 2013, 139, 184113 and J. Chem. Theory Comput. 2014, 10, 5426] suggests that harmonic frequencies computed with this approach will lie within a few cm(-1) of the canonical CCSD(T) CBS limit. These data are used as reference values to gauge the performance of harmonic frequencies obtained with other ab initio methods (e.g., LCCSD(T) and MP2) and water potentials (e.g., TTM3-F and WHBB). This comparison reveals that it is far more challenging to converge harmonic vibrational frequencies for the bound OH stretching modes in these (H2O)n clusters to the CCSD(T) CBS limit than the free OH stretches, the n intramonomer HOH bending modes and even the 6n - 6 intermonomer modes. Deviations associated with the bound OH stretching harmonic frequencies increase rapidly with the size of the cluster for all methods and potentials examined, as do the corresponding frequency shifts relative to the monomer OH stretches.

Melting of the water hexamer

Results are presented for a molecular dynamics simulation of melting of the water hexamer from three-dimensional configurations (the book, cage, and prism isomers). The interactions between water molecules are described using two rigid potentials-TIP4P and TIP5P. The isomer geometry is determined using the connectivity matrix of the H-bond network graph. It is found in a numerical experiment that the Lindemann index, which shows the fluctuations in the length of intermolecular bonds, behaves similarly for all the studied initial isomers of the cluster. The melting of the cluster is found to depend on the interaction potential used in the study. For the cluster’s full energy, regions are identified which can be attributed to the premelting phase, which is dominated by the book isomer for TIP4P and ring isomer for TIP5P. The energy region attributed to the quasi-liquid phase of the cluster is shown to be characterized by isomers with ring and linear structures for both of the potentials.

Characteristics of hydrogen bond revealed from water clusters

The hydrogen bond network is responsible for the exceptional physical and chemical properties of water, however, the description of hydrogen bond remains a challenge for the studies of condensed water. The investigation of structural and binding properties of water clusters provides a key for understanding the H-bonds in bulk water. In this paper, a new set of geometric parameters are defined to describe the extent of the overlap between the bonding orbital of the donor OH and the nonbonding orbital of the lone-pair of the acceptor molecule. This orbital overlap plays a dominant role for the strength of H-bonds. The dependences of the binding energy of the water dimer on these parameters are studied. The results show that these parameters properly describe the H-bond strength. The ring, book, cage and prism isomers of water hexamer form 6, 7, 8 and 9 H-bonds, and the strength of the bonding in these isomers changes markedly. The internally-solvated and the all-surface structures of (H2O)n for n = 17, 19 and 21 are nearly isoenergetic. The internally-solvated isomers form fewer but stronger H-bonds. The hydrogen bonding in the above clusters are investigated in detail. The geometric parameters can well describe the characters of the H-bonds, and they correlate well with the H-bond strength. For the structures forming stronger H-bonds, the H-bond lengths are shorter, the angle parameters are closer to the optimum values, and their rms deviations are smaller. The H-bonds emanating from DDAA and DDA molecules as H-donor are relatively weak. The vibrational spectra of (H2O)n (n = 17, 19 and 21) are studied as well. The stretching vibration of the intramolecular OH bond is sensitive to its bonding environment. The H-bond strength judged from the geometric parameters is in good agreement with the bonding strength judged from the stretching frequencies.

Local-monomer calculations of the intramolecular IR spectra of the cage and prism isomers of HOD(D2O)5 and HOD and D2O ice Ih.

Dilute mixtures of HOD in pure H2O and D2O ices and liquid have been used by experimentalists to focus on the spectrum and vibrational dynamics of the local OH and OD stretches and bend of HOD in these complex and highly heterogeneous environments. The hexamer version of the mixture is HOD(D2O)5. The cage isomer of this cluster was recently studied and analyzed theoretically using local-mode calculations of the IR spectrum by Skinner and co-workers. This and the further possibility of experimental investigation of this cluster have stimulated us to study HOD(D2O)5 using the three-mode, local-monomer model, with the ab initio WHBB dipole moment and potential energy surfaces. Both the cage and prism isomers of this cluster are considered. In addition to providing additional insight into the HOD portion of the spectrum, the spectral signatures of each D2O are also presented in the range of 1000-4000 cm(-1). The OH stretch bands of both the prism and cage isotopomers exhibit rich structures in the range of 3100-3700 cm(-1) that are indicative of the position of the HOD in these hexamers. A preliminary investigation of the site preference of the HOD is also reported for both cage and prism HOD(D2O)5 using an harmonic zero-point energy analysis of the entire cluster. This indicates that the energies of free-OH sites are lower than the ones of H-bonded OH sites. Finally, following our earlier work on the IR spectra of H2O ice models, we present IR spectra of pure D2O and HOD.

Cage versus Prism: Electronic Energies of the Water Hexamer

Recent experiments show that the cage isomer of the water hexamer is lower in energy than the prism isomer near 0 K, and yet state-of-the-art electronic structure calculations predict the prism to be lower in energy than the cage at 0 K. Here, we study the relative energies of the water hexamers from the parametric two-electron reduced density matrix (2-RDM) method in which the 2-RDM rather than the wave function is the basic variable of the calculations. In agreement with experiment and in contrast with traditional wave function methods, the 2-RDM calculations predict the cage to be more stable than the prism after vibrational zero-point correction. Multiple configurations from the hydrogen bonding are captured by the method. More generally, the results are consistent with our previous 2-RDM applications in that they reveal how multireference correlation in molecular systems is important for resolving small energy differences from hydrogen bonding as well as other types of intermolecular forces, even in systems that are not conventionally considered strongly correlated.

The Water Hexamer: Cage, Prism, or Both. Full Dimensional Quantum Simulations Say Both

State-of-the-art quantum simulations on a full-dimensional ab initio potential energy surface are used to characterize the properties of the water hexamer. The relative populations of the different isomers are determined over a wide range of temperatures. While the prism isomer is identified as the global minimum-energy structure, the quantum simulations, which explicitly include zero-point energy and quantum thermal motion, predict that both the cage and prism isomers are present at low temperature down to almost 0 K. This is largely consistent with the available experimental data and, in particular, with very recent measurements of broadband rotational spectra of the water hexamer recorded in supersonic expansions.

Towards an ab initio flexible potential for water, and post-harmonic quantum vibrational analysis of water clusters

We describe a full-dimensional, flexible potential energy surface for arbitrary numbers of water monomers built from ab initio 2- and 3-body potentials. These potentials are each permutationally invariant fits to roughly 30 000 electronic energies. Tests of these potentials are made against direct high-level ab initio results for the water dimer and trimer. An application to the hexamer (prism isomer) is made along with a review of recent ab initio calculations of harmonic frequencies. One-mode quantum calculations of the OH-stretch fundamentals are done with these potentials and their accuracy is assessed by comparisons with coupled MULTIMODE calculations for the water dimer and trimer. Local-mode calculations of OH-stretch fundamentals are presented for the hexamer and decamer, where it appears that the bulk limit is being approached. A new local-monomer model is presented and tested for the water dimer and trimer.

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

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.

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