Lattice QCD approach to the nuclear force

Slides:



Advertisements
Similar presentations
QCD Sum rule による中性子電気双極子 モーメントの再評価 永田 夏海(名古屋大学) 2012 年 3 月 27 日 日本物理学会第 67 回年次大会 共同研究者:久野純治,李廷容,清水康弘 関西学院大学.
Advertisements

相対論的場の理論における 散逸モードの微視的同定 斎藤陽平( KEK ) 共同研究者:藤井宏次、板倉数記、森松治.
Generalized Form Factors of the Nucleon in the Chiral Quark Soliton Model カイラルクォークソリトン模型に基づく 核子の一般化形状 大阪大学 原子核理論研究室 D 1 中小路 義彦.
第 5 章 2 次元モデル Chapter 5 2-dimensional model. Contents 1.2 次元モデル 2-dimensional model 2. 弱形式 Weak form 3.FEM 近似 FEM approximation 4. まとめ Summary.
今まで行ってきた研究 2003年11月 九州工業大学 工学部 物理  鎌田 裕之  .
HBT干渉法における 平均場の効果の準古典理論
Direct Search of Dark Matter
ニュートリノ干渉・回折 飛田 豊 (北海道大学) Collaborators 石川 健三、千徳 仁 (北海道大学)
格子QCDによる核力ポテンシャル 石井理修 (東京大学) 2007年3月28日日本物理学会(都立大)45分講演、日本語.
---核力プロジェクトの近未来の拡張の形(私見)---
木村 匡志 極限ブラックホール近傍の 高速粒子衝突における “バックリアクション“の影響について (YITP 元OCU)
山崎祐司(神戸大) 粒子の物質中でのふるまい.
クラスター変分法による 超新星爆発用 核物質状態方程式の作成
BRST対称性とカイラル対称性の破れの 格子QCDシミュレーションによる検証 scirqcd 帝京大学理工学部 古井貞隆
Pion mass difference from vacuum polarization in lattice QCD E. Shintani, H. Fukaya, S. Hashimoto, J. Noaki, T. Onogi, N. Yamada (for JLQCD Collaboration)
Shell model study of p-shell X hypernuclei (12XBe)
中性子過剰核での N = 8 魔法数の破れと一粒子描像
埼玉大学大学院理工学研究科 物理機能系専攻 物理学コース 06MP111 吉竹 利織
to Scattering of Unstable Nuclei
原子核物理学 第4講 原子核の液滴模型.
RHIC-PHENIX実験での 直接光子測定
クォーク模型による バリオン間相互作用 鈴鹿高専 仲本朝基 1.Introduction 2.クォーク模型の特徴
現実の有限密度QCDの定性的な振る舞いに
中性子星の超流動に対する ハイペロン混在の効果
J-PARCにおけるチャーム原子核生成 須藤 和敬 (二松学舎大学) ストレンジネス核物理2010, KEK
The Nuclear force from lattice QCD
原子核物理学 第8講 核力.
反核子のオフシェルエネルギーでの振舞いおよび ガモフテーラー和則における中間子生成強度
電子 e 光子 g 電磁相互 作用を媒介 陽子 中性子 中間子 p n ハドロン 核力を  媒介 物質の 究極構造 原子 原子核 基本粒子
Λハイパー核の少数系における荷電対称性の破れ
ストレンジネスが拓く エキゾチックな原子核の世界
K凝縮ハイパー核の構造と安定性 1. Introduction 2. Formulation 3. Numerical results
非局所クォーク模型Gaussianバリオン間相互作用とその応用
The Effect of Dirac Sea in the chiral model
Muonic atom and anti-nucleonic atom
全国粒子物理会 桂林 2019/1/14 Implications of the scalar meson structure from B SP decays within PQCD approach Yuelong Shen IHEP, CAS In collaboration with.
クォーク模型バリオン間相互作用 fss2 による 低エネルギー nd 弾性散乱
格子QCDからハドロン間相互作用へ 石井理修(筑波大) 連携の位置関係 現在の計算法 今後5年間での到達点・目標 研究者間での連携
First measurement of interference fragmentation function on longitudinally polarized deuteron target at HERMES 小林 知洋、 Gunar Schnell、 大須賀 弘、 田中 秀和、 長谷川.
Azimuthal distribution (方位角分布)
原子核物理学 第2講 原子核の電荷密度分布.
Y. Fujiwara, Y. Suzuki and C. N., to be published in PPNP;
Charmonium Production in Pb-Pb Interactions at 158 GeV/c per Nucleon
4体離散化チャネル結合法 による6He分解反応解析
チャネル結合AMDによる sd殻Ξハイパー核の研究
「核力についてどんなことができるか」展望
原子核の殻構造の相対論的記述 n n n σ ω n σ ω n 柴田研究室 石倉 徹也 1.Introduction n n
Effective π-N Lagrangian for 0νββ decay in R-parity violating SUSY
G-matrix folding modelによる Ξ核の構造・生成
ストレンジネスで探る原子核 -ハイパー核の世界-
中性子過剰F同位体における αクラスター相関と N=20魔法数の破れ
Perturbative Vacua from IIB Matrix Model
QCDの有効理論とハドロン物理 原田正康 (名古屋大学) at 東京大学 (2006年11月6日~8日)
 中小路義彦 (大阪大学) PHYCAL REVIEW D 77, 074011(2008) M.WAKAMATSU & Y.N
? 格子QCDシミュレーションによる南部-ゴールドストン粒子の 質量生成機構の研究 質量の起源 ドメインウォールフェルミオン作用
格子ゲージ理論によるダークマターの研究 ダークマター(DM)とは ダークマターの正体を探れ!
計画研究代表者 京都大学基礎物理学研究所 大野木 哲也
Marco Ruggieri 、大西 明(京大基研)
非等方格子上での クォーク作用の非摂動繰り込み
原子核物理学 第7講 殻模型.
格子QCDの理論的進展と フレーバー物理への応用II
課題研究 P4 原子核とハドロンの物理 (理論)延與 佳子 原子核理論研究室 5号館514号室(x3857)
ハイパー核物理分野から見た K原子核物理へのコメント
Brueckner-AMDの軽い原子核への適用
低エネルギー3核子分裂反応について 法政大学 石川壮一 1.はじめに 2.3体クーロン問題の定式化 p-p-n系
媒質中でのカイラル摂動論を用いた カイラル凝縮の解析
(K-,K+)反応によるΞハイパー核の生成スペクトル
現実的核力を用いた4Heの励起と電弱遷移強度分布の解析
実数および純虚数化学ポテンシャル領域における 2+1フレーバーPNJL模型を用いた QCD相構造の研究
軽い原子核の3粒子状態 N = 11 核 一粒子エネルギー と モノポール a大阪電気通信大学 b東京工業大学
Presentation transcript:

Lattice QCD approach to the nuclear force 石井理修 (筑波大) 青木慎也 (筑波大) 初田哲男 (東京大) Plan of the talk: Introduction and Backgroud Formalism Lattice QCD results Summary For detail, see nucl-th/0611096.

Introduction 核力は核物理においてもっとも重要な基本概念の一つである。 中間距離から長距離にかけての引力は核子を原子核に束縛するため重要。 近距離に存在する斥力芯は、heavy nucleiや核物質が安定に存在するために重要。 核力は原子核物理以外でも重要な役割を果たす。 中性子星の上限質量や超新星爆発の物理においても重要な影響力をもつ。 湯川秀樹の最初の論文以降 H. Yukawa, Proc. Math. Phys. Japan 17, 48 (1935) 非常に多くの研究が精力的になされ続けている重要なターゲット。

Introduction (cont’d) 核力の性質は空間的に3つの領域に分割されて理解される。 長距離領域 (r > 2fm) OPEP(one pion exchange)が支配的。 中間距離領域 (1fm < r < 2fm) 中間子論の立場からアクセス可能。 heavier meson exchanges such as “σ”, ρ, ω,... multi pion exchange 近距離領域 (r < 1fm) この領域の理解は遅れている。 現象論的アプローチ (1) phenomenological repulsive core model (2) vector meson exchange model (3) quark exchange (constituent quark model) この領域の核力は次を反映していることが期待されている。 (1) the internal structure of the nucleon (2) quark や gluonの自由度の存在 このため、QCDに基づいた核力の理解が長い間望まれていた。. The nuclear force

Various approaches to QCD ★ 高エネルギーではαS は弱くなり、 QCDは摂動論的取り扱いがうまくいく。 ★ 低エネルギーではαS は成長し、摂動論的取り扱いが破綻する。 ⇒ 様々な非摂動論的アプローチが試みられている。. quark models constituent quark models bag model non-linear sigma model Skymeon chiral effective models Nambu-Jona-Lasinio(NJL) model QCD sum rule Schwinger-Dyson(Bethe-Salpeter) equation Light front approach to QCD lattice QCD Monte Carlo calculation 有効理論的アプローチ QCD的アプローチ

Lattice QCD has been applied to many targets. inter-quark potential hadron mass coupling constants and form factors (moments of) structure function/parton distribution function 有限温度/密度QCD、 QCD phase diagram multi-hadron properties 散乱長 phase shift nuclei [bound state deuteron has not yet been obtained] ★現在の技術では次を同時に実現するのはまだ難しい: 物理的に軽い quark mass 大きな空間体積 dynamical quarks (i.e., unquenched QCD) カイラル対称性と相性のよい chiral quark action 細かい格子分割 ★ 状況は次の発展によって徐々に改善されている。 スーパーコンピュータの進歩 Monte Carlo 計算のアルゴリズムの進歩

Standard lattice QCD approach to the inter-quark potential The static qqbar potential qqbarの位置を固定するため、 m=∞のstatic quark とstatic anti-quark を導入する。 the qqbar potentials は total energy として計算される。 (massは ~∞より、運動エネルギーは~0) この方法はいろいろ拡張されて使われている。 (1) unquenched QCD with dynamical quarks (2) 有限温度/有限密度 (3) multi-quark 系: qqq, qqqq, qqqqq, ... (4) 違った color charge をもった仮想的なクォーク from G.S.Bali et al., PRD46,2636(’92)

ハドロン間ポテンシャルへの拡張の試み 次の文献で試みられている: A.Mihaly et al., PRD55, 3077(’97). C.Stewart et al., PRD57, 581(’98). C.Michael et al., UKQCD Collab., PRD60,054012(’99) T.T.Takahashi et al., AIP Conf.Proc.842,246(’06) T.Doi et al., AIP Conf.Proc.842,249(’06). (例) NN ポテンシャル 二つの核子の位置を固定するため、それぞれの核子中にstatic quarkが導入される。. light quarkを4つ導入して、2N系にする。 [(Qqq) + (Qqq)] この light quarkは二つの static quarkに巻きつくと仮定する。 核子の運動エネルギーを無視することで、NNのポテンシャルエネルギーを計算する。. よく考えられた方法なのだが今のところうまくいっていない。 一つには近距離では右のような再配置がおこる可能性があり、 期待した結果に結びつかないらしい [(QQq) + (qqq)]

我々は全く違う方法を用いる。 ここでは、この方法を核子系に拡張し核力を計算する。 この方法はCP-PACS collaborationによってππの散乱長を研究するために使われた。 CP-PACS collab., S. Aoki et al., PRD71,094504(2005) この方法によれば、 (1) 格子 QCDを使って、NNの波動関数を計算する。 (2) NNポテンシャルは、このNNの波動関数が Schroedinger 方程式を満たすことを要請して逆算する。 ※ A GOOD FEATURE: 先ほどの static quark potential を拡張した方法にくらべ、 この方法では核子の位置を固定するために fictitious static quarkを導入する必要がない。 ⇒ 実験データを解析して得られるNN potentialにより近いものが得られる。 ここでは、この方法を核子系に拡張し核力を計算する。

The Formalism The (effective) non-relativistic Schroedinger eq. for NN system. NN のBethe-Salpeter(BS) 方程式から 3D reduction をつかって導かれる。 相互作用 kernel U(r-r’) は一般には non-local でtotal energy E に依存する。 相互作用kernel U(r-r’)の可能な形は様々な対称性により制限される。 低エネルギーで相互作用 kernel U(r-r’) は次のようにparameterizeされる。 ⇒ the familiar NN Schroedinger equation: 波動関数φ(r) があらかじめ分かっている場合, ポテンシャルは次のように書ける。 ←schematical

1S0 channel (計算が一番単純化される) ★ L=0, S=0 ⇒ central force しか残らない ! ★ 通常のNNのSchroedinger 方程式に帰着する。 ★ VC(r) は通常の関数である。 ⇒ central force VC(r) は次のように表せる。 あらかじめ波動関数φ(r)が知られている場合、 これをつかってポテンシャルを計算できる。 ※ 同じことを 3S1 channel でやると、 テンソル力とLS力がeffeciveに繰り込まれた effective central force VCeff(r) が得られる。 ※まじめにやるには例えば3S1, 3D1, 3D2を 使って coupled eq を解く。 ⇒ 同時にVLS(r), VT(r)も求まる。

Bethe-Salpeter (BS) 波動関数 QCDでは、non-rela. のNN 波動関数は本当は近似的な概念である。 この概念に最も近いのは、 Bethe-Salpeter (BS) 波動関数と呼ばれるものである。 このamplitudeは non-rela. のp-n の波動関数と解釈できる。 QCDのpn状態。非常に直感的には次のようなものを想像するればよい。 点xにいる proton を standard proton operator を用いて eliminate する。 点yにいる neutron を standard neutron operatorを用いて eliminateする。 最終的に その結果と、QCD真空 <0| との overlapを計算することで、 non-rela の波動関数 φ(x,y) を取り出せる。 non-rela 波動関数 neutron creation op. proton creation op.

格子QCDによるBS波動関数の求め方: 非弾性散乱に対応する成分 ★ 十分大きな t では、励起状態からの寄与は指数関数的に suppress される。 基底状態のNN系のBS波動関数を取り出せる。 これを先ほどの式に代入して NNポテンシャルを求める。

Lattice QCD parameters Quenched QCD is used. (virtual qqbar 対生成/消滅が起こらない近似) Gauge action: the standard Wilson gauge action. β= 6/g2 = 5.7(⇔ゲージ結合定数) 格子間隔: a ~0.14 fm (from ρmass in the chiral limit) 体積: 324 lattice (L~4.5 fm) 統計: 500 gauge configs (3000 sweeps for thermalization. The gauge config is separated by 200 sweeps) Quark action: Standard Wilson quark action κ=0.1665(⇔クォーク mass): mπ ~0.53 GeV, mρ~0.88 GeV, mN~1.34 GeV (Monte Carlo 計算は quark massが軽いほどどんどん難しくなる) Dirichlet BC along temporal direction Wall source on the time-slice t=5 NN wave function is measured on the time-slice t=11 Blue Gene/L at KEK has been used for the Monte Carlo calculations.

The Lattice QCD result for NN wave function r > 0.7 fm では、座標軸のnearest neighbor のみで計算が行われている。 (クラゲのような形をしているのはその為。 全ての点を計算すると60倍大変になる!) 波動関数の回転対称性が見える。 ★ 近距離でのshrink ⇒ 斥力の存在 convergence against t

The Lattice QCD result of NN potential The effective central force for 3S1. t に関する収束性 ※ 長距離で dominant になる diagram は現在の計算に 取り込まれている。 長距離でone π exch に対応。 中間距離の引力 ~ 30 MeV 重い mπ のせいで引力が弱くなっている。. ⇒ mπ が軽くなると、virtual pion は 長距離をpropagateできるようになる。 ⇒ stronger attraction. 3S1 の有効中心力は1S0より若干つよい。 (deuteron に都合がよい) 近距離の斥力芯 ~ 500~600MeV constituent quark model によると、 軽いquark mass 領域で enhance される。

quark mass 依存性 中間距離の引力が若干成長? PRELIMINARY mπ~530 MeV mπ~370 MeV 軽い quark mass 中間距離の引力が若干成長? 近距離の斥力芯が急激に成長。 This behavior is consistent with the constituent quark picture --- color magnetic interaction is propotional to 1/(m1 m2). 軽い quark mass でMonte Carlo 計算することが、 実験と比較できる核力を得るために必至。(今後の課題) いろいろなquark mass

This NN potential is attractive. ポテンシャルのrangeの外のs-wave波動関数は、低エネルギー(E → 0)で次のように近似される。 散乱長 a0 は rφasy(r)のlinear fitの r-切片になる。 r-切片は原点よりも左側にある。 a0 = 0.066(22) fm (attractive) the net interaction is attraction. r-intercept ※ comments: このポテンシャルは 自動的に正しい散乱長を出す 我々のNNポテンシャルに通常の散乱理 論を適用して求めた散乱長と格子QCDで 散乱長を求める標準的な方法(Luescher の方法)は L→∞で等価であることが示せ る。 a0~0.066 fm は実験値 a0~20 fm と比べ て非常に小さい。これはpionが重たいため mπ~530 MeV. 1S0 ★ This is understood by the formula of the Born approx: The relation tells us The volume factor r2dr hides the repulsive core at short distance a0 is a subtle quantity obtained as a result of a big cancellation between the repulston and the attraction.

波動関数の不定性 はprotonやneutronに対するわりあいよいoperatorであるが、 パーフェクトではない。 proton や neutron 状態とだけ couple すればよいが、 それ以外に、余計な状態とも結合する。 たとえば、 ここから余計な相関が生じる危険性がある。 etc このinelastic contributionは低エネルギーでは抑えられることを期待している。 しかしながら、近い将来その大きさを estimate する必要がある。

Summary 格子QCDを用いて NN potentialを計算した。 最近 CP-PACS collaborationによってππ散乱長の研究に用いられた方法を 拡張して使った。 核力の定性的な部分は再現されていた。 近距離(r < 0.5 fm)の斥力芯(~600 MeV). 中間距離 (0.5 fm < r < 1.2 fm) の引力(~30 MeV). (The attraction is weak due to the heavy pion (mπ~530 MeV).) 将来計画: 斥力芯の物理的期限 (dependences on the quark mass, the flavor structure, ...) ハイペロン相互作用 (ΣN, ΛN, ΛΛ, ΞN, etc.) mson-baryon 系。 LS力、テンソル力 (そして、3体力) ポテンシャルのE依存性、非局所性。 unquenched QCD, physical quark mass, finer discretization, larger spatial volume, ...