Search

Due to their potential role in the peculiar geophysical properties of the ice giants Neptune and Uranus, there has been a growing interest in superionic (SI) phases of water ice. So far, however, little attention has been given to their mechanical properties, even though plastic deformation processes in the interiors of planets are known to affect long-term processes, such as plate tectonics and mantle convection. Here, using density functional theory calculations and machine learning techniques, we assess the mechanical response of high-pressure/temperature solid phases of water in terms of their ideal shear strength (ISS) and dislocation behavior. The ISS results are well described by the renormalized Frenkel model of ideal strength and indicate that the SI ices are expected to be highly ductile. This is further supported by deep neural network molecular dynamics simulations for the behavior of lattice dislocations for the SI face-centered cubic (fcc) phase. Dislocation velocity data indicate effective shear viscosities that are orders of magnitude smaller than that of Earth’s lower mantle, suggesting that the plastic flow of the internal icy layers in Neptune and Uranus may be significantly faster than previously foreseen.

The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κ

This article describes nonequilibrium techniques for the calculation of free energies of solids using molecular dynamics (MD) simulations. These methods provide an alternative to standard equilibrium thermodynamic integration methods and often present superior efficiency. Here we describe the implementation in the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) code of two specific nonequilibrium processes that allow the calculation of the free-energy difference between two different system Hamiltonians as well as the free-energy temperature dependence of a given Hamiltonian, respectively. The theory behind the methods is summarized, and we describe (including fragments of LAMMPS scripts) how the process parameters should be selected to obtain the best-possible efficiency in the calculations of free energies using nonequilibrium MD simulations. As an example of the application of the methods we present results related to polymorphic transitions for a classical potential model of iron.