少数粒子系严密计算-KEK.ppt

少数粒子系厳密計算 肥山詠美子（理研） * 特定領域「ストレンジネスを含むクオーク多体系」公募研究 ～H２１年度）連携研究者？？ 世話人から（A０２：初田さんから） １．現在の計算法について説明 ２．今後５年間でどこまで到達できるか、目指すか ３．研究者間のどのような連携が必要か ４．計算資源やアルゴリズム開発に関する希望 QCDによる核力との関連（A01班とのつながり）、 天体現象との関連（A0３班とのつながり）について 適宜入れるように これを２０分で話すことが要求 原子核だけでなく、素粒子?宇宙の方にも分かりやすく?? １．計算法 　「無限小変位ガウス?ローブ法」を大学院の時に提唱 　　 量子力学的３体系、４体系のシュレディンガー方程式を 　　　厳密に（近似的ではなく）解く方法　 ３体問題 ３体問題シュレーディンガー方程式　：　６変数２階偏微分方程式 ４体問題 汎用性が高いので、物理学のさまざまな分野に適用。 　精密に解くことによって初めて解明できる物理がある 　ことを明らかにしてきた。 ? 構成粒子は何でもよい、 　 質量、電荷を問わない。　 少数多体系のシュレーディンガー方程式を精密に解く 「無限小変位ガウス?ローブ法」の利点 ３体問題 ４体問題 ５体問題 現在は、 さらに クーロン３体問題は １０桁の精度　 0 R ? 粒子間に強い相関がある場合 にも精密に適用できる。 強い相関 (核力など） （電子、陽子、中性子、クオーク、?????） V(R） ３体系のシュレーディンガー方程式の一般形　 ６変数の２階偏微分方程式の固有値問題 H: ハミルトニアン 　（ 固有値 E と 固有関数　　 を求める）。 r (x,y,z) R (X,Y,Z) 量子力学の教科書にある一般的な解き方： 未知関数　　　を、基底関数系 ハミルトニアンをこの基底で対角化して、 固有値 　　 と 係数 　　 を求める。 （できるだけ完全形に近い） で展開し、 ４体系：９変数、　 ５体系：１２変数 (非直交系でもよい) Hi n = 行列要素を計算（多変数積分）する。 行列の一般化固有値問題となる。これを解いて 、 Cn Ni n Cn エネルギー固有値のセット ： 1 , 2 , 3 , ????? 波動関数のセット ： 1 , 2 , 3 , ????? を得る。 　≡ ?ΦiΦn dr dR * N i n = Φi | 1　| Φｎ H i n= Φi | H | Φn 　ここで、厳密解を得るために最も重要なことは、 　良い基底関数を用いること 　良い基底関数とは、　例えば、 　１）すべての粒子間の（強い）相関を取り入れられる。 　２）波動関数の遠方の漸近の形をよく記述できる。 　３）行列要素の計算が容易に、 　　できるだけ解析的に行える。 N i n = Φi | 1　| Φｎ H i n= Φi | H | Φn 基底関数 ４体問題?５体問題に適用するために?? r R 1 1 r R r R 2 2 3 3 C=1 C=2 C=3 それぞれの成分関数を、各座標の基底関数の積で展開する： 粒子間の相関を取り入れるために、3通りのヤコビ座標の組み方を使う。 波動関数を、３通りの座標の成分関数の和で表わす。 未知係数 基底関数 九大流 ガウス型基底関数 （１９８８～） 　　　　　　　　　ガウスの広がり　　　　を等比級数にするのが重要point。 　この基底関数を用いる利点は、 　　（１）短距離の「強い相関」も 　　（２）長距離の「遠方 tail 」も、 　同時によく記述できること。 　九大グループは、さまざまな分野の　３体系の課題に適用して、 　大きな成果を挙げていた。 　　しかしながら、　４体系になると、 　　この基底関数では、計算が非常に困難。 　　球面調和関数の角度積分が大変だから。 九大流 ガウス型基底関数 （１９８８～） X 原点の周りの ガウス関数 色々な方向に 無限小シフトさせ、 固有の重みで　重ね合わせることに より、球面調和関数と同じ角度依存性を表す。 左辺：従来の基底関数 右辺：無限小変位ガウス?ローブ基底関数 シフト係数 シフト?ベクトル シフトのパラメタ（C,D）は容易に決められる 無限小変位ガウス?ローブ法　　(Infinitesimally-Shifted Gaussian Lobe) 　３体?４体系の束縛状態関しては、汎用性が高く、 　非常に精度の高い方法が確立した。 　すべての積分は単純なガウス積分なので、手計算は非常に簡単。 　面倒な球面調和関数は全く現れない。 　シフトを無限小にする極限は、行列要