Search results for "Monte Carlo method."
showing 10 items of 1217 documents
Monte Carlo Study of the Isotropic-Nematic Interface in Suspensions of Spherocylinders
2007
The isotropic to nematic transition in suspensions of anisotropic colloids is studied by means of grand canonical Monte Carlo simulation. From measurements of the grand canonical probability distribution of the particle density, the coexistence densities of the isotropic and the nematic phase are determined, as well as the interfacial tension.
How Do Droplets Depend on the System Size? Droplet Condensation and Nucleation in Small Simulation Cells
2003
Using large scale grandcanonical Monte Carlo simulations in junction with a multicanonical reweighting scheme we investigate the liquid-vapor transition of a Lennard—Jones fluid. Particular attention is focused on the free energy of droplets and the transition between different system configurations as the system tunnels between the vapor and the liquid state as a function of system size. The results highlight the finite size dependence of droplet properties in the canonical ensemble and free energy barriers along the path from the vapor to the liquid in the grandcanonical ensemble.
How do droplets on a surface depend on the system size?
2002
Abstract We investigate the thermodynamics of inhomogeneous polymer melts in the framework of a coarse grained off-lattice model. Properties of the liquid–vapour interface and the packing of the melt in contact with an attractive wall are considered. We employ Monte Carlo simulations in the grand canonical ensemble to determine excess free energies, the wetting temperature and the pre-wetting line, as well as the pre-wetting critical point. Having determined the wetting properties and the phase diagram of the model polymer, we perform canonical Monte Carlo simulations of small droplets on a surface. This allows us to study the dependence of droplet size on the wetting properties. It is foun…
Simulation of Models for Isotropic and Anisotropic Orientational Glasses
1992
“Orientational glass” behavior is found when molecular crystals are randomly diluted, and quadrupole moments get frozen by random alignment of the molecules, similar to “spin glass” behavior of randomly diluted magnets. Monte Carlo simulation of lattice models where quadrupole moments interact with nearest neighbor Gaussian coupling is a unique tool to study this behavior. The time-dependent glass order parameter exhibits anomalously slow relaxation, compatible with the Kohlrausch-Williams-Watts (KWW) stretched exponential function. Both isotropic and anisotropic models exhibit in d=2 and d=3 spatial dimensions glass transitions at zero temperature only. While the glass correlation length a…
The ensemble switch method for computing interfacial tensions
2015
We present a systematic thermodynamic integration approach to compute interfacial tensions for solid-liquid interfaces, which is based on the ensemble switch method. Applying Monte Carlo simulations and finite-size scaling techniques, we obtain results for hard spheres, which are in agreement with previous computations. The case of solid-liquid interfaces in a variant of the effective Asakura-Oosawa model and of liquid-vapor interfaces in the Lennard-Jones model are discussed as well. We demonstrate that a thorough finite-size analysis of the simulation data is required to obtain precise results for the interfacial tension.
Computer Simulations and Coarse-Grained Molecular Models Predicting the Equation of State of Polymer Solutions
2010
Monte Carlo and molecular dynamics simulations are, in principle, powerful tools for carrying out the basic task of statistical thermodynamics, namely the prediction of macroscopic properties of matter from suitable models of effective interactions between atoms and molecules. The state of the art of this approach is reviewed, with an emphasis on solutions of rather short polymer chains (such as alkanes) in various solvents. Several methods of constructing coarse-grained models of the simple bead–spring type will be mentioned, using input either from atomistic models (considering polybutadiene as an example) or from experiment. Also, the need to have corresponding coarse-grained models of t…
Adatom Island Diffusion on Metal Fcc(100) Surfaces
2001
We study the energetics and atomic mechanisms of diffusion of adatom islands on fcc(100) metal surfaces. For small islands, we perform detailed microscopic calculations using semi-empirical embedded-atom model and glue potentials in the case of Cu and Al, respectively. Combining systematic saddle-point search methods and molecular statics simulations allows us to find all the relevant transition paths for island motion. In particular, we demonstrate that there are novel many-body mechanisms such as internal row shearing which can, in some cases, control the island dynamics. Next, we show how using the master equation formalism, diffusion coefficients for small islands up to about five atoms…
Competition between submonolayer ordering and multilayer adsorption: Studies of simple lattice gas models
1986
Abstract We model condensation of adatoms at a substrate surface by a semi-infinite simple cubic lattice gas system. While in the bulk there is just a nearest-neighbour attractive interaction, in the first layer adjacent to the surface we allow for a periodic potential due to the substrate with a period of two lattice spacings, or for a next-nearest-neighbour repulsive interaction mediated by the substrate. Hence order-disorder phenomena may occur in the first layer, while only gas-liquid condensation transitions can occur in layers further away from the substrate surface. The ground-state phase diagrams of this model are obtained exactly, while the behaviour at nonzero temperatures is obta…
Incommensurate phases in adsorbed monolayers: structure and energy of domain walls
2002
Abstract The properties of incommensurate films of domain-wall structure formed on the (1 0 0) plane of face centered cubic crystals are studied by Monte Carlo simulation. The wall energies, wall structure and the wall–wall interaction are determined for different types of domain walls occurring in films which form the c(2×2) registered structure. The systems characterized by different strength and corrugation of the surface potential and of different misfit between adsorbate and adsorbent are discussed. It is demonstrated that heavy as well as light walls are rather strongly localized. Moreover, it is shown that the incommensurate structure with crossing heavy walls has higher stability th…
Monte Carlo study of cluster-diameter distribution: An observable to estimate correlation lengths
1997
We report numerical simulations of two-dimensional $q$-state Potts models with emphasis on a new quantity for the computation of spatial correlation lengths. This quantity is the cluster-diameter distribution function $G_{diam}(x)$, which measures the distribution of the diameter of stochastically defined cluster. Theoretically it is predicted to fall off exponentially for large diameter $x$, $G_{diam} \propto \exp(-x/\xi)$, where $\xi$ is the correlation length as usually defined through the large-distance behavior of two-point correlation functions. The results of our extensive Monte Carlo study in the disordered phase of the models with $q=10$, 15, and $20$ on large square lattices of si…