The IR spectrum of supercritical water: Combined molecular dynamics/quantum mechanics strategy and force field for cluster sampling

Supercritical water was analyzed recently as a gas of small clusters of waters linked to each other by intermolecular hydrogen‐bonds, but unexpected “linear” conformations of clusters are required to reproduce the infra‐red (IR) spectra of the supercritical state. Aiming at a better understanding of clusters in supercritical water, this work presents a strategy combining classical molecular dynamics to explore the potential energy landscape of water clusters with quantum mechanical calculation of their IR spectra. For this purpose, we have developed an accurate and flexible force field of water based on the TIP5P 5‐site model. Water dimers and trimers obtained with this improved force field compare well with the quantum mechanically optimized clusters. Exploration by simulated annealing of the potential energy surface of the classical force field reveals a new trimer conformation whose IR response determined from quantum calculations could play a role in the IR spectra of supercritical water. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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     

The investigation on the relative stability of different clusters of thiourea dioxide in water using gas phase quantum chemical calculations

The relative stability of different clusters of thiourea dioxide (TDO) in water is examined using gas phase quantum chemical calculations at the MP2 and B3LYP level with 6‐311++G(d,p) basis set. The possible equilibrium structures and other energetic and geometrical data of the thiourea dioxide clusters, TDO‐(H₂O)_n (n is the number of water molecules), are obtained. The calculation results show that a strong interaction exists between thiourea dioxide and water molecules, as indicated by the binding energies of the TDO clusters progressively increased by adding water molecules. PCM model is used to investigate solvent effect of TDO. We obtained a negative hydration energy of ?20.6 kcal mol^?1 and free‐energy change of ?21.0 kcal mol^?1 in hydration process. On the basis of increasing binding energies with adding water molecules and a negative hydration energy by PCM calculation, we conclude thiourea dioxide can dissolve in water molecules. Furthermore, the increases of the C? S bond distance by the addition of water molecules show that the strength of the C? S bonds is attenuated. We find that when the number of water molecules was up to 5, the C? S bonds of the clusters, TDO‐(H₂O)₅ and TDO‐(H₂O)₆ were ruptured. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Molecular dynamics extended for fluctuating networks: Application to water

Molecular simulation models are increasingly important tools in efforts to understand the role that water plays in biochemical processes. However, existing models of water have limited capacity to deal with the characteristics of hydrogen bond networks. This article proposes a new fluctuating network (FN) algorithm as an extension of the standard molecular dynamics algorithm. The new algorithm allows for the simulation of a molecular system based on an underlying network, such as the hydrogen bond network in water. This algorithm distinguishes strong from weak network connections, applying a potential that best describes the specific connection behavior. We model liquid water with this new technique using a single‐site, isotropic, short‐range potential. We successfully reproduce liquid water's signature molecular spacing (as represented by the radial distribution function) and characterize its dynamic properties including the exponential hydrogen bond lifetime distribution, diffusion rate, and average hydrogen bonds per molecule. The FN algorithm allows exploration of the behavior of networked systems where explicit coordination limits are required. As such it could also be used to model covalent interactions, reaction dynamics, and applied to simulation of cellular networks. © 2012 Wiley Periodicals, Inc.     

BLUF Hydrogen network dynamics and UV/Vis spectra: A combined molecular dynamics and quantum chemical study

Blue light sensing using flavin (BLUF) protein photoreceptor domains change their hydrogen bond network after photoexcitation. To explore this phenomenon, BLUF domains from R. sphaeroides were simulated using Amber99 molecular dynamics (MD). Five starting configurations were considered, to study different BLUF proteins (AppA/BlrB), Trp conformations (“W_in”/“W_out”), structure determination (X‐ray/NMR), and finally, His protonation states. We found dependencies of the hydrogen bonds on almost all parameters. Our data show an especially strong correlation of the Trp position and hydrogen bonds involving Gln63. The latter is in some contradiction to earlier results (Obanayama et al., Photochem. Photobiol. 2008, 84 10031010). Possible origins and implications are discussed. Our calculations support conjectures that Gln63 is more flexible with Trp104 in W_in position. Using snapshots from MD and time‐dependent density functional theory, UV/vis spectra for the chromophore were determined, which account for molecular motion of the protein under ambient conditions. In accord with experiment, it is found that the UV/vis spectra of BLUF bound flavin are red‐shifted and thermally broadened for all calculated π → π* transitions, relative to gas phase flavin at T = 0 K. However, differences in the spectra between the various BLUF configurations cannot be resolved with the present approach. © 2012 Wiley Periodicals, Inc.     

An explanation of the unusual strength of the hydrogen bond in small water clusters

We seek to explain why the hydrogen bond possesses unusual strength in small water clusters that account for many of the complex behaviors of water. We have investigated and visualized the donation of covalent character from covalent (sigma) to hydrogen bonds by calculating the eigenvector coupling properties of quantum theory of atoms in molecules (QTAIM), stress tensor σ ( r ), and Ehrenfest Force F ( r ) on the F ( r ) molecular graph. The next-generation three-dimensional (3-D) bond-path framework sets are presented, and only the F ( r ) bond-path framework sets reproduce the earlier finding on the coupling between covalent (sigma) and hydrogen bonds that possess a degree of covalent character. Exploration of the bond-path between the covalent (sigma) and hydrogen bond's critical points provides an explanation for the previously obtained coupling results. The directional character of the covalent (sigma) and hydrogen bonds' 3-D bond-path framework sets for the F ( r ) explains differences found in the earlier results from QTAIM and the stress tensor σ ( r ).     

The water molecule arrangement over the side chain of some aliphatic amino acids: A quantum chemical and bottom-up investigation

The geometry of surrounding water molecules on the side chain of glycine, alanine, α-aminoisobutyric acid, α-aminobutyric acid, valine, and related hydrocarbons has been analyzed combining bottom-up and quantum chemical methodologies. To minimize the cavity size and to prevent water-water hydrogen bonding loss, the water molecules adopt a shape, resembling the one found in crystal structure of gas clathrate hydrates, with water molecules tangentially oriented to the surface of hydrophobic side chain. The cage is directly hydrogen bonded to the backbone's polar groups, thus hydration shells around hydrophobic and hydrophilic groups are folded together in amphiphilic molecules. The hydrophobe enclathration implies a substantial freedom degree reduction which makes it entropically disfavored. This disadvantageous entropic contribution is partially compensated by the favorable van der Waals interactions with guest in stabilizing clathrate hydrate formation. The water shell around the side chain relates intimately with the side-chain rotational isomerism. Present data are correlated with the experimental determined populations of the three rotamers, yielding promising results for both α-aminobutyric acid and valine.     

Theoretical study of hydrogen bonded clusters of water and cyanic acid: Hydrogen bonding in terms of the molecular structure

Ab initio and density functional calculations were used to analyze the interaction between a molecule of cyanic acid (HOCN) and up to 4 molecules of water at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) computational levels. The cooperative effect (CE) is increased with the increasing size of the studied clusters. Red shifts of the H–O stretching frequency for complexes involving HOCN as an H-donor were predicted. The strength of the hydrogen bonds in terms of molecular structures could be deduced from a comparison of HOCN–H₂O with HCNO–H₂O, HONC–H₂O and HNCO–H₂O HB clusters. The atom in molecules (AIM) method was used to analyze the cooperative effects on topological parameters.     

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.     

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.     

Dynamics and structural changes of small water clusters on ionization

Despite utmost importance in understanding water ionization process, reliable theoretical results of structural changes and molecular dynamics (MD) of water clusters on ionization have hardly been reported yet. Here, we investigate the water cations [(H₂O)_{n = 2–6}⁺] with density functional theory (DFT), Möller–Plesset second‐order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The complete basis set limits of interaction energies at the CCSD(T) level are reported, and the geometrical structures, electronic properties, and infrared spectra are investigated. The characteristics of structures and spectra of the water cluster cations reflect the formation of the hydronium cation moiety (H₃O⁺) and the hydroxyl radical. Although most density functionals fail to predict reasonable energetics of the water cations, some functionals are found to be reliable, in reasonable agreement with high‐level ab initio results. To understand the ionization process of water clusters, DFT‐ and MP2‐based Born‐Oppenheimer MD (BOMD) simulations are performed on ionization. On ionization, the water clusters tend to have an Eigen‐like form with the hydronium cation instead of a Zundel‐like form, based on reliable BOMD simulations. For the vertically ionized water hexamer, the relatively stable (H₂O)₅⁺ (5sL4A) cluster tends to form with a detached water molecule (H₂O). © 2013 Wiley Periodicals, Inc.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Energetic Ordering of Hydrogen Bond Strengths in Methanol-Water Clusters: Insights via Molecular Tailoring Approach

In this work, we examine the strength of various types of individual hydrogen bond (HB) in mixed methanol-water M_nW_m, (n+m=2 to 7) clusters, with an aim to understand the relative order of their strength, using our recently proposed molecular tailoring-based approach (MTA). Among all the types of HB, it is observed that the O_M−H…O_W HBs are the strongest (6.9 to 12.4 kcal mol⁻¹). The next ones are O_M−H…O_M HBs (6.5 to 11.6 kcal mol⁻¹). The O_W−H…O_W (0.2 to 10.9 kcal mol⁻¹) and O_W−H…O_M HBs (0.3 to 10.3 kcal mol⁻¹) are the weakest ones. This energetic ordering of HBs is seen to be different from the respective HB energies in the dimer i. e., O_M−H…O_M (5.0 to 6.0 kcal mol⁻¹)>O_W−H…O_M (1.5 to 6.0 kcal mol⁻¹)>O_M−H…O_W (3.8 to 5.6 kcal mol⁻¹)>O_W−H…O_W (1.2 to 5.0 kcal mol⁻¹). The plausible reason for the difference in the HB energy ordering may be attributed to the increase or decrease in HB strengths due to the formation of cooperative or anti-cooperative HB networks. For instance, the cooperativity contribution towards the different types of HB follows: O_M−H…O_W (2.4 to 8.6 kcal mol⁻¹)>O_M−H…O_M (1.3 to 6.3 kcal mol⁻¹)>O_W−H…O_W (−1.0 to 6.5 kcal mol⁻¹)>O_W−H…O_M (−1.2 to 5.3 kcal mol⁻¹). This ordering of cooperativity contribution is similar to the HB energy ordering obtained by the MTA-based method. It is emphasized here that, the interplay between the cooperative and anti-cooperative contributions are indispensable for the correct energetic ordering of these HBs.     

Ab initio study of the structure,cooperativity, and vibrational properties in the mixed hydrogen‐bonded trimers of hydrogen isocyanide and water

The five trimers of H₂O···HNC···H₂O, H₂O···H₂O···HNC, HNC···H₂O···H₂O, H₂O···HNC···HNC, and HNC···HNC···H₂O have been studied with quantum chemical calculations. Their structures, harmonic vibrational frequencies and interaction energies have been calculated at the B3LYP and MP2 levels with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets. The cooperative effect on these properties has also been studied quantitatively. For HNC:(H₂O)₂ systems, the cyclic H₂O···H₂O···HNC trimer is most stable with an interaction energy of ?16.01 kcal/mol and a large cooperative energy of ?3.25 kcal/mol at the MP2/aug‐cc‐pVTZ level. For H₂O:(HNC)₂ systems, the interaction energy and cooperative energy in the H₂O···HNC···HNC trimer are larger than those in the HNC···HNC···H₂O trimer. The NH stretch frequency has a blue shift for the terminal HNC molecule in the HNC···H₂O···H₂O and HNC···HNC···H₂O trimers and a red shift in other cases. A many‐body analysis has also been performed to understand the interaction energies in these hydrogen‐bonded clusters. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A scheme for rapid prediction of cooperativity in hydrogen bond chains of formamides,acetamides, and N‐methylformamides

A scheme is proposed in this article to predict the cooperativity in hydrogen bond chains of formamides, acetamides, and N‐methylformamides. The parameters needed in the scheme are derived from fitting to the hydrogen bonding energies of MP2/6‐31+G** with basis set superposition error (BSSE) correction of the hydrogen bond chains of formamides containing from two to eight monomeric units. The scheme is then used to calculate the individual hydrogen bonding energies in the chains of formamides containing 9 and 12 monomeric units, in the chains of acetamides containing from two to seven monomeric units, in the chains of N‐methylformamides containing from two to seven monomeric units. The calculation results show that the cooperativity predicted by the scheme proposed in this paper is in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction, demonstrating that the scheme proposed in this article is reasonable. Based on our scheme, a cooperativity effect of almost 240% of the dimer hydrogen bonding energy in long hydrogen bond formamide chains, a cooperativity effect of almost 190% of the dimer hydrogen bonding energy in long hydrogen bond acetamide chains, and a cooperativity effect of almost 210% of the dimer hydrogen bonding energy in long hydrogen bond N‐methylformamide chains are predicted. The scheme is further applied to some heterogeneous chains containing formamide, acetamide, and N‐methylformamide. The individual hydrogen bonding energies in these heterogeneous chains predicted by our scheme are also in good agreement with those obtained from Møller‐Plesset calculations including BSSE correction. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Insights into the dynamics of evaporation and proton migration in protonated water clusters from Large‐scale Born–Oppenheimer direct dynamics

Large‐scale on‐the‐fly Born–Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H₂O)_nH⁺, for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well‐known Eigen and Zundel motifs and is characterized by specific low‐frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed. © 2012 Wiley Periodicals, Inc.     

A molecular dynamics study of small water clusters comparing two flexible models for water

We report results from a molecular dynamics study of small water clusters, (H₂O) _n=2,3,4,6,8, comparing the recent polarizable, dissociable (PD) model of Halley et al. with the central force Stillinger-Rahman (SR) model. Dynamics calculations in the microcanonical ensemble quantified short-time averaged temperatures and RMS bond length fluctuations, which, in combination with simulated thermal quenching, identified structures and structural and phase changes of the clusters. Both models generally pointed to similar global minimum energy configurations, but local minima configurations, relative energies, and RMS bond length fluctuations varied more significantly. The PD model, with its incorporation of many-body effects, more accurately reproduced the structural features of clusters predicted from ab initio calculations, although it over-estimated the binding energies. Solid-like, pre-melting, and liquid-like states were identified based upon RMS bond length fluctuations for the O-O, O-H and H-H pairs. The SR and PD models yielded very similar trends in melting temperature as a function of cluster size. This suggests that general trends in the melting phenomenon do not depend strongly upon details of the models. The melting temperatures for the dimer and tetramer closely resembled the bulk melting temperature, while those for other sizes were considerably lower.     

Solvation of tetraalkylammonium chlorides in acetonitrile-water mixtures: mass spectrometry and molecular dynamics simulations.

The solvation of tetramethylammonium chloride (Me4NCl) and tetra-n-butylammonium chloride (Bu4NCl) in water-acetonitrile mixtures was investigated by mass spectrometry of clusters isolated from the solution. As far as the positive ions are concerned, clusters composed of alkylammonium ions and acetonitrile molecules only were observed, even for mixtures with high water content. In contrast, for the negative ions, clusters composed of chloride with both water and/or acetonitrile molecules were observed. For the smaller system (Me4NCl) we performed quantum chemical calculations and molecular dynamics simulations. It was found that even though water is present in the solvation shell of Me4N+, only acetonitrile has a strong electrostatic interaction with the cation. Water molecules around Me4N+ form hydrogen bonds with other water molecules, and they interact with Me4N+ mainly via dispersive interactions. These results indicate that Me4N+ behaves like a hydrophobic solute. On the other hand, the interaction of Cl- with water and acetonitrile is of comparable strength and, in both cases, the electrostatic interaction dominates. Herein we demonstrate experimentally and theoretically that positive and negative ions give rise to characteristic solvation structures in mixed solvents: even a relatively small organic cation, such as Me4N+, exhibits a hydrophobic-like solvation shell.     

Conformational analysis,tautomeric preference,intramolecular hydrogen bonding,and solvent effect on dinitrosamine: A quantum chemical study

HF, B3LYP, and MP2 methods with the standard basis set, 6‐311++G(d,p), were used to study various aspects of dinitrosamine. These results were compared with the outcomes of G2 and CBS‐QB3 methods. First, the conformational analysis and characterization of equilibrium conformations, especially global minima, were performed. On the basis of relative energies, we found that the dinitroso tautomers are more stable than the nitroso‐hydroxy (NH) ones. This preference is well‐interpreted in terms of tautomerization process and nitrosamine resonance. Furthermore, the nature of O? H···O intramolecular hydrogen bond (IMHB), in chelated forms of NH (NH‐11 and NH‐13) was comprehensively studied to evaluate the effect of hetero atoms (N) on the characteristic of IMHB systems. According to the results of isodesmic reaction method, the hydrogen bond energy of NH‐11 is greater than the malonaldehyde (MA) and NH‐13, whereas the electron density analysis and energy‐geometry correlation methods clearly predict that the hydrogen bond of NH‐11 is weaker than the MA. Additionally, the geometrical, atoms in molecules (AIM) and natural bond orbital's (NBO) parameters also emphasize on the MA as a chelated form with the strongest hydrogen bond. Finally, the solvent effects on the relative stability of selected dinitrosamine conformers are evaluated by different continuum (polarizable‐continuum model, isodensity polarizable continuum model, and self‐consistent isodensity polarizable continuum model), discrete and mixed solvent models. Theoretical results readily show that the potential energy surface of dinitrosamine, especially global minima, is strongly affected by the solvent. © 2012 Wiley Periodicals, Inc.