Yi Yao

Role of Charge Transfer in Water Diffusivity in Aqueous Ionic Solutions

We performed molecular dynamics simulations on four types of systems containing ion and solvating water. Two systems contained a cation (Na(+) or K(+)), and two other systems an anion (Cl(-) or I(-)). Classical molecular dynamics simulations were performed using three different force fields: a fixed charge force field, a polarizable force field that includes explicit polarization, and also a recently developed force field that includes polarization and charge transfer. These simulations were then compared to first-principles molecular dynamics simulations. While the first-principles simulations showed that the anions accelerated water translational diffusion, the cations slowed it down. In simulations with the classical force fields, only the force field that incorporates explicit charge transfer reproduced this ion-specific behavior. Additional simulations performed to understand the effect of charge transfer demonstrated that two competitive factors determine the behavior of water translational diffusion: the ions diminished charge accelerates water, while the net charge acquired by water either accelerates or slows down its dynamics. Our results show that charge transfer plays a crucial role in governing the water dynamics in aqueous ionic solutions.

Communication: Modeling of concentration dependent water diffusivity in ionic solutions: Role of intermolecular charge transfer

The translational diffusivity of water in solutions of alkali halide salts depends on the identity of ions, exhibiting dramatically different behavior even in solutions of similar salts of NaCl and KCl. The water diffusion coefficient decreases as the salt concentration increases in NaCl. Yet, in KCl solution, it slightly increases and remains above bulk value as salt concentration increases. Previous classical molecular dynamics simulations have failed to describe this important behavior even when polarizable models were used. Here, we show that inclusion of dynamical charge transfer among water molecules produces results in a quantitative agreement with experiments. Our results indicate that the concentration-dependent diffusivity reflects the importance of many-body effects among the water molecules in aqueous ionic solutions. Comparison with quantum mechanical calculations shows that a heterogeneous and extended distribution of charges on water molecules around the ions due to ion-water and also water-water charge transfer plays a very important role in controlling water diffusivity. Explicit inclusion of the charge transfer allows us to model accurately the difference in the concentration-dependent water diffusivity between Na(+) and K(+) ions in simulations, and it is likely to impact modeling of a wide range of systems for medical and technological applications.

Examining real-time time-dependent density functional theory nonequilibrium simulations for the calculation of electronic stopping power

In ion irradiation processes, electronic stopping power describes the energy transfer rate from the irradiating ion to the target materials electrons. Due to the scarcity and significant uncertainties in experimental electronic stopping power data for materials beyond simple solids, there has been growing interest in the use of first-principles theory for calculating electronic stopping power. In recent years, advances in high-performance computing have opened the door to fully first-principles nonequilibrium simulations based on real-time time-dependent density functional theory (RT-TDDFT). While it has been demonstrated that the RT-TDDFT approach is capable of predicting electronic stopping power for a wide range of condensed matter systems, there has yet to be an exhaustive examination of the physical and numerical approximations involved and their effects on the calculated stopping power. We discuss the results of such a study for crystalline silicon with protons as irradiating ions. We examine the influences of key approximations in RT-TDDFT nonequilibrium simulations on the calculated electronic stopping power, including approximations related to basis sets, finite size effects, exchange-correlation approximation, pseudopotentials, and more. Finally, we propose a simple and efficient correction scheme to account for the contribution from core-electron excitations to the stopping power, as it was found to be significant for large proton velocities.

Free Energy Profile of NaCl in Water: First-Principles Molecular Dynamics with SCAN and ωB97X-V Exchange–Correlation Functionals

Properties of water and aqueous ionic solutions are of great scientific interest because they play a central role in the atmosphere, biological environments, and various industrial processes. Employing two advanced exchange-correlation (XC) approximations, ωB97X-V and SCAN, in first-principles molecular dynamics simulations, we calculate the potential of mean force of NaCl in water as a function of the ion separation distance. Compared to the commonly used GGA-PBE functional, both of these XC functionals perform much better in simulating liquid water at room temperature for obtaining structural properties. The potential of mean force of NaCl in water exhibits two minima corresponding to two distinct types of ion pairing. ωB97X-V predicts that the contact ion pair is energetically more stable than the solvent-separated ion pair. The SCAN functional, however, predicts the opposite stability order, similarly to other XC functionals such as PBE. This is notable especially since classical molecular dynamics simulations with widely used force-field models predict greater stability for the contact ion pair. We also discuss how the electronic structures of water molecules and ions depend on the XC approximations. ωB97X-V and SCAN approximations noticeably improve the description of electron charge on Cl- ion in water while the charge on Na+ ion does not vary appreciably among the three XC functionals.

Plane-wave pseudopotential implementation and performance of SCAN meta-GGA exchange-correlation functional for extended systems

We present the implementation and performance of the strongly constrained and appropriately normed, SCAN, meta-GGA exchange-correlation (XC) approximation in the planewave-pseudopotential (PW-PP) formalism using the Troullier-Martins pseudopotential scheme. We studied its performance by applying the PW-PP implementation to several practical applications of interest in condensed matter sciences: (a) crystalline silicon and germanium, (b) martensitic phase transition energetics of phosphorene, and (c) a single water molecule physisorption on a graphene sheet. Given the much-improved accuracy over the GGA functionals and its relatively low computational cost compared to hybrid XC functionals, the SCAN functional is highly promising for various practical applications of density functional theory calculations for condensed matter systems. At same time, the SCAN meta-GGA functional appears to require more careful attention to numerical details. The meta-GGA functional shows more significant dependence on the fast Fourier transform grid, which is used for evaluating the XC potential in real space in the PW-PP formalism, than other more conventional GGA functionals do. Additionally, using pseudopotentials that are generated at a different/lower level of XC approximation could introduce noticeable errors in calculating some properties such as phase transition energetics.

Diffusion quantum Monte Carlo study of martensitic phase transition energetics: The case of phosphorene

Recent technical advances in dealing with finite-size errors make quantum Monte Carlo methods quite appealing for treating extended systems in electronic structure calculations, especially when commonly used density functional theory (DFT) methods might not be satisfactory. We present a theoretical study of martensitic phase transition energetics of a two-dimensional phosphorene by employing diffusion Monte Carlo (DMC) approach. The DMC calculation supports DFT prediction of having a rather diffusive barrier that is characterized by having two transition states, in addition to confirming that the so-called black and blue phases of phosphorene are essentially degenerate. At the same time, the DFT calculations do not provide the quantitative accuracy in describing the energy changes for the martensitic phase transition even when hybrid exchange-correlation functional is employed. We also discuss how mechanical strain influences the stabilities of the two phases of phosphorene.