单轴应变驱动铁bcc—hcp相转变的微观模拟

用分子动力学方法模拟了沿〈001〉晶向应变加载和卸载情况下单晶铁中体心立方(bcc)与六方密排(hcp)结构的相互转变,分析了相变的可逆性和微结构演化特征.微观应力的变化显示样品具有超弹性性质,而温度变化表明在相变和逆相变过程中均出现放热现象.相变起始于爆发式均匀形核,晶核由块状颗粒迅速生长为沿{011}晶面的片状分层结构; 而卸载逆相变则从形核开始就呈现片状形态,且相界面晶面指数与加载相变完全一致,表现出形态记忆效应.在两hcp晶核生长的交界面易形成面心立方(fcc)堆垛层错. fcc通过在hcp晶粒内     

液态金属急冷过程中微观结构转变的研究

本文对液态金属Ａｌ在溶态及急冷过程中的微观结构及其转变进行了分子动力学模拟研究。发现在由ＨＡ（Ｈｏｎｅｙｃｕｔ－Ａｎｄｅｒｓｅｎ）健型指数所确定的双锥体原子团结构的数目随急冷温度变化的关系曲线上存在着二个明显的转变点。第一个点与熟知的玻璃转变温度Ｔｇ相吻合；而第二个点则为新发现的低温端相交点Ｔｇ２。这一结果为用分子动力学方法从微观层次上研究结构相变过程提供了一条新途径。     

一维应变压缩下TC4钛合金绝热剪切研究

 利用内径为57 mm的压缩气炮，在撞击速度为0.2~1.2 km/s（相应的靶中压力为3~15 GPa）范围内进行对称碰撞实验，以研究TC4（Ti-Al6-V4）钛合金在一维应变冲击压缩条件下的绝热剪切现象。对回收得到的受冲击样品，在扫描电镜（SEM）下进行细观金相分析。结果指出，一维应变冲击压缩条件下，TC4钛合金中绝热剪切带产生的对称碰撞速度阈值为500 m/s（相当于样品中的压力为5.87 GPa）；主剪切带与冲击方向约为45°角，带上有圆形和椭圆形两种孔洞且随碰撞速度的增大而增多和长大，这是典型的韧性损伤特征。随碰撞速度增大，产生与主剪切带成15°角的支剪切带。这些与理论预言相符。X射线能谱分析结果指出，剪切带内材料发生了(α+β)→β相的转变，是典型的相变带。剪切带的温度估算与实验提供的信息吻合。     

单晶Cu冲击加载及卸载下塑性行为的微观模拟

使用分子动力学方法对室温下单晶铜沿［001］和［111］方向冲击加载及卸载下的塑性行为进行了模拟,得到了Hugoniot关系以及冲击熔化压力,与实验基本符合. 加载过程中,较高的初始温度有利于位错的形核与发展. 通过对冲击波在自由表面卸载过程的模拟和分析发现：卸载过程呈现“准弹性卸载行为”;沿［001］方向卸载后大量不全位错环与堆积层错消失,而沿［111］方向卸载后只有少量层错消失,部分层错甚至会发展扩大. 关键词： 分子动力学 冲击波 塑性     

液态金属In凝固过程中微观结构转变的模拟研究

采用分子动力学方法对液态金属In的快速凝固过程进行计算机模拟跟踪研究.运用HA键型指 数法和原子成团类型指数法分析了金属原子In的成键类型和形成的原子团簇结构.发现：与通常的液态金属(如Al)相反，随着温度的降低，二十面体及与二十面体相关的1 551 键越来越少；与四方体，六角密堆积相关的1421，1422和1431键数目总和变化很小；而与菱 面体相关的132l，1311，1301和1201的数目却随着温度的降低而显著增加，逐渐占据优势 .最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc，bcc) 关键词： 液态金属In 微结构转变 团簇结构 分子动力学 计算机模拟     

平面一维应变电炮加载技术研究

 在分析电炮加载装置放电回路的基础上，探讨了电炮装置的设计原则。根据所提出的设计原则，设计和建立了一门储能14.4 kJ的电炮装置。实验测试结果表明，其性能参数诸如回路电感小于40 nH、放电周期7 μs、飞片平面度小于40 ns、初始电流上升陡度达750 GA/s。与国内外的同类型装置的参数相比较，该电炮装置达到了较高水平，具有较好的性能。在此基础上，提出了获得平面一维应变加载对电炮装置的一般要求。     

液态金属Ga急冷凝固中微观结构转变的模拟研究

采用分子动力学方法对液态金属Ga的快速凝固过程进行计算机模拟跟踪研究,运用HA键型指数法和原子团类型指数法分析了金属原子Ga的成键类型和形成的基本原子团结构.发现:与通常的液态金属(如Al)相反,随着温度的降低,二十面体及与二十面体相关的1551, 1541键型数目明显下降;与立方体(fcc,bcc),六角密集结构相关的1421,1422键型数目明显增加;而与菱面体相关的1321,1311键型的数目却显著增加,逐渐占据优势.最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc,bcc)、六角密集(hcp)等团簇结构所组成的非晶态结构.这正是非晶态金属Ga的g(r)曲线分裂的第二峰的顺序为前低后高,而与非晶态金属Al的g(r)曲线(其分裂的第二峰的顺序为前高后低)明显不同的微观结构上的物理根源.     

冷却过程中Pd923 团簇的结构转变研究

本文利用分子动力学模拟方法, 研究了钯团簇在不同冷却条件下形成晶体及非晶的过程. 利用平均原子体积、双体分布函数、键对分析和键序参数方法研究了微观局域结构随温度的变化关系. 研究发现：在50 K/ps冷却过程中,液态Pd923团簇在1000 K发生玻璃化转变,在100 K下形成非晶结构; 而在0.1 K/ps冷却过程中,液态Pd923团簇发生结晶,并最终形成六角密排（hcp）晶体结构。     

冷却过程中Pd_(923)团簇的结构转变研究

本文利用分子动力学模拟方法,研究了钯团簇在不同冷却条件下形成晶体及非晶的过程.利用平均原子体积、双体分布函数、键对分析和键序参数方法研究了微观局域结构随温度的变化关系.研究发现:在50 K/ps冷却过程中,液态Pd923团簇在1000 K发生玻璃化转变,在100 K下形成非晶结构;而在0.1K/ps冷却过程中,液态Pd923团簇发生结晶,并最终形成六角密排(hcp)晶体结构.     

高应变率下材料行为的分子动力学研究

 本文用分子动力学方法研究了高应变率下晶体材料的力学行为。在冲击加载下，晶体材料中产生了位错和塑性变形。在强冲击时还可出现相变变化。在讨论应变率变化时，获得了屈服强度随应变率增大而增高的结果。     

Collective modes of one-dimensional lennard-jones systems

The density fluctuations of one-dimensional Lennard-Jones systems are investigated by molecular dynamics simulation. The full Lennard-Jones potential is compared to the repulsive Lennard-Jones potential. It is found that the behavior of the density fluctuations at small wave vectors is determined by the repulsive portion of the potential. The variation of the fluctuations with density is explained. It is shown that these systems do not display hydrodynamics.     

Orientation dependence of structural transition in fcc Al driven under uniaxial compression by atomistic simulations

By molecular dynamics simulations employing an embedded atom method potential,we have investigated structural transformations in single crystal Al caused by uniaxial strain loading along the [001],[011] and [111] directions. We find that the structural transition is strongly dependent on the crystal orientations. The entire structure phase transition only occurs when loading along the [001] direction,and the increased amplitude of temperature for [001] loading is evidently lower than that for other orientations. The morphology evolutions of the structural transition for [011] and [111] loadings are analysed in detail. The results indicate that only 20% of atoms transit to the hcp phase for [011] and [111] loadings,and the appearance of the hcp phase is due to the partial dislocation moving forward on {111} fcc family. For [011] loading,the hcp phase grows to form laminar morphology in four planes,which belong to the {111} fcc family; while for [111] loading,the hcp phase grows into a laminar structure in three planes,which belong to the {111} fcc family except for the (111) plane. In addition,the phase transition is evaluated by using the radial distribution functions.     

The phase transition in a one-dimensional lattice of axisymmetric bodies

The partition function is studied for an array of axisymmetric, hard bodies (capped cylinders, etc.) with each fixed at the base on a regular one-dimensional lattice. It is shown that if a phase transition occurs in a system ofn molecules, then it also occurs in a system of two molecules for the same value of the spacing parameter. With certain additional technical assumptions the converse is also true. Results are reported specifically for a system of thin, hard rods. Necessary and sufficient conditions are shown for a first-order transition to occur in the thermodynamic limit; there is only one transition and that happens when the spacing parameter is equal to the length of the rod. As expected, there is no phase transition when the rod is contracted to a point.     

Atomistic simulation of fcc-bcc phase transition in single crystal Al under uniform compression

By molecular dynamics simulations employing an embedded atom model potential,we investigate the fcc-to-bcc phase transition in single crystal Al,caused by uniform compression.Results show that the fcc structure is unstable when the pressure is over 250 GPa,in reasonable agreement with the calculated value through the density functional theory.The morphology evolution of the structural transition and the corresponding transition mechanism are analysed in detail.The bcc (011) planes are transited from the fcc (11) plane and the (11) plane.We suggest that the transition mechanism consists mainly of compression,shear,slid and rotation of the lattice.In addition,our radial distribution function analysis explicitly indicates the phase transition of Al from fcc phase to bcc structure.     

一维有机导体的Peierls相变研究

根据建立在电子-声子相互作用基础上的Peierls 相变理论，利用形变势模型和半经验晶体轨道方法计算的能带结构数据，对一维有机导体tetrathiafulvalene-tetracyanoquinodimethane的Peierls相变温度进行了计算.并对其金属-绝缘体相变机制进行了讨论. 关键词： 电子-声子相互作用 Peierls相变温度 一维有机导体     

Atomistic simulation of fccben bcc phase transition in single crystal Al under uniform compression

By molecular dynamics simulations employing an embedded atom model potential, we investigate the fcc-to-bcc phase transition in single crystal Al, caused by uniform compression. Results show that the fcc structure is unstable when the pressure is over 250 GPa, in reasonable agreement with the calculated value through the density functional theory. The morphology evolution of the structural transition and the corresponding transition mechanism are analysed in detail. The bcc (011) planes are transited from the fcc (111) plane and the (111) plane. We suggest that the transition mechanism consists mainly of compression, shear, slid and rotation of the lattice. In addition, our radial distribution function analysis explicitly indicates the phase transition of Al from fcc phase to bcc structure.     

应力作用下GdTiO_3磁序相变的第一性原理研究

采用基于密度泛函理论的第一性原理方法,系统研究了GdTiO_3薄膜在压缩应力和拉伸应力作用下的磁序相变.计算结果表明:1)在LaAlO_3压缩衬底的作用下,GdTiO_3薄膜从铁磁基态转变为G型反铁磁基态.该结果不同于YTiO_3和LaTiO_3在LaAlO_3压缩衬底作用时都呈现A型反铁磁基态的情况.若进一步加大压缩应力,例如在(001)平面施加YTiO_3衬底,此时GdTiO_3薄膜基态才为A型反铁磁态.2)在LaScO_3和BaZrO_3拉伸衬底的作用下,GdTiO_3薄膜的基态仍是铁磁态,但是随着拉伸应力的增大,A型反铁磁态的能量和铁磁态的能量差逐渐缩小,即GdTiO_3薄膜的基态有转变为A型反铁磁态的趋势.3)在外加应力的作用下,GdTiO_3薄膜基态的磁序发生了相变,但是其绝缘性并没有变,说明GdTiO_3薄膜仍为Mott型绝缘体.     

Molecular dynamics simulations of phase transition of n-nonadecane under high pressure

The phase transition behavior of n-nonadecane under high pressure was investigated with molecular dynamics (MD) simulations method. A simplified model with amorphous structure and periodic boundary conditions in constant-temperature, constant-pressure ensemble was used in this study. The results showed that the whirling and molecules motion of n-nonadecane chains were restrained by the high pressure. The simulated phase transition temperature of n-nonadecane under high pressure is higher than that under atmospheric pressure. The order parameter of n-nonadecane decreases with the increase in temperature. The simulations reveal that MD is an effective method to understand the phase transition of alkane-based phase change materials on molecular and atomic scale.     

Molecular dynamic studies on materials under laser shocks

Molecular dynamic (MD) simulations offer a powerful means of understanding the microscopic characteristics of shock-propagation through solids and fluids, especially for the short spatial and temporal scales relevant to laser-driven shocks. First-principles molecular dynamics can be directly compared with time-resolved experimental measurements, and methods based on empirical (embedded-atom) potentials fitted to first-principles quantum-mechanical calculations are effective for MD simulations of shock propagation through many millions of atoms. In comparison, thermodynamic approaches based on free-energy considerations do not provide detailed information about mechanical-relaxation or phase-transformation processes within the shock front. We illustrate these ideas by way of embedded-atom simulations of shock-wave propagation through copper crystals of different orientation.