H+2团簇与固体相互作用的Monte Carlo模拟

讨论了H+2团簇在物质中的作用过程, 包括库仑爆炸过程、尾流势的影响和与靶原子的近碰撞等. 开发成功了一套适用于不同靶材料的Monte Carlo模拟程序, 模拟计算了H+2通过碳膜后出射粒子的能谱和角分布, 并与实验结果进行了比较.     

超薄膜外延生长的Monte Carlo模拟

用MonteCarlo(MC)方法对超薄膜外延生长过程进行了计算机模拟.模型中引入Morse势描述粒子间的相互作用,考虑粒子的沉积、吸附粒子的扩散和蒸发三个过程.研究了粒子间相互作用范围α和允许粒子行走的最大步数对薄膜生长形貌的影响.结果表明:在不同的α值下,随粒子行走步数的增加,薄膜的生长经历了从分散、分形、混合到团聚的过程;其中α=6时,基本观察不到粒子的分散生长过程;α值越小且粒子行走步数越小的情况下,薄膜越易趋向于分散生长 关键词： 超薄膜 MonteCarlo方法 外延生长 Morse势 分形     

非均质基底表面上团簇生长的Monte Carlo模拟

Monte Carlo方法模拟了在非均质基底表面上团簇的生长过程.非均质基底由规则分布的A，B两类不同性质的物质(不同材料、不同温度等)组成.沉积粒子的初动能为E₀，在不同相基底上单位长度扩散的能耗为E_A,E_B，粒子从A经过相界扩散到B的能耗为E_AB.模拟中粒子的初动能E₀.远大于E_A,E_B，粒子从B过相界到A能耗E_BA取为0.结果表明，沉积在此类基底表面的粒子凝聚成具有分形结构的团簇.当E_AB远大于E_A,E_B时，团簇的数目、团簇的回旋半径、团簇的分形维数均随E₀/E_AB呈周期性变化. 关键词： 分形团簇 非均质基底 Monte Carlo模拟 周期性     

磁性纳米颗粒系统偶极相互作用的Monte Carlo研究

采用局域Monte Carlo方法模拟不同易轴分布的简单立方排列单分散单畴Fe纳米颗粒系统的ZFC-FC曲线及磁滞回线.结果表明：随着偶极相互作用的增强，系统的阻塞温度T_B逐渐增大，且ZFC曲线的峰变宽.说明偶极相互作用使得系统的有效能垒提高，分布宽度增加.研究FC曲线磁化强度的倒数与温度关系，发现偶极相互作用系统中存在反铁磁有序.系统的阻塞态及超顺磁态的磁滞回线表明，极低低温下，随着偶极相互作用的增强，系统的矫顽力和剩磁减小，偶极相互作用阻碍系统的磁化；系统处于超顺磁态，各向异性作用及偶极相互作用使得系统的磁化曲线偏离Langevin曲线且偶极相互作用展现出退磁相互作用效应.偶极相互作用增强，系统磁化曲线与Langevin曲线偏差量的最大值向低场移动.在偶极相互作用下，易轴与外场夹角为45°的磁性纳米颗粒系统的平均有效能垒和有效能垒分布宽度较易轴随机分布系统的大.     

旋转磁场作用下磁流变液颗粒运动及结构演化的模拟

通过对旋转磁场作用下磁流变液一定数目颗粒运动的数值模拟,在旋转平面内得到盘状聚集结构,在垂直于旋转平面的平面内得到层状结构,得出旋转磁场作用下磁流变液的结构特点。在模拟过程中颗粒受到磁场作用力简化为磁偶极子间相互作用。     

中子引起单粒子翻转过程的Monte Carlo计算模拟

从中子与硅原子相互作用的物理机理出发,利用Monte Carlo方法编制了中子引起单粒子翻转的计算模拟程序,并对14 MeV中子环境下的16K位静态存储器硅片翻转过程中的物理量进行了计算,同时可为中子引起的单粒子翻转的研究提供截面和描述内部物理过程的参考数据。     

用分子自组装技术制备的单电子器件的Monte Carlo模拟

用分子自组装技术制备出纳米金单电子器件,并测量了其伏安特性,根据单电子系统的半经典理论,用MonteCarlo法对其结果进行了模拟.结果表明,模拟出的伏安曲线与实测的伏安曲线有较好的一致性,反映了模拟方法用于单电子器件研究的合理性,此外发现,虽然单电子器件两电极间含有众多的纳米粒子,但在低压区,其伏安特性只与少数纳米粒子有关 关键词： 单电子器件 MonteCarlo模拟 分子自组装     

Ih对称性点阵上Ising相互作用自旋系统的族团Monte Carlo研究

应用一种有效的族团迭方法，我们在具有Ｉｈ对称性点阵上作了关于Ｉｓｉｎｇ自旋相互作用系统的ＭｏｎｔｅＣａｒｌｏ模拟，热力学量的模拟讲师结果显示，在β＝０．６４温度点上系统的磁化率和比热存在一个明显的峰，表征出系统的磁有序性开始消失，该结果与高温组数展开的解析计算的结果一致。     

平面旋转场诱导下磁流变液的结构模型

针对平面旋转场下磁流变液中诱导产生的层状结构 ,提出了恰当的理论模型 ,该模型成功的解释了层状结构在同步旋转状态下 :层间距随磁场的变化关系 ;旋转层的滞后角随旋转角频率以及其他物理量之间的关系 ;层半径与旋转角频率之间的关系 .理论结果与文献中发表的实验结果相吻合     

铜团簇Cun(n=2-60)的基态能量与几何结构的Monte Carlo模拟

采用蒙特卡洛(MC)方法和嵌入原子势(EAM)函数,研究了铜原子团簇的结构及其结合能.表明铜原子团簇在n值较小时,趋向于立体结构,而非平面结构;n＞13的铜团簇都是在正二十面体的基础上添加原子形成的,通过分析团簇的结构和能量之间的关系发现,如果添加原子后形成的Cun具有较高的对称性,那么这种铜团簇的结构就稳定.     

磁性薄膜自旋重取向行为的Monte Carlo模拟

利用Monte Carlo方法模拟了二维简单立方结构磁性薄膜的自旋重取向行为,重点研究了各向异性和偶极相互作用对系统自旋取向的影响．通过计算,获得了系统的相图以及系统组态、磁分量、比热等随偶极相互作用和温度的变化规律．模拟结果表明,在一定的参数范围内,随着温度的升高,系统的自旋取向将由垂直向平行方向转变。     

磁性单层膜磁学性质的Monte Carlo模拟

采用类Ising模型,利用Monte Carlo方法研究了磁性单层膜中退磁性偶极作用和铁磁性交换作用对系统磁学性质的影响.结果显示,随着偶极相互作用的增加,系统在低温下的磁化出现平台现象,此时磁化曲线可分为2个阶段,在低外场下,温度升高,系统易磁化,在高外场下则反之.这种新奇的磁化行为导致系统的磁熵变在低温低外场下出现大于零的反常行为.在模拟过程中,对长程力作用采用了比较精确的处理方法.     

铜团簇(n=55)结构及能量随温度演变的Monte Carlo 模拟研究

采用Monte Carlo方法及Embeded-Atom Method作用势，以55个铜原子体系为对象，研究了高 温铜团簇的结构及能量.首先计算温度为200?K时Cu₅₅的能量及结构，其结果 与土耳 其Süleyman ?zcelik 等人的研究小组的研究结果符合很好，说明此方法是可靠的.然后， 从12000?K开始逐渐降温，间隔1000?K，每个温度点MCS为10⁶，在达到平衡时 得到体系的结构及能量.结果表明：在10000?K以上时，体系只有少数原子结合在一起成核，成核数很少，且很不稳定，单原子数较多；在10000?K到6000?K温度段，随着温度的降低，体系 结构由单原子双原子数较多渐变到有较多原子结合在一起，形成一个大的成核中心，剩余单 原子、双原子逐渐与成核中心集结；在温度5000?K达到平衡时，体系原子已经完全结合在 一起，没有单个原子存在，体系能量为-11206?eV，在进一步降温的过程中，体系结构由 5000?K时的不对称、不规则渐变到600?K时的规则、对称的双二十面体结构，能量为-168 50?eV；在可视化的显示下，得到一幅十分清晰的关于体系随温度降低的结构演变图. 关键词： 铜团簇 结构及能量 Monte Carlo 模拟 高温     

聚合物中多重光散射传导的Monte Carlo数值模拟

基于米氏散射(Mie scattering)理论,建立填充分散粒子群的聚合物对光散射传导的Monte Carlo数学模型.在此基础上,编写了一套仿真模拟程序.通过模拟单个光子在聚合物中的多重散射运动过程,把问题扩展到以激光束或线状光为入射源,得到在聚合物板块内的光传导情况,并且在计算机上图像化地重现整个物理过程,对输出光强的分布情况进行模拟统计分析.模拟结果表明,利用体散射机制,可以将点光源和线光源转换为平面光输出,输出光的状态可以通过对比计算结果实施有效控制.     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

团簇淀积纳米结构薄膜的计算机模拟

对团簇的淀积方式作了简介。着重对团簇的淀积过程进行了计算机模拟,从中得出了淀积的尺寸分布与团簇在衬底上的徙动长度有关,而团簇的聚合长大存在一临界尺寸。并给出低能Ｐｂ团簇在碳膜表面的徙动长度和聚合的临界尺寸。     

酞菁钴/铁纳米填充母粒组成的磁流变液性能

采用有机／无机原位（ｉｎ ｓｉｔｕ）复合方法制备得出本菁钴／铁纳米填充母粒，与甲基硅油组成磁流变液（ＭＲＳ）。ＭＲＳ的附加动态剪切应力（△τ）与分散介质浓度、外加磁场强度呈正比例关系；剪切速率对△τ的影响表明磁致流变为链状结构特征；△τ的对温度不敏感；ＭＲＳ对外加磁场有可逆的开／关变化特征，无记忆效应，磁流变响应时间小于０．１ｓ。     

模拟颗粒凝并过程的快速Monte Carlo方法

对常规异权值Monte Carlo(MC)方法进行改进,基于强核函数思想,通过颗粒群单重遍历即可求得强核函数最大值,采用接受-拒绝法随机搜寻凝并对,并利用搜寻过程中拒绝和接受的所有凝并对的信息来估计凝并事件的等待时间(时间步长),从而避免颗粒群的双重遍历,以提高MC的效率.对典型工况的模拟结果显示该快速方法计算代价仅为O(N_s),能够显著提高计算效率,同时保持足够的计算精度,较好地协调计算代价与计算精度之间的矛盾.     

铜原子纳米团簇热力学性质的分子动力学模拟研究

本文利用分子动力学模拟方法,研究了CuN(N=80、140、216、312、408、500、628和736)纳米团簇在热化和冷凝过程中结构和热力学性质的变化,模型采用的是Johnson的EAM作用势.模拟结果表明:铜团簇的熔点和凝固点随其尺寸线性增加,并逐渐向大块晶体靠拢;所有团簇的凝固点都低于熔点,出现凝固过程中的滞后现象;在熔点和凝固点附近,团簇都具有负热容特性,负热容是由相变前后团簇内部结构突变引起的.     

铜原子纳米团簇热力学性质的分子动力学模拟研究

利用分子动力学模拟方法,研究了CuN(N=80,140,216,312,408,500,628和736)纳米团簇在热化和冷凝过程中结构和热力学性质的变化,模型采用的是Johnson的EAM作用势.模拟结果表明:铜团簇的熔点和凝固点随其尺寸线性增加,并逐渐向大块晶体靠拢;所有团簇的凝固点都低于熔点,出现凝固过程中的滞后现象;在熔点和凝固点附近,团簇都具有负热容特性,负热容是由相变前后团簇内部结构突变引起的.