Search results for "Monte Carlo method"
showing 10 items of 1234 documents
A dosimetric study of Leipzig applicators.
2005
Purpose: To obtain the absolute dose-rate distribution in liquid water for all six cup-shaped Leipzig applicators by means of an experimentally validated Monte Carlo (MC) code. These six applicators were used in high-dose-rate (HDR) afterloaders with the “classic” and v2 192Ir sources. The applicators have an inner diameter of 1, 2, and 3 cm, with the source traveling parallel or perpendicular to the contact surface. Methods and materials: The MC GEANT4 code was used to obtain the dose-rate distribution in liquid water for the six applicators and the two HDR source models. To normalize the applicator output factors, a MC simulation for the “classic” and v2 sources in air was performed to es…
Monte Carlo Simulation of Alloy Phase Diagrams and Short-Range Order
1986
As a prototype model for order-disorder phenomena in binary alloys, a face-centered cubic lattice is considered,the sites of which can be taken by either A-atoms or B-atoms, assuming pair-wise interactions between nearest (J) and next nearest neighbours (J). The phase diagram is constructed from Monte Carlo calculations. Some technical aspects essential for the success of such calculations are briefly mentioned (use of grand-canonical rather than canonical ensemble, how to obtain the free energy needed to locate first-order phase transitions, etc.). It is shown that the topology of the phase diagram changes when the ratio R = Jnnn/Jnn is varied, and this behaviour is discussed in the contex…
Replica-exchange molecular dynamics simulation for supercooled liquids
2000
We investigate to what extend the replica-exchange Monte Carlo method is able to equilibrate a simple liquid in its supercooled state. We find that this method does indeed allow to generate accurately the canonical distribution function even at low temperatures and that its efficiency is about 10-100 times higher than the usual canonical molecular dynamics simulation.
Classical and Quantum Two-Dimensional Fluids in the Gibbs Ensemble
1994
We study the properties of model fluids in two spatial dimensions with Gibbs ensemble Monte Carlo (GEMC) techniques. In particular in the first part of the paper we study the entropy driven phase separation in case of a nonadditive symmetric hard disc fluid and locate by a combination of GEMC with finite size scaling techniques the critical line of nonadditivities as a function of the system density, which separates the mixing/demixing regions, we compare with a simple approximation. In the second part we successfully combine path integral Monte Carlo (PIMC) and GEMC techniques in order to locate the gas-liquid coexistence densities for a fluid with classical degrees of freedom and internal…
Anomalous size-dependence of interfacial profiles between coexisting phases of polymer mixtures in thin-film geometry: A Monte Carlo simulation
1997
The interfacial profile between coexisting phases of a binary mixture (A,B) in a thin film of thickness D and lateral linear dimensions L depends sensitively on both linear dimensions and on the nature of boundary conditions and statistical ensembles applied. These phenomena generic for systems in confined geometry are demonstrated by Monte-Carlo simulations of the bond fluctuation model of symmetric polymer mixtures. Both the canonical and semi-grand-canonical ensemble are studied. In the canonical ensemble, the interfacial width w increases (from small values which are of the same order as the intrinsic profile) like sqrt{D}, before a crossover to a saturation value w_max (w_max^2 proport…
Numerical and experimental verification of a technique for locating a fatigue crack on beams vibrating under Gaussian excitation
2007
The stationary vibrations of a beam excited by Gaussian noise are strongly affected by the presence of a fatigue crack. Indeed, as soon as the crack arises the system response becomes non-linear due to crack breathing and a non-Gaussian behaviour is encountered. The paper presents both numerical and experimental investigations in order to assess the capability of the non-Gaussianity measures to detect crack presence and position. Monte Carlo method is applied to evaluate in time domain the higher order statistics of a cantilever beam modelled by finite elements. The skewness coefficient of the rotational degrees of freedom appears the most suitable quantity for identification purpose being …
Isotropic-nematic interface in suspensions of hard rods: Mean-field properties and capillary waves
2006
We present a study of the isotropic-nematic interface in a system of hard spherocylinders. First we compare results from Monte Carlo simulations and Onsager density functional theory for the interfacial profiles of the orientational order parameter and the density. Those interfacial properties that are not affected by capillary waves are in good agreement, despite the fact that Onsager theory overestimates the coexistence densities. Then we show results of a Monte Carlo study of the capillary waves of the interface. In agreement with recent theoretical investigations (Eur.Phys.J. E {\bf 18} 407 (2005)) we find a strongly anistropic capillary wave spectrum. For the wave-numbers accessed in o…
Capillary Waves in a Colloid-Polymer Interface
2004
The structure and the statistical fluctuations of interfaces between coexisting phases in the Asakura-Oosawa (AO) model for a colloid--polymer mixture are analyzed by extensive Monte Carlo simulations. We make use of a recently developed grand canonical cluster move with an additional constraint stabilizing the existence of two interfaces in the (rectangular) box that is simulated. Choosing very large systems, of size LxLxD with L=60 and D=120, measured in units of the colloid radius, the spectrum of capillary wave-type interfacial excitations is analyzed in detail. The local position of the interface is defined in terms of a (local) Gibbs surface concept. For small wavevectors capillary wa…
Computer simulation studies of finite-size broadening of solid–liquid interfaces: from hard spheres to nickel
2009
Using Molecular Dynamics (MD) and Monte Carlo (MC) simulations interfacial properties of crystal-fluid interfaces are investigated for the hard sphere system and the one-component metallic system Ni (the latter modeled by a potential of the embedded atom type). Different local order parameters are considered to obtain order parameter profiles for systems where the crystal phase is in coexistence with the fluid phase, separated by interfaces with (100) orientation of the crystal. From these profiles, the mean-squared interfacial width w^2 is extracted as a function of system size. We rationalize the prediction of capillary wave theory that w^2 diverges logarithmically with the lateral size o…
Critical point and coexistence curve properties of the Lennard-Jones fluid: A finite-size scaling study
1995
Monte Carlo simulations within the grand canonical ensemble are used to explore the liquid-vapour coexistence curve and critical point properties of the Lennard-Jones fluid. Attention is focused on the joint distribution of density and energy fluctuations at coexistence. In the vicinity of the critical point, this distribution is analysed using mixed-field finite-size scaling techniques aided by histogram reweighting methods. The analysis yields highly accurate estimates of the critical point parameters, as well as exposing the size and character of corrections to scaling. In the sub-critical coexistence region the density distribution is obtained by combining multicanonical simulations wit…