Quantum Monte Carlo calculations of neutron and nuclear matter Stefano Gandolfi Los Alamos National Laboratory (LANL) Advances and perspectives in computational nuclear physics, Hilton Waikoloa Village, Kona, October 5-7, 214 www.computingnuclei.org
Nuclei SciDAC UNEDF/NUCLEI
Homogeneous neutron matter
Outline The model and the method Neutron matter EOS, symmetry energy and neutron stars Nuclei and nuclear matter: role of the NN Hamiltonian Conclusions
Nuclear Hamiltonian Model: non-relativistic nucleons interacting with an effective nucleon-nucleon force (NN) and three-nucleon interaction (TNI). H = 2 2m A i=1 2 i + i<j v ij + i<j<k v ij NN fitted on scattering data. Sum of operators: V ijk v ij = O p=1,8 ij v p (r ij ), O p ij = (1, σ i σ j, S ij, L ij S ij ) (1, τ i τ j ) Argonne AV6 (no LS), AV7 (no LS-τ), AV8. Local chiral forces up to N 2 LO has the same spin/isospin operatorial structure than AV7 - Gezerlis, Tews, et al. PRL (213).
Nuclear Hamiltonian Phase shifts, Argonne AV6, AV7 and AV8 6 4 1 S Argonne V6 Argonne V7 Argonne V8 SAID 12 1 8 3 S1 Argonne V6 Argonne V7 Argonne V8 SAID 1 P1 Argonne V6 Argonne V7 Argonne V8 SAID δ (deg) 2 δ (deg) 6 4 δ (deg) -2 2-2 δ (deg) -4 1 2 3 4 5 6 E lab (MeV) 2-2 Argonne V6 Argonne V7 Argonne V8 SAID 3 P δ (deg) -2 1 2 3 4 5 6 E lab (MeV) -1-2 -3-4 3 P1 Argonne V6 Argonne V7 Argonne V8 SAID δ (deg) -4 1 2 3 4 5 6 E lab (MeV) 3 2 1 Argonne V6 Argonne V7 Argonne V8 SAID 3 P2-4 1 2 3 4 5 6 E lab (MeV) -1 3 D1 Argonne V6 Argonne V7 Argonne V8 SAID -5 1 2 3 4 5 6 E lab (MeV) 6 5 4 Argonne V6 Argonne V7 Argonne V8 SAID 3 D2-1 1 2 3 4 E lab (MeV) 8 6 ε 1 δ (deg) -2 δ (deg) 3 δ (deg) 4-3 -4 2 1 2 Argonne V6 Argonne V7 Argonne V8 SAID 1 2 3 4 5 6 E lab (MeV) 1 2 3 4 5 6 E lab (MeV) 1 2 3 4 5 6 E lab (MeV)
Nuclear Hamiltonian Chiral interactions permit to understand the evolution of theoretical uncertainties with the increasing of A. Slide by Joel Lynn, Scidac NUCLEI meeting 214.
Three-body forces Urbana Illinois V ijk models processes like π π π π π π π π π π Chiral forces at N 2 LO: + short-range correlations (spin/isospin independent).
Nuclear Hamiltonians Advantages: Argonne interactions fit phase shifts up to high energies. At ρ = ρ, k F 33 MeV. Two neutrons have E CM 12 MeV, E LAB 24 MeV. accurate up to (at least) 2-3ρ. Provide a very good description of several observables in light nuclei. Interactions derived from chiral EFT can be systematically improved. Changing the cutoff probes the physics and energy scales entering into observables. They are generally softer, and make most of the calculations easier to converge. Disadvantages: Phenomenological interactions are phenomenological, not clear how to improve their quality. Systematic uncertainties hard to quantify. Chiral interactions describe low-energy (momentum) physics. How do they work at large momenta, (i.e. e and ν scattering)? Important to consider both and compare predictions
Quantum Monte Carlo The goal is to solve the many-body Schrödinger equation: H ψ( r 1... r N ) = E ψ( r 1... r N ) ψ(t) = e (H E T )t ψ() Ground-state extracted in the limit of t. Propagation performed by ψ(r, t) = R ψ(t) = dr G(R, R, t)ψ(r, ) Importance sampling: G(R, R, t) G(R, R, t) Ψ I (R )/Ψ I (R) Constrained-path approximation to control the sign problem. Trial wave function includes spin/isospin dependent two- and three-body correlations. GFMC includes all spin-states of nucleons in the w.f., nuclei up to A=12 AFDMC samples spin states, A 1, limitations to the form of H Ground state obtained in a non-perturbative way. Systematic uncertainties within 1-2 %.
Neutron matter equation of state Why do we care of neutron matter? EOS of neutron matter useful to study the symmetry energy and its slope at saturation. Main input for neutron stars. The three-neutron force (T = 3/2) very weak in light nuclei. The dominant T = 1/2 is zero in neutron matter. No direct T = 3/2 experimental probe available! Theory Esym, L Neutron stars Experiments
What is the Symmetry energy? symmetric nuclear matter pure neutron matter Symmetry energy E = -16 MeV ρ =.16 fm -3 Nuclear saturation Assumption from experiments: E SNM (ρ ) = 16MeV, ρ =.16fm 3, E sym = E PNM (ρ ) + 16 At ρ we access E sym by studying PNM.
Neutron matter We consider different forms of three-neutron interaction by only requiring a particular value of E sym at saturation. 12 Energy per Neutron (MeV) 1 8 6 4 2 E sym = 33.7 MeV different 3N.1.2.3.4.5 Neutron Density (fm -3 ) Gandolfi, et al., EPJA (214).
Neutron matter Equation of state of neutron matter using Argonne forces: Energy per Neutron (MeV) 1 8 6 4 2 E sym = 35.1 MeV (AV8 +UIX) E sym = 33.7 MeV E sym = 32 MeV E sym = 3.5 MeV (AV8 ).1.2.3.4.5 Neutron Density (fm -3 ) Gandolfi, Carlson, Reddy, PRC (212)
Neutron matter and symmetry energy From the EOS, we can fit the symmetry energy around ρ using E sym (ρ) = E sym + L 3 ρ.16.16 + 7 65 AV8 +UIX L (MeV) 6 55 5 45 4 E sym =32. MeV E sym =33.7 MeV 35 AV8 3 3 31 32 33 34 35 36 E sym (MeV) Gandolfi et al., EPJ (214) Tsang et al., PRC (212) Very weak dependence to the model of 3N force for a given E sym. Chiral Hamiltonians give compatible results.
Neutron star structure EOS used to solve the TOV equations. M (M O ) 3 2.5 2 1.5 1.5 Causality: R>2.9 (GM/c 2 ) 32 error associated with E sym ρ central =3ρ E sym = 3.5 MeV (NN) ρ central =2ρ 33.7 35.1 1.97(4) M O 1.4 M O 8 9 1 11 12 13 14 15 16 R (km) Gandolfi, Carlson, Reddy, PRC (212). Strong interplay between E sym and neutron star radii. more details Thur. at 11.15, session DL
Light nuclei spectrum computed with GFMC Energy (MeV) -2-3 -4-5 -6-7 -8 4 He + 6 He 2 + + 1 + 2 + 3 + 1 + 6 Li 7/2 5/2 5/2 7/2 1/2 3/2 7 Li Argonne v 18 with UIX or Illinois-7 GFMC Calculations 1 June 211 2 + + 4 + + 2 + 8 He 2 + 8 Li 1 + 3 + 1 + 2 + 8 Be + 1 + 4 + 3 + 1 + 2 + 4 + 2 + + 7/2 + 5/2 + 4 + 7/2 2 + 7/2 3 + 3/2 2 + 3/2 + 3 + 5/2 + 1 + 1/2 1 + 5/2 1/2 + 3/2 3 + 4 + 2 + 3 + 2 + 1 + 1 + 9 Be 3 + 1 B -9-1 AV18 AV18 +UIX AV18 +IL7 Expt. 12 C + Carlson, Pieper, Wiringa, many papers Also radii, densities, matrix elements,...
4 He energy with chiral two-body interactions. It s important to explore: the role of i) the chiral expansion and ii) the cutoff. E b (MeV) -16-18 -2-22 -24-26 -28-3 -42-44 -46-48 Exp. R = 1. fm R = 1.1 fm R = 1.2 fm LO NLO N 2 LO v 8 Chiral Order Lynn, Carlson, Epelbaum, Gandolfi, Gezerlis, Schwenk, arxiv:146.2787. Three-body terms not yet included. AFDMC gives results within 1% See Joel lynn talk Wed. at 1. pm, session CK
Nuclei AV6 AV7 exp 4 He -27.9(3) -25.7(2) -28.295 16 O -115.6(3) -9.6(4) -127.619 4 Ca -322(2) -29(1) -342.51 Gandolfi, Lovato, Carlson, Schmidt, arxiv:146.3388 About 1 to 4% of binding energy missing. Preliminary: 16 O binding energy with chiral potentials: R =1.2 fm R =1. fm LO -121.1(3) -269.9(4) NLO -87.4(2) -35.9(7) N 2 LO -116.1(2) -91.1(3) Different orders do not even overlap for both 4 He and 16 O.
Nuclear matter EOS of symmetric nuclear matter using Argonne AV6 and AV7 : energy per nucleon (MeV) -2-4 -6-8 -1-12 -14-16 -18-2 AV7 AV6-22 -24.5.1.15.2.25.3.35 ρ (fm -3 ) Gandolfi, Lovato, Carlson, Schmidt, arxiv:146.3388
Nuclear matter Preliminary results! EOS of symmetric nuclear matter with chiral forces: energy per nucleon (MeV) 4-4 -8-12 -16-2 R =1. fm NLO R =1.2 fm N 2 LO (no V3) R =1.2 fm R =1. fm Same qualitative behavior of nuclei..1.15.2.25.3 ρ (fm -3 )
Neutron matter Equation of state of neutron matter using NN chiral forces: 15 AFDMC LO AFDMC NLO AFDMC N 2 LO R =1. fm E/N [MeV] 1 5 R =1.2 fm.5.1.15 n [fm -3 ] Gezerlis, Tews, Epelbaum, SG, Hebeler, Nogga, Schwenk, PRL (213), and arxiv:146.454 (214)
Nuclear/neutron matter E/N [MeV] 15 1 5 AFDMC LO AFDMC NLO AFDMC N 2 LO R =1. fm R=1.2 fm energy per nucleon (MeV) 4-4 -8-12 -16 R =1. fm NLO R =1.2 fm N 2 LO (no V3) R =1.2 fm.5.1.15 n [fm -3 ] -2 R =1. fm.1.15.2.25.3 ρ (fm -3 ) Nuclear matter is more correlated than neutron matter. Open questions: Perturbative vs non-perturbative? Importance of phase shifts? Importance of the (not included yet) V3?
Nuclear/neutron matter Energy per neutron (MeV) 4 3 2 1 AV8 +UIX AV8 /AV7 AV6.5.1.15.2.25 ρ (fm -3 ) energy per nucleon (MeV) -2-4 -6-8 -1-12 -14-16 -18-2 AV7 AV6-22 -24.5.1.15.2.25.3.35 ρ (fm -3 ) Same behavior with Argonne forces.
Summary AFDMC/GFMC methods useful to study nuclear systems in a coherent framework: same Hamiltonians and same many-body machinery Three-neutron force is the (strong) bridge between E sym and neutron star structure. Chiral potentials show some convergence in light nuclei (but three-body terms not included yet). Error bars given by varying the cutoff grow quickly with the size of the system. Nuclear and neutron matter have very different behavior with both Argonne and chiral forces. Further investigation needed. Due to the phase shifts?
Acknowledgments J. Carlson, J. Lynn (LANL) A. Gezerlis (U. Guelph) A. Lovato, S. Pieper, R. Wiringa (ANL) A. Steiner (INT) R. Schiavilla (Jlab/ODU) K. Schmidt (ASU) I. Tews, A. Schwenk (TU Darmstadt) www.computingnuclei.org