Study of thermal properties of the metastable supersaturated vapor with the integral equation method

Pressure, excess chemical potential, and excess free energy data for different densities of the supersaturated argon vapor at reduced temperatures from 0.7 to 1.2 are obtained by solving the integral equation with perturbation correction to the radial distribution function [F. Lado, Phys. Rev. 135, A1013 (1964)]. For those state points where there is no solution, the integral equation is solved with the interaction between argon atoms modeled by Lennard-Jones potential plus a repulsive potential with one controlling parameter, alpha exp(-rsigma) and in the end, all the thermal properties are mapped back to the alpha=0 case. Our pressure data and the spinodal obtained from the current method are compared with a molecular dynamics simulation study [A. Linhart et al., J. Chem. Phys. 122, 144506 (2005)] of the same system.     

The free energy of the metastable supersaturated vapor via restricted ensemble simulations. II. Effects of constraints and comparison with molecular dynamics simulations

Extensive restricted canonical ensemble Monte Carlo simulations [D. S. Corti and P. Debenedetti, Chem. Eng. Sci. 49, 2717 (1994)] were performed. Pressure, excess chemical potential, and excess free energy with respect to ideal gas data were obtained at different densities of the supersaturated Lennard-Jones (LJ) vapor at reduced temperatures from 0.7 to 1.0. Among different constraints imposed on the system studied, the one with the local minimum of the excess free energy was taken to be the approximated equilibrium state of the metastable LJ vapor. Also, a comparison of our results with molecular dynamic simulations [A. Linhart et al., J. Chem. Phys. 122, 144506 (2005)] was made.     

Chemical potential perturbation: a method to predict chemical potentials in periodic molecular simulations

A new method, called chemical potential perturbation (CPP), has been developed to predict the chemical potential as a function of density in periodic molecular simulations. The CPP method applies a spatially varying external force field to the simulation, causing the density to depend upon position in the simulation cell. Following equilibration the homogeneous (uniform or bulk) chemical potential as a function of density can be determined relative to some reference state after correcting for the effects of the inhomogeneity of the system. We compare three different methods of approximating this correction. The first method uses the van der Waals density gradient theory to approximate the inhomogeneous Helmholtz free energy density. The second method uses the local pressure tensor to approximate the homogeneous pressure. The third method uses the Triezenberg-Zwanzig definition of surface tension to approximate the inhomogeneous free energy density. If desired, the homogeneous pressure and Helmholtz free energy can also be predicted by the new method, as well as binodal and spinodal densities of a two-phase fluid region. The CPP method is tested using a Lennard-Jones (LJ) fluid at vapor, liquid, two-phase, and supercritical conditions. Satisfactory agreement is found between the CPP method and an LJ equation of state. The efficiency of the CPP method is compared to that for Widom's method under the tested conditions. In particular, the new method works well for dense fluids where Widom's method starts to fail.     

Investigation of excess adsorption, solvation force, and plate-fluid interfacial tension for Lennard-Jones fluid confined in slit pores

The excess Helmholtz free energy functional is formulated in terms of a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)] for a short ranged repulsion and a first-order mean-spherical approximation theory [Y. P. Tang, J. Chem. Phys. 118, 4140 (2003)] for a long ranged attraction. Within the framework of the density functional theory, the density profile, excess adsorption, solvation force, and plate-fluid interfacial tension of a Lennard-Jones fluid confined in slit pores are predicted, and the results agree well with the simulation data. The phase equilibria inside the slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal, and the plate-fluid interfacial tensions at equilibrium states are predicted consequently.     

The evaporation/condensation transition of liquid droplets

The condensation of a supersaturated vapor enclosed in a finite system is considered. A phenomenological analysis reveals that the vapor is found to be stable at densities well above coexistence. The system size at which the supersaturated vapor condenses into a droplet is found to be governed by a typical length scale which depends on the coexistence densities, temperature and surface tension. When fluctuations are neglected, the chemical potential is seen to show a discontinuity at an effective spinodal point, where the inhomogeneous state becomes more stable than the homogeneous state. If fluctuations are taken into account, the transition is rounded, but the slope of the chemical potential versus density isotherm develops a discontinuity in the thermodynamic limit. In order to test the theoretical predictions, we perform a simulation study of droplet condensation for a Lennard-Jones fluid and obtain loops in the chemical potential versus density and pressure. By computing probability distributions for the cluster size, chemical potential, and internal energy, we confirm that the effective spinodal point may be identified with the occurrence of a first order phase transition, resulting in the condensation of a droplet. An accurate equation of state is employed in order to estimate the droplet size and the coexisting vapor density and good quantitative agreement with the simulation data is obtained. The results highlight the need of an accurate equation of state data for the Laplace equation to have predictive power.     

Spinodal decomposition and droplets entrapment in monotectic solidification

In this article, we present two models to simulate solidification morphologies in monotectic alloys. With the first model, we investigate the morphological evolution under the influence of spinodal decomposition. The model requires that a gradient energy contribution for the concentration field should be incorporated, in order to stabilize phase separation when the liquid concentration is inside the region of miscibility gap. The free energy of the system in this model is derived from direct interpolation of the bulk energy densities. This, however, results in simulation regions in nanometer scale due to contributions from the chemical free energy of the system to the total surface excess. With the second model, our purpose is to develop a phase-field model to simulate scales that are larger than nanometer, where the departures from equilibrium are very small resulting in phase concentrations outside the spinodal region. In view of this, we exclude the concentration gradient contribution to the grand chemical potential functional, and develop a model based on [M. Plapp, Phys. Rev. E 84, 031601 (2011); A. Choudhury and B. Nestler, Phys. Rev. E 85, 021602 (2012)]. The advantage is that the free energy excess across the interface at equilibrium disappears, and hence it is easier to derive the required surface energies with higher interface widths. Due to this benefit, we employ the method to simulate the dynamic entrapment process in the monotectic reaction and study the influence of liquid(1) - liquid(2) surface energy and undercooling on the entrapment process.     

Prediction of phase behavior of nanoconfined Lennard-Jones fluids with density functional theory based on the first-order mean spherical approximation

The recently proposed first-order mean spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to study the phase behavior of nanoconfined Lennard-Jones fluids, which is consistent with the phase equilibria calculation of the corresponding bulk fluid. With a combination of fundamental measure theory, FMSA provides Helmholtz free energy and direct correlation function to formulate density functional theory, which implementation is as easy as the mean-field theory. Following previous success in predicting density profiles inside slit pores, this work is focused specially on the vapor-liquid equilibrium of the Lennard-Jones fluids inside these pores. It is found that outside the critical region FMSA predicts well the equilibrium diagram of slit pores with the sizes of 5.0, 7.5, and 10 molecular diameters by comparing with available computer simulation data. As a quantitative method, FMSA can be treated as an extension from its bulk calculation, while the mean-field theory is only qualitative, as its bulk version.     

Monte Carlo simulation study of droplet nucleation

A new rigorous Monte Carlo simulation approach is employed to study nucleation barriers for droplets in Lennard-Jones fluid. Using the gauge cell method we generate the excess isotherm of critical clusters in the size range from two to six molecular diameters. The ghost field method is employed to compute the cluster free energy and the nucleation barrier with desired precision of (1-2)kT. Based on quantitative results obtained by Monte Carlo simulations, we access the limits of applicability of the capillarity approximation of the classical nucleation theory and the Tolman equation. We show that the capillarity approximation corrected for vapor nonideality and liquid compressibility provides a reasonable assessment for the size of critical clusters in Lennard-Jones fluid; however, its accuracy is not sufficient to predict the nucleation barriers for making practical estimates of the rate of nucleation. The established dependence of the droplet surface tension on the droplet size cannot be approximated by the Tolman equation for small droplets of radius less than four molecular diameters. We confirm the conclusion of ten Wolde and Frenkel [J. Chem. Phys. 109, 9901 (1998)] that integration of the normal component of the Irving-Kirkwood pressure tensor severely underestimates the nucleation barriers for small clusters.     

Molecular dynamics simulation in the grand canonical ensemble

An extended system Hamiltonian is proposed to perform molecular dynamics (MD) simulation in the grand canonical ensemble. The Hamiltonian is similar to the one proposed by Lynch and Pettitt (Lynch and Pettitt, J Chem Phys 1997, 107, 8594), which consists of the kinetic and potential energies for real and fractional particles as well as the kinetic and potential energy terms for material and heat reservoirs interacting with the system. We perform a nonlinear scaling of the potential energy parameters of the fractional particle, as well as its mass to vary the number of particles dynamically. On the basis of the equations of motion derived from this Hamiltonian, an algorithm has been proposed for MD simulation at constant chemical potential. The algorithm has been tested for the ideal gas, for the Lennard-Jones fluid over a wide range of temperatures and densities, and for water. The results for the low-density Lennard-Jones fluid are compared with the predictions from a truncated virial equation of state. In the case of the dense Lennard-Jones fluid and water our predicted results are compared with the results reported using other available methods for the calculation of the chemical potential. The method is also applied to the case of vapor-liquid coexistence point predictions.     

Theoretical and computational investigations on thermodynamic properties, effective site diameters, and molecular free volume of carbon disulfide fluid

A newly proposed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was extended to polyatomics and applied to compute the density and temperature dependence of the effective site diameters of carbon disulfide fluids. The generic van der Waals (GvdW) theory was also extended to polyatomics in order to calculate the GvdW parameters and the molecular free volume using the effective site diameters as the repulsion-attraction separation distance. A three-site Lennard-Jones potential available in the literature was slightly modified and used in Monte Carlo simulations to obtain the functions appearing in the effective site diameter and GvdW expressions. The interaction potential was examined to reproduce the fluid phase thermodynamic properties using Gibbs ensemble Monte Carlo simulations and also the equation of state in the liquid phase using NVT Monte Carlo (NVT-MC) simulations. Comparison between the simulation results and experimental data shows excellent agreement for the densities of the coexisting phases, the vapor pressure, properties of the predicted critical point, and the equation of state. NVT-MC simulations were performed over a wide range of densities and temperatures in sub- and supercritical regions to compute the effective site diameters, the GvdW parameters, and the molecular free volume. The molecular structure in terms of the site-site pair correlation functions, the density dependence of the effective site diameters, and the density and temperature dependence of the GvdW parameters and molecular free volume were studied and discussed. The GvdW parameters were fitted to empirical expressions as a function of density and temperature. The computed molecular free volume will be used in future investigations to study the transport properties of carbon disulfide.     

Wall-liquid and wall-crystal interfacial free energies via thermodynamic integration: A molecular dynamics simulation study

A method is proposed to compute the interfacial free energy of a Lennard-Jones system in contact with a structured wall by molecular dynamics simulation. Both the bulk liquid and bulk face-centered-cubic crystal phase along the (111) orientation are considered. Our approach is based on a thermodynamic integration scheme where first the bulk Lennard-Jones system is reversibly transformed to a state where it interacts with a structureless flat wall. In a second step, the flat structureless wall is reversibly transformed into an atomistic wall with crystalline structure. The dependence of the interfacial free energy on various parameters such as the wall potential, the density and orientation of the wall is investigated. The conditions are indicated under which a Lennard-Jones crystal partially wets a flat wall.     

Determination of the solid-fluid coexistence of the n - 6 Lennard-Jones system from free energy calculations

The solid-fluid coexistence properties of the n - 6 Lennard-Jones system, n from 7 to 12, are reported. The procedure relies on determining Helmholtz free energy curves as a function of volume for each phase independently, from several NVT simulations, and then connecting it to points of known absolute free energy. For n = 12 this requires connecting the simulated points to states of very low densities on the liquid phase, and to a harmonic crystal for the solid phase, which involves many extra simulations for each temperature. For the reference points of the remaining systems, however, the free energy at a given density and temperature can be calculated relative to the n = 12 system. The method presented here involves a generalization of the multiple histogram method to combine simulations performed with different potentials, provided they visit overlapping regions of the phase space, and allows for a precise calculation of relative free energies. The densities, free energies, average potential energies, pressure, and chemical potential at coexistence are presented for up to T? = 5.0 and new estimations of the triple points are given for the n - 6 Lennard-Jones system.     

Multiple free energies from a single simulation: extending enveloping distribution sampling to nonoverlapping phase-space distributions

A recently proposed method to obtain free energy differences for multiple end states from a single simulation of a reference state which was called enveloping distribution sampling (EDS) [J. Chem. Phys. 126, 184110 (2007)] is expanded to situations where the end state configuration space densities do not show overlap. It uses a reference state Hamiltonian suggested by Han in 1992 [Phys. Lett. A 165, 28 (1992)] in a molecular dynamics implementation. The method allows us to calculate multiple free energy differences "on the fly" from a single molecular dynamics simulation. The influence of the parameters on the accuracy and precision of the obtained free energy differences is investigated. A connection is established between the presented method and the Bennett acceptance ratio method. The method is applied to four two-state test systems (dipole inversion, van der Waals perturbation, charge inversion, and water to methanol conversion) and two multiple-state test systems [dipole inversion with five charging states and five (dis-)appearing water molecules]. Accurate results could be obtained for all test applications if the parameters of the reference state Hamiltonian were optimized according to a given algorithm. The deviations from the exact result or from an independent calculation were at most 0.6 kJ/mol. An accurate estimation of the free energy difference is always possible, independent of how different the end states are. However, the convergence times of the free energy differences are longer in cases where the end state configuration space densities do not show overlap [charge inversion, water to methanol conversion, (dis-)appearing water molecules] than in cases where the configuration space densities do show some overlap [(multiple) dipole inversion and van der Waals perturbation].     

Regarding convergence curve of virial expansion for the Lennard-Jones system

We calculate the convergence curve of the virial expansion for a Lennard-Jones system in the density-temperature plane using the approximate method based on the density expansion of the Ornstein-Zernike equation and the condition of thermodynamic consistency [J. Chem. Phys. 106, 6095 (1997)]. At subcritical temperatures, this curve is close to the binodal. At supercritical temperatures, the curve does not coincide with the freezing curve. In the latter case, the densities along the convergence line are distinctly smaller than the densities corresponding to the condensed state.     

Calculation of surface tension via area sampling

We examine the performance of several molecular simulation techniques aimed at evaluation of the surface tension through its thermodynamic definition. For all methods explored, the surface tension is calculated by approximating the change in Helmholtz free energy associated with a change in interfacial area through simulation of a liquid slab at constant particle number, volume, and temperature. The methods explored fall within three general classes: free-energy perturbation, the Bennett acceptance-ratio scheme, and the expanded ensemble technique. Calculations are performed for both the truncated Lennard-Jones and square-well fluids at select temperatures spaced along their respective liquid-vapor saturation lines. Overall, we find that Bennett and expanded ensemble approaches provide the best combination of accuracy and precision. All of the methods, when applied using sufficiently small area perturbation, generate equivalent results for the Lennard-Jones fluid. However, single-stage free-energy-perturbation methods and the closely related test-area technique recently introduced by Gloor et al. [J. Chem. Phys. 123, 134703 (2005)] generate surface tension values for the square-well fluid that are not consistent with those obtained from the more robust expanded ensemble and Bennett approaches, regardless of the size of the area perturbation. Single-stage perturbation methods fail also for the Lennard-Jones system when applied using large area perturbations. Here an analysis of phase-space overlap produces a quantitative explanation of the observed inaccuracy and shows that the satisfactory results obtained in these cases from the test-area method arise from a cancellation of errors that cannot be expected in general. We also briefly analyze the variation in method performance with respect to the adjustable parameters inherent to the techniques.     

Molecular dynamics study of structural and thermodynamic characteristics of nanodroplets of simple fluid

Equilibrium configurations of Lennard-Jones nanodroplets composed of 10–15000 spherically symmetric molecules placed in the center of a spherical container are studied at constant temperature by the molecular dynamics method. The distribution of local density is found and size dependences of density in the center of droplet, first coordination number, and energy surface tension coinciding for equimolecular dividing surface with specific excess free energy of droplet are studied. Radial distribution function is also determined. It is established that the passage of structural characteristics to their macroscopic values is observed for droplets containing as little as about 300 molecules, while, for energy surface tension, analogous passage for energy surface tension occurs for droplets containing 700–6000 molecules.     

Competing interactions in two dimensional Coulomb systems: surface charge heterogeneities in coassembled cationic-anionic incompatible mixtures

A binary mixture of oppositely charged components confined to a plane such as cationic and anionic lipid bilayers may exhibit local segregation. The relative strengths of the net short range interactions, which favors macroscopic segregation, and the long range electrostatic interactions, which favors mixing, determine the length scale of the finite size or microphase segregation. The free energy of the system can be examined analytically in two separate regimes, when considering small density fluctuations at high temperatures and when considering the periodic ordering of the system at low temperatures [F. J. Solis, S. I. Stupp, and M. Olvera de la Cruz, J. Chem. Phys. 122, 054905 (2005)]. A simple molecular dynamics simulation of oppositely charged monomers, interacting with a short range Lennard-Jones potential and confined to a two dimensional plane, is examined at different strengths of short and long range interactions. The system exhibits well-defined domains that can be characterized by their periodic length scale as well as the orientational ordering of their interfaces. By adding salt, the ordering of the domains disappears and the mixture macroscopically phase segregates in agreement with analytical predictions.     

Vapor-liquid interfacial properties of rigid-linear Lennard-Jones chains

We have obtained the interfacial properties of short rigid-linear chains formed from tangentially bonded Lennard-Jones monomeric units from direct simulation of the vapour-liquid interface. The full long-range tails of the potential are accounted for by means of an improved version of the inhomogeneous long-range corrections of Janec?ek [J. Phys. Chem. B 110, 6264-6269 (2006)] proposed recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] valid for spherical as well as for rigid and flexible molecular systems. Three different model systems comprising of 3, 4, and 5 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test-area method. In addition to the surface tension, we also obtain density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length and rigidity on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The surface tension has been scaled by critical properties and represented as a function of the difference between coexistence densities relative to the critical density.     

Thermodynamic behavior and structural properties of an aqueous sodium chloride solution upon supercooling

We present the results of a molecular dynamics simulation study of thermodynamic and structural properties upon supercooling of a low concentration sodium chloride solution in TIP4P water and the comparison with the corresponding bulk quantities. We study the isotherms and the isochores for both the aqueous solution and bulk water. The comparison of the phase diagrams shows that thermodynamic properties of the solution are not merely shifted with respect to the bulk. Moreover, from the analysis of the thermodynamic curves, both the spinodal line and the temperatures of maximum density curve can be calculated. The spinodal line appears not to be influenced by the presence of ions at the chosen concentration, while the temperatures of maximum density curve displays both a mild shift in temperature and a shape modification with respect to bulk. Signatures of the presence of a liquid-liquid critical point are found in the aqueous solution. By analyzing the water-ion radial distribution functions of the aqueous solution, we observe that upon changing density, structural modifications appear close to the spinodal. For low temperatures, additional modifications appear also for densities close to that corresponding to a low density configurational energy minimum.