PubMed Health⌕ Search

Biomedical subjects

I-Feng W Kuo

Publications and source records attributed to I-Feng W Kuo.

7 recordsLinked to original sources

Molecular dynamics simulation of liquid water: hybrid density functionals.

The structure, dynamical, and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first-principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta-functional, four gradient-corrected functionals, and the local density and Hartree-Fock approximations. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self-diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and understructured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller average numbers of hydrogen bonds than pure density functionals but similar hydrogen bond populations. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than those of the corresponding pure density functionals.

Algorithms↗

Simulating fluid-phase equilibria of water from first principles.

Efficient Monte Carlo algorithms and a mixed-basis set electronic structure program were used to compute from first principles the vapor-liquid coexistence curve of water. A water representation based on the Becke-Lee-Yang-Parr exchange and correlation functionals yields a saturated liquid density of 900 kg/m3 at 323 K and normal boiling and critical temperatures of 350 and 550 K, respectively. An analysis of the structural and electronic properties of the saturated liquid phase shows an increase of the asymmetry of the local hydrogen-bonded structure despite the persistence of a 4-fold coordination and decreases of the molecular dipole moment and of the spread of the lowest unoccupied molecular orbital with increasing temperature.

Journal Article↗

Isobaric-isothermal monte carlo simulations from first principles: application to liquid water at ambient conditions.

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T=298 K and p=1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2 K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids can result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm-3 resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm-3, showing the sensitivity of the simulation outcome to this parameter.

Journal Article↗

Bonding in the superionic phase of water.

The predicted superionic phase of water is investigated via ab initio molecular dynamics at densities of 2.0--3.0 g/cc (34-115 GPa) along the 2000 K isotherm. We find that extremely rapid (superionic) diffusion of protons occurs in a fluid phase at pressures between 34 and 58 GPa. A transition to a stable body-centered cubic O lattice with superionic proton conductivity is observed between 70 and 75 GPa, a much higher pressure than suggested in prior work. We find that all molecular species at pressures greater than 75 GPa are too short lived to be classified as bound states. Up to 95 GPa, we find a solid superionic phase characterized by covalent O-H bonding. Above 95 GPa, a transient network phase is found characterized by symmetric O-H hydrogen bonding with nearly 50% covalent character. In addition, we describe a metastable superionic phase with quenched O disorder.

Biophysics↗

Dynamic ionization of water under extreme conditions.

Raman spectroscopy in a laser heated diamond anvil cell and first principles molecular dynamics simulations have been used to study water in the temperature range 300 to 1500 K and at pressures to 56 GPa. We find a substantial decrease in the intensity of the O-H stretch mode in the liquid phase with pressure, and a change in slope of the melting line at 47 GPa and 1000 K. Consistent with these observations, theoretical calculations show that water beyond 50 GPa is "dynamically ionized" in that it consists of very short-lived (<10 fs) H2O, H3O+, and OH- species, and also that the mobility of the oxygen ions decreases abruptly with pressure, while hydrogen ions remain very mobile. We suggest that this regime corresponds to a superionic state.

Journal Article↗

An ab initio molecular dynamics study of the aqueous liquid-vapor interface.

We present an ab initio molecular dynamics simulation of the aqueous liquid-vapor interface. Having successfully stabilized a region of bulk water in the center of a water slab, we were able to reproduce and further quantify the experimentally observed abundance of surface "acceptor-only"(19%) and "single-donor"(66%) moieties as well as substantial surface relaxation approaching the liquid-vapor interface. Examination of the orientational dynamics points to a faster relaxation in the interfacial region. Furthermore, the average value of the dipole decreases and the average value of the highest occupied molecular orbital for each water molecule increases approaching the liquid-vapor interface. Our results support the idea that the surface contains, on average, far more reactive states than the bulk.

Journal Article↗