RMD)シミュレーション実施のためのプログラム開発を実施した。学(九州大学薬学部)らにより,今年度新規に合成されたアクリジニウムエステル誘導体であるDCMACも含め,フェニル基における置換基の種類の差異および位置により発光強度が大きく変化することが報告されている。誘導体の種類により燐光の強度が変わり,定量的理解には至っていない。そこで量子化学計算を用い,遷移状態および安定構造の探索などを実施した。溶媒効果を暗黙分極連続媒体モデル法により静的に考慮した。これまでの成果と本結果より一つの仮説が見出された。反応経路は,大まかに遷移状態であるジオキセタン構造が現れた後,発光種である三重項状態である励起アクリドンが生成され,この励起種が強い燐光を示す。従って,ジオセキタン構造が現れる前後の反応過程が,燐光の強度を決定している。置換基により明らかに二つの異なる経路に分類されることが判明した。この結果でも定量的発光強度の再現は出来なかった。RMDシミュレーションのプログラム開発と予備的計算を実施した。約1ピコ秒時間発展させた古典軌道の結果には,ジオキセタン構造が現れなかった。計算は安定し,系のポテンシャルエネルギーが下がっていることが判明した。従って,脱プロトン化が起こり,ジオキセタン構造を取るものと思われる。planned.Reason基本的には予定通り進んでいる。RMD)シミュレーションプログラムの開発および予備的計算の実施において変更が必要となってきている。主な理由は,以下の通りである。古典軌道の時間発展において毎時間ステップに実施する量子化学計算により,on-the-flyでポテンシャルエネルギーとその傾きを求め,分子動力学シミュレーションを実施する。この量子化学計算において並列計算効率が圧倒的に高い,石村和也(分子科学研究所)が開発したSMASHプログラムをこれまで利用して,開発および予備的計算を実施してきた。本研究で扱う分子のジオセキタン構造が現れる前後の反応過程までは,一重項状態の基底状態で計算を進めるため,1ピコ秒の時間発展も約12日間程度で達成できた。因みに時間ステップ(Δt)は0.05フェムト秒であることを考えると約2000ステップで1ピコ秒に達する。2.40GHzを用い,28スレッド並列計算である。RMDシミュレーションプログラムの全計算コストの9割が,量子化学計算の部分であることから,明らかに素晴らしい性能を示していた。16を利用し,テストを行っている。現在,量子化学計算プログラムの再検討を必要としている。Activity現在までの進捗状況にて述べた通り,量子化学計算における並列化効率が,本研究の推進に大きく依存している。16以外のプログラムパッケージで性能を評価し,最適な解を探し出す。溶媒にジメチルスルホキシド(DMSO)を用い,溶媒を露わに考慮した量子古典混合(QM/MM)モデルを用い,計算を進める。RMD)シミュレーションを実施し,自由エネルギー評価のもとに定量的議論を行う予定である。)により生物発光に匹敵する発光をしめす。本年度は,この反応メカニズム解明に着手する。ヒトアルカリフォスファターゼ(1EW2)のPDBデータより,分子動力学シミュレーションを実施し,タンパク質の構造をコンピュータ上に再現する。水系の溶媒中では発光が激減することが実験では観測されていることから,水中でのシミュレーションを実施する。RMDシミュレーションの予備的計算を進める。