Neutron- and muon-induced background in underground physics experiments

Background induced by neutrons in deep underground laboratories is a critical issue for all experiments looking for rare events, such as dark matter interactions or neutrinoless ββ decay. Neutrons can be produced either by natural radioactivity, via spontaneous fission or (α, n) reactions, or by interactions initiated by high-energy cosmic rays. In all underground experiments, Monte Carlo simulations of neutron background play a crucial role for the evaluation of the total background rate and for the optimization of rejection strategies. The Monte Carlo methods that are commonly employed to evaluate neutron-induced background and to optimize the experimental setup, are reviewed and discussed. Focus is given to the issue of reliability of Monte Carlo background estimates. We dedicate this work to the memory of our friend and colleague Nicola Ferrari, who prematurely passed away in July 2006.     

Stochastic Order and Generalized Weighted Mean Invariance

In this paper, we present order invariance theoretical results for weighted quasi-arithmetic means of a monotonic series of numbers. The quasi-arithmetic mean, or Kolmogorov–Nagumo mean, generalizes the classical mean and appears in many disciplines, from information theory to physics, from economics to traffic flow. Stochastic orders are defined on weights (or equivalently, discrete probability distributions). They were introduced to study risk in economics and decision theory, and recently have found utility in Monte Carlo techniques and in image processing. We show in this paper that, if two distributions of weights are ordered under first stochastic order, then for any monotonic series of numbers their weighted quasi-arithmetic means share the same order. This means for instance that arithmetic and harmonic mean for two different distributions of weights always have to be aligned if the weights are stochastically ordered, this is, either both means increase or both decrease. We explore the invariance properties when convex (concave) functions define both the quasi-arithmetic mean and the series of numbers, we show its relationship with increasing concave order and increasing convex order, and we observe the important role played by a new defined mirror property of stochastic orders. We also give some applications to entropy and cross-entropy and present an example of multiple importance sampling Monte Carlo technique that illustrates the usefulness and transversality of our approach. Invariance theorems are useful when a system is represented by a set of quasi-arithmetic means and we want to change the distribution of weights so that all means evolve in the same direction.     

Stochastic resonance for adhesion of membranes with active stickers

The behavior of two membranes that interact by active adhesion molecules or stickers is studied theoretically using mean-field theory and Monte Carlo simulations. The stickers are anchored in one of the membranes and undergo conformational transitions between on and off states. In their on states, the stickers can bind to ligands that are anchored in the other membrane. The transitions between the on and off states arise from the coupling of the stickers to some active, energy-releasing process, which keeps the system out of equilibrium. As one varies the transition rates of this active process, the membrane separation undergoes a stochastic resonance: this separation is maximal at intermediate rates of the sticker transitions and considerably smaller both at high and at low transition rates. This implies that the effective, fluctuation-induced repulsion between the membranes contains a rate-dependent contribution that arises from the switching of the active stickers.     

Distributions of switching times of single-domain particles using a time quantified Monte Carlo method

Using a time quantified Monte Carlo scheme we performed simulations of the switching time distribution of single mono-domain particles in the Stoner–Wohlfarth approximation. We considered uniaxial anisotropy and different conditions for the external applied field. The results obtained show the switching time distribution can be well described by two relaxation times, either when the applied field is parallel to the easy axis or for an oblique external field and a larger damping constant. We found that in the low barrier limit these relaxation times are in very good agreement with analytical results obtained from solutions of the Fokker–Planck equation related to this problem. When the damping is small and the applied field is oblique the shape of the distribution curves shows several peaks and resonance effects.     

Differences between fixed time step and kinetic Monte Carlo methods for biased diffusion

We consider biased diffusion in a one-dimensional lattice and compare results obtained with fixed time step and kinetic Monte Carlo methods. Spurious dispersion and particle position correlation appear with the fixed time step Monte Carlo approach. The mentioned correlation increases with time. We demonstrate that the correct results, that correspond to a time step that tends to zero, are obtained using the kinetic Monte Carlo method. The conclusions also apply to biased diffusion in two or more dimensions and to random deposition.     

Simulation study of random sequential adsorption of mixtures on a triangular lattice

Random sequential adsorption of binary mixtures of extended objects on a two-dimensional triangular lattice is studied numerically by means of Monte Carlo simulations. The depositing objects are formed by self-avoiding random walks on the lattice. We concentrate here on the influence of the symmetry properties of the shapes on the kinetics of the deposition processes in two-component mixtures. Approach to the jamming limit in the case of mixtures is found to be exponential, of the form: θ(t) ∼ θ_jam - Δθ exp(- t/σ), and the values of the parameter σ are determined by the order of symmetry of the less symmetric object in the mixture. Depending on the local geometry of the objects making the mixture, jamming coverage of a mixture can be either greater than both single-component jamming coverages or it can be in between these values. Results of the simulations for various fractional concentrations of the objects in the mixture are also presented.     

JMCT与子通道程序耦合方法研究及验证

蒙特卡罗与热工水力的耦合计算是目前反应堆数值模拟的重要研究方向,在蒙特卡罗方法连续能量点截面的基础上结合热工程序的温度反馈,反应堆中子计算的准确性得到大幅提高。为了提高计算精度,堆芯模型分辨率也需进一步提高,相比于组件均匀化模型,pin-by-pin的建模方式能够获得更好的结果。利用蒙特卡罗程序JMCT与子通道程序COBRA-EN实现了蒙特卡罗-热工的内耦合,内耦合方式通过内存进行数据传递,其计算效率及安全性均优于外耦合方法。随后利用NURISP项目迷你堆的pin-by-pin模型对耦合程序进行验证。计算结果与同类耦合程序相似,验证了程序的准确性。同时,对耦合过程的收敛性问题进行了初步分析。     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.     

Phase transition of compartmentalized surface models

Two types of surface models have been investigated by Monte Carlo simulations on triangulated spheres with compartmentalized domains. Both models are found to undergo a first-order collapsing transition and a first-order surface fluctuation transition. The first model is a fluid surface one. The vertices can freely diffuse only inside the compartments, and they are prohibited from the free diffusion over the surface due to the domain boundaries. The second is a skeleton model. The surface shape of the skeleton model is maintained only by the domain boundaries, which are linear chains with rigid junctions. Therefore, we can conclude that the first-order transitions occur independent of whether the shape of surface is mechanically maintained by the skeleton (=the domain boundary) or by the surface itself.     

Multilevel coarse graining and nano-pattern discovery in many particle stochastic systems

In this work we propose a hierarchy of Markov chain Monte Carlo methods for sampling equilibrium properties of stochastic lattice systems with competing short and long range interactions. Each Monte Carlo step is composed by two or more sub-steps efficiently coupling coarse and finer state spaces. The method can be designed to sample the exact or controlled-error approximations of the target distribution, providing information on levels of different resolutions, as well as at the microscopic level. In both strategies the method achieves significant reduction of the computational cost compared to conventional Markov chain Monte Carlo methods. Applications in phase transition and pattern formation problems confirm the efficiency of the proposed methods.     

Positron backscattering from solid targets: Modeling of scattering processes via various approaches

Monte Carlo simulation of 1–4 keV positron backscattering from semi-infinite solid targets ranging from Be (z = 4) to Au (z = 79) with normal angle of incidence is here reported. In our study, the elastic and inelastic scattering cross sections are modeled by using various approaches based on either a classical or a quantum mechanical treatment. Calculations of positron backscattering coefficient are then reported for the solid targets of interest. The results obtained show a fairly good agreement with the data available in the literature. The dependence of the positron backscattering coefficient versus the atomic number of the solid target of interest has been investigated. In this respect, polynomial functions are proposed which does not require any recourse to Monte Carlo calculations.     

An algorithm to determine cirrus properties from analysis of multiple-scattering influence on lidar signals

A Monte Carlo-based evaluation of the multiple-scattering influence on ground-based Raman lidar measurements is presented. The lidar returns from cirrus clouds are analyzed in order to evaluate vertical profiles of the extinction and backscattering coefficients. Results show that for the typical cirrus cloud, the presence of the multiple scattering can lead to an underestimation of the extinction coefficient by as large as 200% whereas the backscattering coefficient is almost unaffected for the Raman lidar technique. An algorithm to select one or a set of phase functions which fit to the lidar data is also presented. It is an iterative procedure based on Monte Carlo scattering simulation. By comparison of the experimental value of the lidar ratio, corrected for the multiple scattering influence, and the phase function used in the Monte Carlo simulation, one can determine a suitable phase function. The validity and sensitivity of the algorithm have been demonstrated by applying it to simulated cases. The application to some real cases indicates that our procedure allows for the establishing of a practical scattering model for the cirrus clouds.     

Shape transformation transitions of a tethered surface model

A surface model of Nambu and Goto is studied statistical mechanically by using the canonical Monte Carlo simulation technique on a spherical meshwork. The model is defined by the area energy term and a one-dimensional bending energy term in the Hamiltonian. We find that the model has a large variety of phases; the spherical phase, the planar phase, the long linear phase, the short linear phase, the wormlike phase, and the collapsed phase. Almost all two neighboring phases are separated by discontinuous transitions. It is also remarkable that no surface fluctuation can be seen in the surfaces both in the spherical phase and in the planar phase.     

基于禁忌搜索的车联网蒙特卡洛定位算法

在蒙特卡洛定位算法中引入禁忌搜索算法以提高车联网中快速定位的性能。自组织车联网高速移动的车辆和快速变化的网络拓扑结构,使用传统的蒙特卡洛定位算法,不能迅速的收敛位置信息。在滤波阶段引入禁忌搜索算法对传统蒙特卡洛定位算法进行改进, 优化滤波排除可能性较小的位置点,获得近似最优估计位置采样集。仿真结果表明,改进后的算法在样本采集数、计算时间、定位精度等方面有了显著提升,改进后的算法能更好地解决车联网的定位问题。     

Monte Carlo Simulation of Optical Properties of Wake Bubbles

Based on Mie scattering theory and the theory of multiple light scattering, the light scattering properties of air bubbles in a wake are analysed by Monte Carlo simulation. The results show that backscattering is enhanced obviously due to the existence of bubbles, especially with the increase of bubble density, and that it is feasible to use the Monte Carlo method to study the properties of light scattering by air bubbles.     

Cold collisions in strong laser fields: partial wave analysis of magnesium collisions

We have developed Monte Carlo wave function simulation schemes to study cold collisions between magnesium atoms in a strong red-detuned laser field. In order to address the strong-field problem, we extend the Monte Carlo wave function framework to include the partial wave structure of the three-dimensional system. The average heating rate due to radiative collisions is calculated with two different simulation schemes which are described in detail. We show that the results of the two methods agree and give estimates for the radiative collision heating rate for ²⁴Mg atoms in a magneto-optical trap based on the ¹S₀–¹P₁ atomic laser cooling transition.     

Investigation of nucleation and grain growth in 2-dimensional systems by using generalized Monte Carlo simulations

In this study, nucleation and grain growth was studied by using 2-dimensional generalized Monte Carlo simulations and experiments. As an attempt to improve the JMAK model, we proposed a new differential equation to be able to model nucleation and growth phenomena using nonextensive thermostatistics. One of the reasons that we would like to perform generalized Monte Carlo simulations in studying of nucleation and grain growth phenomena is that the generalized Monte Carlo algorithm was shown to be more effective than the standard Monte Carlo algorithm and also than the standard Molecular Dynamic algorithm in locating the minimum energy configuration. Therefore, for a given temperature, the fact that a configuration of the system with lower energy could be obtained by using the generalized Monte Carlo simulation means that a different textural configuration of grain growth could be also expected. In this respect, it is possible to say that the nonextensive statistics might be an appropriate tool in studying of nucleation and growth phenomena.     

应用似然函数比以及Fourier变换计算置信度的解析方法

新粒子寻找的实验结果要通过计算置信度来解释,置信度的定义或是基于新粒子与本底并存的假设,或是基于唯有本底存在的假设.通常的计算方法是产生大量的toy Monte Carlo实验,这些实验按照某种定义的估计值的大小以顺序排列,然后将真正观测到的实验的估计值与之相比较,从而得到置信度.本文则介绍一种新的计算方法,通过定义似然函数比为实验的估计值,并应用Fourier变换,用解析的方法计算出置信度,与toy Monte Carlo方法相比,解析方法可极大地提高计算的速度和精度.