Simulating vapor-liquid nucleation of water: A combined histogram-reweighting and aggregation-volume-bias Monte Carlo investigation for fixed-charge and polarizable models

The method of histogram-reweighting was integrated with a recently developed approach using aggregation-volume-bias Monte Carlo and self-adaptive umbrella sampling to develop the AVUS-HR algorithm that allows for exceedingly efficient calculations of nucleation properties over a wide range of thermodynamic conditions. Simulations were carried out for water using both fixed-charge and polarizable force fields belonging to the TIP4P family (namely, TIP4P, TIP4P-FQ, TIP4P-pol2, and TIP4P-pol3) to investigate the nucleation of water over a temperature range from 200 to 300 K and the concentration of water clusters in the atmosphere at elevations up to 15 km. It was found that the nucleation free energy barriers and atmospheric concentrations are extremely sensitive to the force field, albeit all of the models investigated in this study support the following general conclusions: (i) ice nucleation is not present under the thermodynamic conditions and cluster-size range investigated here, i.e., the critical nuclei possess liquidlike structures, and (ii) the atmospheric concentrations of water clusters under homogeneous conditions are very low with the mole fraction of hexamers being about 10(-10), a number probably too low to influence the solar radiation balance. Compared to the experimental data, the TIP4P-pol3 model yields the most accurate nucleation results, consistent with its excellent performance for the second virial coefficient and the minimum cluster energies.     

Molecular modelling of liquid crystal systems: An internal coordinate Monte Carlo approach

A Monte Carlo scheme is presented which is designed to provide a convenient mechanism to model accurately the internal molecular structure of liquid crystalline molecules. The technique stores atomic positions in terms of bond lengths, bond angles and dihedral angles within a Z-matrix, and the Monte Carlo scheme involves generating trial configurations from changes to the Z-matrix using the MM2 molecular mechanics potential to describe energy changes between different molecular conformations. The technique is applied to the liquid crystal molecule 4-n-pentyl-4'-cyanobiphenyl (5CB), and results are presented for the conformational populations and dihedral angle distributions of 5CB in the gas phase at 300 K. The effect of a nematic mean field on the distribution of molecular conformations is also examined via the addition of a conformation-dependent potential of mean torque to the internal energy.     

Freezing line of the Lennard-Jones fluid: a phase switch Monte Carlo study

We report a phase switch Monte Carlo (PSMC) method study of the freezing line of the Lennard-Jones (LJ) fluid. Our work generalizes to soft potentials the original application of the method to hard-sphere freezing and builds on a previous PSMC study of the LJ system by Errington [J. Chem. Phys. 120, 3130 (2004)]. The latter work is extended by tracing a large section of the Lennard-Jones freezing curve, the results of which we compare with a previous Gibbs-Duhem integration study. Additionally, we provide new background material regarding the statistical-mechanical basis of the PSMC method and extensive implementation details.     

Homogeneous nucleation in vapor-liquid phase transition of Lennard-Jones fluids: a density functional theory approach

Density functional theory (DFT) with square gradient approximation for the free energy functional and a model density profile are used to obtain an analytical expression for the size-dependent free energy of formation of a liquid drop from the vapor through the process of homogeneous nucleation, without invoking the approximations used in classical nucleation theory (CNT). The density of the liquid drop in this work is not the same as the bulk liquid density but it corresponds to minimum free energy of formation of the liquid drop. The theory is applied to study the nucleation phenomena from supersaturated vapor of Lennard-Jones fluid. The barrier height predicted by this theory is significantly lower than the same in CNT which is rather high. The density at the center of the small liquid drop as obtained through optimization is less than the bulk density which is in agreement with other earlier works. Also proposed is a sharp interface limit of the proposed DFT of nucleation, which is as simple as CNT but with a modified barrier height and this modified classical nucleation theory, as we call it, is shown to lead to improved results.     

Monte Carlo simulations of a liquid crystal copolymer in the solid state

The condensed phase of the alternating copolyester of p-hydroxybenzoic acid (HBA) and 2-hydroxy-6-naphthoic acid (HNA) is investigated by studying the room temperature packing arrangement of the copolymer chains. A molecular modeling methodology is employed with a Monte Carlo sampling of the configurational phase space. Realistic poly(HBA-alt-HNA) polymer chains are represented by an explicit atom representation of the HBA/HNA dimers. States are sampled from the NVT ensemble using a sampling scheme consisting of (1) valence and torsional variations, (2) rigid body rotations of the chain about the chain axis, and (3) rigid body translations of the chain. The effect of chain packing on the conformation of chains, as well as the relative intra- and intermolecular orientations of aromatic rings, is investigated. Correlation of chain positioning along the chain axis is dominated by aromatic rings maintaining a center-to-center plane of registry. These layers of aromatic units pack with a preference for edge-to-face orientations in a herringbone-type pattern and have an intermolecular ring angle between the pairs of aromatic rings in the unit cell that is ca. 68°. The aromatic rings, on average, are rotated 38° out from the b–c plane. The phenylene rings of these copolyesters are less restricted in their relative orientation in comparison to the naphthalene rings. Intramolecular orientational probability density distributions indicate a preference for staggering the successive aromatic rings along the chain, with a staggering angle of ca. 66°. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 727–741, 1998     

Self-similar theory of nonisothermal vapor condensation on a droplet growing in a vapor-gas medium

Rigorous self-similar solutions to the joint problems of vapor diffusion toward a droplet growing in a vapor-gas medium and the removal of heat released during vapor condensation are found. An equation for the temperature of a droplet ensuring the existence of a self-similar solution is derived. This equation sets the constancy of the temperature of a droplet throughout the time of its growth and unambiguously determines this temperature. In the case of the strong heat effects, when the rate of droplet growth decreases substantially, the analytical solution to this equation is obtained. This temperature coincides precisely with the temperature, which is established in the droplet at the diffusion regime of its growth. At the found droplet temperature, interconnected fields of vapor concentration and temperature of vapor-gas medium around the droplet are expressed through the initial (prior to the droplet nucleation) parameters of a vapor-gas medium. These parameters are used to express the dependence of the radius of a droplet on the time at the diffusion regime of its growth and the time required to establish the diffusion regime of droplet growth. The case of weak heat effects is also studied.     

Coagulation of nanoparticles in reverse micellar systems: a Monte Carlo model

The process of formation of nanoparticles obtained by mixing two micellized, aqueous solutions has been simulated using the Monte Carlo technique. The model includes the phenomena of finite nucleation, growth via intermicellar exchange, and coagulation of nanoparticles after their formation. Using the model, an exploratory study has been conducted to analyze whether the coagulation of nanoparticles is the reason for the formation of nanoparticles whose sizes are comparable to the size of the reverse micelles. The model explains the possible mechanism of coagulation of semiconductor nanoparticles formed within reverse micelles and its effect on the evolution of their size with time. The model is predictive in nature, and the simulation results compare well with those observed experimentally.     

Water vapor nucleation on an infinite substrate with complementary structure: 2. The influence of crystal defects

The simulation of the adsorption of water vapor on the infinite surface of a silver iodide crystal with regularly arranged crystal defects is performed by the Monte Carlo method. Long-range electrostatic and polarization interactions are taken into account by the Ewald method. The work of nucleus formation in the field of crystal defects is calculated by the method of bicanonical statistical ensemble. The interaction with a defect decreases the internal energy, the Gibbs energy, and the entropy of the nucleus. The rough surface of the silver iodide crystal exerts a larger stimulating effect on vapor nucleation than an ideally smooth surface. Saturation vapor pressure above the condensate layer formed on the rough surface can be several times lower than the pressure above the smooth surface. This effect is caused by the cooperative action of surface crystal defects. It may be expected that the surface with a characteristic size of rough elements of 15–20 Å is the most efficient for the formation of condensed phase nuclei. Single point crystal defects with extremely small sizes on the substrate surface do not exert a stimulating effect on the formation of the macroscopic condensed phase. Single defects with moderate sizes can lower the barrier of monomolecular film formation; however, such a route of the stimulation of ice-forming activity of the surface is less efficient than the cooperative action of the defects of rough surface.     

Can Monte Carlo detect the absence of ordering in a model liquid crystal?

We investigate the possibility of inferring the absence of an ordered phase using Monte Carlo simulations. The example we have chosen is that of a one dimensional Lebwohl-Lasher model, where an analytic solution is available. We argue that Monte Carlo can be of help even in this delicate sector notwithstanding the complications created by periodic boundary conditions.     

Hydration free energy of a Model Lennard-Jones solute particle: microscopic Monte Carlo simulation studies, and interpretation based on mesoscopic models

In this study, the hydration of a model Lennard-Jones solute particle and the analytical approximations of the free energy of hydration as functions of solute microscopic parameters are analyzed. The control parameters of the solute particle are the charge, the Lennard-Jones diameter, and also the potential well depth. The obtained multivariate free energy functions of hydration were parametrized based on Metropolis Monte Carlo simulations in the extended NpT ensemble, and interpreted based on mesoscopic solvation models proposed by Gallicchio and Levy [J. Comput. Chem. 25, 479 (2004)], and Wagoner and Baker [Proc. Natl. Acad. Sci. U.S.A. 103, 8331 (2006)]. Regarding the charge and the solute diameter, the dependence of the free energy on these parameters is in qualitative agreement with former studies. The role of the third parameter, the potential well depth not previously considered, appeared to be significant for sufficiently precise bivariate solvation free energy fits. The free energy fits for cations and neutral solute particles were merged, resulting in a compact manifold of the free energy of solvation. The free energy of hydration for anions forms two separate manifolds, which most likely results from an abrupt change of the coordination number when changing the size of the anion particle.     

Monte Carlo simulation for the formation of a mixed crystal from two solids in contact

The study focuses on nucleation and growth of a binary mixed crystal phase from two pure crystals in contact. Monte Carlo simulations of this process are conducted, with the dynamics proceeding via activated atom-vacancy exchanges. Intermolecular interactions, ranging up to next-nearest neighbors, are of size typical of hydrogen bonded systems. The process is driven by the formation of strong AB bonds at the expense of weaker AA and BB bonds. In the resulting model, the material is channeled and transported through the mixed phase crust along antiphase boundaries. The flow of molecules through the channels is directed, due to molecular energy lowering via gradual acquisition of an increasing number of nearest neighbors of the second species. On the other hand, defect motion is quasirandom. The model accounts partially for the t(1/alpha) (alpha>3) time dependence observed for conversion of nanoparticles of HBr dihydrate to monohydrate, by exposure to acid adsorbate.     

Monte Carlo versus molecular dynamics simulations in heterogeneous systems: an application to the n-pentane liquid-vapor interface

The Monte Carlo (MC) and molecular dynamics (MD) methodologies are now well established for computing equilibrium properties in homogeneous fluids. This is not yet the case for the direct simulation of two-phase systems, which exhibit nonuniformity of the density distribution across the interface. We have performed direct MC and MD simulations of the liquid-gas interface of n-pentane using a standard force-field model. We obtained density and pressure components profiles along the direction normal to the interface that can be very different, depending on the truncation and long range correction strategies. We discuss the influence on predicted properties of different potential truncation schemes implemented in both MC and MD simulations. We show that the MD and MC profiles can be made in agreement by using a Lennard-Jones potential truncated via a polynomial function that makes the first and second derivatives of the potential continuous at the cutoff distance. In this case however, the predicted thermodynamic properties (phase envelope, surface tension) deviate from experiments, because of the changes made in the potential. A further readjustment of the potential parameters is needed if one wants to use this method. We conclude that a straightforward use of bulk phase force fields in MD simulations may lead to some physical inconsistencies when computing interfacial properties.     

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.     

The ground state tunneling splitting and the zero point energy of malonaldehyde: a quantum Monte Carlo determination

Quantum dynamics calculations of the ground state tunneling splitting and of the zero point energy of malonaldehyde on the full dimensional potential energy surface proposed by Yagi et al. [J. Chem. Phys. 1154, 10647 (2001)] are reported. The exact diffusion Monte Carlo and the projection operator imaginary time spectral evolution methods are used to compute accurate benchmark results for this 21-dimensional ab initio potential energy surface. A tunneling splitting of 25.7+/-0.3 cm-1 is obtained, and the vibrational ground state energy is found to be 15 122+/-4 cm-1. Isotopic substitution of the tunneling hydrogen modifies the tunneling splitting down to 3.21+/-0.09 cm-1 and the vibrational ground state energy to 14 385+/-2 cm-1. The computed tunneling splittings are slightly higher than the experimental values as expected from the potential energy surface which slightly underestimates the barrier height, and they are slightly lower than the results from the instanton theory obtained using the same potential energy surface.     

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.     

2-D islanding of dodecane on an Au(111) surface: an investigation using he beam reflectivity and Monte Carlo modeling

Dodecane is deposited at submonolayer coverages onto an Au(111) surface forming two-dimensional (2-D) islands. The islands sublimate to a 2-D gas at higher substrate temperatures. We observe island formation and subsequent 2-D sublimation between substrate temperatures of 40 and 350 K, using low-energy helium reflectivity. A computer model of the submonolayer islanding process using Monte Carlo simulations shows significant agreement with experimental data and yields an intermolecular potential of 0.10+/-0.03 eV (about half that of the bulk substance) and a significantly higher corrugation potential of 0.3+/-0.1 eV.     

Kinetic model for binary homogeneous nucleation in the H2O-H2SO4 system: comparison with experiments and classical theory of nucleation

A kinetic model to predict nucleation rates in the sulfuric acid-water system is presented. It allows calculating steady-state nucleation rates and the corresponding time lag, using a direct solution of a system of kinetic equations that describe the populations of sub- and near-critical clusters. This kinetic model takes into account cluster-cluster collisions and decay of clusters into smaller clusters. The model results are compared with some predictions obtained with the classical nucleation theory (CNT) and also with available measurement data obtained in smog chambers or flow tubes. It is shown that in the case of slow nucleation processes, the kinetic model and the CNT as used by Shugard et al. [J. Chem. Phys. 75, 5298 (1974)] give the same results. However, in the case of intensive nucleation, a large part of the nucleation flux is due to cluster-cluster collisions and the CNT underestimates the nucleation rates.     

A grand canonical Monte Carlo study of capillary condensation in mesoporous media: effect of the pore morphology and topology

We study by means of Grand Canonical Monte Carlo simulations the condensation and evaporation of argon at 77 K in nanoporous silica media of different morphology or topology. For each porous material, our results are compared with data obtained for regular cylindrical pores. We show that both the filling and emptying mechanisms are significantly affected by the presence of a constriction. The simulation data for a constricted pore closed at one end reproduces the asymmetrical shape of the hysteresis loop that is observed for many real disordered porous materials. The adsorption process is a quasicontinuous mechanism that corresponds to the filling of the different parts of the porous material, cavity, and constriction. In contrast, the desorption branch for this pore closed at one end is brutal because the evaporation of Ar atoms confined in the largest cavity is triggered by the evaporation of the fluid confined in the constriction (which isolates the cavity from the gas reservoir). This evaporation process conforms to the classical picture of "pore blocking effect" proposed by Everett many years ago. We also simulate Ar adsorption in a disordered porous medium, which mimics a Vycor mesoporous silica glass. The adsorption isotherm for this disordered porous material having both topological and morphological defects presents the same features as that for the constricted pore (quasicontinuous adsorption and steep desorption process). However, the larger degree of disorder of the Vycor surface enhances these main characteristics. Finally, we show that the effect of the disorder, topological and/or morphological, leads to a significant lowering of the capillary condensation pressure compared to that for regular cylindrical nanopores. Also, our results suggest that confined fluids isolated from the bulk reservoir evaporate at a pressure driven by the smallest size of the pore.