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

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

液态金属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 凝固过程 微结构转变 分子动力学模拟     

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

采用分子动力学方法对六种不同冷却速率对液态金属Pb凝固过程中微观团簇结构演变的影响进行了模拟跟踪研究.采用双体分布函数、Honeycutt-Andersen键型指数法、原子团类型指数法(CTIM-2)、平均配位数等方法对凝固过程中微观团簇结构的演变进行了分析.结果表明：系统存在一个形成非晶态或晶态结构的临界冷速(介于5×10¹²与1×10¹²之间),大于这个临界冷速时,系统将形成以1551,1541和1431键型为主的非晶态结构;当小于这个临界冷速时,系统将先形成以1441和1661键型或以bcc基本原子团(14 6 0 8 0 0)为主的晶态结构,并稳定存在一段时间,然后又迅速转变为以1421和1422键型为主或以fcc基本原子团(12 0 0 0 12 0)和hcp基本原子团(12 0 0 0 6 6)以一定比例并存的部分晶态结构;同时发现,冷速对系统中fcc结构和hcp结构的相对比例有明显的影响,冷速越低,fcc结构所占的比例越多,越倾向于形成完美的fcc晶态结构. 关键词： 液态金属铅 凝固过程 团簇结构演变 分子动力学模拟     

冷速对液态金属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 凝固过程 分子动力学模拟 原子团类型指数法     

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

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

冷速对液态金属Mg凝固过程中微观结构演变的影响

采用分子动力学方法对不同冷速下液态金属镁(Mg)快速凝固过程中的微观结构演变进行了模拟研究.并采用能量-温度(E-T)曲线、双体分布函数、Honeycutt-Andersen键型指数法、原子团簇类型指数法(CTIM-3)以及三维可视化等方法系统地考察了凝固过程中微观结构演变与相转变过程.结果发现:在以冷速为1×10~(11)K/s的凝固过程中,亚稳态bcc相优先形成,随后大量解体,其变化规律符合Ostwald规则,系统最终形成以hcp结构为主体与fcc结构共存,中间还夹杂部分bcc结构的致密晶体结构.在1×10~(12)K/s冷速下,结晶过程呈现迟缓现象,形成bcc结构的初始温度降低,系统形成以hcp居多、与bcc和fcc三相共存的结构,且因相互竞争、相互制约而导致不易形成粗大的晶粒结构.而在1×10~(13)K/s冷速下,系统则形成以1551,1541,1431键型为主的多种非晶态基本原子团组成的非晶态结构.此外,在冷速1×10~(12)与1×10~(13)K/s之间的确存在一个形成非晶态结构的临界冷速.     

熔体初始温度对液态金属Na凝固过程中微观结构影响的模拟研究

采用分子动力学方法对熔体初始温度热历史条件对液态金属Na凝固过程中微观结构的影响，进行了模拟研究，并采用双体分布函数g(r)曲线、键型指数法和原子团类型指数法对凝固过程中的微观结构进行了分析.结果表明：液态金属Na在不同熔体初始温度条件下以1×10¹¹K/s冷速凝固时，均形成晶化结构，其中1661和1441键型或体心立方基本原子团(14 6 0 8)在凝固过程中对微观结构的转变起决定性作用.同时发现：熔体初始温度对凝固微结构有显著影响，而对液态和过冷态的微观结构影响并不明显，只有在晶化起始温度T_c附近才充分地展现出来.不同熔体初始温度对凝固结构的晶化程度有不同的影响，虽其影响程度是随着熔体初始温度的下降呈非线性变化关系的，但仍表明是可以通过改变熔体初始温度来加以控制的.原子团类型指数法(比键型指数法)更进一步表征了晶化体系中原子团的结构特征，将有利于对液态金属凝固过程中微观结构的转变机理进行更为深入的研究.     

冷却速率对液态金属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 凝固过程 分子动力学模拟     

熔体初始温度对液态金属Ni凝固过程中微观结构演变影响的模拟研究

采用量子 Sutton-Chen多体势, 对熔体初始温度热历史条件对液态金属Ni快速凝固过程中微观结构演变的影响进行了分子动力学模拟研究. 采用双体分布函数g(r)曲线、键型指数法、原子团类型指数法和三维可视化等分析方法对凝固过程中微观结构的演变进行了分析. 结果表明: 熔体初始温度对凝固微结构有显著影响, 但在液态和过冷态时的影响并不明显, 只有在结晶转变温度T_c附近才开始充分显现出来. 体系在1×10¹² K/s的冷速下, 最终均形成以1421和1422键型或面心立方(12 0 0 0 12 0)与六角密集(12 0 0 0 6 6) 基本原子团为主的晶态结构. 末态时, 不同初始温度体系中的主要键型和团簇的数目有很大的变化范围, 且与熔体初始温度的高低呈非线性变化关系. 然而, 体系能量随初始温度呈线性变化关系, 初始温度越高, 末态能量越低, 其晶化程度越高. 通过三维可视化分析进一步发现, 在初始温度较高的体系中, 同类团簇结构的原子出现明显的分层聚集现象, 随着初始温度的下降, 这种分层现象将被弥散开去. 可视化分析将更有助于对凝固过程中微观结构演变进行更为深入的研究. 关键词： 液态金属Ni 熔体初始温度 微观结构 分子动力学模拟     

金属Zn液态结构变化的研究

利用ＴＢ模型给出的原子间相互作用势详细计算了不同温度下Ｚｎ的双体分布函数ｇ（ｒ），结果发现随着温度的不断降低，液态金属Ｚｎ的ｇ（ｒ）第一峰变得高而尖，第二峰由弱变强，说明了液态金属Ｚｎ的有序度随温度降低而不断增强；利用键对分析技术统计了液态金属Ｚｎ在不同温度下的键取向序参数、键对数。键取向序参数及键对数随温度的变化，进一步证明了低温液态的有序度高于高温液态，从而充分说明液态金属在不同温度下有不同的结构形式，而不像人们想象得那样杂乱无章。     

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

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

冷却速度对Ga凝固过程的影响的模拟研究

采用分子动力学方法对Ga在相同初始状态下以不同速度冷却的凝固过程中进行了模拟研究。发现：以3.38×1013、3.38×1012K/s的速度冷却,得到非晶态结构;以2.01×1011K/s的速度冷却,发生明显晶化,结晶转变温度约为133K。这一结果,对于如何正确选择冷却速度获得优良材料性能,将具有重要的实际意义。     

HL—1装置边缘分子过程对中性粒子输运的影响

本文采用新版的一维中性粒子输运程序,在平板几何下模拟了HL-1装置边缘8cm区域的分子过程及中性粒子与壁相互作用对中性粒子输运的影响,得到了电离速率的径向分布、稳态中性氢原子径向分布、壁与边缘等离子体区之间、芯部等离子体与边缘等离子体之间中性氢各种粒子通量、能量通量并与只含原子过程的类似计算进行了比较,结果表明边缘区的分子过程对巾性氢原子输运有显著的影响。有些结果得到实验的证实。     

固体中冲击波的分子动力学研究

 文中用分子动力学方法研究了具有面心结构的晶体在冲击波的作用下温度、压力、粒子速度、物质密度及冲击波速度诸参数的变化,并与宏观规律作了比较。文中采用Fincham的计算格式。原子间的相互作用势用Morse势,其参数用最近邻相互作用近似获得。     

二维超导薄膜中拓扑缺陷数目与冷却速度关系研究

基于含时Ginzburg-Landau（TDGL）方程建立二维超导薄膜经历快速冷却后的拓扑缺陷数目与冷却速度的关系模型.采用差分法对此模型进行了数值模拟,计算出不同冷却速度下拓扑缺陷数目随着时间的变化关系.计算结果发现,在不同冷却速度下,缺陷数目随着时间的变化可以分为两个阶段：第一阶段缺陷数目较多,拓扑缺陷数目随时间急...     

GaAs中过渡金属深杂质能级的计算

本文采用41个原子的原子族模型来模拟晶体,用电荷-组态自洽(SCCC)的EHMO方法和在边界上用"类Ga"和"类As"原子来饱和悬挂键,计算了中性过渡金属原子Cr、Mn、Fe、Co、Ni在半导体GaAs中的深杂质态,所得杂质能级在带隙中的位置与实验结果相符。     

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

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

二元液体中动力学不均匀性与脆性的关系

过冷液体和玻璃有着复杂的结构和动力学性质.目前我们对他的研究是很不够的.当液体的温度远低于它的熔点时,它呈现非指数的驰豫行为.大量的实验和分子动力学模拟的研究表明这种非指数的驰豫行为是由动力学不均匀性引起的.计算机模拟已经提供了一些过冷液体的动力学不均匀性的相关信息.而Agell认为,液体可以粘度随温度的变化分为强,弱液体两种.脆性表示的是随着玻璃态的接近液体性质(如黏度)的改变速度.在这篇文章中我们通过模拟了三种不同二元Lennard-Jones液体,计算了它们的非高斯参数,脆性系数等一些参数,研究动了力学不均匀性和脆性的关系.得出了系统的脆性系数越小他的结构就越稳定,动力学不均匀性就越强的结论.     

微小化毛细泵吸环路(miniature CPL)应用于笔记本计算机传热之研究

本文针对目前笔记本计算机CPU在小空间内的高发热量趋势,提出采用毛细泵吸环路（CPL）取代传统热管,实现小空间内的更高效传热。传统大型CPL环路的研究与应用已经相当完备,但当尺寸缩小时会有许多新的问题产生。本文通过实验观察及理论分析,望能发展出适用未来电子产品高发热需求的微小化CPL环路。