自由能方法与零压下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熔体 结构与黏度相关性     

金属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 晶体生长的全过程。     

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

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

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

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

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

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

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

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

硅熔化特性的分子动力学模拟–-不同势函数的对比研究

分别采用Stillinger-Weber (SW)势、修正的成熟原子嵌入模型(MEAM)势、 Tersoff势和HOEP (highly optimized empirical potential)势来描述硅原子间相互作用, 运用分子动力学方法对比模拟研究了四种势函数的硅晶体的体熔化和表面熔化特性. 结果表明: 四种势函数均能反映出硅的热膨胀、高温熔化和熔化时吸热收缩等基本物理规律. 但综合对比发现, Tersoff势和MEAM势相对更适合描述硅的熔化和凝固过程, SW势次之, HOEP势则不适合描述硅的熔化和凝固过程. 关键词： 硅 势函数 熔化 分子动力学     

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

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

Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid Al and Cu

We investigate the application of embedded atom method (EAM) interatomic potentials in the study of crystallization kinetics from deeply undercooled melts, focusing on the fcc metals Al and Cu. For this application, it is important that the EAM potential accurately reproduces melting properties and liquid structure, in addition to the crystalline properties most commonly fit in its development. To test the accuracy of previously published EAM potentials and to guide the development of new potential in this work, first-principles calculations have been performed and new experimental measurements of the Al and Cu liquid structure factors have been undertaken by X-ray diffraction. We demonstrate that the previously published EAM potentials predict a liquid structure that is too strongly ordered relative to measured diffraction data. We develop new EAM potentials for Al and Cu to improve the agreement with the first-principles and measured liquid diffraction data. Furthermore, we calculate liquid-phase diffusivities and find that this quantity correlates well with the liquid structure. Finally, we perform molecular dynamics simulations of crystal nucleation from the melt during quenching at constant cooling rate. We find that EAM potentials, which predict the same zero-temperature crystal properties but different liquid structures, can lead to quite different crystallization kinetics. More interestingly, we find that two potentials predicting very similar equilibrium solid and liquid properties can still produce very different crystallization kinetics under far-from-equilibrium conditions characteristic of the rapid quenching simulations employed here.     

Experimental and theoretical investigations of the brittle fracture of Ti<Subscript>3</Subscript>Al

The formation of microcracks upon dislocation interactions was studied by means of TEM analysis. The types of microcracks were classified depending on the orientation of deformation axes in single crystals. The Griffith surface energy was simulated, and the energy of unstable stacking faults was calculated by the molecular dynamics method using the N-particle EAM potentials of the interatomic interaction for Ti₃Al. The Rice—Thompson model was used to study the relationship between the tendency toward cleavage and the plastic relaxation of stresses near a crack tip by the emission of dislocations in Ti₃Al single crystals.     

A molecular dynamics study on iridium

In this study, molecular dynamics simulations are performed by using a modified form of Morse potential function in the framework of the Embedded Atom Method (EAM). Temperature-and pressure-dependent behaviours of bulk modulus, second-order elastic constants (SOEC), and the linear-thermal expansion coefficient is calculated and compared with the available experimental data. The melting temperature is estimated from 3 different plots. The obtained results are in agreement with the available experimental findings for iridium.      

Melting curves of FCC-metals by cell-theory

The cell model is developed to account for many body interaction effects, as occurring in the embedded atom method (EAM) potentials, and crystal lattice features to determine the free energy of FCC metals. The well known smearing approximation, which is generally used in conjunction with pair potentials, is also developed for EAM potentials. The free energy so obtained is used to determine melting curves of FCC metals. For this purpose, the liquid phase free energy is calculated using the corrected rigid spheres model. The advantage of our scheme, which is verified with data on Lennard-Jones solids, is that both free energies are based on the same interaction potential. Results match well with the available experimental/theoretical data for Al and Cu. A good agreement is also found for the transition metals, Ni and Pt, for which molecular dynamics as well as theories like Lindemann's law do not give accurate results. It is also found that smearing approximation, which neglects interactions beyond nearest neighbors, also provides good estimates for free energy if many body effects are accounted with EAM potential.     

Molecular dynamics simulations of the melting of Al–Ni nanowires

Molecular dynamics (MD) simulations were performed to investigate the influence of nickel (Ni) composition and nanowire thickness on the thermal properties of Al-x%Ni (at%) nanowires using the embedded atom model (EAM) potential. The melting of the nanowire was characterised by studying the temperature dependence of the cohesive energy and mean square displacement. The effect of the nanowire thickness on the cohesive energy, melting temperature, heat capacity as well as latent heat was studied in canonical ensemble. Moreover, the crystal stability of Al, Al-20%Ni, Al-40%Ni, Al-60%Ni, Al-80%Ni, Al₃Ni, Ni₃Al and Ni nanowires was studied at different temperatures using mean square displacement and cohesive energy.     

A molecular dynamics study on melting point and specific heat of Ni<Subscript>3</Subscript>Al alloy

Using the Embedding Atom Method (EAM) for highly undercooled Ni₃Al alloy, the melting point and the specific heat were studied by a molecular dynamics simulation. The simulation of melting point was carried out by means of the sandwich method and the NVE ensemble method, and the results show a good agreement, whereas are larger than the experimental value of 1663 K. This difference is attributed to the influence of surface melting on experimental results, which causes the smaller measurements compared with the thermodynamic melting point. The simulated specific heat of Ni₃Al alloy weakly and linearly increases with the increase of undercooling in the temperature range from 800 K to 2000 K. Supported by the National Natural Science Foundation of China (Grant No. 50395101)     

Thermodynamic and physical properties of FeAl and Fe3Al: an atomistic study by EAM simulation

With this work we present a newly developed potential for the Fe–Al system, which is based on the analytical embedded atom method (EAM) with long range atomic interactions. The potential yields for the two most relevant phases B2-FeAl and D0₃-Fe₃Al lattice constants, elastic constants, as well as bulk and point defect formation enthalpies, which are in good agreement with experimental and other theoretical data. In addition, the phonon dispersions for B2-FeAl and D0₃-Fe₃Al show a good agreement with available experiments. The calculated lattice constants and formation enthalpy for disordered Fe–Al alloys are in good agreement with experimental data or other theoretical calculations. This indicates that the present EAM potentials of Fe–Al system is suitable for atomistic simulations of structural and kinetic properties for the Fe–Al system.     

On the suitability of the ILJ function to match different formulations of the electrostatic potential for water-water interactions

Energetic and structural properties of liquid water have been calculated using molecular dynamics simulations in order to investigate the effect of different formulations of the van der Waals (vdW) interaction on the behaviour of liquid water. In particular, two model potentials, the SPC/E using a Lennard Jones (LJ) function and the AMPF using an Improved Lennard Jones function (ILJ) have been considered. The flexibility of the ILJ function in the AMPF model has also been analysed, proving that its vdW component can match different parametrizations of the electrostatic component.     

The fundamental vibrational shift of HCl diluted in dense Ar and Kr

The fundamental infrared band vibrational shifts of HCl diluted in dense Ar and Kr have been calculated by means of molecular dynamics simulation techniques and compared with the experimental data reported by Pérez et al. (J. Chem. Phys. 122, 194507 (2005)). The results have been analysed along the liquid–vapour coexistence line as in several supercritical states, with good agreement between the theoretical and experimental values.     

Influence of Defects and Crystallographic Orientation on Mechanical Behavior of Nanocrystalline Aluminium

Simulation of molecular dynamics using Embedded Atom Method (EAM) potentials is performed to investigate the mechanical properties of single crystal Al along various crystallographic orientations under tensile loading. The specimens are provided with one or two embedded circular voids to analyze the damage evolution by void growth and coalescence. The simulation result shows that the Young's modulus, yielding stress and ultimate stress decrease with the emergence of the voids. Besides, the simulations show that the single-crystal Al in different crystallographic orientations behaves differently in elongation deformations. The single-crystal Al with <100> crystallographic orientations has greater ductility than other orientated specimens. The incipient plastic deformation and the stress-strain curves are presented and discussed for further understanding of the mechanical properties of single-crystal Al.