Urania-Postępy Astronomii
Wykłady z astrofizyki   
Urania-Postępy Astronomii Urania-Postępy Astronomii



Spis
treści





Spis
wykładów


Budowa i Ewolucja
Gwiazd
prof. Bohdan Paczyński
ver. 1.00


XXII. Białe karły (karły zdegenerowane).
    Białe karły są gwiazdami podtrzymywanymi przez ciśnienie zdegenerowanego gazu elektronowego, tj. energia kinetyczna kT w ich wnętrzach jest dużo mniejsza aniżeli energia Fermiego EF. Wyprowadzimy równania struktury białych karłów, nazywanych czasami karłami zdegenerowanymi, w przypadku granicznym, gdy ich ciśnienie termiczne może być pominięte, ale zdegenerowany gaz elektronowy może być albo nierelatywistycz-ny, albo nieco relatywistyczny, albo też ultrarelatywistyczny.
    Wprowadzimy zmienną x zdefiniowaną jako bezwymiarowy pęd elektronu:
x ≡ p ⁄ mc ,       xF ≡ pF ⁄ mc .
(wd.1)
Korzystając z wyprowadzeń poczynionych w rozdziale
„Równanie stanu”, gęstość i ciśnienie możemy zapisać jako funkcję bezwymiarowego pędu Fermiego xF gdzie X to obfitość wodoru jako udział masowy, μe jest średnią liczbą nukleonów na elektron, natomiast
A ≡ 
3
  (
mc
h
)3 H = 0.981 × 106        [g cm –3] ,
(wd.3a)
B ≡ 
3
  (
mc
h
)3 mc2 = 4.80 × 1023        [erg cm –3] .
(wd.3b)
    Równanie równowagi hydrostatycznej może być zapisane następująco
dP
dr
 = 
dP
dx
 
dx
dr
 = B 
x4
(1 + x2) 1 ⁄ 2
 
dx
dr
 = – 
GMr
r2
 ρ = – 
GMr
r2
 Aμex3 ,
(wd.4)
gdzie, dla uproszczenia, napisaliśmy x zamiast xF. Równanie zachowania masy, może być napisane jako
dMr
dr
 = 4π r2ρ = 4π r2ex3 .
(wd.5)
Wprowadzimy bezwymiarowe zmienne x1, x2 i x3, zdefiniowane poprzez związki:
r ≡ αrx1 ,       Mr ≡ αmx2 ,       1 + x2 ≡ x3 .
(wd.6)
Łącząc równania
(wd.4), (wd.5) oraz (wd.6), otrzymamy równania wyrażone w zmiennych bezwymiarowych
dx3
dx1
 = –  [
2GA
B
 
αm
αr
]  
x2
x12
 x3 1 ⁄ 2 ,
(wd.7)
dx2
dx1
 =  [ 4πAμe 
αr3
αm
]  x12 (x3 – 1)1.5 .
(wd.8)
Dwa parametry skalujące — αr i αm — mogą być dobrane tak, by stałe w nawiasach kwadratowych w równaniach
(wd.7) i (wd.8) równały się jedności. Osiągamy to kładąc
αr =  (
B
8πG
)1 ⁄ 2  
1
e
 = 5.455 × 108 μe– 1 (cm) = 0.00784RSłońcaμe– 1 ,
(wd.9)
αm = 
1
(2π)1 ⁄ 2
  (
B
G
)1.5  
1
(2Aμe)2
 = 2.00 × 1033 μe– 2 (g) = 1.005MSłońcaμe– 2 .
(wd.10)
Przy takich parametrach skalujących, równania
(wd.7) i (wd.8) mogą być zapisane następująco
dx3
dx1
 – 
x2
x12
 x31 ⁄ 2 ,
(wd.11)
dx2
dx1
 – x12(x3 – 1)1.5 .
(wd.12)
    Bezwymiarowe równania
(wd.11) i (wd.12) trzeba uzupełnić o dwa warunki brzegowe. Są one następujące:
x3 = x3,c ,       x2 = 0 ,       dla       x1 = 0 ,       (wewnętrzny warunek brzegowy) ,
(wd.13)
x3 = 1 ,       x2 = x2,s ,       dla       x1 = x1,s ,       (zewnętrzny warunek brzegowy) ,
(wd.14)
gdzie x3,c jest wartością w centrum zmiennej x3, a x2,s i x1,s wartościami powierzchnio-wymi bezwymiarowych odpowiednio masy i promienia. Całkowita masa i całkowity promień białego karła są dane jako
M = αmx2,s ,       R = αrx1,s .
(wd.15)
    Ponieważ mamy dwa wewnętrzne warunki brzegowe
(wd.13), jak również dwa równania różniczkowe zwyczajne (wd.11) i (wd.12), dlatego też możemy traktować (wd.13) jako warunki początkowe dla całkowań, w których x3,c jest parametrem swobodnym. Dla zadanej wartości x3,c możemy policzyć gęstość centralną korzystając z równań (wd.2a) i (wd.6). Następnie całkujemy numerycznie równania (wd.11) i (wd.12), obliczamy x1,s oraz x2,s odpowiadające x3 = 1, tj. ρ = 0 na powierzchni białego karła, by ostatecznie policzyć całkowitą masę i promień z równań (wd.15). W ten oto sposób możemy otrzymać zależność masa – promień dla białych karłów.

    Zależność masa – promień dla białych karłów może być oszacowana poprzez zastosowanie zwyczajowych przybliżeń algebraicznych do równań różniczkowych struktury gwiazdowej oraz przybliżenia analitycznego równania stanu zdegenerowanego gazu elektronowego. Równania struktury gwiazdy mogą być przybliżone następująco:
M
R
 ≈ R2ρ ,       tj.       ρ ≈ 
M
R3
,
(wd.16)
P
R
 ≈ 
GM
R2
 ρ ≈ 
GM2
R5
,       tj.       P ≈ 
GM2
R4
.
(wd.17)
Równanie stanu może być przybliżone następująco
P ≈  [  (K1ρ 5 ⁄ 3) – 2 + (K2ρ 4 ⁄ 3) – 2 ] – 1 ⁄ 2 .
(wd.18)
    Równania
(wd.16), (wd.17) i (wd.18) można połączyć, by otrzymać
P – 2 ≈ 
R8
G2M4
 ≈ 
R10
K12M 10 ⁄ 3
 + 
R8
K22M 8 ⁄ 3
.
(wd.19)
Powyższe można przearanżować do postaci
R ≈ 
K1
GM 1 ⁄ 3
  [ 1 – 
G2M 4 ⁄ 3
K22
] 1 ⁄ 2 .
(wd.20)
Ostatnie równanie powinno mieć poprawną, asymptotyczną postać, może jednakże zawierać bezwymiarowe współczynniki rzędu jedności, których nasza przybliżona analiza nie może dostarczyć. Współczynniki te możemy wszelako odzyskać zauważając, iż w dwóch przypadkach granicznych, ρ « 106 g cm – 3 i ρ » 106 g cm – 3, równanie stanu
(wd.18) jest bardzo dobrze przybliżone politropą o indeksie odpowiednio n = 1.5 i n = 3. W owych dwóch przypadkach granicznych mamy dokładne zależności masa – promień:
R = 
K1
0.4242GM 1 ⁄ 3
,       dla n = 1.5 ,
(wd.21)
M =  (
K2
0.3639G
)1.5  = 1.142 × 1034 μe – 2 (g) = 5.74MSłońcaμe – 2 ,       dla n = 3 .
(wd.22)
Łącząc równania
(wd.21) i (wd.22), możemy zapisać (wd.20) jako
R ≈ 
K1
0.4242GM 1 ⁄ 3
  [ 1 –  (
M
MCh
)4 ⁄ 3 ] 1 ⁄ 2  =
(wd.23)
0.0126RSłońca  (
2
μe
)5 ⁄ 3   (
M
MSłońca
)– 1 ⁄ 3   [ 1 –  (
M
MCh
)4 ⁄ 3 ] 1 ⁄ 2 ,
z masą Chandrasekhara
MCh = 1.435MSłońca   (
2
μe
)2 .
(wd.24)
Analityczna formuła
(wd.23) przybliża dokładną, numeryczną zależność masa – promień dla białych karłów, z błędem mniejszym niż 15% dla mas w pobliżu granicy Chandrasekhara i znacznie lepszą dokładnością dla mniejszych mas.

    Poniższa tabela zawiera porównanie numerycznych i analitycznych wartości promieni białych karłów dla μe = 2. Pierwsza kolumna to logarytm gęstości centralnej, druga — masa białego karła w jednostkach MSłońca, trzecia i czwarta — odpowiednio numeryczne i analityczne promienie białych karłów, w jednostkach RSłońca, a piąta — błąd ułamkowy promienia analitycznego.
  log ρc   M ⁄ MSłońca R ⁄ RSłońca błąd
  numeryczny     analityczny  
4 .04811 .3448 .3446 .0008
5 .14600 .2339 .2335 .0015
6 .39366 .01566 .01558 .0048
7 .80146 .01013 .00997 .0158
8 1.16176 .00619 .00593 .0411
9 1.34619 .00353 .00325 .0803
10   1.41096   .00188 .00165   .1230  



Autor: prof. Bohdan Paczyński

Tłumaczenie, opr. i wersja HTML:
Marek Gołębiewski

© "Urania-Postępy Astronomii"
webmaster: Marek Gołębiewski