Phase-space path-integral representation of the quantum density of states: Monte Carlo simulation of strongly correlated soft-sphere fermions

The Wigner formulation of quantum mechanics is used to derive a path-integral representation of the quantum density of states (DOS) of strongly correlated fermions in the canonical ensemble. A path-integral Monte Carlo approach for the simulation of DOS and other thermodynamic functions is suggested. The derived Wigner function in the phase space resembles the Maxwell-Boltzmann distribution but allows for quantum effects. We consider a three-dimensional quantum system of strongly correlated soft-sphere fermions at different densities and temperatures. The calculated properties include the DOS, momentum distribution functions, spin-resolved radial distribution functions, potentials of mean force, and related energy levels obtained from the Bohr-Sommerfeld condition. We observe sharp peaks on DOS and momentum distribution curves, which are explained by the appearance of fermionic bound states.

Exchange-correlation bound states of the triplet soft-sphere fermions by path-integral Monte Carlo simulations

Path-integral Monte Carlo simulations in the Wigner approach to quantum mechanics has been applied to calculate momentum and spin-resolved radial distribution functions of the strongly correlated soft-sphere quantum fermions. The obtained spin-resolved radial distribution functions demonstrate arising triplet clusters of fermions, that is the consequence of the interference of exchange and interparticle interactions. The semiclassical analysis in the framework of the Bohr-Sommerfeld quantization condition, applied to the potential of the mean force corresponding to the same-spin radial distribution functions, allows to detect exchange-correlation bound states in triplet clusters and to estimate corresponding averaged energy levels. The obtained momentum distribution functions demonstrate the narrow sharp separated peaks corresponding to bound states and disturbing the Maxwellian distribution.

Phase equilibria of symmetric Lennard-Jones mixtures and a look at the transport properties near the upper critical solution temperature

This study investigates phase equilibria and transport properties of five symmetric binary Lennard-Jones mixtures using molecular simulation and equation of state models. The mixtures are selected for their representation of different types of phase behavior and the research contributes to the development of simulation techniques, mixture theories and understanding of thermophysical mixture properties. A novel method is introduced for determining the critical end point (CEP) and critical azeotropic end point (CAEP) by molecular simulation. The van der Waals one-fluid theory is assessed for its performance in conjunction with Lennard-Jones equation of state models, while addressing different phase equilibrium types simultaneously. An empirical correlation is introduced to account for deviations between the equation of state and simulation that arise when using the same binary interaction parameter. This study also investigates the influence of the liquid-liquid critical point on thermophysical properties, which are found to exhibit no significant anomalies or singularities. System-size effects of diffusion coefficients are addressed by extrapolating simulation data to the thermodynamic limit and applying analytical finite-size corrections.

Interfacial properties of binary azeotropic mixtures of simple fluids: Molecular dynamics simulation and density gradient theory

Interfacial properties of binary azeotropic mixtures of Lennard-Jones truncated and shifted fluids were studied by molecular dynamics (MD) simulation and density gradient theory (DGT) in combination with an equation of state. Three binary mixtures were investigated, which differ in the energetic cross interaction parameter that yields different types of azeotropic behavior. This study covers a wide temperature and composition range. Mixture A exhibits a heteroazeotrope at low temperatures, which changes to a low-boiling azeotrope at high temperatures, mixture B exhibits a low-boiling azeotrope, and mixture C exhibits a high-boiling azeotrope. The phase behavior and fluid interfacial properties as well as their relation were studied. Vapor-liquid, liquid-liquid, and vapor-liquid-liquid equilibria and interfaces were considered. Density profiles, the surface tension, the interfacial thickness, as well as the relative adsorption and enrichment of the components at the interface were studied. The results obtained from the two independent methods (MD and DGT) are overall in good agreement. The results provide insights into the relation of the phase behavior, particularly the azeotropic behavior, of simple fluid mixtures and the corresponding interfacial properties. Strong enrichment was found for the mixture with a heteroazeotrope in the vicinity of the three-phase equilibrium, which is related to a wetting transition.

Bound States of the Exchange—Correlation Excitons in the Uniform Electron Gas by the Monte Carlo Simulations

The modified path integral representation of Wigner functions and the new Monte Carlo approach has been suggested to account for the impact of the interparticle interaction on the Pauli exclusion principle of fermions. This approach also allows to calculate the momentum distribution functions and to reduce the “sign problem” that is inaccessible to the standard path integral Monte Carlo methods. The obtained pair electron–electron distribution functions for the “uniform electron gas” demonstrate the short-range quantum ordering of electrons associated with exchange–correlation excitons. The exchange–correlation exciton is caused by the interaction of electrons with positively charged exchange holes and the excluded volume effect. The developed approach allows one to study the density–temperature range of the exciton arising, existence, and decay. Using the potential of the mean force and semiclassical Bohr–Sommerfeld quantization condition, we have demonstrated the existence of bound states disturbing the Maxwellian distribution and estimated their average energy levels. The exchange–correlation excitons have not been observed earlier in the standard path integral Monte Carlo simulations.

Characteristic Curves of the Lennard-Jones Fluid

Equations of state based on intermolecular potentials are often developed about the Lennard-Jones (LJ) potential. Many of such EOS have been proposed in the past. In this work, 20 LJ EOS were examined regarding their performance on Brown's characteristic curves and characteristic state points. Brown's characteristic curves are directly related to the virial coefficients at specific state points, which can be computed exactly from the intermolecular potential. Therefore, also the second and third virial coefficient of the LJ fluid were investigated. This approach allows a comparison of available LJ EOS at extreme conditions. Physically based, empirical, and semi-theoretical LJ EOS were examined. Most investigated LJ EOS exhibit some unphysical artifacts.

Phase Equilibria of Polydisperse Square-Well Chain Fluid Confined in Random Porous Media: TPT of Wertheim and Scaled Particle Theory

Extension of Wertheim's thermodynamic perturbation theory and its combination with scaled particle theory is proposed and applied to study the liquid-gas phase behavior of polydisperse hard-sphere square-well chain fluid confined in the random porous media. Thermodynamic properties of the reference system, represented by the hard-sphere square-well fluid in the matrix, are calculated using corresponding extension of the second-order Barker-Henderson perturbation theory. We study effects of polydispersity and confinement on the phase behavior of the system. While polydispersity causes increase of the region of phase coexistence due to the critical temperature increase, confinement decreases the values of both critical temperature and critical density making the region of phase coexistence smaller. This effect is enhanced with the increase of the size ratio of the fluid and matrix particles. The increase of the average chain length at fixed values of polydispersity and matrix density shifts the critical point to a higher temperature and a slightly lower density.

Self-Averaging Fluctuations in the Chaoticity of Simple Fluids

Bulk properties of equilibrium liquids are a manifestation of intermolecular forces. Here, we show how these forces imprint on dynamical fluctuations in the Lyapunov exponents for simple fluids with and without attractive forces. While the bulk of the spectrum is strongly self-averaging, the first Lyapunov exponent self-averages only weakly and at a rate that depends on the length scale of the intermolecular forces; short-range repulsive forces quantitatively dominate longer-range attractive forces, which act as a weak perturbation that slows the convergence to the thermodynamic limit. Regardless of intermolecular forces, the fluctuations in the Kolmogorov-Sinai entropy rate diverge, as one expects for an extensive quantity, and the spontaneous fluctuations of these dynamical observables obey fluctuation-dissipation-like relationships. Together, these results are a representation of the van der Waals picture of fluids and another lens through which we can view the liquid state.

Theoretical Rationale for a Thermodynamic Glass State

Using a chemical potential route, the square-well (SW) fluid model is solved in a quasichemical approximation (QCSW). At low temperatures, the liquid reaches a limiting density greater than the triple point density of a SW fluid but less than the equilibrium liquid-to-solid transition density for hard spheres. As this unique density is approached with decreasing temperature, the liquid entropy also approaches an asymptotic value, thus averting the "entropy catastrophe". Mean-field models in the van der Waals (VDW) genre fail to predict this type of behavior. In VDW models, attractive force contributions to the equation of state incorrectly diverge with decreasing temperature as 1/T, whereas those for the QCSW model asymptote to a fixed value. The QCSW model posits the intuitively pleasing idea that, at high densities, attractive contributions to the configurational energy begin to saturate well before zero temperature is reached. As a consequence, the force balance between repulsive and attractive forces stabilizes the liquid density, which thereafter becomes effectively independent of temperature. This fixed density in turn fixes all other density-dependent thermodynamic properties. These low-temperature, force-stabilized states are identified as glass states.

Exchange symmetry, fluctuation-compressibility relation, and thermodynamic potentials of quantum liquids

Liquid helium does not obey the Gibbs fluctuation-compressibility relation, which was noted more than six decades ago. However, still missing is a clear explanation of the reason for the deviation or the correct fluctuation-compressibility relation for the quantum liquid. Here we present the fluctuation-compressibility relation valid for any grand canonical system. Our result shows that the deviation from the Gibbs formula arises from a nonextensive part of thermodynamic potentials. The particle-exchange symmetry of many-body wave function of a strongly degenerate quantum gas is related to the thermodynamic extensivity of the system; a Bose gas does not always obey the Gibbs formula, while a Fermi gas does. Our fluctuation-compressibility relation works for classical systems as well as quantum systems. This work demonstrates that the application range of the Gibbs-Boltzmann statistical thermodynamics can be extended to encompass nonextensive open systems without introducing any postulate other than the principle of equal a priori probability.

TOWARDS A MOLECULAR THEORY FOR THE VAN DER WAALS–MAXWELL DESCRIPTION OF FLUID PHASE TRANSITIONS

We briefly review developments of theories for phase transitions of molecular fluids and mixtures, from semi-phenomenological approaches providing equations of state with adjustable parameters to first-principles microscopic methods qualitatively correct for a variety of molecular models with realistic interaction potentials. We further present the generalization of the van der Waals–Maxwell description of fluid phase diagrams to account for chemical specificities of polar molecular fluids, such as hydrogen bonding. Our theory uses the reference interaction site model (RISM) integral equation approach complemented with the new closure we have proposed (KH approximation), successful over a wide range of density from gas to liquid. The RISM/KH theory is applied to the known three-site models of water, methanol, and hydrogen fluoride. It qualitatively reproduces their vapor-liquid phase diagrams and the structure in the gas as well as liquid phases, including hydrogen bonding. Furthermore, phase transitions of water and methanol sorbed in nanoporous carbon aerogel are described by means of the replica generalization of the RISM approach we have developed. The changes as compared to the bulk fluids are in agreement with simulations and experiment. The RISM/KH theory is promising for description of phase transitions in various associating fluids, in particular, electrolyte as well as non-electrolyte solutions and ionic liquids.

Simulation of two- and three-dimensional dense-fluid shear flows via nonequilibrium molecular dynamics: comparison of time-and-space-averaged stresses from homogeneous Doll's and Sllod shear algorithms with those from boundary-driven shear

Homogeneous shear flows (with constant strainrate dv(x)/dy) are generated with the Doll's and Sllod algorithms and compared to corresponding inhomogeneous boundary-driven flows. We use one-, two-, and three-dimensional smooth-particle weight functions for computing instantaneous spatial averages. The nonlinear normal-stress differences are small, but significant, in both two and three space dimensions. In homogeneous systems the sign and magnitude of the shearplane stress difference, Pxx-Pyy, depend on both the thermostat type and the chosen shearflow algorithm. The Doll's and Sllod algorithms predict opposite signs for this normal-stress difference, with the Sllod approach definitely wrong, but somewhat closer to the (boundary-driven) truth. Neither of the homogeneous shear algorithms predicts the correct ordering of the kinetic temperatures: Txx > Tzz > Tyy.

Thermodynamic properties of model solids with short-ranged potentials from Monte Carlo simulations and perturbation theory

Monte Carlo simulations have been performed to determine the excess energy and the equation of state of fcc solids with Sutherland potentials for wide ranges of temperatures, densities, and effective potential ranges. The same quantities have been determined within a perturbative scheme by means of two procedures: (i) Monte Carlo simulations performed on the reference hard-sphere system and (ii) second-order Barker-Henderson perturbation theory. The aim was twofold: on the one hand, to test the capability of the "exact" MC-perturbation theory of reproducing the direct MC simulations and, on the other hand, the reliability of the Barker-Henderson perturbation theory, as compared with direct MC simulations and MC-perturbation theory, to determine the thermodynamic properties of these solids depending on temperature, density, and potential range. We have found that the simulation data for the excess energy obtained from the two procedures are in close agreement with each other. For the equation of state, the results from the MC-perturbation procedure also agree well with the direct MC simulations except for very low temperatures and extremely short-ranged potentials. Regarding the Barker-Henderson perturbation theory, we have found that in general the second-order approximation does not provide significant improvement over the first-order one.

Improvement on macroscopic compressibility approximation and beyond

A numerical procedure is proposed to extend the thermodynamic perturbation expansion (TPE) to a higher order. It is shown that the present second order term is superior to that due to a macroscopic compressibility approximation (MCA), a local compressibility approximation, and a superposition approximation by Barker and Henderson [Rev. Mod. Phys. 48, 587 (1976)]. Extensive model calculation and comparison with simulation data available in literature and supplied in the present report indicate that the present third order TPE is superior to a previous second order TPE based on the MCA, two previous perturbation theories, which are respectively based on an analytical mean spherical approximation for an Ornstein-Zernike equation, and an assumed explicit functional form for the Laplace transform of radial distribution function multiplied by radial distance, and a recent generalized van der Waals theory. The present critical temperature for a hard core attractive Yukawa fluid of varying range is in very good agreement with that due to a hierarchical reference theory. The present third order TPE is computationally far more modest than the self-consistent integral equation theory, and therefore is a viable alternative to use of the latter.

Thermodynamic perturbation theory in fluid statistical mechanics

A methodology is proposed that pushes the thermodynamic perturbation theory (TPT) from first order to higher order. The second-order correction is superior to a macroscopic compressibility (MC) approximation of Barker and Henderson. The present third-order TPT performs far better than the original first-order TPT and second-order TPT based on the MC approximation for many subfields in fluid statistical mechanics, such as predicting excess Helmholtz free energy, excess chemical potential, bulk pressure, gas-liquid coexistence, and solid-liquid equilibrium of very short-range potential fluids. A nonuniform version of the TPT is proposed; it is also shown that the nonuniform third-order TPT performs far better than the nonuniform first-order TPT in predicting density profile of fluids in critical region. The present report indicates that the TPT still can be a "universal" and accurate theoretical tool that has general applicability in fluid statistical mechanics, especially in soft-matter physics.

A general perturbation approach for equation of state development: Applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.

Vapor-to-droplet transition in a Lennard-Jones fluid: simulation study of nucleation barriers using the ghost field method

We report a comprehensive Monte Carlo (MC) simulation study of the vapor-to-droplet transition in Lennard-Jones fluid confined to a spherical container with repulsive walls, which is a case study system to investigate homogeneous nucleation. The focus is made on the application of a modified version of the ghost field method (Vishnyakov, A.; Neimark, A. V. J. Chem. Phys. 2003, 119, 9755) to calculate the nucleation barrier. This method allows one to build up a continuous trajectory of equilibrium states stabilized by the ghost field potential, which connects a reference droplet with a reference vapor state. Two computation schemes are employed for free energy calculations, direct thermodynamic integration along the constructed trajectory and umbrella sampling. The nucleation barriers and the size dependence of the surface tension are reported for droplets containing from 260 to 2000 molecules. The MC simulation study is complemented by a review of the simulation methods applied to computing the nucleation barriers and a detailed analysis of the vapor-to-droplet transition by means of the classical nucleation theory.