Freezing and melting properties of the Lennard‐Jones system

Melting Is Well-Known, but Is It Also Well-Understood?

Contrary to continuous phase transitions, where renormalization group theory provides a general framework, for discontinuous phase transitions such a framework seems to be absent. Although the thermodynamics of the latter type of transitions is well-known and requires input from two phases, for melting a variety of one-phase theories and models based on solids has been proposed, as a generally accepted theory for liquids is (yet) missing. Each theory or model deals with a specific mechanism using typically one of the various defects (vacancies, interstitials, dislocations, interstitialcies) present in solids. Furthermore, recognizing that surfaces are often present, one distinguishes between mechanical or bulk melting and thermodynamic or surface-mediated melting. After providing the necessary preliminaries, we discuss both types of melting in relation to the various defects. Thereafter we deal with the effect of pressure on the melting process, followed by a discussion along the line of type of materials. Subsequently, some other aspects and approaches are dealt with. An attempt to put melting in perspective concludes this review.

Measuring many-body distribution functions in fluids using test-particle insertion

We derive a hierarchy of equations, which allow a general n-body distribution function to be measured by test-particle insertion of between 1 and n particles. We apply it to measure the pair and three-body distribution functions in a simple fluid using snapshots from Monte Carlo simulations in the grand canonical ensemble. The resulting distribution functions obtained from insertion methods are compared with the conventional distance-histogram method: the insertion approach is shown to overcome the drawbacks of the histogram method, offering enhanced structural resolution and a more straightforward normalization. At high particle densities, the insertion method starts breaking down, which can be delayed by utilizing the underlying hierarchical structure of the insertion method. Our method will be especially useful in characterizing the structure of inhomogeneous fluids and investigating closure approximations in liquid state theory.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.

Coincidence of the freezing and the onset of caging in hard sphere and Lennard-Jones fluids

In this article, we examine the collective particle dynamics, as expressed by the time correlation function of the longitudinal particle current density, of several different fluids in the vicinity of their freezing points/lines. We consider and compare results obtained by dynamic light scattering for a suspension of hard spheres and by molecular dynamics for fluids with hard sphere and Lennard-Jones interactions. The latter are performed along both an isotherm and an isochore. In all cases, we find a qualitative change in the collective dynamics, within the resolution of the data, when their respective freezing lines are crossed. We associate this change with the onset of caging. The new results for the Lennard-Jones fluid reported here confirm that the occurrence of caging, found previously for systems of hard spheres, is a more general feature that distinguishes a metastable fluid from one in thermodynamic equilibrium.

Collective modes in simple melts: Transition from soft spheres to the hard sphere limit

We study collective modes in a classical system of particles with repulsive inverse-power-law (IPL) interactions in the fluid phase, near the fluid-solid coexistence (IPL melts). The IPL exponent is varied from n = 10 to n = 100 to mimic the transition from moderately soft to hard-sphere-like interactions. We compare the longitudinal dispersion relations obtained using molecular dynamic (MD) simulations with those calculated using the quasi-crystalline approximation (QCA) and find that this simple theoretical approach becomes grossly inaccurate for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\,\gtrsim 20$$\end{document}n≳20. Similarly, conventional expressions for high-frequency (instantaneous) elastic moduli, predicting their divergence as n increases, are meaningless in this regime. Relations of the longitudinal and transverse elastic velocities of the QCA model to the adiabatic sound velocity, measured in MD simulations, are discussed for the regime where QCA is applicable. Two potentially useful freezing indicators for classical particle systems with steep repulsive interactions are discussed.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

The Lennard-Jones melting line and isomorphism

The location of the melting line (ML) of the Lennard-Jones (LJ) system and its associated physical properties are investigated using molecular dynamics computer simulation. The radial distribution function and the behavior of the repulsive and attractive parts of the potential energy indicate that the ML is not a single isomorph, but the isomorphic state evolves gradually with temperature, i.e., it is only "locally isomorphic." The state point dependence of the unitless isomorphic number, X̃, for a range of static and dynamical properties of the LJ system in the solid and fluid states, and for fluid argon, are also reported. The quantity X̃ typically varies most with state point in the vicinity of the triple point and approaches a plateau in the high density (temperature) limit along the ML.

Crystal growth kinetics in Lennard-Jones and Weeks-Chandler-Andersen systems along the solid-liquid coexistence line

Kinetics of crystal-growth is investigated along the solid-liquid coexistence line for the (100), (110), and (111) orientations of the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) fcc crystal-liquid interface, using non-equilibrium molecular dynamics simulations. A slowing down of the growth kinetics along the coexistence line is observed, which is due to the decrease of the melting enthalpy with increasing coexistence temperature and pressure. Other quantities such as the melting pressure and liquid self-diffusion coefficient have a comparatively lesser impact on the kinetic growth coefficient. Growth kinetics of the LJ and WCA potentials become similar at large values of the melting temperature and pressure, when both resemble a purely repulsive soft-sphere potential. Classical models of crystallization from the melt are in reasonable qualitative agreement with our simulation data. Finally, several one-phase empirical melting/freezing rules are studied with respect to their validity along the coexistence line.

Soft mean spherical approximation for dusty plasma liquids: One-component Yukawa systems with plasma shielding

The structure and thermodynamics of strongly coupled dusty plasmas are investigated with the soft mean spherical approximation. This integral theory approach is analytically solvable for Yukawa pair interactions yielding a closed-form solution for the direct correlation function. The pair correlation function, the structure factor, and basic thermodynamic quantities are calculated for a wide range of parameters. Exact consistency between the "energy"-"virial" thermodynamic routes and approximate consistency between the "energy"-"compressibility" paths is demonstrated. Comparison with extensive molecular dynamics results is carried out and a remarkable agreement from the Coulomb limit to the strongly screened limit is revealed. The soft mean spherical approximation is concluded to be particularly well suited for the study of dusty plasma liquids, uniquely combining simplicity and accuracy.

Thermodynamic properties and entropy scaling law for diffusivity in soft spheres

The purely repulsive soft-sphere system, where the interaction potential is inversely proportional to the pair separation raised to the power n, is considered. The Laplace transform technique is used to derive its thermodynamic properties in terms of the potential energy and its density derivative obtained from molecular dynamics simulations. The derived expressions provide an analytic framework with which to explore soft-sphere thermodynamics across the whole softness-density fluid domain. The trends in the isochoric and isobaric heat capacity, thermal expansion coefficient, isothermal and adiabatic bulk moduli, Grüneisen parameter, isothermal pressure, and the Joule-Thomson coefficient as a function of fluid density and potential softness are described using these formulas supplemented by the simulation-derived equation of state. At low densities a minimum in the isobaric heat capacity with density is found, which is a new feature for a purely repulsive pair interaction. The hard-sphere and n = 3 limits are obtained, and the low density limit specified analytically for any n is discussed. The softness dependence of calculated quantities indicates freezing criteria based on features of the radial distribution function or derived functions of it are not expected to be universal. A new and accurate formula linking the self-diffusion coefficient to the excess entropy for the entire fluid softness-density domain is proposed, which incorporates the kinetic theory solution for the low density limit and an entropy-dependent function in an exponential form. The thermodynamic properties (or their derivatives), structural quantities, and diffusion coefficient indicate that three regions specified by a convex, concave, and intermediate density dependence can be expected as a function of n, with a narrow transition region within the range 5 < n < 8.

Continuous distributions of charges: extensions of the one component plasma

The electrostatic interaction between finite charge distributions, ρ(r), in a neutralizing background is considered as an extension of the one component plasma (OCP) model of point charges. A general form for the interaction potential is obtained which can be applied to molecular theories of many simple charged fluids and mixtures and to the molecular dynamics (MD) simulation of such systems. The formalism is applied to the study of a fluid of Gaussian charges in a neutralizing background by MD simulation and using hypernetted-chain integral equation theory. The treatment of these interactions is extended to a periodic system using a Fourier Transform formulation and, for a rapidly decaying charge distribution, an application of the Ewald method. The contributions of the self-energy and neutralizing background to the system's energy are explicitly included in the formulation. Calculations reveal differences in behavior from the OCP model when the Wigner-Seitz radius is of order and less than the Gaussian charge density decay length. For certain parameter values these systems can exhibit a multiple occupancy crystalline phase at high density which undergoes re-entrant melting at higher density. An exploration of the effects of the various length scales of the system on the equation of state and radial distribution function is made.

NVU perspective on simple liquids' quasiuniversality

The last half-century of research into the structure, dynamics, and thermodynamics of simple liquids has revealed a number of approximate universalities. This paper argues that simple liquids' reduced-coordinate constant-potential-energy hypersurfaces constitute a quasiuniversal family of compact Riemannian manifolds parametrized by a single number; from this follows the quasiuniversalities.

Fluid-solid phase transitions in three-dimensional complex plasmas under microgravity conditions

Phase behavior of large three-dimensional (3D) complex plasma systems under microgravity conditions onboard the International Space Station is investigated. The neutral gas pressure is used as a control parameter to trigger phase changes. Detailed analysis of structural properties and evaluation of three different melting-freezing indicators reveal that complex plasmas can exhibit melting by increasing the gas pressure. Theoretical estimates of complex plasma parameters allow us to identify main factors responsible for the observed behavior. The location of phase states of the investigated systems on a relevant equilibrium phase diagram is estimated. Important differences between the melting process of 3D complex plasmas under microgravity conditions and that of flat 2D complex plasma crystals in ground based experiments are discussed.

Freezing and melting of 3D complex plasma structures under microgravity conditions driven by neutral gas pressure manipulation

Freezing and melting of large three-dimensional complex plasmas under microgravity conditions is investigated. The neutral gas pressure is used as a control parameter to trigger the phase changes: Complex plasma freezes (melts) by decreasing (increasing) the pressure. The evolution of complex plasma structural properties upon pressure variation is studied. Theoretical estimates allow us to identify the main factors responsible for the observed behavior.

Phase diagram of the Weeks-Chandler-Andersen potential from very low to high temperatures and pressures

A combination of two molecular simulation algorithms has been used to determine the solid-liquid coexistence of the Weeks-Chandler-Andersen (WCA) fluid from low temperatures up to very high temperatures. Values are reported for the coexistence pressure, temperature, energy, enthalpy change, and densities of both the liquid and solid phases. At very high temperatures, the coexistence pressure approaches the same 12th-power soft-sphere asymptote as the 12-6 Lennard-Jones potential. However, in contrast to the Lennard-Jones potential, which shows a discontinuity of pressure at low temperatures, the coexistence pressure of the WCA potential approaches the zero-temperature limit. Empirical relationships are determined to accurately reproduce the coexistence pressure and both solid and liquid phase densities from near zero temperature to very high temperatures. The simulation data are used to improve the accuracy of a WCA equation of state. The validity of common melting and freezing rules is tested.

Evaluation of phenomenological one-phase criteria for the melting and freezing of softly repulsive particles

We test the validity of some widely used phenomenological criteria for the localization of the fluid-solid transition thresholds against the phase diagrams of particles interacting through the exp-6, inverse-power-law, and Gaussian potentials. We find that one-phase rules give, on the whole, reliable estimates of freezing/melting points. The agreement is ordinarily better for a face-centered-cubic solid than for a body-centered-cubic crystal, even more so in the presence of a pressure-driven reentrant transition of the solid into a denser fluid phase, as found in the Gaussian-core model.

Phase diagram of complex fluids using an efficient integral equation method

We present an adaptive technique for the determination of the phase diagram of fluids within the integral equation theory. It enables an efficient and accurate systematic mapping of the thermodynamic space in order to construct the binodal and spinodal lines. Results are obtained with the thermodynamically consistent integral equation proposed by Sarkisov [J. Chem. Phys. 114, 9496 (2001)] within the tangent linear technique that yields an exact differentiation of correlation functions. The generality of the numerical approach is assessed by determining both the liquid-vapor coexistence and the critical parameters of the generalized Lennard-Jones (n,6) potentials with varying repulsive part, including the hard-sphere limit.

Three-particle correlations in simple liquids

We use videomicroscopy to follow the phase-space trajectory of a two-dimensional colloidal model liquid and calculate three-point correlation functions from the measured particle configurations. Approaching the fluid-solid transition by increasing the strength of the pair-interaction potential, one observes the gradual formation of a crystal-like local order due to triplet correlations, while being still deep inside the fluid phase. Furthermore, we show that in a strongly interacting system the Born-Green equation can be satisfied only with the full triplet correlation function but not with three-body distribution functions obtained from superposing pair correlations (Kirkwood superposition approximation).