BaTiO3铁电体中辐射位移效应的分子动力学模拟

运用基于壳模型的分子动力学方法研究了BaTiO₃铁电体中的辐射位移效应.采用O原子作为初级击出原子,模拟了当初级击出原子能量为1 keV时体系内缺陷的产生和演化.模拟结果表明,当入射方向为[001]时,体系内产生的缺陷最多.在所有缺陷中,以O缺陷的含量为最高,达80%以上.同时,这些缺陷的产生并不显著改变体系的自发极化强度,对体系的极化翻转过程也基本没有影响.在外电场作用下,观察到了显著的缺陷迁移. 关键词： 分子动力学 3铁电体')" href="#">BaTiO₃铁电体 辐射位移效应     

功能化单壁碳纳米管的热导率

采用基于BrennerⅡ势的非平衡态分子动力学方法,模拟研究了300K温度下经氢化学修饰的(10,0)单壁碳纳米管的热导率.研究显示功能化后碳管的热导率有明显减小,当有一列碳原子被氢化后（功能化程度为5%）,碳管的热导率减小了大约1/3,为了进一步解释这种功能化对碳纳米管热导率的影响,计算了不同功能化程度下碳纳米管的声子谱.     

分子动力学模拟纯镁熔体凝固过程中结构演化

采用分子动力学模拟方法模拟了在周期性边界条件下由500个原子构成的液态Mg模型系统的凝固过程,分别考察了在5×10^14 K/s、5×10^13 K/s、1×10^13 K/s 、1×10^12 K/s的冷却速率下液态Mg熔体的凝固行为。模拟结果很好地重现了实验值。模拟中原子间作用势采用FS势,结构分析采用径向分布函数、均方位移、系统总能量分析、H-A键对分析技术等方法。结果表明,当冷却速率为5×10^14 K/s时,系统形成以1541键对为主的非晶态结构;当冷却速率分别为5×10^13 K/s、1×10^13 K/s、1×10^12 K/s时,系统形成以1421、1422键对为主的hcp晶态结构;另外,在快速冷却形成非晶的过程中,大部分bcc结构被保留下来,而在慢冷形成晶态的过程中,大部分bcc结构最终演化形成了hcp结构。     

电子关联对聚乙炔中激子类型的影响

在SSH哈密顿基础上引进电子关联,对反式聚乙炔链中光致激子的产生和演化过程实施分子动力学模拟。弱关联强度U取值0～1.250 eV,V =U/2取值0～0.625 eV .计算结果表明,关联强度的大小影响链中元激发类型,U＜0.555eV时产生的元激发为孤子-反孤子对,U ＞0.555eV 时产生的元激发为正负荷电极化子对。为进一步讨论该类型一维系统中不同类型激子的产生、输运、衰减等动态过程,关联强度U的选择提供参考。     

胶体铁磁流体的光子特性研究

这篇综述回顾了我们今年来在铁磁流体方面的一些研究工作。铁磁流体是一种纳米磁颗粒或亚铁磁颗粒悬浮在基液（如：水或煤油）中形成的悬浮液。铁磁流体之所以受到广泛关注,是因为它们在许多方面都有着潜在应用,例如在机械工程和生物医学等领域。再次综述中,我们首先介绍了铁磁流体的场感应各向异性结构,然后,介绍了几种铁磁流体基软物质材料的光学性质,即：光学负折射、磁控光子待隙和非线性光学响应。我们采用的研究方法主要是分子动力学模拟,有效煤质近似和有限元模拟.     

固体微/纳米尺度传热理论研究进展

随着半导体技术的飞速发展,器件的尺寸已进入到微/纳米尺度。由于量子效应、表面及界面效应,使得微尺度下的热物性与宏观尺度下有了明显的区别。人们针对微观传热领域的特点,发展了声子玻尔兹曼传输方程、分子动力学等方法,取得了一定的成果,但仍存在不少问题。一些基础概念问题,特别是非平衡态下的局域温度的定义,有待澄清。本文主要回顾近年来微/纳米尺度传热在理论和数值模拟方面的进展,以及目前所面临的挑战和问题。     

高剪切率的Couette流下液氩流体热导率的分子动力学模拟

基于平衡分子动力学(EMD)模拟方法,对Lees-Edwards周期性边界条件作一修正并成功地应用其形成具有高恒定剪切率的Couette流动,最后通过Green-Kubo公式计算得到了剪切流下液氩流体的热导率。计算结果,表明在剪切流下液氩流体的热导率仍然表现为各向同性,同时还发现流动方向的热导率不随流动剪切率的增加而发生变化,并且该值跟静态下液氩的热导率(λ=0.132 Wm~(-1)K~(-1))基本一致。     

单壁碳纳米管内流体流动的手性效应

纳米尺度范围内流道分子结构形成的粗糙度将会影响其中的流体流动。本文采用分子动力学方法,以氩为工质,模拟了不同手性的单壁碳纳米管（SWCNT）内流体的流动。结果表明：由于范德瓦耳斯力的作用,流体分子与壁而间有一定的距离。可以看出在壁面附近有两个明显的分层,这是纳米流动中持有的流体密度分层现象。由于不同手性碳纳米管的分子排...     

A new probe of neutron skin thickness

The correlation between neutron-to-proton yield ratio (R_np) and neutron skin thickness (δ_np) in neutron-rich projectile induced reactions is investigated within the framework of the Isospin-Dependent Quantum Molecular Dynamics (IQMD) model. The density distribution of the Droplet model is embedded in the initialization of the neutron and proton densities in the present IQMD model. By adjusting the diffuseness parameter of neutron density in the Droplet model for the projectile, the relationship between the neutron skin thickness and the corresponding R_np is obtained. The results show strong linear correlation between R_np and δ_np for neutron-rich Ca and Ni isotopes. It is suggested that R_np may be used as an experimental observable to extract δ_np for neutron-rich nuclei, which is very interesting in the study of the nuclear structure of exotic nuclei, the equation of state (EOS) of asymmetric nuclear matter and neutron-rich matter in astrophysics, etc.     

Zn-Mg合金的长程Finnis-Sinclair势

通过拟合Mg的晶格能、晶格常数、弹性常数,并将其与前人的结果相比较后获得了描述Mg的最优长程Finnis-Sinclair(F-S)势函数参数,使用同样方法并引入修正因子后得到了Zn的长程F-S势参数.基于单质Zn,Mg的F-S势参数,进一步拟合合金Mg₂₁Zn₂₅,MgZn₂,Mg₂Zn₁₁的晶格常数、晶格能获得Zn-Mg原子对的F-S势参数,构建了整套描述Zn-Mg合金的长程F-S势参数.在此 关键词： 长程Finnis-Sinclair势 Zn-Mg合金 分子动力学模拟