Computation of partial molar properties using continuous fractional component Monte Carlo

ABSTRACT

An alternative method for calculating partial molar excess enthalpies and partial molar volumes of components in Monte Carlo (MC) simulations is developed. This method combines the original idea of Frenkel, Ciccotti, and co-workers with the recent continuous fractional component Monte Carlo (CFCMC) technique. The method is tested for a system of Lennard–Jones particles at different densities. As an example of a realistic system, partial molar properties of a [NH₃, N₂, H₂] mixture at chemical equilibrium are computed at different pressures ranging from P = 10 to 80 MPa. Results obtained from MC simulations are compared to those obtained from the PC-SAFT Equation of State (EoS) and the Peng–Robinson EoS. Excellent agreement is found between the results obtained from MC simulations and PC-SAFT EoS, and significant differences were found for PR EoS modelling. We find that the reaction is much more exothermic at higher pressures.     

气泡幕后向散射光信号特性的蒙特卡罗方法研究

利用蒙特卡罗方法对水介质中的气泡幕的后向光散射回波信号进行了系统仿真。通过对计算结果与实验结果的比较,证明了蒙特卡罗方法的有效性。根据仿真结果,结合理论分析的方法,对气泡幕位置、厚度和接收器视场角等参数对回波信号的影响进行了分析,结果表明:回波信号出现时间与气泡幕位置存在一一对应的关系;在气泡幕的衰减系数ρσt不变的情况下,存在一个有效气泡幕厚度,当大于该厚度时,气泡幕的后向光散射回波信号基本不变;在系统各项参数不变的情况下,适当增加接收器的视场角,可有效地提高回波信号的信噪比。     

Thermodynamic properties of ternary alloys from Monte Carlo simulations

Methods are presented for obtaining both the integral and partial molar energies, entropies and free energies of the components in rigid ternary substitutional alloy systems from the one Monte Carlo simulation. Tests of the methods for some model systems are presented and discussed.     

Magnetic properties for cobalt nanorings: Monte Carlo simulation

In this paper, two structure models of cobalt nanoring cells (double-nanorings and four-nanorings, named as D-rings and F-rings, respectively) have been considered. Base on Monte Carlo simulation, the magnetic properties of the D-rings and F-rings, such as hysteresis loops, spin configuration, coercivity, etc., have been studied. The simulated results indicate that both D-rings and F-rings with different inner radius (r) and separation of ring centers (d) display interesting magnetization behavior and spin configurations (onion-, vortex- and crescent shape vortex-type states) in magnetization process. Moreover, it is found that the overlap between the nearest single nanorings connect can result in the deviation of the vortex-type states in the connected regions. Therefore, the appropriate d should be well considered in the design of nanoring device. The simulated results can be explained by the competition between exchange energy and dipolar energy in Co nanorings system. Furthermore, it is found that the simulated temperature dependence of the coercivity for the D-rings with different d can be well described by H_c=H₀ exp[−(T/T₀)^p].     

R32热力学性质的蒙特卡罗模拟

建立在统计热力学和分子力学理论基础上的分子模拟方法逐渐运用于计算制冷剂的热力学性质。文中首先在NVT系综条件下,采用吉布斯蒙特卡罗模拟方法(GEMC),模拟了R32的气液相平衡的密度、饱和蒸汽压及蒸发焓;其次在NPT系综条件下,采用蒙特卡罗模拟方法(MC),模拟了R32在4MPa条件下的过冷液态密度。模拟结果同美国国家标准研究所(NIST)的Refprop 8.0相比,有很好的一致性。结果表明,运用该方法预测单一制冷剂的热力学性质是可行的。     

Performance characteristics of single photon pulse ranging system using Monte Carlo simulation

Single photon pulse ranging system with extremely high sensitivity has been widely used in distance measurement and 3D imaging. To analyze the factors that affect the measurement precision and accuracy will help to improve system performance. According to system structure and principle, we mainly discussed the following factors: laser intensity, pulse width, detection efficiency and time jitter. A simulation model based on Monte Carlo stochastic method was constructed in this paper, and we get the specific influence of factors on measurement precision and accuracy by simulation. Finally, we set up laboratory experiment system and took effective experiments on ranging precision and accuracy.     

Investigation of scatter from out of the field of view and multiple scatter in PET using Monte Carlo simulations

In fully three-dimensional (3D) positron emission tomography (PET) imaging, the scatter fraction (SF) is about 40%-60%, which may degrade the imaging quality severely. Scatter correction is important for high quality image reconstruction. Model-based scatter correction has been proved to be accurate and available in clinical PET. However, it does not correct the scatter from out of the field of view (OFOV) and multiple scatters. In this study, we demonstrate the radial and axial distribution of scatters from OFOV when the source is located in different radial positions. In order to apply the above conclusions to different PET systems, we characterize the scatters from OFOV as a function of the ratio of the scanner diameter to the length of the axial field of view (AFOV) by modeling several typical whole-body and micro PET systems. The proportions of true events (S_0-0), single scatter of one photon (S_1-0) , single scatter of both photons (S_1-1) , double scatter of one photon (S_2-0) and multiple scatter (S_m) are also calculated and compared. Here the 3D-PET Monte Carlo simulations are performed with the Geant4 Application for Tomography Emission (GATE). In summary, the scatters from OFOV tend to be recorded on the lines of response (LOR) far away from the source. They have a much more serious impact on whole-body PET than micro PET depending on the ratio of scanner diameter to the length of AFOV. In whole-body PET, twice scatters including single scatter of both photons (S_1-1) and double scatter of one photon (S_2-0) add up to about 12% so that twice scatter correction must be taken into account to acquire a high quality reconstruction image.     

Magnetic phase diagram simulation of La1−xCaxMnO3 system by using Monte Carlo, Metropolis algorithm and Heisenberg model

This work presents a critical temperature study of La_1−xCa_xMnO₃ manganites in bulk by means of Monte Carlo method thermal activated magnetic properties. The analysis was carried out for stoichiometries in the range of 0≤x≤1. The model is based on a three-dimensional classical Heisenberg-Hamiltonian involving the presence of Mn^3+eg, Mn^3+eg′ and Mn⁴⁺ ions, and their nearest neighbor interaction. For this modeling, simple cubic lattice samples of size L³, with L=6, 15 and 30 were used. The values of exchange parameters were determined by using LaMnO₃ (x=0), La_0.5Ca_0.5MnO₃ and CaMnO₃ (x=1) phases. Relationships between exchange parameters and anisotropy constants for different hole densities were found. Results of transition temperatures for each phase showed good agreement with experimental reports, especially for L=30 and L→∞.     

铜原子纳米团簇热力学性质的分子动力学模拟研究

利用分子动力学模拟方法,研究了CuN(N=80,140,216,312,408,500,628和736)纳米团簇在热化和冷凝过程中结构和热力学性质的变化,模型采用的是Johnson的EAM作用势.模拟结果表明:铜团簇的熔点和凝固点随其尺寸线性增加,并逐渐向大块晶体靠拢;所有团簇的凝固点都低于熔点,出现凝固过程中的滞后现象;在熔点和凝固点附近,团簇都具有负热容特性,负热容是由相变前后团簇内部结构突变引起的.     

Thermodynamic properties of the Lennard-Jones FCC solid: perturbation theory parameterisation and Monte Carlo simulation

This work is concerned with a valid representation of the solid-phase equation of state (EOS), the validity of which is evaluated by comparing to Monte Carlo (MC) simulation results. The proposed EOS has been developed by employing an optimal division of the Lannard-Jones (LJ) potential and an effective temperature- and density-dependent diameter into the framework of the simplified perturbation theory. Then, with the aim of extending to the chain systems, the conventional chain contribution (i.e. TPT1) is added to the proposed model (i.e. the atomic LJ system). Finally, the solid-state EOS based on Helmholtz free energy will be introduced for low temperature and high density conditions. To verify the accuracy of the proposed model, its performance is compared with the results of MC simulation. The comparison between the obtained results from the proposed model and the MC simulations shows that the EOS can satisfactorily predict the properties of the solid LJ system, both for the atomic system and for the chains.     

用分子自组装技术制备的单电子器件的Monte Carlo模拟

用分子自组装技术制备出纳米金单电子器件,并测量了其伏安特性,根据单电子系统的半经典理论,用MonteCarlo法对其结果进行了模拟.结果表明,模拟出的伏安曲线与实测的伏安曲线有较好的一致性,反映了模拟方法用于单电子器件研究的合理性,此外发现,虽然单电子器件两电极间含有众多的纳米粒子,但在低压区,其伏安特性只与少数纳米粒子有关 关键词： 单电子器件 MonteCarlo模拟 分子自组装     

A geometrical model for the double octahedral group based on a genus two Riemann surface

Geometrical models of double point groups are provided by double-sheeted Riemann surfaces of nonzero genus. This avoids the geometrically confusing picture of the doubling operation R as a rotation by 2π (i.e. 360°), which should instead be the identity operation E. In this manner a Riemann surface of genus 2 doubly covering a sphere and platonically tessellated b 16 equilateral triangles provides a geometrical model for the double octahedral group ²O_h. Stretching this Riemann surface along one axis to convert the 16 equilateral triangles to isosceles triangles and the underlying sphere to a prolate ellipsoid provides a model for the ²D_4h double group arising from the Jahn-Teller elongation of the regular octahedron into a prolate tetragonal bipyramid.     

Simulation of adsorbate-induced faceting on curved surfaces

A simple solid-on-solid model, proposed earlier to describe overlayer-induced faceting of bcc(1 1 1) surface, is applied to faceting of curved surfaces covered by an adsorbate monolayer. Surfaces studied in this paper are formed by a part of sphere around the [1 1 1] pole. Results of Monte Carlo simulations show that the morphology of a faceted surface depends on the annealing temperature. At an initial stage the surface around the [1 1 1] pole consists of 3-sided pyramids and step-like facets, then step-like facets dominate and their number decreases with temperature, finally a single big pyramid is formed. It is shown that there is a reversible phase transition at which a faceted surface transforms to an almost spherical one. It is found that the temperature of this phase transition is an increasing function of the surface curvature. Simulation results show that measurements of high temperature properties performed directly and after fast cooling down to a low temperature lead to different results.     

Monte Carlo simulation of supercooled liquids using a self-consistent local temperature

We combine Creutz energy conservation with Kawasaki spin exchange to simulate the microcanonical dynamics of a system of interacting particles. Relaxation occurs via Glauber spin-flip activation using a self-consistent temperature. Heterogeneity in the dynamics comes from finite-size constraints on the spin exchange that yield a distribution of correlated regions. The simulation produces a high-frequency response that can be identified with the boson peak, and a lower-frequency peak that contains non-Debye relaxation and non-Arrhenius activation, similar to the primary response of supercooled liquids.     

Molecular simulation of chemical reaction equilibria by Kinetic Monte Carlo

ABSTRACT

We describe a new algorithm for the molecular simulation of chemical reaction equilibria, which we call the Reactive Kinetic Monte Carlo (ReKMC) algorithm. It is based on the use of the equilibrium Kinetic Monte Carlo (eKMC) method (Ustinov et al., J. Colloid Interface Sci., 2012, 366, 216–223) to generate configurations in the underlying nonreacting system and to calculate the species chemical potentials at essentially zero marginal computational cost. We consider in detail the typical case of specified temperature, T and pressure, P, but extensions to other thermodynamic constraints are straightforward in principle. In the course of this work, we also demonstrate an alternative method for calculating simulation box volume changes in NPT ensemble simulations to achieve the specified P. We consider two sets of example reacting systems previously considered in the literature, and compare the ReKMC results and computational efficiencies with those of different implementations of the REMC algorithm (Turner et al., Molec. Simulation, 2008, 34, 119–146).     

溶液中晶体生长动力学的MC计算机模拟研究

用蒙特卡罗计算机模拟方法研究了溶液中晶体表面的生长台阶和生长动力学,着重讨论了表面尺寸、表面粗糙度以及溶液过饱和度对晶体生长速率的影响以及相应的生长机制.     

Fully parallel code for Monte Carlo simulation

We present a fully parallel version of Monte Carlo simulation of the Ising model using the Metropolis algorithm. In the 3-dimensional version the performance can be enhanced by a factor >20 in 16-bit word processors relative to other multispin codes. This factor could be further increased if implemented in 64-bit word computers.     

A simulation study of the e~+/e~- sensitivity of a low resistive phosphate glass electrode in a RPC using GEANT MC

The resistivity of conventional glass is quite high and is unacceptable in a high rate environment. Low resistive glass-electrodes could be a solution for this problem. The present study reports the e+/e-simulation results of an RPC detector made from low resistive phosphate glass electrodes. The detailed geometrical configuration of the content materials which are the essential components of the glass of the RPC detector have been created with the GEANT4 simulation code. Two different types of particle sources, i.e. for e+/e-, have been located on the detectors surface to evaluate the performance of the phosphate glass RPC. Both of the particles have been simulated as a function of their respective energies in the range of 0.1 MeV – 1.0 GeV. The present simulation work has shown that the resistive electrode plays an important role for the particle production in the RPC configuration.     

Computer simulations of critical phenomena and phase behaviour of fluids

Computer simulation techniques such as Monte Carlo (MC) and Molecular Dynamics (MD) methods yield numerically exact information (apart from statistical errors) on model systems of classical statistical mechanics. However, a systematic limitation is the restriction to a finite (and often rather small) particle number N (or box linear dimension L, respectively). This limitation is particularly restrictive near critical points (due to the divergence of the correlation length of the order parameter) and for the study of phase equilibria (possibly involving interfaces, droplets, etc.). Starting out with simple lattice gas (Ising) models, finite size scaling analyses have been developed to overcome this limitation. These techniques work for both simple Lennard-Jones fluids and their mixtures, including generalizations to approximate models for quadrupolar fluids such as carbon dioxide, benzene etc. and various mixtures, whose phase behaviour can be predicted. A combination of MC and MD allows the study of dynamic critical phenomena, and specialised techniques (umbrella sampling plus thermodynamic integration) yield the surface free energy of droplets as function of droplet size. Thus, computer simulation has become a versatile and widely applicable tool for the study of fluids.     

Monte Carlo simulation of the enantioseparation process

By means of Monte Carlo simulation, a study of enantioseparation by capillary electrophoresis has been carried out. A simplified system consisting of two enantiomers S (R) and a selector chiral C, which reacts with the enantiomers to form complexes RC (SC), has been considered. The dependence of Δμ$\Delta \mu$ (enantioseparation) with the concentration of chiral selector and with temperature have been analyzed by simulation. The effect of the binding constant and the charge of the complexes are also analyzed. The results are qualitatively satisfactory, despite the simplicity of the model.