Spectroscopic characterization of microscopic hydrogen-bonding disparities in supercritical water

The local hydrogen-bonding environment in supercritical water (380 degrees C, 300 bars, density 0.54 gcm3) was studied by x-ray Raman scattering at the oxygen K edge. The spectra are compared to those of the gas phase, liquid surface, bulk liquid, and bulk ice, as well as to calculated spectra. The experimental model systems are used to assign spectral features and to quantify specific local hydrogen-bonding situations in supercritical water. The first coordination shell of the molecules is characterized in more detail with the aid of the calculations. Our analysis suggests that approximately 65% of the molecules in supercritical water are hydrogen bonded in configurations that are distinctly different from those in liquid water and ice. In contrast to liquid water the bonded molecules in supercritical water have four intact hydrogen bonds and in contrast to ice large variations of bond angles and distances are observed. The remaining approximately 35% of the molecules exhibit two free O-H bonds and are thus either not involved in hydrogen bonding at all or have one or two hydrogen bonds on the oxygen side. We determine an average O-O distance of 3.1+/-0.1 A in supercritical water for the H bonded molecules at the conditions studied here. This and the corresponding hydrogen bond lengths are shown to agree with neutron- and x-ray-diffraction data at similar conditions. Our results on the local hydrogen-bonding environment with mainly two disparate hydrogen-bonding configurations are consistent with an extended structural model of supercritical water as a heterogeneous system with small patches of bonded molecules in various tetrahedral configurations and surrounding nonbonded gas-phase-like molecules.     

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.     

Cooperativity in ordinary ice and breaking of hydrogen bonds

The total interaction energy between two H-bonded water molecules in a condensed phase is composed of a binding energy between them and an energy due to a cooperative effect. An approximate simple expression is suggested for the dependence of the interaction energy between two H-bonded water molecules on the number of neighboring water molecules with which they are H-bonded. Using this expression, the probabilities of breaking a H bond with various numbers of H-bonded neighbors are estimated. These probabilities are used in computer simulations of the breaking of specified fractions of H bonds in an ordinary (hexagonal) ice. A large "piece" of hexagonal ice (up to 8 millions molecules) is built up, and various percentages of H bonds are considered broken. It is shown that 62-63% of H bonds must be broken in order to disintegrate the "piece" of ice into disconnected clusters. This value is only a little larger than the percolation threshold (61%) predicted both by the percolation theory for tetrahedral ice and by simulations in which all H bonds were considered equally probable to be broken. When the percentage of broken bonds is smaller than 62-63%, there is a network of H-bonded molecules which contains the overwhelming majority of water molecules. This result contradicts some models of water which consider that water consists of a mixture of water clusters of various sizes. The distribution of water molecules with unequal probabilities for breaking is compared with the simulation involving equal probabilities for breaking. It was found that in the former case, there is an enhanced number of water monomers without H bonds, that the numbers of 2- and 3-bonded molecules are smaller, and the number of 4-bonded molecules is larger than in the latter case.     

Hydrogen‐Bond Free Energy of Local Biological Water

Here, we propose an experimental methodology based on femtosecond‐resolved fluorescence spectroscopy to measure the hydrogen (H)‐bond free energy of water at protein surfaces under isothermal conditions. A demonstration was conducted by installing a non‐canonical isostere of tryptophan (7‐azatryptophan) at the surface of a coiled‐coil protein to exploit the photoinduced proton transfer of its chromophoric moiety, 7‐azaindole. The H‐bond free energy of this biological water was evaluated by comparing the rates of proton transfer, sensitive to the hydration environment, at the protein surface and in bulk water, and it was found to be higher than that of bulk water by 0.4 kcal mol^?1. The free‐energy difference is dominated by the entropic cost in the H‐bond network among water molecules at the hydrophilic and charged protein surface. Our study opens a door to accessing the energetics and dynamics of local biological water to give insight into its roles in protein structure and function.     

Interaction between NH2NO and H2O2: A quantum chemistry study

The H‐bonded complexes formed from interaction between NH₂NO (NA) and H₂O₂ (HP) have been investigated by using B3LYP and MP2 methods with a wide range of basis sets. We found six H‐bonded complexes in which three of them have cyclic structure. Calculations carried out at various levels show that the seven‐membered cyclic structure with O···HO and O···HN hydrogen bonding interactions is the most stable complex. The large binding energy obtained for A1 complex probably results from a more linear arrangement of the O···H N and O H···OH‐bonds in the seven‐membered structure A1. The natural bond orbital (NBO) analysis and the Bader's quantum theory of atoms in molecules have been used to elucidate the interaction characteristics of the NA‐HP complexes. The NBO results reveal that the charge transfer energy corresponds to the H‐bond interactions for A1 complex is grater than other complexes. The electrostatic nature of H‐bond interactions is predicted from QTAIM analysis. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Nature of one-dimensional short hydrogen bonding: bond distances, bond energies, and solvent effects

On the basis of recently synthesized calix[4]hydroquinone (CHQ) nanotubes which were self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (SHB), we have investigated the nature of 1-D SHB using first-principles calculations for all the systems including the solvent water. The H-bonds relay (i.e., contiguous H-bonds) effect in CHQs shortens the H...O bond distances significantly (by more than 0.2 A) and increases the bond dissociation energy to a large extent (by more than approximately 4 kcal/mol) due to the highly enhanced polarization effect along the H-bond relay chain. The H-bonds relay effect shows a large increase in the chemical shift associated with the SHB. The average binding energies for the infinite 1-D H-bond arrays of dioles and dions increase by approximately 4 and approximately 9 kcal/mol per H-bond, respectively. The solvent effect (due to nonbridging water molecules) has been studied by explicitly adding water molecules in the CHQ tube crystals. This effect is found to be small with slight weakening of the SHB strength; the H...O bond distance increases only by 0.02 A, and the average binding energy decreases by approximately 1 kcal/mol per H-bond. All these results based on the first-principles calculations are the first detailed analysis of energy gain by SHB and energy loss by solvent effect, based on a partitioning scheme of the interaction energy components. These reliable results elucidate not only the self-assembly phenomena based on the H-bond relay but also the solvent effect on the SHB strength.     

Spectroscopic Monitoring of the Acidity of Water Films on Ru(0001): Orientation‐Specific Acidity of Adsorbed Water

We examined the acid–base properties of water films adsorbed onto a Ru(0001) substrate by using surface spectroscopic methods in vacuum environments. Ammonia adsorption experiments combined with low‐energy sputtering (LES), reactive ion scattering (RIS), reflection–absorption infrared spectroscopy (RAIRS) and temperature‐programmed desorption (TPD) measurements showed that the adsorbed water is acidic enough to transfer protons to ammonia. Only the water molecules in an intact water monolayer and water clusters larger than the hexamer exhibit such acidity, whereas small clusters, a thick ice film or a partially dissociated water monolayer that contains OH, H₂O and H species are not acidic. The observations indicate the orientation‐specific acidity of adsorbed water. The acidity stems from water molecules with H‐down adsorption geometry present in the monolayer. However, the dissociation of water into H and OH on the surface does not promote but rather suppresses the proton transfer to ammonia.     

Surface Charges at the CaF2/Water Interface Allow Very Fast Intermolecular Vibrational‐Energy Transfer

We investigate the dynamics of water in contact with solid calcium fluoride, where at low pH, localized charges can develop upon fluorite dissolution. We use 2D surface‐specific vibrational spectroscopy to quantify the heterogeneity of the interfacial water (D₂O) molecules and provide information about the sub‐picosecond vibrational‐energy‐relaxation dynamics at the buried solid/liquid interface. We find that strongly H‐bonded OD groups, with a vibrational frequency below 2500 cm^?1, display very rapid spectral diffusion and vibrational relaxation; for weakly H‐bonded OD groups, above 2500 cm^?1, the dynamics slows down substantially. Atomistic simulations based on electronic‐structure theory reveal the molecular origin of energy transport through the local H‐bond network. We conclude that strongly oriented H‐bonded water molecules in the adsorbed layer, whose orientation is pinned by the localized charge defects, can exchange vibrational energy very rapidly due to the strong collective dipole, compensating for a partially missing solvation shell.     

Adducts of acylphloroglucinols with explicit water molecules: Similarities and differences across a sufficiently representative number of structures

Acylphloroglucinols constitute a broad class of compounds, derivatives of 1,3,5‐trihydroxybenzene, characterized by at least one COR group and exhibiting a variety of biological activities. The presence of several hydrogen bond donor or acceptor sites (the three phenol OH of the phloroglucinol moiety and the sp² O of the COR group), and their comparatively close spacing, makes the study of adducts with explicit water molecules particularly interesting, because it is possible to consider adducts in which water molecules surround the entire acylphloroglucinol molecule, or a large part of it, providing expectedly realistic images of possible arrangements of water molecules in the close vicinity of the acylphloroglucinol molecule in the aqueous medium. This work considers a number of different monomeric structures sufficiently representative of the broad structural variety of acylphloroglucinols and considers adducts of all the relevant conformers for each structure. Calculations use the HF/6‐31G(d,p) level because of affordability reasons in view of the adducts' size. The results: show that the intramolecular hydrogen bond (IHB) between the sp² O of COR and an ortho OH does not break on competition with solute–solvent intermolecular H‐bonding; highlight general trends and trends related to specific geometry features of the conformers; enable an interpretation of the additional solvent stabilization of the conformers without IHBs, observed from polarizable continuum model results in water solution; and highlight the significance, for this class of compounds, of considering adducts in whichthe water molecules directly H‐bonded to the central molecule are bridged by other water molecules, to approximate a continuous layer. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2378–2390, 2010     

X-ray absorption spectra of hexagonal ice and liquid water by all-electron Gaussian and augmented plane wave calculations

Full potential x-ray spectroscopy simulations of hexagonal ice and liquid water are performed by means of the newly implemented methodology based on the Gaussian augmented plane waves formalism. The computed spectra obtained within the supercell approach are compared to experimental data. The variations of the spectral distribution determined by the quality of the basis set, the size of the sample, and the choice of the core-hole potential are extensively discussed. The second part of this work is focused on the understanding of the connections between specific configurations of the hydrogen bond network and the corresponding contributions to the x-ray absorption spectrum in liquid water. Our results confirm that asymmetrically coordinated molecules, in particular, those donating only one or no hydrogen bond, are associated with well identified spectral signatures that differ significantly from the ice spectral profile. However, transient local structures, with half formed hydrogen bonds, may still give rise to spectra with dominant postedge contributions and relatively weaker oscillator strengths at lower energy. This explains why by averaging the spectra over all the O atoms of liquid instantaneous configurations extracted from ab initio molecular dynamics trajectories, the spectral features indicating the presence of weak or broken hydrogen bonds turn out to be attenuated and sometimes not clearly distinguishable.     

On the Debye-Waller factor of hexagonal ice: a computer simulation study

We investigate by molecular dynamics (MD) simulations the temperature dependence of the Debye-Waller (DW) factor of hexagonal ice with 25 different proton-disordered configurations. Each initial configuration is composed of 288 water molecules with no net dipole moment. The intermolecular interaction of water is described by TIP4P potential. Each production run of the simulation is 15 ns or longer. We observe a change in slope of the DW factor around 200 K, which cannot be explained within the framework of either classical or quantum harmonic approximation. Configurations generated by MD simulations are subjected to the steepest descent energy minimization. Analysis of the local energy minimum structures reveals that water molecules above 200 K jump to other lattice sites via some local energy minimum structures which contain some water molecules sitting on the locations other than the lattice sites. As time evolves, these defect molecules move back and forth to the lattice sites yielding defect-free structures. Those motions are responsible for the unusual increase in the DW factor at high temperatures. In making a transition from an energy-minimum structure to another one, a small number of water molecules are involved in a highly cooperative fashion. The larger DW factor at higher temperature arises from jump-like motions of water molecules among these locally stable configurations which may or may not be a family of the proton-disordered ice forms satisfying the "ice rule".     

Parallel Water/Aromatic Interactions of Non‐Coordinated and Coordinated Water

The parallel interactions of non‐coordinated and coordinated water molecules with an aromatic ring were studied by analyzing data in the Cambridge structural database (CSD) and by using quantum chemical calculations. The CSD data show that water/aromatic contacts prefer parallel to OH/π interactions, which indicates the importance of parallel interactions. The results reveal the influence of water coordination to a metal ion; the interactions of aqua complexes are stronger. Coordinated water molecules prefer a parallel‐down orientation in which one O?H bond is parallel to the aromatic ring, whereas the other O?H bond points to the plane of the ring. The interactions of aqua complexes with parallel‐down water/benzene orientation are as strong as the much better known OH/π orientations. The strongest calculated interaction energy is ?14.89 kcal mol^?1. The large number of parallel contacts in crystal structures and the quite strong interactions indicate the importance of parallel orientation in water/benzene interactions.     

DFT Study of Hydrogen-Bonded 1,3,5-Triazine-Water Complexes

The 1,3,5-triazine-water hydrogen bonding interactions have been investigated using the density functional theory B3LYP method and 6-31 ＋＋G^＊＊ basis, obtaining one, two and seven energy minima of the ground states for the 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively. The fully optimized geometries and binding energies were reported for the various stationary points. The global minima of 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes have a hydrogen bond N…H-O and a chain of water molecules, terminated by a hydrogen bond O…H-C. The binding energies are 13.38, 39.52 and 67.79 kJ/mol for the most stable 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively, after the basis set superposition error and zero point energy corrections. The H-O symmetric stretching modes of water in the complexes are red-shifted relative to those of the monomer water. In addition, the NBO analysis indicates that inter-molecule charge transfer is 0.02145 e, 0.02501 e and 0.02777 e for the most stable 1 ： 1, 1 ： 2 and 1 ： 3 complexes between 1,3,5-triazine and water, respectively.     

Use of molecular dynamics simulations with ab initio SCF calculations for the determination of the deuterium quadrupole coupling constant in liquid water and bond lengths in ice

The deuterium quadrupole coupling constant and asymmetry parameter in heavy water were determined using ab initio SCF calculations. Snapshots from a molecular dynamics simulation were used to give liquid water cluster configurations and the influence of simulation parameters on the quadrupole coupling constant was investigated. The electronic potential model and the number of molecules in the molecular dynamics simulation and the pressure of the system were found to have only a small influence on the quadrupole coupling constant. The average value of the quadrupole coupling constant at room temperature, corrected for the known deficiency of the ab initio calculation in the gas phase, yields a quadrupole coupling constant of 253 kHz, in perfect agreement with the most recent experiments. The oxygen—deuterium bond lengths in ice Ih, ice II, and ice IX were determined using experimental quadrupole coupling constants and a model equation. An averaged bond length of 98.9 pm was obtained for the Ih form, which is approximately 2 pm shorter than that determined by neutron diffraction studies, whereas the bond lengths for the four deuterium sites in ice II and the three sites in ice IX are in fair agreement with experiment. © John Wiley & Sons, Inc.     

Transfer of electron density as a result of hydrogen bond formation— A case of charged complexes

Previous investigation of transfer of electron density accompanying hydrogen bond formation has been extended to complexes between positively charged donors and neutral acceptors, as well as to the complexes between a neutral donor and a negatively charged acceptor molecules. The amount of transferred electron density from acceptor to donor for the charged complexes may be adequately described by the same exponential dependence on the equilibrium distance between the hydrogen atom and the nearest atom of the acceptor molecule as it was found for neutral complexes. Relation of the H‐bond energy to electron density at the H‐bond critical point was dependent on the sign of Laplacian of the electron density. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

The effect of boron nitride nanotubes size on the HArF interaction by NBO and AIM analysis

Three molecules HArF@BNNT(5,0), HArF@BNNT(6,0), and HArF@BNNT(7,0) have been formed by HArF encapsulated in boron nitride nanotubes (BNNTs) with different sizes. Due to the interaction between the HArF and the BNNTs, the H? Ar bond lengths are in a decrease trend, while the Ar? F bond lengths are in an increase trend compared with those of HArF. To investigate the nature of the interaction between H and Ar and the interaction between Ar and F, the quantum theory of “atoms in molecules” was carried out. The Laplacian (?²ρ_b) values of H? Ar suggest that the covalent interaction plays a key role in the H‐Ar interaction. For Ar? F, the results indicate that the Ar‐F interaction has a dominant noncovalent character. Moreover, the results obtained from the ratio of the kinetic‐energy density (G_b) and the potential‐energy density (V_b) (?G_b/V_b) and the total energy density (H_b) are in good agreement with that of ?²ρ_b values. In addition, the results of natural bond orbital charge and electron density difference between the HArF and BNNTs show that less electrons transfer away from the HArF to BNNTs with the gradual increase in the diameters of the BNNTs. © 2014 Wiley Periodicals, Inc.     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

Hydrogen atom formation from the photodissociation of water ice at 193 nm

The TOF spectra of photofragment hydrogen atoms from the 193 nm photodissociation of amorphous ice at 90-140 K have been measured. The spectra consist of both a fast and a slow components that are characterized by average translational energies of 2k(B)T(trans)=0.39+/-0.04 eV (2300+/-200 K) and 0.02 eV (120+/-20 K), respectively. The incident laser power dependency of the hydrogen atom production suggests one-photon process. The electronic excitation energy of a branched cluster, (H(2)O)(6+1), has been theoretically calculated, where (H(2)O)(6+1) is a (H(2)O)(6) cyclic cluster attached by a water molecule with the hydrogen bond. The photoabsorption of this branched cluster is expected to appear at around 200 nm. The source of the hydrogen atoms is attributed to the photodissociation of the ice surface that is attached by water molecules with the hydrogen bond. Atmospheric implications are estimated for the photodissociation of the ice particles (Noctilucent clouds) at 190-230 nm in the region between 80 and 85 km altitude.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Structures and electronic transport of water molecular nanotubes embedded in carbon nanotubes

In this paper, ice nanotubes confined in carbon nanotubes are investigated by molecular dynamics. The trigonal, square, pentagonal, and hexagonal water tubes are obtained, respectively. The current-voltage (I-V) curves of water nanotubes are found to be nonlinear, and fluctuations of conductance spectra of these ice nanotubes show that the transport properties of ice nanotubes are quite different from those of bulk materials. Our studies indicate that the conductance gap of ice nanotube is related to the difference value from the Fermi energy EF to the nearest molecular energy level E0. Increasing the diameter of a water molecular nanostructure results in the increase of the conductance.