「バリオン間相互作用に基づく核物質の構造」 2009年6月25－27日盛岡つなぎ温泉ホテル紫苑密度汎関数理論の基礎と原子核への応用筑波大学数理物質科学研究科計算科学研究センター矢花一浩自然界を構成するミクロな物質の世界原子核原子・分子・固体・液体強い相互作用で束縛した陽子と中性子の多体系クーロン力で束縛した電子とイオンの多体系基底状態を記述する理論の分類 Few-body approach Green function Monte Carlo (GFMC) Configuration mixing (CI) Hartree-Fock theory MBPT Many-body perturbation theory 生の相互作用か、有効相互作用か DFT （Density functional theory）密度汎関数理論 Hartree-Fock理論と密度汎関数理論の違い解いているのは似たような方程式なのだけれど・・・       r1, r2 ,, rN   1 r1 2 r2 N rN  （反対称化＝Slater行列式）   p2 i  v r i   ii 1   2  ... 2m   2  r    i r  i Hartree-Fock理論密度汎関数理論エネルギー期待値を最小にする１粒子軌道を探す。原理的には密度の変分により基底状態を厳密に求めることができる。（エネルギー汎関数が存在定理。）  H 0 i   ・無限系に対しても近似解となる。・変分原理を基礎に持つ。・摂動の高次項を計算すれば、系統的に解の精度を改善できる。  min  H   min  H   min    r    r   E r  ・「厳密な」エネルギー密度を求める手続きは与えられていない。・系統的に解の精度を上げる手段がない。 Kohn-Shamの理論 Kohn-Sham (1965) 原理的に基底状態を、ある自己無撞着な平均場方程式を解くことにより求めることが可能     E r   TS  r   E r   TS  r   V r   i   p2 i i  V  r  2m 2体相互作用のない仮想的な系の  エネルギー汎関数 TS  r  を導入   2     r    i r  , i     i  j   ij  一様無限系のエネルギー密度を出発点にとる（局所密度近似）・無限系は定義からexact (Hartree-Fockでは近似）。・運動エネルギーの非局所性はKohn-Sham理論で考慮。・局所密度近似のもとでは、ポテンシャルはlocal (NonlocalなFockポテンシャルがないので、解くのが楽）物質科学の世界で、密度汎関数理論が受け入れられた経緯物理の世界（物性理論）で密度汎関数理論が受け入れられた訳（1965頃）（某年長の物性理論の人から聞いた話） Slaterが原子のHartree-Fock計算を簡単化するため、交換ポテンシャルに対して局所密度近似を行っていたが、0.7倍程度に弱めることが必要だった。（Xα法）密度汎関数理論で、エネルギー汎関数を局所密度近似してから変分すると、ポテンシャルを局所密度近似する場合に比べて(2/3)の因子が現れた。化学の世界（量子化学）で密度汎関数理論が受け入れられた訳（1990頃）（某中堅の物性理論研究者から聞いた話）当初、局所密度近似による計算は、「chemical accuracy」が無かった。勾配補正を取り入れることにより、化学分野で意味のある答えになった。（W. Kohnはノーベル化学賞を受賞）物質科学のエネルギー密度（電子ガス）・相互作用：クーロン力（２体力のみ、３体力なし）   ・一様物質のエネルギー密度 E  無限物質計算（電子ガス）：済み（1980）・「局所密度近似」は1980年頃に完成・1990年頃に、勾配補正。量子化学研究者が参画。・精度向上の探究が続く最も成功しているのは、「ハイブリッド」交換項に対する局所密度近似と非局所Fockポテンシャルをブレンド（７：３）B3YLP 原子核のSkyrme Hartree-Fock計算は、「平均場理論」か、「密度汎関数理論」か？ Skyrme Hartree-Fock法では、 Skyrme力（２体力＋３体力）に対するSlater行列式の期待値 (3) E   vij (2)  vijk  ij ijk を取るが、・生の核力とSkyrme力は直接つながらず、Skyrme力のパラメータは、エネルギー汎関数を表現するパラメータと見ることができる。・実際に使われているエネルギー汎関数は、期待値の形に表せない項を含む。例えばρα（αは非整数）密度汎関数理論と考えた方が良いだろう。 Skyrme-Hatree-Fock法で用いるエネルギー汎関数   E   dr H r     q r    i r , , q 2 i ,   2   q r    i r , , q i ,       *  J q r    i  i r , , q i r , ' , q     ' i , , '      2  1  1  2   1 1 1 H r    r   t0 1  x0     x0   n2   2p   t1  t2   t2  t1 n n   p p  2m 2  2  2 8   4   1 1 1 1  t2  3t1   3t1  t2 n n   p  p  t2  t1  J n2  J p2  t3 n  p   HC 16 32 16 4       1  W0   J  n  J n   p  J p 2    スピン軌道力、∇ρを含む項は、核物質の性質だけでは決まらない。  原子核のDFT計算の目標全核図表にわたる記述と予言（エネルギー、密度、形、・・・）信頼できるエネルギー汎関数の構築 ⇔核物質の状態方程式（ρ0付近） r-process核（不安定核）実験・弱いビーム強度、短い実験時間、限られた核種計算の強み・安定核と不安定核は（計算の観点からは）差がない・計算機能力の増大様々な理論による原子核の質量予測 Model r.m.s. error # of nuclei # of parameters Bethe-Weizsäcker formula 2.97 MeV 1768 5 Lunney et al., RMP 75 (2003) 1021 Finite-Range Droplet Model 0.609 MeV 1654 30 Möller, Mix et al., At. Data Nucl. Data Tables 59, 185(1995) KUTY 0.68 MeV 2921 34 Koura, Uno, Tachibana, Yamada, NPA 674 (2000) 47 DZ 0.375 MeV 1571 28 Duflo, Zucker, PRC52 (1995) R23 Extended Thomas-Fermi + Strutinsky integral 0.709 MeV 1719 Hartree-Fock + BCS 2 - 3 MeV 1029 2.6 MeV 1315 HFB-8 0.635 MeV 2149 About 15 Samyn and Goriely et al., PRC70, 044309 (2004) HFB-15 0.678 MeV 2149 About 15 Gorirly and Samyn et al., PRC77, 031301 (2008) Relativistic Mean-Field References Goriely, AIP Conf. Prec. No. 529 (2000), p. 287 About 10 Tajima et al., NPA603 (1996) 23 Lalazissis et al., At. Data Nucl. Data Tables 71, 1 (1999) Bethe-Weizsaker mass formulaからのずれ（平均偏差2.97MeV) Skyrme-HFB mass formulaからのずれ（平均0.729MeV) 密度汎関数理論の拡張による励起状態の記述・時間依存密度汎関数理論・配位混合、射影、相関エネルギー励起状態の記述：時間依存密度汎関数理論（TDDFT=Time Dependent Density Functional Theory) バネで結ばれた質点（のつもり）基底状態＝つりあいの位置  V x1, x2 ,, xN   0 xi Hartree-Fock or Kohn-Sham   p2 i  v r i   ii 2m   2  r    i r  i 励起状態＝基準振動モード  mi xi   V x1, x2 ,, xN  xi xi t   Ai eit RPA or TDDFT   p2   i t   v r , t  i t   i  i t  2m t   2  r , t     i r , t  i 時間依存密度汎関数理論＋線形応答＝RPA 2  p   i t   v r , t  i t   i  i t  2m t   2  r , t     i r , t  i 原子核の応答関数と核物質の性質例：圧縮率とGMR（モノポール振動）元素合成との関わり (n,γ）反応（光応答）分子の光吸収スペクトル（電子の応答） N2 H2O 川下、中務、矢花 O2 C6H6 原子核の光応答計算 (Continuum RPA) 稲倉、中務、矢花 Expt. vs. Calc. 原子核の応答関数を計算してみると、原子分子（電子応答）に比べてまだまだ粗が目立つ。分子 C6H6 ・軽い核、高エネルギー側で顕著な差エネルギー汎関数の問題？テンソル相関の重要性？いくつかの球形核の巨大双極共鳴（光吸収）密度汎関数理論殻模型 AMD 少数系厳密計算 GFMC 原子核の大きさにより違う理論が用いられているようにみえるけど・・・核力の精緻化：有効核力から現実的核力、さらにQCD核力へ配位空間の拡大：模型から非経験へ殻模型球形ポテンシャル中のいくつかの軌道の中で完全対角化を目指す。 ⇒軌道を増やす。平均場（密度汎関数）変形したポテンシャル中の一粒子軌道からなる単一のSlater行列式 ⇒重ね合わせる「模型空間」を広げていくとどの模型もAb-initioになる。（有効核力が難問）クラスター模型いくつかのクラスター自由度で状態を記述する。 ⇒多様な分解を取り入れる微視的クラスター計算（延与、京大基研） Excitation energy (MeV) GFMCによる原子核の第一原理計算（Argonne National Lab.の理論グループ） 10 2+2 0+3,4 ? α粒子の直鎖分子構造？３ 7.65 0+2 α粒子のボーズ凝縮状態 0 0+1 ground states 基底状態重元素生成の鍵になる３α→12C反応が経由する共鳴状態虚時間法による基底状態探索の途中に様々な構造が現れる。任意のSlater行列式 16O Local minimum Hartree-Fock 16O原子核：計算で得られる波動関数（Slater行列式）様々なLocal minimum や shoulder が現れる • 決まったパリティ・角運動量状態を抜き出す（射影計算） • 波動関数を重ね合わせる（配位混合）軌道関数の表現：３次元格子303程度 Input parameter: N=6, Z=6、相互作用は現象論（Skyrme力） Energy spectrum in 12C SGII second 0+: 42.4% ground 0+: 94.8% 3-: 78.7% third 0+: 62.1% Energy spectrum in 20Ne SGII (0p1/2)-1(sd)5 Skyrme力で配位混合計算をするときに困ること・密度汎関数理論で配位混合をしてもいい？・ゼロレンジ力で本当の基底状態はない？・３体力やρα依存力の行列要素はどうしたら良い？まとめ・密度汎関数理論（DFT）は、一様無限系のエネルギー密度を出発点として有限・非一様物質の性質を探るアプローチである。・CIやMBPTに比較すると、系統的に精度を上げる手段に欠ける。・精度と労力のバランスが良く、物質科学・原子核ともに非常に良く用いられている。・密度汎関数法の対象は、基底状態に加え、励起状態も。・時間依存密度汎関数法＋小振幅近似＝RPA 巨大共鳴⇔核物質情報、エネルギー汎関数のアセスメント、光応答など元素合成にも重要な情報をもたらす。・不安定核構造研究の醍醐味の一つは、 N,Z,E３次元で眺めた構造の多様性。模型研究⇒非経験計算へ DFT＋CI+Projection,、N,Z,E３次元空間での核構造予測