高压下固相硝基甲烷分解的分子动力学计算

基于ReaxFF, 采用NVT系综和Berendsen方法对0–7 GPa时和2500 K时固相硝基甲烷的 分解过程进行分子动力学计算, 通过分析硝基甲烷发生分解反应生成的碎片数量随时间的变化, 对不同压强下硝基甲烷的分解机理进行研究. 计算结果表明在0–3 GPa时, 初始分解路径为C–N键断裂和硝基甲烷的异构化; 在4–7 GPa 时, 初始分解路径为分子间质子转移和C–N, N–O键的断裂; 在硝基甲烷的第二阶段反应中存在H₂O, NO, NO₂, HONO, 硝基甲烷分子自身的催化反应. 硝基甲烷在高温高压下发生热分解反应生成碳团簇, 且团簇中碳原子的数量和碳团簇的空间构型随着压强的变化而变化. 关键词： ReaxFF 分子动力学 热分解 压强效应 碳团簇     

撞击诱导硝基甲烷分子第一步分解反应的动力学计算

本文基于分子温度与压强的关系,计算在不同压强下基态和最低三态硝基甲烷的分子温度,对应计算其沿着CN键裂解反应的热化学和动力学参数.发现基态的硝基甲烷沿着CN键的分解反应是吸热反应,不具自发性,反应转换温度为1550.2 K,平衡常数在80-1202 K温度范围内很低.最低三态的硝基甲烷沿着CN键的裂解是放热反应,反应的Gibbs自由能在80-2558.5 K范围内为负,有好的自发性,且反应较为彻底.298.15-2558.5 K温度范围内反应活化能随着温度的升高而改变,使反应速率随着温度的升高而急剧增大.对应硝基甲烷爆压15 GPa,其分子温度为4617.6 K,该温度下三态分子分解反应的反应速率为1.088×10~8cm~3·mol~(-1)·s~(-1).推算硝基甲烷沿着CN键分解反应混合物的终态温度,当混合物为硝基、甲基和基态的硝基甲烷分子时,反应的终温为1611.37 K,等效能为1676.47 cm~(-1).当混合物为硝基、甲基、基态和最低三态的硝基甲烷分子时,反应的终温为1184.79 K,等效能为1232.65 cm~(-1).两种情况下终态等效能都足以维持硝基甲烷分子沿C-N键裂解反应的发生.这个能量也足以导致混合物中的NO_2分解为NO和O,这与实验检测的结论相一致.     

高压下NH_4ClO_4结构、电子及弹性性质的第一性原理研究

研究高压下NH_4ClO_4的结构和性质对于NH_4ClO_4在固体推进剂和炸药的安全应用具有重要意义.采用基于色散校正密度泛函理论的第一性原理方法,研究了0—15 GPa静水压力下NH_4ClO_4的晶体结构、分子结构、电子性质和弹性性质,计算结果与实验值具有较好的一致性.在压强为1,4和9 GPa时,NH_4ClO_4的晶体参数、键长和分子构型等均出现不连续变化,说明了在压强作用下结构发生变化.随着压强增加,氢键增多且作用增强,由分子内氢键向分子内和分子间的氢键转变;导带态密度峰值增加,电子局域性增强,晶体内N-H和Cl-O共价键作用增强,带隙增大,不同相变区域内带隙呈线性关系.0—15 GPa条件下NH_4ClO_4的弹性常数满足力学稳定性标准,采用Voigt-Reuss-Hill方法计算了体积模量B,剪切模量G和杨氏模量E,根据Cauchy压力和B/G值,说明NH_4ClO_4属于韧性材料,随着压强增加韧性增强.     

LLM-105的分子间相互作用和热力学性质

应用第一性原理可以计算含能材料0 K下的结构和物理性质,但温度效应的缺失通常会导致计算数据与实验结果产生偏差.同时,与温度相关的热力学参数是含能材料在宏观和介观尺度下建模的关键输入.为此,本文以高能低感炸药1-氧-2, 6-二氨基-3, 5-二硝基吡嗪(LLM-105)为研究体系,基于准简谐近似,采用色散修正的密度泛函理论研究温度加载下LLM-105的分子间相互作用和热力学性质.晶格参数和热膨胀系数的演化表明LLM-105分子间相互作用具有强烈的各向异性,其中b轴方向(分子层间)的膨胀率远高于ac平面(分子层内). Hirshfeld表面及其指纹图分析进一步证实LLM-105的分子间相互作用主要取决于O···H构成的氢键.结合Mulliken布居数和结构分析,温度加载下氢键相互作用的变化可诱发硝基旋转,并使得C—NO₂键的强度明显减弱,为高温分解反应的触发键提供了理论依据.此外,本文计算了等容和等压条件下的热容、熵以及等温和绝热条件下的体模量等基础热力学参数.其中绝热条件下的体模量与实验值吻合,同时体模量随温度的演化反映了LLM-105在温度加载下的软化行为.上述...     

压强对冰Ⅶ相结构影响的拉曼光谱研究

测量了室温下2.5～23GPa压强下冰Ⅶ相的拉曼光谱。测量结果表明:随压强增加,冰的氧原子间距离dO—O减小,使氢键长度变短,导致O—H化学键键长增加,力常数减小,拉曼光谱发生红移。质子(氢核)取向有序性随压强先增加而后减小致使拉曼光谱强度先增加而后减小;拉曼光谱线宽先减小而后增加,当压强约为13GPa时呈现最小值。     

高压下硝基甲烷分子温度的密度泛函计算

用密度泛函理论在B3LYP/6-311++G(2d,2P)计算水平上对硝基甲烷分子进行了结构优化、频率和热化学分析.发现:在相同温度条件下改变压强,分子熵函数产生了改变,当温度和压强条件相同时,对于不同物质熵函数的改变是相同的.以热力学理论中麦克斯韦关系为基础,通过计算等温过程中分子的熵函数对压强的变化率,用数值拟合方法得到不同压强条件下分子温度的表达式:T=T0+(1-B)[18.3858+0.5392P]V0,式中T0、V0分别表示分子系统初态的温度和体积,T、V分别表示系统在末态的温度和体积,B是体积的压缩比.在选定参数的情况下该表达式可以计算不同压强条件下CHNO含能材料的分子温度.同时,以硝基甲烷为验证,选取基本参数V0和B,计算其在C-J条件对应的爆压14GPa下,分子温度为3461K,对应爱因斯坦温度,相当于3228cm-1的能量,在实验中该能量足以激发硝基甲烷分子内振动能量重新分配过程,有可能激发C-N键的红外振动而引起单分子分解反应的发生.因此,此表达式可用于预测含能材料撞击点火过程单分子分解可能的反应通道.     

高温下季戊四醇晶体的结构和振动性质的分子动力学模拟(英文)

通过分子动力学研究了季戊四醇晶体在温度高达500K范围内的结构和振动性质.考察了季戊四醇晶体的晶格参数和分子构型等随温度变化情况,相对晶格常数a来说,晶格常数c容易被压缩.发现C-C,C-H和C-O键随温度的变化较大,这意味着在分解过程中这些键可能较敏感,即最先断裂可能性较大.研究了季戊四醇晶体在不同温度下的振动,发现大部分频率随温度的升高而降低.     

BaF2高压相变行为的第一性原理研究

通过理论计算研究了BaF₂在高压下的晶体结构及物理性质.结果表明,在3.5和18.3 GPa,BaF₂依次经历了Fm■m-Pnma-P6₃/mmc两次结构相变,相变过程伴随着体积的塌缩,均为一级相变.约15 GP时,Pnma相晶轴压缩性出现异常,表现为随压强增大,晶轴b_o轻微增加,a_o略微减小.对其电子态密度进行分析发现,在16 GPa以后,由于F1原子的p_y+p_z与p_x轨道电子离域,导致其带隙随压强增加而降低.在约20 GPa时,Pnma相完全转变为P6₃/mmc相,相变完成.对BaF₂的拉曼峰位随压强变化进行了计算,为其高压拉曼光谱行为提供了相应的理论依据.计算了P6₃/mmc相在不同压强下的声子色散曲线,揭示了其卸压过程中的滞后机制,计算结果还预测该物相至少可以稳定到80 GPa.     

碳酸铅在不同传压介质下的高压拉曼研究

为研究PbCO₃在高压下的稳定性,利用金刚石对顶砧技术,采用NaCl固体、甲醇-乙醇-水混合液体(16∶3∶1)和甲醇-乙醇混合液体(4∶1)做传压介质,开展了PbCO₃的高压拉曼实验,最高压强分别达到24.5、25.0和67.0 GPa。研究发现,PbCO₃在10、15和30 GPa左右发生相变,在静水压强条件下CO^2-₃基团的ν^2-外弯曲振动模出现了软化现象。通过对比得到不同传压介质中PbCO₃的Grüneisen参数γ,发现相变机制略有不同,并且压强对晶格振动的影响CO^2-₃基团的影响大,这是由Pb²⁺―O键的键长较大造成的。在所研究的压强范围内,PbCO₃没有发生分解或非晶化,30.0 GPa以上出现的PbCO₃-Ⅳ相直至67.0 GPa都很稳定。     

C轴压力下ZnO晶体结构及电子性质的第一性原理研究

本文采用基于密度泛函理论(DFT)的第一性原理方法对ZnO晶体在c轴取向压力作用下的晶体结构、电子结构的变化进行了研究. 结果表明,当压力在0到6 GPa区间时,晶格参数呈线性变化,带隙随压力增大而增大,显示弹性应变特征;当压力从6 GPa增大到10 GPa的过程中,晶体结构有了较大变化,出现了介于常压下纤锌矿结构和等静压高压下NaCl结构之间的类石墨结构(Graphitelike structure). 伴随着这一结构相变,ZnO的晶格参数,能隙和态密度等电子结构出现了较大跃变.     

Intramolecular and intermolecular hydrogen bonds in aminophenols

IR-Fourier spectroscopy methods are adopted to study intramolecular and intermolecular hydrogen bonds that form in CCl₄ solutions of aminophenol derivatives and in a solid phase of these compounds pressed in KBr. If a hydroxyl group is present in the molecule in the ortho-position to an amino group, then intramolecular interactions between OH and NH groups will take place in aminophenol solutions. Intramolecular O–H⋅⋅⋅O=S=O and N–H⋅⋅⋅O=S=O hydrogen bonds are found in solutions of compounds containing a sulfonamide fragment. Additional acylation of the amino group causes an intramolecular O–H⋅⋅⋅O=C hydrogen bond to form in solutions. Functional groups OH, NH, SO₂, and C=O interact with one another in various ways in the solid phase to form intermolecular hydrogen bonds in aminophenols.     

Theoretical and Experimental Study of the Conformation and Vibrational Frequencies of N‐Acetyl‐L‐alanine and N‐Acetyl‐L‐alaninate

The vibrational frequencies of N‐acetyl‐L‐alanine (NAAL), its potassium salt (NAALK) and its free anionic form (NAAL^?) are calculated using density functional theory (B3LYP) combined with the 6‐311 + + G(d,p) basis set. The experimental Raman spectrum of solid NAALK and the scaling factors for calculated values are discussed as well. The three species are characterized by intramolecular NH…O hydrogen bonds leading to the formation of a five‐membered ring. As indicated by the intramolecular (N)H…O distances and by the ν(NH) frequencies, the strength of the intramolecular hydrogen bond is ordered as follows: NAAL^? < NAALK < NAAL^?. Owing to their difference in the coupling with other vibrational modes, the in‐plane and out‐of‐plane vibrations do not reflect the strength of the hydrogen bond.     

Intramolecular hydrogen bonds in sulfur-containing aminophenols

IR Fourier spectroscopy methods have been adopted to study intramolecular interactions that occur in CCl₄ solutions of antiviral derivatives of aminophenol. Analysis of the IR spectra showed that intramolecular bonds O–H···N, O–H···O=C, N–H···O=S=O, and O–H···O=S=O can occur in these compounds depending on the substituent on the amino group. Not only the presence of intramolecular O–H···N, O–H···O=S=O, and N– H···O=S=O hydrogen bonds in 2-amino-4,6-di-tert-butylphenol derivatives containing a sulfonamide fragment but also conformational equilibrium among these types of intramolecular interactions are essential for the manifestation of high efficiency in suppressing HIV-infection in cell culture.     

Vibrational-electronic coupling in IrF6, OsCl6 2− and IrCl6 2−

The results obtained in a previous paper [1] for an unsymmetric regular model are applied to solid hydrogen at low temperatures when the lattice structure and the quadrupole-quadrupole interaction make the interaction energies non-isotropic. It is shown that, as for the lattice structures with isotropic interactions, no second-order transition is likely to occur through a cooperative rotational effect. The possible occurrence of spatial ordering on sub-lattices is discussed. It is found that such ordering is unlikely to occur at low temperatures on the face-centred cubic lattice, but that on the hexagonal close-packed lattice (which is the probable crystal structure) there is a second-order transition to an ordered state. Using a zeroth-order approximation the temperature at which this transition occurs in pure ortho-hydrogen is found to be 5·8°k.     

A theoretical study on the ESPT mechanism for a novel Bis‐HPBT fluorophore

In this work, based on the density functional theory and time‐dependent density functional theory methods, the properties of the 2 intramolecular hydrogen bonds (O1‐H2···N3 and O4‐H5···N6) of a new photochemical sensor 4‐(3‐(benzo[d]thiazol‐2‐yl)‐5‐tert‐butyl‐4‐hydroxybenzyl)‐2‐(benzo[d]thiazol‐2‐yl)‐6‐tert‐butyl phenol (Bis‐HPBT) have been investigated in detail. The calculated dominating bond lengths and bond angles about these 2 hydrogen bonds (O1‐H2···N3 and O4‐H5···N6) demonstrate that the intramolecular hydrogen bonds should be strengthened in the S₁ state. In addition, the variations of hydrogen bonds of Bis‐HPBT have been also testified based on infrared vibrational spectra. Our theoretical results reproduced absorption and emission spectra of the experiment, which verifies that the theoretical level we used is reasonable and effective in this work. Further, hydrogen bonding strengthening manifests the tendency of excited state intramolecular proton transfer (ESIPT) process. Frontier molecular orbitals depict the nature of electronically excited state and support the ESIPT reaction. According to the calculated results of potential energy curves along stepwise and synergetic O1‐H2 and O4‐H5 coordinates, the potential energy barrier of approximately 1.399 kcal/mol is discovered in the S₁ state, which supports the single ESIPT process along with 1 hydrogen bond of Bis‐HPBT. In other words, the proton transfer reaction can be facilitated based on the electronic excitation effectively. In turn, through the process of radiative transition, the proton‐transfer Bis‐HPBT‐SPT form regresses to the ground state with the fluorescence of 539 nm.     

How universal are hydrogen bond correlations? A density functional study of intramolecular hydrogen bonding in low-energy conformers of α-amino acids

Hydrogen bonding is one of the most important and ubiquitous interactions present in Nature. Several studies have attempted to characterise and understand the nature of this very basic interaction. These include both experimental and theoretical investigations of different types of chemical compounds, as well as systems subjected to high pressure. The O–H..O bond is of course the best studied hydrogen bond, and most studies have concentrated on intermolecular hydrogen bonding in solids and liquids. In this paper, we analyse and characterise normal hydrogen bonding of the general type, D–H...A, in intramolecular hydrogen bonding interactions. Using a first-principles density functional theory approach, we investigate low energy conformers of the twenty α-amino acids. Within these conformers, several different types of intramolecular hydrogen bonds are identified. The hydrogen bond within a given conformer occurs between two molecular groups, either both within the backbone itself, or one in the backbone and one in the side chain. In a few conformers, more than one (type of) hydrogen bond is seen to occur.

Interestingly, the strength of the hydrogen bonds in the amino acids spans quite a large range, from weak to strong. The signature of hydrogen bonding in these molecules, as reflected in their theoretical vibrational spectra, is analysed. With the new first-principles data from 51 hydrogen bonds, various parameters relating to the hydrogen bond, such as hydrogen bond length, hydrogen bond angle, bond length and vibrational frequencies are studied. Interestingly, the correlation between these parameters in these bonds is found to be in consonance with those obtained in earlier experimental studies of normal hydrogen bonds on vastly different systems. Our study provides some of the most detailed first-principles support, and the first involving vibrational frequencies, for the universality of hydrogen bond correlations in materials.     

Intramolecular hydrogen bonds and antioxidant activity of aminophenols

We have used IR Fourier spectroscopy to study intramolecular interactions in solutions of aminophenols in n-hexane. When the hydroxyl group in the molecule is ortho to the amino group, O-H⋯N and N-H⋯O intramolecular hydrogen bonds are formed in the aminophenols. Adding two tert-butyl groups to the benzene ring of ortho-aminophenols strengthens the O-H⋯N bond in the molecules, and prevents formation of an N-H⋯O bond. Additional acylation of the amino group in ortho-aminophenols leads to formation of an O-H⋯O=C intramolecular hydrogen bond. Formation of the above-indicated intramolecular hydrogen bonds in aminophenols affects the course of radiation-induced reactions occurring in n-hexane with participation of these compounds. The antioxidant properties of the aminophenols are enhanced when the hydroxyl groups in the molecules are found in the free state, and are diminished when strong O-H⋯N or O-H⋯O=C intramolecular hydrogen bonds are formed. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 74, No. 5, pp. 577–582, September–October, 2007.     

氨的半水合物的高压拉曼光谱研究

在室温条件下, 利用金刚石对顶砧超高压技术, 对氨的半水合物(2NH3·H2O) 进行了原位高压拉曼光谱研究, 采用红宝石荧光压标测压, 实验的最高压力为41.0 GPa。装入金刚石对顶砧样品腔的初始样品为液态的氨的半水合物, 当压力达到3.5 GPa时, 显微镜下观察到整个样品腔内均匀的出现块状晶体, 同时, 测量到的拉曼谱上出现许多新的拉曼峰。因此, 我们判断在此压力下液态的2NH3·H2O发生了液固相变。当压力增加到19.0 GPa左右时, 2NH3·H2O的拉曼频移随压力变化的曲线有拐点, 并且具有软化特性的N-H伸缩振动模式消失。我们分析这是因为在高压下, 通过O-H…N成键的II型氨分子发生了旋转, 所以2NH3·H2O在此压力下发生了一次固固相变。     

Manifestation of intramolecular and intermolecular interactions in infrared absorption spectra of biologically active aminophenols

FTIR methods were used to study intramolecular and intermolecular interactions in solutions and the solid state of the biologically active aminophenols 2-anilino-4,6-di-tert-butylphenol, N-(3,5-di-tert-butyl-2-methoxyphenyl)aniline, and 2,4-di-tert-butyl-10H-1-phenothiazinol. An analysis of the IR spectra has shown that intramolecular interactions between the OH and NH groups occur in solutions of 2-anilino-4,6-di-tertbutylphenol in CCl₄ to form O-H⋯N hydrogen bonds. The NH groups in solutions of N-(3,5-di-tert-butyl-2-methoxyphenyl)aniline in CCl₄ are present in the non-associated state because of the absence of hydroxyl groups in the molecular structure. The OH and NH groups in solutions of 2,4-di-tert-butyl-10H-1-phenothiazinol in CCl₄ do not interact within the molecule due to a decrease in the conformational mobility of the molecular fragments due to the presence of the rigid C_Ar-S-C_Ar molecular bond in this compound. Intermolecular interactions involving the NH groups occur in the solid state of these compounds. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 75, No. 1, pp. 68–73, January–February, 2008.     

A first-principles study of the structural, electronic and elastic properties of solid nitromethane under pressure

The structural,electronic and elastic properties of solid nitromethane are investigated under pressure by performing first-principles density functional theory(DFT)calculations within the generalized gradient approximation(GGA)and the local density approximation(LDA).The obtained ground state structure properties are found to be consistent with existing experimental and theoretical results.The pressure-induced variations of structure parameters(a,b,c and V)indicate that the solid nitromethane has an anisotropic compressibility,and the compression along the c direction is more difficult than along a and b directions.From the vibration curves of intermolecular bond length and bond angle,we find that the C—N bond is the most sensitive among these bonds under pressure,suggesting that the C—N bonds may be broken first under external loading.The influence of pressure on the electronic properties of solid NM has been studied,indicating that solid NM is an insulating compound with a large indirect band gap and tends to be a semiconductor with increasing pressure.Finally,we predict the elastic constants and their pressure dependence for the solid NM with the bulk modulus,Young’s modulus,shear modulus and the Poisson’s ratio derived.