Effect of confinement on the solid-liquid coexistence of Lennard-Jones fluid

Theoretical investigation on the solid-liquid phase transition of gallium through free energy analysis

CONTEXT

Gallium, renowned for its notably low melting point and unique property of becoming liquid at room temperature, is a valuable constituent in phase change materials. In this study, we investigate the solid-liquid phase transition of gallium using the modified embedded atom method (MEAM) potential. It addresses the technique to compute the free energy difference between the solid and liquid without using a reference state. We examine various thermodynamic and dynamic properties, including density, specific heat capacity, diffusivity, and radial distribution functions. We compute the coexistence temperature of the solid-liquid phase transitions of gallium from free energy analysis. This information is crucial for understanding the behavior of the material under different pressure conditions and can be valuable for various applications, such as materials processing and high-pressure studies. The analysis, findings, and insights of the present work will be of great significance to the broad scientific and engineering communities in the field of phase transformation of materials.

METHODS

A series of molecular dynamics(MD) simulations were conducted using the LAMMPS software packages. The gallium atoms are modeled using the modified embedded atom method (MEAM) potential. To accurately predict the solid-liquid phase transitions of gallium, we calculated free energy by employing the "constrained λ integration" method, coupled with multiple histogram reweighting (MHR). The solid-liquid coexistence line is determined through the Gibbs-Duhem integration technique.

Confined Quantum Hard Spheres

We present computer simulation and theoretical results for a system of N Quantum Hard Spheres (QHS) particles of diameter σ and mass m at temperature T, confined between parallel hard walls separated by a distance Hσ, within the range 1≤H≤∞. Semiclassical Monte Carlo computer simulations were performed adapted to a confined space, considering effects in terms of the density of particles ρ*=N/V, where V is the accessible volume, the inverse length H−1 and the de Broglie’s thermal wavelength λB=h/2πmkT, where k and h are the Boltzmann’s and Planck’s constants, respectively. For the case of extreme and maximum confinement, 0.5<H−1<1 and H−1=1, respectively, analytical results can be given based on an extension for quantum systems of the Helmholtz free energies for the corresponding classical systems.

On the Gibbs-Thomson equation for the crystallization of confined fluids

The Gibbs-Thomson (GT) equation describes the shift of the crystallization temperature for a confined fluid with respect to the bulk as a function of pore size. While this century old relation is successfully used to analyze experiments, its derivations found in the literature often rely on nucleation theory arguments (i.e., kinetics instead of thermodynamics) or fail to state their assumptions, therefore leading to similar but different expressions. Here, we revisit the derivation of the GT equation to clarify the system definition, corresponding thermodynamic ensemble, and assumptions made along the way. We also discuss the role of the thermodynamic conditions in the external reservoir on the final result. We then turn to numerical simulations of a model system to compute independently the various terms entering in the GT equation and compare the predictions of the latter with the melting temperatures determined under confinement by means of hyper-parallel tempering grand canonical Monte Carlo simulations. We highlight some difficulties related to the sampling of crystallization under confinement in simulations. Overall, despite its limitations, the GT equation may provide an interesting alternative route to predict the melting temperature in large pores using molecular simulations to evaluate the relevant quantities entering in this equation. This approach could, for example, be used to investigate the nanoscale capillary freezing of ionic liquids recently observed experimentally between the tip of an atomic force microscope and a substrate.

Liquid-solid and solid-solid phase transition of monolayer water: high-density rhombic monolayer ice

Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.