铁磁性物质Ni的液态结构的分子动力学模拟

用分子动力不模拟方法在１８７３－３００Ｋ对液态Ｎｉ的微正则系统进行了模拟研究。模拟采用ＥＡＭ相互作用势对时间和空间的平均,得到了不同温度下Ｎｉ的双体分布函数及原子组态变化的重要信息,并利用键对分析技术对模拟结果作了深入讨论。     

冷却速率对液态金属Cu凝固过程中微观结构演变影响的模拟研究

采用Quantum Sutton-Chen(Q-SC)多体势对液态金属Cu在四个不同冷却速率下的凝固过程进行了分子动力学模拟研究. 通过双体分布函数、键型指数、配位数、均方位移及可视化分析, 结果表明：冷却速率对液态金属Cu的微观结构演变有决定性影响. 当冷却速率为1.0×10¹⁴K/s时得到非晶态结构；当冷速分别为1.0×10¹³K/s，1.0×10¹²K/s和1.3×10¹¹K/s时，系统形成以1421键型为主体的面心立方(fcc)与六角密集(hcp)共存的混合晶体结构；且其结晶温度分别为373K，773K和873K，即冷速越慢，其结晶温度越高，结晶程度也越高；且冷速越慢，1421键型越多，混合晶体中面心立方(fcc)结构所占的比例越高. 同时发现，原子的平均配位数的变化与1551，1441，1661键型的变化密切相关, 反映出体系对称性结构的变化规律与配位数的变化有关. 在可视化分析中，进一步采用中心原子法展现出非晶态与晶体结构的2D截面，及在3D下混合晶体中两个基本原子团分别为面心立方(fcc)与六角密集(hcp)基本原子团的具体结构. 关键词： Q-SC多体势 液态金属Cu 凝固过程 分子动力学模拟     

液态合金NiAl凝固过程中微观结构转变的分子动力学模拟

采用分子动力学模拟方法对液态NiAl凝固过程进行了研究 ,考察了不同冷却速度下液态NiAl结构变化特点 ,原子间相互作用势采用F S多体势 ,结构分析采用键取向序和对分析技术 .计算结果表明 ,冷却速度对液态NiAl结构转变有重要影响 ,在不同的冷却速度下 ,NiAl凝固过程出现了明显不同 ,冷速为 4× 10 13 和4× 10 12 K/s时 ,NiAl快速凝固为无序的非晶体结构 ;而在较慢的 8× 10 11K/s冷速下 ,NiAl凝固为晶态结构 .给出了不同冷却速度下液态NiAl结构转变的微观信息 .     

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

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

液态金属Cu快速凝固过程中原子团簇及非晶态结构转变特性的模拟研究

本文采用Quantum Sutton-Chen多体势,对由5万个液态金属Cu原子组成系统的原子团簇的形成与演变特性进行了分子动力学模似研究.我们采用原子团类型指数法(CTIM)来描述复杂的微观结构转变.研究发现:系统形成以1551、1541和1431三种键型为主的非晶态结构;二十面体原子团(12 0 12 0)和(12 2 8 2)、(13 1 10 2)、(13 3 6 4)、(14 1 10 3)、(14 2 8 4)、(14 3 6 5)缺陷多面体基本原子团在液态转变为非晶体过程中起着关键性的作用.系统所形成的纳米团簇是由一些基本团簇和由这些基本团簇相互连接而成的中等团簇所组成,这正是与由气相沉积法和离子溅射法所获得的团簇结构的本质差别所在.通过双体分布函数g(r)、HA键型、基本原子团、平均原子体积和比值g_(min)/g_(max)的分析,还得到液态金属Cu在冷却速率为1.0×10~(14)K/S时的非晶转变温度T_8约为673 K.同时还发现,1551、1441、1661三种键型随温度有相同的变化趋势,这反映出体系对称性结构有相同的变化规律.     

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

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

冷速对液态金属Ga凝固过程中微观结构演变影响的模拟研究

采用分子动力学方法对液态金属Ga凝固过程中不同冷却速率对微观结构演变的影响进行了模拟跟踪研究. 运用HA键型指数法和原子成团类型指数法(CTIM)分析了金属原子Ga的成键类型和形成的基本原子团结构. 结果发现，冷却速率对凝固过程中的微观结构起着非常重要的作用. 在以1.0×10¹⁴，1.0×10¹³，1.0×10¹²K/s的速率冷却时，系统形成以与1311，1301键型相关的菱面体结构为主，夹杂着立方体、六角密集等其他团簇结构所构成的非晶态结构；在以1.0×10¹¹K/s的速率冷却时，系统明显发生结晶，其结晶转变温度T_c约为198K，且冷却速率越慢，结晶转变温度T_c越高，形成以与1421键型相关的斜方晶体(经可视化分析确认)为主要构型的晶态结构. 这将为研究液态金属的结晶转变过程提供一种新方法. 关键词： 液态金属Ga 凝固过程 微结构转变 分子动力学模拟     

冷速对液态金属Na凝固过程中微观结构影响的模拟研究

采用分子动力学方法对液态金属Na在四种不同冷速下的快速凝固过程进行了模拟跟踪研究.采用双体分布函数g(r)曲线、Honeycutt-Andersen键型指数法和原子团类型指数法对凝固过程中微观结构的变化进行了分析.结果表明：冷却速率对微结构的转变有决定性影响，当冷速为1.0×10¹⁴和1.0×10¹³K/s时，系统形成以1551和1541键型或以缺陷多面体基本原子团(13 1 10 2)和二十面体基本原子团(12 0 12 0)为主体的非晶态结构；当冷速为1.0×10¹²和1.0×10¹¹K/s时，系统则形成以1441和1661键型或以体心立方基本原子团(14 6 0 8)为主体的晶态结构.同时发现：不同冷速对液态金属Na在液态和过冷态时微观结构的影响甚小；但不同冷速对其固态(非晶态利晶态)时的微观结构有显著的影响，且要在液-固转变点(分别在玻璃转变温度T_g和晶化起始温度T_c)附近或以后才能充分展现出来.根据这一特点，有可能建立另一种确定液态金属T_g和T_c的新方法.原子团类型指数法比键型指数法更有利于研究液态、非晶态等无序体系和一些晶化体系的具体结构特征. 关键词： 液态金属Na 凝固过程 分子动力学模拟 原子团类型指数法     

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

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

单壁碳纳米管端口结构的热稳定性

利用紧束缚势分子动力学模拟方法,研究了温度在1000 K-3000 K之间单壁碳纳米管端口结构的变化趋势.计算表明,温度对整个管端口结构有重要影响,温度升高容易使理想单壁碳管端口封闭.温度在3000 K下碳管端口达到封闭,而且端口封闭导致碳管系统能量的降低.碳纳米管长度越长,端口封闭越快,且扶手型碳纳米管比锯齿型碳纳米管更容易形成端口封闭的结构.     

液态金属的初始状态对凝固微结构影响的模拟研究

采用分子动力学方法对液态金属Ａｌ在不同的初始状态下，以相同急冷速率凝固的过程进行模拟跟踪研究，发现：在玻璃转变温度Ｔｇ以上（即过冷液态）时，系统的微结构组态情况基本一致，相差甚微；但在Ｔｇ以下时，不同的初始液态微结构对其固态微结构有明显的影响；且在Ｔｇ处各种固态微结构之间的差别发生突发性的变化。这一结果对于深入理解液－固微结构之间的转变关系，具有一定的重要意义。     

热镀铜层冷凝过程的分子动力学模拟

采用常温、常压分子动力学模拟方法和FS(Finnis Sinclair)势 ,研究了在周期性边界条件下由 5 0 0个原子构成的液态Cu模型系统的凝固过程 ,考察了不同降温速率下Cu的凝固行为 ,得到了不同温度、不同冷却速率下Cu的双体分布函数 ;采用HA键型指数法统计了各种小原子团在不同温度下所占比例 ,采用键取向序分析了体系降温全过程的局域取向对称性 ,得到原子体系微观结构组态变化的重要信息 ;最后 ,利用能量分析的方法对体系微观结构的变化进行了说明 ,给出了液态Cu冷凝过程中微观结构转变的重要信息 .     

第一过渡金属酞菁分子的电子结构的第一性原理计算

对第一过渡金属酞菁化合物(Metal Phthalocyanine,MPc,M=Sc,Ti,V,Cr,Mn,Fe,Co,Ni和Cu)的电子结构和基本物理化学性质进行了第一性原理计算.理论模拟出来的STM图像表现出亚分子结构,与已有的实验观察结果相当吻合,且跟金属原子的d电子组态明显有关.在小偏压条件下,第一过渡金属首尾端ScPc,NiPc和CuPc分子的中央金属离子在STM图像表现为空洞,其他所有金属酞菁分子的中央金属离子均为亮斑.同时还研究了ScPc和NiPc分子的STM图像与偏压的关系,当针尖偏压分别 关键词： STM图像模拟 金属酞菁 电子结构     

液态Au和Ag冷却过程的分子动力学模拟

采用反映原子多体相互作用的FS势模型,通过分子动力学方法研究了贵金属Au、Ag在温度2000~300 K之间的冷却过程,考察其结构变化特点,给出了冷却过程微观结构转变的重要信息.     

贵金属Au冷却过程中结构及能量变化的分子动力学计算机模拟

通过分子动力学方法,研究了不同冷速下贵金属Ａｕ在温度２０００～３００Ｋ的冷却过程中微观结构的变化特点,结果发现,冷却速度对Ａｕ的微观结构产生重要影响,采用偶关联函数和键对分析技术对原子局域团簇结果者分析,并考察了冷却过程中原子势能随温度的变化,比较了Ａｕ的微观结构转变与能量变化的对应关系,从能量转化的角度对冷却过程中Ａｕ的结构变化进行了说明。     

耐蚀合金Au3Cu高温冷却过程中能量及结构转变的分子动力学模拟

用分子动学模拟方法对液态Au3Cu冷却过程进行了研究,考察了不同冷却速度下Au3Cu结构变化特点,原子间相互作用势采用F-S多体势,结构分析采用键取向序和对分析技术。计算结果表明,冷却速度对液态Au3Cu能量及结构转变有重要影响,给出了不同冷却速度下液态Au3Cu结构转变的微观信息。     

金属线性热膨胀系数的研究

 从温度T时的压强表达式出发，结合普适状态方程（UEOS）、Debye模型和热力学关系式，给出了计算金属线性热膨胀系数的普遍公式。并对Al、Cu、Pb的这一物理量进行了具体计算。结果表明，理论与实验符合较好。     

微通道内电解质溶液离子分布动电学效应模拟

动电学效应对微通道内流体流动特性影响很大,其对通道内粒子分布特性的影响使得通道近壁面流体流动特性极不稳定。本文采用分子动力学方法模拟了二维矩形微通道内NaCl稀电解质溶液的流动特性,考虑存在于不同粒子间的Lennard-Jones势能、静电力、以及带电离子与水分子间的相互作用,得到了粒子在通道内的分布特征。结果显示在动电学效应主要作用于通道壁面附近,而主流区域影响极小。Na~+离子在无量纲通道高度达到0.08和0.91时其浓度达到最大值,沿远离壁面其浓度逐渐降低,与壁面电性相反的Cl~-离子则在无量纲通道高度达到0.15和0.84附近浓度最高。其结果与基于连续介质解理论的Boltzamnn统计分布一致。水分子的浓度在壁面附近也较通道中心处高。