BaTiO3晶体结构及弹性的分子动力学模拟

有效的势函数是分子动力学模拟的关键. 引入了一种势函数,该势函数的特点是运用参数r_eff计算原子间的静电作用. 通过分子动力学方法模拟得到了BaTiO₃晶体立方相、四方相结构的对关联函数和X射线衍射谱,计算得出了它们的晶格常数及弹性常数. 模拟结果与实验结果符合较好.该势函数可以有效地模拟BaTiO₃晶体的热学和力学性能. 关键词： 分子动力学模拟 势函数 3铁电晶体')" href="#">BaTiO₃铁电晶体     

利用电子结构数据拟合H2分子的电子壳模型势函数参数

电子壳模型势函数在离子晶体的原子级计算机模拟中有广泛应用,其势参数主要通过拟合晶体的实验数据或电子结构数据得到.提出了通过拟合双原子分子的量子化学从头计算电子结构数据来获得该势函数的方法,并由H2分子的电子结构数据建立了H原子间的电子壳模型势函数.此外,还应用该势函数对H+2分子离子进行了计算.该势函数拟合方案更适合于共价键型的分子.     

低温下液氦、液氖热物性的分子动力学模拟

本文基于平衡分子动力学模拟方法,采用了含量子修正项的QFH势函数和WK势函数对低温下液氦和液氖进行计算,并通过Green-Kubo公式统计得到热导率和黏性系数。结果表明:液氦和液氖的热导率、黏性系数随温度升高呈线性增大;采用QFH势函数和WK势函数比采用LJ势函数计算得到的热导率更接近实验值,且WK势函数计算结果比采用QFH势获得的结果更准确;WK势函数计算的液氦和液氖的黏性系数也与实验值非常接近。     

利用电子结构数据拟合H2分子的电子壳模型势函数参数

电子壳模型势函数在离子晶体的原子级计算机模拟中有广泛应用,其势参数主要通过拟合晶体的实验数据或电子结构数据得到.提出了通过拟合双原子分子的量子化学从头计算电子结构数据来获得该势函数的方法,并由H₂分子的电子结构数据建立了H原子间的电子壳模型势函数.此外,还应用该势函数对H⁺₂分子离子进行了计算.该势函数拟合方案更适合于共价键型的分子. 关键词： 电子壳模型势 参数拟合 共价键 2分子')" href="#">H₂分子     

α-B晶体的Lennard-Jones对势和对势型多体势构建

精确原子间相互作用势函数的建立是分子动力学模拟的核心. 针对α-B晶体(R3m群), 分别构建Lennard-Jones (L-J)型对势和对势型多体势, 这两类势函数的构建仅需考虑晶体的原子平均结合能和几何构型参数. 前者取最近邻原子间距的势函数为最小值; 后者采用负指数(L-J型)函数和分段插值多项式形式来分别构造势函数, 在近邻处为势能极小值, 相邻极小值位置的中点引入势垒. 针对完整的α-B 晶体和偏离完整α-B晶体结构, 首先利用第一性原理计算结果来确定势函数的参数, 然后借助分子动力学的能量最小化方法优化结构, 并对这两类势函数以及Tersoff 势函数进行比较. 结果表明: L-J对势和Tersoff势的结果与α-B晶体构型有较大偏离; 而对势型多体势, 无论初始晶体构型完整与否, 其结果与完整α-B晶体构型比较一致.     

银纳米团簇负热容现象的分子动力学模拟

采用分子动力学方法和嵌入原子法(EAM)多体势函数,模拟研究了银纳米团簇在不同温度直到熔化过程中的结构变化,得到了体系能量和热容量随温度的变化关系.结果显示:银纳米团簇在临近熔点附近出现了负热容现象.研究了弛豫后银纳米团簇的稳态结构变化及其在不同时刻结构的演变过程.结果表明:产生负热容现象的主要原因是纳米团簇在熔点附近,结构发生了巨大的变化,形成由{111}和{100}面围成的结构十分稳定和能量更低的多面体.     

氧化硅层中的锗纳米晶体团簇量子点

采用氧化和析出的方法在氧化硅中凝聚生成锗纳米晶体量子点结构. 其形成的锗晶体团簇没有突出的棱角和支晶结构，锗晶体团簇的轮廓较圆混，故可以用球形量子点模型来模拟实际的锗晶体团簇. 对比了在长时间退火氧化条件下和在短时间退火用激光照射氧化条件下所生成的锗纳米晶体结构的PL光谱和对应的锗纳米晶体团簇的尺寸分布. 短时间退火氧化条件下生成的锗纳米晶体较小(3.28—3.96nm)，长时间退火用激光照射氧化条件下所生成的锗纳米晶体较大(3.72—4.98nm)；其分布结构显示某些尺寸的锗纳米晶体团簇较稳定，适当的氧化条件可以得到尺寸分布范围较窄的锗纳米晶体团簇. 用量子点受限模型计算了锗纳米晶体团簇的能隙结构，用Monte Carlo方法模拟了PL光谱和对应的锗纳米晶体团簇的尺寸分布，分别与实验结果符合较好. 关键词： 锗晶体团簇 纳米晶体 量子点 激光照射     

红外区二维光子晶体偏振分光镜透射谱分析

采用时域有限差分（Finite-Difference Time-Domain,FDTD）法模拟了垂直入射红外波段二维光子晶体偏振分光镜的透射谱,从而验证了用平面波展开方法计算的二维光子晶体能带结构.用电子束蒸发方法在光栅基板上制备了Ge/BaF2多层膜,实现了二维光子晶体偏振分光镜.实验初步证明了这种光子晶体在垂直入射时具有偏振特性,分析了实验结果与理论计算之间偏差的产生原因.     

纳米晶体结构与性能的模拟研究

通过分子动力学方法模拟纳米晶体（1—3nm）的结构．利用模拟的结果，进行了X射线衍射点阵常数、晶粒尺寸及点阵畸变的模拟计算，还计算了结合能及弹性模量等．结果表明纳米晶体无论是晶界和晶粒都与传统的粗晶粒晶体材料没有本质的区别，只是由于晶粒尺寸变小，以及晶界的体积分数等的作用，导致一系列的性能差异 关键词：     

金属铜Cu凝固过程中晶体生长与团簇微观结构演变特性的模拟研究

采用分子动力学方法和Quantum Sutton-Chen（QS-C）多体势,对液态金属铜（Cu）凝固过程中的晶体生长规律及纳米团簇微观结构转变特性进行了模拟跟踪研究．运用Honeycutt-Andersen（HA）键型指数法和新的原子团类型指数法（CTIM-2）分析了金属Cu原子的成键类型和原子团簇结构演变特性.结果发现：当以1.0×1013K/s速率凝固时,系统最终形成晶体和非晶体混合共存结构;在以4.0×1012K/s速度冷却时,系统从673K就开始结晶,并形成以1421和1422二种键型为主的晶体结构;面心立方（FCC）和六角密集（HCP）结构在形成晶体铜时起着非常重要的作用,尤其是由1421键型构成的面心立方（12 0 0 0 12 0）基本原子团在晶体生长和纳米团簇结构形成过程中占主导地位.     

The peculiarities of crystallization of low-temperature eutectic in magnesium-lead system

Density of low-temperature eutectic in magnesium-lead system (83.06 at. % Pb) has been measured by gamma-raying of the samples with narrow beam from cesium-137 isotope over the temperature range 293–1000 K of solid and liquid states. Approximation density dependences have been obtained and data of this work and other authors have been compared. Reference tables of temperature dependences of the alloy thermal properties have been compiled for the entire range of measurements and their errors estimated. It is shown that differences in the values of volumetric changes obtained during melting and crystallization are bound with the metastable β′-phase formation.     

Characteristic features of the melting of two-dimensional mesoscopic Wigner clusters

Two-dimensional Wigner microclusters in a semiconductor dot are studied. Their melting is investigated in detail and it is shown that, for typical mesoscopic clusters possessing a shell structure, melting occurs in two stages: orientational melting (rotation of the shells relative to one another) and total melting, where the shells start to overlap with one another and exchange particles. An example of a “magic” microstructure which has a triangular structure and melts in a single stage is presented. For this, the temperature dependences of various quantities characterizing cluster structure are investigated. The change in the distribution of cluster configurations over local minima of the potential energy with increasing temperature is investigated. At temperatures below the temperature of total melting, a cluster is always located near the configuration of a global minimum and, at temperatures above the temperature of complete melting, a cluster can be located with finite probability near configurations corresponding to various local minima of the potential energy. Fiz. Tverd. Tela (St. Petersburg) 41, 1499–1504 (August 1999)     

Generalized random phase approximation and gauge theories

Mean-field treatments of Yang-Mills theory face the problem of how to treat the Gauss law constraint. In this paper we try to face this problem by studying the excited states instead of the ground state. For this purpose we extend the operator approach to the Random Phase Approximation (RPA) well-known from nuclear physics and recently also employed in pion physics to general bosonic theories with a standard kinetic term. We focus especially on conservation laws, and how they are translated from the full to the approximated theories, demonstrate that the operator approach has the same spectrum as the RPA derived from the time-dependent variational principle, and give—for Yang-Mills theory—a discussion of the moment of inertia connected to the energy contribution of the zero modes to the RPA ground state energy. We also indicate a line of thought that might be useful to improve the results of the Random Phase Approximation.     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

The effect of the topography of heterogeneous surfaces on the adsorption energy distribution,evaluated by the asymptotically correct,condensation approximation

In our previous paper [1], we used the well-known CA(Condensation Approximation) method, to evaluate the adsorption energy distribution from experimental isotherms, assuming that adsorption sites of different adsorption energies are distributed on adsorbent surface completely at randon.This note is an extension of our theory, using the ACCA, (Asymptotically Correct Condensation Approximation), which, in the hitherto investigation, was used only to adsorption on surfaces, with the patchwise topographical distribution of adsorption sites.     

碳纳米管-硅纳米线复合结构的形成和热稳定性

通过分子动力学方法模拟了在碳纳米管内填充一定数目的半导体元素硅形成碳纳米管-硅纳米线复合结构的过程,并采用Lindemann指数研究了这种复合结构的热稳定性．计算结果表明,当考虑碳纳米管和硅纳米线轴向方向的周期性边界条件之后,在C（13,0）和C （14,0）碳纳米管内能够形成亚稳结构的硅纳米线Si₁₆NW和Si₂₀NW,从而获得一种碳纳米管-硅纳米线的新型复合结构．通过计算这种复合结构的Lindemann指数,可以看到由于碳纳米管的保护作用,碳纳米管包裹的硅纳 关键词： 复合结构 纳米线 碳纳米管 分子动力学     

Thermal conductivity of liquid rubidium in the interval of 312–873 K

The thermal conductivity λ and the thermal diffusivity a of liquid rubidium were measured by the laser flash method in the temperature interval from the melting point up to 873 K. The measurement error was 4–6%. The data of this paper were compared with the results of other authors. Approximation equations and the table of reference values for the temperature dependence of λ and a have been obtained. The dependence of the Lorentz number on temperature has been calculated.     

Empirical correlation between melting temperature and cohesive energy of binary Laves phases

A list of 143 binary Laves phases with their melting temperature and melting type is collected, and used to study a correlation between melting temperature and cohesive energy. It is found that the melting temperature of Laves phases is roughly proportional to its cohesive energy calculated by Miedema's empirical model from their intrinsic atomic properties. The average predicted error of melting temperature of compounds is as low as 8.0%. This empirical rule is consistent with the result of the universal binding energy theory of solids.     

The melting and solidification of nanowires

A mathematical model is developed to describe the melting of nanowires. The first section of the paper deals with a standard theoretical situation, where the wire melts due to a fixed boundary temperature. This analysis allows us to compare with existing results for the phase change of nanospheres. The equivalent solidification problem is also examined. This shows that solidification is a faster process than melting; this is because the energy transfer occurs primarily through the solid rather than the liquid which is a poorer conductor of heat. This effect competes with the energy required to create new solid surface which acts to slow down the process, but overall conduction dominates. In the second section, we consider a more physically realistic boundary condition, where the phase change occurs due to a heat flux from surrounding material. This removes the singularity in initial melt velocity predicted in previous models of nanoparticle melting. It is shown that even with the highest possible flux the melting time is significantly slower than with a fixed boundary temperature condition.     

Correlations between thermodynamic characteristics of alkali halide crystals and the binding energy of dipolons

The thermodynamic properties of alkali halide crystals are considered. Correlations between the thermodynamic characteristics (such as the melting temperature, the melting energy, and the jump in the entropy upon melting) and the binding energy of dipolons are established for 16 alkali halide crystals.