On the Sampling Method for Grand Canonical Monte Carlo Simulations

Abstract

A bulk Lennard-Jones fluid was simulated using the grand canonical Monte Carlo method. Three different sampling methods were used in the transition matrix, namely the Metropolis, Barker and a third novel method. While it can be shown that the Metropolis method will give the most accurate ensemble averages in the limit of an infinitely long run, the new method termed “Modified Barker Sampling” (MBS), is shown to be superior for the runs of practical length for the particular system studied.     

Calculation of Vapour-Liquid Equilibria of Lennard-Jones Binary Systems by the Gibbs Ensemble Monte Carlo Simulation

Abstract

The Gibbs ensemble Monte Carlo simulation has been used to calculate vapour-liquid equilibria of a Lennard-Jones (LJ) binary mixture. The mixture studied is the LB-2-1 model which has been used in our previous calculations on PVT relation and density-dependent local composition. The P-x-y relation has been established at two different temperatures and used to determine vapour-liquid coexistence region in the PVTx space.     

Gibbs Ensemble Calculations with an Equation of State: An Application to Vapor-Liquid Equilibria

Abstract

A new modification of the Gibbs ensemble Monte Carlo computer simulation method for fluid phase equilibria is described. The modification is based on a thermodynamic model for the vapor phase, and uses an equation of state to account for the weak interactions between the vapor phase molecules. Reductions in the computational time by 30–40% as compared to the original Gibbs ensemble method are obtained. The algorithm is applied to Lennard-Jones - (12,6) fluids and their mixtures and the results are in good agreement with results obtained from simulations using the full Gibbs ensemble method.     

Smart Monte Carlo Algorithm for the Adsorption of Molecules at a Surface

Abstract

A modified grand canonical ensemble Monte Carlo (GCMC) technique has been developed to simulate adsorption isotherms for molecules on or near a surface. The speed and accuracy of the simulation is increased by using a non-uniform distribution function, related to the force field exerted by the surface and the current configuration, to generate coordinates for the creation of new particles in the simulation. With this method, isotherms are generated more efficiently than with current techniques in which the creation step relies on a uniform distribution to generate the coordinates of a new molecule. This is shown by comparing the calculation of an isotherm for a simple molecule adsorbed on a graphite substrate from a traditional GCMC simulation with that calculated using this new technique.     

The Detailed Balance Energy-scaled Displacement Monte Carlo Algorithm

Abstract

The Detailed Balance Energy-scaled Displacement Monte Carlo method that stems from the previously published Energy Scaled Displacement Monte Carlo method is presented. The results of tests performed on a dense Lennard-Jones liquid and on two particles in one dimension are reported.     

Test of the Inverse Monte Carlo Method for the Calculation of Interatomic Potential Energies in Atomic Liquids

Abstract

The principle purpose of this paper is to demonstrate the use of the Inverse Monte Carlo technique for calculating pair interaction energies in monoatomic liquids from a given equilibrium property. This method is based on the mathematical relation between transition probability and pair potential given by the fundamental equation of the “importance sampling” Monte Carlo method. In order to have well defined conditions for the test of the Inverse Monte Carlo method a Metropolis Monte Carlo simulation of a Lennard Jones liquid is carried out to give the equilibrium pair correlation function determined by the assumed potential. Because an equilibrium configuration is prerequisite for an Inverse Monte Carlo simulation a model system is generated reproducing the pair correlation function, which has been calculated by the Metropolis Monte Carlo simulation and therefore representing the system in thermal equilibrium. This configuration is used to simulate virtual atom displacements. The resulting changes in atom distribution for each single simulation step are inserted in a set of non-linear equations defining the transition probability for the virtual change of configuration. The solution of the set of equations for pair interaction energies yields the Lennard Jones potential by which the equilibrium configuration has been determined.     

On computational efficiency of the hybrid Monte Carlo method applied to the multicanonical ensemble

Abstract

In the present paper, computational efficiency of the hybrid Monte Carlo (HMC) method applied to the multicanonical ensemble is studied; the HMC is an equation of motion guided Monte Carlo method. As in the standard HMC for the canonical ensemble, the multicanonical HMC calculations with high acceptance ratio show better efficiency; about 60% acceptance yields the best performance for the system examined.     

Gibbs-Ensemble Molecular Dynamics: Liquid-Gas Equilibria for Lennard-Jones Spheres and n-Hexane

Abstract

We present a novel method to simulate phase equilibria in atomic and molecular systems. The method is a Molecular Dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles (or molecules) in a thermodynamically consistent fashion. Here we pres the derivation of the generalized equations of motion and discuss the relation of the resulting trajectory averages to the relevant ensemble. We test this Gibbs-Ensemble Molecular Dynamics algorithm by applying it to an atomic and a molecular system, i.e. to the liquid-gas coexistence in a Lennard-Jones fluid and in n-hexane. In both cases our results are in good accord with previous mean field and Gibbs-Ensemble Monte Carlo results as well as with the experimental data in the case of hexane. We also show that our Gibbs-Ensemble Molecular Dynamics algorithm like other Molecular Dynamics techniques can be used to study the dynamics of the system. Self-diffusion coefficients calculated with this method are in agreement with the result of conventional constant temperature Molecular Dynamics.     

A Modified Valley Restrained Monte Carlo Method to Efficiently Search the Low Energy Structures of Peptides

Abstract

A new Monte Carlo sampling scheme, namely the Modified Valley Restrained Monte Carlo procedure, is used to obtain the global energy minimum conformations for polypeptides, such as Met-enkephalin and Melittin. For each peptide, we found close agreement with previous results from both theoretical and experimental studies. The simple idea for controlling the step size according to the Valley Function, provides useful suggestions in searching the global energy minimum structures, and furthermore helps solve the multiple minima problem.     

Monte Carlo simulations in the isothermal-isobaric ensemble: use of a ‘shell’ particle for simulating polyatomic fluids

The recent reformulation of the isothermal-isobaric ensemble requires the use of a ‘shell’ particle to define uniquely the volume of the system, thereby avoiding the redundant counting of configurations. A previous modification of the Monte Carlo method, in which trial moves are generated and accepted consistent with the correct constant pressure partition function, is extended here to the case of polyatomic fluids. With a ‘shell’ molecule, either the centre of mass of the molecule or the location of any one of the atoms within the molecule can be chosen to define the system volume. Ensemble averages obtained with the use of the shell molecule differ from ensemble averages determined with the old (i.e. no shell particle) Monte Carlo algorithm, specifically for small system sizes, although both sets of averages become equal, as they must, in the thermodynamic limit. Monte Carlo simulations in the constant pressure ensemble for various Lennard-Jones polyatomic fluids, both for pure component and binary mixtures, demonstrate these differences for small systems. For mixtures, Monte Carlo simulations may include attempted identity swaps for the shell molecule, as the choice of which component serves as the shell molecule is arbitrary when periodic boundary conditions are applied.     

The Lennard-Jones Fluid Revisited: Computer Simulation Results

Abstract

Pseudoexperimental data of high accuracy on the pressure and the internal energy of the Lennard-Jones fluid have been generated both by the Monte Carlo and molecular dynamics methods for five subcritical and three supercritical isotherms. Values of the chemical potential of the Lennard-Jones fluid computed by a new version of the gradual insertion particle method for two isotherms up to very high densities are also reported and discussed, and compared with existing data.     

Simulations on the Primitive Electrolyte Environment of a High Charge-density Polyelectrolyte. A Sampling Problem and its Solution

Abstract

We show that the classical Metropolis Monte Carlo (MMC) algorithm converges very slowly when applied to the primitive electrolyte environment for a high charge-density polyelectrolyte. This slowness of convergence, which is due to the large density inhomogeneity around the polyelectrolyte, produces noticeable errors in the ion distribution functions for MMC runs of 1.3 × 10⁶ trial steps started from nonequilibrium distributions. We report that an algorithm which we call DSMC (for density-scaled Monte Carlo) overcomes this problem and provides relatively rapid convergence in this application. We suggest that DSMC should be well-suited for other Monte Carlo simulations on physical systems where large density inhomogeneities occur.     

Thermal Properties of Supercritical Carbon Dioxide by Monte Carlo Simulations

We present simulation results for the volume expansivity, isothermal compressibility, isobaric heat capacity, Joule-Thomson coefficient and speed of sound for carbon dioxide (CO 2 ) in the supercritical region, using the fluctuation method based on Monte Carlo simulations in the isothermal-isobaric ensemble. We model CO 2 as a quadrupolar two-center Lennard-Jones fluid with potential parameters reported in the literature, derived from vapor-liquid equilibria (VLE) of CO 2 . We compare simulation results with an equation of state (EOS) for the two-center Lennard-Jones plus point quadrupole (2CLJQ) fluid and with a multiparametric EOS adjusted to represent CO 2 experimental data. It is concluded that the VLE-based parameters used to model CO 2 as a quadrupolar two-center Lennard-Jones fluid (both simulations and EOS) can be used with confidence for the prediction of thermodynamic properties, including those of industrial interest such as the speed of sound or Joule-Thomson coefficient, for CO 2 in the supercritical region, except in the extended critical region.     

Application of Enhanced Sampling Monte Carlo Methods for High-Resolution Protein-Protein Docking in Rosetta

The high-resolution refinement of docked protein-protein complexes can provide valuable structural and mechanistic insight into protein complex formation complementing experiment. Monte Carlo (MC) based approaches are frequently applied to sample putative interaction geometries of proteins including also possible conformational changes of the binding partners. In order to explore efficiency improvements of the MC sampling, several enhanced sampling techniques, including temperature or Hamiltonian replica exchange and well-tempered ensemble approaches, have been combined with the MC method and were evaluated on 20 protein complexes using unbound partner structures. The well-tempered ensemble method combined with a 2-dimensional temperature and Hamiltonian replica exchange scheme (WTE-H-REMC) was identified as the most efficient search strategy. Comparison with prolonged MC searches indicates that the WTE-H-REMC approach requires approximately 5 times fewer MC steps to identify near native docking geometries compared to conventional MC searches.     

Direct Evaluation of Vapour-Liquid Equilibria of Mixtures by Molecular Dynamics Using Gibbs-Duhem Integration

Abstract

We present an extension of the Gibbs-Duhem integration method that permits direct evaluation of vapour-liquid equilibria of mixtures by molecular dynamics. The Gibbs-Duhem integration combines the best elements of the Gibbs ensemble Monte Carlo technique and thermodynamic integration. Given conditions of coexistence of pure substances, simultaneous but independent molecular dynamics simulations of each phase at constant number of particles, constant pressure, constant temperature and constant fugacity fraction of species 2 are carried out in succession along coexistence lines. In each simulation, the coexistence pressure is adjusted to satisfy the Clapeyron-type equation. The Clapeyron-type equation is a first-order nonlinear differential equation that prescribes how the pressure must change with the fugacity fraction of species 2 to maintain coexistence at constant temperature. The Clapeyron-type equation is solved by the predictor-corrector method. Running averages of mole fraction and compressibility factor for the two phases are used to evaluate the right-hand side of the Clapeyron-type equation. The Gibbs-Duhem integration method is applied to three prototypes of binary mixtures of the two-centre Lennard-Jones fluid having various elongations. The starting points on the coexistence curve were taken from published data.     

Monte Carlo Simulation Study of Adsorption Characteristics in Slit-Like Micropores Under Supercritical Conditions

Abstract

Adsorption characteristics of a solute diluted in supercritical fluids has been investigated by using the Monte Carlo simulation techniques. The Lennard-Jones potential function is used for describing interactions for a model system of CO₂ + benzene in slit-like micropores with infinite graphitic carbon walls. A modified μVT ensemble method with particle exchange proposed by Cracknell, Nicholson and Quirke (1993) is found to be much superior to the conventional μVT ensemble method especially for dense mixtures in a pore. Adsorption isotherms of CO₂ and benzene, in equilibrium with a dilute benzene mixture in CO₂ (mole fraction of benzene = 0.001), are computed by varying pressure, temperature, the benzene–surface interaction potential, and the slitwidth. Adsorption isotherm curve of CO₂ increases with an increase in pressure while that of benzene shows a maximum at a pressure far below the critical pressure of CO₂ and then it decreases with increasing pressure. The decrease in benzene adsorption with increasing pressure is attributable to both the enhanced solubility in supercritical CO₂ and the competitive adsorption of CO₂. The isotherm curves of each component at two temperatures, 313.2 K and 323.2 K, show to cross at a pressure near the critical pressure due to the “density effect” on the chemical potentials of a solute at supercritical fluid conditions. When the interaction between a solute and a surface increases, the adsorption isotherm increases. Narrowing the slitwidth results in the increase in the adsorption of solute since the external potential from two walls becomes deeper.     

Transition Path Sampling Study of the Local Molecular Structure in the Aqueous Solvation of Sodium Chloride

Abstract

The aqueous solvation of sodium chloride has been investigated using the recently introduced technique of the transition path sampling. We performed a series of Monte Carlo simulations for each element of an ensemble of chains of states. The evolution of the local solvent structure during the dissociation process has been observed. The incoming of a couple of waters to the first coordination shell is responsible for the structural changes which allow the dissociation occur: waters which leave the second coordination shell produce voids and a local molecular reorganization in order to allocate the dissociated ion pair.     

A Modified Gibbs Ensemble Method for Calculating Fluid Phase Equilibria

Abstract

A modification of the Gibbs ensemble Monte Carlo computer simulation method for fluid phase equilibrium is described. The modification, which is based on the assumption of a thermodynamic model for the vapor phase, reduces the computational time for the simulation as compared to the original Gibbs ensemble methods. Since the computational time is largely proportional to the number of particle-particle interactions, avoiding the direct simulation of the vapor phase typically leads to a thirty to forty percent reduction in computational time. For a pure Leonard-Jones-(12,6) fluid the results obtained at moderate reduced temperatures, T/T_c < 0.8, are in good agreement with the full Gibbs ensemble method.     

Local Order in Some Aprotic Dipolar Liquids

Abstract

Analysis of the local order in four aprotic, dipolar liquids (acetonitrile, acetone, N,N- dimethylformamide, and pyridine) has been carried out by computer simulation methods. The effect of the dipole-dipole-as well as of the Lennard-Jones interactions on the pair distribution functions and some orientational characteristics are discussed. A general observation for this group of liquids is that the local order is accompanied by a favorable orientation of the nearest neighbor molecules with dipole vectors anti-parallel to the central one. This orientation decays rapidly and extends to slightly more than to the first neighbor molecules only. Differences in local order due to the differences in shapes and symmetries of the molecules are discussed.

The present study partly employs the recently developed Reverse Monte Carlo method, special features of which in contrast to the traditional simulation methods are also discussed.     

Some Recent Developments in Computational Chemistry

Abstract

Some recent developments in the use of computational methods to predict the properties of condensed phases are discussed; the use of Gibbs ensemble Monte Carlo to predict the phase equilibria of bulk phases, the use of molecular dynamics to elucidate Atomic Force Microscopy experiments on organic films, and the use of combined Monte Carlo/molecular dynamics techniques to enable the direct prediction of particle fluxes along pressure gradients in model microporous materials.