scholarly journals A 3D particle Monte Carlo approach to studying nucleation

2017 ◽  
Author(s):  
Christoph Köhn ◽  
Martin Bødker Enghoff ◽  
Henrik Svensmark

Abstract. The nucleation of sulphuric acid molecules plays a key role in the formation of aerosols. We here present a three dimensional particle Monte Carlo model to study the growth of sulphuric acid clusters as well as its dependence on the ambient temperature and the initial particle density.We initiate a swarm of sulphuric acid molecules with a size of 0.15 nm with densities between 107 and 108 cm−3 at temperatures of 200 and 300 K. After every time step, we update the position and velocity of particles as a function of size-dependent diffusion coefficients. If two particles encounter, we merge them and add their volumes and masses. Inversely, we check after every time step whether a polymer evaporates liberating a molecule.We present the spatial distribution as well as the size distribution calculated from individual clusters. We also calculate the nucleation rate of clusters with a radius of 0.85 nm as a function of time, initial particle density and temperature. For 200 K, the nucleation rate increases as a function of time; for 300 K we observe an interplay between clustering and evaporation and thus the oscillation of the nucleation rate around the mean nucleation rate. The nucleation rates obtained from the presented model agree well with experimentally obtained values which serves as a benchmark of our code. In contrast to previous nucleation models, we here present for the first time a code capable of tracing individual particles and thus of capturing the physics related to the discrete nature of particles.

1988 ◽  
Vol 141 ◽  
Author(s):  
James P. Lavine ◽  
Gilbert A. Hawkins

AbstractA three-dimensional Monte Carlo computer program has been developed to study the heterogeneous nucleation and growth of oxide precipitates during the thermal treatment of crystalline silicon. In the simulations, oxygen atoms move on a lattice with randomly selected lattice points serving as nucleation sites. The change in free energy that the oxygen cluster would experience in gaining or losing one oxygen atom is used to govern growth or dissolution of the cluster. All the oxygen atoms undergo a jump or a growth decision during each time step of the anneal. The growth and decay kinetics of each nucleation site display interesting fluctuation phenomena. The time dependence of the cluster size generally differs from the expected 3/2 power law due to the fluctuations in oxygen arrival at and incorporation in a precipitate. Competition between growing sites and coarsening are observed.


2021 ◽  
Vol 36 (25) ◽  
pp. 2150182
Author(s):  
Khusniddin K. Olimov ◽  
Vladimir V. Lugovoi ◽  
Kosim Olimov ◽  
Maratbek Shodmonov ◽  
Kadyr G. Gulamov ◽  
...  

To describe [Formula: see text] interactions with production of three [Formula: see text]-particles at incident neutron kinetic energy of 14 MeV in a nuclear (photo) emulsion, a Monte Carlo model is proposed for four channels of decay of an excited carbon-12 nucleus into three [Formula: see text]-particles. The Monte Carlo calculation results describe well the experimental data on the distribution of the angle between the three-dimensional momenta of all pairs of [Formula: see text]-particles in a collision event, on the distribution of the angle between the projections of the momentum vectors of all pairs of [Formula: see text]-particles in collision event on each of the coordinate planes, on the distribution of the sum of the kinetic energies of all pairs of [Formula: see text]-particles in a collision event, and the distribution of projections of the momenta of [Formula: see text]-particles on the coordinate planes. The best agreement of the Monte Carlo model results with the experimental data is achieved if the direct decay [Formula: see text] and decay through the formation of an intermediate beryllium nucleus [Formula: see text] are generated with equal probabilities, while the excitation energies of 3.04 MeV, 1.04 MeV, and 0.1 MeV for the beryllium nucleus are generated with relative weights of 75%, 15%, and 10%, respectively.


2013 ◽  
Vol 1559 ◽  
Author(s):  
Andreas Latz ◽  
Lothar Brendel ◽  
Dietrich E. Wolf

ABSTRACTWhile the self-learning kinetic Monte Carlo (SLKMC) method enables the calculation of transition rates from a realistic potential, implementations of it were usually limited to one specific surface orientation. An example is the fcc (111) surface in Latz et al. 2012, J. Phys.: Condens. Matter 24, 485005. This work provides an extension by means of detecting the local orientation, and thus allows for the accurate simulation of arbitrarily shaped surfaces. We applied the model to the diffusion of Ag monolayer islands and voids on a Ag(111) and Ag(001) surface, as well as the relaxation of a three-dimensional spherical particle.


2014 ◽  
Vol 1704 ◽  
Author(s):  
Hayley R. Osman ◽  
Saibal Mitra

ABSTRACTInitiated Chemical Vapor Deposition (iCVD) is a well-known method for depositing polymers that are used in chemical, biological, and electrical applications. It is a variation of hot filament deposition and can used to produce conformal coatings of polymer films at relatively low reaction temperatures. It is also a solventless technique in which thin polymeric films are deposited by introducing controlled ratios of monomer and initiator gasses into the reaction chamber. Low temperatures in the reaction chamber allow the deposition of polymer films on a wide variety of substrates that include biological substrates.We have simulated the growth of a monolayer of polymer films on two-dimensional surfaces using Monte Carlo simulation. We saw the formation of polymer chains over a time scale on the order of microseconds. We have assumed the substrate to be at room temperature while the reactor pressure close of 800 mTorr.The grid on which we have simulated this polymer growth is represented by a 100x100 matrix, on which a series of specialized functions are executed in each time-step, or iteration. These functions can be divided into three categories: population, translation, and polymerization.The goal of this simulation is to observe the initial growth of the iCVD surface reaction. We have obtained favorable results with the simulation and we are now looking to compare these results with experimental results for initiation growth.


1997 ◽  
Vol 490 ◽  
Author(s):  
Myung-Sik Son ◽  
Ho-Jung Hwang

ABSTRACTAn alternative three-dimensional (3D) Monte Carlo (MC) dynamic simulation model for phosphorus implant into (100) single-crystal silicon has been developed which incorporates the effects of channeling and damage. This model calculates the trajectories of both implanted ions and recoiled silicons and concurrently and explicitly affects both ions and recoils due to the presence of accumulative damage. In addition, the model for room-temperature implant accounts for the self-annealing effect using our defined recombination probabilities for vacancies and interstitials saved on the unit volumes. Our model has been verified by the comparison with the previously published SIMS data over commonly used energy range between 10 and 180 keV, using our proposed empirical electronic energy loss model. The 3D formations of the amorphous region and the ultra-shallow junction around the implanted region could be predicted by using our model, TRICSI (TRansport Ions into Crystal-Silicon).


1982 ◽  
Vol 49 (18) ◽  
pp. 1299-1302 ◽  
Author(s):  
Dietrich Stauffer ◽  
Antonio Coniglio ◽  
Dieter W. Heermann

2009 ◽  
Vol 66 (10) ◽  
pp. 3131-3146 ◽  
Author(s):  
Robert Pincus ◽  
K. Franklin Evans

Abstract This paper examines the tradeoffs between computational cost and accuracy for two new state-of-the-art codes for computing three-dimensional radiative transfer: a community Monte Carlo model and a parallel implementation of the Spherical Harmonics Discrete Ordinate Method (SHDOM). Both codes are described and algorithmic choices are elaborated. Two prototype problems are considered: a domain filled with stratocumulus clouds and another containing scattered shallow cumulus, absorbing aerosols, and molecular scatterers. Calculations are performed for a range of resolutions and the relationships between accuracy and computational cost, measured by memory use and time to solution, are compared. Monte Carlo accuracy depends primarily on the number of trajectories used in the integration. Monte Carlo estimates of intensity are computationally expensive and may be subject to large sampling noise from highly peaked phase functions. This noise can be decreased using a range of variance reduction techniques, but these techniques can compromise the excellent agreement between the true error and estimates obtained from unbiased calculations. SHDOM accuracy is controlled by both spatial and angular resolution; different output fields are sensitive to different aspects of this resolution, so the optimum accuracy parameters depend on which quantities are desired as well as on the characteristics of the problem being solved. The accuracy of SHDOM must be assessed through convergence tests and all results from unconverged solutions may be biased. SHDOM is more efficient (i.e., has lower error for a given computational cost) than Monte Carlo when computing pixel-by-pixel upwelling fluxes in the cumulus scene, whereas Monte Carlo is more efficient in computing flux divergence and downwelling flux in the stratocumulus scene, especially at higher accuracies. The two models are comparable for downwelling flux and flux divergence in cumulus and upwelling flux in stratocumulus. SHDOM is substantially more efficient when computing pixel-by-pixel intensity in multiple directions; the models are comparable when computing domain-average intensities. In some cases memory use, rather than computation time, may limit the resolution of SHDOM calculations.


2020 ◽  
Vol 24 ◽  
pp. 472-485 ◽  
Author(s):  
Vigneshwaran Thangavel ◽  
Oscar Xavier Guerrero ◽  
Matias Quiroga ◽  
Adelphe Matsiegui Mikala ◽  
Alexis Rucci ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document