Quantum Monte Carlo calculations of two neutrons in finite volume

Transcription

1 Quantum Monte Carlo calculations of two neutrons in finite volume Philipp Klos with J. E. Lynn, I. Tews, S. Gandolfi, A. Gezerlis, H.-W. Hammer, M. Hoferichter, and A. Schwenk Nuclear Physics from Lattice QCD INT, Seattle, April 6, 2016 April 6, 2016 Institut für Kernphysik Philipp Klos 1

2 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 2

3 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 2

4 Two neutrons in a box April 6, 2016 Institut für Kernphysik Philipp Klos 3

5 Motivation Lattice QCD to describe nuclei from QCD April 6, 2016 Institut für Kernphysik Philipp Klos 4

6 Motivation I Lattice QCD to describe nuclei from QCD I marekich/wikimedia However, constrained to finite volume and very few particles April 6, 2016 Institut für Kernphysik Philipp Klos 4

7 Motivation I Lattice QCD to describe nuclei from QCD I Strategy: I Lattice QCD calculations for few nucleon systems I Matching of nuclear forces (EFT) in finite volume I EFT calculations of nuclear properties using advanced many-body methods April 6, 2016 Institut für Kernphysik Philipp Klos 4

8 Motivation Why two neutrons? Comparison to Lüscher formula possible "proof of principle" Low-energy constants (LEC) in chiral EFT are fitted to experimental data Neutron-neutron scattering length cannot be measured directly Why QMC? Exact method to solve Schrödinger equation Very successful for light nuclei Naturally extendable to more particles (3n, 4n,...) April 6, 2016 Institut für Kernphysik Philipp Klos 5

9 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 6

10 Two neutrons in finite volume Lüscher formula Consider interaction of nucleons N through contact interactions ( ) L = Ñ i t + 2 Ñ + 2M S-wave scattering amplitude A Now use A = 4π M 1 p cot δ ip. ( ) 4 D µ [C 0 (Ñ Ñ) 2 + C (1) 2 2 (Ñ Ñ) 2 + C (2) 2 [i (Ñ Ñ)] ] and periodic boundary conditions of the box to obtain energies. April 6, 2016 Institut für Kernphysik Philipp Klos 7

11 Two neutrons in finite volume Lüscher formula Energy E = p2 of two nucleons in a box predicted by infinite volume phase shift δ(p) m M. Lüscher, Commun. Math. Phys. 105, 153 (1986). p cot δ(p) = 1 πl S S(η) = lim Λ (( ) 2 ) Lp 2π ( Λ Beane et al., PLB 585, 106 (2006). j 1 j 2 η 4πΛ ) pionless EFT Possible to infer scattering length, effective range,..., from finite volume calculations (p cot δ = 1 a r ep ) Finite volume energies scattering length a, effective range r e,... April 6, 2016 Institut für Kernphysik Philipp Klos 8

12 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 9

13 Quantum Monte Carlo QMC in three lines: Ground state: H Ψ 0 = E 0 Ψ 0 Trial state: Ψ T = α i Ψ i i Propagate: lim τ e Hτ Ψ T Ψ 0 QMC in more than three lines: J. Carlson et al., RMP 87, 1067 (2015). April 6, 2016 Institut für Kernphysik Philipp Klos 10

14 Quantum Monte Carlo Trial wave function The trial wave function written in terms of eigenstates of H Ψ T = α i Ψ i Propagate in imaginary time to project out the ground state Ψ 0 Ψ(τ) = e (H ET )τ Ψ ( T ) = e (E0 ET )τ α 0 Ψ 0 + α i e (Ei E0)τ Ψ i Ψ(τ) τ Ψ 0 i i 0 April 6, 2016 Institut für Kernphysik Philipp Klos 11

15 Quantum Monte Carlo Example Particle in a{ box 0, 0 < x < L V =, otherwise Trial wave function Ψ T (x) = Mixed estimator c n Ψ n (x) n=1 Ψ n (x) = 2 sin(nπx) E(τ) = Ψ T He (H ET )τ Ψ T Ψ T e (H ET )τ Ψ T c n ψ n (x) E(τ) n= n=5 n= x τ [1/E sep. ] E 1 April 6, 2016 Institut für Kernphysik Philipp Klos 12

16 Quantum Monte Carlo Example Particle in a{ box 0, 0 < x < L V =, otherwise Trial wave function Ψ T (x) = Mixed estimator c n Ψ n (x) n=1 Ψ n (x) = 2 sin(nπx) E(τ) = Ψ T He (H ET )τ Ψ T Ψ T e (H ET )τ Ψ T c n ψ n (x) E(τ) n= n=5 n= x τ [1/E sep. ] E 1 April 6, 2016 Institut für Kernphysik Philipp Klos 12

17 Quantum Monte Carlo Example Particle in a{ box 0, 0 < x < L V =, otherwise Trial wave function Ψ T (x) = Mixed estimator c n Ψ n (x) n=1 Ψ n (x) = 2 sin(nπx) E(τ) = Ψ T He (H ET )τ Ψ T Ψ T e (H ET )τ Ψ T c n ψ n (x) E(τ) n= n=5 n= x τ [1/E sep. ] E 1 April 6, 2016 Institut für Kernphysik Philipp Klos 12

18 Quantum Monte Carlo Trial wave function Trial wave function for N-body system: Ψ T (R, S) = f J (r ij )Ψ SD (R, S) i<j R = (r 1,..., r N ) and S = (s 1,..., s N ) Jastrow function f J (r ij ) Slater determinant Ψ SD (R, S) = det({φ α (r i, s i )}) φ α (r i, s i ) = e ikα ri χ s,ms (s i ) April 6, 2016 Institut für Kernphysik Philipp Klos 13

19 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 14

20 Two neutrons in finite volume Contact interaction Comparison: AFDMC vs. Lüscher prediction with contact interaction [ ( r ) ] 4 V (r) = C 0 exp R ground state pionless EFT E = q 2 4π2 ML 2 q AFDMC, C 0 physical a AFDMC, C 0 a= fm Lüscher, physical a Lüscher, a= fm L [fm] April 6, 2016 Institut für Kernphysik Philipp Klos 15

21 Two neutrons in finite volume Excited state Challenging problem: Excited state has nodal surface ψ(r node ) = 0 QMC requires nodal surface of wave function as input (fixed node approximation, plus small release possible) spherical nodal surface April 6, 2016 Institut für Kernphysik Philipp Klos 16

22 Two neutrons in finite volume Excited state Introduce node in Jastrow function f J (r 12 ) in Ψ T (R, S) = f J (r 12 )Ψ SD (R, S) f J (r 12 ) Ground state jastrow GS r 12 [fm] f J (r 12 ) Excited state jastrow EX r 12 [fm] April 6, 2016 Institut für Kernphysik Philipp Klos 17

23 Two neutrons in finite volume Excited state How to determine radius of the nodal surface? For local potential Schrödinger equation Hψ(R) = Eψ(R) yields the same energy E independent of coordinates R 0 r node r max Perform separate simulations in the two pockets until E left = E right April 6, 2016 Institut für Kernphysik Philipp Klos 18

24 Two neutrons in finite volume Excited state Adjust rnode such that Eleft = Eright 0 rnode rmax E [MeV] counts left pocket right pocket r [fm] April 6, 2016 Institut für Kernphysik Philipp Klos t [MeV ] 4

25 Two neutrons in finite volume Excited state Comparison: AFDMC vs. Lüscher prediction with contact interaction [ ( r ) ] 4 V (r) = C 0 exp R st excited state 0.66 q AFDMC, C 0, physical a Lüscher, physical a L [fm] April 6, 2016 Institut für Kernphysik Philipp Klos 20

26 Nodal surface Determine exact nodal surface through diagonalization of H ψ i = E i ψ i st excited state q AFDMC, C 0 physical a Diagonalization Lüscher, physical a Lüscher fit, a=-18.9(2) fm L [fm] April 6, 2016 Institut für Kernphysik Philipp Klos 21

27 Nodal surface Extract nodal surface r node (θ, ϕ) from first excited state ψ ex (r node, θ, ϕ) = 0 Nodal surface not spherical! April 6, 2016 Institut für Kernphysik Philipp Klos 22

28 Nodal surface April 6, 2016 Institut für Kernphysik Philipp Klos 23

29 Decomposition in cubic harmonics Nodal surface can be decomposed in spherical harmonics r node (θ, φ) = c lm Y lm (θ, φ) Rotation symmetry group is broken down to the cubic symmetry group O h J. Muggli, Z. Angew. Math. Mech. 23, 311 (1972). cubic harmonics Y c l = l m=0,4,8,... c m Y lm Express nodal surface in terms of cubic harmonics r node (θ, φ) = c l Y c l (θ, φ). l April 6, 2016 Institut für Kernphysik Philipp Klos 24

30 Decomposition in cubic harmonics Nodal surface decomposed in cubic harmonics r node (θ, φ) = c l Y c l (θ, φ). l cubic harmonics Y c l = c m Y lm m=0,4,8,... We find c l ( r node ) l L Two-body scattering in infinite volume R(r) j l (pr) pr 1 (pr) l c l L=10 fm L=20 fm L=30 fm l April 6, 2016 Institut für Kernphysik Philipp Klos 25

31 Decomposition in cubic harmonics Nodal surface decomposed in cubic harmonics r node (θ, φ) = c l Y c l (θ, φ). l cubic harmonics Y c l = c m Y lm m=0,4,8,... Estimate correction due to l = 4 contribution E = ψ l=0 H ψ l=0 + c 2 4 ψ l=4 H ψ l= L=10 fm L=20 fm L=30 fm c l E E (c 4) 2 1% 0.01 in agreement with results l April 6, 2016 Institut für Kernphysik Philipp Klos 25

32 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 26

33 Chiral effective field theory Separation of scales: Momentum Q breakdown scale Λ b Most general Lagrangian consistent with QCD symmetries Expansion in powers of Q Λ b Local forces up to N 2 LO by using linearly indep. operators Gezerlis et al., PRL 111, (2013) [Weinberg, van Kolck, Kaplan, Savage, Wise, Epelbaum, Hammer, Kaiser, Meißner,...] April 6, 2016 Institut für Kernphysik Philipp Klos 27

34 Limitations of Lüscher for pionful theory Chiral leading-order (LO) potential q ground state AFDMC, LO, R 0 =1.0 fm AFDMC, LO, R 0 =1.2 fm Lüscher, a=-18.9 fm, r=2.01 fm Lüscher, a=-18.9 fm, r=2.15 fm L [fm] q st excited state AFDMC, LO R 0 =1.0 fm Lüscher, a= 18.9 fm, r=2.01 fm L [fm] Lüscher formula assumes zero-range interaction (pionless EFT, p < m π /2) Not applicable for nuclear interactions at small box sizes L Direct matching of lattice QCD and chiral EFT necessary for small L April 6, 2016 Institut für Kernphysik Philipp Klos 28

35 Limitations of Lüscher for pionful theory Chiral leading-order (LO) potential q st excited state AFDMC, LO R 0 =1.0 fm Lüscher, a= 18.9 fm, r=2.01 fm L [fm] δ 0 (p) AFDMC, LO R 0 =1.0 fm, excited state δ 0 (p) chiral LO ERE with a= 18.9 fm, r e =2.01 fm p [fm -1 ] Compare directly to phase shift Shows challenges for p > m π /2 p cot δ(p) = 1 (( ) 2 ) Lp πl S 2π April 6, 2016 Institut für Kernphysik Philipp Klos 29

36 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 30

37 Extraction of scattering parameters from finite volume results Energy E = p2 of two nucleons in a box predicted by infinite volume phase shift δ(p) m Beane et al., PLB 585, 106 (2006). p cot δ(p) = 1 πl S (( ) 2 ) Lp 2π Effective range expansion: 1 a r ep = 1 (( ) 2 ) Lp πl S 2π pionless EFT Finite volume energies scattering length a, effective range r e,... April 6, 2016 Institut für Kernphysik Philipp Klos 31

38 Fit to AFDMC data Chiral NLO and N 2 LO potentials Fit to ground + excited state results st excited state st excited state q q AFDMC, C 0 physical a Diagonalization Lüscher, physical a Lüscher fit, a=-18.9(2) fm L [fm] Contact interaction C AFDMC, LO R 0 =1.0 fm 0.60 Lüscher, a= 18.9 fm, r=2.01 fm Lüscher fit, a=-18.5(4) fm L [fm] Chiral LO interaction Fit: a = 18.9(2) fm, r e = 1.10(1) fm Exact: a = 18.9 fm, r e = fm Fit: a = 18.5(4) fm, r e = 2.00(7) fm Exact: a = 18.9 fm, r e = 2.01 fm April 6, 2016 Institut für Kernphysik Philipp Klos 32

39 Motivation Lüscher formula Quantum Monte Carlo Ground and excited states in AFDMC Chiral effective field theory Extraction of scattering parameters from finite volume Summary and outlook April 6, 2016 Institut für Kernphysik Philipp Klos 33

40 Summary First results for two-neutron finite-volume ground and excited states in AFDMC Approximate construction of excited state vs. exact diagonalization Extraction of scattering parameters from AFDMC simulations yields accurate results QMC techniques can serve to match chiral EFT and lattice results beyond limitations of the Lüscher formula Outlook Generalizable to more particles (3n, 4n,...) where extensions of Lüscher s formula are only partially available Extraction of resonance properties through calculations of excited states arxiv: Thank you! April 6, 2016 Institut für Kernphysik Philipp Klos 34