You are using a browser not supported by the Wolfram Cloud
Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.
I understand and wish to continue anyway »
WOLFRAM NOTEBOOK
WDs: Lead in
WDs: Lead in
CopyRight@Zheng SHEN
E-mail: kantopenhauer@whu.edu.cn
E-mail: kantopenhauer@whu.edu.cn
In[]:=
Clear["Global`*"]
Basic Formulae
Newtonian Structure Equations
Newtonian Structure Equations
[r]
p
x
r
Gϵ
x
M
x
2
c
2
r
[r]
M
x
r
ϵ
x
2
c
EoS under Fermi Gas Model of Degenerated Electrons
EoS under Fermi Gas Model of Degenerated Electrons
Energy Density
Energy Density
n=
3
k
F
3
2
π
3
ℏ
ϵ
electron
4
m
e
5
c
2
π
3
ℏ
k
F
m
e
∫
0
1/2
(+1)
2
u
2
u
ϵ
total
ϵ
electron
nA
m
N
2
c
Z
ϵ
electron
m
N
Z
3
m
e
5
c
3
2
π
3
ℏ
3
k
F
3
m
e
3
c
[Energy density is dominated by baryon rest mass term]
[Energy density is dominated by baryon rest mass term]
it can be shown, for WDs composed of degenerated electrons and baryons, the total energy density is always dominated my the rest mass term (1)the factor in front of the integral in is about , while the integral itself is obviously dimensionless
nA
m
N
2
c
Z
ϵ
electron
1.4*10^23J/
3
m
In[]:=
Module{m=9.1093837015*,c=3*,ℏ=1.054517871*},
-31
10
8
10
-34
10
4
m
5
c
2
π
3
ℏ
Out[]=
1.44577×
23
10
(2) the integral take values as plotted if upper limit of integral set from 0 to 2 (when taking upper limit 2, the Fermi momentum is sufficiently relativistic
In[]:=
Plotu,{v,0,2}
v
∫
0
1/2
(+1)
2
u
2
u
Out[]=
(3)while the scale of (the factor with dimension of energy density) is about , is over 100 times larger than degenerated electrons’ contribution
m
N
3
m
e
3
c
3
2
π
3
ℏ
8.8×J/
25
10
3
m
In[]:=
Module{m=9.1093837015*,c=3*,ℏ=1.054517871*,M=1.67*},M
-31
10
8
10
-34
10
-27
10
3
m
5
c
3
2
π
3
ℏ
Out[]=
8.83501×
25
10
[#Domination? A More Precise Evaluation]
[#Domination? A More Precise Evaluation]
The portion of degenerated-electron energy density in total energy density, i.e. + when ctakes value from 0 to 100
ϵ
electron
ϵ
electron
ϵ
baryon
k
F
m
e
In[]:=
Module{m=9.1093837015*,c=3*,ℏ=1.054517871*,M=1.67*,A=56,Z=26},Plotu+u,{v,0,10000}
-31
10
8
10
-34
10
-27
10
4
m
5
c
2
π
3
ℏ
v
∫
0
1/2
(+1)
2
u
2
u
MA
3
m
5
c
3Z
2
π
3
ℏ
3
v
4
m
5
c
2
π
3
ℏ
v
∫
0
1/2
(+1)
2
u
2
u
Out[]=
In[]:=
Module{m=9.1093837015*,c=3*,ℏ=1.054517871*,M=1.67*,A=56,Z=26},Plotu+u,{v,0,1000}
-31
10
8
10
-34
10
-27
10
4
m
5
c
2
π
3
ℏ
v
∫
0
1/2
(+1)
2
u
2
u
MA
3
m
5
c
3Z
2
π
3
ℏ
3
v
4
m
5
c
2
π
3
ℏ
v
∫
0
1/2
(+1)
2
u
2
u
Out[]=
Thus, is relevantly small, compared with the contribution of rest mass of baryons, even when Fermi momentum takes extremely relativistic values.But we shall still be aware that when ctakes values over this approximation will be no longer practical.
ϵ
electron
k
F
m
e
3
10
Pressure
Pressure
p=cu
4
m
e
5
c
3
2
π
3
ℏ
k
F
m
e
∫
0
-1/2
(+1)
2
u
4
u
[ParametricPlot: a Panoramic on p~ϵ Relation]
[ParametricPlot: a Panoramic on p~ϵ Relation]
Integral p and ϵ numerically with Fermi momentum as parameter, we can estimate the relation between p and ϵ under the limit ; by observation the curve has a shape like
k
F
0<c<2
k
F
m
e
p~
4/3
ϵ
In[]:=
Module{m=9.1093837015*,c=3*,ℏ=1.054517871*,M=1.67*,Z=26,A=56},ParametricPlotM,u,{v,0,2}
-31
10
8
10
-34
10
-27
10
A
Z
3
m
5
c
3
2
π
3
ℏ
3
v
4
m
5
c
3
2
π
3
ℏ
v
∫
0
-1/2
(+1)
2
u
4
u
Out[]=
Polytrope EoS
Extreme Relativistic ( kF>>me)
Extreme Relativistic ( >>)
k
F
m
e
Pressure:[prel≈c4kF122π3ℏ=4me5c122π3ℏ4kFmec]
Pressure:[≈=c]
p
rel
c
4
k
F
12
2
π
3
ℏ
4
m
e
5
c
12
2
π
3
ℏ
4
k
F
m
e
Dimensionless Equations
Dimensionless Equations
[Failed Attempt]
[Failed Attempt]
Dimensionless Equations
Dimensionless Equations
Numerical Integration
Numerical Integration
the result of R is unexpectedly small, that is because the choice of initial condition of pressure p is not realistic
[[another attempt]]
[[another attempt]]
[[Another Attempt: Try Something New]]
[[Another Attempt: Try Something New]]
For WDs, once EoS has been determined, pressure~radius and mass~radius relation depend only on the choice of initial condition (central pressure). Thus in this subsubsection I use Plot3D function to show the pressure~(radius, central pressure) and mass~(radius, central pressure) dependence graphically.
*The results are beyond my expectation ( I didn’t expect MMA can do this) but still primitive. Since the results are of no importance in our research into NS, the adjustment of arguments is postponed until we deal with NS models.
**what inspired me to use Plot3D is that peaks of some diagram might have some relationship with the behavior of system under disturbance, and this might be represented more explicitly on 3D diagrams than 2Ds.
**what inspired me to use Plot3D is that peaks of some diagram might have some relationship with the behavior of system under disturbance, and this might be represented more explicitly on 3D diagrams than 2Ds.
Dimensionless Equations
Dimensionless Equations
Dimensionless Equations
Dimensionless Equations
Numerical Integration
Numerical Integration
[realistic initial condition]
[realistic initial condition]
Rel vs NR Approximations: Conclusions
Rel vs NR Approximations: Conclusions