A Study of Surfactant-Induced Growth of Ag on Ag （111）

A new growth model is introduced to describe surfactant-induced growth of Ag on Ag （111） with realistic physical parameters. In this model, the A-S exchange mechanism is considered for the first time. Using the Monte Carlo simulations, the influence of exchange mechanism, surface temperature T, the exchange barrier Eεx, and the coverage of surfactant θM on the growth mode and morphology during multilayer film growth of Ag/Ag （111） are studied in detail Both the referenced value of surfactant coverage and the method to obtain perfect layer-by-layer film in surfactantinduced Ag/Ag （111） system are provided. Our simulation results are consistent with many experimental observations for surfactant-induced growth of Ag on Ag （111）.     

Coverage Dependence of Growth Mode in Heteroepitaxy of Ni on Cu（100） Surface

A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.     

Study of neutron response for two hybrid RPC setups using the GEANT4 MC simulation approach

The present article describes a detailed neutron simulation study in the energy range 10^-10 MeV to 1.0 GeV for two different RPC configurations. The simulation studies were taken by using the GEANT4 MC code. Aluminum was utilized on the GND and readout strips for the （a） Bakelite-based and （b） glass-based RPCs. For the former type of RPC setup the neutron sensitivity for the isotropic source was Sn = 2.702 × 10^-2 at En = 1.0 GeV, while for the latter type of RPC, the neutron sensitivity for the same source was evaluated as Sn = 4.049 × 10^-2 at En = 1.0 GeV. These results were further compared with the previous RPC configuration in which copper was used for ground and pickup pads. Additionally A1 was employed at （GND＋strips） of the phosphate glass RPC setup and compared with the copper-based phosphate glass RPC. Good agreement with sensitivity values was obtained with the current and previous simulation results.     

Monte Carlo simulation of the property of a scintillation bar in the multi-neutron correlation spectrometer

To perform a kinematically complete measurement of the dissociation reaction for neutron-rich nuclei, a multi-neutron correlation spectrometer is proposed at Peking University. A Monte Carlo simulation code based on GEANT4 is developed for a single scintillation bar which processes not only the energy deposition but also the light propagation in the scintillator and the light collection and conversion to signal at the end of the bar in a realistic way. The simulating method is described in detail in this paper, and the timing and position resolutions and detector efficiency are studied based on the simulation and compared with the experimental results. A new method of crosstalk rejection has been demonstrated to be important for the design of the whole spectrometer.     

Speckle intensity images of target based on Monte Carlo method

Speckle intensity in the detector plane is deduced in the free-space optical system and imaging system based on Van Cittert-Zemike theorem. The speckle intensity images of plane target and conical target are obtained by using the Monte Carlo method and measured experimentally. The results show that when the range extent of target is smaller, the speckle size along the same direction become longer, and the speckle size increase with increasing incident light wavelengths. The speckle size increases and the speckle intensity images of target is closer to the actual object when the aperture scale augments. These findings are useful to access the target information by speckle in laser radar systems.     

Theoretical study of molecular hydrogen and spiltover hydrogen storage on two-dimensional covalent-organic frameworks

Molecular hydrogen and spiltover hydrogen storages on five two-dimensional （2D） covalent-organic frameworks （COFs） （PPy-COF, TP-COF, BTP-COF, COF-18 A, and HHTP-DPB COF） are investigated using the grand canonical Monte Carlo （GCMC） simulations and the density functional theory （DFT）, respectively. The GCMC simulated results show that HHTP-DPB COF has the best performance for hydrogen storage, followed by BTP-COF, TP-COF, COF-18 A, and PPy-COE However, their adsorption amounts at room temperature are all too low to meet the uptake target set by US Department of Energy （US-DOE） and enable practical applications. The effects of pore size, surface area, and isosteric heat of hydrogen on adsorption amount are considered, which indicate that these three factors are all the important factors for determining the H2 adsorption amount. The chemisorptions of spiltover hydrogen atoms on these five COFs represented by the cluster models are investigated using the DFT method. The saturation cluster models are constructed by considering all possible adsorption sites for these cluster models. The average binding energy of a hydrogen atom and the saturation hydrogen storage density are calculated. The large average binding energy indicates that the spillover process may pro- ceed smoothly and reversibly. The saturation hydrogen storage density is much larger than the physisorption uptake of H2 molecules at 298 K and 100 bar （1 bar = 105 Pa）, and is close to or exceeds the 2010 US-DOE target of 6 wt% for hydrogen storage. This suggests that the hydrogen storage capacities of these COFs by spillover may be significantly enhanced. Thus 2D COFs studied in this paper are suitable hydrogen storage media by spillover.     

Finite Size Effect in Path Integral Monte Carlo Simulations of ^4He Systems

Path integral Monte Carlo （PIMC） simulations are a powerful computational method to study interacting quantum systems at finite temperatures. In this work, PIMC has been applied to study the finite size effect of the simulated systems of ^4He. We determine the energy as a function of temperature at saturated-vapor-pressure （SVP） conditions in the temperature range of T ∈ [1.0 K,4.0 K], and the equation of state （EOS） in the grmmd state For systems consisted of 32, 64 and 128 ^4He atoms, respectively, We find that the energy at SVP is influenced significantly by the size of the simulated system in the temperature range of T ∈ [2.1 K, 3.0 K] and the larger the system is, the better results are obtained in comparison with the experimental values; while the EOS appeared to be unrelated to it.     

Energy of Long-Lifetime Configurations in Zero-Temperature Dynamics

Using Monte Carlo method with zero-temperature dynamics, we investigate energy evolution of Ising spin configuration on a square lattice. The energies of some configurations exhibit long duration before those configurations reach the final state -- ground state or frozen stripe state. For ground-state dynamical realization, the duration occurs when the energy per spin is 4/L, where L is the lattice size. For stripe-state dynamical realization, the energy is slightly higher than 2/L when the duration appears in the last evolution stage. In addition, it is found that the average energy per spin in final state is approximately 2/3L.     

Strong and electromagnetic J/ψ and ψ（2S） decays into pion and kaon pairs

We discuss interference effects important for the form factors extraction in the vicinity of J/ψ and ψ（2S） resonances in combination with resonance parameters determination. The implementation to the Monte Carlo event generator PHOKHARA of the J/ψ and ψ（2S） contributions to the muon, pion and kaon pairs production associated with a photon at next-to-leading order is also described.     

Design and simulations for the detector based on DSSSD

The present paper describes the design and simulation results of a position-sensitive charged particle detector based on the Double Sided Silicon Strip Detector （DSSSD）. Also, the characteristics of the DSSSD and its testing result were are discussed. With the application of the DSSSD, the position-sensitive charged particle detector can not only give particle flux and energy spectra information and identify different types of charged particles, but also measure the location and angle of incident particles. As the detector can make multi-parameter measurements of charged particles, it is widely used in space detection and exploration missions, such as charged particle detection related to earthquakes, space environment monitoring and solar activity inspection.     

Reverse electric field Monte Carlo simulation for vector radiative transfer in the atmosphere

In this paper, a reverse electric field Monte Carlo （REMC） method is proposed to study the vector radiation transfer in the atmosphere. The REMC is based on tracing the multiply scattered electric field to simulate the vector transmitted radiance. The reflected intensities with different total optical depth values are obtained, which accord well with the results in the previous research. Stokes vector and the degree of polarization are numerically investigated. The simulation result shows that when the solar zenith angle is determined, the zenith angle of detector has two points, of which the degree of polarization does not change with the ground albedo and the optical depth. The two points change regularly with the solar zenith angle. Moreover, our REMC method can be applied to the vector radiative transfer in the atmosphere-ocean system.     

Influence of thyroid volume reduction on absorbed dose in 131I therapy studied by using Geant4 Monte Carlo simulation

A simulation study has been performed to quantify the effect of volume reduction on the thyroid absorbed dose per decay and to investigate the variation of energy deposition per decay due to β- and γ-activity of 131I with volume/mass of thyroid, for water, ICRP- and ICRU-soft tissue taken as thyroid material. A Monte Carlo model of the thyroid, in the Geant4 radiation transport simulation toolkit was constructed to compute the β- and γ-absorbed dose in the simulated thyroid phantom for various values of its volume. The effect of the size and shape of the thyroid on energy deposition per decay has also been studied by using spherical, ellipsoidal and cylindrical models for the thyroid and varying its volume in 1-25 cm3 range. The relative differences of Geant4 results for different models with each other and MCNP results lie well below 1.870%. The maximum relative difference among the Geant4 estimated results for water with ICRP and ICRU soft tissues is not more than 0.225%. S-values for ellipsoidal, spherical and cylindrical thyroid models were estimated and the relative difference with published results lies within 3.095%. The absorbed fraction values for beta particles show a good agreement with published values within 2.105% deviation. The Geant4 based simulation results of absorbed fractions for gammas again show a good agreement with the corresponding MCNP and EGS4 results （±6.667%） but have 29.032% higher values than that of MIRD calculated values. Consistent with previous studies, the reduction of the thyroid volume is found to have a substantial effect on the absorbed dose. Geant4 simulations confirm dose dependence on the volume/mass of thyroid in agreement with MCNP and EGS4 computed values but are substantially different from MIRD8 data. Therefore, inclusion of size/mass dependence is indicated for 131I radiotherapy of the thyroid.     

Microwave propagation with the gas breakdown

By combining the microwave propagation theory and the gas breakdown theory, the microwave propagation with the gas breakdown is analyzed theoretically. Particle-in-cell/Monte Carlo collision （PIC/MCC） simulations are carried out to verify the theoretical results. Based on this theoretical method, the breakdown phenomenon of the pulse microwave is analyzed. The results show that the product values of the initial electron density and the propagation length are the criterion to distinguish the pulse peak decline breakdown and the pulse width reduction breakdown. Furthermore, the energy transmission is also studied, which shows that the total output energy is approximately independent of the input electric field if the electric field is not extremely large.     

Supersolid Phase in One-Dimensional Hard-Core Boson Hubbard Model with a Superlattice Potential

The ground state of the one-dimensional hard-core boson Hubbard model with a superlattice potential is studied by quantum Monte Carlo methods. We demonstrate that besides the CDW phase and the Mott insulator phase, the supersolid phase emerges due to the presence of the superlattice potential, which reflects the competition with the hopping term. We also study the densities of sublattices and have a clear idea about the distribution of the bosons on the lattice.     

Potential-Decomposition Strategy in Markov Chain Monte Carlo Sampling Algorithms

We introduce the potential-decomposition strategy （PDS）, which can be used in Markov chain Monte Carlo sampling algorithms. PDS can be designed to make particles move in a modified potential that favors diffusion in phase space, then, by rejecting some trial samples, the target distributions can be sampled in an unbiased manner. Furthermore, if the accepted trial samples are insumcient, they can be recycled as initial states to form more unbiased samples. This strategy can greatly improve efficiency when the original potential has multiple metastable states separated by large barriers. We apply PDS to the 2d Ising model and a double-well potential model with a large barrier, demonstrating in these two representative examples that convergence is accelerated by orders of magnitude.     

Status and accuracy of the Monte Carlo generators for luminosity measurements

The status and accuracy of the precision Monte Carlo generators used for luminosity measurements at flavour factories is reviewed. It is shown that, thanks to a considerable, long-term effort in tuned comparisons between the predictions of independent programs, as well as in the validation of the generators against the presently available calculations of the next-to-next-to-leading order QED corrections to Bhabha scattering, the theoretical accuracy reached by the most precise tools is of about one per mille. This error estimate is valid for realistic experimental cuts, appears to be quite robust and is already sufficient for very accurate luminosity measurements. However, recent progress and possible advances to further improve it are also discussed.     

Quantum Monte Carlo Study of Frustrated Spin Chain under External Field

In this paper a stochastic series expansion quantum Monte Carlo algorithm is used to study a frustrated spin chain with diagonal next-nearest-neighbor interactions. The detailed balance conditions are carefully analyzed to improve the efficiency of simulation process. As an application of this algothrim, the total magnetization, the static structure factor and spin-stiffness are calculated for a certain set of system parameters as a function of external field strength.     

A simulation study of the e+/e－ sensitivity of a low resistive phosphate glasselectrode 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.     

News on PHOTOS Monte Carlo： γ＊→π^＋π^-（γ） and K^±→π^＋π^-e^±υ（γ）

PHOTOS Monte Carlo is widely used for simulating QED effects in decay of intermediate particles and resonances. It can be easily connected to other main process generators. In this paper we consider decaying processes γ＊→π^＋π^-（γ） and K^±→π^＋π^-e^±υ（γ） in the framework of Scalar QED. These two processes are interesting not only for the technical aspect of PHOTOS Monte Carlo, but also for precision measurement of αED（Mz）, g--2, as well as ππ scattering lengths.     

Estimation of pulsed laser-induced single event transient in a partially depleted silicon-on-insulator 0.18-μm MOSFET

In this paper, we investigate the single event transient(SET) occurring in partially depleted silicon-on-insulator(PDSOI) metal–oxide–semiconductor(MOS) devices irradiated by pulsed laser beams. Transient signal characteristics of a 0.18-μm single MOS device, such as SET pulse width, pulse maximum, and collected charge, are measured and analyzed at wafer level. We analyze in detail the influences of supply voltage and pulse energy on the SET characteristics of the device under test(DUT). The dependences of SET characteristics on drain-induced barrier lowering(DIBL) and the parasitic bipolar junction transistor(PBJT) are also discussed. These results provide a guide for radiation-hardened deep sub-micrometer PDSOI technology for space electronics applications.