在使用量化軟件BDF完成Azulene的基態(tài)和激發(fā)態(tài)的結(jié)構(gòu)優(yōu)化、頻率計算以及非絕熱耦合的計算，并計算完MOMAP輸入文件中需要的參數(shù)后，接下來將使用MOMAP軟件對azulene的S1→S0的輻射速率和內(nèi)轉(zhuǎn)換速率進(jìn)行計算，通過對比輻射速率和內(nèi)轉(zhuǎn)換速率來解釋azulene的S1→S0的熒光難以被觀測到的原因。 首先計算S1→S0的熒光輻射速率，第一步為做電子振動耦合（electron-vibration coupling, EVC）計算，該計算基于量化計算輸出的分子振動頻率、力常數(shù)矩陣，同時在內(nèi)坐標(biāo)以及直角坐標(biāo)系下，計算分子躍遷發(fā)生初末態(tài)間的模式位移、黃昆因子、重整能以及Duschinsky轉(zhuǎn)動矩陣。 將bdf的S0優(yōu)化頻率計算結(jié)果改為azulene-s0.out，將S1的計算結(jié)果改為azulene-s1.out，同時放在EVC計算文件夾中。

EVC計算的輸入文件momap.inp為：

do_evc = 1

&evc

ffreq(1) = "azulene-s0.out"

ffreq(2) = "azulene-s1.out"

proj_reorg = .t.

/

提交腳本文件momap.slurm，任務(wù)運行結(jié)束后，生成如下文件

其中 evc.cart.dat 為利用直角坐標(biāo)系計算得到的模式位移、黃昆因子、重整能以及 Duschinsky 轉(zhuǎn)動矩陣的結(jié)果；而evc.dint.dat 為內(nèi)坐標(biāo)計算模式位移、黃昆因子、重整能，直角坐標(biāo)計算Duschinsky轉(zhuǎn)動矩陣的結(jié)果。

打開evc.cart.dat文件，查看總重整能的數(shù)值：

并與evc.dint.dat文件中的數(shù)值進(jìn)行比較：

比較 evc.cart.dat 以及 evc.dint.dat 文件內(nèi)的重整能，若重整能相差不大（< 1000 cm-1），使用evc.cart.dat文件進(jìn)行后續(xù)計算，若重整能相差較大（一般情況evc.cart.dat大于evc.dint.dat），使用evc.dint.dat文件進(jìn)行后續(xù)計算。這里兩者較為接近，且數(shù)值較為合理（<10000 cm-1），可使用evc.cart.dat進(jìn)行下一步S1→S0的熒光輻射速率計算。 此外，還可以根據(jù)EVC計算的結(jié)果文件做更多的后處理。 evc.dx.x.com和evc.dx.x.xyz為兩個電子態(tài)分子疊加圖，其中evc.dx.x.com可用GaussView打開，在View-Display Format-Molecule中選擇Tube類型，顯示如下：

evc.dx.x.xyz可用Jmol軟件打開，顯示如下：

evc.dx.v.xyz也是分子態(tài)疊加圖，用Jmol軟件打開，選擇顯示-矢量-0.1A，顯示如下：

evc.cart.abs為Duschinsky矩陣文件，可用來畫Duschinsky矩陣二維圖。可以在Device Studio中，選擇Simulator-momap-analysis，打開evc.cart.abs文件，顯示如下：

同樣的方法打開evc.dint.abs文件，可在ColorType下拉選框選擇顯示顏色，顯示如下：

Duschinsky矩陣都是用直角坐標(biāo)計算的，二者選一即可。 將evc.vib1.xyz和evc.vib2.xyz文件與evc.cart.dat文件放在同一路徑下，在Device Studio中，選擇Simulator-momap-analysis，打開evc.cart.dat文件，出現(xiàn)基態(tài)和激發(fā)態(tài)下各個振動模式對應(yīng)的重整能和黃昆因子圖，顯示如下：

點擊Choose Color可以改變線的顏色，改變Set Width中的數(shù)值可改變線的粗細(xì)。點擊圖中的線，右側(cè)結(jié)構(gòu)將展示相應(yīng)的振動模式。振動的快慢可通過Animation Frequency調(diào)整，振動的幅度可根據(jù)Displacement Amplitude顯示。