HTPB与A1不同晶面结合能和力学性能的分子动力学模拟

采用分子力学（MM）和分子动力学（MD）方法,在250、300、350、400、450K,对固体推进剂端羟基聚丁二烯（HTPB）和铝晶胞不同晶面结构所组成的层模型在COMPASS力场下,进行模拟计算,求得结合能和静态力学性能（弹性系数、模量和泊松比）．模拟结果表明,在400K时HTPB与A1（011）面的结合能最大,从综合力学性能优劣上看,各个面从优到劣的排序为（011）〉（221）〉（001）,HTPB与A1的结合能与力学性能具有对应关系,结合能大的力学性能优异,结合能小的力学性能较差．     

HTPB/增塑剂玻璃化转变温度及力学性能的分子动力学模拟

为了预测高分子粘结剂端羟基聚丁二烯(HTPB)与增塑剂癸二酸二辛酯(DOS)、硝化甘油(NG)的相容性及HTPB/增塑剂共混物的玻璃化转变温度(Tg)和力学性能,在COMPASS力场条件下采用分子动力学(MD)模拟方法对相容体系(HTPB-DOS)和不相容体系(HTPB-NG)进行了研究.结果表明,通过比较溶度参数差值(Δδ)的大小可以预测HTPB与增塑剂的相容性,即HTPB与DOS属于相容体系,而HTPB与NG不相容.通过温度-比容曲线可以得到HTPB、HTPB/DOS与HTPB/NG的Tg分别为197.54,176.30和200.03K.力学性能分析结果表明,添加DOS增塑剂后使HTPB的弹性模量(E),体积模量(K)和剪切模量(G)下降,材料刚性减弱,柔性增强,力学性能得到改善.本模拟方法可以作为预测聚合物/增塑剂共混物性能的有利工具,也可以为固体推进剂和高聚物粘结炸药的配方设计提供理论指导.     

CL-20/HMX共晶及其为基PBX界面作用和力学性能的MD模拟研究

为提高共晶炸药的实际使用价值, 改善其安全性和力学性能, 以CL-20/HMX共晶炸药为基, 分别添加2种高聚物粘结剂Estane 5703(聚氨基甲酸乙酯)和HTPB(端羟基聚丁二烯), 共构建两种共晶基高聚物粘结炸药(PBX)模型, 进行细致的295 K NPT分子动力学(MD)模拟研究. 通过两种PBX模型及其与该共晶炸药的MD模拟结果比较表明, 与基炸药之间的结合能Estane 5703大于HTPB, 预示含少量Estane 5703的PBX稳定性和相容性更佳; 对相关函数g(r)揭示粘结剂与基炸药界面相互作用的方式, 以基炸药中H分别与Estane 5703中羰基O和HTPB中端羟基O之间的氢键较强. 与CL-20/HMX共晶炸药相比, 少量粘结剂Estane 5703或HTPB的加入, 使弹性系数C_ij下降, 拉伸模量(E)、体积模量(K)和剪切模量(G)均显著减小, 而泊松比(ν), 柯西压(C₁₂－C₄₄)和K/G值明显增大, 表明PBXs体系刚性减小, 延展性增强, 力学性能大为改善. 少量粘结剂包覆使PBXs致钝, 主要归因于其隔热、吸热并使体系变“软”的缓冲作用, 而界面作用造成的分子结构引发键键长变化变为次要因素.     

端羟基聚丁二烯/增塑剂共混物相容性的分子动力学模拟和介观模拟

应用分子动力学(MD)和介观动力学(MesoDyn)模拟方法对固体推进剂中端羟基聚丁二烯(HTPB)与增塑剂癸二酸二辛酯(DOS)、硝化甘油(NG)的相容性进行了研究. 采用MD模拟方法在COMPASS力场下, 对纯物质、HTPB/增塑剂共混物的密度、内聚能密度、溶度参数和共混物分子间的Flory-Huggins作用参数及结合能等进行了模拟计算, 通过比较溶度参数差值(Δδ)的大小、模拟前后体系密度变化情况均可以预测HTPB与增塑剂的相容性, 结合能的分析揭示了HTPB/增塑剂共混物组分间的相互作用及本质. 将Flory-Huggins作用参数转化为MesoDyn模拟的输入参数, 采用MesoDyn模拟方法对HTPB/增塑剂共混体系的介观形貌与动力学演变过程进行了研究, 通过模拟得到的等密度图、自由能密度和有序度参数等可以判断共混体系的相容性. MD和MesoDyn模拟结果均表明: HTPB/DOS属于相容体系, 而HTPB/NG属于不相容体系, 其结论与实验结果一致.     

端羟基聚丁二烯(HTPB)的~(13)C-NMR研究

<正> 以过氧化氢催化丁二烯聚合而得到的端羟基聚丁二烯(HTPB)在工业上具有广泛的用途,当它们发生交联时,即能形成结构体,用作推进剂、涂料和其它材料的基体,HTPB的质量对其最终产物的力学性能起着决定性的作用。然而人们已经发现,在某些情况下,使用相同牌号的HTPB,甚至用常规的分析方法进行鉴定认为是相同的HTPB,其交联体     

温度对Al2O3吸附诺氟沙星的影响及其分子动力学模拟

采用批量吸附实验及分子动力学(MD)模拟,分别从实验和理论角度研究了温度对氧化铝(Al2O3)吸附诺氟沙星(NOR)性能的影响。结果表明,NOR分子与Al2O3晶体表面存在相互作用,使得Al2O3可有效吸附NOR;对于不同的Al2O3初始投加浓度,在所设定的三个实验温度(293K,298K,303K)中,298K下NOR的吸附去除效率r298K最高,不同温度条件下的吸附去除效率依次为:r298Kr293Kr303K;MD模拟发现NOR分子与Al2O3的(012)晶面间的结合能(Eb)及NOR分子形变能(ΔEdeform)均有以下顺序:NOR-Al2O3(298K)NOR-Al2O3(293K)NOR-Al2O3(303K),该规律与吸附实验结果一致。通过对MD模拟的对关联函数分析,理论上揭示了Al2O3对NOR的吸附主要为非键相互作用。通过对该体系293~303K间各温度的MD模拟,得出Al2O3吸附NOR的理论最佳温度为298K。研究结果为明晰土壤中抗生素在Al2O3界面的微观吸附过程提供了理论依据。     

RhCo双金属催化剂的研究 I. RhCo双金属催化剂的金属-金属及金属-载体相互作用

利用XPS对以金属盐共浸渍制备的Rh+Co/Al2O3和Rh2Co2(CO)12为前体制备的Rh2Co2/Al2O3催化剂的金属-金属及金属-载体相互作用进行了详细研究。发现Rh+Co/Al2O3经400℃H2还原后, Rh的结合能与Rh^0的结合能接近, Co基本上以CoAl2O4状态存在, Co^0的谱峰很弱; 而Rh2Co2/Al2O3经400℃H2还原后, Rh的结合能与Rh^0接近, Co除以CoAl2O4状态存在外, 还有相当一部分以Co^0状态存在。上述结果揭示出两个样品的Rh-Al2O3的相互作用弱, Co-Al2O3的相互作用强。Rh+Co/Al2O3上Rh-Co相互作用弱, 而Rh2Co2/Al2O3在H2还原后仍保持RhCo簇合物的强的Rh-Co相互作用, 导致显著量的Co^0存在。     

RhCo双金属催化剂的研究 I. RhCo双金属催化剂的金属-金属及金属-载体相互作用

利用XPS对以金属盐共浸渍制备的Rh+Co/Al2O3和Rh2Co2(CO)12为前体制备的Rh2Co2/Al2O3催化剂的金属-金属及金属-载体相互作用进行了详细研究。发现Rh+Co/Al2O3经400℃H2还原后, Rh的结合能与Rh^0的结合能接近, Co基本上以CoAl2O4状态存在, Co^0的谱峰很弱; 而Rh2Co2/Al2O3经400℃H2还原后, Rh的结合能与Rh^0接近, Co除以CoAl2O4状态存在外, 还有相当一部分以Co^0状态存在。上述结果揭示出两个样品的Rh-Al2O3的相互作用弱, Co-Al2O3的相互作用强。Rh+Co/Al2O3上Rh-Co相互作用弱, 而Rh2Co2/Al2O3在H2还原后仍保持RhCo簇合物的强的Rh-Co相互作用, 导致显著量的Co^0存在。     

端羟基聚丁二烯/增塑剂共混物相容性的分子动力学模拟

固体推进剂和炸药的力学性能在很大程度上依赖于配方中高分子粘结剂与增塑剂的相容性. 本文对相容和非相容两种体系进行了分子动力学(MD)模拟, 以考察分子模拟方法的实用性. 为预测固体推进剂中端羟基聚丁二烯(HTPB)与增塑剂癸二酸二辛酯(DOS)、硝化甘油(NG)的相容性, 采用MD模拟方法在COMPASS力场下, 对HTPB、DOS、NG和共混物HTPB/DOS、HTPB/NG的密度、内聚能密度及溶度参数等进行了模拟计算. 通过比较溶度参数差值(△δ)的大小、分子间径向分布函数和模拟前后体系密度变化情况均可以预测HTPB/DOS属于相容体系,而HTPB/NG属于不相容体系, 与实验结果一致. 径向分布函数分析同时揭示了HTPB/增塑剂组分之间的相互作用及本质. 本文的模拟方法可以作为预测聚合物与增塑剂相容性的有利工具, 也可以为固体推进剂和炸药的配方设计提供理论指导.     

ε-CL-20基PBX结构和性能的分子动力学模拟——HEDM理论配方设计初探

以高能量密度化合物ε-CL-20(六硝基六氮杂异伍兹烷)为主体,分别添加5种高聚物黏结剂(Estane5703、GAP、HTPB、PEG和F2314)构成高聚物黏结炸药(PBXs).用分子动力学(MD)方法模拟研究各PBX的结合能、相容性、安全性、力学性能和能量性质,通过比较和分析,为优选黏结剂、指导HEDMs配方设计提供信息和规律.由结合能预测各PBX的相容性和稳定性排序为:ε-CL-20/PEG>ε-CL-20/Estane5703≈ε-CL-20/GAP>ε-CL-20/HTPB>ε-CL-20/F2314.以对相关函数g(r)描述了组分之间相互作用的方式.5种黏结剂的少量加入均能显著改善ε-CL-20的弹性力学性能,增强各向同性.各黏结剂并非通过改变ε-CL-20的分子结构影响PBX的感度.它们主要通过自身的热容(C°p)和密度(ρ)影响PBX的安全性和能量性质.     

First-principles study of the atomic structures,electronic properties,and surface stability of BaTiO3 (001) and (011) surfaces

The atomic structures, electronic properties, and surface stability of (001) and (011) surfaces of BaTiO₃ are studied by first-principles calculations. Four differently terminated BaTiO₃ surfaces are considered in this study, including (001)-BaO, (001)-TiO₂, (011)-BaTiO, and (011)-O₂ terminations. The relaxations and rumplings are calculated and discussed, finding that the first layer relaxes inwards, while the second layer relaxes outwards for (001) and (110) surfaces. The data obtained for electronic properties show that O_2p states in (001)-BaO/(001)-TiO₂ termination shift to the lower/higher energy region, leading to a wide/narrow band gap. And the new produced surface states are observed in (011) surface terminations, which is mainly attributed to the supplied electrons from outermost surface atoms, even O atoms are oxidized. Furthermore, the (001) surface of BaTiO₃ is found to be more stable than the (011) surface according to the predicted surface energy which is 0.86 and 2.92 J/m² for (001) and (011) surfaces, respectively. Of which, BaO termination is predicted to be more likely to cleavage from the (001) direction than the TiO₂ termination is.     

Fe2C晶体及低指数晶面的结构与稳定性研究

采用自旋极化的密度泛函理论(DFT)对正交与六方的Fe₂C晶体体相与表面性质进行了研究,计算了晶胞的聚合能、磁矩以及低指数晶面的表面能。研究结果表明,两种晶型Fe₂C 的磁性质相似,但正交堆积的Fe₂C比六方堆积的Fe₂C更稳定。正交晶系Fe₂C低指数晶面的稳定性以 (011) > (110) > (100) > (101) > (001) 顺序降低。对一系列碳化程度不同的碳化铁最稳定表面(Fe₂C(011)、Fe₃C(001)和Fe₄C(100))表面能的比较显示,碳化铁表面的相对稳定性与碳化度非线性相关。另外,与面心立方(BCC)铁最稳定表面(110)相比,Fe₂C、Fe₃C及Fe₄C晶体最稳定表面具有较低的表面能,表明铁表面碳化在热力学上是有利的。     

The relaxation of HCN(011) by V-T,R and V-V energy transfer

Tuned output from an optical parametric oscillator has been used to excite HCN directly to its (011) level. By careful use of a “cold-gas filter”, it has proved possible to distinguish between the time-resolved fluorescence from HCN(011) and that from HCN(001) formed during collisional relaxation. Rate constants for relaxation from both levels have been obtained for the partners: (i) He, Ne, Ar and Kr, and (ii) HCN, CO₂, N₂O, OCS, CS₂, C₂H₂ and C₂D₂. With the rare gases, HCN(011) is relaxed to (001) by V-T,R energy transfer, with rate constants (cm³ molecule^?1 s^?1) at 298 ± 4 K of: k_He⁰¹¹ = (7.9 ± 1.05) × 10^?13; k_Ne⁰¹¹ = (1.56 ± 0.12) × 10^?13; k_Ar⁰¹¹ = (1.20 ± 0.17) × 10^?13; k_Kr⁰¹¹ = (6.7 ± 0.65) × 10^?14. The molecular collision partners also transfer HCN(011) to (001). The rates are much greater and clearly near-resonant V-V energy exchange is important. The results are compared to first-order Sharma-Brau theory, with fair agreement where near-resonant channels exist.     

Crystal nucleation, growth, and morphology of the synthetic malaria pigment beta-hematin and the effect thereon by quinoline additives: the malaria pigment as a target of various antimalarial drugs

The morphology of micrometer-sized beta-hematin crystals (synthetic malaria pigment) was determined by TEM images and diffraction, and by grazing incidence synchrotron X-ray diffraction at the air-water interface. The needle-like crystals are bounded by sharp {100} and {010} side faces, and capped by {011} and, to a lesser extent, by {001} end faces, in agreement with hemozoin (malaria pigment) crystals. The beta-hematin crystals grown in the presence of 10% chloroquine or quinine took appreciably longer to precipitate and tended to be symmetrically tapered toward both ends of the needle, due to stereoselective additive binding to {001} or {011} ledges. Evidence, but marginal, is presented that additives reduce crystal mosaic domain size along the needle axis, based on X-ray powder diffraction data. Coherent grazing exit X-ray diffraction suggests that the mosaic domains are smaller and less structurally stable than in pure crystals. IR-ATR and Raman spectra indicate molecular based differences due to a modification of surface and bulk propionic acid groups, following additive binding and a molecular rearrangement in the environment of the bulk sites poisoned by occluded quinoline. These results provided incentive to examine computationally whether hemozoin may be a target of antimalarial drugs diethylamino-alkoxyxanthones and artemisinin. A variation in activity of the former as a function of the alkoxy chain length is correlated with computed binding energy to {001} and {011} faces of beta-hematin. A model is proposed for artemisinin activity involving hemozoin nucleation inhibition via artemisinin-beta-hematin adducts bound to the principal crystal faces. Regarding nucleation of hemozoin inside the digestive vacuole of the malaria parasite, nucleation via the vacuole's membranous surface is proposed, based on a reported hemozoin alignment. As a test, a dibehenoyl-phosphatidylcholine monolayer transferred onto OTS-Si wafer nucleated far more beta-hematin crystals, albeit randomly oriented, than a reference OTS-Si.     

Molecular dynamics simulations for pure epsilon-CL-20 and epsilon-CL-20-based PBXs

Molecular dynamics has been employed to simulate the well-known high energy density compound epsilon-CL-20 (hexanitrohexaazaisowurtzitane) crystal and 12 epsilon-CL-20-based PBXs (polymer bonded explosives) with four kinds of typical fluorine polymers, i.e., polyvinylidenedifluoride, polychlorotrifluoroethylene, fluorine rubber (F(2311)), and fluorine resin (F(2314)) individually. The elastic coefficients, isotropic mechanical properties (tensile moduli, bulk moduli, shear moduli, and poission's ratios), and bonding energy are first reported for epsilon-CL-20 crystal and epsilon-CL-20-based polymer bonded explosives (PBXs). The mechanical properties of epsilon-CL-20 can be effectively improved by blending with a small amount of fluorine polymers, and the whole effect of the adding fluorine polymers to improve mechanical properties of PBXs along the three crystalline surfaces of epsilon-CL-20 is found to be (100) approximately (001) > (010). The interaction between each of the crystalline surfaces and each of the fluorine polymers is different, and the ordering of binding energy for the three surfaces is (001) > (100) > (010); F(2314) always has the strongest binding ability with the three different surfaces. F(2314) can best improve the ductibility and tenacity of PBX when it is positioned on epsilon-CL-20 (001) crystal surface. The calculations on detonation performances for pure epsilon-CL-20 crystal and the four epsilon-CL-20-based PBXs show that adding a small amount of fluorine polymer into pure epsilon-CL-20 will lower detonation performance, but each detonation parameter of the obtained PBXs is still excellent.     

冲击载荷下TATB晶体滑移和各向异性的分子动力学研究

采用ReaxFF反应力场和分子动力学方法,研究了1,3,5-三氨基-2,4,6-三硝基苯(TATB)炸药晶体在沿不同方向冲击载荷下的滑移和各向异性。冲击方向分别垂直于(101)、(111)、(011)、(110)、(010)、(100)和(001)晶面,冲击强度为10 GPa。研究结果表明,各冲击方向下可能被激发的滑移系均在{001}面,而其它滑移系均因很大的剪切阻力不容易被激发,这与TATB晶体沿c轴的层状结构和平面分子结构相符。预测了七个冲击方向下最容易被激发的滑移系,分别为(101)/{001}100、(111)/{001}010、(011)/{001}010、(110)/{001}010、(010)/{001}110、(100)/{001}120和(001)/{001}010。TATB晶体的冲击响应具有各向异性,动力学过程中体系的应力、能量、温度和化学反应都依赖于冲击方向。对垂直于(100)和(001)晶面的冲击,体系在滑移过程中遭遇的剪切阻力较高、持续时间较长,使得能量和温度较快升高,化学反应较容易发生;对垂直于(101)和(111)晶面的冲击,体系在滑移过程中遭遇的阻力较小且出现次数少,使得能量和温度缓慢升高,化学反应不易发生;对其余冲击方向,体系的响应居中。据此评价了7个冲击方向的相对敏感程度:(101)、(111)(011)、(110)、(010)(100)、(001)。本研究有助于在微观层次深入认识动载荷下TATB的响应机制、结构与性能的关系,为高能低感炸药的设计和研制提供理论参考。     

Hot bands in jet-cooled and ambient temperature spectra of chloromethylene

Rotationally resolved spectra of several bands lying to the red of the origin of the A(1)A" - X (1)A' band system of chloromethylene (HCCl), were recorded by laser absorption spectroscopy in ambient temperature and jet-cooled samples. The radical was made by excimer laser photolysis of dibromochloromethane, diluted in inert gas, at 193 nm. The jet-cooled sample showed efficient rotational but less vibrational cooling. Analysis showed that the observed bands originate in the (upsilon(1),upsilon(2),upsilon(3)) = (010), (001), and (011) vibrational levels of the ground electronic state of the radical, while the upper-state levels involved were (000), (010), (001), and (011). Vibrational energies and rotational constants describing the rotational levels in the lower-state vibrational levels were determined by fitting to combination differences. The analysis also resulted in a reevaluation of the C-Cl stretching frequency in the excited state and we find E(001)' = 13 206.57 or 926.17 cm(-1) above the A(1)A" (000) rotationless level for HC(35)Cl. Scaled ab initio potential energy surfaces for the A and X states were used to compute the transition moment surface and thereby the relative intensities of different vibronic transitions, providing additional support for the assignments and permitting the prediction of the shorter wavelength spectrum. All the observed upper state levels showed some degree of perturbation in their rotational energy levels, particularly in K(a) = 1, presumably due to coupling with near-resonant vibrationally excited levels of the ground electronic state. Transitions originating in the low-lying a(3)A" were also predicted to occur in the same wavelength region, but could not be identified in the spectra.     

NiO表面磁性对分子吸附的影响的第一性原理研究

过渡金属氧化物广泛应用在当今能源与环境相关的催化领域,理解其表面化学性质以及结构-反应活性之间的关系对于先进催化材料的进一步发展以至理性设计至关重要.3d后过渡系金属(Mn,Fe,Co,Ni)的氧化物以其中金属离子独特的自旋状态和由此产生的铁磁/反铁磁性为典型特征.研究过渡金属氧化物的自旋状态以及磁性对表面化学的影响将使我们更加完整了解这些材料的表面化学.以NiO为代表的后过渡系金属岩盐结构一元氧化物具有反铁磁性,被经常作为反铁磁研究的模型体系.尽管在低温(低于其Neel温度)下NiO体相的完整晶体具有确定的反铁磁序,但是一系列最新研究表明,在条件变化时NiO表面的Ni离子可以产生不同的磁序.以此为背景,本工作以NiO为模型体系,采用DFT+U的第一性原理方法研究了NiO表面磁序对表面的小分子吸附活性的影响,包括表面吸附活性对各磁性相的表面取向以及吸附物种磁性的依赖关系.我们考察了NiO的5种反铁磁相和一种铁磁相,两个晶面NiO(001)和NiO(011),顺磁性分子NO和非顺磁性分子CO.我们发现表面能受磁性的影响较轻微,NiO(001)面上从49到54 meV/?2,NiO(011)面上从162到172 meV/?2.在NiO(001)面上,CO与NO都倾向于在Ni离子的顶位吸附.对于不同的体相磁序与表面取向,CO吸附能的变化范围为-0.33~-0.37 eV,NO吸附能的变化范围为-0.42~-0.46 eV.在NiO(011)表面,两种分子都倾向于吸附在由两个Ni离子构成的桥位.我们发现相对于NiO不同磁性相的体相长程磁序,吸附位点处构成桥位的两个Ni离子的局部磁矩相对取向对于分子的吸附具有更加显著的影响.计算得到NO在局部磁矩相对取向反平行(↑↓)吸附位点处的吸附能为-0.99~-1.05 eV,在局部磁矩相对取向平行(↑↑)吸附位点处吸附会增强,吸附能为-1.21~-1.30 eV.对于CO,尽管计算的吸附能在(↑↓)吸附位点(-0.73~-0.75 eV)与在(↑↑)吸附位点(-0.71~-0.72 eV)非常接近,两种吸附位点处的CO吸附时分子轨道杂化方式以及吸附后CO的局域电子态密度却具有明显不同的特征.本工作突出揭示了分子在过渡金属氧化物表面的多重吸附位点上吸附时吸附位点的局域磁矩相对取向对吸附性能的影响.     

MgF2表面结构稳定性及电子特性的密度泛函理论研究

用密度泛函理论(DFT)研究了MgF₂(010)、MgF₂(001)、MgF₂(011)及MgF₂(110)四种表面10种构型的稳定性和电子特性. 结果表明: 四种表面的邻近表面几层原子均出现了明显的驰豫现象, 终止于单层F原子的表面相对稳定; 进一步对比分析四种表面(终止于单层F原子的稳定构型)的表面能发现, 稳定性依次减弱排列为MgF₂(110)、MgF₂(011)、MgF₂(010)、MgF₂(001); 最稳定的MgF₂(110)表面的态密度显示在费米能级以下较多的成键电子处于低能级区, 同时由于表面的影响, 导致表面F原子电荷聚集显负电性, 促使表面活性增加.     

Static and dynamical properties of heavy water at ambient conditions from first-principles molecular dynamics

The static and dynamical properties of heavy water have been studied at ambient conditions with extensive Car-Parrinello molecular-dynamics simulations in the canonical ensemble, with temperatures ranging between 325 and 400 K. Density-functional theory, paired with a modern exchange-correlation functional (Perdew-Burke-Ernzerhof), provides an excellent agreement for the structural properties and binding energy of the water monomer and dimer. On the other hand, the structural and dynamical properties of the bulk liquid show a clear enhancement of the local structure compared to experimental results; a distinctive transition to liquidlike diffusion occurs in the simulations only at the elevated temperature of 400 K. Extensive runs of up to 50 ps are needed to obtain well-converged thermal averages; the use of ultrasoft or norm-conserving pseudopotentials and the larger plane-wave sets associated with the latter choice had, as expected, only negligible effects on the final result. Finite-size effects in the liquid state are found to be mostly negligible for systems as small as 32 molecules per unit cell.