利用Bayesian-MCMC方法从雷达回波反演海洋波导

应用贝叶斯-蒙特卡罗（Bayesian-MCMC）方法将海洋波导参数的先验信息描述为先验概率密度,结合雷达回波资料（电磁波传播损耗）,得到待反演海洋波导参数的后验概率密度,用马尔可夫链蒙特卡罗（MCMC）-Gibbs采样器采样后验概率密度分布,并用样本最大似然估计值作为对海洋波导参数分布的估计.数值实验结果表明,该方法对先验信息进行了有效利用,反演精度高于遗传算法的反演精度.该方法较为充分利用先验信息,得到解的概率分布,即解的不确定性分析,这在实际应用中有一定的参考价值. 关键词： 波导 电磁波传播损耗 贝叶斯-蒙特卡罗 概率分布     

Critical properties of Rule 22 elementary cellular automata

The one-dimensional elementary cellular automaton Rule 22 is studied by means of Monte Carlo simulation on the dedicated K2 high-speed computer. If one considers random initialization with probability p for one-initialization per site, it is shown that the system behaves like a normal one-dimensional statistical ensemble with critical points atp=0 andp=1. Critical slowing down is exhibited, with a dynamical exponent of 1.0. The standard initialization ofp=0.5 is too far away from the critical point to allow similar observations.     

模拟回火马尔可夫链蒙特卡罗全波形分析方法

针对传统的全波形分析方法不能快速自动处理全波形数据的缺点,提出了一种模拟回火马尔可夫链蒙特卡罗全波形分析法,用于求解全波形数据中的波峰数和峰值位置等参量.该方法采用Metropolis更新策略求解波峰数量和噪声两个参量,以达到快速求解的目的;而峰值位置和波峰幅值则采用改进的模拟回火策略求解,通过添加的主动干预回火步骤实现对参量更新过程的有效探测,以满足对速度或运算收敛性的要求.模拟回火马尔可夫链蒙特卡罗全波形分析方法以马尔可夫算法为基础,仍保持马氏链的收敛性,从而保证本方法具有良好的鲁棒性,实现对全波形数据的自动化处理.     

On the formation of spin dependent double layer structure in astrophysical plasmas

Double layers (DLs) are nonlinear structures, and are thought to be responsible for particle acceleration in laboratory plasmas and astrophysical plasmas. In this paper, the existence of spin dependent DLs structure is explored using separate spin evolution quantum hydrodynamic model. Based on reductive perturbation method, we derived an extended Korteweg–de Vries (eKdV) equation to demonstrate the existence and nature of DLs. We found that spin polarization significantly enhanced the amplitude of the electrostatic potential associated with DLs. Further, spin polarization also increases the depth and width of the Sagdeev potential. It is noted that the contribution of Bohm potential effect is essential for the formation of DLs structure. Our results may be helpful to explain the phenomenon of particle acceleration in dense astrophysical environments specifically in a white dwarf.     

蒙特卡罗模拟单阻止柱双散射体质子束流扩展方法

束流扩展是研究质子治疗大体积深度肿瘤的基础,利用Monte Carlo多粒子模拟软件Fluka2006.3b.10,计算质子经过单阻止柱双散射体后束流横向分布特性,与相同条件下由Highland公式经解析计算相比较.结果表明:在束流分布均匀性大于95%的范围内,Fluka计算得到的质子通量小于由Highland公式经解析计算得到的结果,膜薄的条件下二者符合好于膜厚的条件.束流扩展半径随着两个散射膜厚度的增加而增大,束流利用率先增加后减小,最大值是16.9%.     

Evidence for a second,local contact condition for a planar electric double layer containing size and valency asymmetric ions

A generalised form of a local contact condition for the charge profile in a primitive model planar double layer [Bhuiyan, Outhwaite, and Henderson, Mol. Phys. 107, 343 (2009)] at low electrode charge is examined for completely asymmetric, binary electrolytes. The cation and anion sizes are taken to be different from each other with the valencies being 2⁺:1^? or 1⁺:2^?, while the electrode surface charge density is varied from being negative through zero to being positive. Monte Carlo simulation data obtained for such double layer systems at varying ionic radius ratios and electrolyte concentrations suggest the generalised contact relation to be valid at low charge on the electrode.     

降雨融化层后向散射的蒙特卡罗仿真

云层上端冰雪粒子的融化形成了降雨融化层，随着粒子的下落融化过程开始，因此融化层的微观特性在垂直方向上是连续变化的.建立了降雨融化层的仿真模型，应用更为切合实际的三层球形粒子代替了融化层中的粒子.根据降雨过程中雨滴的尺寸分布推导了融化层中融化粒子的尺寸分布.根据Mie理论，计算了降雨率小于12.5 mm/h时，降雨融化层对5，10，35，94 GHz电磁波的雷达反射率和特征衰减因子的垂直廓线.计算结果表明，当电磁波频率高于20 GHz时，将无法观测到“雷达亮带”，这与实验结果相符.根据辐射传输理论，应用蒙特卡罗方法计算了垂直方向微观特性连续变化的降雨融化层对不同频率电磁波的反射率，比较了两种不同尺寸分布(Gamma分布和Marshall-Palmer分布)融化层反射率的差别，这为利用高频电磁波对降雨融化层进行遥感提供了理论和数值依据. 关键词： 降雨融化层 蒙特卡罗方法 后向散射 电磁波     

分子核间距对非时序双电离的影响

利用半经典再散射模型研究了激光强度在一定范围内时分子核间距对非时序双电离过程的影响.分别计算了非时序双电离率、两电子的电离能、两电子的动量相关、Coulomb和激光场的复合势随分子核间距的变化关系.研究表明,分子核间距在1.0—6.0 a.u.范围内时,非时序双电离率和两电子动量和为零的双电离事件数随着分子核间距的增大而增大.当分子核间距继续增大(大于6.0 a.u.)时,非时序双电离率和两电子动量和为零的双电离事件数却减小.     

Monte Carlo study of the critical and the compensation behaviors of the double perovskite Sr2CrIrO6

The phase diagrams and magnetic properties of double perovskite $S{r}_{2}CrIr{O}_6$ have been studied by using Monte Carlo simulation based on the heat bath algorithm. The ground-state diagrams of the compound $S{r}_{2}CrIr{O}_6$ have been calculated for different combinations of system parameters. The diagrams obtained are very rich and they give an idea of all the most stable configurations. The effects of the exchange interactions and the crystal field on the phase diagrams and magnetic properties of the system have been examined. A number of interesting phenomena have been observed such as the compensation temperature, the first and second order phase transitions, the critical triple point and the terminal critical point.     

Modeling of Terahertz Radiation from InSb and InAs

In this paper the THz radiation dynamics of InSb and InAs, two typical narrow band semiconductors, was investigated using the ensemble Monte Carlo method. Our simulations indicated that the single mechanism of current surge only can result in small difference of THz emission efficiency for InSb and InAs. The great advantage of InAs over InSb in THz emission efficiency that was found in a published work is possibly due to the mechanism of optical rectification. In addition, under low excitation level, we found the emission efficiency of InSb is advantage over that of InAs, but under high excitation level, the result is reversed. On the other hand, through the Fourier transforms of temporal THz waveforms we found that the main frequency of THz pulses from InAs is always higher than that of InSb.     

Spin frustration-related behavior of mixed spin diamond chain

The magnetization of mixed spin-1/2, 1 Ising diamond chain with antiferromagnetic spin exchange interaction has been investigated by using the Monte Carlo simulation. Magnetization behaviors with different antiferromagnetic interactions are simulated under magnetic field at low temperature. The behaviors are studied in two different cases. In one case the absolute value of the interaction between two spin-1/2 atoms is the largest and in the other case it is not the largest. Multiple plateaus of magnetization M, such as plateaus at M=±1, M=0 and M=±0.566 are obtained. The plateaus at M=±0.566 correspond to two metastable states. The M=0 plateau may disappear at some specific values of spin exchange interactions. The spin configurations were checked when the plateaus appear. The dimmer and trimmer states are found.     

Electrostatic double layers in a warm negative ion plasma with nonextensive electrons

The electrostatic double layer (DL) structures are studied in negative ion plasma with nonextensive electrons q-distribution. The extended Korteweg–de Vries (EKdV) equation is derived using a reductive perturbation method. It is found that both fast (compressive) and slow (rarefactive) ion acoustic (IA) DLs can propagate in such type of plasmas. The effects of various plasma physical parameters; such as nonextensivity of electrons, presence of negative ions, temperature of both positive and negative ions and different mass ratios of positive to negative ions on the formation of DL structures are discussed in detail with numerical illustrations.     

Effect of ionic size on the structure of spherical double layers: a Monte Carlo simulation and density functional theory study

The effect of ionic size on the diffuse layer characteristics of a spherical double layer is studied using Monte Carlo simulation and density functional theory within the restricted primitive model. The macroion is modelled as an impenetrable charged hard sphere carrying a uniform surface charge density, surrounded by the small ions represented as charged hard spheres and the solvent is taken as a dielectric continuum. The density functional theory uses a partially perturbative scheme, where the hard sphere contribution to the one particle correlation function is evaluated using weighted density approximation and the ionic interactions are calculated using a second-order functional Taylor expansion with respect to a bulk electrolyte. The Monte Carlo simulations have been performed in the canonical ensemble. The detailed comparison is made in terms of zeta potentials for a wide range of physical conditions including different ionic diameters. The zeta potentials show a maximum or a minimum with respect to the polyion surface charge density for a divalent counterion. The ionic distribution profiles show considerable variations with the concentration of the electrolyte, the valency of the ions constituting the electrolyte, and the ionic size. This model study shows clear manipulations of ionic size and charge correlations in dictating the overall structure of the diffuse layer.     

Monte Carlo simulation of overdamped motion in a double well potential under the influence of coloured noise

We present a Monte Carlo simulation algorithm for evaluating the stationary probability and the mean and the variance of first passage times in any dynamical system under the influence of additive coloured Gaussian and Marcovian noise (Ornstein-Uhlenbeck process). Our algorithm generates the Ornstein-Uhlenbeck process by a superposition of a finite number of random telegraph processes. We obtain our results from a direct evaluation of the trajectories. We apply our method to the overdamped motion of a particle in a double well potential. We compare our simulation results with various analytic approximations for the stationary probability and the mean first passage times.     

Application of a probabilistic double-fibre structure model to diffusion-weighted MR images of the human brain

A Markov chain Monte Carlo (MCMC) algorithm has been reported which is capable of determining the probabilistic orientation of two-fibre populations from high angular resolution diffusion-weighted data (HARDI). We present and critically discuss the application of this algorithm to in vivo human datasets acquired in clinically realistic times. We show that by appropriate model selection areas of multiple fibre populations can be identified that correspond with those predicted from known anatomy. Quantitative maps of fibre orientation probability are derived and shown for one- and two-fibre models of neural architecture. Fibre crossings in the pons, the internal capsule and the corona radiata are shown. In addition, we demonstrate that the relative proportion of anisotropic signal may be a more appropriate measure of anisotropy than summary measures derived from the tensor model such as fractional anisotropy in areas with multi-fibre populations.     

Measuring preferential attachment in growing networks with missing-timelines using Markov chain Monte Carlo

Preferential attachment is widely recognised as the principal driving force behind the evolution of many growing networks, and measuring the extent to which it occurs during the growth of a network is important for explaining its overall structure. Conventional methods require that the timeline of a growing network is known, that is, the order in which the nodes of the network appeared in time is available. But growing network datasets are commonly accompanied by missing-timelines, in which instance the order of the nodes in time cannot be readily ascertained from the data. To address this shortcoming, we propose a Markov chain Monte Carlo algorithm for measuring preferential attachment in growing networks with missing-timelines. Key to our approach is that any growing network model gives rise to a probability distribution over the space of networks. This enables a growing network model to be fitted to a growing network dataset with missing-timeline, allowing not only for the prevalence of preferential attachment to be estimated as a model parameter, but the timeline also. Parameter estimation is achieved by implementing a novel Metropolis–Hastings sampling scheme for updating both the preferential attachment parameter and timeline. A simulation study demonstrates that our method accurately measures the occurrence of preferential attachment in networks generated according to the underlying model. What is more, our approach is illustrated on a small sub-network of the United States patent citation network. Since the timeline for this example is in fact known, we are able to validate our approach against the conventional methods, showing that they give mutually consistent estimates.     

Complex formation between a nanoparticle and a weak polyelectrolyte chain: Monte Carlo simulations

Understanding the complexation processes between nanoparticles and polyelectrolytes is an essential aspect in many branches of nanotechnology, nanoscience, chemistry, and biology to describe processes such as nanoparticle stabilization/destabilization and dispersion, water treatment, microencapsulation, complexation with biomolecules for example, and evolution of the interface of many natural and synthetic systems. In view of the complexity of such processes, applications are often based on empirical or semiempirical observations rather than on predictions based on theoretical or analytical models. In this study, the complex formation between an isolated weak polyelectrolyte and an oppositely charged nanoparticle is investigated using Monte Carlo simulations with screened Coulomb potentials in the grand canonical ensemble. The roles of the nanoparticle surface charge density , solution pH and ionic concentration C_i are systematically investigated. The phase diagrams of complex conformations are also presented. It is shown that the polyelectrolyte conformation at the surface of the nanoparticle is controlled by the attractive interactions with the nanoparticle but also by the repulsive interactions between the monomers. To bridge the gap with experiments titration curves are calculated. We clearly demonstrate that an oppositely charged nanoparticle can significantly modify the acid/base properties of a weak polyelectrolyte.     

Effective-field and Monte Carlo studies of mixed spin-2 and spin-1/2 Ising diamond chain

The magnetic properties of a mixed spin-2 and spin-1/2 ferromagnetic diamond chain are studied by effective-field theory and Monte Carlo(MC) simulation based on the Ising model.The temperature dependences of magnetization,magnetic susceptibility,internal energy,and specific heat are studied,respectively.The exchange interaction dependences of magnetization and the critical temperature are calculated by MC simulation.The changes of magnetization depending on the field increasing and then the field decreasing under steady-static conditions are also given.     

Measurement of the dead layer thickness in a p-type point contact germanium detector

A 994 g mass p-type PCGe detector has been deployed during the first phase of the China Dark matter EXperiment, aiming at direct searches for light weakly interacting massive particles. Measuring the thickness of the dead layer of a p-type germanium detector is an issue of major importance since it determines the fiducial mass of the detector. This work reports a method using an uncollimated ~(133)Ba source to determine the dead layer thickness. The experimental design, data analysis and Monte Carlo simulation processes, as well as the statistical and systematic uncertainties are described. A dead layer thickness of 1.02 mm was obtained based on a comparison between the experimental data and the simulated results.