BNCT优化网格设计及相关算法研究

用MCNP蒙特卡罗程序模拟了硼中子俘获治疗(BNCT)3种国际基准网格模型, 并与 修正的Snyder椭球模型进行了比较. 在此基础上, 给出了一种保质量守恒、内存量少、易于产生输入文件的4种基本材料成分的BNCT网格模型. 计算结果表明, 在4mm网格下, 新模型可以达到基准模型的精度; 根据解析模型剂量随深度的变化规律, 研究构造了多网格组合模型, 在重要区域计算精度不损失的条件下, 计算时间大大缩短. 最后研究给出了一个既保证精度、又在可接受的时间内完成剂量计算的模型、样本数和相应的算法, 它基本上满足临床BNCT的要求.     

Combinational simulations of free energy of polymer solution

It is conceptually proposed that the total entropy of polymer solution is contributed from two distinct parts: the positional and the oomformational. The former can be represented analytically, while the latter can be simulated with the random self-avoiding walk model on the simple cubic lattice for multichain systems. The obtained results indicated that both the conformational entropy and the mixing heat are consistent with the scaling laws wry well.     

High-precision Monte Carlo test of the conformai-invariance predictions for two-dimensional mutually avoiding walks

Let _l be the critical exponent associated with the probability thatl independentN-step ordinary random walks, starting at nearby points, are mutually avoiding. Using Monte Carlo methods combined with a maximum-likelihood data analysis, we find that in two dimensions ₂=0.6240±0.0005±0.0011 and ₃=1.4575±0.0030±0.0052, where the first error bar represents systematic error due to corrections to scaling (subjective 95% confidence limits) and the second error bar represents statistical error (classical 95% confidence limits). These results are in good agreement with the conformal-invariance predictions ₂=5/8 and ₃=35/24.     

Branched polymers in a wedge geometry in three dimensions

We investigate the statistical and dimensional properties of uniform star polymers attached by the branching vertex of degreef in a wedge geometry in three dimensions and described by the wedge angles and. We show that the growth constant is equal to ^f , where is the self-avoiding walk limit. Thef and (, ) dependences of the corresponding critical exponent _f (, ) are studied using Monte Carlo techniques. In the casef=1, our results are compared with existing predictions obtained from series expansion and renormalization group methods. We have also estimated the amplitudes for the mean square radius of gyration and the mean square end-to-end branch length. Our results for the ratio of the mean square radius of gyration of anf-star to that of a linear polymer of the same degree of polymerization attached in a similar wedge, and the analogous ratio for the mean square end-to-end branch length, are consistent with these ratios being lattice-independent quantities.     

Superdegenerate point in FCC phase diagram: CVM and Monte Carlo investigations

We investigate the topology of the phase diagram of binary alloys on the fee lattice with first-neighbor antiferromagnetic interactions around the superdegenerate point, where the L1₀ and L1₂ phases meet. We treat the system as a hard-constraint lattice gas, following a procedure previously described by Lebowitzet al. We perform cluster variation method calculations in theT0 limit and Monte Carlo simulations directly atT=0 K on the ground states of the superdegenerate point. We find that: (i) there is no disordered phase in the neighborhood of this point; (ii) a phase L for which two of the four cubic sublattices have the same average occupancy and each of the two others are different appears between L1₀ and L1₂; (iii) the transition L/L1₂ is of first order.     

Likelihood estimation of soft-core interaction potentials for Gibbsian point patterns

The likelihood method is developed for the analysis of socalled regular point patterns. Approximating the normalizing factor of Gibbs canonical distribution, we simultaneously estimate two parameters, one for the scale and the other which measures the softness (or hardness), of repulsive interactions between points. The approximations are useful up to a considerably high density. Some real data are analyzed to illustrate the utility of the parameters for characterizing the regular point pattern.     

On the spectrum of the Heisenberg Hamiltonian

The quantum, antiferromagnetic, spin-1/2 Heisenberg Hamiltonian on thed-dimensional cubic lattice ^d is considered for any dimensiond. First the anisotropic case is considered for small transversal coupling and a convergent expansion is given for a family of eigenprojections which is complete in all finite-volume truncations. Then the general case is considered, for which an upper bound to the ground-state energy is given which is optimal for strong enough anisotropy. This bound is expressed through a functional involving the statistical expectation value at finite temperature of a certain correlation function of an Ising model defined on the lattice ^d itself.     

MONTE CARLO STUDY ON THE CRITICAL ADSORPTION POINT OF BOND-FLUCTUATED POLYMER CHAINS TETHERED ON ADSORBING SURFACES

 The behavior of three-dimensional bond fluctuation model chains tethered on an adsorbing flat surface was simulated by the Monte Carlo method. The dependence of the number of surface contacts M on the interaction strength e and the chain length N was investigated by a finite-size scaling law M = N[a₀ + a₁N^1/v_k + O((N ^1/v_k)²)] for e near the critical adsorption point e_c, i.e., k ≡(e-e_c)/e_c closes to 0. The critical adsorption point was estimated to be e_c = 0.93, and the exponents  = 0.49 and 1/v= 0.57.     

模板剂与ＭＣＭ-22分子筛匹配作用的分子模拟计算

本文用Ｍｏｎｔｅ Ｃａｒｌｏ方法和分子力学方法研究了模板剂分子六次甲基亚胺(ＨＭＩ)在ＭＣＭ-22分子筛孔道中的排布与取向。模拟结果表明模板剂分子能够稳定地嵌入在ＭＣＭ-22分子筛的超笼、窗口和正弦孔道中。对这三种位置上模板剂与分子筛骨架之间的主客体非键相互作用能(包括范德华能和库仑能)进行了计算,结果表明模板剂分子与分子筛的骨架之间作用力主要以范德华吸引力为主。本文还进一步讨论了模板剂分子的结构导向作用和对骨架Ａｌ的靶向作用。     

大分子固液界面吸附行为的数字模拟研究方法

介绍了随机数字模拟方法在研究链状大分子吸附行为的一般原理,以立方格子模型的自避行走为例,分析了模型的建立、初始位形的产生以及分子链运动的处理方法.回顾了随机数字模拟方法在研究大分子界面吸附行为方面的发展状况.预测了随机数字模拟方法在大分子吸附行为研究方面的应用前景.