QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Structural and thermodynamic properties of different phases of supercooled liquid water.

Computer simulation results are reported for a realistic polarizable potential model of water in the supercooled region. Three states, corresponding to the low density amorphous ice, high density amorphous ice, and very high density amorphous ice phases are chosen for the analyses. These states are located close to the liquid-liquid coexistence lines already shown to exist for the considered model. Thermodynamic and structural quantities are calculated, in order to characterize the properties of the three phases. The results point out the increasing relevance of the interstitial neighbors, which clearly appear in going from the low to the very high density amorphous phases. The interstitial neighbors are found to be, at the same time, also distant neighbors along the hydrogen bonded network of the molecules. The role of these interstitial neighbors has been discussed in connection with the interpretation of recent neutron scattering measurements. The structural properties of the systems are characterized by looking at the angular distribution of neighboring molecules, volume and face area distribution of the Voronoi polyhedra, and order parameters. The cumulative analysis of all the corresponding results confirms the assumption that a close similarity between the structural arrangement of molecules in the three explored amorphous phases and that of the ice polymorphs I(h), III, and VI exists.

Communication: energy benchmarking with quantum Monte Carlo for water nano-droplets and bulk liquid water.

We show the feasibility of using quantum Monte Carlo (QMC) to compute benchmark energies for configuration samples of thermal-equilibrium water clusters and the bulk liquid containing up to 64 molecules. Evidence that the accuracy of these benchmarks approaches that of basis-set converged coupled-cluster calculations is noted. We illustrate the usefulness of the benchmarks by using them to analyze the errors of the popular BLYP approximation of density functional theory (DFT). The results indicate the possibility of using QMC as a routine tool for analyzing DFT errors for non-covalent bonding in many types of condensed-phase molecular system.

Structure of coexisting liquid phases of supercooled water: analogy with ice polymorphs.

The structural changes occurring in supercooled liquid water upon moving from one coexisting liquid phase to the other have been investigated by computer simulation using a polarizable interaction potential model. The obtained results favorably compare with recent neutron scattering data of high and low density water. In order to assess the physical origin of the observed structural changes, computer simulation of several ice polymorphs has also been carried out. Our results show that there is a strict analogy between the structure of various disordered (supercooled) and ordered (ice) phases of water, suggesting that the occurrence of several different phases of supercooled water is rooted in the same physical origin that is responsible for ice polymorphism.

Movement of deposited water on turbomachinery rotor blade surfaces

A theoretical approach for calculating the movement of liquid water following deposition onto a turbomachine rotor blade is described. Such a situation can occur during operation of an aero-engine in rain. The equation of motion of the deposited water is developed on an arbitrarily oriented plane triangular surface facet. By dividing the blade surface into a large number of facets and calculating the water trajectory over each one crossed in turn, the overall trajectory can be constructed. Apart from the centrifugal and Coriolis inertia effects, the forces acting on the water arise from the blade surface friction, and the aerodynamic shear and pressure gradient. Non- dimensionalisation of the equations of motion provides considerable insight and a detailed study of water flow on a flat rotating plate set at different stagger angles demonstrates the paramount importance of blade surface friction. The extreme cases of low and high blade friction are examined and it is concluded that the latter (which allows considerable mathematical generalisation) is the most likely in practice. It is also shown that the aerodynamic shear force, but not the pressure force, may influence the water motion. Calculations of water movement on a low-speed compressor blade and the fan blade of a high bypass ratio aero-engine suggest that in low rotational speed situations most of the deposited water is centrifuged rapidly to the blade tip region. Copyright © 2006 by ASME.

A comparison of the value of viscosity for several water models using Poiseuille flow in a nano-channel.

The viscosity-temperature relation is determined for the water models SPC/E, TIP4P, TIP4P/Ew, and TIP4P/2005 by considering Poiseuille flow inside a nano-channel using molecular dynamics. The viscosity is determined by fitting the resulting velocity profile (away from the walls) to the continuum solution for a Newtonian fluid and then compared to experimental values. The results show that the TIP4P/2005 model gives the best prediction of the viscosity for the complete range of temperatures for liquid water, and thus it is the preferred water model of these considered here for simulations where the magnitude of viscosity is crucial. On the other hand, with the TIP4P model, the viscosity is severely underpredicted, and overall the model performed worst, whereas the SPC/E and TIP4P/Ew models perform moderately.

Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng-Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes. © 2013.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding.

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.

Machine learning for predictive condensed-phase simulation

We show how machine learning techniques based on Bayesian inference can be used to reach new levels of realism in the computer simulation of molecular materials, focusing here on water. We train our machine-learning algorithm using accurate, correlated quantum chemistry, and predict energies and forces in molecular aggregates ranging from clusters to solid and liquid phases. The widely used electronic-structure methods based on density-functional theory (DFT) give poor accuracy for molecular materials like water, and we show how our techniques can be used to generate systematically improvable corrections to DFT. The resulting corrected DFT scheme gives remarkably accurate predictions for the relative energies of small water clusters and of different ice structures, and greatly improves the description of the structure and dynamics of liquid water.

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number W e = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20 % of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number We = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20% of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

Measurement of liquid film thickness by optical fluorescence and its application to an oscillating piston positive displacement flowmeter

The movement of the circular piston in an oscillating piston positive displacement flowmeter is important in understanding the operation of the flowmeter, and the leakage of liquid past the piston plays a key role in the performance of the meter. The clearances between the piston and the chamber are small, typically less than 60 νm. In order to measure this film thickness a fluorescent dye was added to the water passing through the meter, which was illuminated with UV light. Visible light images were captured with a digital camera and analysed to give a measure of the film thickness with an uncertainty of less than 7%. It is known that this method lacks precision unless careful calibration is undertaken. Methods to achieve this are discussed in the paper. The grey level values for a range of film thicknesses were calibrated in situ with six dye concentrations to select the most appropriate one for the range of liquid film thickness. Data obtained for the oscillating piston flowmeter demonstrate the value of the fluorescence technique. The method is useful, inexpensive and straightforward and can be extended to other applications where measurement of liquid film thickness is required. © 2011 IOP Publishing Ltd.

First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.