A Monte Carlo Study on the Identification of Quark and Gluon Jets

Three-jet events in e^ e^- collisions at 91.2 GeV are investigated using both HERWIG and JETSET Monte Carlo generators.The angles among the three jets are used to identify the quark and gluon jets.An analysis at the parton level is carried out to ensure the reasonableness of this method and an angular cut is utilized to improve the purity of this identification.The multiplicity inside the identified quark or gluon jets agrees with the quantum chromodynamic predictions qualitatively.     

Monte Carlo Study of the Neutron Detection Wall at RIBLLⅡ in HIRFL-CSR

With the construction of the new Radioactive Ion Beam Line in Lanzhou (RIBLLⅡ) which connects the CSRm and the CSRe, an experimental setup for physics research is highly required. A large area neutron detection wall is the main part of the setup. This paper introduced the detection principle of the neutron detection wall and the Monte-Carlo simulation of its design under the environment of the Geant4 toolkit. We presented the final design with the optimized parameters and the performance of the wall.     

Monte Carlo Study of CO-NO Catalytic Surface Reaction Including CO-CO Repulsion

The CO-NO reaction on a catalytic surface is studied by using Langmuir-Hinshclwood thermal mechanism with Monte Carlo computer simulation. In this model, a novel concept of CO CO repulsion is introduced, which has experimental evidence due to the formation of dipoles when these molecules are chemisorbed on the surface. The system is investigated by applying two approaches of NO dissociation. In the first ca.se, NO always decomposes into N and O before adsorption on the surface, In the second case, NO adsorbs on the surface molecularly and then dissociates into N and O if a vacancy is present in its adjacent neighbourhood. The steady state reactive window （i.e. the continuous production of CO2 and N2） is obtained only with the diffusion of N-atoms on the surface, which extends with CO-CO repulsion in the first, case. Itowever, in the second case, reactive window is obtained with CO-CO repulsion alone, The reactive window width in this case is reasonably large. The first-order phase transition is eliminated in both the cases with CO-CO repulsion.     

Monte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky-Moriya Interaction

With the help of an improvement Monte Carlo method, the Berezinskii-Kosterlitz-Thouless phase transition arising in two-dimensional planar rotator model with weak Dzyaloshinsky-Moriya （DM） interaction is investigated. The effects of the DM interaction on specific heat, susceptibility, and magnetization are simulated. The critical temperature of transitions is determined by the so-called Binder cumulant and the susceptibility of finite-size scaling. We find that the chiral Z2 symmetry reduced by the DM interactions plays an important role in a two-dimensional XY spin system, typically, the critical temperature is sensitive to weak DM spin couplings.     

Monte Carlo Investigation of Adsorption Stage of O on Ru（0001）：a Study of a Lattice Model

The adsorption of O on Ru(0001) is studied by means of Monte Carlo simulation of lattice gas model on a triangular lattice. A recent STM study shows that at low coverage the p(2×2) structure grows via island formation but the p(2×1) structure is abruptly formed at a critical coverage. Moreover, it also shows that there is a coexistence of the p(2×2) and p(2×1) structures. The above results seem not to coincide with the former studies of the system by both the LEED and Monte Carlo simulation. We therefore carried out the Monte Carlo study for the system again in the present paper and found that our simulation almost agrees with the results of the STM.     

A Monte Carlo generator of nucleon configurations in complex nuclei including nucleon–nucleon correlations

We developed a Monte Carlo event generator for production of nucleon configurations in complex nuclei consistently including effects of nucleon–nucleon (NN) correlations. Our approach is based on the Metropolis search for configurations satisfying essential constraints imposed by short- and long-range NN correlations, guided by the findings of realistic calculations of one- and two-body densities for medium-heavy nuclei. The produced event generator can be used for Monte Carlo (MC) studies of pA and AA collisions. We perform several tests of consistency of the code and comparison with previous models, in the case of high energy proton–nucleus scattering on an event-by-event basis, using nucleus configurations produced by our code and Glauber multiple scattering theory both for the uncorrelated and the correlated configurations; fluctuations of the average number of collisions are shown to be affected considerably by the introduction of NN correlations in the target nucleus. We also use the generator to estimate maximal possible gluon nuclear shadowing in a simple geometric model.     

Investigation of the energy dependence of the $$p_{\text {T0}}$$ parameter in the Pythia 8 Monte Carlo event generator

By using predictions from the Pythia 8 Monte Carlo event generator, we determine the energy-dependent turn-over values (\(p_{\text {T0}}\)) of the partonic cross section for the simulation of multiparton interactions. Since the observed energy dependence of the \(p_{\text {T0}}\) values is not well described by a power-law function, we introduce an additional energy-dependent term, to better describe the experimental observations in the energy range \(\sqrt{s} = 0.3-13 \, \hbox {TeV}\). We obtain a similar level of agreement for predictions using various parton densities.     

Multi-step magnetization of the Ising model on a Shastry-Sutherland lattice: a Monte?Carlo simulation

The magnetization behaviors and spin configurations of the classical Ising model on a Shastry-Sutherland lattice are investigated using Monte Carlo simulations, in order to understand the fascinating magnetization plateaus observed in TmB(4) and other rare-earth tetraborides. The simulations reproduce the 1/2 magnetization plateau by taking into account the dipole-dipole interaction. In addition, a narrow 2/3 magnetization step at low temperature is predicted in our simulation. The multi-step magnetization can be understood as the consequence of the competitions among the spin-exchange interaction, the dipole-dipole interaction, and the static magnetic energy.     

Monte Carlo Simulation of Kinesin Movement with a Lattice Model

Kinesin is a processive double-headed molecular motor that moves along a microtubule by taking about 8nm steps. It generally hydrolyzes one ATP molecule for taking each forward step. The processive movement of the kinesin molecular motors is numerically simulated with a lattice model. The motors are considered as Brownian particles and the ATPase processes of both heads are taken into account. The Monte Carlo simulation results agree well with recent experimental observations, especially on the relation of velocity versus ATP and ADP concentrations.     

Monte Carlo Study on Focus Properties of Portable Ultrabright Microfocus X-Ray Sources

The construct and electrode potential of emitting systems are very important for the portable ultrahigh brightness microfocus x-ray sources. The ratio of Dw/H （Dw is the diameter of Wehnelt grid aperture and H is the setting height of the cathode） and the grid bias are determinative parameters for the emission current and focus properties of an electron beam. The Monte Carlo method is used to numerically calculate the electron beam trajectories within a microfocus x-ray tube. The calculated results indicate that an optimum source system with the electron beam focal spot diameter smaller than 10μm at the anode, reasonable short focal length of about 76mm and the cathode emission current of more than 30mA can be attained under conditions of Dw/H = 5 and grid bias Vg = -0.6 k V.     

Monte Carlo Study of Griffiths Phase in Randomly Site Diluted Ising Magnetic System

We study the dynamics of spin autocorrelation function and the non-equilibrium decay of magnetization in the Criffiths phase of randomly site diluted Ising magnetic systems by Monte Carlo simulation. The results show that the relaxations are non-exponential for both two- and three-dimenslonal diluted Ising magnetic systems. At a long timescale, the asymptotic dynamics of the spin autocorrelation function C（t） behaves in the form of C（t） ～ exp{-A（lnt）^d/d-1}, expected by the theoretical prediction based on clustering arguments. Furthermore we show that the evolution of the spin autocorrelation function satisfies rather well a stretched exponential form of C（t） ～ exp{-（t/t）^β} even at quite small timescales.     

A Monte Carlo Study of Low-Energy Electron Transport in Liquid Water： Influence of the Rutherford Formula and the Mott Model

Influence of the Rutherford formula and the Mott model on the Monte Carlo simulation of low-energy electron （〈10 keV） transport in liquid water is investigated. One of the features of the constructed Monte Carlo code is the implementation of the new optical-data model from Emfietzoglou et al. in inelastic cross section based on the dielectric response theory. In addition, a novel mean elastic cross section by means of the Mott model is proposed to calculate the electron elastic scattering for high simulation efficiency. The systematical calculations of both the distribution of energy depositions and penetration parameter for low-energy electrons in liquid water are performed by using the Rutherford formula and the Mott model, respectively, for the elastic collisions in the simulations. The calculated results show that the effect of the two models is evident at energies below about 1 keV.     

Development and Application of Fast Monte Carlo Simulation Method for Beijing τ-Charm Factory (BTCF)

Based on the existing concept design and the expected performance of the detector at the future τ-charm factory,a fast Monte Carlo simulation program package using the BES software frame and data structure has been devcloped.Vadous important and/or interesting physics issues have been simulated by using this method,and a series of important physics results have been got.The design idea,the program frame,and the advantages of this method have been described.One can conclude that this Monte Carlo simulation method can greatly save CPU and disk space,and is specially suitable for BES users.The physics results obtained from the simulation demonstrate that it is essential to build the BTCF.     

Chaos and Monte Carlo Approximations of the Flip-Annihilation process

The flip-annihilation process is a random particle process with one-dimensional local interaction in discrete time, initially presented by one of us, namely Toom in 2004. Its components are enumerated by integer numbers and every component has two states, “minus” and “plus”. At every time step two transformations occur. The first one, called “flip”, independently turns every minus into plus with probability β. The second one, called “annihilation”, acts thus: whenever a plus is a left neighbor of a minus, both disappear with probability α independently from other components. What is interesting about this process is that it is ergodic for β>α/2 and non-ergodic for β<α ²/250. It is natural to conjecture that there is some transition curve, which we call the true curve and denote by , which separates the areas of ergodicity and non-ergodicity of this process from each other. The estimates, mentioned above, albeit rigorous, leave a large gap between them and the present article’s purpose is to obtain some closer, albeit non-rigorous, approximations of the true curve. We do it in two ways, one of which is a chaos approximation and the other is a Monte Carlo simulation. Thus we obtain two curves, which are much closer to each other than the rigorous estimations. Also we fill in, albeit only numerically, another shortcoming of the rigorous estimation β<α ²/250, namely that it leaves us uncertain whether the true curve has a zero or positive slope at the point α=β=0. Both approximate curves have a positive slope at α=0, as we hoped.     

Energy distribution of cosmic rays in the Earth’s atmosphere and avionic area using Monte Carlo codes

Cosmic rays cause significant damage to the electronic equipments of the aircrafts. In this paper, we have investigated the accumulation of the deposited energy of cosmic rays on the Earth’s atmosphere, especially in the aircraft area. In fact, if a high-energy neutron or proton interacts with a nanodevice having only a few atoms, this neutron or proton particle can change the nature of this device and destroy it. Our simulation based on Monte Carlo using Geant4 code shows that the deposited energy of neutron particles ranging between 200 MeV and 5 GeV are strongly concentrated in the region between 10 and 15 km from the sea level which is exactly the avionic area. However, the Bragg peak energy of proton particle is slightly localized above the avionic area.     

Study of the zero modes of the Faddeev–Popov operator in the maximal Abelian gauge

A study of the zero modes of the Faddeev–Popov operator in the maximal Abelian gauge is presented in the case of the gauge group SU(2)$SU\left(2\right)$ and for different Euclidean space–time dimensions. Explicit examples of classes of normalizable zero modes and corresponding gauge field configurations are constructed by taking into account two boundary conditions, namely: (i) the finite Euclidean Yang–Mills action, (ii) the finite Hilbert norm.