Solving, we haveα = 9 5 + 5B + 1 ( 3T o ǫ 5 + 5B − K )o= 18 T o 3T o 5 + 20B − K o,T o T oβ ′′ = 5K/T o −3−9ǫ−45(1+ǫ)B/T o10−15ǫη = 18−10K/T o +150B/T o15ǫ(2−3ǫ)= K oT o− 6 5 − 12BT o,= 18 5 + 30BT o− 2K T o,K o ′ = 143 K o −15B − 3 (5 T o +ǫ K o −15B − 9 )5 T o(3.20)= 5K o −20B − 6 5 T o.The last expressions are obtained with ǫ = 1/3. We also find, <strong>for</strong> ǫ = 1/3,[ 1S v = T o3 + α 2 −α L + 5 3 β′′ L − 5 9 β′′ +η L − η ],2[ 2S v ′ = T o9 + α 2 −α L + 259 β′′ L − 2527 β′′ + 4 3 η L − 2 ]3 η .Pure neutron matter has pressure[ 2p N = T o ρ o u5 (2u)2/3 −α L u+ 5 ]3 22/3 β Lu ′′ 5/3 +(1+ǫ)η L u 4/3(3.21)(3.22)(3.23)At ρ o , p N has a value about half <strong>of</strong> the non-interacting Fermi pressure. Combining this resultwith reasonable values S v = 30 MeV <strong>and</strong> S v ′ = 15 MeV, <strong>and</strong> assuming ǫ = 1/3, one findsα L ≃ −4.4, β L ′′ ≃ 3.6, <strong>and</strong> η L ≃ −10.7. In the case S v ′ = 30 MeV, one finds α L ≃ −8.4,β L ′′ ≃ 6.7, <strong>and</strong> η L ≃ −19.9. In both cases, since β L ′′ > 0, the neutron matter pressure willcontinuously rise with u.However, we now observe that[ ( 3Q = T o5 u2/3 2 2/3 − 14 ) [ β′′+(29 18 + 2/3 − 5 ) ]β L]u ′′ 5/3 , (3.24)3evaluated at the saturation density <strong>and</strong> using ǫ = 1/3, S v = 30 MeV <strong>and</strong> S v ′ = 15 MeV, isQ 1 ≃ −10.3 MeV. Obviously, the neutron matter properties <strong>and</strong> the symmetry propertiescannot be separately adjusted.Most Skyrme <strong>for</strong>ces have the properties that both Q 1 <strong>and</strong> Q ′ 1 ≡ (dQ(u)/dlnu) 1 = [p n −p s −S v ′ ] have small magnitudes: if they don’t, their neutron matter has problems. Setting44
Q 1 to a small number also results in small values <strong>for</strong> Q ′ 1. since <strong>for</strong> ǫ = 1/3, the value <strong>of</strong> β ′′is nearly zero, setting β L ′′ = 0 will ensure that Q 1 is small. The coupled set <strong>of</strong> equations todetermine α L <strong>and</strong> η L in this case areThese result inS vT o= 1 3 + α 2 − η 2 −α L +η L ,S ′ vT o= 2 9 + α 2 − 2 3 η −α L + 4 3 η L.α L = 2 3 + α 2 − η 2 − 4S v −3S ′ vT o≃ 1.6,η L = 1 3 + η 2 −3S v −S ′ vT o≃ 1.1.(3.25)(3.26)The neutron matter energy <strong>and</strong> pressure are now implicitly closely coupled to S v <strong>and</strong> S ′ vsuch that e N,1 ≃ S v −B <strong>and</strong> p N,1 /ρ o ≃ S ′ v . For Hebler & Schwen [63] estimates <strong>of</strong> e N,1 ≃ 14MeV <strong>and</strong> p N,1 /ρ o ≃ 2 MeV, we find S v ≃ 30 MeV <strong>and</strong> S ′ v ≃ 10 MeV. In this case, α L ≃ 1.2<strong>and</strong> η L ≃ 0.7.Alternatively, one could use a value <strong>of</strong> ǫ < 1/3, resulting in β ′′ > 0. To ensure Q 1 ≡ 0,a positive value <strong>for</strong> β L ′′ is also obtained. This results in different values <strong>for</strong> α L <strong>and</strong> η L <strong>and</strong> adifferent relation between S v , S v ′ <strong>and</strong> e N,1, p N,1 .Thus, it would be interesting to fit nuclei to determine the correlation between S v <strong>and</strong> S v ′ <strong>for</strong>a couple different values <strong>of</strong> ǫ, such as 1/3 or 1/6. One could then also add the correlations<strong>for</strong> S v <strong>and</strong> S v ′ obtained from the relations reduced from theoretical studies <strong>of</strong> neutron matterenergy <strong>and</strong> pressure at the saturation density.3.3 Optimized Parameter SetAddingnewdensitydependentinteractionssolvesnuclearincompressibility, opticalpotential,<strong>and</strong> pure neutron matter properties. To find an optimized parameter set, α L,U , β L,U , <strong>and</strong>,η L,U , we may compare the experimental binding energy <strong>of</strong> single nucleus with theoreticalcalculation from the modified model. For this set, we fix β ′′ = 0, which gives proper nuclearincompressibility (K ≃ 235 MeV) <strong>and</strong> varies S v <strong>and</strong> S v ′ (= L/3). Since the experimentalerror in the binding energy is extremely small, so that it is meaningless if we define χ 2 inconventional way, instead we define the χ 2 asχ 2 = 1 N∑ ( ) 2B ex,i −B th,i(3.27)σwhere N(= 2336) is the total number <strong>of</strong> nuclei in the calculation <strong>and</strong> σ is a mean error(arbitrary). The B ex,i is the data [45] <strong>and</strong> B th,i is the numerical calculation.The χσ has a unit <strong>of</strong> MeV so we can estimate the mean difference between B ex <strong>and</strong> B th ineach calculation. Even though our code is fast to calculate single nucleus’s properties, thenumber <strong>of</strong> calculation <strong>for</strong> this contour plot is 2336×41×41 ∼ 4·10 6 so we wrote a simple45
- Page 1 and 2:
Theory of nuclear matter of neutron
- Page 3 and 4:
Abstract of the DissertationTheory
- Page 5 and 6:
To Jaenyeong, Jonghyun, and Jongbum
- Page 7 and 8:
3.3 Optimized Parameter Set . . . .
- Page 9 and 10: List of Figures1.1 Mass and radius
- Page 11 and 12: SymbolsNuclear Physics SymbolE ener
- Page 13 and 14: ω δt 90−10QtT cTaylor expansion
- Page 15 and 16: List of Tables1.1 Range of Tables .
- Page 17 and 18: Chapter 1IntroductionStars give us
- Page 19 and 20: Since the mass of neutron star is d
- Page 21 and 22: shown in Fig. 1.2, LS220 is the onl
- Page 23 and 24: density:p = ρ2ρ o[ K9µ n = −B
- Page 25 and 26: nuclear energy density functional i
- Page 27 and 28: 2.3.2 Thermodynamic QuantitiesThe t
- Page 29 and 30: point r and with momentum p isU n (
- Page 31 and 32: Thus, in the case of zero temperatu
- Page 33 and 34: For standard nuclear matter, u = 1
- Page 35 and 36: a α L α U β L β U σ0.5882 1.11
- Page 37 and 38: Figure 2.2: Bulk equilibrium of T =
- Page 39 and 40: Figure 2.4: Coexistence curves for
- Page 41 and 42: proton fraction x are irrelevant fo
- Page 43 and 44: potential model analogue for the Ha
- Page 45 and 46: with the convention that z II is th
- Page 47 and 48: Figure 2.10: The surface tension of
- Page 49 and 50: Figure 2.12: Contour plots of surfa
- Page 51 and 52: 2.6.2 Dense Matter and the Wigner-S
- Page 53 and 54: Figure 2.15: Proton number per unit
- Page 55 and 56: Chapter 3Modified model† The trun
- Page 57 and 58: We choose the value ǫ = 1/3. One m
- Page 59: where∆V t = T oρ o[( ρρ o) ǫ[
- Page 63 and 64: where p = (S v ,L) and M ij = 1 ∂
- Page 65 and 66: 3.4 Nuclear Surface TensionCompared
- Page 67 and 68: 1.41.2FRTF II ω 0 , T=0 MeVFRTF II
- Page 69 and 70: For the finite-range term, we use a
- Page 71 and 72: Landau’s quasi-particle formula g
- Page 73 and 74: The symmetry energy in nuclear matt
- Page 75 and 76: Figure 4.2: The left figure shows t
- Page 77 and 78: Figure 4.4: This figure shows the b
- Page 79 and 80: Figure 4.5: In a neutron star, heav
- Page 81 and 82: The energy exchange from the gradie
- Page 83 and 84: Figure 4.8: The left panel shows th
- Page 85 and 86: 4.5 Astrophysical application4.5.1
- Page 87 and 88: 0.550.5GaussianFRTF IIFRTF I0.450.4
- Page 89 and 90: neutron matter respectively. Those
- Page 91 and 92: Table 5.3: Nuclear properties of si
- Page 93 and 94: which aref Rc = a 0 +a 1 x+a 2 x 2
- Page 95 and 96: parametrized[ht]504540EDF, Z nucFit
- Page 97 and 98: Chapter 6Nuclear Equation of State
- Page 99 and 100: Coulomb, and, symmetry part [24], i
- Page 101 and 102: We can eliminate N s , leading to t
- Page 103 and 104: pressure from pure neutron matter,
- Page 105 and 106: 2.52GSRsSGISIVSkaSkI1SkI2SkI3SkI5St
- Page 107 and 108: HNsn, p, α,eFigure 6.7: In the Wig
- Page 109 and 110: Table 6.1: Surface tension analytic
- Page 111 and 112:
free energy density for the alpha p
- Page 113 and 114:
Here s ′ = ∂s(u)/∂u.The minim
- Page 115 and 116:
6.5.2 Determination of Coulomb surf
- Page 117 and 118:
gives analytic derivatives of therm
- Page 119 and 120:
200180160Atomic number, Y p =0.45,
- Page 121 and 122:
for Fermi integral. This fitting fu
- Page 123 and 124:
A.1 Taylor expansion integrationThe
- Page 125 and 126:
thenwhere0u − i+1 = fu− i +J
- Page 127 and 128:
We separate the two integrals as in
- Page 129 and 130:
Each w i can be obtained[ae −∆/
- Page 131 and 132:
By expanding e r 0/a −e −r 0/a
- Page 133 and 134:
M 3 2a 0.433 1p 0 (e 2 /a) √ π/3
- Page 135 and 136:
of find exact polynomial expression
- Page 137 and 138:
Appendix CPhase coexistenceBulk equ
- Page 139 and 140:
To summarize, we need to solve 5 eq
- Page 141 and 142:
D.1 Non-uniform electron density ap
- Page 143 and 144:
Appendix ENuclear Quantities in Non
- Page 145 and 146:
Table E.1: Non-relativistic Skyrme
- Page 147 and 148:
[20] C.J. Pethick and D.G. Ravenhal