2026-03-13T04:20:55Z https://ir.soken.ac.jp/oai

oai:ir.soken.ac.jp:00000175 2023-06-20T14:59:48Z 2:427:9

Electronic Structure and Nonadiabatic Dynamics in the Photoisomerization of Polyenes Electronic Structure and Nonadiabatic Dynamics in the Photoisomerization of Polyenes 伊藤, 正勝 イトウ, マサカツ ITO, Masakatsu ポリエンの光異性化は、視覚過程とも関係し、生物的にも化学的にも最も重要な反応の1つであり多くの実験的、理論的研究が行われてきた。特に最近のフェムト、ピコ秒時間分解レーザースペクトロスコピーの発展にともない、反応ダイナミックスの詳細な情報が得られるようになってきた。また、大規模MCSCF法、CI法により電子相関を考慮して、異性化経路に沿ったポテンシャルエネルギーを計算することが可能になってきている。しかし、反応に関与するπ結合は複数あり、多次元的な取り扱いが、光異性化グイナミクスを理解するために不可欠となる。本研究では、このダイナミックスを取り扱うために、多次元での電子状態変化のモデルにもとづくハミルトニアンを開発し、半古典動力学シュミレーションにより状態遷移、エネルギー緩和を解析した。 　光異性化反応を理解するために、関与する3つの電子状態1 1Bu、2 1Agと1 1Agの性質を知る必要がある。実験と電子状態計算の結果から、光学許容の1 1Bu状態へ励起されたポリエン分子は、数十フェムト秒という極めて短い時間のうちに、ほぼ平面構造を保ったままで、最低励起状態の2 1Ag状態に内部転換を起こし、この状態で異性化反応を経て基底状態1 1Agに失活するまでの比較的長い時間を過ごすと考えられている。従って、最低励起状態（2 1Ag）、基底状態（1 1Ag）の2つのポテンシャルエネルギー面、その間の非断熱カップリングが異性化のダイナミックスを主にコントロールしていると考えられる。 　彼らは、これらの電子状態を原子価結合（VB）法の描像により記述するモデルハミルトニアンを構築した。まず、π、σ電子を分離し、π電子はVB型のモデルハミルトニアンで、σ結合は分子力学法に似たポテンシャル関数で記述できると仮定した。次に、異性化に伴う2 1Ag、1 1Ag状態の変化は、CASSCF波動関数によって扱うと、イオン的VB構造の重み変化も無視できないが、Heisenberg型ハミルトニアンでのスピン間相互作用がイオン的VB構造からの寄与を繰り込むと考え、中性VB構造のみで記述した。つまり、VB型波動関数のスピンカップリングで電子状態を表現し、スピン間相互作用として表現したπ結合が異性化により変化するため、状態変化が起こるとした。 　隣接スピン間の相互作用のみを考えた単純なHeisenberg型モデルとab-initio計算で、ポリエン分子CnH2n（n=2，3，4）のCASSCF計算によるポテンシャル面の特徴を定性的に再現できる。すなわち、2 1Ag状態において、平面構造はC4H6、C6H8は不安定であるのに対し、C8H10は安定である。2つのポテンシャル面はCC結合が2つまでのCC結合がどんなやり方でねじれても非断熱遷移が起こるのに充分なほど近づかないが、特定の組合せの3つのCC結合が90℃までねじれると、2 1Ag、1 1Ag状態は同じ規約表現に属するにもかかわらず縮遇する。その結果、蛍光寿命C4H6（＜測定限界）　C6H8（～1ps）で短く、C8H10（100ns）で長い事を説明している。 　特に、s-trans butadieneについては、種々の分子配座に対してCASSCF/DZ+d計算（4軌道、4電子）を行い、非線形最小自乗法によりHeisenberg型モデルのパラメーターを決定したところ、様々な異性化経路に沿ったポテンシャルエネルギー面の起伏、非断熱カップリングの変化が、このHeisenberg型モデルで充分に定量的な精度で再現できることが明らかになった。またこれらポテンシャルエネルギー面、非断熱カップリングの特徴は、Heisenberg型2x2ハミルトニアンを介して、π電子のスピン間相互作用の異性化経路に沿った変化で容易に説明する事ができる。 　このHeisenberg型モデルを電子状態の記述に用い、半古典動力学シュミレーションにより非断熱遷移を含むトラジェクトリ一計算を行い、異性化の速度、非断熱遷移の確率、異性化に伴う分子内エネルギー緩和の機構を解析した。このシュミレーションでは、電子状態は量子論的に、分子配座は古典座標として扱い、系の時間発展はsurface hopping、断熱状態のスイッチと、その時の速度補正はTullyのアルゴリズムを用いた。 　トラジェクトリーの解析から、2 1Agから1 1Agへの非断熱遷移は3つのCC結合がねじれてそのエネルギー差が小さくなり、さらにCCC結合角の大きな減少することにより引き起こされていることが示された。更に、swarm dynamicsを解析した結果、2 1Ag→1 1Agの内部転換は、（1）励起直後の0.1ps： 励起直後のπ電子エネルギーの急激な緩和、（2）0.1ps～0.7ps： 二重結合まわりのねじれが緩和し、一重結合がねじれ始めるまでの誘導期間、（3）0.7ps～： 2 1Ag状態が減衰し始める、といった、3つのステージに分かれることも示唆された。基底状態に遷移してからは、π電子エネルギーが多くのaccepting modeに分散し、promoting modeであるCC結合ねじれ、CC伸縮、CCC変角のモードに充分に集中しないため状態遷移は不可逆となる。 　今までに、モデルハミルトニアンの構築と非断熱動力学シュミレーションの組合せにより、butadieneの内部転換メカニズムをあきらかにすることができた。一方で、もつと長いポリエンについても、モデルハミルトニアンによる定性的記述が可能となったので、同様の手法の組合せにより光異性化グイナミクスの一般的な描像に近づくことができると考えている。 application/pdf 総研大甲第250号 eng thesis AM https://ir.soken.ac.jp/records/175 博士（理学） 1997-03-24 総合研究大学院大学 https://ir.soken.ac.jp/record/175/files/甲250_要旨.pdf application/pdf 306.0 kB 2016-02-17 https://ir.soken.ac.jp/record/175/files/甲250_本文.pdf application/pdf 9.4 MB 2016-02-17