基于含時(shí)密度泛函理論的紫外光對(duì)太安的定性影響

張寶森，張樹(shù)海，茍瑞君，陳亞紅，朱雙飛，馬坤

基于含時(shí)密度泛函理論的紫外光對(duì)太安的定性影響

張寶森1，張樹(shù)海1，茍瑞君1，陳亞紅1，朱雙飛1，馬坤2

（1.中北大學(xué) 環(huán)境與安全工程學(xué)院，太原 030051；2.陜西應(yīng)用物理化學(xué)研究所，西安 710061）

研究太安炸藥在紫外光作用下的穩(wěn)定性退化機(jī)制。基于TDDFT（含時(shí)密度泛函）理論，在pbe1pbe/6-311G**水平下對(duì)太安分子50個(gè)激發(fā)態(tài)進(jìn)行計(jì)算，依據(jù)計(jì)算結(jié)果，繪制吸收光譜，使用空穴–電子方法對(duì)最大吸收峰3個(gè)激發(fā)態(tài)（9、10和11）的激發(fā)特征進(jìn)行分析，這3個(gè)激發(fā)態(tài)對(duì)最大吸收峰的總貢獻(xiàn)率達(dá)97.31%。將此3個(gè)激發(fā)態(tài)定為研究對(duì)象，對(duì)太安分子被特定的紫外光激發(fā)至激發(fā)態(tài)后弱鍵的Mayer和Laplace鍵級(jí)進(jìn)行分析，并基于IFCT（Interfragment Charge Transfer）方法對(duì)太安分子激發(fā)至激發(fā)態(tài)過(guò)程中的電子轉(zhuǎn)移情況進(jìn)行描述。太安紫外吸收光譜的最大吸收峰位置為186.6 nm，小于實(shí)驗(yàn)測(cè)試的吸收峰位置8.4 nm。對(duì)此吸收峰的強(qiáng)度貢獻(xiàn)最大的3個(gè)激發(fā)態(tài)中，9態(tài)貢獻(xiàn)最高，為48.27%，其余2個(gè)激發(fā)態(tài)10、11為簡(jiǎn)并態(tài)，貢獻(xiàn)相同，為24.52%。通過(guò)空穴–電子的電荷轉(zhuǎn)移分析結(jié)果可知，3個(gè)激發(fā)態(tài)均存在整體激發(fā)并帶有局域電荷轉(zhuǎn)移的特征。PETN分子在吸收特定波長(zhǎng)（187.00、186.92 nm）紫外光并激發(fā)至激發(fā)態(tài)時(shí)，O—NO2鍵的Mayer與Laplace鍵級(jí)均有所降低。通過(guò)IFCT分析可知，引發(fā)鍵鍵級(jí)變化由O—NO2上的n→Pi*躍遷主導(dǎo)，這種效應(yīng)會(huì)促使太安分子的穩(wěn)定性降低。

太安；紫外光；含時(shí)密度泛函；吸收光譜；分子穩(wěn)定性；激發(fā)態(tài)

季戊四醇四硝酸酯（太安，PETN）因其具有較高的化學(xué)穩(wěn)定性和爆轟性能，被廣泛用于雷管、導(dǎo)爆索以及傳爆藥的制造之中[1]。近年來(lái)，激光起爆技術(shù)以其安全、可靠、低成本和對(duì)靜電不敏感等傳統(tǒng)技術(shù)無(wú)法代替的特點(diǎn)在國(guó)內(nèi)外火工品中得到廣泛應(yīng)用[2-5]。

1996年，Tarzhanov等[6]對(duì)1.06 μm激光起爆PETN的實(shí)驗(yàn)結(jié)果進(jìn)行了分析，分析結(jié)果表明，1.06 μm激光作用下的PETN響應(yīng)實(shí)驗(yàn)結(jié)果符合光熱起爆理論。2008年，陜西應(yīng)用物理化學(xué)研究所[7]對(duì)多種配合物起爆藥、炸藥的吸收光譜及其不同波長(zhǎng)下的激光感度進(jìn)行了測(cè)試，測(cè)試結(jié)果表明，可見(jiàn)紫外波段激光或?qū)⑵茐幕衔锏娜蹑I，并引發(fā)化學(xué)反應(yīng)。同時(shí)，激光敏感藥劑的感度選擇性調(diào)控，可以通過(guò)化合物輻射躍遷的選擇定律和弱化合鍵的分解機(jī)理來(lái)實(shí)現(xiàn)。2011年，Aluker等[8]對(duì)釹玻璃激光作用添加光散射劑后的PETN樣品進(jìn)行了研究，研究表明，MgO光散射劑的添加會(huì)顯著提升樣品對(duì)激光能量的轉(zhuǎn)換概率。2013年，中國(guó)工程物理研究院化工材料研究所[9]采用實(shí)驗(yàn)結(jié)合量子化學(xué)計(jì)算的方法，對(duì)多種炸藥的激光敏感波長(zhǎng)進(jìn)行了分析，其中PETN最弱化學(xué)鍵為O—NO2鍵，實(shí)驗(yàn)測(cè)得PETN最敏感的紫外激光波長(zhǎng)值為195 nm。2014年，Tsyshevsky等[10]對(duì)PETN光學(xué)吸收光譜中的能帶進(jìn)行了分析，強(qiáng)度較高能帶對(duì)應(yīng)于一系列的單線態(tài)–單線態(tài)躍遷，而強(qiáng)度較低的能帶歸因于緊密結(jié)合的單線態(tài)和三線態(tài)激子，所有躍遷過(guò)程都位于PETN的O—NO2基團(tuán)。2017年，Tarzhanov等[11]研究了Al粉與PETN混合物的激光起爆特性，研究表明，Al粉與PETN粉混合形成的沉淀物的起爆難度較高。2019年，Aduev等[12-13]對(duì)PETN/Al復(fù)合物的激光起爆機(jī)理進(jìn)行了研究，結(jié)果表明，樣品中Al的添加會(huì)降低樣品激光起爆的臨界能量密度，當(dāng)復(fù)合物中含0.1%（質(zhì)量分?jǐn)?shù)）鋁時(shí)，樣品的臨界能量密度最低。2020年，北京理工大學(xué)嚴(yán)楠團(tuán)隊(duì)[14]使用仿真模擬與實(shí)驗(yàn)相結(jié)合的方法，對(duì)半導(dǎo)體激光直接起爆炸藥的可行性進(jìn)行了研究，驗(yàn)證了半導(dǎo)體激光直接起爆炸藥的可行性。

近紅外激光起爆炸藥，通過(guò)添加摻雜物的方法提升炸藥的激光感度，以激光熱效應(yīng)完成起爆過(guò)程。可見(jiàn)紫外激光是通過(guò)選擇性激發(fā)炸藥分子的電子能級(jí)，降低化合物中弱化學(xué)鍵的穩(wěn)定性，誘導(dǎo)化學(xué)反應(yīng)的發(fā)生。因此，研究激光器紫外–可見(jiàn)波長(zhǎng)與藥劑吸收光譜之間的匹配技術(shù)是提高激光特征感度的有效途徑之一[7]。本文將借助TDDFT理論，對(duì)炸藥分子激發(fā)態(tài)進(jìn)行計(jì)算，對(duì)比激發(fā)前后分子引發(fā)鍵強(qiáng)度，以研究紫外光作用下炸藥的穩(wěn)定性退化機(jī)制。在分子層面，研究PETN在紫外光下的定性響應(yīng)將為特征感度尋找提供數(shù)據(jù)支撐，為含能材料在激光環(huán)境下的防護(hù)理論依據(jù)。

1 計(jì)算細(xì)節(jié)及結(jié)構(gòu)來(lái)源

分子對(duì)光的吸收過(guò)程涉及振動(dòng)態(tài)和電子態(tài)的改變，在實(shí)際的計(jì)算過(guò)程中，往往只考慮電子的勢(shì)能面，而忽視核振動(dòng)的量子效應(yīng)，其目的是簡(jiǎn)化計(jì)算。本文計(jì)算的PETN分子結(jié)構(gòu)來(lái)源于CCDC晶體數(shù)據(jù)庫(kù)[15]，使用Gaussian09[16]程序在pbe1pbe/6-311G**水平下對(duì)PETN分子結(jié)構(gòu)進(jìn)行優(yōu)化，結(jié)果無(wú)虛頻。借助TDDFT理論，在pbe1pbe/6-311G**水平下基于優(yōu)化后結(jié)構(gòu)（如圖1所示），對(duì)PETN的50個(gè)激發(fā)態(tài)進(jìn)行計(jì)算。使用Multiwfn[17]軟件，借助展寬函數(shù)和FWHM（Full Width at Half Maximum）對(duì)PETN 50個(gè)激發(fā)態(tài)下的紫外吸收光譜進(jìn)行繪制，并與實(shí)驗(yàn)測(cè)試的PETN吸收光譜最大吸收峰位置[9]進(jìn)行對(duì)比。

為明確PETN分子的躍遷特征，分析了激發(fā)態(tài)對(duì)吸收峰的貢獻(xiàn)值。使用空穴–電子[18]的方法對(duì)貢獻(xiàn)最高的3個(gè)激發(fā)態(tài)的電子和空穴分布進(jìn)行了繪制，空穴–電子分析方法已多次用于含能材料的電子激發(fā)態(tài)研究[19-20]。為研究PETN激發(fā)過(guò)程中的電荷轉(zhuǎn)移情況，進(jìn)一步使用IFCT[17]方法對(duì)PETN片段間電子轉(zhuǎn)移量進(jìn)行計(jì)算，并使用Gaussian09軟件對(duì)3個(gè)激發(fā)態(tài)密度進(jìn)行計(jì)算。通過(guò)Multiwfn軟件生成自然軌道，以對(duì)不同態(tài)下O—NO2鍵的Mayer與Laplace鍵級(jí)進(jìn)行計(jì)算。

圖1 PETN分子優(yōu)化后結(jié)構(gòu)及原子編號(hào)

2 結(jié)果及分析

2.1 PETN分子紫外吸收光譜與激發(fā)特征分析

使用Gaussian09軟件，基于優(yōu)化后的PETN基態(tài)結(jié)構(gòu)，依據(jù)含時(shí)密度泛函理論，對(duì)其50個(gè)激發(fā)態(tài)進(jìn)行了計(jì)算。借助Multiwfn，依據(jù)激發(fā)態(tài)計(jì)算結(jié)果，繪制PETN紫外吸收光譜（如圖2所示），其中FWHM 參數(shù)設(shè)定為0.7 eV，計(jì)算水平為pbe1pbe/6-311G**。

圖2 基于TDDFT方法計(jì)算的PETN分子的紫外吸收光譜

由圖2可知，PETN分子計(jì)算光譜的最大吸收峰位置為186.6 nm。在以往的研究中，實(shí)驗(yàn)測(cè)得PETN炸藥在195 nm處[9]的吸光度最高，與本文所計(jì)算的紫外吸收光譜最大吸收峰位置僅相差8.4 nm。由此可知，本文計(jì)算的紫外吸收光譜與實(shí)驗(yàn)測(cè)試的吸收光譜符合較好。

PETN分子吸收光譜中，對(duì)最大吸收峰起到主要貢獻(xiàn)的3個(gè)激發(fā)過(guò)程為0→9、0→10、0→11，貢獻(xiàn)率分別為48.27%、24.52%、24.52%，振子強(qiáng)度分別為0.144、0.073、0.073。其中，振子強(qiáng)度表示分子吸收某一波長(zhǎng)的光的概率[21]。在3種激發(fā)態(tài)中，10與11對(duì)波峰的貢獻(xiàn)率相同，并且2個(gè)態(tài)的激發(fā)能相等，均為6.633 0 eV，由此可以確定10和11為簡(jiǎn)并態(tài)。

計(jì)算3個(gè)激發(fā)態(tài)的軌道躍遷貢獻(xiàn)（大于5%），結(jié)果見(jiàn)表1。其中10與11態(tài)的2類(lèi)軌道躍遷貢獻(xiàn)率相近，2個(gè)激發(fā)態(tài)的激發(fā)能相等。LUMO+1與LUMO+2能量分別為–2.130 1、–2.130 0 eV，2個(gè)分子軌道能量數(shù)值上近似相等，由此可將這2軌道視為簡(jiǎn)并軌道。

表1 3個(gè)主要貢獻(xiàn)的激發(fā)過(guò)程的主要軌道躍遷貢獻(xiàn)率

Tab.1 Main orbital transition contribution rates in excitation processes of excited states S9, S11 and S10

空穴–電子方法適用于描述大部分體系的激發(fā)過(guò)程。為明確以上激發(fā)態(tài)的激發(fā)特征，本文將采用空穴–電子方法來(lái)分析PETN的電子激發(fā)特征?？昭èC電子分析是基于空穴和電子的密度分布，“空穴”是指電子供體，是電子被激發(fā)和轉(zhuǎn)移的地方；而“電子”表示電子受體。

表2中指數(shù)表示空穴和電子質(zhì)心之間的距離；r指數(shù)表示電子和空穴的重疊程度；指數(shù)表示電子空穴的總體平均分布廣度；指數(shù)衡量空穴和電子的分離程度；coul表示空穴與電子間的庫(kù)侖吸引能；空穴離域指數(shù)（HDI）和電子離域指數(shù)（EDI）表示空穴或電子的分布均勻程度。由表2可知，3種激發(fā)態(tài)指數(shù)、r指數(shù)數(shù)值均較小，而指數(shù)大于0.25 nm。通過(guò)以上3種參數(shù)可以確定，3種激發(fā)過(guò)程具有整體激發(fā)特征。3種激發(fā)態(tài)的指數(shù)為負(fù)值，表示在電荷轉(zhuǎn)移激發(fā)方向上，空穴與電子的分布未發(fā)生顯著的分離。HDI與EDI數(shù)值相近，表示空穴和電子的分布均勻程度類(lèi)似，這與以上分析相符合。觀察3種激發(fā)態(tài)的空穴–電子分布圖（圖3）可知，3種激發(fā)過(guò)程均伴隨局部的電荷轉(zhuǎn)移（與硝基相連的O和與中心C原子連接的亞甲基上電子主要向硝基處移動(dòng)）特征。

表2 基態(tài)到9、11、10態(tài)的、r、、、coul、HDI和EDI指數(shù)數(shù)值

Tab.2 Indexes of D, Sr, H, t, Ecoul, HDI and EDI from the ground state to excited states S9, S11, and S10

圖3 在等值面數(shù)值為0.002（a. u.）時(shí)基態(tài)到S9、S10、S11的空穴（深色）和電子（淺色）分布

從圖3可以看出，大部分電子位于PETN的硝基部分，而空穴位于C原子的兩側(cè)，0→9態(tài)的整體激發(fā)特征更為明顯，并且10、11態(tài)的空穴–電子分布具有明顯的對(duì)稱(chēng)特征，這與上文所分析的2個(gè)激發(fā)態(tài)為簡(jiǎn)并態(tài)相對(duì)應(yīng)。在指數(shù)分析中發(fā)現(xiàn)，除EDI參數(shù)外，其他參數(shù)均相等，同樣體現(xiàn)了2個(gè)激發(fā)態(tài)簡(jiǎn)并。同時(shí)，空穴–電子分布圖明顯可見(jiàn)PETN單分子的3種激發(fā)態(tài)都具有從中心向四周電荷轉(zhuǎn)移激發(fā)的特征，并且空穴集中于與硝基相連的O原子和與中心C原子相連的亞甲基。如圖3所示的局部電荷轉(zhuǎn)移可能會(huì)引發(fā)PETN內(nèi)O—NO2鍵強(qiáng)度的變化。

2.2 PETN主要激發(fā)過(guò)程分析

炸藥中最弱鍵斷裂通常被認(rèn)為是炸藥起爆的一個(gè)重要特征[22]。一般認(rèn)為，PETN的引發(fā)鍵為O—NO2鍵[23-25]，研究PETN激發(fā)態(tài)中O—NO2鍵的穩(wěn)定性具有重要意義，同時(shí)引發(fā)鍵穩(wěn)定性的變化也是PETN吸收特定波長(zhǎng)并激發(fā)至激發(fā)態(tài)的重要響應(yīng)體現(xiàn)。2.1節(jié)中對(duì)PETN主要貢獻(xiàn)的3個(gè)激發(fā)態(tài)的電子激發(fā)特征進(jìn)行了分析，其中PETN中局域電荷轉(zhuǎn)移或許會(huì)改變O—NO2的強(qiáng)度。在以往的研究中認(rèn)為，在特征波長(zhǎng)的作用下，會(huì)降低炸藥的起爆能量[8]。因此，本文對(duì)PETN分子3種激發(fā)態(tài)O—NO2鍵的Mayer鍵級(jí)與Laplace鍵級(jí)進(jìn)行了計(jì)算，計(jì)算結(jié)果見(jiàn)表3。

PETN基態(tài)O—NO2的Mayer鍵級(jí)與Laplace鍵級(jí)分別為0.884和0.285。由表3可知，幾種激發(fā)態(tài)下，PETN分子最弱鍵的鍵級(jí)均有所降低，即PETN分子在吸收波長(zhǎng)為187.00 nm（9）和186.92 nm（10、11）的光后，激發(fā)至激發(fā)態(tài)時(shí)，分子內(nèi)O—NO2鍵的穩(wěn)定性被削弱。

表39、10、11激發(fā)態(tài)下PETN分子O—NO2鍵的Mayer鍵級(jí)、Laplace鍵級(jí)

Tab.3 O—NO2 Mayer and Laplace bond orders of PETN molecules in excited states S9, S11, and S10

2.1節(jié)的分析中，3種激發(fā)態(tài)均包含一定程度的整體激發(fā)和電荷轉(zhuǎn)移激發(fā)特征，3種主要激發(fā)過(guò)程中，PETN分子內(nèi)電子的轉(zhuǎn)移可能影響PETN最弱鍵的穩(wěn)定性。通過(guò)分析O—NO2鍵級(jí)發(fā)現(xiàn)，3種主要貢獻(xiàn)的激發(fā)過(guò)程都會(huì)降低O—NO2的穩(wěn)定性。為對(duì)此過(guò)程進(jìn)行進(jìn)一步分析，使用IFCT方法對(duì)PETN的激發(fā)過(guò)程中片段間的電子轉(zhuǎn)移進(jìn)行了計(jì)算。IFCT分析基于空穴與電子貢獻(xiàn)來(lái)衡量體系中片段間電子的轉(zhuǎn)移量。將PETN分子內(nèi)的4個(gè)硝基和與硝基相連的4個(gè)O原子，以及與1號(hào)C原子相連的4個(gè)亞甲基視為研究對(duì)象（原子編號(hào)如圖1所示），IFCT計(jì)算結(jié)果見(jiàn)表4。

表4 3種激發(fā)過(guò)程下PETN分子內(nèi)片段上電子的變化

Tab.4 Variation of electrons on fragments in PETN molecules under three excitation processes

表4中片段1、4、7、10表示PETN的硝基，2、5、8、11表示亞甲基，3、6、9、12表示與硝基相連的O原子。在9態(tài)中，PETN硝基基團(tuán)均得到0.194e，與9態(tài)明顯由中心向四周CT的特征相吻合，除硝基外，PETN其余基團(tuán)都發(fā)生了電子的轉(zhuǎn)移，并且與硝基相連的O原子的電子轉(zhuǎn)移占主導(dǎo)（61.7%）。在10態(tài)中，PETN分子4、10號(hào)硝基得到0.368e，6、12號(hào)氧原子各轉(zhuǎn)移0.196e，5、11號(hào)亞甲基各轉(zhuǎn)移0.119e。11態(tài)與10態(tài)類(lèi)似，1、7號(hào)硝基基團(tuán)得到0.368e，而3、9號(hào)氧原子各轉(zhuǎn)移0.196e，2、8號(hào)PETN亞甲基各轉(zhuǎn)移0.119e。PETN其他片段如5、11號(hào)亞甲基各轉(zhuǎn)移0.029e，6、12號(hào)氧原子各轉(zhuǎn)移0.044e，電荷轉(zhuǎn)移數(shù)值較小，對(duì)分子內(nèi)總的電子轉(zhuǎn)移貢獻(xiàn)程度較低。3種激發(fā)態(tài)中，電子主要由與硝基相連的O原子提供，而電子主要轉(zhuǎn)移至硝基上。

以上對(duì)PETN各片段間的電子轉(zhuǎn)移分析表明，電子主要來(lái)自于PETN分子內(nèi)與硝基相連的氧原子。進(jìn)一步使用IFCT的方法對(duì)3個(gè)激發(fā)態(tài)中與硝基相連的O原子與硝基之間的電子凈轉(zhuǎn)移進(jìn)行了分析。在9態(tài)中，每個(gè)與硝基相連的O原子向此硝基凈轉(zhuǎn)移的電荷為0.092 1e，這與空穴–電子的分析結(jié)果相吻合；在10態(tài)中，3、9號(hào)O原子分別向4、10號(hào)硝基基團(tuán)凈轉(zhuǎn)移的電子均為0.019 6e；在11態(tài)中，6、12號(hào)O原子分別向1、7號(hào)硝基轉(zhuǎn)移0.019 6e，都與空穴–電子分布圖分析相吻合。觀察圖3可知，與硝基相連的O原子上，空穴分布在氧原子兩側(cè)，同時(shí)硝基上的電子分布具有明顯的對(duì)稱(chēng)面。結(jié)合IFCT分析可知，氧原子與硝基的n→pi*躍遷，促使與硝基相連的O原子發(fā)生電子轉(zhuǎn)移，并導(dǎo)致O原子攜帶的電荷降低，同時(shí)這也是PETN內(nèi)O—NO2鍵削弱的重要原因。

綜上所述，PETN分子內(nèi)0→9、0→10和0→11這3個(gè)激發(fā)過(guò)程，以n→pi*躍遷特征為主導(dǎo)的電子轉(zhuǎn)移會(huì)導(dǎo)致PETN引發(fā)鍵穩(wěn)定性的降低，即PETN在吸收特定的紫外光并激發(fā)至激發(fā)態(tài)（9、10、11）會(huì)促使分子的穩(wěn)定性降低。

3 結(jié)論

本文基于TDDFT理論對(duì)PETN在紫外光下的穩(wěn)定性退化機(jī)制進(jìn)行了研究，首先基于PETN在pbe1pbe/6-311G**水平下的50個(gè)激發(fā)態(tài)計(jì)算繪制吸收光譜，并與實(shí)驗(yàn)測(cè)試數(shù)據(jù)對(duì)比，其中理論計(jì)算峰位置小于實(shí)驗(yàn)吸收峰位置8.4 nm。

使用空穴–電子方法，對(duì)PETN紫外吸收光譜貢獻(xiàn)較高的3種激發(fā)態(tài)的激發(fā)特征進(jìn)行了分析。3種激發(fā)態(tài)的空穴集中于與硝基相連的O原子和與中心C原子相連的亞甲基，激發(fā)過(guò)程可能會(huì)促使PETN內(nèi)O—NO2鍵強(qiáng)度的變化。3種激發(fā)態(tài)均存在一定程度整體激發(fā)和由中心向四周電荷轉(zhuǎn)移激發(fā)的特征。

對(duì)PETN分子激發(fā)態(tài)與基態(tài)的引發(fā)鍵O—NO2的Mayer與Laplace鍵級(jí)進(jìn)行了計(jì)算，計(jì)算結(jié)果表明，3種起到主要貢獻(xiàn)的激發(fā)過(guò)程都會(huì)降O—NO2鍵的強(qiáng)度，即PETN在吸收特定波長(zhǎng)的紫外光后，激發(fā)至激發(fā)態(tài)，這一過(guò)程會(huì)削弱引發(fā)鍵的穩(wěn)定性。

使用IFCT方法對(duì)PETN分子片段間的電子轉(zhuǎn)進(jìn)行了分析，分析結(jié)果表明，與硝基相連的氧原子和硝基之間的n→pi*躍遷是PETN引發(fā)鍵穩(wěn)定性削弱的重要原因。

[1] 董海山, 周芬芬. 高能炸藥及相關(guān)物性能[M]. 北京: 科學(xué)出版社, 1989.

DONG Hai-shan, ZHOU Fen-fen. Performance of High Explosives and Correlates[M]. Beijing: Science Press, 1989.

[2] 王宇銳, 范新中, 王曉鵬, 等. 激光火工品系統(tǒng)在戰(zhàn)術(shù)導(dǎo)彈領(lǐng)域的應(yīng)用前景與展望[J]. 導(dǎo)彈與航天運(yùn)載技術(shù), 2021(6): 15-20.

WANG Yu-rui, FAN Xin-zhong, WANG Xiao-peng, et al. Application Prospect of Laser Ignition System in Tactical Missile Field[J]. Missiles and Space Vehicles, 2021(6): 15-20.

[3] YANG L C, MENICHELLI V J. Optically Detonated Explosive Device: US3812783[P]. 1974-05-28.

[4] 王宇, 魏超, 張嵩. 鈍感火工品中新技術(shù)、新含能材料研究進(jìn)展[J]. 科技導(dǎo)報(bào), 2013, 31(10): 74-79.

WANG Yu, WEI Chao, ZHANG Song. Developing Status of New Techniques and New Energetic Materials in Insensitive Pyrotechnics[J]. Science ＆ Technology Review, 2013, 31(10): 74-79.

[5] 梅崴, 徐宇新, 邢朝洋, 等. 激光點(diǎn)火系統(tǒng)用1×2 MEMS光開(kāi)關(guān)研究[J]. 導(dǎo)航與控制, 2017, 16(1): 85-90.

MEI Wei, XU Yu-xin, XING Chao-yang, et al. Research of 1×2 MEMS Optical Switch for Laser Ignition System[J]. Navigation and Control, 2017, 16(1): 85-90.

[6] TARZHANOV V I, ZINCHENKO A D, SDOBNOV V I, et al. Laser Initiation of PETN[J]. Combustion, Explosion and Shock Waves, 1996, 32(4): 454-459.

[7] 盛滌倫, 朱雅紅, 陳利魁, 等. 激光與含能化合物相互作用機(jī)理研究[J]. 含能材料, 2008, 16(5): 481-486.

SHENG Di-lun, ZHU Ya-hong, CHEN Li-kui, et al. Interactional Mechanism between Laser and Energetic Compound[J]. Chinese Journal of Energetic Materials, 2008, 16(5): 481-486.

[8] ALUKER E D, KRECHETOV A G, MITROFANOV A Y, et al. Laser Initiation of Energetic Materials: Selective Photoinitiation Regime in Pentaerythritol Tetranitrate[J]. The Journal of Physical Chemistry C, 2011, 115(14): 6893-6901.

[9] 劉建, 熊鷹, 蔣小華, 等. 含能材料激光誘導(dǎo)起爆特征波長(zhǎng)解析[J]. 激光技術(shù), 2013, 37(6): 816-819.

LIU Jian, XIONG Ying, JIANG Xiao-hua, et al. Characteristic Wavelength Analysis for Laser-Induced Initiation in Energetic Material[J]. Laser Technology, 2013, 37(6): 816-819.

[10] TSYSHEVSKY R V, SHARIA O, KUKLJA M M. Energies of Electronic Transitions of Pentaerythritol Tetranitrate Molecules and Crystals[J]. The Journal of Physical Chemistry C, 2014, 118(18): 9324-9335.

[11] TARZHANOV V I, SDOBNOV V I, ZINCHENKO A D, et al. Laser Initiation of Mixtures of PETN and Aluminum by a Deposit[J]. Combustion, Explosion, and Shock Waves, 2017, 53(6): 724-729.

[12] ADUEV B P, NURMUKHAMETOV D R, NELYUBINA N V, et al. Initiation Thresholds and Dynamic Characteristics of an Explosion for Thin Samples of PETN-Al Composites under Laser Irradiation[J]. Technical Physics, 2019, 64(6): 858-864.

[13] ADUEV B P, NURMUKHAMETOV D R, ZVEKOV A A, et al. Laser Initiation of PETN-Based Composites with Additives of Ultrafine Aluminium Particles[J]. Combustion, Explosion, and Shock Waves, 2016, 52(6): 713-718.

[14] 湛贊, 譚明, 嚴(yán)楠, 等. 半導(dǎo)體激光起爆炸藥試驗(yàn)與仿真研究[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2020, 40(5): 10-14.

ZHAN Zan, TAN Ming, YAN Nan, et al. Experimental and Simulation Study of Laser Diode Detonation[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2020, 40(5): 10-14.

[15] CADY H H, LARSON A C. Pentaerythritol Tetranitrate II: Its Crystal Structure and Transformation to PETN I; An Algorithm for Refinement of Crystal Structures with Poor Data[J]. Acta Crystallographica Section B Structural Crystallography and Crystal Chemistry, 1975, 31(7): 1864-1869.

[16] FRISCH M J, TRUCKS G W, SCHLEGEL H B, et al. 2009. Gaussian 09[R]. Inc. USA: Wallingford CT.

[17] LU Tian, CHEN Fei-wu. Multiwfn: A Multifunctional Wavefunction Analyzer[J]. Journal of Computational Chemistry, 2012, 33(5): 580-592.

[18] LIU Ze-yu, LU Tian, CHEN Qin-xue. An Sp-Hybridized All-Carboatomic Ring, Cyclo[18]Carbon: Electronic Structure, Electronic Spectrum, and Optical Nonlinearity[J]. Carbon, 2020, 165: 461-467.

[19] YANG Jun-qing, ZHANG Gu-dan, ZHANG Jian-guo, et al. New Perspectives on the Laser Initiation for Metal Tetrazine Complexes: A Theoretical Study[J]. Physical Chemistry Chemical Physics: PCCP, 2021, 24(1): 305-312.

[20] LIU Run-ze, SHANG Fang-jian, XIONG Ying, et al. A Combined Experimental and Theoretical Investigation of the Excited-State Dynamics of 2, 4, 6-Trinitrotoluene (TNT) in DMSO Solvent[J]. Physical Chemistry Chemical Physics: PCCP, 2021, 23(36): 20718-20723.

[21] HILBORN R C. Einstein Coefficients, Cross Sections,Values, Dipole Moments, and all that[J]. American Journal of Physics, 1982, 50(11): 982-986.

[22] 曹霞, 向斌, 張朝陽(yáng). 炸藥分子和晶體結(jié)構(gòu)與其感度的關(guān)系[J]. 含能材料, 2012, 20(5): 643-649.

CAO Xia, XIANG Bin, ZHANG Chao-yang. Review on Relationships between the Molecular and Crystal Structure of Explosives and Their Sensitivities[J]. Chinese Journal of Energetic Materials, 2012, 20(5): 643-649.

[23] LANDERVILLE A C, OLEYNIK I I, WHITE C T, et al. Reactive Molecular Dynamics of Detonating Petn[J]. AIP Conference Proceedings, 2009, 1195(1): 813-816.

[24] WU C, REE F, YOO C S. A Quantum Mechanical Molecular Dynamics Study of Binary Collisions of Pentaerythritol Tetranitrate (PETN): Its Correlation to Shock Sensitivity[J]. Propellants, Explosives, Pyrotechnics, 2004, 29(5): 296-303.

[25] RICE B M, SAHU S, OWENS F J, et al. Density Functional Calculations of Bond Dissociation Energies for NO2Scission in some Nitroaromatic Molecules[J]. Journal of Molecular Structure: Theochem, 2002, 583(1/2/3): 69-72.

Qualitative Effects of Ultraviolet Light on PETN Based on TDDFT

ZHANG Bao-sen1, ZHANG Shu-hai1, GOU Rui-jun1, CHEN Ya-hong1, ZHU Shuang-fei1, MA Kun2

(1. School of Environmental and Safety Engineering, North University of China, Taiyuan, 030051, China; 2. Shaanxi Applied Physics-chemistry Research Institute, Xi'an 710061, China)

The work aims to study the stability degradation mechanism of PETN explosives under the action of ultraviolet light. Based on the TDDFT (time-dependent density functional) theory, 50 excited states of PETN molecules at the pbe1pbe/6-311G** level were calculated and the ultraviolet absorption spectrum was drawn according to the calculation results, and the excitation characteristics of the three excited states of the maximum absorption peak (the total contribution rate of the three excited states of S9, S10, and S11to the maximum absorption peak was 97.31%) were analyzed by the hole-electron method. These three excited states were set as the research objects to study the Mayer and Laplace bond orders of PETN molecules' weak bond after absorbing specific ultraviolet light. Finally, the electron transfer of PETN molecules excited to the excited states were described based on the IFCT (Interfragment Charge Transfer) method. The maximum absorption peak position of the PETN ultraviolet absorption spectrum calculated in this work was 186.6 nm, which was less than 8.4 nm of the absorption peak measured in experiment. Among the three excited states that contributed the most to this absorption peak, the S9state had the highest contribution of 48.27%, and the other two excited states S10and S11were degenerate states with the same contribution of 24.52%. According to the results of the hole-electron analysis, all the three excited states had the characteristics of overall excitation and partial charge transfer. When the PETN molecules absorbed ultraviolet light of a specific wavelength (187.00 nm and 186.92 nm) and then was excited to the corresponding excited state, the Mayer and Laplace bond orders at O—NO2bonds decreased. The IFCT analysis shows that the change of Mayer bond order is guided by the change from n→Pi* on O—NO2. This effect may reduce the stability of PETN molecules.

PETN; ultraviolet; TDDFT; absorption spectrum; molecular stability; excited state

TJ450

A

1672-9242(2022)12-0048-06

10.7643/ issn.1672-9242.2022.12.008

2022–08–16；

2022–11–21

2022-08-16；

2022-11-21

陜西應(yīng)用物理化學(xué)研究所應(yīng)用物理化學(xué)重點(diǎn)實(shí)驗(yàn)室基金（6142602200304）

Science and Technology on Applied Physicial Chemistry Laboratory, Shaanxi Applied Physics-Chemistry Research Institute (6142602200304)

張寶森（1998—），男，碩士研究生，主要研究方向?yàn)楹懿牧舷嚓P(guān)計(jì)算化學(xué)。

ZHANG Bao-sen (1998-)，Male, Postgraduate, Research focus: related computational chemistry of energetic materials.

張樹(shù)海（1969—），男，博士，教授，主要研究方向?yàn)樾滦突鸸に巹┘夹g(shù)及工程爆破。

ZHANG Shu-hai (1969-), Male, Doctor, Professor, Research focus: new initiating explosive technology and engineering blasting.

張寶森,張樹(shù)海,茍瑞君, 等. 基于含時(shí)密度泛函理論的紫外光對(duì)太安的定性影響[J]. 裝備環(huán)境工程, 2022, 19(12): 048-053.

ZHANG Bao-sen, ZHANG Shu-hai, GOU Rui-jun, et al. Qualitative Effects of Ultraviolet Light on PETN Based on TDDFT[J]. Equipment Environmental Engineering, 2022, 19(12): 048-053.

責(zé)任編輯：劉世忠