Modeling desorption of fluids from disordered mesoporous materials

The desorption mechanism of fluids in disordered mesoporous glasses is studied by Monte Carlo simulations of a coarse-grained lattice model with realistic matrix configurations representative of Vycor. Two methods of simulation are considered: grand canonical ensemble Monte Carlo simulations and dynamic Monte Carlo simulations which mimic the diffusion of the fluid in and out of the material using Kawasaki dynamics. In the grand canonical simulations, cavitation via nucleation of bubbles inside the pores plays the dominant role in determining the fluid configurations along the desorption isotherm. The Kawasaki dynamics simulations indicate that such configurations are achieved dynamically via the gradual advancement of macroscopic front interfaces toward the interior. This is made possible by the bubble nucleation mechanism operating on a length scale that is determined by both the typical pore size and the strength of the solid-fluid interaction.     

Reactive Monte Carlo and grand-canonical Monte Carlo simulations of the propene metathesis reaction system

The influence of silicalite-1 pores on the reaction equilibria and the selectivity of the propene metathesis reaction system in the temperature range between 300 and 600 K and the pressure range from 0.5 to 7 bars has been investigated with molecular simulations. The reactive Monte Carlo (RxMC) technique was applied for bulk-phase simulations in the isobaric-isothermal ensemble and for two phase systems in the Gibbs ensemble. Additionally, Monte Carlo simulations in the grand-canonical ensemble (GCMC) have been carried out with and without using the RxMC technique. The various simulation procedures were combined with the configurational-bias Monte Carlo approach. It was found that the GCMC simulations are superior to the Gibbs ensemble simulations for reactions where the bulk-phase equilibrium can be calculated in advance and does not have to be simulated simultaneously with the molecules inside the pore. The confined environment can increase the conversion significantly. A large change in selectivity between the bulk phase and the pore phase is observed. Pressure and temperature have strong influences on both conversion and selectivity. At low pressure and temperature both conversion and selectivity have the highest values. The effect of confinement decreases as the temperature increases.     

Path integral Monte Carlo method for ab initio calculation

The Feynman path integral method is applied to the many-electron problem. We first give new closure relations in terms of ordinary complex and real numbers, which could be derived from an arbitrary complete set of state vectors. Then, in the path integral form, the partition function of the system and the ensemble average of energy are explicitly expressed in terms of these closure relations. It is impossible to evaluate the path integral by direct numerical integrations because of its huge amount of integration variables. Therefore, we develop an algorithm by the Monte Carlo method with constraints corresponding to the normalization condition of states to calculate the required integral. Finally, the ensemble average of energy for the hydrogen molecule is explicitly evaluated by the quantum Monte Carlo method and results are compared with the result obtained by the ordinary full configuration interaction (CI) method. © 1996 John Wiley & Sons, Inc.     

A rapidly convergent simulation method: Mixed Monte Carlo/stochastic dynamics

Although Monte Carlo and molecular dynamics are the primary methods used for free energy simulations of molecular systems, their application to molecules that have multiple conformations separated by energy barriers of ≥ 3 kcal/mol is problematic because of slow rates of convergence. In this article we introduce a hybrid simulation method termed MC-SD which mixes Monte Carlo (MC) and stochastic dynamics (SD). This new method generates a canonical ensemble via alternating MC and SD steps and combines the local exploration strengths of dynamics with the barrier-crossing ability of large-step Monte Carlo. Using calculations on double-well potentials and long simulations (10⁸ steps of MC and 1 μs of SD) of the simple, conformationally flexible molecule n-pentane, we find that MC-SD simulations converage faster than either MC or SD alone and generate ensembles which are equivalent to those created by classical MC or SD. Using pure SD at 300 K, the conformational populations of n-pentane are shown to be poorly converged even after a full microsecond of simulation. © 1994 by John Wiley & Sons, Inc.     

Computational studies on thermodynamic properties, effective diameters, and free volume of argon using an ab initio potential

A quantum mechanical derived ab initio interaction potential for the argon dimer was tested in molecular simulations to reproduce the thermophysical properties of the vapor-liquid phase equilibria using the Gibbs ensemble Monte Carlo simulations as well as the liquid and supercritical equation of state using the NVT Monte Carlo simulations. The ab initio interaction potential was taken from the literature. A recently developed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was used to compute the effective diameters of argon in fluid phases and the results were subsequently applied in the generic van der Waals theory to compute the free volume of argon. The calculated densities of the coexisting phases, the vapor pressure, and the equation of state show excellent agreement with experimental values. The effective diameters and free volumes of argon are given over a wide range of densities and temperatures. An empirical formula was used to fit the effective diameters as a function of density and temperature. The computed free volume will be used in future investigations to calculate the transport properties of argon.     

Monte Carlo simulations of solid nitrogen in the isothermal–isobaric ensemble with an ab initioSCF–CI potential

Monte Carlo simulations in the isothermal–isobaric ensemble for the α phase of solid N₂ have been carried out with two different pair potentials obtained from ab initio quantum chemical calculations. Comparison is made with data obtained from empirical potentials as well as with experimental data, and in general, the agreement with experiment is found to be good. It is also found that the differences between experiment and theory can largely be attributed to correlation effects, not considered in the quantum mechanical calculations of the pair potential.     

Rigid quantum Monte Carlo simulations of condensed molecular matter: water clusters in the n=2-->8 range

The numerical advantage of quantum Monte Carlo simulations of rigid bodies relative to the flexible simulations is investigated for some simple systems. The results show that if high frequency modes in molecular condensed matter are predominantly in the ground state, the convergence of path integral simulations becomes nonuniform. Rigid body quantum parallel tempering simulations are necessary to accurately capture thermodynamic phenomena in the temperature range where the dynamics are influenced by intermolecular degrees of freedom; the stereographic projection path integral adapted for quantum simulations of asymmetric tops is a significantly more efficient strategy compared with Cartesian coordinate simulations for molecular condensed matter under these conditions. The reweighted random series approach for stereographic path integral Monte Carlo is refined and implemented for the quantum simulation of water clusters treated as an assembly of rigid asymmetric tops.     

Prediction of the thermophysical properties of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton by Monte Carlo simulation using ab initio potentials

Gibbs ensemble Monte Carlo simulations were used to test the ability of intermolecular pair potentials derived ab initio from quantum mechanical principles, enhanced by Axilrod-Teller triple-dipole interactions, to predict the vapor-liquid phase equilibria of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton. The interaction potentials for Ne-Ne, Ar-Ar, Kr-Kr, and Ne-Ar were taken from literature; for Ar-Kr a different potential has been developed. In all cases the quantum mechanical calculations had been carried out with the coupled-cluster approach [CCSD(T) level of theory] and with correlation consistent basis sets; furthermore an extrapolation scheme had been applied to obtain the basis set limit of the interaction energies. The ab initio pair potentials as well as the thermodynamic data based on them are found to be in excellent agreement with experimental data; the only exception is neon. It is shown, however, that in this case the deviations can be quantitatively explained by quantum effects. The interaction potentials that have been developed permit quantitative predictions of high-pressure phase equilibria of noble-gas mixtures.     

Statistical thermodynamics of internal rotation in a hindering potential of mean force obtained from computer simulations

A method of statistical estimation is applied to the problem of one-dimensional internal rotation in a hindering potential of mean force. The hindering potential, which may have a completely general shape, is expanded in a Fourier series, the coefficients of which are estimated by fitting an appropriate statistical-mechanical distribution to the random variable of internal rotation angle. The function of reduced moment of inertia of an internal rotation is averaged over the thermodynamic ensemble of atomic configurations of the molecule obtained in stochastic simulations. When quantum effects are not important, an accurate estimate of the absolute internal rotation entropy of a molecule with a single rotatable bond is obtained. When there is more than one rotatable bond, the "marginal" statistical-mechanical properties corresponding to a given internal rotational degree of freedom are reduced. The method is illustrated using Monte Carlo simulations of two public health relevant halocarbon molecules, each having a single internal-rotation degree of freedom, and a molecular dynamics simulation of an immunologically relevant polypeptide, in which several dihedral angles are analyzed.     

Thermodynamics and partitioning of homopolymers into a slit-A grand canonical Monte Carlo simulation study

Grand canonical ensemble Monte Carlo simulation (GCMC) combined with the histogram reweighting technique was used to study the thermodynamic equilibrium of a homopolymer solution between a bulk and a slit pore. GCMC gives the partition coefficients that agree with those from canonical ensemble Monte Carlo simulations in a twin box, and it also gives results that are not accessible through the regular canonical ensemble simulation such as the osmotic pressure of the solution. In a bulk polymer solution, the calculated osmotic pressure agrees very well with the scaling theory predictions both for the athermal polymer solution and the theta solution. However, one cannot obtain the osmotic pressure of the confined solution in the same way since the osmotic pressure of the confined solution is anisotropic. The chemical potentials in GCMC simulations were found to differ by a translational term from the chemical potentials obtained from canonical ensemble Monte Carlo simulations with the chain insertion method. This confirms the equilibrium condition of a polymer solution partition between the bulk and a slit pore: the chemical potentials of the polymer chain including the translational term are equal at equilibrium. The histogram reweighting method enables us to obtain the partition coefficients in the whole range of concentrations based on a limited set of simulations. Those predicted bulk-pore partition coefficient data enable us to perform further theoretical analysis. Scaling predictions of the partition coefficient at different regimes were given and were confirmed by the simulation data.     

Computational study of the structures of gaseous helium-3 at low temperature

A computational study of gaseous helium-3 at T=5.23 K, for number densities rho N<0.0021 A(-3), analyzing the different pair and triplet structures in both the r and the k spaces, is presented. Structures in r space (i.e., instantaneous, total continuous linear response, and centroids) are determined via path-integral Monte Carlo simulations in the canonical ensemble by utilizing the Aziz-Slaman and the SAPT2 interatomic potentials. Additional results obtained with the application of two closures for triplets in r space, the Kirkwood superposition approximation and the Jackson-Feenberg convolution, are also reported. Besides, an analysis of the nonsuitability of quantum hard spheres for describing this system is included. The pair structures in k space are fixed via Ornstein-Zernike schemes appropriate for dealing with quantum diffraction effects in fluids. The effect on the quality of the computed isothermal compressibilities brought about by increasing the sample size in the simulations and by the subsequent application of a grand ensemble correction to the asymptotic behavior of the canonical pair radial correlation functions is also investigated. Furthermore, it is demonstrated analytically that the methods of classical statistical mechanics for dealing with the higher-order direct correlation functions remain fully valid for studying the higher-order correlations of path-integral centroids. By taking advantage of this result, the triplet structure factors for the centroid (also for the instantaneous) correlations are computed by following a number of distinct closures and methods that involve triplet direct correlation functions. The latter computations are intended to explore an alternative scheme to the very expensive fixing of triplet structure factors through direct path integral simulations, an alternative which is expected to yield the main features of these triplet quantities for this gas. Comparison with experiment is made wherever possible, and the results presented allow one to explain the substantial structural features existing in gaseous helium-3.     

Monte Carlo simulation of polarizable systems: Early rejection scheme for improving the performance of adiabatic nuclear and electronic sampling Monte Carlo simulations

An early rejection scheme for trial moves in adiabatic nuclear and electronic sampling Monte Carlo simulation (ANES-MC) of polarizable intermolecular potential models is presented. The proposed algorithm is based on Swendsen–Wang filter functions for prediction of success or failure of trial moves in Monte Carlo simulations. The goal was to reduce the amount of calculations involved in ANES-MC electronic moves, by foreseeing the success of an attempt before making those moves. The new method was employed in Gibbs ensemble Monte Carlo (GEMC) simulations of the polarizable simple point charge-fluctuating charge (SPC-FQ) model of water. The overall improvement in GEMC depends on the number of swap attempts (transfer molecules between phases) in one Monte Carlo cycle. The proposed method allows this number to increase, enhancing the chemical potential equalization. For a system with 300 SPC-FQ water molecules, for example, the fractions of early rejected transfers were about 0.9998 and 0.9994 at 373 and 423 K, respectively. This means that the transfer moves consume only a very small part of the overall computing effort, making GEMC almost equivalent to a simulation in the canonical ensemble.     

DFT-LDA pseudopotentials in quantum Monte Carlo

We investigate the portability of standard norm-conserving pseudopotentials outside the density functional theory-local density approximation (DFT-LDA) framework, i.e., their use and interpretation as electron-ion effective potentials in valence-only diffusion Monte Carlo simulations. While first-principles many-body pseudopotentials are not available in the literature yet, the use of approximate pseudopotentials in quantum Monte Carlo simulations is becoming widespread. Here we attempt a systematic analysis of the portability of norm-conserving pseudopotentials generated within DFT-LDA, focusing on a model many-body system, the two-electron valence-only ion. Our results indicate that the portability is good in most cases, hence the use of pseudopotentials in quantum Monte Carlo simulations is in general a reasonable approximation but suggest that in some cases this approximation may be relevant. © 1997 John Wiley & Sons, Inc.     

Microscopic approaches to liquid nitromethane detonation properties

In this paper, thermodynamic and chemical properties of nitromethane are investigated using microscopic simulations. The Hugoniot curve of the inert explosive is computed using Monte Carlo simulations with a modified version of the adaptative Erpenbeck equation of state and a recently developed intermolecular potential. Molecular dynamic simulations of nitromethane decomposition have been performed using a reactive potential, allowing the calculation of kinetic rate constants and activation energies. Finally, the Crussard curve of detonation products as well as thermodynamic properties at the Chapman-Jouguet (CJ) point are computed using reactive ensemble Monte Carlo simulations. Results are in good agreement with both thermochemical calculations and experimental measurements.     

Phase behavior of elemental aluminum using monte carlo simulations

Monte Carlo simulations are presented for two models of aluminum: an embedded-atom model and an explicit many-body model. Vapor/liquid coexistence curves are determined using Gibbs ensemble Monte Carlo simulations. The normal boiling points predicted by both models are somewhat higher (by about 10%) than the experimental value. Isothermal constant-stress simulations are used to simulate solid Al from 300 K to the triple point. The solid structures are at least metastable in the face-centered cubic configuration, and the specific heat is determined to be lower than the experimental value. The melting point for the embedded-atom model determined via thermodynamic integration along a pseudo-supercritical path is approximately 20% higher than the experimental value.     

Multiparticle moves in acceptance rate optimized monte carlo

Molecular Dynamics (MD) and Monte Carlo (MC) based simulation methods are widely used to investigate molecular and nanoscale structures and processes. While the investigation of systems in MD simulations is limited by very small time steps, MC methods are often stifled by low acceptance rates for moves that significantly perturb the system. In many Metropolis MC methods with hard potentials, the acceptance rate drops exponentially with the number of uncorrelated, simultaneously proposed moves. In this work, we discuss a multiparticle Acceptance Rate Optimized Monte Carlo approach (AROMoCa) to construct collective moves with near unit acceptance probability, while preserving detailed balance even for large step sizes. After an illustration of the protocol, we demonstrate that AROMoCa significantly accelerates MC simulations in four model systems in comparison to standard MC methods. AROMoCa can be applied to all MC simulations where a gradient of the potential is available and can help to significantly speed up molecular simulations. © 2015 Wiley Periodicals, Inc.     

Spatially adaptive lattice coarse-grained Monte Carlo simulations for diffusion of interacting molecules

While lattice kinetic Monte Carlo (KMC) methods provide insight into numerous complex physical systems governed by interatomic interactions, they are limited to relatively short length and time scales. Recently introduced coarse-grained Monte Carlo (CGMC) simulations can reach much larger length and time scales at considerably lower computational cost. In this paper we extend the CGMC methods to spatially adaptive meshes for the case of surface diffusion (canonical ensemble). We introduce a systematic methodology to derive the transition probabilities for the coarse-grained diffusion process that ensure the correct dynamics and noise, give the correct continuum mesoscopic equations, and satisfy detailed balance. Substantial savings in CPU time are demonstrated compared to microscopic KMC while retaining high accuracy.