Thermodynamic properties of van der Waals fluids from Monte Carlo simulations and perturbative Monte Carlo theory

Computer simulations have been performed for fluids with van der Waals potential, that is, hard spheres with attractive inverse power tails, to determine the equation of state and the excess energy. On the other hand, the first- and second-order perturbative contributions to the energy and the zero- and first-order perturbative contributions to the compressibility factor have been determined too from Monte Carlo simulations performed on the reference hard-sphere system. The aim was to test the reliability of this "exact" perturbation theory. It has been found that the results obtained from the Monte Carlo perturbation theory for these two thermodynamic properties agree well with the direct Monte Carlo simulations. Moreover, it has been found that results from the Barker-Henderson [J. Chem. Phys. 47, 2856 (1967)] perturbation theory are in good agreement with those from the exact perturbation theory.     

Molecular organization in liquids consisting of anisometric molecules: Monte Carlo simulations and comparison with theories

Monte Carlo computer simulations of systems of rigid anisometric molecules with axial ratio x 3, 4 and 5 have been performed at densities comparable to those of the mesogenic groups in liquid crystalline materials. The agreement between computer experiments and Flory-Ronca-Irvine theory is quite good. In particular, the nematic-isotropic transition temperatures found in the simulations are within 20 K of those predicted by the theory in all cases examined with x 3 and 4, and only slightly higher when x 5. Furthermore, the orientational distributions of the molecular long axes with respect to the nematic director are always nearly coincident with those evaluated according to the FRI theory, provided that the order parameter is the same. However, the order parameters at the transition points are smaller than expected on the basis of the theory, and closer to experimental values. Calculations performed by changing the size of the simulated systems prove that the reported results are size independent, while size dependent results are obtained with smaller systems.     

Monte Carlo simulations of thermodynamic and structural properties of Mie(14,7) fluids

The vapor-liquid phase envelope of Mie(14,7) fluids is determined by the Gibbs ensemble Monte Carlo (MC) simulation technique. The NVT-MC simulation method is then utilized to compute the equation of state and the pair correlation function over a wide range of densities and temperatures. The effective diameters are calculated via the virial minimization method and the results are applied as the repulsion-attraction splitting distance within the generic van der Waals (GvdW) theory to compute the mean free volume. The density and temperature dependence of these parameters are studied and discussed. The results for the effective diameter, and the GvdW parameters are fitted to analytical functions of density and temperature. An examination of the results for the fluid phase equilibria of argon shows excellent agreement with empirical data for the densities of the coexisting phases, the vapor pressure, and the critical point. The computed free volumes are used to compute the diffusion coefficient of argon and the results are compared with experimental data.     

Surface tension of associating fluids by Monte Carlo simulations

Canonical Monte Carlo (NVT-MC) simulations were performed to obtain surface tension and coexistence densities at the liquid-vapor interface of one-site associating Lennard-Jones and hard-core Yukawa fluids, as functions of association strength and temperature. The method to obtain the components of the pressure tensor from NVT-MC simulations was validated by comparing the equation of state of the associative hard sphere system with that coming from isothermal-isobaric Monte Carlo simulations. Surface tension of the associative Lennard-Jones fluid determined from NVT-MC is compared with previously reported results obtained by molecular dynamics simulations of a pseudomixture model of monomers and dimers. A good agreement was found between both methods. Values of surface tension of associative hard-core Yukawa fluids are presented here for the first time.     

Optimized ensemble Monte Carlo simulations of dense Lennard-Jones fluids

We apply the recently developed adaptive ensemble optimization technique to simulate dense Lennard-Jones fluids and a particle-solvent model by broad-histogram Monte Carlo techniques. Equilibration of the simulated fluid is improved by sampling an optimized histogram in radial coordinates that shifts statistical weight towards the entropic barriers between the shells of the liquid. Interstitial states in the vicinity of these barriers are identified with unprecedented accuracy by sharp signatures in the quickly converging histogram and measurements of the local diffusivity. The radial distribution function and potential of mean force are calculated to high precision.     

Heat capacities of Stockmayer fluids from Monte Carlo simulations and perturbation theory

The heat capacities of dipolar fluids are investigated using a thermodynamic perturbation theory approach and the NVT and NpT Monte Carlo simulation methods. The theoretical results are compared to corresponding simulation data. The comparison shows that the applied perturbation theory is appropriate for the heat capacity calculations. As an application, the isobaric heat capacity of ammonia is also studied by the Stockmayer fluid model.     

Inhomogeneous fluids of colloidal hard dumbbells: fundamental measure theory and Monte Carlo simulations

Recently, a density functional theory for hard particles with shape anisotropy was developed, the extended deconvolution fundamental measure theory (edFMT). We apply edFMT to hard dumbbells, arguably the simplest non-convex shape and readily available experimentally in the form of colloids. We obtain good agreement between edFMT and Monte Carlo simulations for fluids of dumbbells in a slit and for the same system under gravity. This indicates that edFMT can be successfully applied to nearly all colloidal shapes, not just for the convex shapes for which edFMT was originally derived. A theory, such as edFMT, that allows a fast and general way of mapping the phase behavior of anisotropic colloids, can act as a useful guide for the design of colloidal shapes for various applications.     

Ion binding to polyelectrolytes: Monte Carlo simulations versus classical mean field theories

The influence of ion size and surface charge model in titrations of ionizable polyelectrolytes is studied by means of the Semi Grand Canonical Monte Carlo simulation method in the context of the primitive model. Three models describing a discrete distribution of charged functional groups on the polyelectrolyte and different values for the radius of the background electrolyte spanning from ionic to hydrated radii values were analyzed. The polyelectrolyte titrations were simulated by calculating the degree of ionization versus pH curves at two ionic strengths. The results allow us to quantify the impact of the sizes of the background salt ions and surface functional groups of the polyelectrolyte on the dissociation degree. This influence is explained in terms of the effectiveness of the screening of the charged surface sites. Finally, by comparison with the Non-Linear Poisson–Boltzmann model, the influence of ionic correlations and finite size of the solution ions is assessed.     

The use of Monte Carlo simulations to evaluate kinetic data and analytic approximations

Experimental kineticists are always faced with the problem of reducing kinetic data to extract physically meaningful information. A particularly vexing problem arises when different models reproduce the data but yield different values for the physical parameters. For over forty-five years Monte Carlo simulation techniques have been used to study the statistical behavior of parameters extracted from data. Not only do these simulations provide realistic uncertainties, correlation coefficients, and confidence envelopes, but they also provide insight into the nature of the model. These insights may be obtained by viewing two-dimensional scatter plots of the fractional changes of the parameters and one-dimensional histograms of the distributions of the changes in the parameters. Monte Carlo simulations are illustrated with examples from OH+CH₄ → CH₃+H₂O and the high-pressure rate coefficient for methyl-methyl association. A more complex problem involves models for pressure-dependent rate coefficients in the falloff region. We have modeled methyl-methyl association with five of the most current analytic approximations for behavior in the falloff region. All of these reproduce the data to within their uncertainties. However, when Monte Carlo techniques are applied the correlations between the parameters and the nonlinear nature of their behavior become evident. We postulate that the statistical behavior of the parameters of a model may be used to distinguish one model from another and, thereby, identify those analytic approximations that hold promise for further investigation and utilization. Finally, the recent advent of high-speed workstations implies that Monte Carlo simulations should become a routine part of the analysis of kinetic data. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 803–817, 1997     

Monte Carlo study of interfacial properties of associating fluids

Canonical Monte Carlo Simulations have been performed to calculate liquid-vapor properties of the associating square well and Lennard-Jones fluids with one and two sites. Simulations were carried out by using several values of reduced temperatures and association energies. The orthobaric densities, as well as the surface tension of associating square well fluids, were calculated and compared with those reported previously in literature; a good agreement was found among them. Results of surface tension of two-sites associating Lennard-Jones fluids are presented here for the first time.     

Monte Carlo simulations of fluids whose particles interact with a logarithmic potential

Monte Carlo simulations of a model fluid in which the particles interact via a continuous potential that has a logarithmic divergence at a pair separation of sigma, which we introduced in J. G. Powles et al., Proc. R. Soc. London, Ser. A 455, 3725 (1999), have been carried out. The potential has the form, phi(r)= -epsilon ln(fr), where epsilon sets the energy scale and fr=1-(sigma/r)m. The value of m chosen was 12 but the qualitative trends depend only weakly on the value of m, providing it is greater than 3. The potential is entirely repulsive and has a logarithmic divergence as approximately -ln(r/sigma-1) in the r-->sigma limit. Predictions of the previous paper that the internal energy can be computed at all temperatures using the standard statistical mechanics formula for continuous potentials are verified here. The pressure can be calculated using the usual virial expression for continuous potentials, although there are practical limitations in resolving the increasingly important contribution from the r-->sigma limit at reduced temperatures greater than approximately 5. The mean square force F2 and infinite frequency shear Ginfinity and bulk Kinfinity moduli are only finite for T*=kBT/epsilon<1. The logarithmic fluid's physical properties become increasingly more like that of the hard sphere fluid with increasing temperature, showing a sharp transition in the behavior of the mean square force and infinite frequency elastic constants at T*=1. The logarithmic fluid is shown to exhibit a solid-fluid phase transition.     

Mercedes-Benz water molecules near hydrophobic wall: integral equation theories vs Monte Carlo simulations

Associative version of Henderson-Abraham-Barker theory is applied for the study of Mercedes-Benz model of water near hydrophobic surface. We calculated density profiles and adsorption coefficients using Percus-Yevick and soft mean spherical associative approximations. The results are compared with Monte Carlo simulation data. It is shown that at higher temperatures both approximations satisfactory reproduce the simulation data. For lower temperatures, soft mean spherical approximation gives good agreement at low and at high densities while in at mid range densities, the prediction is only qualitative. The formation of a depletion layer between water and hydrophobic surface was also demonstrated and studied.     

Monte Carlo simulations of liquid crystals between microstructured substrates

The structure of a model liquid crystalline fluid confined between two microstructured substrates is studied through Monte Carlo simulations. A simple model for a structured substrate, similar in spirit to those used for rough walls of walls with grafted polymers, is introduced. It is found that varying the structure of the substrate, a transition in the alignment of the confined fluid, from parallel to perpendicular, is induced. For particular substrate structures, it is possible to induce tilted alignment in the confined fluid, the tilt angle being temperature dependent.     

Efficient multiparticle sampling in Monte Carlo simulations on fluids: application to polarizable models

A novel Monte Carlo simulation scheme based on biased simultaneous displacements of all particles of the system has been developed. The method is particularly suited for systems with nonadditive interactions and its efficiency is demonstrated by its implementation for the polarizable Stockmayer fluid. Performance of the method is compared with both the standard one-particle move method and an unbiased multiparticle scheme by computing the mean squared displacements, rotation relaxation, and the speed of equilibration (translational order parameter). It is shown that the proposed biased method is about a factor of 10 faster, for the system considered, when compared with the other schemes.     

Elastic properties of soft sphere crystal from Monte Carlo simulations

Elastic properties of faced centered cubic (fcc) crystals composed of soft spheres, interacting through potentials of the form u(r) ~ r(-n), have been investigated by Monte Carlo (MC) simulations. It is shown that both the softness parameter (n(-1)) and temperature strongly influence the elastic properties of the studied system. The simulations show explicitly that when T > 0 the elastic constants of the hard sphere crystal can be obtained by taking the limit n --> infinity of soft spheres. When T --> 0 for any finite n, the elastic constants of the soft spheres tend to those of the static model. At all temperatures and softness parameters studied here, n, the Poisson's ratio in [110] (perpendicular direction) is negative.     

Generalizations for the potential of mean force between two isolated colloidal particles from Monte Carlo simulations

A substantial amount of experimental and numerical evidence has shown that the Derjaguin-Landau-Verwey-Overbeek theory is not suitable for describing those colloidal solutions that contain multivalent counterions. Toward improved understanding of such solutions, the authors report Monte Carlo calculations wherein, following Rouzina and Bloomfield, they postulate that, in the absence of van der Waals forces, the overall force between two isolated charged colloidal particles in electrolyte solutions is determined by a dimensionless parameter Gamma=z(2)l(B)/a, which measures the electrostatic repulsion between counterions adsorbed on the macroion surface, where z = counterion valence, l(B)=Bjerrum length, and a = average separation between counterions on the macroion surface calculated as if the macroion were fully neutralized. The authors find, first, that the maximum repulsion between like-charged macroions occurs at Gamma approximately 0.5 and, second, that onset of attraction occurs at Gamma approximately 1.8, essentially independent of the valence and concentration of the surrounding electrolyte. These observations might provide new understanding of interactions between electrostatic double layers and perhaps offer explanations for some electrostatic phenomena related to interactions between DNA molecules or proteins.     

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.     

Monte Carlo simulations of the polymer glass transition: From the test of theories to material modeling

We present results on the glass transition in polymer melts using Monte Carlo simulations of the bond fluctuation lattice model. There are two questions we address in this work. What is the temperature dependence of the entropy density in such a model polymer melt and how well is it described by theories like the Gibbs-DiMarzio theory of the glass transition? And to what degree is one able to map the Hamiltonian of such an abstract lattice model onto a specific polymer material and use it to model the large scale and long time properties of a realistic polymer melt?