Validity of The Wigner-Seitz Approximation in Neutron Star Crust
Validity of The Wigner-Seitz Approximation in Neutron Star Crust
Validity of The Wigner-Seitz Approximation in Neutron Star Crust
N. Chamel
Institut dAstronomie et dAstrophysique, Universite Libre de Bruxelles,
CP226, Boulevard du Triomphe, 1050 Brussels, Belgium
Institut de Physique Nucleaire, Universite Paris-Sud, IN2 P3 -CNRS, F-91406 Orsay Cedex, France
(Dated: February 5, 2008)
Since the seminal work of Negele and Vautherin, the Wigner-Seitz approximation has been widely
applied to study the inner crust of neutron stars formed of nuclear clusters immersed in a neutron
sea. In this article, the validity of this approximation is discussed in the framework of the band
theory of solids. For a typical cell of 200 Zr, present in the external layers of the inner crust, it is shown
that the ground state properties of the neutron gas are rather well reproduced by the Wigner-Seitz
approximation, while its dynamical properties depend on the energy scale of the process of interest or
on the temperature. It is concluded that the Wigner-Seitz approximation is well suited for describing
the inner crust of young neutron stars and the collapsing core of massive stars during supernovae
explosions. However the band theory is required for low temperature transport properties as, for
instance, the effective neutron mass giving rise to entrainment effects.
PACS numbers: 97.60.Jd,26.60.+c,21.10.Ma,21.60.Jz,71.15.Ap,71.18.+y,21.10.-k
Keywords: neutron star crust, band theory, Wigner-Seitz approximation, Hartree-Fock
(1)
2
electrons as protons. The effects of the ion lattice on
the electrons, which give rise to complicated band structures in ordinary terrestrial matter are negligible in the
inner crust of a neutron star due to the very high densities [16]. Nevertheless, the neutron band effects due to
nuclear inhomogeneities cannot be ignored.
In the present study, we consider the outermost layers
of the inner crust where pairing effects are negligible [17].
In the Hartree-Fock approximation with Skyrme forces
which we shall consider in the following, the occupied
nucleon single particle wave functions are obtained by
solving the self-consistent equations (q = n, p for neutrons and protons respectively)
(q) (q)
(q)
(q)
h0 kk (rr ) = kk kk (rr )
(2)
h0
2
h
Wq (rr ) ,
+ Uq (rr ) iW
r)
2m
q (r
(3)
(q)
(4)
(q)
,kk +K
K = k
k .
(5)
(q)
(q)
(q)
(q)
(6)
(q)
hk
2 k2
h
hk ,
+ vq
r)
2m
q (r
(7)
B.
Wigner-Seitz approximation
The spherical W-S approximation is a computationally very efficient method with the advantage of reducing the 3D partial differential Eqs. (2) to ordinary differential radial equations. This approximation is twofold.
(q)
First of all, the Hamiltonian hk in equation (6) is neglected. Consequently the wave functions and the energies are independent of k and approximated by the
solutions at k = 0. Only the band index remains.
Secondly, the W-S polyhedron is replaced by a sphere
of equal volume. The equations are then usually solved
with the Dirichlet-Neumann mixed boundary conditions
which yield a nearly constant neutron density outside the
cluster.
The W-S approximation turns out to be very good if
the boundary conditions play a minor role. For instance,
bound states whose associated wave functions are vanishingly small outside the clusters are very well treated
provided that the spatial extent of these states is smaller
than the lattice spacing. This condition is fulfilled almost everywhere in the crust except in the bottom layers
where the clusters nearly touch. The aim of this paper is
to investigate the validity of the W-S approximation for
the outermost layers of the inner crust where the bound
neutron states are not altered by the boundary conditions.
Let us emphasize that in the W-S approximation, the
nuclear clusters are supposed to be spherical while in
the full band theory, no assumption is made about their
shape. For the low densities of interest in this study, the
nuclear clusters can still be considered as spherical. It
should mentioned that in a recent development of the
W-S approximation [18], the W-S cell is replaced by a
cube with strictly periodic boundary conditions. Possible deformations of the nuclear clusters are thus included
but at the price of unphysical boundary conditions because the W-S cell of the body centered cubic lattice is a
truncated octahedron and not a cube (the cube being the
W-S cell of the simple cubic lattice). This is the reason
why we still consider the spherical W-S approximation
closer to the physical situation than the cubic one at low
density.
0.14
0.12
10
-2
WS
BT
0.1
10
-3
-3
(fm )
(fm-3)
0.08
0.06
0.04
10
-4
0.02
0
10
12
10
-5
0 5 10 15 20 25 30 35 40 45 50
r (fm)
r (fm)
A.
10 2
10
WS
BT
10
-1
0 5 10 15 20 25 30 35 40 45 50
r (fm)
FIG. 3: Integrated unbound neutron number (see text) calculated with the W-S approximation (WS, full lines) and the
band theory (BT, dotted lines)
B.
0.5
0.4
E (MeV)
0.3
0.2
0.1
0
N
5
0.6
0.5
400
4p3/2 ,4s1/2
2h11/2
,2h9/2
350
0.3
0.2
-1
g(E) (MeV )
E (MeV)
0.4
4p1/2
2g9/2 ,2g7/2
3d5/2 ,3d3/2
3p1/2
2f7/2 ,2f5/2
0.1
0
3p3/2
300
250
200
150
100
50
XZ
BZ
d3k
(E kk )
(2)3
(11)
(12)
E (MeV)
2m
h2
3/2 p
E Ev ,
(14)
10 2
10
WS
BT
Fermi
10
-1
E (MeV)
CONCLUSION
[1] G.A. Baym, H.A. Bethe and C.J. Pethick, Nucl. Phys. A
175 (1971), 225.
[2] C. Kittel, Introduction to solid state physics, 7th edition,
Wiley&Sons (1996).
[3] E.P. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933).
[4] E.P. Wigner and F. Seitz, Phys. Rev. 46, 509 (1934).
[5] J.W. Negele and D. Vautherin, Nucl. Phys. A 207
(1973), 298.
[6] P. Bonche, D. Vautherin, Nucl. Phys. A 372 (1981), 496.
Acknowledgments
8
[12] M. Baldo, E.E. Saperstein and S.V. Tolokonnikov,
Nucl.Phys. A 775 (2006), 235.
[13] B. Carter, N. Chamel, P. Haensel, Nucl. Phys. A 748
(2005), 675.
[14] N. Chamel, Nucl.Phys. A 747 (2005), 109.
[15] N. Chamel, Nucl.Phys. A 773 (2006), 263.
[16] C.J. Pethick and D.G. Ravenhall, Ann. Rev. Nucl. Part.
Sci. 45 (1995), 429.
[17] M. Baldo, E.E. Saperstein and S.V. Tolokonnikov, Nucl.
Phys. A 749 (2005), 42.
[18] P. Magierski, A. Bulgac, P.-H. Heenen, Nucl. Phys. A
719 (2003), 217.
[19] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. A 635 (1998), 231 ; Erratum, Nucl.
Phys. A 643 (1998), 441.
[20] G.F. Koster, Solid State Physics 5, edited by F. Seitz, D.
Turnbull, Academic Press, New York (1957), 173.
[21] M. Baldo, E.E. Saperstein, S.V. Tolokonnikov, arxiv
preprint nucl-th/0609031.
[22] J. Margueron, J. Navarro, and P. Blottiau, Phys. Rev. C
70, 028801 (2004).
[23] D. Page, U. Geppert, F. Weber, Nucl. Phys. A 777
(2006), 497.