自從磷光材料被應用到OLED以來，基于磷光機制的有機發光器件研究快速發展，其中代表性的有紅、綠、藍單基色磷光器件和全磷光白色磷光器件。由于自旋禁阻的原因，磷光量子產率要比熒光低得多，因此通常采用重原子效應等方法來提高磷光量子產率。當引入重原子后，自旋軌道耦合作用加強，使禁阻的三重態向單重態的躍遷變為允許，分子停滯在三線態的時間大幅縮短，極大地提高了器件的內量子效率，常用的重金屬原子有Ir、Pt、Re、Os、Cu等。

磷光材料一般應具備良好的光熱穩定性、較大的分子吸光截面、較高的系間竄越能力、室溫下有較高的磷光量子產率和較短的三線態壽命等特征。由于有機材料中磷光的綠光材料最容易獲得，因而研究最早也是從磷光綠光OLED開始，其中Ir(ppy)?是最具代表性的材料，被廣泛地應用到磷光器件中。

磷光發射速率是發光機理研究中的重要參數，本文將以Ir(ppy)?為例，使用BDF和MOMAP軟件計算其磷光發射速率。首先需要通過BDF量化軟件進行結構優化、頻率計算和旋軌耦合計算，之后基于BDF的結構優化和頻率計算結果文件、旋軌耦合計算結果文件使用MOMAP軟件進行磷光輻射速率的計算。

首先使用量子化學計算軟件BDF進行Ir(ppy)?的基態S0和第一激發三重態T1的結構優化和頻率計算。

準備Ir(ppy)?分子結構的xyz文件如下：

61

Ir -2.606160230000 -0.262817540000 0.032585640000

C -3.837298770000 2.407777490000 0.243683290000

C -1.553000180000 2.622608110000 0.521271830000

C -3.991643180000 3.786525490000 0.422094600000

C -4.929719550000 1.476909230000 -0.003856060000

C -1.634158680000 3.988444110000 0.709067180000

H -0.591261780000 2.126459120000 0.554253960000

C -2.889293280000 4.581787990000 0.656041260000

H -4.972006580000 4.231616950000 0.376099570000

C -6.263398690000 1.888616610000 -0.076321100000

H -0.744068550000 4.570116630000 0.889344200000

H -2.999939670000 5.647195600000 0.795286530000

C -5.623129880000 -0.778047400000 -0.424041550000

C -7.264449980000 0.973962040000 -0.319493370000

H -6.527094820000 2.928025340000 0.057085040000

C -6.940359500000 -0.364034500000 -0.495248330000

H -5.397358100000 -1.824987250000 -0.570011260000

H -8.293855190000 1.298147200000 -0.375111620000

H -7.723296650000 -1.084118350000 -0.689505650000

C -2.780095460000 -2.271307610000 -0.073978550000

C -2.962145630000 -2.905938610000 1.179649240000

C -2.704720750000 -3.101798200000 -1.194651040000

C -3.053849400000 -4.297470620000 1.273150410000

C -3.052375550000 -2.037406190000 2.345262930000

C -2.797310890000 -4.477732940000 -1.095358170000

H -2.574810200000 -2.657655530000 -2.171485550000

C -2.970830350000 -5.080470100000 0.142706420000

H -3.190782260000 -4.777448690000 2.231582660000

C -3.251325250000 -2.483356810000 3.656088400000

H -2.735191030000 -5.089028800000 -1.985144810000

H -3.042943740000 -6.155896270000 0.221114210000

C -2.995026490000 0.148733040000 3.092507280000

C -3.319187740000 -1.576319800000 4.692855730000

H -3.353725050000 -3.536644150000 3.859562860000

C -3.186631940000 -0.222801710000 4.408749290000

H -2.888815760000 1.194668600000 2.833837140000

H -3.472767910000 -1.912494740000 5.707724440000

H -3.232610730000 0.522480480000 5.186904590000

N -2.927015180000 -0.717593640000 2.090156960000

C 0.125501920000 -0.309636820000 -1.075330180000

C 0.242318970000 -0.710559210000 1.197687130000

C 1.518232510000 -0.401785790000 -1.166132500000

C -0.762473880000 -0.044853840000 -2.198985150000

C 1.619464440000 -0.811598610000 1.179643800000

H -0.298309850000 -0.828666640000 2.128258480000

C 2.270161400000 -0.653514600000 -0.037593100000

H 2.007628100000 -0.277927520000 -2.118201950000

C -2.150460220000 -0.005742870000 -1.917121780000

C -0.287950380000 0.157363030000 -3.497950870000

H 2.165728000000 -1.009392350000 2.088259610000

H 3.346017750000 -0.726918300000 -0.099292420000

C -3.004998770000 0.237880450000 -2.994864780000

C -1.165335710000 0.397865570000 -4.532504760000

H 0.771521710000 0.127613130000 -3.708268180000

C -2.529020200000 0.436535570000 -4.277665910000

H -4.071982740000 0.267717660000 -2.824596670000

H -0.792761550000 0.552584980000 -5.535050450000

H -3.220629300000 0.622039860000 -5.087916480000

N -0.487689310000 -0.467248190000 0.117089220000

N -2.608631830000 1.850704260000 0.298598520000

C -4.578169480000 0.115167240000 -0.176155400000

打開Device studio，點擊File-new project，命名為phosphorescence.hpf，將Ir(ppy)?_s0.xyz拖入Project中，雙擊Ir(ppy)?_s0.hzw。接下來進行Ir(ppy)?基態S0的結構優化和頻率計算。選中Simulator → BDF → BDF，在界面中設置參數。在Basic Settings界面中的Calculation Type選擇Opt+Freq，方法選擇PBE0泛函，基組在Basis中的All Electron類型中，選擇Def2-SVP。Basic Settings界面中的其它參數使用推薦的默認值，不需要做修改。

在SCF Settings界面中，DFT Integral Grid選擇Coarse，Convergence Threshold選擇Tight。SCF Settings界面中的其它參數使用推薦的默認值，不需要做修改。

在OPT Settings界面中，Convergence Threshold選擇Tight。OPT Settings界面中的其它參數使用推薦的默認值，不需要做修改。

Freq Settings界面中的參數使用推薦的默認值，不需要做修改。之后點擊 Generate files 即可生成對應計算的輸入文件。選中生成的bdf.inp文件，右擊選擇open with，打開該文件，如下所示：

$compass

Title

C33H24IrN3

Geometry

Ir -2.60616023 -0.26281754 0.03258564

C -3.83729877 2.40777749 0.24368329

C -1.55300018 2.62260811 0.52127183

C -3.99164318 3.78652549 0.42209460

C -4.92971955 1.47690923 -0.00385606

C -1.63415868 3.98844411 0.70906718

H -0.59126178 2.12645912 0.55425396

C -2.88929328 4.58178799 0.65604126

H -4.97200658 4.23161695 0.37609957

C -6.26339869 1.88861661 -0.07632110

H -0.74406855 4.57011663 0.88934420

H -2.99993967 5.64719560 0.79528653

C -5.62312988 -0.77804740 -0.42404155

C -7.26444998 0.97396204 -0.31949337

H -6.52709482 2.92802534 0.05708504

C -6.94035950 -0.36403450 -0.49524833

H -5.39735810 -1.82498725 -0.57001126

H -8.29385519 1.29814720 -0.37511162

H -7.72329665 -1.08411835 -0.68950565

C -2.78009546 -2.27130761 -0.07397855

C -2.96214563 -2.90593861 1.17964924

C -2.70472075 -3.10179820 -1.19465104

C -3.05384940 -4.29747062 1.27315041

C -3.05237555 -2.03740619 2.34526293

C -2.79731089 -4.47773294 -1.09535817

H -2.57481020 -2.65765553 -2.17148555

C -2.97083035 -5.08047010 0.14270642

H -3.19078226 -4.77744869 2.23158266

C -3.25132525 -2.48335681 3.65608840

H -2.73519103 -5.08902880 -1.98514481

H -3.04294374 -6.15589627 0.22111421

C -2.99502649 0.14873304 3.09250728

C -3.31918774 -1.57631980 4.69285573

H -3.35372505 -3.53664415 3.85956286

C -3.18663194 -0.22280171 4.40874929

H -2.88881576 1.19466860 2.83383714

H -3.47276791 -1.91249474 5.70772444

H -3.23261073 0.52248048 5.18690459

N -2.92701518 -0.71759364 2.09015696

C 0.12550192 -0.30963682 -1.07533018

C 0.24231897 -0.71055921 1.19768713

C 1.51823251 -0.40178579 -1.16613250

C -0.76247388 -0.04485384 -2.19898515

C 1.61946444 -0.81159861 1.17964380

H -0.29830985 -0.82866664 2.12825848

C 2.27016140 -0.65351460 -0.03759310

H 2.00762810 -0.27792752 -2.11820195

C -2.15046022 -0.00574287 -1.91712178

C -0.28795038 0.15736303 -3.49795087

H 2.16572800 -1.00939235 2.08825961

H 3.34601775 -0.72691830 -0.09929242

C -3.00499877 0.23788045 -2.99486478

C -1.16533571 0.39786557 -4.53250476

H 0.77152171 0.12761313 -3.70826818

C -2.52902020 0.43653557 -4.27766591

H -4.07198274 0.26771766 -2.82459667

H -0.79276155 0.55258498 -5.53505045

H -3.22062930 0.62203986 -5.08791648

N -0.48768931 -0.46724819 0.11708922

N -2.60863183 1.85070426 0.29859852

C -4.57816948 0.11516724 -0.17615540

End Geometry

Basis

Def2-SVP

Skeleton

Group

C(1)

$end

$bdfopt

Solver

1

MaxCycle

366

TolGrad

3.0E-5

TolStep

1.2E-4

IOpt

3

Hess

final

$end

$xuanyuan

Direct

$end

$scf

RKS

Charge

0

SpinMulti

1

DFT

PBE0

D3

Grid

Coarse

ThreshConv

1.0D-9 5.0D-7

MPEC+COSX

Molden

$end

$resp

Geom

$end

選中bdf.inp文件，右擊選擇Run，彈出如下服務器提交的界面：

點擊Run提交作業。任務結束后bdf.out，bdf.out.tmp，bdf.scf.molden三個結果文件會出現在Project中。

選擇bdf.out，右擊show view，在Optimization對話框中，顯示結構已經達到收斂標準。

在Frequency對話框中，檢查頻率，若不存在虛頻證明結構已經優化到極小點。

選擇bdf.out文件，右擊open with containing folder打開所在文件夾，在bdf.out文件中查找‘converged in’，緊接著輸出的‘Molecular Cartesian Coordinates (X,Y,Z) in Angstrom : ’下的結構即為優化好的Ir(ppy)?的S0結構。

編輯：黃飛