蒙特卡罗模拟中相关变量随机数序列的产生方法

蒙特卡罗模拟有时需要对多维相关随机变量进行模拟抽样．本文介绍基于Choesky因子线性变换一非线性变换产生具有指定边缘分布和相关系数的多维相关随机变量抽样序列的一般方法,给出一种简单易行的高效数值实现途径和一些模拟结果．模拟结果表明,该方法产生的各随机变量抽样序列间具有预期要求的相关性,并能通过指定边缘分布Kolmogorov—Smimov非参数假设检验．对该方法应用中的一些限制问题进行了讨论．     

局部Klein-Nishina分布抽样

局部Klein-Nishina分布是指，光子与物质发生非相干散射后对散射方向的抽样限定在局部范围内。近一些年来，一种通过确定对探测器所张立体角计算探测器通量的方法被提了出来，并受到了普遍重视，局部Klein-Nishina分布抽样问题就是由此产生的。较全面地研究了局部Klein-Nishina分布抽样问题，给出了一些可行的抽样方法，并在计算机上进行了模拟抽样比较。     

非理想表面对声子态密度的影响——一维简谐链模型和格林函数方法

本文讨论非理想表面对一维简谐链单声子态密度函数的影响,发现非理想表面可以导致声子的软化或硬化。     

振幅衰减模型中弱相干光场的退相干

利用热纠缠态表象求解密度矩阵主方程的方法,给出了任意时刻振幅衰减模型中弱相干场密度算符的表示式。采用数值计算方法计算了弱相干光场的压缩效应、统计性质和Wigner函数。讨论了耗散对这些量子特性的影响。研究结果表明:随衰减时间的增大,弱相干场的压缩效应、统计性质和Wigner函数的负性等非经典效应均减弱。     

使用非正交曲线座标与速度分量S_1流面正问题流场矩阵解

本文基于吴仲华提出的使用非正交曲线座标与相应的非正交速度分量的叶轮机械三元流动气动热力学基本方程组,引入流函数,得出求解的主要方程:流函数的二阶拟线性偏微分方程.除了与密度有关的项以外,流函数的各阶导数都置在方程的左端.这样加快了收敛的速度.用中心九点差分格式,将微分方程离散化后,所得的线性代数方程组用矩阵[L][u]分解直接求解.这种解法收敛速度较快.系数矩阵为对角线带状稀疏矩阵.采用了:(1)非零元素按对角线编号;(2)增设虚点两项办法.大量减少了计算机内存量.由流函数求密度时采用了内存密度函数表插值方法.简单地讨论了松弛因子的选取.用此程序对一些压气机和透平的叶型进行了计算,同实验结果及理论解析解进行了比较,相互是一致的.     

局部Klein-Nishina分布抽样

局部Ｋｌｅｉｎ－Ｎｉｓｈｉｎａ分布是指，光子与物质发生非相干散射后对散射方向的抽样限定在局部范围内。近一年来，一种通过确定对探测器所张立体角计算探测器通量的方法被提了出来，并受到了普遍重视，局部Ｋｌｅｉｎ－Ｎｉｓｈｉｎａ分布抽样问题就是由此产生的。较全面地研究了局部Ｋｌｅｉｎ－Ｎｉｓｈｉｎａ分布抽样问题，给出了一些可行的抽样方法，并在计算机上进行了模拟抽样比较。     

超强短脉冲激光在不同密度分布等离子体中的自聚焦

讨论高斯型强激光束在具有初始柱对称密度分布的低密度冷等离子体中传播时，等离子体密度分布的不同对激光自聚焦的影响.推导出可以判断更有利于自聚焦发生的评价函数，这样通过比较不同密度分布的评价函数值就可以判断哪种密度分布更有利于自聚焦的发生.为了说明这种方法的有效性，对评价函数进行分析得出：在相同的激光场中等离子体柱的轴心密度给定时(以激光的光轴为轴)，离轴越远的地方密度越大及密度变化越陡，自聚焦越容易发生；相对论效应与有质动力共同作用比相对论的单独作用，自聚焦更容易发生.数值模拟证实了评价函数能准确的预测在不 关键词： 自聚焦 相对论效应 有质动力 评价函数     

蒙特卡罗技术在放射诊断剂量快速计算中的应用

采用蒙特卡罗（MC）技术研究放射性核素131I的放疗特性。给出了131I衰变过程的MC抽样表达式，实现了含电子和光子的混合抽样模拟。对水客体中131I衰变的MC模拟结果表明:131I能够有效地治疗甲状腺肿瘤而几乎不影响正常组织。变密度客体的MC模拟结果表明:沉积能量与面密度之间的关系与密度无关，满足卷积核自适应的要求；并且MC模拟能够提供剂量卷积核函数，为快速剂量计算提供数据。     

一维介观系统的隧道电流零频散粒噪声谱密度

导出了一维介观系统隧道结的以非平衡格林函数(闭路格林函数)表示的电流散粒噪声谱密度表达式,对隧道结中无杂质及含杂质两种情形下的电流零频散粒噪声谱密度随系统各物理参数的关系做了详细的计算及讨论. 关键词：     

表面对V_3Si电子态密度的影响

本文从V_3Si的准一维链模型出发,用双时格林函数方法,正确地计算了链的理想中断所引起的自由表面、非理想中断的表面以及含杂的自由表面对电子态密度的影响,分析了可能提高V_3Si费米面附近态密度的物理条件,讨论了超细微颗粒化提高V_3Si超导转变温度的可能性.     

卷积与误差分布密度函数的合成

实验误差并不是单一的，它有可能是几个独立产生的误差之和，本文讨论了其和的误差分布密度函数，并给出了几种典型的误差分布密度函数的合成结果．     

Si4团簇电子输运性质的第一性原理计算

采用密度泛函理论和非平衡格林函数相结合的方法对Si₄ 团簇与Au(100)电极空位相连的纳米结点的电子输运性质进行了理论模拟计算, 得到了纳米结点在不同距离下的几何结构、电子结构、电导、透射谱、电荷转移量; 讨论了当距离d_z=12.004 Å时纳米结点的电导、电流随电压的变化关系. 关键词： 密度泛函理论 非平衡格林函数 4团簇')" href="#">Si₄团簇 电子输运     

导体表面的电荷分布

从对称性角度出发讨论了导体表面的电荷分布问题，并对“电荷密度分布函数”的理论作了评述。     

聚合物侧基NMR弛豫相关函数

本文导出了聚合物为链段运动时,侧基进行非等价二位置和三位置跃迁内旋转、随机自由扩散内旋转以及受阻内旋转时的NMR弛豫相关函数,同时给出了相应的谱密度函数,并讨论了两个应用实例。     

对称性关系与研究强子物理的自洽非微扰途径

讨论了对称性导致的格林函数间的关系——横向及通常的纵向Ward-Takahashi关系. 自洽求解这组WT关系，导出了非微扰的费米子－玻色子项角函数. 应用非微扰顶角函数，得到了自洽截断的对传播子封闭的Dyson-Schwinger方程，构成了研究强子物理的自洽非微扰途径.     

非均匀大气中γ射线输运的蒙特卡罗模拟方法

γ射线在高空大气中的输运问题常采用蒙特卡罗方法进行模拟,但由于大气是非均匀连续介质,通常需要根据高度上的密度变化将大气进行分层处理.高空大气密度变化非常剧烈,精细的分层会带来很大的计算量.针对非均匀连续介质中的粒子输运问题,利用质量距离抽样方法代替传统的步长抽样方法,提出一种可直接模拟粒子在非均匀连续介质中输运的蒙特卡罗方法,实现γ射线在非均匀大气中输运的蒙特卡罗模拟,并与MCNP程序的计算结果进行比较,结果符合较好.     

C74分子光谱特性及电子结构的密度泛函研究

利用第一性原理的密度泛函理论和非平衡格林函数的方法研究了C74分子的分子轨道密度、振动光谱和态密度。结果显示C74分子的能隙0.6eV。比较C74分子的电子云分布发现，在HOMO和LUMO状态下的电子只分布在C74分子的半球表面侧，电子在成键过程中都表现出了较好的离域特性。C74分子的红外光谱和拉曼光谱显示，碳碳键的振动在有些频率处只有红外活性，有些频率处只有拉曼活性，还有一些频率处既有红外又有拉曼活性。另外，考虑电子自旋作用时C74分子态密度的极大值数量有所增加，并使简并的能级发生分裂。     

跨声速流函数方程强隐式解及确定密度场的新方案

基于人工可压缩性对密度的修正,本文用强隐式格式快速求解了非正交曲线坐标系下跨声速流函数方程;在流函数场解出后,通过求解一个由动量方程、能量方程和连续方程组合而成的关于密度的一阶偏微分方程来获得密度场,因此流函数解法中常遇到的密度双值问题在这里已不存在;通常所讲的完成流函数场{Ψ}与密度场{ρ}间的迭代在本文便体现在流函数主方程与这个新推出的一阶微分方程间的迭代计算上;几个典型算例表明了本方法的有效性。     

CFHH—GLF方法直接求解氢分子体系中相关函数的选取及相应矩阵元 …

给出了在非正则超球坐标系下用相关函数超球谐－广大拉盖尔函数方法（ＣＦＨＨ－ＧＬＦ）对氢分子薛定谔方程进行直接求解时相关函数的选取原则，引入相关函数，进行了方程的推导，并列出所涉及的相应矩阵元的解析表达式，在推导过程中，利用了合流超几何函数、拉卡系数和维格纳系数的一些性质。     

光子增加混沌场的退相干和非经典效应

将产生算符作用在混沌场上,构造了光子增加混沌场.利用有序算符内的积分技术和热纠缠态表象求解密度矩阵主方程的方法,研究了振幅衰减模型中光子增加混沌场的退相干和非经典效应.通过解振幅衰减模型中的密度算符的主方程,得到了初态为光子增加混沌场的密度算符的演化公式.计算了终态密度算符的P表示和Wigner函数,并数值计算了耗散对其P表示和Wigner函数的影响.结果表明:随耗散时间的增长,光子增加混沌场的非经典效应减弱.另一方面,随光子增加数的增加,其非经典效应也减弱.