Hajime Togashi - International Workshop on Neutrino Physics and

Variational study of nuclear equation of state
for core-collapse supernovae
and hyperonic neutron stars
H. Togashi (RIKEN)
Outline
1:Introduction
2:Variational EOS for core-collapse supernovae
3:Variational EOS for hyperonic neutron stars
4:Summary
International Workshop on Neutrino Physics and Astrophysics
@ Mimar Sinan Fine Arts University, Istanbul, Turkey, March 16, 2015
1. Introduction
Nuclear Equation of State (EOS) plays important roles
for astrophysical studies.
(Neutron stars, core-collapse supernovae (SNe), black hole formations)
Nuclear EOS available for SN simulation
Thermodynamic quantities in a wide range of ρ, T, Yp
SN matter contains uniform and non-uniform nuclear matter.
1. Lattimer-Swesty EOS : The Skyrme-type interaction (NPA 535 (1991) 331)
2. Shen EOS : The Relativistic Mean Field Theory (NPA 637 (1998) 435)
Those EOSs are based on phenomenological models for uniform matter.
There is no nuclear EOS based on the microscopic many-body theory.
We construct a new SN-EOS with the variational method
starting from the realistic nuclear forces.
EOS with Variational Method
Fermi Hypernetted Chain (FHNC) method
Zero temperature : APR
A. Akmal et al., PRC58(1998)1804
Finite temperature : AM
A. Mukherjee, PRC 79(2009) 045811
Potential: AV18+UIX
Wave function: Jastrow wave function
For Symmetric Nuclear Matter and Pure Neutron Matter
SN-EOS table must cover in the following wide ranges.
Density ρ : 105.1 ≤ ρm ≤ 1016.0g/cm3
110 point
Temperature T : 0 ≤ T ≤ 400 MeV
92 point
Proton fraction Yp : 0 ≤ Yp ≤ 0.65
66 point
We have to treat Asymmetric Nuclear Matter directly.
SN-EOS for Uniform Nuclear Matter is constructed
with the simplified cluster variational method.
Our Plan to Construct the EOS for SN Simulations
Collaborators : M. Takano (Waseda Uniersity) ,
Y. Takehara, S. Yamamuro, K. Nakazato, H. Suzuki (Tokyo Univ. of Science)
K. Sumiyoshi (Numazu College of Tech.)
Uniform Nuclear Matter
EOS is constructed with the cluster variational method
Non-uniform Nuclear Matter
EOS is constructed with the Thomas-Fermi calculation
★We are here.★
Completion of a Nuclear EOS table for SN simulations
Density ρ : 105.1 ≤ ρm ≤ 1016.0g/cm3
110 point
Temperature T : 0 ≤ T ≤ 400 MeV
92 point
Proton fraction Yp : 0 ≤ Yp ≤ 0.65
66 point
2. Variational EOS for core-collapse supernovae
EOS for uniform nuclear matter at zero temperature Two-body Hamiltonian
Three-body Hamiltonian
The AV18 two-body nuclear potential
Jastrow wave function
ΦF: The Fermi-gas wave function
fij : Correlation function
Central
The UIX three-body nuclear potential
Ptsµ: Spin-isospin projection operators�
Tensor
Spin-orbit
E2/N is the expectation value of H2 in the two-body cluster approximation.
E2/N is minimized with respect to fCtsµ(r), fTtµ(r) and fSOtµ(r)
with the appropriate constraints. Two-Body Energy E2/N
APR : PRC58(1998)1804
Our results are in good agreement with the results by
APR (FHNC method).
Three-Body Energy
UIX potential
Three-Body Energy
: 2π exchange part
: Repulsive part
Correction term
Expectation value with the Fermi-gas wave function
α, β, γ, δ : Parameters in E3/N
・ Total energy per nucleon E/N = E2/N + E3/N reproduce
the empirical saturation properties.
・ Thomas-Fermi calculation of isolated atomic nuclei with E/N
reproduces the gross feature of the experimental data.
13/32
Total Energy per Nucleon E/N
Our EOS : NPA902 (2013) 53
APR : PRC58(1998)1804
n0[fm-3]
E0[MeV]
K [MeV]
Esym[MeV]
0.16
-16.1
245
30.0
Application of EOS to Neutron Stars
J0348+0432: Science 340 (2013) 1233232
J1614-2230: Nature 467 (2010) 1081
Shaded region: Steiner et al.,
Astrophys. J. 722 (2010) 33
The NS mass-radius relation is consistent with observational data.
EOS of Uniform Nuclear Matter at Finite Temperature
We follow the prescription proposed by Schmidt and Pandharipande.
(Phys. Lett. 87B(1979) 11)
(A. Mukherjee et al., PRC 75(2007) 035802)
Free energy per nucleon at T=20MeV
Our EOS : NPA902 (2013) 53
Free energy per nucleon at T=30MeV
FHNC : A. Mukherjee, PRC 79(2009) 045811
EOS for Non-uniform Nuclear Matter
1013)
We follow the Thomas-Fermi (TF) method by Shen et. al. (PTP100(1998)
(APJS 197(2011) 20)
Free energy in the Wigner-Seitz (WS) cell
Bulk energy
Gradient energy
Coulomb energy
Local free energy density : f = fN + fα
F0 = 68.00 MeV fm5 a : Lattice constant
f N: Free energy density of uniform nuclear matter
Parameter
Minimum
log10(T) [MeV]
-1.08
Maximum
1.52
Number
66 + 1
Yp
0.0
0.5
213
nB [fm-3]
0.000001
0.18
1980
67×213×1980
= 28256580 points
F/Vcell is minimized with respect to parameters in np (r), nn (r), nα (r)
for various densities, temperatures and proton fractions.
Thermodynamic Quantities
Mass Number and Proton Number
Mass number A
Proton number Z
Shen EOS (APJ Suppl.197 (2011) 20)
Application to supernova simulation
1.  Fully GR spherically symmetric simulation�(1D)
2. �Adiabatic collapse�(Ye is fixed during simulation.)
3. Progenitor: Woosley Weaver 1995, 15M◉
Astrophys. J. Suppl. 101 (1995) 181
SN simulation numerical code: K. Sumiyoshi, et al., NPA 730 (2004) 227
Density profiles at the bounce
Adiabatic indices at the bounce
Variational EOS is softer than Shen EOS in the adiabatic simulation.
Compositions of Central Core
Mass, neutron and proton numbers of nuclei
The mass fraction of particles
Xα and Xp at the initial profile are larger than those with the Shen EOS.
Our new SN-EOS table will be available soon!
3. Variational EOS for hyperonic neutron stars
GR
2.5
Hyperons (Λ, Σ, Ξ ) mixing soften
an EOS of neutron star matter.
<∞
AP3
y
it
sal
au
C
2.0
AP4
2
MS2
J1614-2230
SQM1
1.5
1
PAL1
ENG
SQM3
FSU
J1903+0327
Mass (M()
HYPERON PUZZLE
(HYPERON CRISIS?)
P
R
MS0
MPA1
4
GM3
PAL6
J1909-3744
3
MS1
5
GS1
Double neutron star
s sy
systems
6
1.0
7
8
9
tion
0.5
Rota
1
0.0
7
8
9
10
11
12
Radius (km)
13
14
15
(P. B. Demorest et al., NATURE 467 (2010)
Figure 3 | Neutron star mass–radius diagram. The plot shows non-rotating
mass versus physical radius for several typical EOSs27: blue, nucleons; pink,
nucleons plus exotic matter; green, strange quark matter. The horizontal bands
show the observational constraint from our J1614-2230 mass measurement of
(1.97 6 0.04)M[, similar measurements for two other millisecond pulsars8,28
and the range of observed masses for double neutron star binaries2. Any EOS
line that does not intersect the J1614-2230 band is ruled out by this
measurement. In particular, most EOS curves involving exotic matter, such as
kaon condensates or hyperons, tend to predict maximum masses well below
2.0M[ and are therefore ruled out. Including the effect of neutron star rotation
increases the maximum possible mass for each EOS. For a 3.15-ms spin period,
this is a =2% correction29 and does not significantly alter our conclusions. The
grey regions show parameter space that is ruled out by other theoretical or
observational constraints2. GR, general relativity; P, spin period.
There is no modern hyperon EOS with the variational method.
•  Our cluster variational method for asymmetric nuclear matter
is extended to that for hyperonic nuclear matter. common feature of models that include the appearance of ‘exotic’
hadronic matter such as hyperons4,5 or kaon condensates3 at densities
of a few times the nuclear saturation density (ns), for example models
GS1 and GM3 in Fig. 3. Almost all such EOSs are ruled out by our
results. Our mass measurement does not rule out condensed quark
matter as a component of the neutron star interior6,21, but it strongly
constrains quark matter model parameters12. For the range of allowed
EOS lines presented in Fig. 3, typical values for the physical parameters
of J1614-2230 are a central baryon density of between 2ns and 5ns and a
radius of between 11 and 15 km, which is only 2–3 times the
Schwarzschild radius for a 1.97M[ star. It has been proposed that
the Tolman VII EOS-independent analytic solution of Einstein’s
equations marks an upper limit on the ultimate density of observable
cold matter22. If this argument is correct, it follows that our mass mea-
•  We study contributions from bare hyperon interactions
to neutron star structure.
Collaborators : E. Hiyama (RIKEN) ,
M. Takano (Waseda University), Y. Yamamoto (RIKEN)
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
Cluster Variational Method for hyperon matter
Two-body Hamiltonian
Two-body potential YN ・YY potential: Central potential
(µ = Λp, Λn, ΛΛ, Σ-p, Σ-n)
Jastrow wave
function
Two-body correlation function
YN・YY correlation : Central component
(µ = Λp, Λn, ΛΛ, Σ-p, Σ-n)
<H2> is calculated in the two-body cluster approximation.
ΛΛ Interaction and EOS of Neutron Star Matter
ΛΛ potential(Triplet-odd state)
Pressure of neutron star matter
As the odd-state part of ΛΛ interaction becomes repulsive,
the pressure of neutron star matter increases.
Particle Fraction in Neutron Star Matter
Λ threshold density nthΛ : 0.46 fm-3
Σ- threshold density nthΣ : 0.92 fm-3 → 0.82 fm-3 → 0.77 fm-3 → 0.72 fm-3
Gravitational Mass of Neutron Stars
As the odd-state part of ΛΛ interaction becomes repulsive,
maximum mass of neutron star increases about 10%.
1.45M◉ → 1.50M◉→ 1.55M◉→ 1.60M◉ EOS of Neutron Star Matter with Three Baryon Force
We supplement an phenomenological Three Baryon Forces (TBF)
�as a density dependent two-body effective potential. Y. Yamamoto et al., PRC 90 (2014) 045805
Λ threshold density nthΛ : 0.59 fm-3
Maximum mass of NS : 2.12 M◉
4. Summary
1. New nuclear EOS for SN simulations is constructed
with the variational method.
•  Nuclear Hamiltonian composed of realistic Nuclear forces (AV18 + UIX).
•  Variational EOS is softer than Shen EOS in the adiabatic simulation.
Our EOS is advantageous for SN explosion! Future Plan
•  SN simulations with the neutrino transfer
2. EOS for hyperonic neutron star is constructed
with the extended variational method.
•  Cluster variational method for hyperonic nuclear matter is developped .
•  Odd-state part of ΛΛ interaction affect the structure of neutron stars.
Future Plan
•  Modern NY and YY interaction (e.g. Nijmegen)