公募研究
本研究で採用する分割統治法では、全系をドメイン(小部分)に分割し、各ドメインに対して隣接するドメインとは独立に電子状態を計算する。ドメイン内では正しく波動関数が得られるように、各ドメインの周りにはバッファー領域を設定して電子状態計算を行う。今年度は、非断熱・非平衡ダイナミクスを扱う理論手法と本研究で提案するオーダーN電子状態計算法を結合するために、線形応答時間依存密度汎関数法で電子の励起状態を扱い、Fewest-Switches Surface-Hopping法により励起子の非断熱的な動的過程を追跡する分子動力学(MD)法を開発した。これにより、アモルファス有機半導体におけるsinglet fissionの機構を解明し、励起子ダイナミクスに複数の時間スケールがある理由を説明した。具体的には、128個の5,12-diphenyl tetracene分子(C30H20)(6400原子)を用いて、古典的な力場によるMD計算を行うことにより原子の初期配置を得た。分割統治計算を行う際は、各分子をドメインに選び、最近接の数分子をバッファー領域として選んだ。ドメイン+バッファーの系に対して、バッファーの外側からの分子の影響を取り入れて非断熱第一原理MD計算を行い、(1) 各分子間の励起子の遷移確率、(2) 各分子上で励起子が消滅する確率、(3) singlet fissionにより二つのスピン三重項の励起子が生じる確率を求めた。観測された励起子ダイナミクスを再現するためには1 ns程度の計算が必要であるため、得られた遷移確率等を用いてkinetic Monte Carlo計算を行った。この計算により、分子のバックボーンプレーンが平行に配向した分子間がsinglet fissionのhot spotとなり、1 ps程度での急激な三重項励起子の増加に寄与していることが明らかになった。
24年度が最終年度であるため、記入しない。
すべて 2013 2012
すべて 雑誌論文 (5件) (うち査読あり 5件) 学会発表 (13件) (うち招待講演 1件)
Applied Physics Letters
巻: 100 ページ: 203306 (5 pp.)
10.1063/1.4719206
Journal of Chemical Physics
巻: 136 ページ: 184705 (8 pp.)
10.1063/1.4712616
Physical Review B
巻: 86 ページ: 104202 (9 pp.)
10.1103/PhysRevB.86.104202
ournal of Physics: Conference Series
巻: 402 ページ: 012012 (7 pp.)
doi:10.1088/1742-6596/402/1/012012
巻: 402 ページ: 012044 (9 pp.)
doi:10.1088/1742-6596/402/1/012044