金纳米颗粒周围水的结构和动力学性质的分子动力学模拟

采用经典的分子动力学模拟方法系统地研究了在常温条件下金纳米颗粒周围水的结构与动力学性质. 结果表明, 水分子在纳米颗粒附近形成了明显的多层结构. 同时随着径向距离的减小, 水分子的空间取向也从无序排列趋向于有序排列. 通过分析界面处不同水层中的均方位移及停留时间分布, 发现紧贴颗粒表面的第一和第二水层中的水分子表现出很低的扩散系数, 而第三和第四水层中的水分子则能够轻易地离开界面区域而进入主体相区域. 此外, 在界面处的每个水分子的氢键平均数要高于在主体相的平均值.     

离子液体1-乙基-3-甲基咪唑三氟甲基磺酸盐与水之间的氢键作用(英文)

通过衰减全反射红外(ATR-IR)光谱、二维红外相关谱结合量子化学计算研究了1-乙基-3-甲基咪唑三氟甲基磺酸盐([emim][OTf])和水之间的氢键作用.结果表明,在[emim][OTf]-水体系中,当水的浓度较低时(0.1x(D2O)0.3),水分子的主要存在形式是包裹在离子液体中的没有缔合的单体.水分子优先填充到[emim][OTf]的空隙中,并且与[emim][OTf]的阴离子形成"[OTf]-…HOH…[OTf]-"结构,水分子与[emim][OTf]的阳离子的相互作用位点是烷基氢而不是芳香氢;当水分子浓度较高时,水分子倾向于自身缔合形成小团簇结构,水分子与[emim][OTf]的阳离子的相互作用位点是芳香氢而不是烷基氢.     

聚乙烯在石墨表面结晶的分子动力学模拟

采用分子动力学方法模拟了聚乙烯在石墨(001)表面的吸附和结晶过程;直观的给出了聚乙烯链被石墨(001)面吸附并诱导形成有序的片层晶体的过程;发现结晶温度对得到的有序结构中的聚乙烯链相对石墨表面的特定取向有影响(300 K和600 K时的取向方向不同);表面覆盖率影响聚乙烯吸附层的厚度,对取向的方向无影响.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附.设置不同的水层厚度,观察同液界面和气液界面吸附的差异.模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面,受碳链和固体表面之间相互作用的影响形成表面活性剂分子层,并依据吸附量的大小形成不同的聚集结构;在水层足够厚的情况下,由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构;计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关,解离能垒远大于结合能垒,引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中;无论气液还是固液界面,极性头均伸向水相,与水分子形成不同类型的氢键.模拟表明,分子动力学方法可以作为实验的一种补充,为实验提供必要的微观结构信息.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附. 设置不同的水层厚度, 观察固液界面和气液界面吸附的差异. 模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面, 受碳链和固体表面之间相互作用的影响形成表面活性剂分子层, 并依据吸附量的大小形成不同的聚集结构; 在水层足够厚的情况下, 由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构; 计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关, 解离能垒远大于结合能垒, 引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中; 无论气液还是固液界面, 极性头均伸向水相, 与水分子形成不同类型的氢键. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观结构信息.     

含有二维有序水层的癸二酸锌(Ⅱ)、钴(Ⅱ)超分子配合物的构筑

室温下,在醇水溶液中合成了2个新型配合物,[Zn(phen)3]2·[Zn(C10H16O4)·(H2O)3]·(C10H16O4)2·20H2O(C10H18O4=癸二酸)(1)和[Co(phen)3]2·[Co(H2O)6]·(C10H16O4)3·30H2O(2),并对配合物进行了元素分析分析、红外光谱分析、热重分析以及晶体结构研究.配合物(1)的基本结构单元中含有一个电中性配位单元[Zn(C10H16O4)·(H2O)3]、二个配位阳离子[Zn(phen)3]2+、二个游离的癸二酸根和20个晶格水.Zn原子有两种配位模式,在[Zn(phen)3]2+配位单元中Zn原子与三个邻菲啰啉的六个N原子配位,构成略有畸变的八面体,Zn原子位于八面体的对称中心;在[Zn(C10H16O4)·(H2O)3]配位单元中Zn原子采取五配位的三角双锥构型,两个羧基均采取单齿配位,同一个癸二酸根与相邻的不同Zn原子配位,将相邻的[Zn(C10H16O4)·(H2O)3]配位单元连接起来,自组装得到了无限链状结构.配合物以癸二酸根为模板在ab平面形成了有序水层,该水层由5元,6元水簇,以及其他由羧基参与的各元环组成.与一般常见的6元水簇不同,配合物1中的6元水簇采取了高能量的类似苯环的平面构象.配合物(2)的基本结构单元中含有一个[Co(H2O)6]2+配位阳离子、二个[Co(phen)3]2+配位阳离子、三个游离的癸二酸根和30个晶格水.Co原子也有两种配位模式.在[Co(phen)3]2+配位单元中,Co原子与三个邻菲啰啉的六个N原子配位,构成略有畸变的八面体,Co原子位于八面体的对称中心.在[Co(H2O)6]2+配位单元中Co(II)与6个配位水的氧原子配位,6个配位水中的6个原子和Co(Ⅱ)形成八面体结构.通过水分子之间及水分子和羧基阴离子间的氢键形成了二维有序水层,二维水层中最小的构筑基元为6元环的水,最大的构筑基元为16元环水,每个水环中分别含有6个和16个游离的没有配位的水分子,水分子之间通过很强的氢键作用形成环形的超分子水簇.这些基元片断在二维空间扩展开来,形成二维有序水层.有趣的是,在该水层上通过水分子之间的分子间氢键形成了16元大环水簇,它与两端的六元环共用六条边,每一个16元水环中含有16个游离的没有配位的水分子,水分子之间通过氢键作用形成环形超分子水簇,每个水分子同时作为氢键的给体和受体.[Co(H2O)6]2+和四个羧酸根离子位于十六元环水的中心,与其周围的游离水分子形成了十二个氢键,同时与环外的游离水和羧酸根形成了十个氢键,对大环水聚集体起到稳定作用.     

高岭石-水体系中水分子结构的分子动力学模拟

以Hendricks模型为初始结构, 利用CLAYFF力场对高岭石-水体系进行无晶体学限制的分子动力学模拟. 结果表明, 层间水有三种类型: I型类似于Costanzo提出的“洞水”分子, 其HH矢量(水分子中从一个氢原子位置指向另一个氢原子位置的方向矢量)平行于(001)平面, 而C2轴稍微倾斜于(001)面法线; II型类似于“连接水”, 一个氢氧键指向临近的层间四面体氧形成氢键, 另一个氢氧键与(001)面近似平行; III型水分子在层间近似保持为竖直状, 一个氢与层间四面体氧形成氢键, 而另一个氢与对面层的羟基氧形成氢键. 高岭石羟基氢沿(001)晶面法线的浓度曲线显示一部分羟基指向变为近似平行于(001)面, 羟基氧因此能够暴露出来与层间水分子氢形成氢键. 此外, 模拟中还观察到部分II型水分子氧偏离于层间的平均位置而更靠近四面体层, 这和Costanzo的实验结果一致, 可能是X射线谱图中(002)弱衍射峰出现的原因.     

氧化镧表面反应性对甲烷和氧分子活化的影响

以氧化镧催化剂在甲烷氧化偶联(OCM)反应中的结构敏感性实验研究为基础, 采用周期性密度泛函理论(DFT)计算研究氧化镧(001), (110)和(100)3个晶面及OCM反应物分子甲烷和氧在其上的吸附、 活化和解离. 结果表明, 氧化镧(001), (110)和(100)3个晶面的表面能大小顺序为(110)>(100)>(001), 3个晶面的价带和导带间隙大小顺序为(110)<(100)<(001), 即(001)是3个晶面中最稳定的晶面, 而(110)则是最活泼的晶面. 甲烷分子在氧化镧(001), (110)和(100)晶面上的吸附很弱(0.03 eV), H—CH₃解离吸附能分别为2.16, 0.68和0.90 eV, 解离反应的难易性与晶面的活性顺序一致; 而氧分子在氧化镧(001), (110)和(100)晶面上的分子吸附能分别为-0.04, -0.31和-0.12 eV, 解离吸附能分别为1.22, 0.53和1.52 eV, 即氧化镧晶面结构对氧分子吸附具有明显的影响, 其中, (001)晶面上吸附最弱, (110)晶面上吸附最强, 以致O—O在(110)晶面上可以较低能垒(0.53 eV)解离, 形成亲电的过氧物种. 由于氧分子在氧化镧表面的吸附较甲烷分子强, 因此, 氧化镧在OCM反应中结构敏感性应与氧分子的吸附和活化密切相关. 甲烷和氧分子在氧化镧表面上活化的本质源自于电子自表面流向甲烷和氧分子的反键轨道, 且表面结构的改变会导致不同强度的电子流动驱动.     

含水离子液体/金属界面结构的SERS研究

利用表面增强拉曼光谱(SERS)研究了不同含水量下离子液体及水分子在银电极上随电位变化吸附方式的改变,通过水的O-H伸缩振动谱峰频率变化特征,详细探究了水在离子液体/电极界面上的存在形式及作用方式以及体系零电荷电位与水含量的关系.水含量较低时O-H伸缩振动的Stark系数值较低,随水含量的增加O-H伸缩振动的谱峰位置逐渐向高波数方向移动,同时O-H伸缩振动的Stark系数也逐渐增大,1molL-1[BMIM]Br水溶液中达到76cm-1V-1,且体系的零电荷电位正移,这些差异与水在离子液体中所形成氢键的程度及水分子的存在形式密切相关,在水的含量较低时水与离子液体阳离子通过氢键作用而存在于界面层中,当水的含量增加时,水分子间氢键的作用增强,水与电极表面直接作用的可能性增大.     

Co3O4纳米晶表面晶格氧催化氧化CO的理论研究

近年来,纳米科学技术的迅速发展给催化领域,特别是多相催化带来了新的机遇和挑战.科学家们开始着眼于在纳米尺度上对催化剂结构和催化性能进行表征、控制和设计.Co3O4作为一种重要的半导体金属氧化物材料,由于其优异的氧化还原性质,在锂电池、气体传感器以及多相催化领域得到了十分广泛的应用.最近,研究者发现Co3O4纳米晶在催化CO低温氧化和CH4活化等一系列重要反应中表现出显著的反应活性和晶面效应,表明有效设计和合成特定的高活性、高选择性的纳米晶面,对催化领域的发展将具有十分重要的意义.因而,从原子层面对纳米晶所表现出的这种高活性和晶面效应进行深入解释,将为高效催化剂设计提供重要指导.低温CO氧化作为一种重要的催化反应在燃料电池、空气净化与汽车尾气处理中具有重要的应用价值.本文采用密度泛函理论对Co3O4纳米晶催化CO氧化反应的机理、晶面效应以及结构敏感性进行了理论研究.首先,研究了CO在Co3O4(001)和(011)表面Co,CoOo和Co-Ot三种不同位点的吸附扩散行为,发现CO在Co位点表现出较强的吸附行为,但这种吸附构型需要克服很高的能垒(～1 eV)才能转变到Co-O离子对位点,在低温下这种转变将不可能发生,因此我们推断CO在Co位点的吸附对Co3O4催化CO氧化的晶面效应没有显著影响.接着,对CO在Co-O离子对位点抽提晶格氧生成CO2的反应机理进行了研究.我们发现,(011)表面Co-Ot位点可以较强地吸附CO(吸附能-1.15 eV),并十分容易夺取晶格氧离子(能垒0.26 eV),具有很低的势能面,因而其CO氧化活性明显大于(001)面.为了更清楚地理解这种晶面效应和结构敏感效应的本质,我们提出将CO2形成步的过渡态在反应路径上的能级作为反应活性指标.这种活性指标兼顾考虑了CO在Co-O氧位点的吸附覆盖度和CO2形成步的反应能垒,可以近似理解为反应的表观活化能.据此我们得出,Co3O4不同表面不同品格位点催化CO氧化的反应活性顺序为:(011)-Co-Ot＞＞(001)-Co-Oo＞(011)-Co-Oo＞(001)-Co-Ot.由于CO吸附和CO2形成步都涉及到表面被还原的过程,我们因此发现CO催化氧化活性的高低与表面晶格氧位点的可还原性具有正相关性.这种表面不同位点的还原性可以直接通过对空穴形成能的计算获得,降低表面氧空穴的生成能将有利于提高CO氧化的活性.催化设计的终极目标是在对催化活性位点的本质及反应机理深入认识的基础上在原子层面上对催化剂进行可控设计,从而实现催化剂材料的高效、经济的利用.本文研究表明离子对活性位点是Co3O4纳米晶催化CO氧化反应的活性位点,其中阳离子负责对CO的吸附,阴离子则负责CO2的形成过程,这种协同作用实现了Co3O4纳米晶的高反应活性.我们相信,寻找有效的方法在催化剂表面增加离子对位点活性中心的数目是一种实现高性能催化剂设计的途径.     

Molecular simulation of water behaviors on crystal faces of hydroxyapatite

The water behavior on (001) and (100) crystal faces of hydroxyapatite (HAP) were studied using molecular dynamics (MD) simulations. The study showed that the water molecules between the HAP faces were under conditions of strong electrical field and high pressure, and hence formed 2–3 well-organized water layers on the crystal surfaces. These structured water layers had ice-like features. Compared with the crystallographic [100] direction of HAP, the polarity along the [001] direction was stronger, which resulted in more structured water layers on the surface. The interaction of water molecules with the calcium and phosphate sites at the HAP-water interface was also studied. The results indicated the multiple pathways of water adsorption onto the HAP surfaces. This study revealed the formation and the detailed structure of water layers on HAP surfaces and suggested that the interfacial water played an important role in stabilizing the HAP particles in aqueous solutions. Translated from Chinese Journal of Inorganic Chemistry, 2006, 22(8): 1392–1400 (in Chinese)     

Small molecule‐mediated control of hydroxyapatite growth: Free energy calculations benchmarked to density functional theory

The unique, plate‐like morphology of hydroxyapatite (HAP) nanocrystals in bone lends to the hierarchical structure and functions of bone. Proteins enriched in phosphoserine (Ser‐OPO₃) and glutamic acid (Glu) residues have been proposed to regulate crystal morphology; however, the atomic‐level mechanisms remain unclear. Previous molecular dynamics studies addressing biomineralization have used force fields with limited benchmarking, especially at the water/mineral interface, and often limited sampling for the binding free energy profile. Here, we use the umbrella sampling/weighted histogram analysis method to obtain the adsorption free energy of Ser‐OPO₃ and Glu on HAP (100) and (001) surfaces to understand organic‐mediated crystal growth. The calculated organic‐water–mineral interfacial energies are carefully benchmarked to density functional theory calculations, with explicit inclusion of solvating water molecules around the adsorbate plus the Poisson–Boltzmann continuum model for long‐range solvation effects. Both amino acids adsorb more strongly on the HAP (100) face than the (001) face. Growth rate along the [100] direction should then be slower than in the [001] direction, resulting in plate‐like crystal morphology with greater surface area for the (100) than the (001) face, consistent with bone HAP crystal morphology. Thus, even small molecules are capable of regulating bone crystal growth by preferential adsorption in specific directions. Furthermore, Ser‐OPO₃ is a more effective growth modifier by adsorbing more strongly than Glu on the (100) face, providing one possible explanation for the energetically expensive process of phosphorylation of some proteins involved in bone biomineralization. The current results have broader implications for designing routes for biomimetic crystal synthesis. © 2013 Wiley Periodicals, Inc.     

Adsorption processes of Gly and Glu amino acids on hydroxyapatite surfaces at the atomic level

The regulation mechanism of organic additives on the crystallization of inorganic crystal is fundamentally important in biomineralization. Experimentally, it was found that the amino acids glycine (Gly) and glutamic acid (Glu) could lead to the formation of rod- and plate-like hydroxyapatite (HAP) crystallites, respectively. The detailed adsorption behavior of Gly and Glu on HAP crystal faces was studied by molecular dynamics (MD) simulation. The specific adsorption sites and patterns of Gly and Glu on the (100) and (001) faces of HAP crystals were revealed at the atomic level. The amino acids adsorbed on the HAP (001) and (100) faces with their positive amino groups occupied vacant calcium sites, and their negative carboxylate groups occupied vacant P or OH sites precisely and formed an ordered adsorption layer. The atomic force microscopy pulling simulation and free energy calculation showed that Glu was much more difficult to depart from the HAP (001) face than that from the (100) face. This result indicated that Glu preferred to adsorb strongly onto the HAP (001) face, which resulted in the formation of plate-like HAP. However, Gly did not show any significantly preferential adsorption between these two HAP faces. Thus, the habits of HAP, rod-like crystallites, were not altered during the HAP crystallization in the presence of Gly. Combined with experimental results, our study demonstrated that the MD simulation of interfacial structures could improve our understanding of biological regulation in mineralization processes at the atomic level.     

Atomic force microscope images of the surfaces of aspirin crystals at submolecular resolution.

The atomic force microscope has been developed and used to image arrays of molecules at the (001) and (100) faces of aspirin crystals in water and in air. Lattice spacings composed of methyl groups and the part of the phenyl groups on the surface of the (001) in water, are consistent with X-ray diffraction data. The surface of (100) face which shows most perfect cleavage in bulk, is more difficult to image. This initial success in imaging at drug crystal surfaces clarified the different structural behavior at the submolecular level for three crystal faces, and the close relationship to the differences in the dissolution velocity.     

Chemisorption and dissociation of single oxygen molecules on Ag110

The chemisorption of single oxygen molecules on Ag110 and the dissociation of the adsorbed molecules induced by tunneling electrons were studied at 13 K using a variable-low-temperature scanning tunneling microscope. Two predominant types of chemisorbed O2 molecules were identified, one with the O2 molecular axis aligned along the [001] direction of the substrate [O2(001)], and the other with the molecular axis aligned along the [110] direction [O2(110)]. Tunneling of electrons between the scanning tunneling microscope tip and O2(001) caused the molecule either to rotate or dissociate, depending on the direction of electron tunneling. In contrast, electron tunneling caused O2(110) to dissociate regardless of tunneling direction. In addition to O2(001) and O2(110), several other oxygen species and their dynamical behaviors were observed.     

Crystal structure of {aquaimidazole-[2-(2-carbamoylhydrazone)-propionato]} copper(II) nitrate

The crystal structure of {aquaimidazole[2-(2-carbamoylhydrazone)-propionato]}copper(II) nitrate [Cu(L)Im(H₂O)]NO₃ (I), where HL is the semicarbazone of pyruvic acid, Im is imidazole, is dtermined. The crystal structure of I contains two independent complexes IA and IB in which copper atoms coordinate once deprotonated tridentate HL, imidazole, and water molecules. Outer spheres of the complexes contain nitrate ions. In the compounds studied the coordination polyhedron of the copper atom is a distorted tetragonal pyramid. Its base is composed of carboxyl and carbamide oxygen atoms, azomethine nitrogen of monodeprotonated HL molecules, and the imidazole nitrogen atom. In the crystal, nitrate ions and imidazole molecules link the complexes via hydrogen bonds into 2D networks parallel to the (010) plane. These networks in turn are in pairs arranged into layers along the [010] direction due to hydrogen bonds between water molecules and oxygen atoms of nitrate ions, and also by water molecules and O3 atoms of the neighboring 2D networks. In the crystal, the π-π stacking interaction is observed between the imidazole rings from different layers and there is also a N-O…Cg (π ring) interaction inside the layers.     

Solventless synthesis of monodisperse Cu2S nanorods, nanodisks, and nanoplatelets

Cu(2)S nanocrystals with disklike morphologies were synthesized by the solventless thermolysis of a copper alkylthiolate molecular precursor. The nanodisks ranged from circular to hexagonal prisms from 3 to 150 nm in diameter and 3 to 12 nm in thickness depending on the growth conditions. High resolution transmission electron microscopy (HRTEM) revealed the high chalcocite (hexagonal) crystal structure oriented with the c-axis ([001] direction) orthogonal to the favored growth direction. This disk morphology is thermodynamically favored as it allows the extension of the higher energy [100] and [110] surfaces with respect to the [001] planes. The hexagonal prism morphology also appears to relate to increased C-S bond cleavage of adsorbed dodecanethiol along the more energetic [100] facets relative to [001] facets. Monodisperse Cu(2)S nanodisks self-assemble into ribbons of stacked platelets. This solventless approach provides a new technique to synthesize anisotropic metal chalcogenide nanostructures with shapes that depend on both the face-sensitive thermodynamic surface energy and the surface reactivity.     

Cooperative effects at water-crystalline silica interfaces strengthen surface silanol hydrogen bonding. An ab initio molecular dynamics study

Silica and silica based materials are widely used in chemistry and materials science due to their importance in many technological fields. The properties of these materials, which are crucial for their applications, are mainly determined by the presence of hydrogen bonding between surface silanols. Here, we present ab initio molecular dynamics simulations (AIMD) on different surfaces derived from the crystallographic α-quartz (100) and the α-cristobalite (001) and (101) faces, both free and at the interface with liquid water. The focus was on studying whether water adsorption can disrupt the H-bond pattern at the pristine free silica surface and how deep the perturbation due to the contact with the surface affects the structure of the water multilayer. Results highlight that the water phase is over structured at the interface with silica, as compared to water bulk. Furthermore, an apparent counterintuitive behavior has been observed for quartz (100) and cristobalite (001) surfaces: the interaction with water does not cleave the pre-existent H-bonds between the surface silanol groups. On the contrary, in several cases, it is observed that SiOH···OHSi H-bonds are even strengthened, as the result of a mutual cooperative H-donor/H-acceptor enhancement between silanols and water molecules, which may alter the adsorption capability of these silica surfaces.     

Texture development in polyethylene II. Unidirectional rolling

The effect of degree of strain on texture development in high-density polyethylene has been studied by pole figure analysis for unidirectional rolling. The crystallite orientation distribution in rolling textures has been quantified with an efficient technique which fits three-parameter, two-dimensional Gaussian-type distributions to pole figure intensity data around ideal single crystal orientations. During flat rolling of polyethylene a texture consisting of a strong (100) [001] component and a weak (110) [001] component develops continuously from the lowest true strain of 0.24 (21% reduction) up to the highest true strain of 1.36 (74% reduction). The peak intensity of the Gaussian distributions of both (100) [001] and (110) [001] components increase continuously to the highest strain. The maximum angular breadth of both component distributions, which are roughly perpendicular to the strain direction, remains constant with increasing strain. The minimum angular breadth of both component distributions, which are roughly parallel to the strain direction, decreases continuously owing to gradual alignment of the covalently bonded chain backbone parallel to the strain direction. The development of the (100) [001] component is explained by slip on (100) planes while the weak (110) [001] component is explained by slip on (110) planes. Although the latter component was previously attributed to (110) or (310) relaxation twinning, this seems unlikely because of the lateral constraint during plane strain deformation conditions used in this study.     

The Anisotropic Chemical Reaction Mechanism of 1,3,3-trinitroazetidine (TNAZ) under Different Shock Wave Directions by ReaxFF Reactive Molecular Dynamics Simulations

1,3,3-Trinitroazetidine (TNAZ) has good thermal stability and low shock sensitivity, among other properties, and it has broad prospects in insensitive ammunition applications. In this study, a molecular dynamics calculation based on the ReaxFF-lg force field and multiscale shock technique (MSST) was used to simulate the shock-induced chemical reaction of TNAZ with different shock wave directions. The results showed that the shock sensitivity of TNAZ was in the order of [100] > [010] > [001]. There were significant differences in molecular arrangements in different shock directions, which affected the reaction rate and reaction path in different directions. The molecular arrangement in the [010] and [001] directions formed a “buffer” effect. The formation and cleavage of bonds, formation of small molecules and growth of clusters were analyzed to show the effect of the “buffer”. The polymerization reactions in the [010] and [001] directions appeared later than that in the [100] direction, and the cluster growth in the [010] and [001] directions was slower than that in the [100] direction. In different shock loading directions, the formation and cleavage mechanisms of the N-O bonds of the TNAZ molecules were different, which resulted in differences in the initial reaction path and reaction rate in the three directions