研究実績の概要 |
光合成システム(PSII)における太陽光を用いる水分解反応は、Mn錯体と高度にカップリングしており、それを触媒するのが「CaMn4O5クラスター」である。これは強相関電子系に属しており、その電子物性・機能発現・反応性の解析には高精度量子化学計算が必須である。しかし、分子軌道法での高精度計算の王道と言われるCC法やMR-CI法は、その適用範囲が比較的少数原子系に限定されており、この巨大分子系に適用することは、計算機資源や計算時間などの制約から極めて困難である。その理由で、電子相関を交換相関汎関数に繰り込んだDFT法、特にHybrid-DFT法が、現在までに数多く行なわれた。しかし、手法内パラメータを順次変更すると、中間体の相対安定性が大幅に変動するため、水分解反応機構の解明には予見性を備えたアプローチが必須であり、その限界をこの研究を遂行することで突破する。
太陽光を利用する水分解反応を触媒するCaMn4O4XYZWクラスターの分子構造は、実験により解明されており、ここでのX, Y, Z, Wはクラスターに配位している水(H2O)あるいはその酸化誘導体(OH-, O2-)などである。しかし、実験では水素原子の位置は見えないため、理論計算で詳細な物性データを得る。
本年度では、特にKokサイクル(S0, S1, S2, S3, S4)の最高酸化状態であるS3状態の解析を行った。UB3LYP法(およびHF成分を15%,10%に変化させたHybrid-DFT)を構造最適化に用いて、複数の安定構造を決定した。加えて、高精度計算としてDLPNO-CCSD(T)法を実行し、相対安定性を厳密に求めた。一方、DMRG-CAS, QMC-CAS法をMn二核・Mn三核錯体に適用して、その磁気的相互作用の解析を進めた。
|
現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
現在までのUB3LYP法とDLPNO-CCSD(T)法を用いたCaMn4O5クラスターの解析を進めてきた。このKokサイクル(S1→S2→S3→(S4)→S0)に関する理論的解析では、基本S1状態(暗所で安定)の前後にあるS0・S2状態と、酸化数が最高度に進んだS3状態に関して、初年度と昨年度に学術論文にまとめることができた。それらの結果、UB3LYP法による構造最適化計算にて、L-open, R-open(およびその派生)という複数の安定構造が発見された。加えて、各酸素原子(O)の状態がO2-,OH-, H2Oであるかを理論的に考察すると(実験測定で水素原子(H)は観測できないから)、さらにこれらの安定構造は複雑になることが判明した。特にS3状態の分子構造は酸素発生の最終段階であり、続く化学反応経路を解析するために重要であると期待された。
計算手法論的には、UB3LYP法(およびHFとDFTの混合割合を15%,10%に変化させたHybrid-DFT法)を上記の構造最適化計算に活用してきた。ここでは、(2)蛋白場を連続体モデルで考慮した計算を実行した。さらに、これら判明した安定構造の全てにおいてDLPNO-CCSD(T)計算を実行した。この計算実行には非常にコストがかかり、高性能計算機(80core, 6TBメモリ)を2週間から1ヶ月間連続運転することで可能となった。その結果各安定構造間の相対安定性を詳細に解析できた。
|