11.07.2015 Views

Theory of Nuclear Matter for Neutron Stars and ... - Graduate Physics

Theory of Nuclear Matter for Neutron Stars and ... - Graduate Physics

Theory of Nuclear Matter for Neutron Stars and ... - Graduate Physics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

The energy exchange from the gradient has the <strong>for</strong>mβ = D pp +2D np ξ +D nn ξ 2 , ξ = − ∂µ p/∂ρ n∂µ n /∂ρ n(4.35)where the coefficients <strong>of</strong> the gradient terms are given by D pp = D np = D nn = 132MeV·fm 5[18]. The k TF in the Coulomb interaction represents the inverse Thomas-Fermi screeninglength <strong>of</strong> the electrons. When we see the change in the sign <strong>of</strong> v, the uni<strong>for</strong>m matter phaseis more stable than the periodic structure <strong>of</strong> the nuclei. The v has a minimum atwhen q 2 min = (4πe2 /β) 1/2 −k 2 TF .v min = v 0 +2(4πe 2 β) 1/2 −βk 2 TF , (4.36)Another way to see the phase transition is to use the thermodynamic instability. Thethermodynamic stability condition can be described using the inequalities [18, 19],( ) ∂P− >0,∂vµ( (4.37)∂µ− >0.∂q c)vwhere P = P b + P e is the total pressure from electrons <strong>and</strong> baryons <strong>and</strong> µ = µ n − µ p isthe difference between the neutron <strong>and</strong> proton chemical potentials, which is the electronchemical potential in beta-stable matter. q c is defined as q c = x p −ρ e /ρ. Mathematically,the inequalities in Eq. (4.37) show that the energy per baryon is convex. Eq. (4.37) can beverified to be [18, 19]( ) [ ∂P− = ρ 2 2ρ ∂E(ρ,x (p)+ρ 2∂2 E(ρ,x p ) ∂ 2 E(ρ,x p )− ρ∂vµ∂ρ ∂ρ 2 ∂ρ∂x p( ∂µ− =∂q c)v( ∂ 2 E(ρ,x p )∂x 2 p) 2/ ] ∂ 2 E(ρ,x p )> 0,∂x 2 p) −1+ µ2 eπ 2 3 ρ > 0. (4.38)The second <strong>of</strong> Eq. (4.38) always holds, so the first will determine the phase transition in theneutron star crust. Xu et al.[18], use a simple equation to determine the instability usingthe thermodynamic relation,2∂E∂ 2 Eρ ∂ρ ∂x 2 p+ ∂2 E∂ρ 2 ∂ 2 E∂x 2 p−( ) ∂ 2 2E= ∂µ ( ) 2n ∂µ p ∂µn− . (4.39)∂ρ∂x p ∂ρ n ∂ρ p ∂ρ pEq. (4.39) is equivalent to the volume part <strong>of</strong> the thermodynamic perturbation Eq. (4.34)method except that there is a ∂µ n /∂ρ n difference. Comparing the two methods (pertubation<strong>and</strong> thermodynamic instability) shows the effects <strong>of</strong> the gradient <strong>and</strong> Coulomb terms in theperturbation method on the transition densities. Fig. 4.7 shows transition densities using65

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!