Melting conditions and entropies of superionic water ice: Free-energy calculations based on hybrid solid/liquid reference systems

Superionic (SI) water ices-high-temperature, high-pressure phases of water in which oxygen ions occupy a regular crystal lattice whereas the protons flow in a liquid-like manner-have attracted a growing amount of attention over the past few years, in particular due to their possible role in the magnetic anomalies of the ice giants Neptune and Uranus. In this paper, we consider the calculation of the free energies of such phases, exploring hybrid reference systems consisting of a combination of an Einstein solid for the oxygen ions occupying a crystal lattice and a Uhlenbeck-Ford potential for the protonic fluid that avoids irregularities associated with possible particle overlaps. Applying this approach to a recent neural-network potential-energy landscape for SI water ice, we compute Gibbs free energies as a function of temperature for the SI fcc and liquid phases to determine the melting temperature T_m at 340 GPa. The results are consistent with previous estimates and indicate that the entropy difference between both phases is comparatively small, in particular due to the large amplitude of vibration of the oxygen ions in the fcc phase at the melting temperature.

Optimizing the driving function for nonequilibrium free-energy calculations in the linear regime: a variational approach

We consider the issue of optimizing linear-regime nonequilibrium simulations to estimate free-energy differences. In particular, we focus on the problem of finding the best-possible driving function lambda(t) that, for a given thermodynamic path, simulation algorithm, and amount of computational effort, minimizes dissipation. From the fluctuation-dissipation theorem it follows that, in the linear-response regime, the dissipation is controlled by the magnitude and characteristic correlation time of the equilibrium fluctuations in the driving force. As a result, the problem of finding the optimal switching scheme involves the solution of a standard problem in variational calculus: the minimization of a functional with respect to the switching function. In practice, the minimization involves solving the associated Euler-Lagrange equation subject to a set of boundary conditions. As a demonstration we apply the approach to the simple, yet illustrative problem of computing the free-energy difference between two classical harmonic oscillators with very different characteristic frequencies.

Transitions between disordered phases in supercooled liquid silicon

We have investigated the transitions between disordered phases in supercooled liquid silicon using computer simulations. The thermodynamic properties were directly obtained from the free energy, which was computed using the recently proposed reversible scaling method. The calculated free energies of the crystalline and liquid phases of silicon at zero pressure, obtained using the environment dependent interatomic potential, are in excellent agreement with the available experimental data. The results show that, at zero pressure, a weak first-order liquid-liquid transition occurs at 1135 K and a continuous liquid-amorphous transition takes place at 843 K. These results are consistent with the existence of a second critical point for the liquid-liquid transition at a negative pressure.