1,1-二氨基-2,2-二硝基乙烯晶体热膨胀各向异性的模拟研究

对新型钝感高能炸药组分1,1-二氨基-2,2-二硝基乙烯(FOX-7),考察并选择了适用的分子力场和超晶胞初始构型,通过303K-378K温度范围在NPT系综下的分子动力学模拟(MD),得到常压下各温度的FOX-7晶体平衡构型,发现随温度升高晶体内分子的取向和堆积方式不变;通过线性拟合求算出线膨胀系数,发现αb明显较大而αc与αa接近但不相等;接着采用密度泛函理论(DFT)方法随各向棱长膨胀的能量变化进行了模拟计算,发现得到的能量变化成各向异性并与各向热膨胀系数相关,从而建立FOX-7晶体特定分子堆积模式与热膨胀各向异性的关联.     

1,1-二氨基二硝基乙烯晶体的密度泛函理论研究

对 1,1 二氨基二硝基乙烯晶体进行了DFT B3LYP水平计算 .计算所得晶格能为 -10 5 .81kJ/mol,与文献值相近 .晶体的前线能带较为平坦 ,表明分子轨道能态受分子晶体场的影响较小 .电荷的分布决定了晶体中DADNE以“头 尾”方式通过分子间氢键相连形成层状结构 ,而层与层之间相互作用较弱 .从带隙 4.0eV推知DADNE的导电性介于半导体和绝缘体之间 .前线轨道由C -NO2 的原子轨道所组成 ,说明分子的强共轭性 ,也表明C -NO2 为化学反应活性部位 .C -NO2 键的布居数远小于其它键 ,提示该键为起爆引发键     

β-1,3,5,7-四硝基-1,3,5,7-四氮杂环辛烷热膨胀各向异性的量子化学计算与分子动力学模拟

在303～383 K和NPT系综和COMPASS力场下对β-1,3,5,7-四硝基-1,3,5,7-四氮杂环辛烷(HMX)超晶胞初始结构的分子动力学模拟,得到常压下各温度的晶体平衡构型并发现分子的堆积方式不变;通过线性拟合求算出线膨胀系数与实验值相近,体现出明显的各向异性. 采用密度泛函理论方法对沿各晶轴方向膨胀率变化(100%～105%)的HMX单胞模型进行了总能计算,得到的能量变化率体现各向异性并与热膨胀系数值关联,建立了关联方程. 由此阐 释了HMX晶体热膨胀各向异性的本质即特定的分子堆积模式.     

4-氨基-3,5-二硝基吡唑二聚体分子间相互作用的理论研究

在密度泛函（DFT）B3LYP/6-311++G**水平下,求得4-氨基-3,5-二硝基吡唑（LLM-116）二聚体势能面上的4种全优化几何构型。经基组叠加误差（BSSE）和零点能（ZPE）校正, 求得分子间最大相互作用能为 –29.56 kJ×mol-1。研究表明,其相互作用是放热过程,并随着温度的升高,相互作用减弱;自然键轨道（NBO）分析揭示了相互作用的本质;对优化构型进行振动分析,得其红外光谱,并基于统计热力学求得200.0-800.0 K 温度范围从单分子形成二聚体的热力学性质变化。     

双(2,2-二硝基丙基)甲缩醛的结构和爆轰性能的量子化学研究

利用B3LYP/6-311+G(2d,p)方法对一种新型含能增塑剂双(2,2-二硝基丙基)甲缩醛进行几何优化,计算了其红外光谱、生成焓和爆轰特性. 分析了最弱键的键离解能和键级并预测了目标化合物的热稳定性. 结果表明双(2,2-二硝基丙基)甲缩醛中的四个N-NO₂键的键离解能都为164.38 kJ/mol. 表明目标化合物是一个热力学性能稳定的化合物. 以凝聚相生成焓和分子密度为基础,采用Kamlet-Jacobs方法预测其爆速和爆压. 目标化合物的晶体结构属于P21空间群.     

1，1-二氨基-2，2-二硝基乙烯的晶型及相变

人们一直尝试合成可以替代现今高能炸药的新型炸药，期望新型炸药在相同的感度水平有更高的性能，或者在相同的性能水平下有更低的感度。1，1-二氨基-2，2-二硝基乙烯（简称FOX-7）就是在这种要求下研究较为活跃的一种新型低感炸药。文献报道FOX-7有3种或4种相态，在室温下存在α相，在113℃转变成β相，下一个相变开始于160℃或170℃。文中利用四圆衍射仪测定FOX-7的晶体结构，变温X射线衍射仪（XRD）研究FOX-7在加热时发生的相变。     

紫外吸收光谱研究4,4′-二硝基二苯乙烯-2,2′-二磺酸稀溶液的光致顺反异构

首次发现并证实了在4,4′-二硝基二苯乙烯-2,2′-二磺酸稀溶液中存在顺反异构现象。采用紫外吸收光谱对其稀溶液的光致顺反异构进行了研究。结果表明，溶液的浓度越稀，光照时间越长，反式异构体转化成顺式异构体越快、越多。     

2-氨基吡啶气相氟化取代反应的密度泛函理论研究

采用量子化学方法,运用密度泛函理论(DFT)B3LYP方法,在6-311G*基组水平上,以2-氨基吡啶为起始反应物,全优化计算了亚硝酰阳离子NO+对2-氨基吡啶中氨基进攻的重氮化反应及氟硼酸根阴离子BF4-对吡啶重氮正离子进攻的席曼反应的机理,搜索到整个反应过程中的过渡态结构,得到了反应路径并通过IRC验证,并对各反应中间体、过渡态和产物进行分子几何构型优化、振动频率分析和标准热力学函数计算．通过计算,搜索到的过渡态均为环状结构,与配位插入反应类似,利于原子间化学键的断裂及新键的形成;通过热力学及动力学计算,重氮化-席曼反应是剧烈的放热反应,与实验结果相吻合．     

2-氨基吡啶气相氟代反应的密度泛函理论研究

采用量子化学方法,运用密度泛函理论(DFT)B3LYP方法,在6-311G*基组水平上,以2-氨基吡啶为起始反应物,全优化计算了亚硝酰阳离子NO~+对2-氨基吡啶中氨基进攻的重氮化反应及氟硼酸根阴离子BF_4~-对吡啶重氮正离子进攻的席曼反应的机理,搜索到整个反应过程中的过渡态结构,得到了反应路径并通过IRC验证,并对各反应中间体、过渡态和产物进行分子几何构型优化、振动频率分析和标准热力学函数计算.通过计算,搜索到的过渡态均为环状结构,与配位插入反应类似,利于原子间化学键的断裂及新键的形成;通过热力学及动力学计算,重氮化-席曼反应是剧烈的放热反应,与实验结果相吻合.     

二芳基乙烯光致各向异性动态特性研究

建立了光致变色二芳基乙烯样品光致各向异性动力学的理论模型.计算了在线偏振光激发下,二芳基乙烯分子数密度、光致二向色性和光致双折射随时间(曝光量)的动态变化曲线.结果表明,随着激发曝光量的增加,二芳基乙烯呈色态分子数密度不断下降直至饱和,光致各向异性随曝光量增加先增加,达到最大之后开始下降,最佳曝光量为260J/cm².利用正交偏振检测方法测量了反映材料各向异性的探测光透过率动力学曲线,并与理论计算结果进行了比较,二者基本相符. 关键词： 光致各向异性 二芳基乙烯 动态特性 最佳曝光量     

Ab initio study of the structure and dynamical properties of crystalline ice

We investigated the structural and dynamical properties of a tetrahedrally coordinated crystalline ice from first principles based on density functional theory within the generalized gradient approximation with the projected augmented wave method. First, we report the structural behaviour of ice at finite temperatures based on the analysis of radial distribution functions obtained by molecular dynamics simulations. The results show how the ordering of the hydrogen bonding breaks down in the tetrahedral network of ice with entropy increase, in agreement with the neutron diffraction data. We also calculated the phonon spectra of ice in a 3× 1× 1 supercell using the direct method. So far, due to the direct method used in this calculation, the phonon spectra are obtained without taking into account the effect of polarization arising from dipole–dipole interactions of water molecules, which is expected to yield the splitting of longitudinal and transverse optic modes at the Γ point. The calculated longitudinal acoustic velocities from the initial slopes of the acoustic mode are in reasonable agreement with the neutron scattering data. Analysis of the vibrational density of states shows the existence of a boson peak at low energy of the translational region, a characteristic common to amorphous systems.     

采用密度泛函理论与分子动力学对聚甲基丙烯酸甲酯双折射性的理论计算

双折射性是各种光学材料的重要性能之一,具有高双折射率的光学材料在诸多研究及工业领域的应用越来越广泛.然而,作为常用的光学薄膜及光波导材料之一的聚合物材料的双折射性通常却很弱,只能通过实验对其双折射率进行大致的表征,缺乏对其双折射率的系统性理论计算,从而限制了提高聚合物双折射性的研究.本文建立了从聚合物的单体分子结构到多分子链的系统性的双折射率理论计算方法,并借助此方法研究了导致聚合物弱双折射性的限制因素.以聚甲基丙烯酸甲酯(PMMA)为研究对象,运用密度泛函理论研究了其本征双折射率,这里的本征双折射率是指分子链完全取向时其单体单元的双折射率.计算结果表明其本征双折射率高达0.0738左右,并且通过计算给出了PMMA单体单元的平均双折射率色散曲线.采用分子动力学方法研究了该聚合物(包含20个分子链)的材料双折射率.理论计算结果表明,尽管该聚合物具有较大的本征双折射率,但是由于聚合物中分子链取向度极低,聚合物材料最终表现出来的双折射率只有0.00052.本文建立的研究方法及研究结果为研究增强聚合物材料双折射性提供了理论依据.     

Vibrational analysis of L-serine using the density functional theory

In this paper, we present a computational study of L-serine using ab initio molecular dynamics simulation based on density functional theory (DFT) within the ultrasoft pseudopotentials and generalized-gradient approximation. Taking into account the intermolecular interactions, we can indeed simulate the features of the experimental results very well for L-serine zwitterions in its solid state. The vibrational spectrum of L-serine performed by DFT was in excellent agreement with our previous inelastic incoherent neutron scattering spectra measured at 20K for L-serine in the 10--200meV region on HET spectrometers at ISIS, Rutherford Appleton Laboratory.     

质子与羟基碰撞的含时密度泛函理论研究

采用将含时密度泛函理论和分子动力学非绝热耦合的方法,研究了不同入射速度下质子与羟基碰撞的反应动力学.计算了碰撞前后质子动能和羟基动能的变化及羟基电子和质子的运动.计算结果表明,质子沿垂直羟基分子轴方向入射时,质子与羟基碰撞后,质子被反弹且动能损失并俘获了羟基中氧的一部分电子,而丢失部分电子的羟基则获得动能以伸缩振动的形式向计算边界平动.随着入射质子的初动能增加,质子从羟基中俘获的电子增多,碰撞后羟基的键长变长,羟基振动变强而伸缩振动频率降低.此外,还发现质子的入射方向对碰撞过程的激发动力学有很大的影响.质子从不同的方向入射时,质子的入射初动能越大,其损失的动能越多且损失的动能与入射初动能呈线性关系,而入射方向对质子动能损失的影响很小.在质子入射初动能较低(小于25 eV)的情况下,羟基获得的动能与质子入射初动能呈线性关系且与入射方向无关;在质子入射初动能较高(大于25 eV)时,当质子沿羟基分子轴方向入射时,羟基动能的增量远大于质子沿垂直于羟基分子轴方向入射时羟基动能的增量.     

单层二硫化钼多相性质及相变的第一性原理研究

本文基于密度泛函的第一性原理,并引入范德瓦耳斯力修正,研究了单层二硫化钼2H,1T,ZT三种相的电学性质及相变原理.首先通过结构弛豫确定了三种相的几何结构,能带和态密度计算证实1T相具有金属性质,ZT相具有半导体性质,带隙为0.01 eV.然后结合变形势理论计算了2H和ZT相的迁移率,ZT相的迁移率高达104cm~2·V~(-1)·s~(-1),进一步拓展了单层二硫化钼的应用范围.最后通过对比三种相吸附锂原子结合能,计算2H-1T相变能量曲线,解释了引起二硫化钼相变的原因.本文的研究结果将对单层二硫化钼实验制备表征以及相关光电器件性能分析提供重要参考.     

First principles study of oxygen adsorption on nickel-doped graphite

Density functional theory is used in a spin-polarized plane wave pseudopotential implementation to investigate molecular oxygen adsorption and dissociation on graphite and nickel-doped graphite surfaces. Molecular oxygen physisorbs on graphite surface retaining its magnetic property. The calculated adsorption energy is consistent with the experimental value of ?0.1?eV. It is found that substituting a carbon atom of the graphite surface by a single doping nickel atom (2.8% content) makes the surface active for oxygen chemisorption. It is found that the molecular oxygen never adsorbs on doping nickel atom while it adsorbs and dissociates spontaneously into atomic oxygens on the carbon atoms which are bound to the nickel. The adsorption energy of ?1.4?eV and zero activation energy barrier indicate that O₂ dissociative adsorption is both thermodynamically and kinetically favoured over the surface. The large electric field near the doping nickel atom along with the excess electrons on the neighbouring carbon atoms, which are bound to the nickel induce molecular oxygen to adsorb and dissociate favourably.     

U_(1-x)Pu_xO_2热膨胀性质分子动力学模拟研究

基于部分离子势函数的分子动力学方法,研究了U_(1-x)Pu_xO_2(x=0.07,0.15,0.25和0.5)在300—3000 K温度范围和0—1.5 GPa下的热力学性质,研究发现,在等压和不同温度下,随着钚(Pu)比率的增加,U_(1-x)Pu_xO_2的晶格常数线性减小,线性膨胀系数有所增大,等温压缩系数减小,通过拟合得到了零压下U_(1-x)Pu_xO_2晶格常数和线膨胀系数随Pu比率变化的经验表达式。     

A theoretical study on the mechanism and dynamics of reactions (CF3)2CHOCH2F/(CF3)2CHOCHF2 with OH radical

The information related with the mechanism of reactions (CF₃)₂CHOCH₂F + OH (R1) and (CF₃)₂CHOCHF₂ + OH (R2) was explored theoretically at the BMC-CCSD//BMK/6-311 + G(d,p) level. Based on the optimised structures, energies, and other information, the rate constants were evaluated by the canonical variational transition-state theory with small curvature tunneling contributions in a temperature range of 220–2000 K. For each reaction, there are both hydrogen-abstraction and displacement channels. In addition, more than one hydrogen atom can be abstracted. The relationship between hydrogen abstraction and displacement, between different hydrogen-abstraction channels, and between reactions R1 and R2 are elucidated.     

First principles determination of structural,electronic and lattice dynamical properties of a model dipeptide molecular crystal

We present the results of ab initio calculations of the structural, electronic and lattice dynamical properties of the solid-state crystal of the glycyl-l-alanine dipeptide. Intramolecular bond lengths are found to be in good agreement with experimental values; lattice constants are in reasonable agreement, although it is found that discrepancies do exist. A hierarchy of hydrogen bond strengths is found, with those between (oppositely-charged) amine and carboxy functional groups being strongest. The crystal is found to be an indirect-bandgap material, with indirect bandgaps ≈4.95 eV, compared to a direct bandgap of 5.00 eV. Analysis of the electronic structure reveals that the electronic states in the near vicinity of the energy gap arise from carboxylate and amide oxygen atoms. The arrangement of both molecules and hydrogen bonds in the unit cell is found to manifest itself in increased bandwidth along specific reciprocal space directions, reflecting coupling brought about by hydrogen bonds. Determination of the zone-centre lattice dynamical behaviour permits the IR absorption spectrum to be explained. Intermolecular hydrogen bonds are found to couple intramolecular motions in adjacent moelcules, revealing the importance of an accurate treatment of intermolecular interactions, even for high-frequency vibronic modes.     

Density function study transition metal chromium-doped alkali clusters: the finding of magnetic superatom

The structures, stabilities and magnetic properties of CrX_n (X = Na, Rb and Cs; n up to 9) clusters are studied using density functional theory to search for the stable magnetic superatoms. The geometrical optimisations indicate the ground-state structures of CrX_n evolve toward a close packed structure with an interior Cr atom surrounded by X atoms as the cluster size increase. Their stabilities are analysed by the relative energy, gain in energy (ΔE(n)) and the highest unoccupied molecular orbital and lowest unoccupied molecular orbital gaps. Furthermore, the magnetic moments of CrX_n clusters show an odd–even oscillation. Here, we mainly focus on the CrX₇ (X = Na, Rb and Cs) clusters due to the same valence count as the known stable magnetic superatoms VNa₈, VCs₈ and TiNa₉. Although these clusters all have a filled electronic configuration 1S²1P⁶ and large magnetic moment 5 μ_B, our studies indicate that only CrNa₇ is highly stable compared to its nearest neighbours, while CrRb₇ and CrCs₇ clusters are less stable. This suggests that Cr-doped Na₇ is most appropriate for filled electronic configuration and CrNa₇ is shown to be a stable magnetic superatom. More interesting, we find CrRb₈ and CrCs₈ with the filled electronic configuration 1S²1P⁶ have higher stability and large magnetic moment 6 μ_B in their respective series.