A Monte Carlo study of liquid benzene

A Monte Carlo study for liquid benzene, where the pair potential is described as a sum of Lennard-Jones (12, 6) interactions acting between six equivalent centres of each molecule, has been carried out. The potential parameters have been chosen in such a way that the internal energy and virial pressure were in agreement with the experimental values. The liquid structure has been examined in detail by computing several distribution functions: the radial, the angular correlations and the radial-angular distribution between molecular planes.     

A Monte Carlo study of water clusters

The free energy, internal energy and entropy for water clusters of 8 and 64 molecules have been calculated using the Monte Carlo method and both the Lie-Clementi (C-XII) and Stillinger (ST2) potential functions. Detailed structural data for a cluster of 64 water molecules was also obtained by averaging over a large number of configurations. Average dipole orientations and radial density distributions indicated that individual water molecules were not particularly restricted to fixed orientations and that general ordering tendencies were small.     

A hybrid fluctuating hydrodynamics and kinetic Monte Carlo method for modeling chemically-powered nanoscale motion

We present a stochastic multiscale method for modeling heterogeneous catalysis at the nanoscale. The system is decomposed into the fluid domain and the catalyst-fluid interface. We implemented the fluctuating hydrodynamics framework to model the diffusion of the chemical species in the fluid domain, and the chemical master equation to describe the catalytic activity at the interface. The coupling between the domains occurs simultaneously. Using a simple one-dimensional (1D) linear model, we showed that the predictions of our scheme are in excellent agreement with deterministic simulations. The method was specifically developed to model the spatially asymmetric catalysis on the surface of self-propelled nanoswimmers. Numerical simulations showed that our approach can estimate the uncertainty in the swimming velocity resulting from inherent stochastic nature of the chemical reactions at the catalytic interface. Although the method has been applied to simple 1D and 2D models, it can be generalized to handle different geometries and more sophisticated chemical reactions. Therefore, it can serve as a practical mathematical tool to study how the efficiency of chemically powered nanomachines is affected by the interplay between structural complexity, nonlinear reactivity, and nonequilibrium fluctuations.     

A Monte Carlo simulation study of branched polymers

Monte Carlo simulations are presented for the static properties of highly branched polymer molecules. The molecules consist of a semiflexible backbone of hard-sphere monomers with semiflexible side chains, also composed of hard-sphere monomers, attached to either every backbone bead or every other backbone bead. The conformational properties and structure factor of this model are investigated as a function of the stiffness of the backbone and side chains. The average conformations of the side chains are similar to self-avoiding random walks. The simulations show that there is a stiffening of the backbone as degree of crowding is increased, for example, if the branch spacing is decreased or side chain length is increased. The persistence length of the backbone is relatively insensitive to the stiffness of the side chains over the range investigated. The simulations reproduce most of the qualitative features of the structure factor observed in experiment, although the magnitude of the stiffening of the backbone is smaller than in experiment.     

Mode excitation Monte Carlo simulations of mesoscopically large membranes

Solvent-free coarse grained models represent one of the most promising approaches for molecular simulations of mesoscopically large membranes. In these models, the size of the simulated membrane is limited by the slow relaxation time of the longest bending mode. Here, we present a Monte Carlo algorithm with update moves in which all the lipids are simultaneously displaced. These collective moves result in fast excitation and relaxation of the long wavelength thermal fluctuations. We apply the method to simulations of a bilayer membrane with linear size of approximately 50 nm and show reduction in the relaxation time by two orders of magnitudes when compared to conventional Monte Carlo.     

Monte Carlo study of polyelectrolyte gels

The discontinuous transition between dense and dilute phases in polyelectrolyte gels is observed in Bond-Fluctuation Method Monte Carlo simulations of gels. The transition is driven by the competition between local attractive interactions of a poor-quality solvent and global repulsive interactions from counter-ion pressure. A procedure is introduced that prevents local attractive interactions from destroying ergodicity. Under good solvent conditions, lengths and volumes of gels are found to follow self-avoiding random walk scaling. © 1996 John Wiley & Sons, Inc.     

Quantum Monte Carlo study of first-row atoms using transcorrelated variational Monte Carlo trial functions

The effect of using the transcorrelated variational Monte Carlo (TC-VMC) approach to construct a trial function for fixed node diffusion Monte Carlo (DMC) energy calculations has been investigated for the first-row atoms, Li to Ne. The computed energies are compared with fixed node DMC energies obtained using trial functions constructed from Hartree-Fock and density functional levels of theory. Despite major VMC energy improvement with TC-VMC trial functions, no improvement in DMC energy was observed using these trial functions for the first-row atoms studied. The implications of these results on the nodes of the trial wave functions are discussed.     

Dynamics of grafted star-branched polymers. A Monte Carlo study

Monte Carlo simulations of simple models of star-branched polymers were carried out. The model chains were confined to simple cubic lattice and consisted of f = 3 branches of equal length and the total number of polymer segments as well as the density of grafted chains on the surface were varied. The chains have had one arm end attached to an impenetrable plate. The simulations were performed by employing the set of local micromodifications of the chain conformations. The model chains were athermal, i.e. good solvent conditions were modeled, the excluded volume effect was present at the model. The density of grafted chains on the surface was varied from a single chain up to 0.3. The static and dynamic properties of the system were studied. The influence of polymer concentration as well as the polymer length on static and dynamic properties of the system studied was shown. The relation between the structure and short-time dynamics (relaxation times) was discussed.     

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.     

Monte Carlo study of lipid membranes: Simulation of diparmitoylphosphatidylcholine bilayers in gel and liquid-crystalline phases

The statistical properties of the bilayer membranes of diparmitoylphosphatidylcholine (DPPC) in the gel and liquid-crystal phases were studied by Monte Carlo (MC) simulation using potential functions of the Lennard-Jones, the simple Coulomb, and the bond torsion. The simulation was undertaken on a two-dimensional periodic condition imposed on the bilayer model consisting of faithfully described molecules. The structure and ordering of the model bilayers accorded well with experiments, and the segment order parameters were in agreement with those of the nuclear magnetic resonance (NMR) experiments. The two kinds of lipid chains of DPPC do not equivalently behave in the bilayers, and chain 2 has lower ordering than chain 1. The order parameters of the first eight segments of chain 2 in the liquid-crystal model are particularly small and are roughly constant. From electron density analysis, it has been observed that the liquid-crystal bilayer has about one excess water molecule per one lipid molecule in comparison with the gel bilayer. The energy difference between the two bilayer models, taking account of the water contribution, is consistent with the latent heat of the phase transition. © 1995 by John Wiley & Sons, Inc.     

Structure analysis of microporous plane membranes by a Monte Carlo method

Plane microporous PTFE membranes, obtained by stretching, have been modelled by layers of cylinders representing the polymer filaments. The concept of the simulation of Knudsen diffusion in the networks by a Monte Carlo method is used to estimate the complexity of the three-dimensional structure of these membranes. A tortuosity parameter can be calculated from a relation depending only upon the volumic porosity and the mean filament diameter. The results obtained are compared to data obtained for real membranes.     

Molten globule state in homopolymers:A Monte Carlo study

Simple models of polymer chains were based on a simple cubic lattice. The model chains were star‐branched with f = 3 and f = 6 branches. The attractive potential between polymer segments was introduced to study the properties of polymer chains in the different temperature regimes. The computer simulations were carried out by means of the dynamic Monte Carlo method. It was found that contrary to recent real experiments, the ratio of the radius of gyration to the hydrodynamic radius did not exhibit a maximum near the coil‐globule transition but decreased monotonically with the temperature. The distribution of polymer‐polymer contacts and their lifetimes were also studied. It appeared that in homopolymer chains the lifetimes of these contacts were very short. At low temperatures contacts were distributed over the entire chain and at high temperatures only contacts that were close to the chain survived longer times.     

Polyion adsorption onto catanionic surfaces. A Monte Carlo study

The adsorption of a single and negatively charged polyion with varying flexibility onto a surface carrying both negative and positive charges representing a charged membrane surface has been investigated by using a simple model employing Monte Carlo simulations. The polyion was represented by a sequence of negatively charged hard spheres connected with harmonic bonds. The charged surface groups were also represented by charged hard spheres, and they were positioned on a hard surface slightly protruding into the solution. The surface charges were either frozen in a liquidlike structure or laterally mobile. With a large excess of positive surface charges, the classical picture of a strongly adsorbed polyion with an extended and flat configuration emerged. However, adsorption also appeared at a net neutral surface or at a weakly negatively charged surface, and at these conditions the adsorption was stronger with a flexible polyion as compared to a semiflexible one, two features not appearing in simpler models containing homogeneously charged surfaces. The presence of charged surface patches (frozen surface charges) and the ability of polarization of the surface charges (mobile surface charges) are the main reasons for the enhanced adsorption. The stronger adsorption with the flexible chain is caused by its greater ability to spatially correlate with the surface charges.     

Structure, dynamics, and phase transitions of tethered membranes: a Monte Carlo simulation study

A coarse-grained model of a self-avoiding tethered membrane with hexagonal coordination, embedded in three-dimensional space, is studied by means of extensive Monte Carlo computer simulations. The simulations are performed at various temperatures for membranes with linear size 5< or =L< or =50. We find that the membrane undergoes several folding transitions from a high-temperature flat phase to multiple-folded structure as the temperature is steadily decreased. Using a suitable order parameter and finite size scaling analysis, these phase transitions are shown to be of first order. The equilibrium shape of the membranes is analyzed by calculating the eigenvalues lambda(max) (2)> or =lambda(med) (2)> or =lambda(min) (2) of the inertia tensor. We present a systematic finite size scaling analysis of the radius of gyration and the eigenvalues of the inertia tensor at different phases of the observed folding transitions. In the high-temperature flat phase, the radius of gyration R(g) grows with the linear size of the membrane L as R(g) proportional to L(nu), where the exponent nu is approximately equal to 1.0. The eigenvalues of the inertia tensor scale as lambda(max) proportional to lambda(med) proportional to L(nu) and lambda(min) proportional to L(nu(min) ), whereby the roughness exponent nu(min) is approximately equal to 0.7. We also find that the time tau(R) of a self-avoiding membrane to diffuse a distance R(g) scales as tau(R) proportional to L(2nu+2), which is in good agreement with the theoretical predictions.     

Protection of nano-powders by adsorbed surfactants: A Monte Carlo study

The production of high-performance ceramics requires the protection of powder particles against chemical reactions. Hydrolysis and oxidation of nanoscaled non-oxidic powders can be impeded by a coating consisting of a dense adsorbed layer of amphiphilic molecules. Using Monte Carlo simulations for a coarse grained model the adsorption equilibrium of differently shaped amphiphiles in apolar and polar solvents is investigated. For estimating the protection capability of the adsorbed surfactant film in aqueous environment we study the diffusion of small hydrophilic particles through the adsorbed surfactant film. The surfactants considered as coating agents differ in the number of hydrocarbon tails. It is found that amphiphiles with a single hydrocarbon tail or at most two branches are more suitable to protect particle surfaces than amphiphiles with three or four branches, although the adsorption energy of amphiphiles with many branches is higher.     

Interfacial properties of glassy polymer melts: A Monte Carlo study

The properties of the interface between a polymer melt and a solid wall are studied over a wide range of temperatures by dynamic Monte Carlo simulations. It is shown that in the supercooled state near the glass transition of the melt an “interphase” forms, the structure of which is influenced by the wall. The thickness of this interphase is determined from the monomer density profile near the surface and is strongly temperature dependent. At low glass-like temperatures it is larger than the bulk radius of gyration of the chains.     

A Monte Carlo study of the mesophases formed by polar bent-shaped molecules

Liquid crystal phases formed by bent-shaped (or "banana") molecules are currently of great interest. Here we investigate by Monte Carlo computer simulations the phases formed by rigid banana molecules modeled combining three Gay-Berne sites and containing either one central or two lateral and transversal dipoles. We show that changing the dipole position and orientation has a profound effect on the mesophase stability and molecular organization. In particular, we find a uniaxial nematic phase only for off-center dipolar models and tilted phases only for the one with terminal dipoles.     

Grafting of maleic anhydride onto linear polyethylene: A Monte Carlo study

Monte Carlo simulation was used to study the graft of maleic anhydride (MAH) onto linear polyethylene (PE‐g‐MAH) initiated by dicumyl peroxide (DCP). Simulation results revealed that major MAH monomers attached onto PE chains as branched graft at higher MAH content. However, at extremely low MAH content, the fraction of bridged graft was very close to that of branched graft. This conclusion was somewhat different from the conventional viewpoint, namely, the fraction of bridged graft was always much lower than that of branched graft under any condition. Moreover, the results indicated that the grafting degree increased almost linearly to MAH and DCP concentrations. On the other hand, it was found that the amount of grafted MAH dropped sharply with increasing the length of grafted MAH, indicating that MAH monomers were mainly attached onto the PE chain as single MAH groups or very short oligomers. With respect to the crosslink of PE, the results showed that the fraction of PE‐(MAH)_n‐PE crosslink structure increased continuously, and hence the fraction of PE‐PE crosslink decreased with increasing MAH concentration. Finally, quantitative relationship among number average molecular weight of the PE, MAH, and DCP contents was given. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5714–5724, 2004