He-HI复合物势能面及微分散射截面的理论研究

采用超分子单双迭代(包括非迭代三重激发)耦合簇理论CCSD(T)方法和由键函数3s3p2d1f组成的大基组, 计算得到了基态He-HI复合物相互作用的全程势能面. 该势能面上存在2个势阱, 分别对应于线性He-I-H和He-H-I构型, 势阱深度分别为4.473和2.996meV, He原子到HI分子质心的距离R分别为0.363和0.442nm. 使用Barker, Fisher和Watts提出的BFW势函数拟合计算得到的相互作用能数据, 获得了He原子与HI分子相互作用势的解析表达式. 在 关键词： He-HI复合物 势能面 微分散射截面     

有效质量法调控原子玻色-爱因斯坦凝聚体的双阱动力学

操控原子玻色-爱因斯坦凝聚体在双势阱中的动力学通常是通过改变势阱深度来实现,本文提出了一种基于调节原子有效质量的控制方案,可以在不改变双阱势的前提下操控凝聚体的双阱动力学.利用双模近似,本文解析地导出了超冷原子在双阱势中的隧穿强度和相互作用强度对有效质量的依赖关系,并基于平均场近似数值模拟了在有效质量调节下的凝聚体动力学演化,展示了隧穿振荡和自束缚等典型的双阱动力学行为.此外,本文的研究还发现,借助负有效质量效应,这一方案甚至可以等效地实现对负散射长度原子凝聚体双阱动力学行为的操控.     

He-HF体系势能模型对散射截面影响的理论研究

用BFW势函数拟合在CCSD(T)/aug-cc-pVQZ理论水平下计算的He-HF相互作用能，获得了He原子与HF分子相互作用的各向异性势，并与其他势模型进行比较，验证了拟合势的可靠性；然后采用量子密耦(Close-Coupling)方法分别计算了He-HF碰撞体系在五种不同势能模型下的微分散射截面、分波散射截面和总截面，并对计算结果进行了详细的比较和分析.研究表明：势能球平均零点能位置、势阱深度、排斥势的强度以及势能在势阱附近的方向性都对散射截面有较大影响. 关键词： 势能模型 密耦方法 散射截面     

碱金属Li原子与CO分子体系相互作用势从头计算研究

采用三重激发项校正CCSD的CCSD(T)方法和AUG-CC-PVnZ(n=2,3,4)基组,优化了CO基态分子结构,并计算了碱金属Li原子与CO分子的相互作用势,共1010个构型势能点得到体系的势能面.结果表明:同一方法下,不同基组得到的CO基态分子的键长、能量等均与实验符合很好.Li-CO势能面体现较小的各向异性势,存在两个势阱,且都为非严格T型结构;基组重叠误差(BSSE)对相互作用势的影响比较明显,采用CP方法(Counterposie method)消除基组重叠误差,不同基组计算的相互作用势显示较好一致.体系各项异性势势阱值远高于CO基态分子转动常数,使得碰撞产生强烈的非弹性碰撞,这将阻止协同冷却制备超冷分子.     

锂原子与部分惰性气体原子间相互作用势的计算

基于碱金属原子与惰性气体原子间相互作用的交换势主要来自碱金属原子的最外层S电子与惰性气体原子的最外层S电子之间 ,将面积分方法 (thesurfaceintegralmethod)推广 ,计算了Li原子与Ne、Ar、Kr和Xe原子间相互作用的vanderWaals势能曲线。其势阱的位置和深度与实验结果和其它理论计算结果符合得相当好     

超冷铯(60D5/2)2 Rydberg分子的双色光缔合光谱

本文主要从理论和实验上研究超冷铯（60D_5/2）₂ Rydberg分子的双色光缔合光谱.数值计算了铯60D_5/2 Rydberg原子对态的长程电多极相互作用和（60D_5/2）₂ Rydberg分子的绝热势能曲线,获得了（60D_5/2）₂ Rydberg分子的势阱深度和平衡间距.实验上利用双色光缔合超冷铯原子的方法制备了（60D_5/2）₂ Rydberg分子.其中,第一色激光（pulse-A）双光子共振激发种子Rydberg原子A;第二色激光（pulse-B,失谐于分子的束缚能）共振激发第二个Rydberg原子B,原子A与B由分子势阱束缚形成超冷（60D_5/2）₂ Rydberg分子.由脉冲场电离探测技术获得Rydberg分子的光缔合光谱,测量的Rydberg分子的势阱深度与理论计算结果相一致.     

利用原子芯片上Z形磁阱囚禁中性87Rb原子

利用解析和数值方法计算了Z形磁阱的囚禁势,发现当囚禁中心和芯片表面距离较远时(该距离和Z形线中部导线的一半长度相差不超过一个量级),势阱的深度不能近似表示成偏置磁场By对应的能量,而要减去囚禁中心的势能高度;而增加By进行磁阱压缩到一定值时,势阱深度反而会下降.此外介绍了原子芯片的制作方法,以及利用原子芯片上Z形磁阱囚禁中性87Rb原子的实验装置和实验过程.最终有2×106个87Rb原子被转移到Z形磁阱中.     

长程铯里德堡分子的势能曲线

本文介绍了半经典近似下的低能电子-原子散射理论, 引入贋势描述里德堡电子与基态原子的相互作用, 数值计算了铯原子nS (n=30-60)里德堡态与6S基态原子形成的长程里德堡分子的势能曲线. 并对最外层势阱进行分析, 获得长程里德堡分子的势阱深度、平衡距离与主量子数n的关系. 为实验制备里德堡分子并进一步分析其性质提供理论依据. 里德堡分子对外界非常敏感, 可用于微弱信号的检测.     

原子芯片上指数型布居增长的原子输运

超冷原子体系中要观测增益平衡的宇称-时间(PT)对称,需产生受控的原子布居增益/损耗及相干耦合。本文提出了动态控制原子芯片上双阱中原子输运实现原子布居指数增长的方法。采用直接蒙特卡罗方法数值研究了原子系综在双阱间的输运动力学,发现初始原子数和温度会显著影响输运效果,如目标势阱中的原子布居增益速率和转移效率。此外,还细致分析了左侧势阱抬升时间对右侧势阱中原子布居增长趋势的影响效果。该方案为在超冷原子气体中实现有增益/损耗的PT对称量子体系提供了切实可行的方法。     

一个新的方阱流体状态方程

将方阱势能函数视为硬球势的一短程均匀负势的叠加,并由统计力学方法导出了相应的状态方程,这个方程包含两个定标的径向分布函数。利用已知的压缩因子分子动力学模拟数据,回归得到了这两个径向分布函数,从而建立了一个新的方肼流体状态方程。它的吸引项要比Alder等人回归得到的简单得多,与方阱流体压缩因子的MD数据一致性也更好,特别是在高密度区域。     

Self-evolving atomistic kinetic Monte?Carlo: fundamentals and applications

The fundamentals of the framework and the details of each component of the self-evolving atomistic kinetic Monte Carlo (SEAKMC) are presented. The strength of this new technique is the ability to simulate dynamic processes with atomistic fidelity that is comparable to molecular dynamics (MD) but on a much longer time scale. The observation that the dimer method preferentially finds the saddle point (SP) with the lowest energy is investigated and found to be true only for defects with high symmetry. In order to estimate the fidelity of dynamics and accuracy of the simulation time, a general criterion is proposed and applied to two representative problems. Applications of SEAKMC for investigating the diffusion of interstitials and vacancies in bcc iron are presented and compared directly with MD simulations, demonstrating that SEAKMC provides results that formerly could be obtained only through MD. The correlation factor for interstitial diffusion in the dumbbell configuration, which is extremely difficult to obtain using MD, is predicted using SEAKMC. The limitations of SEAKMC are also discussed. The paper presents a comprehensive picture of the SEAKMC method in both its unique predictive capabilities and technically important details.     

Molecular dynamics simulation of atomic hydrogen diffusion in strained amorphous silica

Understanding hydrogen diffusion in amorphous SiO_2(a-SiO_2),especially under strain,is of prominent importance for improving the reliability of semiconducting devices,such as metal-oxide-semiconductor field effect transistors.In this work,the diffusion of hydrogen atom in a-SiO_2 under strain is simulated by using molecular dynamics(MD) with the ReaxFF force field.A defect-free a-SiO_2 atomic model,of which the local structure parameters accord well with the experimental results,is established.Strain is applied by using the uniaxial tensile method,and the values of maximum strain,ultimate strength,and Young's modulus of the a-SiO_2 model under different tensile rates are calculated.The diffusion of hydrogen atom is simulated by MD with the ReaxFF,and its pathway is identified to be a series of hops among local energy minima.Moreover,the calculated diffusivity and activation energy show their dependence on strain.The diffusivity is substantially enhanced by the tensile strain at a low temperature(below 500 K),but reduced at a high temperature(above500 K).The activation energy decreases as strain increases.Our research shows that the tensile strain can have an influence on hydrogen transportation in a-SiO_2,which may be utilized to improve the reliability of semiconducting devices.     

Mg-Al合金熔体中固液界面结构的分子动力学研究

本文采用分子动力学研究了Mg-3%Al合金熔体中固液界面结构及界面附近原子的扩散行为.计算结果表明,该二元合金的固液界面为粗糙界面.垂直于界面方向的数密度分布,表现出复杂波动的特征,这种波动一直延伸到液体中.在界面附近的区域,扩散系数的三个不同方向的分量表现出了明显的各向异性,并且这种各向异性一直持续到液相当中.对界面二维结构的分析表明,界面附近液相原子的二维排列呈现出从长程有序逐渐过渡到短程有序的变化.     

氢及其同位素在聚苯乙烯与Cu界面的扩散

为了探索气体在固体表面高分子链中的扩散,使用分子动力学（MD）的方法,对H2,D2,T2在聚苯乙烯与金属铜（PS-Cu）界面的扩散进行了计算模拟,通过所得到气体的均方位移计算了气体在不同金属表面与聚苯乙烯界面中的扩散系数。结果显示:气体在界面的扩散系数比在聚苯乙烯本体中的扩散系数小,气体在PS-Cu(110)界面的扩散系数最大,在PS-Cu(111)界面的扩散系数最小。计算和分析了PS与金属表面的相互作用,发现其相互作用能越大,气体在此界面的扩散系数越小。同时,金属表面的晶面密度对气体在界面中的扩散也有一定的影响。     

Universal Long-Time Relaxation on Lattices of Classical Spins: Markovian Behavior on Non-Markovian Timescales

The long-time behavior of certain fast-decaying infinite temperature correlation functions on one-, two-, and three-dimensional lattices of classical spins with various kinds of nearest-neighbor interactions is studied numerically, and evidence is presented that the functional form of this behavior is either simple exponential or exponential multiplied by cosine. Due to the fast characteristic timescale of the long-time decay, such a universality cannot be explained on the basis of conventional Markovian assumptions. It is suggested that this behavior is related to the chaotic properties of the spin dynamics.     

Multiscale molecular dynamics using the matched interface and boundary method

The Poisson-Boltzmann (PB) equation is an established multiscale model for electrostatic analysis of biomolecules and other dielectric systems. PB based molecular dynamics (MD) approach has a potential to tackle large biological systems. Obstacles that hinder the current development of PB based MD methods are concerns in accuracy, stability, efficiency and reliability. The presence of complex solvent-solute interface, geometric singularities and charge singularities leads to challenges in the numerical solution of the PB equation and electrostatic force evaluation in PB based MD methods. Recently, the matched interface and boundary (MIB) method has been utilized to develop the first second order accurate PB solver that is numerically stable in dealing with discontinuous dielectric coefficients, complex geometric singularities and singular source charges. The present work develops the PB based MD approach using the MIB method. New formulation of electrostatic forces is derived to allow the use of sharp molecular surfaces. Accurate reaction field forces are obtained by directly differentiating the electrostatic potential. Dielectric boundary forces are evaluated at the solvent-solute interface using an accurate Cartesian-grid surface integration method. The electrostatic forces located at reentrant surfaces are appropriately assigned to related atoms. Extensive numerical tests are carried out to validate the accuracy and stability of the present electrostatic force calculation. The new PB based MD method is implemented in conjunction with the AMBER package. MIB based MD simulations of biomolecules are demonstrated via a few example systems.     

Alchemical molecular dynamics for inverse design

We present a molecular dynamics (MD) implementation of an extended statistical mechanical ensemble that includes ‘alchemical’ degrees of freedom describing particle attributes as thermodynamic variables. We demonstrate the use of this alchemical MD method in inverse design simulations of particles interacting via the Oscillating Pair Potential (OPP) and the Lennard–Jones–Gauss potential (LJG) – two general, previously studied models for which phase diagrams are known. We show that alchemical MD can quickly and efficiently optimise pair potentials for target structures within a specified design space in the low-temperature regime, where internal energy adequately represents the features of the alchemical free energy landscape. We show that alchemical MD can be also used to inversely design pair potentials to achieve target materials properties (here, bulk modulus) directly, without explicit knowledge of the structure–property relationship. Alchemical MD can easily be generalised and applied to any target materials properties or structures and used with any differentiable interaction potential.     

抛物量子阱中束缚极化子的极化势和结合能

利用改进的Lee-Low-Pines(LLP)方法,用变分法计算了无限深抛物量子阱中同时考虑与体纵光学声子和界面纵光学声子相互作用的束缚极化子的极化势和结合能.数值计算得出:阱宽较大时极化势很小,阱宽较小时极化势较大,所以对于较窄的抛物阱必须考虑极化势.对于给定阱宽的抛物阱,随着远离阱中心极化势迅速减小,当到达阱的界面附近极化势又开始增大.阱宽较小时,束缚极化子的结合能随着阱宽L的增大而急剧减小;阱宽较大时,结合能减小的非常缓慢,最后接近体材料中的三维值.     

The exciton luminescence in Zn(Cd)Se/ZnMgSSe quantum wells

We study exciton states in Zn(Cd)Se/ZnMgSSe quantum wells (QWs) with various degrees of diffusion blurring in the interfaces by the methods of optical spectroscopy. We show that at low temperatures the QW emission spectra are determined by free and neutral donor-bound excitons. Blurring of the heterointerfaces leads to the increase in the energy shift between the emission line maxima of free and bound excitons. We explain the nonlinear dependence of the steady-state photoluminescence intensity on the excitation-power density in terms of the neutralization of charged donors at the photoexcitation of heterostructures. We observed a complex long-time dynamics of the reflection coefficient, evoked by the charge-redistribution processes in the heterostructure, near the QW exciton resonances under the irradiation.     

Diffusion of 1,2,3-Benzotriazole as a Volatile Corrosion Inhibitor through Common Polymer Films Using the Molecular Dynamics Simulation Method

Diffusion of 1,2,3-benzotriazole (BT), as one of the volatile corrosion inhibitors (VCI) for copper and steel, through several polymers was investigated using molecular dynamics simulation (MD). MD were performed by employing the COMPASS force field to estimate the diffusivity of BT through polyethylene (PE), polypropylene (PP), polyvinyl chloride (PVC), and Nylon 6 as potential hosts for anticorrosion film packaging purposes. The diffusion coefficients (D) of BT in these polymers were calculated by constructing an amorphous cell, each containing BT and one of these polymers. After constructing the cell, a molecular dynamics simulation was performed to calculate the mean square displacement of the BT molecule. Simulation results showed that BT can diffuse through PVC easier than the other polymers. Among these hosts the Nylon 6 had the lowest D value, implying that this polymer can maintain BT for a long time. The temperature dependence of diffusion through PE, as the most common VCI film, was studied and the activation energy (Q_d) and pre-exponential diffusion coefficient (D₀) in Arrhenius equation were calculated.