Viscosity and stress autocorrelation function in supercooled water: a molecular dynamics study

Assessing the impact of valence asymmetry in ionic solutions and its consequences on the performance of supercapacitors

Molecular dynamics simulations were performed to describe the properties of hypothetical salt electrolytic solutions. The main focus of this work is the valence asymmetry, which in recent years has been considered an important aspect in the physical chemistry of aqueous electrolytes. In general, our results show that the structural, energetic, and dynamic properties respond differently to the asymmetry of ionic solutions, but in all cases, appreciable changes were observed. Graphene supercapacitors based on the investigated electrolytes were studied in light of their electrostatic properties. We observed that the electrode capacitances, positive and negative, were greatly influenced by the presence of cations in the electrical double layer of the negative electrode and by the absence of these cations, in the double layer of the positive electrode. In general, we assess that quantitative variations due to valence asymmetry may indeed be an important factor for the development of new and more efficient electrolytes.

Transport properties of a quasisymmetric binary nitrogen-oxygen mixture in the warm dense regime

Transport properties of mixtures in the warm dense matter (WDM) regime play an important role in natural astrophysics. However, a physical understanding of ionic transport properties in quasisymmetric liquid mixtures has remained elusive. Here, we present extensive ab initio molecular dynamics (AIMD) simulations on the ionic diffusion and viscosity of a quasisymmetric binary nitrogen-oxygen (N-O) mixture in a wide warm dense regime of 8-120 kK and 4.5-8.0 g/cm^{3}. Diffusion and viscosity of N-O mixtures with different compositions are obtained by using the Green-Kubo formula. Unlike asymmetric mixtures, the change of proportions in N-O mixtures slightly affects the viscosity and diffusion in the strong-coupling region. Furthermore, the AIMD results are used to build and verify a global pseudo-ion in jellium (PIJ) model for ionic transport calculations. The PIJ model succeeds in reproducing the transport properties of N-O mixtures where ionization has occurred, and provides a promising alternative approach to obtaining comparable results to AIMD simulations with relatively small computational costs. Our current results highlight the characteristic features of the quasisymmetric binary mixtures and demonstrate the applicability of the PIJ model in the WDM regime.

Salt-in-water and water-in-salt electrolytes: the effects of the asymmetry in cation and anion valence on their properties

We investigated the structural, dynamic, energetic, and electrostatic properties of electrolytes based on the ion pairs LiCl and Li₂SO₄. Atomistic molecular dynamics simulations were used to simulate these aqueous electrolytic solutions at two different concentrations 2 M (normal) and 21 M (superconcentrated, WiSE). The effects of the valence asymmetry of the Li₂SO₄ electrolyte were also discussed for both salt concentrations. Our results differ in the physical aspect of pure electrolytes, showing the drastic effect of high concentration, in particular on the viscosity, which is dramatically increased in WiSE. This is a consequence of their reduced ionic mobility and has a direct effect on ionic conductivity. Also, our results for graphene-based supercapacitors, as indicated by some experimental work, do not indicate any better performance of WiSEs over normal electrolytes. In fact, the differences in the total capacitance, due to the concentration of ions, presented by both electrolytes are negligible. The valence asymmetry can be clearly observed in some properties but for most of them its effects could not be quantified or isolated.

Bulk viscosity of liquid noble gases

An equation of state for the bulk viscosity of liquid noble gases is proposed. On the basis of dedicated equilibrium molecular dynamics simulations, a multi-mode relaxation ansatz is used to obtain precise bulk viscosity data over a wide range of liquid states. From this dataset, the equation of state emerges as a two-parametric power function with both parameters showing a conspicuous saturation behavior over temperature. After passing a temperature threshold, the bulk viscosity is found to vary significantly over density, a behavior that resembles the frequency response of a one pole low-pass filter. The proposed equation of state is in good agreement with available experimental sound attenuation data.

Benchmarking the effective one-component plasma model for warm dense neon and krypton within quantum molecular dynamics simulation

The effective one-component plasma (EOCP) model has provided an efficient approach to obtaining many important thermophysical parameters of hot dense matter [J. Clérouin, et al., Phys. Rev. Lett. 116, 115003 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.115003]. In this paper, we perform extensive quantum molecular dynamics (QMD) simulations to determine the equations of state, ionic structures, and ionic transport properties of neon and krypton within the warm dense matter (WDM) regime where the density (ρ) is up to 12 g/cm^{3} and the temperature (T) is up to 100 kK. The simulated data are then used as a benchmark to explicitly evaluate the EOCP and Yukawa models. It is found that, within present ρ-T regime, the EOCP model can excellently reproduce the diffusion and viscosity coefficients of neon and krypton due to the fact that this model defines a system which nearly reproduces the actual physical states of WDM. Therefore, the EOCP model may be a promising alternative approach to reasonably predicting the transport behaviors of matter in WDM regime at lower QMD computational cost. The evaluation of Yukawa model shows that the consideration of the energy level broadening effect in the average atom model is necessary. Finally, with the help of EOCP model, the Stokes-Einstein relationships about neon and krypton are discussed, and fruitful plasma parameters as well as a practical ρ-T-dependent formula of the effective coupling parameter are obtained. These results not only provide valuable information for future theoretical and experimental studies on dense neon and krypton but also reveal the applicability of the EOCP model and the limitation of the Yukawa model in WDM regime and further support the continuing search for a unified description of ionic transport in dense plasma.

Quantum molecular dynamics study on the proton exchange, ionic structures, and transport properties of warm dense hydrogen-deuterium mixtures

Comprehensive knowledge of physical properties such as equation of state (EOS), proton exchange, dynamic structures, diffusion coefficients, and viscosities of hydrogen-deuterium mixtures with densities from 0.1 to 5 g/cm^{3} and temperatures from 1 to 50 kK has been presented via quantum molecular dynamics (QMD) simulations. The existing multi-shock experimental EOS provides an important benchmark to evaluate exchange-correlation functionals. The comparison of simulations with experiments indicates that a nonlocal van der Waals density functional (vdW-DF1) produces excellent results. Fraction analysis of molecules using a weighted integral over pair distribution functions was performed. A dissociation diagram together with a boundary where the proton exchange (H_{2}+D_{2}⇌2HD) occurs was generated, which shows evidence that the HD molecules form as the H_{2} and D_{2} molecules are almost 50% dissociated. The mechanism of proton exchange can be interpreted as a process of dissociation followed by recombination. The ionic structures at extreme conditions were analyzed by the effective coordination number model. High-order cluster, circle, and chain structures can be founded in the strongly coupled warm dense regime. The present QMD diffusion coefficient and viscosity can be used to benchmark two analytical one-component plasma (OCP) models: the Coulomb and Yukawa OCP models.

Method to manage integration error in the Green-Kubo method

The Green-Kubo method is a commonly used approach for predicting transport properties in a system from equilibrium molecular dynamics simulations. The approach is founded on the fluctuation dissipation theorem and relates the property of interest to the lifetime of fluctuations in its thermodynamic driving potential. For heat transport, the lattice thermal conductivity is related to the integral of the autocorrelation of the instantaneous heat flux. A principal source of error in these calculations is that the autocorrelation function requires a long averaging time to reduce remnant noise. Integrating the noise in the tail of the autocorrelation function becomes conflated with physically important slow relaxation processes. In this paper we present a method to quantify the uncertainty on transport properties computed using the Green-Kubo formulation based on recognizing that the integrated noise is a random walk, with a growing envelope of uncertainty. By characterizing the noise we can choose integration conditions to best trade off systematic truncation error with unbiased integration noise, to minimize uncertainty for a given allocation of computational resources.

Effect of guests on the adsorption interaction between a hydrate cage and guests

A criterion is proposed to judge which guest can enter the cage through which face.

Effect of correlation on viscosity and diffusion in molecular-dynamics simulations

In the warm dense matter (WDM) regime, material properties like diffusion and viscosity can be obtained from lengthy quantum molecular dynamics simulations, where the quantum behavior of the electrons is represented using either Kohn-Sham or orbital-free density functional theory. To reduce the simulation duration, we fit the time dependence of the autocorrelation functions (ACFs) and then use the fit to find values of the diffusion and viscosity. This fitting procedure avoids noise in the long time behavior of the ACFs. We present a detailed analysis of the functional form used to fit the ACFs, which is always a more efficient means to obtain mass transport properties. We use the fits to estimate the statistical error of the transport properties. We apply this methodology to a dense correlated plasma of copper and a mixture of carbon and hydrogen. Both systems show structure in their ACFs and exhibit multiple time scales. The mixture contains different structural forms of the ACFs for each component in the mixture.

Lifetimes of cagelike water clusters immersed in bulk liquid water: A molecular dynamics study on gas hydrate nucleation mechanisms

Molecular dynamics simulations were performed to observe the evolution of cagelike water clusters immersed in bulk liquid water at 250 and 230 K. Totally, we considered four types of clusters--dodecahedral (5(12)) and tetrakaidecahedral (5(12)6(2)) cagelike water clusters filled with or without a methane molecule, respectively. The lifetimes of these clusters were calculated according to their Lindemann index (delta) using the criterion of delta> or =0.07. The lifetimes of the clusters at 230 K are longer than that at 250 K, and their ratios are the same as the ratio of structure relaxation times of bulk water at these temperatures. For both the filled and empty clusters, the lifetimes of 5(12)6(2) cagelike clusters are similar to that of 5(12) cagelike clusters. Although the methane molecules indeed make the filled cagelike water clusters live longer than the empty ones, the empty cagelike water clusters still have the chance of being long lived. These observations support the cluster nucleation hypothesis for the formation mechanisms of gas hydrates.

Comment on "Computation of the viscosity of a liquid from time averages of stress fluctuations"

In a recent paper, Hess and Evans [Phys. Rev. E 64, 011207 (2001)] propose a method different from the conventional Green-Kubo and Einstein methods to calculate viscosity in equilibrium molecular dynamics simulations. For a comparison, we calculate the shear viscosity of SPC/E water at 303 K using these three different methods. We find that the Hess-Evans method is not as good as the other two in practical application, especially for the fluids with high viscosity and complicated relaxation.