量子化学における超大规模半正定値计画问题と并列计算…解析.ppt

量子化学における超大規模半正定値計画問題と並列計算による高速求解 神奈川大学工学部情報システム創成学科 山下真 Makoto.Yamashita@is.kanagawa-u.ac.jp 東京工業大学 Mituhiro Fukuda 東京大学 　　　　　　 中田真秀 Outline 量子化学の基本的で重要な問題に 基底状態のエネルギー(Ground State Energy) そのときの電子構造 を求めることがある この問題は、半正定値計画問題(SemiDefinite Programming, SDP)に緩和できる 並列計算を用いて実用時間内に解くことができる 量子化学＋数理最適化＋並列計算の３つの分野 目次 量子論とSDP 基底状態におけるエネルギーとは SDPへの帰着 並列計算による高速化 数値実験 量子論とSDP(その１：座標回転） まずは、大まかな話から 波動関数と物理量のオペレータの行列表現をそれぞれ　　　　とすると、物理量は　　　　　　　の２次形式で計算できる 適切な回転行列　　を用いると対角行列　　で 本質的な量（固有値）と座標回転を同時に考えることが必要 量子論とSDP（その2：対称行列） 本質的な量（固有値）と座標の回転を同時に考えることが必要 対称行列で変数の表現 電子構造を対称行列で表現すると、（直観的には）固有値が電子の個数となる 変数が半正定値行列に制約される 半正定値行列を変数に持つ線形計画問題が半正定値計画問題 半正定値計画問題（SemiDefinite Programming, SDP) 基底状態におけるエネルギーとは 直観的には、原子?分子系の最も安定した状態が基底状態 化学反応で発生するエネルギーなどを計算するときに必要となるエネルギー これを求めることは、量子化学における基本的なこと 量子力学における定式化 Heisenberg の行列力学 Schr?dinger の波動方程式 これらは、Hilbert 空間の内積が一致するため、物理学的には等価 時間に依存しないSchr?dinger方程式 . 最小固有値がＧＳＥ 水素原子のときには 微分方程式を解かずに最小固有値のみを求める 行列表現 水素原子の場合、全空間での電子の数は 1 であるため、波動関数は Hilbert 空間 　　　に属している Hilbert 空間での基底をとすると、波動関数を展開できる このとき波動関数　　とベクトル　　を同一視できる 行列の最小固有値 Schr?dinger 方程式を行列の固有方程式にただし、 したがって、(Ritz の変分原理ともよばれている） 多電子ハミルトニアン ハミルトニアンに電子-電子間などの項が必要 原子核の位置を固定して近似 (Born-Oppenheimer approximation) . 多電子波動関数の基底 多電子波動関数はPauliの排他原理を満たす必要がある 一電子波動関数の基底をとすると、 電子波動関数の基底全体はである。ただし、 既存手法 一般の多電子波動関数はと展開される GSE計算の既存手法 Hartree-Fock Full CI CCSD,CCSD(T), etc Hartree-Fock Self Consistent Field 法で得られた一電子波動関数 エネルギーの低い順に 基底状態の波動関数を　　　　　　　　　　　　　　　　とする 使用しているのは、のひとつの項のみ GSEと誤差１％程度とされている Full CI 一次結合に現れるすべての項を用いて最小固有値計算 与えられた基底については、厳密解を与える 計算量が膨大　　　　　　　　　の行列の固有値分解が必要 大きな分子に適用できない CCSD(T) Coupled Cluster expansion using Single and Double excitations with perturbation treatment of Triples 一次結合　　　　　　　　　　　　　　の中で効果の高い項のみで計算 量子化学の分野で信頼されている手法 これと同等の精度を出すことがひとつの目安 SDPによるGSE計算 Full CI 同様に最小固有値の計算をベース Full CI では、von Neumann 密度行列と同等 必要な情報量は、２次の縮約密度行列のみ ２次の縮約密度行列を変数とすることで計算量を抑える SDPによるGSE計算の歴史（１） 1960年代：Coleman, Garrodら 半無限線形計画問題へ