Application of shell model in molecular dynamics simulation to MgO

The P－V－T equation of state of MgO has been simulated under high pressure and elevated temperature using the molecular dynamics (MD) method with the breathing shell model (BSM). It is found that the MD simulation with BSM is very successful in reproducing accurately the measured molar volumes of MgO over a wide range of temperature and pressure. In addition, the MD simulation reproduces accurately the measured volume compression data of MgO up to 100GPa at 300K. It is demonstrated that the MD simulated P－V－T equation of state of MgO could be applied as a useful internal pressure calibration standard at elevated temperatures and high pressures.     

两体和三体相互作用下的非线性Landau-Zener隧穿的圈结构研究

在加速光晶格中玻色-爱因斯坦凝聚体（BEC）,当同时考虑两体和三体相互作用时,其能级结构、隧穿率出现了独特的特性.在一维加速光晶格中的BEC,当两体和三体相互作用参数满足一定条件时,非线性两能级体系的能级结构中出现了圈结构,研究得出了圈的宽度随两体和三体相互作用参数变化的关系,并由此分析圈结构的出现及其大小对BEC隧穿率的影响.     

球形复合柱表面波声子晶体的带隙特性仿真

设计了一种由镍球与环氧树脂垫层组成的复合柱沉积在铌酸锂基体上构成的表面波声子晶体结构,采用有限元法计算了其能带结构和位移矢量场.结果表明:与具有相同晶格常数的倒圆锥形表面波声子晶体结构相比,研究结构可以在更低的频率范围打开更宽的声表面波完全带隙,且随着复合柱半径增大,镍球体与压电基体的硬边界之间形成限制腔模,相邻高阶带隙间存在能量的耦合以及振动模式的继承;此外,温度场的引入可以实现带隙的主动调控,带隙频率范围随着温度升高向低频移动;通过增加复合柱体的层数,多振子结构与行波发生多极共振耦合,可在高阶能带间打开完全带隙.本文的研究结果为微米级表面波声子晶体结构在100 MHz以下频率范围的带隙特性优化提供了理论支持.     

含能材料中键离解能的密度泛函理论计算

通过计算一个包括臭氧,硝基甲烷以及1,3,5 -三硝基-1,3,5-三氮杂环已烷(RDX)在内的典型系统的键离解能,对由四种交换/相关函数(BLYP,B3LYP,B3PW91和B3P86)加上不同的基函数组合而成的多种密度泛函方法的准确性进行了比较研究.结果表明:B3P86/6-31G**是计算该系统C-NO2,O-O和N-NO2键离解能的最可靠的方法.     

高温高压下闪锌矿相GaN结构和热力学特性的分子动力学研究

利用分子动力学方法和Buckingham经验势模型对重要半导体材料GaN立方闪锌矿相的晶格常数、相变压力(从闪锌矿到岩盐结构)、热膨胀、等温体模量、定压热容等结构和热力学特性在300—3000K的温度范围和0—65GPa的压力范围内进行了研究.研究表明，闪锌矿相GaN常态下的结构和热力学参数的模拟结果与实验数据及其他理论结果相符.同时在所选作用势模型可靠性检验的基础上，对等温体模量、定压热容诸非谐性参量在高温高压下的热力学行为进行了预测.所得结果在材料科学等领域的研究中具有一定的应用背景和参考价值. 关键词： GaN Buckingham势 分子动力学模拟 高温高压     

硝酸丙酯键离解能和热解机理的密度泛函理论研究

 采用B3LYP方法、在6-31G*基组水平上,优化了硝酸丙酯单体的平衡几何构型,计算了分子中各键的离解能。结果表明该分子中最弱的键为O—NO₂键,次弱键为C—O键。同时研究了热解机理,探索性的预测了断裂O—NO₂键及C—O键后的产物,并采用从头算法、半经验方法和密度泛函理论分别计算了硝酸丙酯的生成热,由半经验方法中的PM3得到的数值和实验结果符合较好。     

Molecular dynamics study for the melting curve of MgO at high pressure

Shell-model molecular dynamics method is used to study the melting temperatures of MgO at elevated temperatures and high pressures using interaction potentials. Equations of state for MgO simulated by molecular dynamics are in good agreement with available experimental data. The pressure dependence of the melting curve of MgO has been calculated. The surface melting and superheating are considered in the correction of experimental data and the calculated values, respectively. The results of corrections are compared with those of previous work. The corrected melting temperature of MgO is consistent with corrected experimental measurements. The melting temperature of MgO up to 140GPa is calculated.     

Isothermal bulk modulus and its first pressure derivative of NaCl at high pressure and high temperature

The isothermal bulk modulus and its first pressure derivative of NaCl are investigated using the classical molecular dynamics method and the quasi-harmonic Debye model. To ensure faithful molecular dynamics simulations, two types of potentials, the shell-model (SM) potential and the two-body rigid-ion Born-Mayer-Huggins-Fumi-Tosi (BMHFT) potential, are fully tested. Compared with the SM potential based simulation, the molecular dynamics simulation with the BMHFT potential is very successful in reproducing accurately the measured bulk modulus of NaCl. Particular attention is paid to the prediction of the isothermal bulk modulus and its first pressure derivative using the reliable potential and to the comparison of the SM and the BMHFT potentials based molecular dynamics simulations with the quasi-harmonic Debye model. The properties of NaCl in the pressure range of 0-30 GPa at temperatures up to the melting temperature of 1050 K are investigated.     

Theoretical study on 2s2p6np Rydberg states of Cu19+ ion

We report on the calculations of transition wavelengths and weighted oscillator strengths for 2s²2p⁶-2s2p⁶ⁿp (4 ≤ n ≤ 20) electric dipole (E1) transitions of Cu¹⁹⁺ ion. The flexible atomic code (FAC) has been adopted for the calculations. Comparisons are made with the experimental data available, showing that the present results for 4 ≤ n ≤ 6 are more accurate than the previous calculated values. Furthermore, combining the quantum defect theory (QDT) with the transition energies of 2s²2p⁶-2s2p⁶np, the quantum defects for 2s2p⁶np Rydberg series of Cu¹⁹⁺ ion are determined. In addition, the energies of any highly excited states (n > 20) for this series can be reliably predicted using the QDT and the given quantum defects. The ionization energies for Cu¹⁹⁺ and Cu²⁰⁺ ions are also calculated and they excellently accord with previous experimental and calculated values.     

高压下钙钛矿结构MgSiO3的分子动力学研究

利用分子动力学方法,研究了高温高压下钙钛矿结构MgSiO3的状态方程.研究表明,分子动力学模拟结果很好地再现了广泛温度和压强范围内钙钛矿结构MgSiO3的摩尔体积.温度300 K压强上升到120 GPa模拟的钙钛矿结构MgSiO3状态方程和有效的实验结果基本一致.在更高温度和更高压强下模拟的钙钛矿结构MgSiO3状态方程和他人的计算值吻合的很好.另外,还分别计算了温度300 K,900 K,1500 K和2500 K压强上升到120 GPa时MgSiO3的体积压缩率.