自由能方法与零压下Al的熔化温度

利用自由能方法的分子动力学模拟,计算了零压下Al的熔化温度.在计算液相自由能的过程中,采用勒纳-琼斯(LJ)液体作为参考系统,同时将计算结果与Mei和Davenport等人的计算结果进行了比较,计算结果表明：1)选用LJ参考系统使液相自由能的计算时间节省一半,并且不影响熔化温度的计算结果;2)采用不同的埋入原子势(EAM)的分子动力学模拟计算得到的熔化温度与实验值都存在偏差,而就金属Al而言,采用Cai等人的EAM势的熔化温度的计算结果比Mei和Davenport及Morris等人采用的势模型的结果略有改 关键词： 熔化温度 自由能方法 分子动力学模拟     

金属熔体黏度与结构相关性的分子动力学模拟

运用EAM(embed atom method)作用势，采用非平衡分子动力学模拟获得Al熔体的偶分布函数与黏度数值随温度的变化曲线，偶分布函数的计算结果与实验值符合得较好.对模拟所得到的黏度数据编程实现黏度的Arrhenius公式拟合，得到激活能E.并利用模拟所得到的黏度值及激活能对Lennard-Jones(L-J)作用势进行修正，获得黏度与偶分布函数及原子间相互作用势之间的关系式，两条黏度拟合曲线与分子动力模拟结果符合得比较好，说明拟合程序的编写是比较成功的，实现了对L-J作用势的修正.该研究为金属及合金原子间相互作用势的建立提供了新的思路. 关键词： 非平衡分子动力学模拟 L-J作用势修正 Al熔体 结构与黏度相关性     

Aln(n=13–32)团簇熔化行为的分子动力学模拟研究

基于半经验的Gupta原子间多体相互作用势, 采用分子动力学方法并结合模拟退火及淬火技术, 系统研究了小尺寸铝团簇Al_n (n=13–32)的熔化行为. 模拟结果表明: 除个别尺寸(Al₁₃ 和Al₁₉) 外, 团簇熔化过程比热曲线普遍呈现杂乱无规(无明显单峰)现象, 这与实验观测小Al团簇比热普遍无规的结果完全一致. 通过分析不同温度点上团簇淬火构型的势能分布图给出了小Al团簇比热呈现无规或有规现象的成因. 对于比热无规团簇, 可以利用原子等价指数判断给出团簇熔点, 所得团簇熔点随团簇尺寸的变化趋势与实验观测结果完全一致. 关键词： Gupta势 团簇 分子动力学 熔化     

LJ模型的热力学特性、相平衡及用于真实分子的研究

结合描述硬球固体Helmholtz自由能的自由体积方法与描述硬球固体径向分布函数的拟合的分析表达式与一阶热力学摄动理论,用于描述Lennard Jones(LJ)固体的Helmholtz自由能.按照一个修正的WCA方法将LJ势分为短程排斥部分与长程吸引部分,将文献中一个用于求取液相的等价的硬球直径的简单的迭代法扩展到固相,用于求取固相的等价的硬球直径.在固体Helmholtz自由能的计算中,使用200壳层,以便获得精确的结果.体相LJ液体的热力学特性由一个最近提出的状态方程求取.该方法很好地描述了LJ固体的过量Helmholtz自由能与状态方程,满意地描述了Lennard Jones模型的相平衡;通过选取合适的LJ势参数,能很好地描述了真实分子的融化曲线.     

液态Co物性的平衡及非平衡分子动力学模拟

本文采用EAM作用势,通过平衡分子动力学(EMD)模拟的方法计算了Co熔体的自扩散系数、剪切粘度等物理性质.同时采用非平衡分子动力学(NEMD)方法计算了Co的剪切粘度.研究表明有关传输性质的计算是可与实验比较的,能够反映出液态Co典型的动力学特性.     

Pt-Ag合金纳米线热力学性质的分子动力学模拟研究

基于原子嵌入势(EAM),采用分子动力学方法,对临界尺寸下的Pt_(0.95)Ag_(0.05)合金纳米线多边形结构的熔化行为进行了计算模拟.结果表明:径向尺寸对Pt_(0.95)Ag_(0.05)合金纳米线的熔点影响较为显著,而长度对其影响较小;引入林德曼因子得到的熔点和用势能-温度变化曲线找到的熔点基本一致;合金纳米线的染色原子由外向内运动;综合分析发现Pt_(0.95)Ag_(0.05)合金纳米线以先外后内的模式进行熔化.     

金属钨熔化过程的分子动力学模拟

本文使用分子动力学方法对金属钨的熔化过程进行了数值模拟,分析了钨在熔化过程中的结构、系统内能变化以及表面熔化过程固-液界面变化情况,初步分析了表面熔化现象的机理。模拟过程采用嵌入原子模型(EAM)描述原子间相互作用,模拟结果表明,嵌入原子模型适合于计算固-液相变过程,表面熔化过程是由表面处最外层原子的不稳定性触发的。对于均匀熔化过程,晶体在4700 K下发生固-液相变;对于表面熔化过程,计算获得了不同温度(3800～4800 K)下的熔化速度,拟合出熔化速度公式,得到的表面熔化热力学熔点与已有实验结果基本符合。     

金属Al表面熔化各向异性的分子动力学模拟

采用嵌入原子势,使用分子动力学方法对金属Al不同低指数晶面的表面熔化现象分别进行了模拟.分析了熔化过程中样品结构组态的变化.模拟结果表明对于不同的自由表面,表面熔化呈现出明显的各向异性行为.Al（110）面在低于熔点的温度之下发生预熔化;（111）与（001）面都出现过热现象.与（111）面不同,（001）面发生过热现象时表面原子层为类液层,而（111）面仍然保持很好的晶格结构,即预熔化的Al（001）面在高于熔点的温度下,仍可以在很长的时间内处于相对稳定的亚稳态.由模拟得到Al的热力学熔点为950 K左右,与实验值基本符合. 关键词： 分子动力学 表面熔化 过热     

Ni_(50)Al_(50)合金的液固转变及晶体生长过程的分子动力学模拟

采用分子动力学模拟技术研究了液态Ni50 Al50 合金在不同冷速下的凝固特点 ,模拟采用EAM作用势 ,计算了不同温度 ,不同冷速下Ni50 Al50 的偶分布函数。结果表明EAM作用势能很好地描述液态Ni50 Al50 的无序结构 ,在快速凝固条件下 ,液态Ni50 Al50 形成非晶 ,当冷速较慢时 ,液态Ni50 Al50 形成晶体 ,分析了不同冷速下体系的相变热力学及相变动力学特点。最后采用液固两层构型法 ,清楚地观察到Ni50 Al50 晶体生长的全过程。     

Ni50Al50合金的液固转变及晶体生长过程的分子动力学模拟

采用分子动力学模拟技术研究了液态Ni50Al50合金在不同 冷速下的凝固特点,模拟采用EAM作用势,计算了不同温度,不同冷速下Ni50Al5 0的偶分布函数。结果表明EAM作用势能很好地描述液态Ni50Al50的无序结构 ,在快速凝固条件下,液态Ni50Al50形成非晶,当冷速较慢时,液态Ni50 Al50形成晶体,分析了不同冷速下体系的相变热力学及相变动力学特点。最后采 用液固两层构型法,清楚地观察到Ni50Al50晶体生长的全过程。     

关于"太阳在银河系中所受引力与到银心距离r的关系"的讨论

讨论了银道面内的引力场强分布和太阳在银河系中所受引力与到银心距离r的关系，指出引力场强g并不是与r的平方成反比，这是由于银河系的大小和形状不能忽略造成的．强调了万有引力定律的适用条件．     

On the predictability of the positions of the convergence zones in the ocean

The disagreement between the experimental and calculated positions of the first convergence zone are known from many publications. The most probable cause for such a disagreement, namely, the incorrect specification of the input data for the calculations, is considered. The lack of simultaneity between the hydrological surveys of the region and the acoustic experiments is emphasized. The experimental data obtained by the author in five ocean regions are presented. These data characterize the diurnal variability of the distance from the source to the nearest boundary of the convergence zone. The relations proposed by different researchers for calculating the sound speed from the temperature, salinity, and hydrostatic pressure are analyzed. It is shown that these relations lead to a substantial difference in the estimated depth dependence of the hydrostatic gradient of the sound speed. The position of the first convergence zone is calculated for the propagation conditions determined by vertical temperature and salinity profiles with the subsequent recalculation of these profiles into sound speed profiles by using eight different formulas known from the literature. It is shown that different formulas lead to different values of the distance to the first zone; this difference is substantially greater than that between the calculations and experiment. The necessity of improving the recalculation relations in view of the experimental data on sound propagation in natural oceanic waveguides, including the data on the actual positions of the convergence zones, is emphasized.     

On the calculation of the bremsstrahlung in the region of the plasma frequency

On the basis of the classical theory of radiation produced in incomplete Coulomb interaction, a formula is derived for the coefficient of the bremsstrahlung of fully ionized plasma. In the region of rectilinear approximation a comparison is made with the results of the quantum formula derived from the Born approximation. Furthermore, the influence of dispersion upon the emission coefficient is taken into account; some relations for the region of dispersion are given.     

Dependence of the Intensity of the Spectral Lines and the Composition of the Sample and the Plasma Parameters of the Unipolar Arc

The investigation was made of the dependence of the intensity of Tl, Ga, Mo, Mg, Mn, Sn, Bi, Ni, Zn, Pt and Au spectral lines and the plasma parameters (temperature T, electron concentration n_e, degree of 6 ionnization α) from concentration of lithium additive.     

Extracting the wavefunction of the LSP at the LHC

We consider associated production of squarks and gluinos with the lightest supersymmetric particle (LSP), or states nearly degenerate in mass with it. Though sub-dominant to pair production of color SU(3)-charged superpartners, these processes are directly sensitive to the wavefunction composition of the lightest neutralinos. Exploiting event-shape variables - including some introduced here for the first time - we are able to identify the composition of the LSP by selecting events involving a single high-pT jet recoiling against missing transverse energy. We illustrate the proposed technique on a set of benchmark cases and propose methods for applying these results in more realistic experimental environments.     

Set the controls for the heart of the Sun

In this paper we describe experiments conducted with high-power lasers that are attempting to replicate, for a very short time and in miniature, conditions found in the Sun. Experiments to date have reached conditions in the outer part of the Sun. To reach the Sun's centre requires compression of material to very much greater than solid density and heating to over ten million degrees. To achieve this, a new class of experiments and a new generation of high-power lasers are required.     

On the Question of the Source of the Apokamp

Technical Physics - The object of this work is the apokamp—a new type of plasma jet, which is formed from a bright offshoot emerging at the bending point of a channel of a high-voltage...