关于负超几何分布期望与方差算法的注记

研究了非还原取样模型中负超几何随机变量的联合分布,得到了若干有用的推论.据此给出了负超几何分布的期望和方差的一种分解算法.     

241Am-Be中子参考辐射的蒙特卡罗模拟研究

研究用于校准场所中子剂量监测仪表的241Am-Be中子参考辐射场计量特性。采用蒙特卡罗方法模拟了空气自由中子参考辐射(FRNR),GB/T 14055规定的最小尺寸中子参考辐射(SRNR)和实际中子参考辐射(ARNR)中不同检验点处中子周围剂量当量率、散射中子占比和能谱分布特征。研究结果表明,空气对FRNR中的剂量率和能谱分布影响小,近似为理想中子参考辐射;采用5%含硼聚乙烯作屏蔽的最小尺寸SRNR可减少热中子,降低散射中子占比,影锥法不适用于小尺寸中子参考辐射中对散射中子的修正;ARNR中的散射中子更少、占比更低,影锥法所得散射中子占比与理论值基本一致。     

均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

Monte Carlo simulations of magnetic particle suspensions with a simple assessment method for the particle overlap between magnetic spheroids

We have developed a simple assessment method for the overlap between spheroidal particles, which neither requires the complex manipulation of vectors and matrices that is indispensable in the ordinary methods, nor is based on a model potential. Moreover, we have developed an evaluation method for the interaction energy arising from the overlap of the steric layer coating spheroidal particles. This is based on a sphere-connected particle model, but some modifications are introduced in order to express an appropriate repulsive interaction energy at the deepest overlapping position. We have investigated the phase change in a magnetic spheroidal particle suspension for a two-dimensional system by means of Monte Carlo simulations. In the case of no external magnetic field, if the magnetic particle-particle interaction is sufficiently strong to favour cluster formation, long raft-like clusters tend to be formed in a dilute situation. With decreasing values of area fraction, a chain-like structure in a dense situation transforms into a raft-like structure within a narrow range of the particle area fraction. Similarly, the raft-like clusters are preferred in a weak applied magnetic field, but an increase in the field strength induces a phase change from a raft-like into a chain-like structure.

Highlights of the present paper:

1. A simple assessment method has been proposed for the overlap between two spheroidal particles.

2. The particle overlap assessment is free from a complex mathematical manipulation regarding vectors and matrices.

3. A modified sphere-connected model has been proposed in order to more accurately evaluate a repulsive interaction due to the overlap of the steric layers coating spheroidal particles.

4. 2D Monte Carlo simulations have been performed to elucidate the phenomenon of a phase change by magnetic spheroidal particles on a material plane surface.

5. A phase change between a raft-like and a chain-like aggregate structure is able to be controlled by the area fraction of particles and an external magnetic field.

     

Prediction of branch on branch and topological characteristics of low‐density polyethylene polymerization by a novel stochastic approach

A novel approach using Monte Carlo method applied to simulation of low‐density polyethylene (LDPE) polymerization in tubular reactor showing topological characteristics, and the comprehensive kinetic mechanism has been taken into consideration. The results show the precise details of the structure of a chain in the three levels of the backbone, the main branches, and branches on branch. The chain types include dead polymer, dead polymer with unsaturated end, and live polymer with primary radical, secondary radical, and tertiary radical. In this work, the branches on branch were identified in terms of number, length, and position of the branch. Sixty percent of branches on branch are 1 to 5 carbons long, and the longest branch on branch is about 50 carbons. Thus, this study provides a tool for more accurately mapping the polymer chains architecture, superior to determine the number, and position of long‐ and short‐chain branches in past researches. Finally, this approach will advance the prediction of microstructure‐related properties of polymer one step further.     

Molecular scale simulations on thermoset polymers: A review

This article reviews the field of molecular simulations of thermoset polymers. This class of polymers is of interest in applications ranging from structural components for aerospace to electronics packaging and predictive simulations of their response is playing an increasing role in understanding the molecular origin of their properties and complementing experiments in the search for tailored materials for specific applications. It focuses on modeling and simulation of the process of curing to predict the molecular structure of these polymers and their thermomechanical response by all-atom molecular dynamics simulations. Results from Monte Carlo and coarse-grained simulations are briefly summarized. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 103–122     

Delineating Protein–Protein Curvilinear Dissociation Pathways and Energetics with Naïve Multiple-Walker Umbrella Sampling Simulations

The protein–protein interaction energetics can be obtained by calculating the potential of mean force (PMF) from umbrella sampling (US) simulations, in which samplings are often enhanced along a predefined vector as the reaction coordinate. However, any slight change in the vector may significantly vary the calculated PMF, and therefore the energetics using a random choice of vector may mislead. A non-predefined curve path-based sampling enhancement approach is a natural alternative, but was relatively less explored for protein–protein systems. In this work, dissociation of the barnase–barstar complex is simulated by implementing non-predefined curvilinear pathways in US simulations. A simple variational principle is applied to determine the lower bound PMF, which could be used to derive the standard free energy of binding. Two major dissociation pathways, which include interactions with the RNA-binding loop and the Val 36 to Gly 40 loop, are observed. Further, the proposed approach was used to discriminate the decoys from protein–protein docking studies. © 2019 Wiley Periodicals, Inc.     

A Theoretical study on the charge and discharge states of Na-ion battery cathode material,Na1+xFePO4F

Na₂FePO₄F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na₂FePO₄F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na_1.5FePO₄F has the P2₁/b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P2₁/b11 Na_1.5FePO₄F and Pbcn Na₂FePO₄F or Na₁FePO₄F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na₂FePO₄F has diffusion paths in the a- and c-axis directions, Na₂FePO₄F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.     

Monte Carlo simulation of electric conductivity for pure and doping fullerene (C60)

In this paper, we have used Monte Carlo (MC) method to simulate and study the temperature and doping effects on the electric conductivity of fullerene (C60). The results show that the band gap has reduced by the doping and the charge carrier transport is facilitated from valence band to conduction band by the temperature where is touched a 300 K. In this case, the conductivity reached a value of $4\times {10}^{-7}\text{S}{\text{cm}}^{-1}$. The electric conductivity of C60 can increase by the triphenylmethane dye crystal violet (CV) alkali metal to reach $4\times {10}^{-3}\text{S}{\text{cm}}^{-1}$ at 303 K. Our results of MC simulation have a good agreement with those extracted from literature [10], [33].