This paper continues the investigation of the exponentially repulsive EXP pair-potential system of Paper I with a focus on isomorphs in the low-temperature gas and liquid phases. As expected from the EXP systems strong virial potential-energy correlations, the systems reduced-unit structure and dynamics are isomorph invariant to a good approximation. Three methods for generating isomorphs are compared: the small-step method that is exact in the limit of small density changes and t
The exponentially repulsive EXP pair potential defines a system of particles in terms of which simple liquids quasiuniversality may be explained [A. K. Bacher et al., Nat. Commun. 5, 5424 (2014); J. C. Dyre, J. Phys. Condens. Matter 28, 323001 (2016)]. This paper and its companion present a detailed simulation study of the EXP system. Here we study how structure monitored via the radial distribution function and dynamics monitored via the mean-square displacement as a function of time evolve along the systems isotherms and isochores. The focus is on the gas and liquid phases, which are distinguished pragmatically by the absence or presence of a minimum in the radial distribution function above its first maximum. An NVU-based proof of quasiuniversality is presented, and quasiuniversality is illustrated by showing that the structure of the Lennard-Jones system at four selected state points is well approximated by those of EXP pair-potential systems with the same reduced diffusion constant. The companion paper studies the EXP systems isomorphs, focusing also on the gas and liquid phases.
To study the possibility of a fluid-fluid phase transition, we analyze a three-dimensional soft-core isotropic potential for a one-component system. We utilize two independent numerical approaches, (i) integral equation in the hypernetted-chain approximation and (ii) molecular dynamics simulations, and find for both approaches a fluid-fluid phase transition as well as the conventional gas-liquid critical point. We also study the possible existence of a triple point in the supercooled fluid phase at which three phases---gas, high-density fluid, and low-density fluid---coexist.
We study in this paper the possible existence of Roskilde-simple liquids and their isomorphs in a rough-wall nanoconfinement. Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant in suitable nondimensionalized units. Two model liquids using molecular dynamics computer simulations are considered: the single-component Lennard-Jones (LJ) liquid and the Kob-Andersen binary LJ mixture, both of which in the bulk phases are known to have isomorphs. Nanoconfinement is implemented by adopting a slit-pore geometry with fcc crystalline walls; this implies inhomogenous density profiles both parallel and perpendicular to the confining walls. Despite this fact and consistent with an earlier study [Ingebrigtsen et. al, Phys. Rev. Lett. 111, 235901 (2013)] we find that these nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good scaling of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm along the isomorphs. From this study, we conjecture that in experiments, Roskilde-simple liquids may exhibit isomorphs if confined in a suitable manner, for example with carbon nanotubes. Our study thus provides an alternative framework for understanding nanoconfined liquids.
An iterative Monte Carlo inversion method for the calculation of particle pair potentials from given particle pair correlations is proposed in this paper. The new method, which is best referred to as Iterative Ornstein-Zernike Inversion, represents a generalization and an improvement of the established Iterative Boltzmann Inversion technique [Reith, P{u}tz & M{u}ller-Plathe, J. Comput. Chem. 24, 1624 (2003)]. Our modification of Iterative Boltzmann Inversion consists of replacing the potential of mean force as an approximant for the pair potential with another, generally more accurate approximant that is based on a trial bridge function in the Ornstein-Zernike integral equation formalism. As an input, the new method requires the particle pair correlations both in real space and in the Fourier conjugate wavenumber space. An accelerated iteration method is included in the discussion, by which the required number of iterations can be greatly reduced below that of the simple Picard iteration that underlies most common implementations of Iterative Boltzmann Inversion. Comprehensive tests with various pair potentials show that the new method generally surpasses the Iterative Boltzmann Inversion method in terms of reliability of the numerical solution for the particle pair potential.
In a companion paper, we derived analytical expressions for the structure factor of the square-shoulder potential in a perturbative way around the high- and low-temperature regimes. Here, various physical properties of these solutions are derived. In particular, we investigate the large wave number sector, and relate it to the contact values of the pair-correlation function. Then, thermoelastic properties of the square-shoulder fluids are discussed.