Isoviscosity lines and the liquid-glass transition in simple liquids

Properties of aqueous 1,4-dioxane solution via molecular dynamics

Polyethers are promising compounds for the creation of electrochemical energy storage systems. The molecular dynamics method can facilitate the search of compounds that have the most potential. However, the application of this method requires verification of the force fields. We perform molecular dynamics calculations of the physical properties of the aqueous 1,4-dioxane solution (density, enthalpy of mixing, and viscosity) and compare them to the available experimental data. In addition, we confirm the idea that the solution structure depends on the dioxane molar fraction, proposed in the experiment of Takamuku et al. [J. Mol. Liq. 83(1-3), 163-177 (1999)]. The hydrogen bonds between dioxane and water are analyzed. The correlation between the excess viscosity and enthalpy of mixing is demonstrated.

Transport coefficients of model lubricants up to 400 MPa from molecular dynamics

In this paper, the predictive power of molecular dynamics methods is demonstrated for the cases of model paraffinic and aromatic lubricant liquids at pressures up to 400 MPa. The shear viscosity and self-diffusion coefficients are calculated for 2,2,4-trimethylpentane (C₈H₁₈) at 298 K and 1,1-diphenylethane (C₁₄H₁₄) at 333 K. Three force fields with different levels of accuracy are compared by the ability to predict the experimental data. The Stokes-Einstein correlation between viscosity and self-diffusion is demonstrated for both compounds.

Contributions of force field interaction forms to Green-Kubo viscosity integral in n-alkane case

Decades of molecular simulation history proved that the Green-Kubo method for shear viscosity converges without any problems in atomic and simple molecular liquids, unlike liquids with high values of viscosity. In the case of highly viscous liquids, the time decomposition method was developed in 2015 by Maginn and co-authors [Y. Zhang, A. Otani, and E. J. Maginn, J. Chem. Theory Comput. 11, 3537-3546 (2015)] which allows us to improve the convergence of the Green-Kubo integral. In this paper, the contributions of intramolecular and intermolecular force field parts to the viscosity integral are discovered to gain the understanding of the Green-Kubo method. The n-alkanes from n-ethane to n-pentane at 330 K in the optimized potentials for liquid simulations-all atom force field are used as reference models. The dependencies of these contributions and decay times of the corresponding correlation functions on the chain length are observed. The nonequilibrium simulations are carried out to verify the Green-Kubo results. The obtained values of viscosity are compared with experimental data.

Properties of liquid iron along the melting line up to Earth-core pressures

We report a molecular dynamics study of the transport coefficients and the infinite frequency shear modulus of liquid iron at high temperatures and high pressures. We observe a simultaneous rise of both the shear viscosity and the diffusion coefficient along the melting line and estimate whether liquid iron can vitrify under Earth-core conditions. We show that in the conditions of the model studied in our work iron demonstrates a moderate increase of viscosity along the melting line. It is also demonstrated that at the limit of high temperatures and high pressures the liquid iron behaves similarly to the soft sphere system with exponent n ≈ 4.6.