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

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

The crystal structure on the set of Mirkovi?-Vilonen polytopes

In an earlier work, we proved that MV polytopes parameterize both Lusztig's canonical basis and the Mirkovi?-Vilonen cycles on the affine Grassmannian. Each of these sets has a crystal structure (due to Kashiwara-Lusztig on the canonical basis side and due to Braverman-Finkelberg-Gaitsgory on the MV cycles side). We show that these two crystal structures agree. As an application, we consider a conjecture of Anderson-Mirkovi? which describes the BFG crystal structure on the level of MV polytopes. We prove their conjecture for sln and give a counterexample for sp₆. Finally we explain how Kashiwara data can be recovered from MV polytopes.     

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.     

A complementary thermodynamic limit for classical Coulomb matter

The canonical equilibrium measure of classical two-component Coulomb matter with regularized interactions is analyzed in a finite volume. It is shown that, in the mean-field regime, the one-particle density is inhomogeneous on a new characteristic length scale _inh. For a system ofN positive andN negative particles, _inh and the characteristic length scale of correlations _corr (=Debye screening length) are related via _inh=(2N)^1/2 _corr. The major conceptual conclusion that is drawn from this is that one needs two nontrivial complementary thermodynamic limits to define the equilibrium thermodynamics of two-component Coulomb systems. One of them is the standard thermodynamic limit (infinite volume), where one takesN, _corr fixed. Its complementary limit is characterized byN, _inh fixed, and is a finite-volume inhomogeneous mean-field limit. The most prominent new feature in the mean-field thermodynamic limit, which is absent in the standard thermodynamic limit, is an anomalous first-order phase transition where the Coulomb system explodes or implodes, respectively. The phase transition is connected with the existence of a metastable plasma phase far below the ionization temperature.     

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.     

Surface properties of finite classical Coulomb systems: Debye-Hückel approximation and computer simulations

We report analytical and numerical studies of surface correlations in finite, homogeneously polarizable, classical Coulomb systems placed in an insulating or conducting environment. Their purpose is to understand the phenomenological, shape-dependent laws of electrostatics, from the point of view of statistical mechanics; we focus on the knowledge of the dielectric susceptibility of the system, a quantity proportional to the equilibrium fluctuation of the system's instantaneous polarization per unit volume. This goal has been achieved for a system in a conducting state. The picture is that the shape-dependent part of the susceptibilities results from the action of unbounded observables (the second moments of the instantaneous polarization of the system) on long-range surface correlations and that the relations of electrostatics are verified by means of shape-dependent thermodynamic limits. This picture is supported (i) by exact solutions and asymptotic analysis of the Debye-Hückel approximation of multicomponent plasmas in disks and spheres with insulating and conducting environment and also in ellipses in a vacuum, and (ii) by computer simulations of a one-component plasma in a disk with different environments, notably a conducting environment with permeable and impermeable wall. These observations have revealed for the first time the reason why the susceptibility of a conducting disk in a conductor with impermeable walls diverges linearly with the radius of the disk: this is due to the occurrence of long-range radial correlations in the conductor. These findings are quantitatively interpreted in terms of a novel canonical Debye-Huckel approximation as contrasted to the ordinary grand canonical version. Lastly a fresh look at the problem of the surface correlations of a conductor in a vacuum, which places the observer close to the surface of the conductor but in the vacuum, is presented and applied to the disk, the ellipse, the cylinder, the sphere, and the wedge.     

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.