Search results for " Monte Carlo method"
showing 10 items of 107 documents
Phase transitions in nonadditive hard disc systems: a Gibbs ensemble Monte Carlo Study
2007
we study the properties of a model fluid in two dimensions with Gibbs ensemble Monte Carlo (GEMC) techniques, in particular we analyze the entropy-driven phase separation in case of a nonadditive symmetric hard disc fluid. By a combination of GEMC with finite size scaling techniques we locate the critical line of nonadditivities as a function of the system density, which separates the mixing/demixing regions and compare with a simple analytical approximation.
Monte Carlo Simulations of Polymer Systems
1988
The impact of Monte Carlo “computer experiments” in polymer physics is described, emphasizing three examples taken from the author’s research group. The first example is a test of the classical Flory—Huggins theory for polymer mixtures, including a discussion of cricital phenomena. Also “technical aspects” of such simulations (“grand-canonical” ensemble, finite—size scaling, etc.) are explained briefly. The second example refers to configurational statistics and dynamics of chains confined to cylindrical tubes; the third example deals with the adsorption of polymers at walls. These simulations check scaling concepts developed along the lines of de Gennes.
Monte Carlo Simulations in Polymer Science
2012
Monte Carlo methods are useful for computing the statistical properties of both single macromolecules of various chemical architectures and systems containing many polymers (solutions, melts, blends, etc.). Starting with simple models (lattice models such as the self-avoiding walk or the bond fluctuation model, as well as coarse-grained or chemically realistic models in the continuum) various algorithms exist to generate conformations typical for thermal equilibrium, but dynamic Monte Carlo methods can also model diffusion and relaxation processes (as described by the Rouse and the reptation models for polymer melt dynamics). Limitations of the method are explained, and also the measures to…
Monte Carlo Simulation of Langmuir Monolayer Models
1998
We study a coarse grained, continuum model for Langmuir monolayers, i.e., monolayers of amphiphilic molecules on a polar substrate. The amphiphiles are represented by stiff chains of beads with one end grafted to a planar surface. Monte Carlo Simulations at constant pressure have been performed, using simulation boxes of variable size and variable shape. A number of techniques have been explored in order to obtain an efficient simulation algorithm. We discuss the resulting phase diagram, characterize the different phases, and analyze the conditions, under which they can be found.
Monte Carlo Simulations of Growth Kinetics and Phase Transitions at Interfaces: Some Recent Results
1991
ABSTRACTIn the first part Monte Carlo studies of the kinetics of multilayer adsorption (without screening) are described. The approach to the jamming coverage in each layer is asymptotically exponential. The jamming coverages approach the infinite-layer limit value according to a power law. In the second part, studies of phase transitions in two dimensional fluids are reviewed. With a combination of Monte Carlo and finite size scaling block analysis techniques, accurate values are obtained for the critical temperatures, coexistence densities and the compressibilities of an adsorbed fluid layer in an NVT ensemble.
Isotropic–isotropic phase separation in mixtures of rods and spheres: Some aspects of Monte Carlo simulation in the grand canonical ensemble
2008
Abstract In this article we consider mixtures of non-adsorbing polymers and rod-like colloids in the isotropic phase, which upon the addition of polymers show an effective attraction via depletion forces. Above a certain concentration, the depletant causes phase separation of the mixture. We performed Monte Carlo simulations to estimate the phase boundaries of isotropic–isotropic coexistence. To determine the phase boundaries we simulated in the grand canonical ensemble using successive umbrella sampling [J. Chem. Phys. 120 (2004) 10925]. The location of the critical point was estimated by a finite size scaling analysis. In order to equilibrate the system efficiently, we used a cluster move…
Monte Carlo study of surface critical behavior in the XY model.
1989
We have used Monte Carlo simulations to study the behavior of $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}D$ slabs containing classical spins which interact via nearest-neighbor $\mathrm{XY}$ coupling. The coupling constant ${J}_{S}$ for spins in the surface layer is fixed at $0.5J$. Finite-size scaling is used to analyze data for $D=59$ and to extract estimates for the surface critical exponents. We find that ${\ensuremath{\beta}}_{1}$ is in good agreement with theoretical predictions.
HIGH-PRECISION MONTE CARLO DETERMINATION OF α/ν IN THE 3D CLASSICAL HEISENBERG MODEL
1994
To study the role of topological defects in the three-dimensional classical Heisenberg model we have simulated this model on simple cubic lattices of size up to 803, using the single-cluster Monte Carlo update. Analysing the specific-heat data of these simulations, we obtain a very accurate estimate for the ratio of the specific-heat exponent with the correlation-length exponent, α/ν, from a usual finite-size scaling analysis at the critical coupling Kc. Moreover, by fitting the energy at Kc, we reduce the error estimates by another factor of two, and get a value of α/ν, which is comparable in accuracy to best field theoretic estimates.
Study of the helium-cooled lithium lead test blanket module nuclear behaviour under irradiation in ITER
2009
Abstract The present paper deals with the detailed investigation of the helium-cooled lithium lead test blanket module (HCLL-TBM) nuclear behaviour under irradiation in ITER, carried out at the Department of Nuclear Engineering of the University of Palermo adopting a numerical approach based on the Monte Carlo method. A realistic 3D heterogeneous model of the HCLL-TBM was set-up and inserted into an ITER 3D semi-heterogeneous model that realistically simulates the reactor lay-out up to the cryostat. A Gaussian-shaped neutron source was adopted for the calculations. The main features of the HCLL-TBM nuclear response were assessed, paying a particular attention to the neutronic and photonic d…
Magic numbers, excitation levels, and other properties of small neutral math clusters (N < 50)
2006
The ground-state energies and the radial and pair distribution functions of neutral math clusters are systematically calculated by the diffusion Monte Carlo method in steps of one math atom from 3 to 50 atoms. In addition the chemical potential and the low-lying excitation levels of each cluster are determined with high precision. These calculations reveal that the “magic numbers” observed in experimental math cluster size distributions, measured for free jet gas expansions by nondestructive matter-wave diffraction, are not caused by enhanced stabilities. Instead they are explained in terms of an enhanced growth due to sharp peaks in the equilibrium concentrations in the early part of the e…