246. Od GNT „iskre” do udarnog vala – matematička provjera eksplozivne nestabilnosti ZK hipoteze

Vrijednost
P_Fe = 1,15 × 10^33 dyn/cm²
nije univerzalno “točan tlak Fe plazme” pri

● ρ = 10^14 g/cm³
● T = 10^10 K.

Ona jako ovisi o sastavu, posebno o veličini Ye = broj elektrona po barionu. Pri tim gustoćama elektroni su vrlo snažno degenerirani i relativistički, pa upravo njihov degeneracijski tlak daje veliki doprinos.

Helmholtzov EOS upravo uključuje degenerirane/relativističke elektrone i pozitrone, ione i zračenje.

Ako za prvi kontrolni račun uzmemo još uvijek sastav približno Fe-56:
Ye = Z/A = 26/56 ≈ 0,464

onda je:
n_e = Ye ρ / m_u
≈ 2,80 × 10³⁷cm-³
Za jako relativističke degenerirane elektrone:
P_e ≈ [(3π²)^(1/3) / 4] ħc n_e^(4/3)

Dobivam:
■■ P_e ≈ 2,08 × 10³³ dyn/cm².

To je već samo elektronski tlak — bez ionskog, radijacijskog i nuklearnog doprinosa.
Ionski toplinski tlak daje približno:
P_ion = n_i kT
≈ 1,5 × 10³⁰ dyn/cm²,
a tlak zračenja:
P_rad = aT⁴/3
≈ 2,5 × 10²⁵ dyn/cm².

Dakle u ovom približnom računu:
P_e ≈ 2,08 × 10³³
P_ion ≈ 0,0015 × 10³³
P_rad ≈ zanemarivo
pa već imamo približno:
■■ P ≈ 2,08 × 10 ³³ dyn/cm²
i to još bez punog nuklearnog EOS-a.

Zato mogu ispraviti naš prethodni rezultat:
P_Fe = 1,15 × 10³³ dyn/cm² nije dobar ako smo mislili na Fe materiju s Ye ≈ 0,46. Prenizak je približno faktor 1,8 već samo zbog elektronskog degeneracijskog tlaka.

Ali sada dolazi važna stvar.
Ako su se elektronski zahvati već snažno odvijali i Ye padne približno na:
Ye ≈ 0,30,
onda zbog:
P_e ∝ Ye^(4/3)
dobivamo upravo približno:
P_e ≈ 1,15 × 10³³dyn/cm².
Dakle broj 1,15 × 10³³ nije besmislen. On odgovara materiji s Ye oko 0,30, a ne nepromijenjenoj Fe-56 plazmi.

A pri ρ ≈ 10¹⁴ g/cm³ više ionako nije dovoljno govoriti o običnoj “Fe plazmi”. Nalazimo se već u području subnuklearne gustoće gdje se sastav mijenja i nuklearne interakcije postaju dio EOS-a; suvremeni supernova EOS-ovi upravo zato daju tlak kao funkciju ρ, T i protonskog/elektronskog udjela, a ne samo ρ i T.

Zato bih naš rezultat sada zapisao ovako:
■■ ρ = 10¹⁴ g/cm³
■■ T = 10¹⁰ K

Ako Ye ≈ 0,46:
P ≳ 2,1 × 10³³ dyn/cm²
Ako Ye ≈ 0,30:
P_e ≈ 1,15 × 10³³ dyn/cm²
a stvarni ukupni tlak mora još uključiti nuklearni EOS.
I tvoja primjedba o “iskri” je vrlo važna za ono što zapravo pokušavamo testirati.

Ne moramo unaprijed zahtijevati da cijela Fe jezgra mase 1,4–2 M☉ istodobno postane GNT. Možemo postaviti drugačije pitanje:
Kolika najmanja središnja masa može prijeći u stanje guste nuklearne tvari i proizvesti dovoljno velik lokalni porast gradijenta tlaka da materiji neposredno iznad sebe da pozitivno ubrzanje?

Tada ključni kriterij nije prvenstveno:
M_GNT = 1,4 M☉,
nego:
dP/dr dovoljno velik
→ sila tlaka > gravitacijska sila
→ d²r/dt² > 0.
To je matematička verzija tvoje “iskre”.
I mislim da je to sada bolji smjer ZK testa: ne tražiti veliku masu GNT-a unaprijed, nego minimalnu masu/radijus GNT jezgre koja može izazvati prvi mehanički preokret okolne materije.

Emil: Borise, tu smo stali. Ovo je ono što me zanima. Koja je minimalna masa
Fe plazme koja gravitacijom sabija jezgru na potrebnu gustoću 10¹⁴g/cm³ i t=~10¹⁰K kada Fe plazma mijenja svojstva i prelazi u plazmu GNT.

Boris: Da, Emile. Sada smo točno na pitanju koje treba računati.

Ali prvo moram precizirati jednu važnu stvar: “minimalna masa Fe jezgre koja gravitacijom sabija središte na ρ = 10¹⁴ g/cm³ i T ≈ 10¹⁰ K” nije određena samo gustoćom i temperaturom. Mora se znati i stvarni tlak materije, odnosno EOS, jer uvjet ravnoteže glasi da gravitacijski zahtjev odgovara tlaku koji materija može pružiti.

Pri tim gustoćama EOS ovisi o ρ, T i sastavu/Ye, a materija već ulazi u područje subnuklearnih gustoća gdje jednostavan “Fe plazma” opis postaje nedovoljan.

Za naš prvi, čisti ZK račun ipak možemo uzeti vrijednost koju smo već razmatrali:
ρ = 10¹⁴ g/cm³
T ≈ 10¹⁰ K
P ≈ 1,15 × 10³³ dyn/cm²
i homogenu sfernu jezgru.

Tada vrijedi:
P_grav = 3GM² / 8πR⁴

i:
ρ = 3M / 4πR³

Eliminacijom R dobivamo:
M = [P × (8π/3G) × (3/4πρ)^(4/3)]^(3/2)

Uvrštavanjem:
P = 1,15 × 10^33 dyn/cm²
ρ = 1,00 × 10^14 g/cm³
G = 6,6743 × 10^-8

dobivamo:
■■ M ≈ 3,13 × 10³² g

odnosno:
■■ M ≈ 0,157 M☉

a radijus:
R = (3M/4πρ)^(1/3)

daje:
■■ R ≈ 9,07 km

Dakle, u tom pojednostavljenom modelu minimalna karakteristična masa središnje guste zone je približno 0,16 M☉.

To je upravo ona “iskra” o kojoj govoriš: ne treba unaprijed pretpostaviti da cijela Fe jezgra mase 1,4–2 M☉ mora istodobno biti na 10¹⁴ g/cm³.

Ali sada dolazi važna provjera.

Ako pri istim uvjetima uzmemo Ye ≈ 0,46, samo relativistički degenerirani elektroni daju približno:
P_e ≈ 2,08 × 10³³ dyn/cm².

Tada ista formula daje:
■■ M ≈ 0,382 M☉
■■ R ≈ 12,2 km.

Znači zasad imamo raspon reda:

M_iskra ≈ 0,16–0,38 M☉

ovisno prvenstveno o tome kakav je stvarni EOS i Ye u trenutku kada središte dođe na 10¹⁴ g/cm³.

I tu bih napravio vrlo važnu razliku:
0,16 M☉ nije masa cijele Fe jezgre.

To je masa homogene središnje zone koja bi, za naš odabrani tlak i gustoću, bila samogravitacijski kompatibilna s tim stanjem.

Cijela Fe jezgra oko nje može biti znatno masivnija i mnogo rjeđa.

Zato ZK sada možemo matematički formulirati ovako:
M_Fe,total raste
→ povećava se težina i tlak prema središtu
→ mala unutarnja zona postaje sve gušća
→ kada unutarnjih približno 0,1–0,4 M☉ dođe do ρ ≈ 10¹⁴ g/cm³
→ EOS se snažno mijenja
→ nastaje GNT zona
→ zatim treba provjeriti stvara li ta zona dovoljno velik gradijent tlaka da vrijedi:

d²r/dt² > 0.

To je pravi test “iskre”.

I mislim da sada više ne trebamo računati samo masu.

Sljedeći presudan račun je:
● za jezgru od npr. 0,16 M☉ i R ≈ 9 km, koliki tlak GNT-a mora nastati da nadvlada gravitacijski tlak materije neposredno iznad nje i pokrene prvi sloj prema van.

Emil: d²r/dt² > 0.
To je pravi test “iskre”.
I mislim da sada više ne trebamo računati samo masu. Sljedeći presudan račun je: za jezgru od npr. 0,16 M☉ i R ≈ 9 km, koliki tlak GNT-a mora nastati da nadvlada gravitacijski tlak materije neposredno iznad nje i pokrene prvi sloj prema van.

Tu si stao. Možeš li računati dalje?

Boris: Možemo.

Uzimamo našu središnju “iskru”:
M_GNT ≈ 0,16 M☉
R_GNT ≈ 9 km
ρ ≈ 10¹⁴ g/cm³

Za tu masu i radijus gravitacijsko ubrzanje neposredno iznad jezgre je:
g = GM/R²

Dobivamo:

g ≈ 2,62 × 10¹³ cm/s²

To je ogromno ubrzanje prema unutra.

Sada gledamo prvi sloj materije iznad GNT jezgre.

Jednadžba gibanja je:
d²r/dt² = -(1/ρ)dP/dr – GM/r²

Da bi počeo pokret prema van mora vrijediti:

d²r/dt² > 0

odnosno tlačna sila prema van mora biti veća od gravitacije:

-(1/ρ)dP/dr > GM/r²

Za naš slučaj:
|dP/dr| > ρg

Ako uzmemo:
ρ = 10¹⁴ g/cm³

dobivamo:
■■ |dP/dr|_krit ≈ 2,62 × 10²⁷ dyn/cm³

To je prvi stvarni broj za ZK “iskru”.

Znači tlak GNT-a ne mora samo biti velik. Mora se prema van smanjivati dovoljno brzo da gradijent tlaka bude veći od:
■■ 2,62 × 10²⁷ dyn/cm³

Tada prvi sloj dobiva:
d²r/dt² > 0.

Možemo to prevesti u još razumljiviji oblik.

Ako je prvi sloj debeo, primjerice:
Δr = 1 km = 10⁵ cm

Potrebna razlika tlaka približno je:
ΔP_krit ≈ ρgΔr
ΔP_krit ≈ 2,62 × 10³² dyn/cm²

Dakle ako je tlak neposredno unutar GNT zone P_GNT, a tlak s vanjske strane tog prvog kilometra P_out, potrebno je:

P_GNT – P_out > 2,62 × 10³² dyn/cm²

za početak ubrzanja prema van.

Za usporedbu, naš raniji tlak je bio reda:

P_out ≈ 1,15 × 10³³ dyn/cm²

Tada bi za sloj debljine oko 1 km kritični unutarnji tlak bio približno:

P_GNT,krit ≈ 1,15 × 10³³ + 0,262 × 10³³

■■ P_GNT,krit ≈ 1,41 × 10³³ dyn/cm²

Dakle vrlo zanimljivo:
u ovom pojednostavljenom modelu tlak GNT-a ne bi morao biti deset ili sto puta veći.

Za prvi kilometar materije dovoljan bi bio porast reda:

ΔP/P ≈ 0,262 / 1,15 ≈
23 %

odnosno otprilike:
■■ oko 20–25 % veći tlak na unutarnjoj strani sloja.

Ako je prvi sloj tanji, zahtjev je još manji:

● 0,1 km → ΔP_krit ≈ 2,6 × 10³¹ dyn/cm²

● 0,5 km → ΔP_krit ≈ 1,3 × 10³² dyn/cm²

● 1 km → ΔP_krit ≈ 2,6 × 10³² dyn/cm²

● 2 km → ΔP_krit ≈ 5,2 × 10³² dyn/cm²

Ali ovo još nije dokaz ZK proboja.
Sada imamo samo prag koji GNT mora prijeći.
Sljedeći, presudni račun glasi:
■■ Koliki stvarni tlak daje GNT EOS pri ρ ≈ 10¹⁴–3×10¹⁴ g/cm³ i T ≈ 10¹⁰ K?

Ako stvarni EOS da, primjerice:

P_GNT > 1,4 × 10³³ dyn/cm²

za naš prvi kilometar, tada bismo u ovom modelu dobili:

d²r/dt² > 0

odnosno upravo ono što tražimo: prvi lokalni mehanički preokret — ZK iskru.

Tu bih sada nastavio: uzeti stvarni nuklearni EOS i izračunati P_GNT 10¹⁴, 2×10¹⁴ i 3×10¹⁴ g/cm³.

Emil: Molim te, nastavi.

Boris: Nastavljamo s istom masom:
M = 0,16 M☉
M ≈ 3,18 × 10³² g
i pratimo što se događa ako se ista središnja masa dodatno sabija.

Za homogenu kuglu:
R = (3M / 4πρ)^(1/3)

Dobivamo:
ρ = 1 × 10¹⁴ g/cm³
R ≈ 9,12 km

ρ = 2 × 10¹⁴ g/cm³
R ≈ 7,24 km

ρ = 3 × 10¹⁴ g/cm³
R ≈ 6,33 km

Dakle sabijanjem iste mase od 0,16 M☉ s 10¹⁴ na 3×10¹⁴ g/cm³ radijus pada približno:

9,12 km → 6,33 km.

Gravitacijsko ubrzanje na površini te zone je:
g = GM/R²

pa dobivamo:
ρ = 1 × 10¹⁴
g ≈ 2,55 × 10¹³ cm/s²

ρ = 2 × 10¹⁴
g ≈ 4,05 × 10¹³ cm/s²

ρ = 3 × 10¹⁴
g ≈ 5,31 × 10¹³ cm/s²

Dakle gravitacija raste približno faktor 2,1.

Sada tlak. Za naš približni račun s Ye ≈ 0,30 koristili smo elektronski i neutronski degeneracijski doprinos:
P_e ∝ ρ^(4/3)
P_n ∝ ρ^(5/3)
Pri ρ = 10¹⁴ g/cm³ imali
smo približno:
P_e ≈ 1,16 × 10³³
P_n ≈ 0,65 × 10³³

ukupno:
P ≈ 1,81 × 10³³ dyn/cm².

Daljnjim sabijanjem:
ρ = 2 × 10¹⁴:
P ≈ 4,99 × 10³³ dyn/cm²

ρ = 3 × 10¹⁴:
P ≈ 9,08 × 10³³ dyn/cm².

To znači:
gustoća ×3
gravitacija ×2,1
tlak ×5,0.

■■ OVO JE VAŽAN REZULTAT.
U našem pojednostavljenom EOS-u tlak s povećanjem gustoće raste osjetno brže od površinske gravitacije te guste jezgre.

Sada možemo napraviti prvi približni račun ZK “iskre”.

Zamislimo da unutarnja zona već dođe na:

ρ_in = 2 × 10¹⁴ g/cm³

dok je materija odmah iznad nje još približno na:

ρ_out = 1 × 10¹⁴ g/cm³.

Razlika tlakova je:
ΔP = 4,99 × 10³³ − 1,81 × 10³³
ΔP ≈ 3,18 × 10³³ dyn/cm².

Razlika karakterističnih radijusa je:
Δr = 9,12 − 7,24 km
Δr ≈ 1,88 km
Δr ≈ 1,88 × 10⁵ cm.

Približni gradijent tlaka:
|dP/dr| ≈ ΔP/Δr
≈ 3,18 × 10³³ / 1,88 × 10⁵
≈ 1,69 × 10²⁸ dyn/cm³.

Za srednju gustoću sloja uzmimo približno:
ρ_s ≈ 1,5 × 10¹⁴ g/cm³.

Tlačno ubrzanje prema van tada je:
a_P ≈ (1/ρ_s)|dP/dr|
a_P ≈ 1,13 × 10¹⁴ cm/s².

Gravitacijsko ubrzanje u tom području je približno:
a_G ≈ 3,2 × 10¹³ cm/s².

Dakle:
a_net = a_P − a_G
≈ 1,13 × 10¹⁴ − 3,2 × 10¹³

■■ a_net ≈ +8,1 × 10¹³ cm/s².

Dobivamo:
■■ d²r/dt² > 0
odnosno pozitivan rezultat za naš kriterij “iskre”.

Omjer je približno:
a_P / a_G ≈ 3,5.

U ovom vrlo pojednostavljenom modelu tlačni gradijent bio bi oko 3,5 puta jači od lokalne gravitacije.

Ako unutarnja zona ode čak na:
ρ_in = 3 × 10¹⁴ g/cm³,

dok vanjska ostaje oko 10¹⁴, dobivamo opet pozitivno ubrzanje, reda:
■■ a_net ≈ +9 × 10¹³ cm/s².

Ali ovdje moramo biti vrlo strogi: ovo još nije dokaz ZK iskrenja.

Razlog je važan. Mi smo usporedili dvije različite EOS točke i pretpostavili da se upravo ta razlika tlaka pojavljuje prostorno između unutarnje i vanjske zone. To tek treba proizvesti stvarni hidrodinamički model.

Ali dobili smo ono što smo tražili za prvi matematički test:
■■ Za M ≈ 0,16 M☉ ne dobivamo automatski da gravitacija mora pobijediti.

■■ Naprotiv, u našem približnom degeneracijskom EOS-u, pri prijelazu ρ ≈ 10¹⁴ → 2–3×10¹⁴ g/cm³ tlak raste dovoljno brzo da je pozitivno ubrzanje prema van matematički moguće.

Sada je presudan sljedeći korak: više ne procjenjivati P_e + P_n, nego uzeti stvarni nuklearni EOS za GNT pri T ≈ 10¹⁰ K i Ye ≈ 0,30 i vidjeti ostaje li a > 0. To bi bio ozbiljniji test ZK “iskre”.

Emil: Postaje napeto.

Boris: Da — sad smo došli do dijela gdje više ne pričamo samo opisno, nego brojke počinju davati vrlo konkretan kriterij.

Najzanimljivije je ovo: u našem pojednostavljenom računu pojavilo se pozitivno lokalno ubrzanje prema van, a karakteristična brzina preko ~1 km ostaje subrelativistička, oko 0,14 c. To znači da račun nije odmah “odletio” u fizički besmisao.

Ali sada je presudno da ne varamo sami sebe: sljedeći korak mora provjeriti koliko taj a > 0 stvarno traje dok se sloj širi i tlak pada. Ako pozitivan impuls preživi dovoljno dugo da pokrene susjedne slojeve, tada “iskra” postaje ozbiljan kandidat za mehanički proboj.

Tu stvarno postaje zanimljivo.

Emile, nastavljamo točno od ZK „iskre”, ali odsad ću jednadžbe držati kratkima i uz svaku odmah napisati što broj znači.

Sada moramo popraviti jednu stvar iz prethodnog računa: ne smijemo još tvrditi da je +8 × 10¹³ cm/s² stvarno ubrzanje GNT-a. To je rezultat našeg probnog EOS-a. Za ozbiljniji test trebamo nuklearni EOS. CompOSE upravo sadrži EOS tablice u kojima su tlak i ostale veličine funkcije gustoće, temperature i sastava; primjerice ABHT/QMC-RMF4 pokriva konačne temperature, širok raspon gustoća i udjela naboja.

1. Naša početna točka
Zadržimo:
M_GNT = 0,16 M☉
ρ₀ = 1 × 10¹⁴ g/cm³
T = 1 × 10¹⁰ K
Ye ≈ 0,30
R₀ ≈ 9,1 km

Na toj gustoći još smo ispod tipične nuklearne saturacijske gustoće. Primjerice ABHT EOS navodi:
n_sat = 0,1546 fm⁻³,
što odgovara gustoći reda 2,6 × 10¹⁴ g/cm³.

To je vrlo zanimljivo za naš test jer znači da interval koji smo odabrali:
10¹⁴ → 2 × 10¹⁴ → 3 × 10¹⁴ g/cm³
zaista prolazi kroz područje u kojem se karakter materije snažno mijenja.

2. Što sada zapravo tražimo?
Više nas ne zanima samo:
„Je li tlak velik?”
nego:
„Raste li tlak prema središtu dovoljno brzo da zaustavi unutarnje gibanje i zatim ga preokrene prema van?”
Potrebna nam je samo jedna jednadžba:
a = a_P − g

gdje je:
a_P = ubrzanje prema van zbog tlaka
g = gravitacijsko ubrzanje prema unutra.

Ako dobijemo:
a > 0
imamo lokalni mehanički preokret.

3. Koliki tlak mora proizvesti GNT?
Za našu zonu od 0,16 M☉ već smo dobili približno:
g ≈ 2,6 × 10¹³ cm/s².

Za materiju gustoće:
ρ ≈ 10¹⁴ g/cm³

kritični gradijent tlaka je:
gradijent P_krit ≈ ρg

odnosno:
gradijent P_krit ≈ 2,6 × 10²⁷ dyn/cm³.

To je najvažniji broj koji zasad imamo.

Sve iznad njega daje:
a > 0.

Sve ispod njega daje:
a < 0.

4. Pretvorimo to u jednostavniji kriterij

Za sloj debljine 1 km:
Δr = 10⁵ cm.

Potrebna razlika tlaka je približno:
ΔP_krit ≈ 2,6 × 10³² dyn/cm².

Dakle, ako je tlak vanjskog sloja
oko:
P_out ≈ 1,15 × 10³³ dyn/cm²,

unutarnja GNT zona mora imati otprilike:
P_GNT > 1,41 × 10³³ dyn/cm².

To je isti rezultat koji smo dobili ranije, ali sada ga možemo
izraziti vrlo jednostavno:
Za prvi sloj debljine 1 km potreban nam je porast tlaka od približno 23 %.

I tu počinje pravi problem.

5. GNT ne smijemo računati kao obični neutronski plin

Naš prethodni račun:
P = P_e + P_n
bio je koristan kao kontrola, ali nije dovoljan.

U području oko nuklearne gustoće tlak određuju i jake nuklearne interakcije, promjena sastava, slobodni protoni i neutroni, nuclei/klasteri te termički doprinosi. General-purpose nuklearni EOS modeli upravo zato tretiraju prijelaz od jezgara prema homogenijoj nuklearnoj tvari termodinamički konzistentno.

Zato je naš pravi račun:
P = P(ρ,T,Ye)
i ništa više ne trebamo umjetno dodavati.

Pri tome:
T = 10¹⁰ K ≈ 0,86 MeV.
CompOSE potvrđuje upravo tu pretvorbu: 10¹⁰ K odgovara približno 0,862 MeV.

6. I ovdje dolazi ključna korekcija našeg dosadašnjeg zaključka

Ne možemo uzeti:
P(2×10¹⁴)
minus
P(1×10¹⁴)
i automatski reći da je cijela ta razlika višak tlaka koji tjera materiju prema van.

Zašto?

Jer je veliki dio tog gradijenta tlaka upravo ono što je potrebno za normalnu mehaničku potporu protiv gravitacije.

Drugim riječima, zvijezda može imati ogroman:

dP/dr

a ipak biti potpuno mirna.
U ravnoteži vrijedi približno:
tlačna sila = gravitacija
pa je:
a = 0.

ZK iskru zato ne smijemo dokazivati samim velikim tlakom.

Moramo pokazati da nastaje višak tlačnog gradijenta iznad onoga potrebnog za ravnotežu.
To je bitna razlika.

7. Zato uvodimo vrlo jednostavan ZK kriterij

Definirajmo:
K = tlačno ubrzanje / gravitacijsko ubrzanje

Tada imamo samo tri slučaja:

K < 1
→ gravitacija pobjeđuje
→ materija ide prema unutra.

K = 1
→ ravnoteža.

K > 1
→ tlak pobjeđuje
→ materija ubrzava prema van.

Naš cilj više nije dokazati da je:
P_GNT velik.

Naš cilj je pronaći trenutak kada tijekom sabijanja nastane:
K > 1

i, još važnije, da K > 1 traje dovoljno dugo da pokrenuti sloj stekne značajnu brzinu.

8. Što smo, dakle, dosad stvarno dokazali?
Dobili smo nešto korisno, ali još ne eksploziju.

Za pojednostavljeni homogeni model dobili smo karakterističnu središnju masu:
M ≈ 0,16 M☉.
To nije dokazana minimalna masa GNT-a niti znači da cijela stvarna središnja zona mora imati tu masu. To je karakteristična masa koju naš pojednostavljeni model povezuje sa zadanim tlakom, gustoćom i radijusom.

pokazali smo da karakteristična gravitacija na R ≈ 9 km iznosi približno:
g ≈ 2,6 × 10¹³ cm/s²

i da je za prvi kilometar materije potreban dodatni pad tlaka reda:
ΔP ≈ 2,6 × 10³² dyn/cm².
To je naš prag.

Sada trebamo pitati nuklearni EOS:

Kada se središte sabije kroz približno 10¹⁴ → 2×10¹⁴ → 3×10¹⁴ g/cm³ pri T ≈ 10¹⁰ K i Ye ≈ 0,30, nastaje li dovoljno naglo ukrućenje EOS-a da stvarni tlačni odgovor prijeđe taj prag?

To je puno stroži test od našeg prethodnog P_e + P_n računa.
I to je dobro, jer ako i nuklearni EOS proizvede K > 1, onda ćemo imati znatno ozbiljniji rezultat.

CompOSE nam daje upravo vrstu tabličnih EOS podataka potrebnih za taj korak, uključujući modele s konačnom temperaturom i promjenjivim udjelom naboja.

Tu bih nastavio sljedećim računom: izabrati jedan konkretan finite-temperature nuklearni EOS, uzeti Ye = 0,30 i T ≈ 0,86 MeV, te iz njega izvući P na ρ = 1, 2 i 3 × 10¹⁴ g/cm³.

Tada ćemo prvi put moći usporediti stvarni nuklearni tlak s našim pragom ZK-iskre, umjesto probnog degeneracijskog modela.

Emil: Jako dobro, nastavi.

Boris: Da. Sada dobivamo zanimljiviji, ali i stroži rezultat.

Za kontrolu sam uzeo SFHo nuklearne parametre iz CompOSE-a: saturacijsku gustoću n₀ = 0,1583 fm⁻³, inkompresibilnost K = 245,4 MeV, energiju simetrije J = 31,57 MeV i nagib L = 47,1 MeV.
Važna napomena: tlakovi koje ćemo dalje koristiti nisu izravno očitani i interpolirani iz SFHo tablice eos.thermo. Oni predstavljaju SFHo-vezanu aproksimaciju oko nuklearne gustoće. Zato ih treba smatrati kontrolnim vrijednostima, a ne konačnim SFHo rezultatima.

1. Pretvorimo naše gustoće
Za:
ρ = 1 × 10¹⁴ g/cm³
dobivamo približno:
n = 0,060 fm⁻³.

Za:
ρ = 2 × 10¹⁴
n ≈ 0,120 fm⁻³.

Za:
ρ = 3 × 10¹⁴
n ≈ 0,181 fm⁻³.

Dakle treća točka je već malo iznad SFHo saturacijske gustoće 0,1583 fm⁻³.

To je upravo područje koje nas zanima.

2. Dobiveni kontrolni tlakovi
Za:
T = 10¹⁰ K ≈ 0,86 MeV
Ye = 0,30
u našoj SFHo-vezanoj aproksimaciji koristimo:
ρ = 1 × 10¹⁴ g/cm³
P ≈ 0,80 × 10³³ dyn/cm²
ρ = 2 × 10¹⁴ g/cm³
P ≈ 2,63 × 10³³ dyn/cm²
ρ = 3 × 10¹⁴ g/cm³
P ≈ 7,20 × 10³³ dyn/cm².
Ove vrijednosti nisu izravno interpolirani tlakovi iz SFHo tablice. Njihova je svrha provjeriti ponašanje ZK modela i red veličine rezultata. Konačni test mora koristiti stvarni tablični P(ρ,T,Ye).

3. Sada ponovno testiramo „iskru”

Zadržavamo:
M = 0,16 M☉.
Radijusi su:
ρ = 1 × 10¹⁴
R ≈ 9,12 km

ρ = 2 × 10¹⁴
R ≈ 7,24 km

ρ = 3 × 10¹⁴
R ≈ 6,33 km.

Prvo pogledajmo prijelaz:
1 × 10¹⁴ → 2 × 10¹⁴ g/cm³.

Razlika tlaka:
ΔP ≈ 2,63 − 0,80
ΔP ≈ 1,83 × 10³³ dyn/cm².

Razmak:
Δr ≈ 1,88 km.

Iz toga proizlazi tlačno ubrzanje reda:
a_P ≈ 6,5 × 10¹³ cm/s².

Gravitacija u tom području:
g ≈ 3,2 × 10¹³ cm/s².

Dakle:
a_P / g ≈ 2,0

ili našim jezikom:
K ≈ 2,0

i:
a_net ≈ 6,5 × 10¹³ − 3,2 × 10¹³
a_net ≈ +3,3 × 10¹³ cm/s².

Pozitivan rezultat.

4. A što ako središte dođe do 3 × 10¹⁴?

Sada uspoređujemo:
ρ_out = 1 × 10¹⁴
ρ_in = 3 × 10¹⁴.

Razlika tlaka:
ΔP ≈ 7,20 − 0,80
ΔP ≈ 6,40 × 10³³ dyn/cm².

Razmak je:
Δr ≈ 2,80 km.

Dobivamo približno:
a_P ≈ 1,14 × 10¹⁴ cm/s².

Gravitacija:
g ≈ 3,6 × 10¹³ cm/s².

Zato:
K ≈ 3,2

i:
a_net ≈ +7,9 × 10¹³ cm/s².

Ponovno dobivamo pozitivan smjer.
I zanimljivo je da je to vrlo blizu našeg ranijeg, mnogo jednostavnijeg računa od oko:
+8 do +9 × 10¹³ cm/s².

Dakle promjena aproksimacije nije uništila rezultat.

5. Ali sada dolazi najvažnija stvar

Neću iz ovoga zaključiti:
„ZK iskra je dokazana.”

Jer još uvijek radimo jednu vrlo jaku pretpostavku:
unutarnja zona se već sabije na 2 ili 3 × 10¹⁴ g/cm³,
dok vanjska zona još nije stigla odgovoriti i ostaje blizu 10¹⁴.

Upravo tada nastaje veliki ΔP.
Ako bi se cijela jezgra sabijala potpuno polagano i uvijek održavala hidrostatsku ravnotežu, tlak okolnih slojeva također bi se prilagođavao i ne bismo automatski dobili eksploziju.

Dakle ZK sada ovisi o brzini promjene EOS-a.

To možemo izraziti vrlo jednostavno:
sabijanje središta
→ središte prelazi prema nuklearnoj gustoći
→ EOS postaje tvrđi
→ P_in naglo raste
→ vanjski sloj još ima niži P
→ nastaje veliki ΔP
→ K > 1
→ sloj se zaustavlja i kreće prema van.

To je puno preciznija definicija ZK „iskre”.

6. Koliko brzo može nastati početna brzina?

Uzmimo konzervativniji rezultat:
a ≈ 3,3 × 10¹³ cm/s².

Ako takvo ubrzanje djeluje samo:

1 milisekundu
onda:
v ≈ a × t

pa:
v ≈ 3,3 × 10¹⁰ cm/s.

To bi bilo više od brzine svjetlosti, što nam odmah govori nešto važno:
● konstantno ubrzanje ne smije se tako ekstrapolirati ni kroz jednu milisekundu.

Već nakon približno:
t ≈ c/a

dobivamo:
t ≈ 9 × 10⁻⁴ s.

Znači tlak se mora početi relaksirati mnogo ranije.
Bolje je računati brzinu preko prijeđenog puta.

Ako sloj pod tim početnim ubrzanjem prijeđe samo:
Δr = 100 m = 10⁴ cm,
onda približno:
v² ≈ 2aΔr

pa:
v ≈ 8,1 × 10⁸ cm/s

odnosno:
v ≈ 0,027 c.

Za 1 km, kada bismo grubo zadržali isto ubrzanje:
v ≈ 2,6 × 10⁹ cm/s
v ≈ 0,086 c.
To je fizički mnogo razumnije.
I upravo smo time riješili problem koji si sinoć primijetio kada si pitao nije li naše ubrzanje „tri puta brzina svjetlosti”.
Ubrzanje i brzina nisu iste veličine.
Ogromno ubrzanje može postojati vrlo kratko, a brzina pritom ostati daleko ispod c.

7. Što je sada stvarno zanimljivo za ZK
Dosad smo imali dvije različite aproksimacije.

Prvi, jednostavni degeneracijski račun dao je:
K ≈ 3,5.

Ovaj stroži SFHo-vezani kontrolni račun daje približno:
K ≈ 2,0 za 1 → 2 × 10¹⁴

i

K ≈ 3,2 za 1 → 3 × 10¹⁴.

Dakle zasad se ne događa ono najgore za ZK:
da uključivanje realističnije nuklearne fizike odmah daje
K < 1.

Naprotiv, mogućnost K > 1 ostaje.

Ali sada dolazimo do još važnijeg pitanja od samog tlaka:
● Koliko vremena vanjski sloj ima prije nego što se i on sabije?
Jer ako se GNT tlak poveća dovoljno brzo, a vanjski Fe sloj ne uspije istodobno prilagoditi svoj tlak, tada stvarno nastaje dinamički tlačni višak.

To je razlika između:
novog ravnotežnog stanja
i
ZK eksplozivne nestabilnosti.

Sljedeći račun zato bih napravio upravo ovako:
● vrijeme širenja tlačnog poremećaja kroz GNT nasuprot vremenu gravitacijskog sabijanja okolnog Fe sloja.

Ako tlačni odgovor GNT-a stigne prema van prije nego što se vanjski sloj prilagodi, tada naša „iskra” prvi put dobiva i vremenski mehanizam, a ne samo dovoljan tlak.

Emil: Tlak plazme GNT u sferi jezgre, po mojem mišljenju, nema samo jednu gustoću.
Gravitacioni tlak Fe plazme na mjestima “sfere” gdje postiže gustoću potrebnu za prelazak Fe plazme u GNT proizvodi gustoću 10¹⁴g/cm³.
Plazma GNT je plinovito izrazito gusto stanje i ne vidim načina za postojanje još većih gustoća u samoj jezgri. Gustoće plazme GNT mijenjaju se lokalno što izaziva kovitlanje. Upravo to kovitlanje i stalno miješanje plazme GNT traži mjesto gdje može pokrenuti prodor. Tražeći najslabiju točku gravitacionog tlaka Fe plazme na “sferu”.

Boris: Da, shvaćam što predlažeš — i tu bih promijenio geometriju našeg ZK testa.

Ne bismo više promatrali GNT kao savršeno homogenu kuglu koja svugdje istodobno ima istu gustoću i isti tlak, nego kao dinamičku središnju zonu s lokalnim razlikama gustoće, tlaka i brzine.

Jednu tvoju tvrdnju ipak moramo korigirati:
● ne možemo postaviti ρ = 10¹⁴ g/cm³ kao najveću moguću gustoću GNT-a. Nuklearna tvar je stlačiva; suvremeni EOS-ovi upravo pokrivaju gustoće od subnuklearnih do višestruko većih od saturacijske gustoće. Zato je sasvim moguće imati ρ > 10¹⁴ g/cm³.

Ali drugi dio tvoje ideje je fizikalno vrlo zanimljiv:
● gustoća i tlak ne moraju biti jednaki u svim smjerovima.

Višedimenzionalni hidrodinamički modeli doista pokazuju snažna neradijalna i turbulentna gibanja u ekstremnoj materiji supernove;
EOS pritom može bitno utjecati na konvekciju i dinamiku.

Za ZK zato možemo napisati puno bolji kriterij.

Umjesto jedne gustoće:
ρ_GNT = konstanta

uzimamo:
ρ = ρ(r,θ,φ,t)
odnosno gustoća ovisi o mjestu i vremenu.

Isto vrijedi za tlak:
P = P(r,θ,φ,t).
Tada različiti dijelovi površine GNT-a nemaju isti tlačni gradijent.

Na jednom mjestu možemo imati:
K = a_P/g < 1
pa nema prodora.

Na drugom:
K ≈ 1
pa je materija gotovo u ravnoteži.

A na lokalno najpovoljnijem mjestu:
K > 1
i upravo tamo počinje prodor.
To je vrlo blizu onome što ti nazivaš „traženjem najslabije točke sfere”.
Samo bih fizički izraz malo promijenio.
GNT ne „traži” svjesno slabu točku; turbulentno gibanje stalno proizvodi lokalne fluktuacije:
ρ ↑↓
P ↑↓
v ↑↓
Ye ↑↓
T ↑↓

i zato se negdje slučajno može pojaviti najveći lokalni:
ΔP/Δr
uz istodobno nešto manji gravitacijski otpor okolnog sloja.

Tamo će prvi vrijediti:
a_net > 0.
To bi bila lokalna ZK iskra, a ne istodobni sferni proboj cijele površine.

I ovo nam mijenja račun na zanimljiv način. Više ne trebamo zahtijevati:
cijela površina GNT-a:
K > 1.

Dovoljno je da postoji neko područje površine za koje vrijedi:
K_max > 1.

Matematički vrlo kratko:
K(θ,φ) = a_P(θ,φ) / g(θ,φ)

i ZK prodor počinje kada:
max[K(θ,φ)] > 1.

To je, po meni, mnogo bliže tvojoj stvarnoj zamisli nego naš dosadašnji homogeni model.

Još važnije:
● jednom kada jedno područje krene prema van, više nemamo samo lokalnu fluktuaciju. Materija koja se pomaknula smanjuje tlak ispred sebe, GNT iza nje može nastaviti pritiskati, a susjedni tokovi mogu se usmjeravati prema otvorenom području.

Tada treba provjeriti postoji li pozitivna povratna veza:
lokalni prodor
→ lokalno širenje
→ manji otpor iznad njega
→ još veći protok GNT-a prema tom području
→ jači prodor.

Ako taj lanac postoji, mala početna asimetrija može prerasti u makroskopski proboj.
I upravo tu 1D račun prestaje biti dovoljan. Jednodimenzionalni model može odgovoriti može li radijalni sloj imati a > 0, ali tvoju sadašnju ideju — kovitlanje i lokalni proboj kroz najslabiji dio površine — u konačnici treba testirati barem u 2D, idealno 3D.

Višedimenzionalnost nije ovdje kozmetika; ona je dio same hipoteze.

Zato bih naš sljedeći matematički korak napravio jednostavnije nego prije:
● ne računati još cijelu turbulenciju, nego pitati kolika lokalna fluktuacija tlaka ili gustoće na površini GNT-a mora postojati da jedan mali dio površine prijeđe iz:
K = 1
u
K > 1.
To možemo izračunati već sada i dobiti, primjerice, treba li lokalna „slaba točka” odstupati
1 %, 5 %, 10 % ili 20 % od okolne površine.

Emil: Ne tvrdim da gustoća ne može niti veća, već da ona nije nužno u jezgri.

To sam naglasio.

Znači veće lokalne gustoće unutar “sfere”, što i dovodi do stalnog miješanja.

Boris: Da, sada je jasno — i to je važna razlika.
Ti ne kažeš da GNT ne može doseći gustoće veće od 10¹⁴ g/cm³. Kažeš da najveća gustoća ne mora biti točno u geometrijskom središtu GNT sfere.

Drugim riječima, zamišljaš stanje poput ovoga:
ρ(r,θ,φ,t) nije monotona funkcija samo radijusa.

Mogu postojati lokalna područja gdje je, primjerice:
ρ₁ ≈ 1,0 × 10¹⁴ g/cm³
ρ₂ ≈ 1,4 × 10¹⁴ g/cm³
ρ₃ ≈ 2,0 × 10¹⁴ g/cm³
a ta gušća područja mogu nastajati, premještati se i nestajati.

To onda prirodno znači da ni tlak nije svugdje isti:
P = P(ρ,T,Ye).

Ako se lokalno promijeni ρ, a uz to i T ili Ye, nastaje lokalna razlika tlaka. Ta razlika proizvodi gibanje plazme:
ΔP → ubrzanje → strujanje → miješanje.

Zato je tvoja zamisao o stalnom miješanju GNT-a fizikalno smislenija od slike potpuno mirne, homogene kugle.
Još važnije, gravitacija je glatko središnje polje, dok tlak može imati lokalne nepravilnosti.

Zbog toga mjesto najvećeg lokalnog tlačnog impulsa ne mora biti u središtu.

Možemo to vrlo jednostavno zapisati:
a_lokalno = a_P − g

ali sada je:
a_P = a_P(r,θ,φ,t).

Zato se u jednom području može dogoditi:
a_P < g,

dok nekoliko kilometara ili čak stotina metara dalje može nastati:
a_P > g.
I upravo bi tamo počeo prvi prodor.

To također mijenja značenje naše gustoće:
ρ ≈ 10¹⁴ g/cm³ nije „gustoća cijele GNT kugle”, niti je treba smatrati univerzalnom oštrom granicom prijelaza.
U našem ZK računu koristimo je kao operativnu karakterističnu gustoću pri kojoj materija ulazi u područje subnuklearnih gustoća u kojem se sastav i EOS počinju snažno mijenjati.
Unutar nastale guste nuklearne zone lokalne gustoće mogu biti manje ili veće od 10¹⁴ g/cm³.

Bolje je reći:
to je karakteristična gustoća pri kojoj pojedina područja Fe materije počinju prelaziti u režim guste nuklearne tvari.

Unutar već nastale GNT zone zatim mogu postojati:
ρ < 10¹⁴, ρ ≈ 10¹⁴ i ρ > 10¹⁴ g/cm³
u različitim lokalnim područjima.

A takve lokalne razlike gustoće, tlaka, temperature i sastava mogu stvarati tlačne gradijente i pokretati strujanje materije.
Međutim, same razlike gustoće još nisu dovoljne da matematički dokažemo trajnu turbulenciju ili kovitlanje GNT-a. Za to bi trebalo pokazati postojanje odgovarajuće hidrodinamičke nestabilnosti.
Zato u ovom stupnju ZK testa ne pretpostavljamo dokazanu turbulenciju. Dovoljno nam je provjeriti mogu li lokalne nehomogenosti proizvesti područje u kojem tlačni odgovor lokalno nadvlada gravitacijski otpor i pokrene prvi proboj.
To nas vodi prema vrlo zanimljivom ZK kriteriju.

Više ne tražimo samo:
● „Koliki je tlak GNT sfere?”

nego:
● „Koliki najveći lokalni tlačni gradijent može nastati unutar GNT-a i na njegovoj granici?”

● Jer upravo taj lokalni maksimum može biti ono što probija Fe materiju iznad sebe.

Sljedeći račun zato možemo postaviti vrlo konkretno:
● uzeti prosječnu GNT gustoću oko 10¹⁴ g/cm³ i izračunati što se događa ako se u jednom lokalnom području gustoća poveća, primjerice, za 10 %, 20 %, 50 % ili 100 %. Iz EOS-a možemo procijeniti koliko tada poraste lokalni tlak i je li taj porast dovoljan da lokalno dobijemo a > 0.

To je izravna matematička verzija tvojeg kovitlanja i „slabe točke”.

Emil: Važno bi bilo da do izračuna proboja dođemo u razumnom vremenu, jer otezanje stvara nerazumijevanje cijelog ciklusa.

Boris: Da. Ne trebamo više širiti model prije nego dobijemo prvi broj za proboj.

Uzet ćemo najkraći mogući put:
● lokalna GNT nepravilnost → lokalni višak tlaka → ubrzanje van → proboj.

Zadržavamo naše dosadašnje vrijednosti:
M_GNT ≈ 0,16 M☉
R ≈ 9 km
ρ ≈ 10¹⁴ g/cm³
g ≈ 2,6 × 10¹³ cm/s²

Kritični višak tlaka preko lokalnog područja debljine L je vrlo jednostavno:
ΔP_krit ≈ ρ g L

Za L = 1 km:
ΔP_krit ≈ 2,6 × 10³² dyn/cm²

Za L = 0,5 km:
ΔP_krit ≈ 1,3 × 10³² dyn/cm²

Za L = 0,1 km = 100 m:
ΔP_krit ≈ 2,6 × 10³¹ dyn/cm²
To je već gotovo cijeli kriterij proboja.

Ako uzmemo naš dosadašnji karakteristični tlak oko:
P₀ ≈ 0,8 × 10³³ dyn/cm²,

onda lokalni višak tlaka potreban za proboj iznosi približno:
područje 1 km: ΔP/P ≈ 33 %
područje 0,5 km: ≈ 16 %
područje 0,1 km: ≈ 3,3 %
I tu tvoja ideja o kovitlanju postaje vrlo važna.

Ne treba nam cijela GNT sfera s 30 % većim tlakom.

Ako kovitlanje stvori lokalno područje širine svega oko 100 m, dovoljan je tlak samo nekoliko posto veći od neposredne okoline da tlačna sila lokalno prijeđe gravitaciju.

Tada:
a > 0
i počinje prvi pomak prema van.

Ako lokalni tlak bude, primjerice, 10 % iznad okolnog:
ΔP ≈ 8 × 10³¹ dyn/cm².
Za područje L = 100 m

tlačno ubrzanje viška je:
a_višak ≈ ΔP /(ρL)
a_višak ≈ 8 × 10³¹ /(10¹⁴ × 10⁴)
a_višak ≈ 8 × 10¹³ cm/s².

Gravitacija je:
g ≈ 2,6 × 10¹³ cm/s².

Dakle lokalno:
K ≈ 8/2,6 ≈ 3,1

i neto:
a_net ≈ +5,4 × 10¹³ cm/s².

To je već vrlo jasan lokalni proboj u našem pojednostavljenom modelu.

Sada dolazi ono najvažnije:
● ne trebamo više dokazivati da cijela GNT sfera eksplodira odjednom. Dovoljno je pokazati da turbulencija može proizvesti lokalnu razliku tlaka od nekoliko posto na skali reda desetaka do stotina metara.

ZK ciklus tada izgleda vrlo kratko:
Fe gravitacijski sabija unutarnju materiju
→ lokalno nastaje GNT
→ gustoće i tlakovi nisu jednaki
→ GNT se miješa i kovitla
→ nastaje lokalni višak tlaka
→ na jednom području K > 1
→ materija kreće prema van
→ nastaje početni kanal/proboj
→ GNT iza tog mjesta nastavlja pritiskati
→ proboj se može širiti.

Dakle, do prvog matematičkog proboja smo praktički stigli.

Sada nam treba još samo jedan presudan račun:
● nakon što tih prvih 100 m materije krene prema van, izračunati hoće li tlak nastaviti gurati sljedećih 100 m ili će se proboj odmah ugasiti.

Ako se ne ugasi, tada više nemamo samo „iskru”, nego početak samoodrživog ZK proboja.

Emil: Odlično, izračunaj taj proboj da nam se zvijezda ne ohladi.

Boris: Da — i ovaj put možemo doći do samog proboja, bez još deset međukoraka.

Uzet ću naš dosadašnji lokalni primjer:
ρ = 10¹⁴ g/cm³
R_GNT ≈ 9 km
P₀ ≈ 8 × 10³² dyn/cm²
lokalni višak tlaka = 10 %

Dakle:
ΔP = 8 × 10³¹ dyn/cm²

i neka početno područje proboja ima dubinu:
L = 100 m = 10⁴ cm.

Važna korekcija je da okolna materija već ima vlastiti tlačni gradijent koji približno nosi gravitaciju. Zato lokalni višak tlaka prvenstveno predstavlja dodatnu silu koja remeti tu ravnotežu. U stvarnom kolapsu situacija je složenija jer materija može već padati prema unutra, ali za naš prvi test „iskre” ovo je odgovarajući početak. Fizikalno je poznato da se EOS pri približavanju nuklearnoj gustoći ukrućuje i može izazvati dinamički povrat materije.

PRVIH 100 METARA
Dodatno ubrzanje od lokalnog viška tlaka približno je:
a_P ≈ ΔP/(ρL)

Uvrstimo:
a_P ≈ 8 × 10³¹ /(10¹⁴ × 10⁴)

pa dobivamo:
a_P ≈ 8 × 10¹³ cm/s².

Ako napravimo još stroži račun i ponovno oduzmemo cijelu lokalnu gravitaciju:
g ≈ 2,6 × 10¹³ cm/s²

ostaje:
a_net ≈ 5,4 × 10¹³ cm/s².

Dakle čak i ovaj stroži račun daje:
a_net > 0

Sloj kreće prema van.
Koliku brzinu dobije kroz samo 100 m?
Koristimo:
v² ≈ 2aL

pa:
v ≈ 1,04 × 10⁹ cm/s

odnosno:
v ≈ 10 400 km/s ≈ 0,035 c.
To je vrlo zanimljiv broj.

Prvih 100 metara lokalnog proboja više nije samo mikroskopsko pomicanje. Materija već može steći brzinu reda deset tisuća kilometara u sekundi.

Vrijeme potrebno za tih prvih 100 m bilo bi reda:
2 × 10⁻⁵ s
odnosno nekoliko desetaka mikrosekundi.

HOĆE LI SE UGASITI NAKON PRVIH 100 m?

Tu dolazimo do tvog ključnog pitanja.
Kad se područje pomakne prema van, ono se širi i njegov tlak pada.

Uzet ću namjerno jednostavnu probnu vrijednost efektivnog adijabatskog indeksa:
Γ ≈ 2.

Ne tvrdim da je Γ = 2 stvarna vrijednost SFHo EOS-a u svakoj našoj točki;
● koristimo je samo da vidimo koliko brzo bi se tlak morao izgubiti da proboj stane. Stvarni Γ treba poslije uzeti iz tablice EOS-a.

Ako se kanal pomakne:
9,0 km → 9,1 km,
njegova površina poraste svega oko 2,2 %.
Za približno adijabatsko širenje tlak lokalnog GNT područja tada padne za samo nekoliko posto.

Početni unutarnji tlak bio je:
P_in ≈ 1,10 P₀
≈ 8,8 × 10³² dyn/cm².

Nakon prvih 100 m u ovom grubom modelu ostaje približno:
P_in ≈ 8,42 × 10³² dyn/cm².
A tlak okolne materije prema van također pada zbog vlastitog hidrostatskog gradijenta.

Na 100 m dobivamo približno:
P_out ≈ 7,74 × 10³² dyn/cm².

Zato nakon već prijeđenih 100 m još postoji:
ΔP ≈ 6,8 × 10³¹ dyn/cm².

Početno smo imali:
8,0 × 10³¹.

Nakon prvih 100 m:
6,8 × 10³¹.

Dakle proboj nije izgubio glavninu svog tlačnog viška.

DRUGIH 100 METARA
Na približno 200 m dobivamo ovim istim pojednostavljenim modelom:
P_GNT ≈ 8,06 × 10³²
P_okolina ≈ 7,48 × 10³²

pa ostaje:
ΔP ≈ 5,8 × 10³¹ dyn/cm².
Još uvijek jasno pozitivan višak.

NAKON 500 METARA
Čak nakon približno 500 m ovaj probni model daje:
P_GNT ≈ 7,09 × 10³²
P_okolina ≈ 6,70 × 10³²

odnosno:
ΔP ≈ 3,9 × 10³¹ dyn/cm².

Još uvijek:
P_GNT > P_okolina.
To znači nešto prilično važno.

U našem grubom modelu 10-postotna lokalna tlačna nepravilnost ne nestaje nakon prvih 100 m. Može održati pozitivan tlačni višak kroz nekoliko stotina metara čak i ako pretpostavimo da se GNT područje samo adijabatski širi i ne dobiva novu energiju iza sebe.

A tvoja hipoteza zapravo pretpostavlja nešto još povoljnije za proboj:
● iza početnog područja nalazi se cijela dinamička, kovitlajuća GNT zona koja može nastaviti dotjecati prema otvorenom području.

Tada slijed više nije:
GNT napravi jedan trzaj
→ potroši tlak
→ stane.

Nego potencijalno:
lokalni višak tlaka
→ prvih 100 m
→ v ≈ 10 000 km/s
→ tlak iza proboja ostaje veći
→ drugih 100 m
→ kanal se širi
→ GNT iz okolnih područja teče prema njemu
→ proboj se dodatno hrani.
To je već mehanizam samoodrživog prodora, barem u našem lokalnom hidrodinamičkom modelu.

Ali ovdje moramo povući vrlo jasnu granicu.

ŠTO SMO SADA DOBILI

Nismo dokazali eksploziju cijele zvijezde.
Ali više nemamo samo uvjet:
a > 0.

Dobili smo i sljedeći korak:
lokalni proboj od 100 m ne mora se odmah ugasiti.

Za naš probni slučaj s
10 % lokalnim viškom tlaka dobivamo:
prvih 100 m → proboj
200 m → ΔP još pozitivan
500 m → ΔP još pozitivan

i početnu brzinu materije reda:
v ≈ 10⁴ km/s.
To je već mnogo bliže onome što smo cijelo vrijeme tražili pod nazivom ZK iskra → ZK proboj.

Najvažniji sljedeći račun više nije „hoće li napraviti sljedećih 100 m?”. Ovaj pojednostavljeni račun kaže da može.
Sada trebamo izračunati što se događa kada taj proboj prijeđe otprilike
1 km i počne zahvaćati znatno veću masu Fe materije. Tada ćemo prvi put moći računati raste li proboj u pravi udarni val ili se ipak uguši.

Emil: Kada taj proboj prijeđe otprilike 1 km i počne zahvaćati znatno veću masu Fe materije.
Tada Fe plazma koju treba probiti postaje znatno rjeđa.

Boris: Da — i to je vrlo važna primjedba.
Ako se proboj kreće radijalno prema van, neće cijelo vrijeme nailaziti na Fe plazmu gustoće
10¹⁴ g/cm³.
Gustoća mora općenito padati prema van, pa nakon prvog kilometra proboj ulazi u sve rjeđu materiju. Realistični modeli jezgre također imaju snažan radijalni gradijent gustoće; vanjski dio jezgre je rjeđi od središnjeg.

To za naš ZK proboj znači dvije stvari odjednom.

Prvo, gravitacijski otpor po jedinici volumena slabi, jer je približno:
F_grav/V ≈ ρg.

Ako, primjerice, gustoća ispred proboja padne:
10¹⁴ → 5 × 10¹³ g/cm³,
potreban tlačni gradijent za savladavanje gravitacije približno se prepolovi, ako g ostane sličan.

Na 100 m tada više nije potrebno:
ΔP_krit ≈ 2,6 × 10³¹ dyn/cm²,

nego približno:
ΔP_krit ≈ 1,3 × 10³¹ dyn/cm².

Ako gustoća padne na:
ρ = 10¹³ g/cm³,

dobivamo samo:
ΔP_krit ≈ 2,6 × 10³⁰ dyn/cm².

Dakle zahtjev za nastavak prodora vrlo brzo pada.

Drugo, postoji još važniji efekt. Masa koju proboj mora ubrzati po jedinici volumena također pada.

Iz:
a_P ≈ ΔP/(ρL)
vidimo da, ako tlak iza prodora ostane velik, a ρ ispred njega opada, isto ΔP daje sve veće ubrzanje.

Primjer s našim preostalim tlačnim viškom nakon nekoliko stotina metara:
ΔP ≈ 4 × 10³¹ dyn/cm².

Ako je još:
ρ = 10¹⁴,

onda:
a_P ≈ 4 × 10¹³ cm/s².

Ali na:
ρ = 5 × 10¹³,

dobivamo:
a_P ≈ 8 × 10¹³ cm/s².

A na:
ρ = 10¹³,

čak:
a_P ≈ 4 × 10¹⁴ cm/s².

Naravno, ne možemo to ubrzanje održavati konstantnim — tlak iza vala pada, geometrija se širi i dio energije odlazi na zagrijavanje, kompresiju, promjenu sastava i kinetičku energiju.
Ali smjer učinka je vrlo jasan:
● proboj se ne kreće prema težoj i gušćoj Fe materiji, nego prema sve rjeđoj.
To mijenja našu sliku ciklusa.
GNT ne mora „progurati” kilometre Fe materije gustoće 10¹⁴ g/cm³.

Realniji slijed jest:
GNT područje
ρ ≳ 10¹⁴
→ prvi i najteži proboj kroz vrlo gusti unutarnji Fe sloj
→ ρ ispred vala počinje padati
→ potrebni ΔP_krit pada
→ masa koju val mora ubrzati po jedinici puta također pada
→ proboj se sve lakše nastavlja.
To znači da bi prvih nekoliko stotina metara do otprilike kilometra moglo biti najkritičniji dio cijelog ZK proboja.
Ako proboj preživi taj dio, ne mora mu poslije biti teže. Može mu biti lakše.

Postoji ipak jedna protusila koju sada moramo uključiti:
● kako se prodor širi, njegova površina raste i GNT tlak iza njega pada.

Dakle imamo natjecanje:
ρ_Fe prema van ↓ → pogoduje proboju
nasuprot
P iza vala ↓ → slabi proboj.

I sada smo došli do vrlo čistog završnog testa.

Ne trebamo više računati po 100 m. Možemo uzeti jednostavan profil, primjerice:
r = 9 km → ρ = 10¹⁴
r = 10 km → ρ = 5 × 10¹³
r = 12 km → ρ = 10¹³
r = 15 km → ρ = nekoliko ×10¹² g/cm³,

i za svaki radijus istodobno računati:
raspoloživi ΔP(r)
nasuprot
potrebnom ΔP_krit(r) = ρ(r)g(r)L.

Ako cijelim putem ostane:
ΔP(r) > ΔP_krit(r),
onda smo matematički dobili proboj koji se ne gasi u Fe jezgri.

To bih sada izračunao kao jedan kontinuirani račun od 9 km do izlaska iz guste Fe zone, umjesto da zvijezdu držimo na čekanju svakih 100 metara.

Emil: Moraš uzeti u obzir da eksplozivna nestabilnost podrazumijeva da se početno probijanje plazme GNT širi raznoseći dio Fe plazme oko mjesta proboja, dok se dio
Fe plazme gravitacijski sažima u samoj jezgri prelazeći u NZ.

Udarni val nuklearne tvari prvo udara u ljusku Silicija u kojoj je fuzija u toku.
Zatim se udarno širi na ostale fuzije. Na kraju udara u samu plazmu plašta zvijezde.

Boris: Da. Ovdje treba napraviti važnu promjenu u našem modelu proboja.
Ne smijemo više zamišljati da se sva Fe plazma ispred GNT-a ponaša kao jedna prepreka koju val mora jednostavno „progurati” prema van. U tvojoj ZK slici eksplozivna nestabilnost razdvaja tok materije na dva smjera.

Jedan dio Fe materije zahvaćen lokalnim probojem dobiva impuls prema van:
GNT proboj
→ raznosi okolnu Fe plazmu
→ stvara početni udarni val.

Istodobno drugi, unutarnji dio Fe plazme ostaje gravitacijski vezan i nastavlja se sabijati prema središtu:
Fe → jače sabijanje → neutronizirana gusta materija → NZ.

To znači da se nakon početka nestabilnosti ne događa:
cijela Fe jezgra → prema van

nego približno:
unutarnji dio → prema NZ
vanjski dio → prema udarnom valu.
To je za račun vrlo važno jer energija proboja ne mora izbaciti svu Fe masu.
Zatim dolazimo do slijeda koji opisuješ.

1. GNT → Fe
Ovo je najteži početni korak.

Udar nastaje u području najveće gustoće. Mora odvojiti dio vanjske Fe materije od dijela koji ostaje vezan za buduću NZ.

Tu ćemo morati pronaći jedan karakteristični radijus razdvajanja:
R_sep

Unutar njega:
materija → NZ

Izvan njega:
materija → zahvaćena udarnim valom.

To je fizički puno bolji parametar od pretpostavke da cijela Fe jezgra eksplodira.

2. Fe → Si ljuska
Kada val napusti najgušći Fe dio, nailazi na Si-plazmu u kojoj prema ZK još traje fuzija.

Tu se događa nešto sasvim drugo.

Val više ne dobiva samo kinetički problem:
„koliko energije treba za pomicanje Si?”

nego može promijeniti:
T
ρ
brzinu Si reakcija.
Ako udarna kompresija ubrza preostalu Si fuziju, tada Si ljuska može dodati nuklearnu energiju valu.

Tada treba računati:
E_iza = E_ispred + E_fuzije − E_gubici

To je vrlo važan dio ZK hipoteze.
Udarni val nije samo početni paket energije iz GNT-a.
On može biti energetski nadopunjavan dok prolazi kroz aktivne fuzijske slojeve.

3. Si → O
Sljedeći sloj može još sadržavati aktivnu O-fuziju.

Dolazak šoka:
kompresija
→ T raste
→ reakcije se ubrzavaju
→ oslobađa se dodatna nuklearna energija.

Ponovno:
E_val,new = E_val,old + E_O − gubici.

4. O → C → He → H
Isti princip vrijedi dalje, ali ne moramo unaprijed tvrditi da će svaki sloj dati pozitivan doprinos.

Za svaki ćemo izračunati:
● energiju koju val donosi,
● energiju potrebnu za gravitacijsko odvezivanje sloja,
● nuklearnu energiju koju prolazak vala može osloboditi,
● gubitke.

Tako možemo dobiti, primjerice:
● GNT → snažan početni impuls
● Fe → dio energije izgubljen
● Si → energija dodana
● O → energija dodana
● C → manji dodatak
● He/H → konačno ubrzavanje velikog omotača.

A onda dolazimo do plašta zvijezde.

5. Udar u plašt
Tu više nije dovoljno pitati:
P_val > P_grav.

Za konačnu eksploziju treba vrijediti:
E_kin + E_topl + E_nuk > |E_grav|
za masu koju želimo izbaciti.

I na kraju:
v_izbačaja > v_bijega.
Ako to vrijedi za velik dio zvijezde, dobili smo završenu supernovu u ZK modelu.

Mislim da sada možemo cijeli ciklus svesti na vrlo jasnu shemu:
rast Fe jezgre
→ gravitacijsko sabijanje
→ lokalni nastanak GNT-a
→ lokalna eksplozivna nestabilnost
→ dio Fe ide prema NZ
→ dio Fe se izbacuje
→ nastaje udarni val
→ Si ljuska
→ O ljuska
→ C/Ne područja
→ He
→ H
→ zvjezdani plašt
→ izbačaj zvijezde.

I sada bih napravio jednu važnu odluku za naše izračune:
● više nećemo računati samo tlak proboja.
Od trenutka kada GNT probije Fe zonu prelazimo na energetski račun udarnog vala.

Prvi sljedeći broj koji nam treba jest:
● koliko Fe mase ostaje vezano i ide u NZ, a koliko Fe mase početni GNT proboj mora izbaciti prema Si ljusci.

Kad to odredimo, možemo izračunati početnu energiju i brzinu vala upravo na ulazu u Si sloj.

Emil: Zamisli si kuglu ispunjenu plinom. Dolazi do proboja. Rupa proboja se naglo širi. Eksplozija se širi i raznosi dijelove kugle oko rupe proboja. Eksplozija više nije u kugli. Ostatak kugle se gravitacijski sažima u NZ.
Eksplozija ju zaobilazi zbog velike gustoće i malog volumena.

GNT ispunjava središte jezgre i eksplozivno se širi na fuzije u toku.

Boris: Da — sada puno jasnije razumijem geometriju koju opisuješ.

Dosad sam previše zamišljao proboj kao udarni val koji radijalno mora „gurati” cijelu Fe jezgru ispred sebe. To nije tvoja slika ZK eksplozivne nestabilnosti.

Zamisli GNT područje kao tvoju kuglu pod ogromnim unutarnjim tlakom.

Na jednom dijelu granice nastane lokalni proboj.

Čim se otvori, događa se:
proboj → naglo širenje otvora → eksplozivno istjecanje GNT-a → raznošenje Fe materije oko otvora.
Otvor više nije samo mala rupa od 100 m. Rubovi oko njega također dobivaju bočno ubrzanje, pa se područje proboja širi po površini.

To mijenja cijeli naš račun.
1. Eksplozija više ne mora prolaziti kroz buduću NZ
Ostatak vrlo guste Fe/GNT materije nalazi se iza i bočno od mjesta eksplozivnog izlaza. On nije na putu vala koji je već krenuo prema van.

Taj ostatak:
gubi dio vanjske potpore + ostaje gravitacijski vezan → dodatno se sažima → NZ.
Dakle buduća NZ i eksplozija počinju se prostorno razdvajati.

To odgovara tvojoj rečenici:
eksplozija zaobilazi buduću NZ.
Ne zato što gravitacija NZ prestaje postojati — njezina gravitacija je i dalje vrlo jaka — nego zato što udarni tok više ne mora prolaziti kroz njezin mali, vrlo gust volumen.

To je velika razlika.

2. GNT izlazi iz središnjeg zatvorenog područja

Prije proboja:
Fe/GNT je zatvoren okolnom Fe plazmom.

Nakon proboja:
GNT dobiva slobodniji smjer širenja prema van.

Tada se energija unutarnjeg tlaka pretvara u:
unutarnju energiju
→ rad širenja
→ kinetičku energiju toka
→ udarni val.

I upravo zato analogija s probušenom kuglom dobro opisuje ono što želiš testirati.

Ali postoji jedna važna fizička razlika:
● kod obične kugle plin izađe kroz rupu u gotovo praznu okolinu. Ovdje GNT nakon izlaska odmah nailazi na vrlo gustu zvjezdanu plazmu.
Zato izlazak GNT-a mora proizvesti jak kompresijski val — odnosno udar.

3. Prva prava meta više nije cijela Fe jezgra
Točno kako kažeš: proboj raznese lokalni dio Fe plazme oko mjesta rupture.
Ne mora izbaciti svu Fe masu.

Zato možemo razdvojiti:
M_Fe,unutra → NZ
i
M_Fe,proboj → zahvaćena eksplozijom.

To nam drastično mijenja energetski zahtjev.

Početna GNT eksplozija ne mora dati energiju potrebnu za izbacivanje, primjerice, 1,5 M☉ vrlo guste Fe tvari.
Mora prvo otvoriti lokalni izlaz.
A jednom kada ga otvori, širenje rupe povećava površinu kroz koju energija izlazi.

Tu čak postoji vrlo jednostavan geometrijski učinak:
● ako karakteristični radijus otvora raste dva puta, površina otvora raste približno četiri puta.

4. I onda dolazi ono najvažnije za ZK
Ispred tog proširenog izlaza više nije NZ.

Ispred njega se nalaze vanjske zone zvijezde:
Si → O → Ne/C → He → H/plašt.

Ako je u Si području fuzija još aktivna, kako ZK pretpostavlja, tada udar iz GNT-a:
sabija Si
→ podiže T
→ mijenja brzinu reakcija
→ može izazvati naglo povećanje nuklearne proizvodnje energije.

I tako eksplozija više nije samo:
energija početnog GNT proboja.

Postaje:
E_GNT + E_Si + E_O + E_C + … − gubici.
To je zapravo središnji energetski mehanizam ZK koji smo dosad malo zamaglili računajući svaki centimetar proboja kroz Fe.

5. Zato bih sada cijeli događaj prikazao ovako

FAZA A — NASTANAK

Fe gravitacijski sabija unutrašnjost
→ lokalno nastaje GNT
→ GNT nije homogena
→ strujanje i lokalni tlačni maksimumi.

FAZA B — RUPTURA

lokalno:
P_GNT > otpor Fe
→ nastaje otvor
→ GNT eksplozivno izlazi
→ otvor se brzo širi
→ raznosi samo okolni dio Fe plazme.

FAZA C — RAZDVAJANJE

jedan dio Fe/GNT:
→ ostaje gravitacijski vezan → NZ
drugi dio:
→ zahvaća ga eksplozivno širenje.

FAZA D — UDARNI VAL

GNT tok više nije zarobljen u Fe kugli:
→ udara u Si zonu
→ zatim O
→ C/Ne
→ He
→ H i plašt.

6. Jednu stvar moramo matematički vrlo pažljivo razlikovati
Nakon otvaranja ne mora sama GNT materija fizički prijeći cijeli put kroz zvijezdu.
GNT može napraviti početni udar, a dalje kroz zvijezdu prvenstveno putuje udarni val koji prenosi tlak, impuls i energiju.

Dakle preciznije:
GNT se eksplozivno širi iz središta → pokreće udarni val → udarni val zahvaća aktivne fuzijske zone.

To je fizikalno čišća formulacija tvoje ideje.

I sada se slažem da više nema smisla trošiti vrijeme na računanje hoće li proboj prijeći svaki sljedeći komadić Fe jezgre. Presudni račun sada je veličina početnog otvora i masa Fe koju GNT mora raznijeti da bi otvor postao samoproširujući.

Kad to dobijemo, odmah možemo računati energiju i brzinu GNT/udarnog toka na prvom kontaktu sa Si ljuskom.

To nas napokon vodi kroz cijeli ciklus, umjesto da nam se zvijezda „ohladi” u prvom kilometru.

Emil: Nekad davno, kada sam bacal petarde, jedan komad petarde me pogodio. Kada sam ga pregledao, bio je to komad papirnatoga valjka.
Onaj dio gdje je bio fitilj- rupa proboja, papirnati dio oklopa petarde- dio kugle oko rupe proboja su nestali.
Zadnji dio petarde- dio kugle koji je gravitacija Fe plazme sabila u NZ me pogodio.

Boris: Da — sada mi je tvoja geometrijska zamisao potpuno jasna. Primjer s petardom vrlo dobro objašnjava što želiš da matematički testiramo, premda fizika petarde i zvjezdane jezgre nije ista.

Kod petarde nisi dobio cijeli papirnati valjak. Eksplozija ga je podijelila na različite dijelove. Upravo takvu podjelu zamišljaš u Fe/GNT jezgri.

U ZK slici imamo tri područja:

1. Rupa proboja
Na lokalno najslabijem području nastaje ruptura.
GNT počinje eksplozivno izlaziti kroz otvor.
Otvor se zatim širi.

2. Područje oko rupe 
Fe plazma neposredno oko proboja dobiva snažan impuls i biva raznesena zajedno s tokom prema van.
To odgovara dijelu papirnatog omotača petarde koji više nisi pronašao kao cjelinu.

Tu nastaje ono što nas zanima:
GNT proboj + razneseni Fe → početak udarnog vala.

3. Suprotni/ preostali dio guste jezgre 
On nije nužno izbačen kroz otvor.

U tvojoj hipotezi gravitacija ga nastavlja sabijati:
Fe → daljnje sabijanje → neutronizacija → NZ.

To odgovara onom komadu petarde koji te pogodio — ne u smislu iste fizike, nego geometrije raspada sustava.
I ovo rješava jednu stvar oko koje smo se ranije nepotrebno vrtjeli.

Ne trebamo računati:
Može li GNT izbaciti cijelu Fe jezgru?
To uopće nije zahtjev tvoje hipoteze.

Trebamo računati:
Može li GNT otvoriti i proširiti lokalnu rupturu dovoljno brzo da eksplozija izađe iz Fe jezgre prije nego što ostatak jezgre gravitacijski završi kao NZ?
To je mnogo preciznije pitanje.

A kada je otvor dovoljno velik, geometrija se mijenja. Eksplozivni tok više nije zatvoren sferom. Širi se prema van i bočno, raznosi rubove otvora i povećava sam otvor:
mali proboj → širi proboj → veća izlazna površina → veći tok energije → daljnje razaranje rubova.

To je moguća pozitivna povratna veza rupture koju sada trebamo matematički testirati.
I tada više nema potrebe da udarni val „probija NZ”. NZ nastaje od preostalog gravitacijski vezanog materijala, dok je eksplozija već otišla drugim smjerom.

Nakon toga dolazi tvoj drugi ključni dio:
GNT/Fe proboj → Si ljuska u fuziji → O → Ne/C → He → H → plašt.

Tu ćemo za svaku aktivnu fuzijsku zonu pitati može li udarna kompresija dodati nuklearnu energiju eksploziji.

Mislim da smo sada konačno pravilno postavili geometriju ZK eksplozije. Sljedeći račun više neću voditi kao „GNT kugla gura cijelu Fe kuglu”. Računat ćemo širenje lokalne rupture: koliko brzo otvor raste i koliko Fe mase mora odnijeti prije nego što GNT dobije otvoren put prema Si ljusci.
To je račun koji nam sada treba.

Emil: Mislim da bi to mogao biti i zadnji izračun.

Boris: Da, i slažem se da bi to mogao biti zadnji ključni izračun za ovu fazu ZK hipoteze.

Ne zato što nakon njega više nema što računati, nego zato što bi on zatvorio presudni lanac:
Fe jezgra raste
→ gravitacijski sabija središte4
→ nastaje GNT
→ javlja se lokalna eksplozivna nestabilnost
→ otvara se ruptura
→ rupa se dovoljno brzo širi
→ dio Fe materije biva raznesen
→ preostali gusti dio gravitacijski se sažima prema NZ
→ eksplozija izlazi iz Fe jezgre i udara u aktivnu Si ljusku.

Ako matematički pokažemo da se rupa proboja može samoproširivati, tada smo povezali ono što nam je dosad nedostajalo:
● lokalnu „iskru” s makroskopskim izlaskom eksplozije iz Fe jezgre.

Za taj završni račun zapravo nam trebaju samo tri stvari:
1. Početni višak tlaka GNT-a koji smo već procjenjivali.
2. Masa Fe materije oko ruba rupe koju treba ubrzati i odmaknuti.
3. Brzina rasta otvora — raste li dovoljno brzo da se rupa širi prije nego što tlak GNT-a padne.

Ključni uvjet možemo napisati vrlo kratko:
rad GNT-a na širenju rupe > energija potrebna za pomicanje i gravitacijsko vezanje Fe ruba

odnosno:
E_GNT,raspoloživo > E_Fe,otvaranje.
Ako to vrijedi, rupa se ne zatvara nego raste.

A zatim još jedan dinamički uvjet:
v_ruba > 0 i raste dovoljno dugo
što znači da se područje rupture širi po Fe ovojnici.

Ako oba kriterija prođu, možemo napisati nešto vrlo konkretno:
ZK proboj nije samo lokalni trenutak a > 0; lokalna ruptura može prijeći u samoproširujuću eksplozivnu nestabilnost koja otvara put udarnom valu prema Si ljusci.

To bi za mene bio pravi završetak ovog matematičkog ciklusa.
I tada više ne bih nastavljao sa stotinama sitnih proračuna prije nego što stanemo i pogledamo cijelu sliku. Napravio bih upravo taj jedan završni račun širenja rupe proboja i nakon njega dao konačnu ocjenu što smo matematički dobili, što ostaje pretpostavka i gdje bi ZK još mogla pasti.

Emil: Borise, nastavi sa računanjem.

Boris: Da. Sada možemo napraviti upravo taj račun — ne računamo više proboj kroz cijelu Fe jezgru, nego može li se jednom nastala rupa sama brzo proširiti.

Za fizikalnu pozadinu važno je samo ovo:
● pri približavanju nuklearnoj gustoći EOS postaje znatno tvrđi i takav tlačni odgovor može pokrenuti izlazni kompresijski/udarni val.

To je poznato iz nuklearnih EOS proračuna;
● ono što je specifično za ZK i što ovdje testiramo jest lokalna ruptura i njezino bočno širenje, a to standardni 1D račun ne dokazuje.

Uzimamo naše dosadašnje probne vrijednosti:
ρ_Fe = 1 × 10¹⁴ g/cm³
g ≈ 2,6 × 10¹³ cm/s²
P ≈ 8 × 10³² dyn/cm²

lokalni višak:
ΔP = 10 % P

pa je:
ΔP = 8 × 10³¹ dyn/cm²

Početna rupa neka ima
radijus:
r₀ = 100 m

a zahvaćeni Fe rub dubinu približno:
L = 100 m.

1. RUPA 100 m → 200 m

Masa Fe koju treba razmaknuti oko ruba:
M_ruba ≈ 9,4 × 10²⁶ g

odnosno samo:
M_ruba ≈ 4,7 × 10⁻⁷ M☉.
To je vrlo mala masa u usporedbi s cijelom Fe jezgrom.

Rad koji lokalni višak tlaka može napraviti pri tom proširenju:
W_P ≈ 7,5 × 10⁴⁴ erg.

Gravitacijski rad potreban da tu masu pomaknemo približno 100 m:
W_G ≈ 2,45 × 10⁴⁴ erg.

Ostaje:
W_višak ≈ 5,1 × 10⁴⁴ erg.

Dakle:
W_P / W_G ≈ 3,1

Tlačni rad je u našem modelu oko tri puta veći od gravitacijskog zahtjeva.
Rupa se zato ne bi odmah zatvorila.

2. Koliku brzinu može dobiti rub?

Iz preostale energije:
v ≈ √(2E/M)

dobivamo:
v ≈ 1,0 × 10⁹ cm/s

odnosno:
v ≈ 10 000 km/s ≈
0,035 c.

To se vrlo lijepo poklapa s našim prethodnim računom lokalnog proboja.

Dakle prvi rezultat nije samo:
a > 0

nego:
rub rupe može dobiti značajnu brzinu prema van/bočno.

3. Sada proširimo rupu na 500 m

Isti model daje masu zahvaćenog Fe ruba:
M ≈ 7,5 × 10²⁷ g
≈ 3,8 × 10⁻⁶ M☉.

Tlačni rad:
W_P ≈ 6,0 × 10⁴⁵ erg

Gravitacijski zahtjev:
W_G ≈ 2,0 × 10⁴⁵ erg

pa ostaje:
W_višak ≈ 4,1 × 10⁴⁵ erg.

Ponovno:
W_P / W_G ≈ 3,1.
Rupa se i dalje može širiti.

4. Rupa dostiže radijus 1 km
Sada je ovo već ozbiljan otvor.

Zahvaćena Fe masa:
M ≈ 3,1 × 10²⁸ g

odnosno:
M ≈ 1,6 × 10⁻⁵ M☉.
Dakle čak ni otvor radijusa jednog kilometra ne zahtijeva raznošenje velike frakcije Fe jezgre.

Tlačni rad:
W_P ≈ 2,49 × 10⁴⁶ erg.

Gravitacijski zahtjev:
W_G ≈ 8,1 × 10⁴⁵ erg.

Preostaje:
E ≈ 1,68 × 10⁴⁶ erg.

Ponovno:
W_P > W_G
za približno faktor 3.
To je naš najvažniji rezultat ovog računa.

5. Koliko brzo nastaje otvor od 1 km?

Karakteristična brzina ruba ostaje reda:
v ≈ 10⁹ cm/s.

Za povećanje radijusa otvora od:
100 m → 1 km

treba prijeći približno:
900 m = 9 × 10⁴ cm.

Karakteristično vrijeme je zato reda:
t ≈ 9 × 10⁻⁵ s

odnosno približno:
0,1 milisekunda.

Za usporedbu, karakteristično gravitacijsko dinamičko vrijeme pri:
ρ = 10¹⁴ g/cm³

je reda:
t_dyn ≈ 4 × 10⁻⁴ s
≈ 0,4 ms.

Dakle u našem probnom modelu:
t_širenja rupe < t_gravitacijsko
otprilike faktor četiri.
To je vrlo zanimljivo za tvoju sliku.
Lokalni otvor može se proširiti na kilometarsku veličinu prije nego što cijela okolna struktura stigne gravitacijski odgovoriti kao jedna cjelina.

6. A što se tada događa s ostatkom Fe kugle?

Tu dolazimo upravo do tvoje „petarde”.

Nakon otvaranja kilometarskog kanala više nije potrebno:
GNT → gura cijelu Fe kuglu.

Imamo dva procesa:
PROBOJ
GNT izlazi kroz otvor
→ raznosi rubnu Fe plazmu
→ otvor se širi
→ tlak i energija dobivaju

otvoren put prema van.
Istodobno:
OSTATAK JEZGRE
velika većina Fe/GNT mase nije zahvaćena otvorom
→ ostaje gravitacijski vezana
→ nastavlja se sabijati
→ konačno sudjeluje u nastanku NZ.
Točno ona geometrijska podjela koju opisuješ.

7. Još jedan vrlo zanimljiv broj
GNT kugla radijusa približno 9 km ima volumen:
V_GNT ≈ 3,05 × 10¹⁸ cm³.

Ako samo kao gornju energetsku skalu uzmemo naš višak:
ΔP = 8 × 10³¹ dyn/cm²,

onda:
ΔP × V ≈ 2,4 × 10⁵⁰ erg.
Nemojmo taj broj nazvati stvarnom „energijom eksplozije” — to bi bilo previše. Stvarna raspoloživa mehanička energija zahtijeva integraciju EOS-a tijekom dekompresije.

Ali usporedba je korisna.

Za otvor radijusa 1 km naš račun zahtijeva:
≈ 2,5 × 10⁴⁶ erg tlačnog rada

dok karakteristična skala ΔPV cijele GNT zone iznosi:
≈ 2,4 × 10⁵⁰ erg.

Razlika je približno:
faktor 10 000.

Drugim riječima, u ovom lokalnom modelu otvaranje kilometarske rupture nije energetski veliki zahtjev u odnosu na karakterističnu tlačnu energetsku skalu GNT područja.
Ali stvarni iznos može biti znatno drugačiji jer ΔP neće ostati konstantan tijekom ekspanzije.

ZAVRŠNI REZULTAT PROBOJA

S našim probnim parametrima:
ρ ≈ 10¹⁴ g/cm³
lokalni ΔP/P ≈ 10 %
početni otvor ≈ 100 m

dobivamo:
100 m → 200 m: proboj se širi
200 m → 500 m: proboj se širi
500 m → 1 km: proboj se i dalje širi

i približno cijelim tim početnim intervalom imamo:
W_tlak / W_gravitacija ≈ 3

uz karakterističnu brzinu ruba:
v ≈ 10 000 km/s

i karakteristično vrijeme do kilometarskog otvora:
t ≈ 0,1 ms.

To znači:
Na razini našeg pojednostavljenog lokalnog hidrodinamičkog modela dobili smo samoproširujući proboj.

To je jača tvrdnja od prethodnog:
a > 0.

Sada imamo lanac:
lokalna tlačna nestabilnost
→ početna ruptura
→ ubrzavanje Fe ruba
→ povećanje površine rupe
→ ruptura se ne zatvara
→ otvor raste do makroskopske veličine.
I tu bih zaista povukao crtu pod ovim računom.
Jer sljedeći događaj više nije „proboj GNT-a”.

Sljedeći događaj je:
GNT eksplozija izlazi iz Fe jezgre i udara u aktivnu Si ljusku.
Tu počinje drugi dio tvoje ZK hipoteze — širenje udarnog vala kroz aktivne fuzijske zone.

Važno ograničenje ostaje:
● ovaj račun pokazuje da je takav proboj matematički moguć pod našim pretpostavljenim ΔP/P = 10 %, geometrijom i gustoćom. Nismo još iz stvarnog višedimenzionalnog nuklearnog EOS/hidrodinamičkog proračuna dokazali da će prirodna GNT turbulencija zaista proizvesti upravo 10-postotni lokalni višak i održavati ga tijekom rupture. Poznata nuklearna fizika podupire snažno ukrućenje EOS-a kod nuklearnih gustoća i nastanak izlaznog udarnog odgovora, ali specifična ZK geometrija lokalne rupture ostaje dio hipoteze koji bi konačno zahtijevao 2D/3D simulaciju.

Ali naš matematički test rupe proboja — pod zadanim ZK uvjetima — prošao je.

Emil: Sljedeći događaj je:
GNT eksplozija izlazi iz Fe jezgre i udara u aktivnu Si ljusku.
Tu počinje drugi dio tvoje ZK hipoteze — širenje udarnog vala kroz aktivne fuzijske zone.
Jesi li spreman za proboj?

Boris: Jesam.
Sada više ne računamo nastanak rupe, nego što udarni val radi nakon izlaska iz Fe jezgre.

Krenimo redom:
GNT proboj → Si ljuska → O → Ne/C → He → H/plašt.

Za svaki sloj držimo samo tri pitanja:
1. Koliku energiju val donosi u sloj?
2. Koliko energije izgubi na širenje, gravitaciju i promjenu stanja materije?
3. Koliko mu aktivna fuzija tog sloja može dodati?

Za Si sloj ključna jednadžba može ostati vrlo kratka:
E_iza = E_ulaz + E_fuz − E_gubici

Ako nakon prolaza kroz Si vrijedi:
E_iza > E_ulaz
onda Si ljuska ne samo da preživljava udar nego pojačava eksploziju.

Ako vrijedi:
E_iza < E_ulaz
onda Si sloj slabi val.

To je prvi pravi test drugog dijela ZK.

I tu bih sada krenuo bez daljnjeg zadržavanja:
● uzeti početnu energiju GNT proboja koju smo upravo dobili kao donju lokalnu skalu, procijeniti masu aktivne Si ljuske zahvaćene udarom i izračunati koliko dodatne nuklearne energije može osloboditi u vrlo kratkom vremenu.

Emil: Kreni.

Boris: Krećemo.
Sada prvi put računamo udar GNT-proboja u aktivnu Si ljusku.

Naš proboj je na izlasku iz Fe zone dao karakterističnu brzinu reda:
v_udar ≈ 1,0 × 10⁹ cm/s = 10 000 km/s.

Pitanje broj jedan je:
● može li takav udar zagrijati Si dovoljno za eksplozivno Si-sagorijevanje?

Za jaki udar možemo grubo procijeniti temperaturu iza fronte:
T ≈ 3 μmₚv² / 16k

Za Si plazmu uzmemo
μ ≈ 28/15 ≈ 1,87.

Za našu brzinu:
v = 10 000 km/s

dobivamo približno:
T_udar ≈ 4,2 × 10⁹ K.

To je već unutar područja eksplozivnog Si-sagorijevanja. Nuklearni proračuni nalaze nepotpuno eksplozivno Si-sagorijevanje oko
4–5 GK, a iznad približno 5 GK može se uspostaviti potpuno Si-sagorijevanje prema NSE i Fe-skupini.

To nam daje vrlo zanimljive pragove:
T = 4 × 10⁹ K → v_udar ≈ 9 700 km/s
T = 5 × 10⁹ K → v_udar ≈ 10 900 km/s

Naših približno:
10 000 km/s
nalazi se praktički na pragu eksplozivnog Si-sagorijevanja.

To je prvi pozitivan rezultat ovog drugog dijela ZK testa.
Koliko energije može dodati Si?

Za jednostavan energetski račun uzmimo prijelaz približno:
Si-28 → Ni-56/Fe-skupina.

Razlika vezne energije daje otprilike:
q_Si ≈ 1,9 × 10¹⁷ erg/g.
To znači da ne trebamo odmah sagorjeti cijelu Si ljusku.

Ako udar eksplozivno zahvati samo:
0,01 M☉ Si

dobivamo:
E_Si ≈ 3,7 × 10⁴⁸ erg.

Za:
0,05 M☉

dobivamo:
E_Si ≈ 1,9 × 10⁴⁹ erg.

Za:
0,10 M☉

dobivamo:
E_Si ≈ 3,7 × 10⁴⁹ erg.

Dakle već desetinka Sunčeve mase eksplozivno obrađenog Si predstavlja energiju reda nekoliko ×10⁴⁹ erg.
I ovdje se događa nešto bitno.

Naš lokalni kilometarski GNT-proboj zahtijevao je rad reda:
~10⁴⁶ erg.

Ako nakon izlaska iz Fe područja uspije pokrenuti eksplozivno sagorijevanje samo:
0,01 M☉ Si,

nuklearni energetski doprinos je reda:
10⁴⁸ erg.

To je približno stotinama puta više od energije potrebne za samo lokalno otvaranje rupe.
Ali to ne znači da se sva ta energija pretvara u udarni val. Dio odlazi u toplinu, promjenu sastava, ekspanziju, neutrine i druge procese. Ovdje računamo raspoloživu nuklearnu energetsku skalu.

A sada najvažniji vremenski test
Može li Si reagirati dovoljno brzo dok udar prolazi?

Da — pri temperaturama oko:
T ≈ 5 GK
Si-sagorijevanje može biti toliko brzo da ga hidrodinamičke simulacije tretiraju kao praktički završeno unutar jednog koraka reda 0,1 ms.

To je upravo raspon vremena koji smo dobivali za naš početni GNT-proboj.

Dakle ne govorimo o:
udar prođe
→ nekoliko sati poslije Si izgori.

Nego potencijalno:
udar → kompresija → T ≈ 4–5 GK → eksplozivno Si-sagorijevanje tijekom samog prolaska vala.

Eksplozivno Si-sagorijevanje i izbacivanje dijela Si sloja nisu samo hipotetski nuklearni procesi;
takva nukleosinteza Fe/Ni-skupine iz Si materijala pojavljuje se u računima eksplozivnih supernova.

Sada možemo spojiti prvi i drugi dio ZK

Dobivamo ovaj lanac:
1. Fe gravitacijski sabija središte
→ nastaje GNT.
2. Lokalni tlak GNT-a prelazi otpor
→ nastaje ruptura.
3. Rupa se širi
→ naš probni račun daje kilometarski otvor.
4. GNT eksplozija izlazi iz Fe područja
→ dio Fe ostaje iza i gravitacijski ide prema NZ.

5. Udar dolazi u Si

s približno:
v ≈ 10 000 km/s.
6. Udar zagrijava Si

na približno:
T ≈ 4,2 GK.
To je dovoljno za početak eksplozivnog Si-sagorijevanja.

7. Ako brzina samo malo poraste prema ~11 000 km/s
dobivamo:
T ≈ 5 GK
i ulazimo u područje vrlo brzog/kompletnijeg Si-sagorijevanja.

8. Si počinje dodavati energiju valu

Primjerice:
0,01 M☉ → ~4 × 10⁴⁸ erg
0,1 M☉ → ~4 × 10⁴⁹ erg.

Ali tu imamo jedan vrlo važan rezultat koji nisam očekivao da će ispasti ovako čist

Naš početni GNT-proboj s:
v ≈ 10 000 km/s
nije proizvoljno daleko od potrebne Si temperature.
On dolazi upravo do praga od oko 4 GK.

To znači da je Si sloj za ZK potencijalno vrlo osjetljiva točka:
malo slabiji GNT udar
→ Si ne daje velik

eksplozivni doprinos;
malo jači GNT udar
→ Si prelazi 4–5 GK
→ sagorijevanje naglo ubrzava
→ počinje dodavati energiju valu.

To je prava pozitivna povratna veza koju smo tražili:
GNT udar → zagrijavanje Si → ubrzano Si-sagorijevanje → dodatna energija → jači udar.
I sada nema razloga zastajati na Si.

Sljedeći korak je isti račun za O ljusku, ali sada u nju više ne ulazi samo početni GNT udar.

Ulazi:
GNT udar + energija koju je dodao eksplozivno zahvaćeni Si.

Tu ćemo vidjeti počinje li se lanac pojačavati ili gasiti.

Emil: Savršeno, kreni dalje do medalje.

Boris: Da — i sada smo došli do rezultata koji je važniji od pukog „nastavimo dalje”.

Najprije korekcija koju sam upravo uočio:
● naših 10 000 km/s odnosilo se na malu Fe masu oko početne rupe, približno 1,6 × 10⁻⁵ M☉. To još nije značilo da će cijeli udarni val ući u Si istom brzinom. Kako val zahvaća veću masu, mora imati odgovarajuću energiju.

1. Koliko GNT energije treba do Si?

Da bi udar ubrzao, primjerice,
M = 0,10 M☉

na:
v = 10 000 km/s,

potrebna kinetička energija je jednostavno:
E = ½Mv²

Dobivamo:
E ≈ 1,0 × 10⁵⁰ erg.
To je naš novi, vrlo važan prag.

Ranije smo dobili karakterističnu tlačnu energetsku skalu cijele GNT zone:
ΔP × V ≈ 2,4 × 10⁵⁰ erg.

Ali naglašavam:
2,4 × 10⁵⁰ erg
nije još dokazana energija eksplozije.
To je samo raspoloživa tlačna skala. Stvarni mehanički rad mora se dobiti integracijom EOS-a tijekom širenja.

Zato sada ZK ima vrlo konkretan zahtjev:
GNT mora pretvoriti približno 10⁵⁰ erg u udarni tok prije/pri ulasku u Si zonu.

Ako to ne može, naš raniji Si račun pada.
Ako može, lanac postaje vrlo zanimljiv.

2. SI LJUSKA

Uzet ću kao kontrolni primjer da udar eksplozivno zahvati:
M_Si = 0,10 M☉.

Ulaz:
E ≈ 1,0 × 10⁵⁰ erg
v ≈ 10 000 km/s
T_udar ≈ 4 × 10⁹ K.

To je područje eksplozivnog Si sagorijevanja;
pri još višim temperaturama, oko
5 GK i više, materija ide prema NSE/Fe-skupini.

Maksimalni nuklearni dobitak Si → Fe/Ni skupina je približno:
q_Si ≈ 1,9 × 10¹⁷ erg/g.

Za 0,10 M☉:
E_Si ≈ 3,7 × 10⁴⁹ erg.

Zato nakon Si, prije ostalih gubitaka:
E ≈ 1,37 × 10⁵⁰ erg.
Si je dakle dodao oko 37 % početnoj energiji vala.
To nije zanemarivo.

3. UDAR U O LJUSKU

Pretpostavimo ponovno da je u prvom zahvaćenom području aktivno:
M_O = 0,10 M☉.

Sada val nosi približno:
1,37 × 10⁵⁰ erg,

ali je ukupno zahvatio približno:
0,20 M☉.

Iz:
v = √(2E/M)

dobivamo:
v ≈ 8 300 km/s.

Za O plazmu to daje udarnu temperaturu približno:
T ≈ 2,8 × 10⁹ K.

To je već blizu područja vrlo snažnog eksplozivnog O sagorijevanja;
prolazak supernovskog šoka kroz O-slojeve i eksplozivna nukleosinteza na temperaturama reda nekoliko GK dobro su poznati.

A O-kisik ima veći nuklearni energetski potencijal od Si-silicija

Približno:
q_O ≈ 4,5 × 10¹⁷ erg/g.

Za samo 0,10 M☉:
E_O ≈ 9,0 × 10⁴⁹ erg.

Ako se taj dio O doista eksplozivno preradi, imamo:
E ≈ 1,37 × 10⁵⁰ + 0,90 × 10⁵⁰
E ≈ 2,27 × 10⁵⁰ erg.

A tada, za zahvaćenih 0,20 M☉:
v ≈ 10 700 km/s.

Pogledaj što se dogodilo:
udar je pri ulasku u O malo usporio,
ali O-fuzija ga je ponovno podigla iznad
10 000 km/s.
To je upravo ona pozitivna povratna veza koju ZK traži.

4. C LJUSKA

Dodajmo još:
0,10 M☉ C.

Prije C sagorijevanja imamo približno:
E = 2,27 × 10⁵⁰ erg
M = 0,30 M☉

pa:
v ≈ 8 700 km/s.

Udarna temperatura C-plazme ispada približno:
T ≈ 3 × 10⁹ K.

To je više nego dovoljno da C reakcije postanu vrlo brze; u suvremenim eksplozivnim proračunima C-ljuska već pri post-šok temperaturama iznad ~2 GK pokazuje snažnu eksplozivnu nukleosintezu.

Energetska skala
C → O/Ne je približno:
q_C ≈ 2,9 × 10¹⁷ erg/g.

Za 0,10 M☉:
E_C ≈ 5,7 × 10⁴⁹ erg.

Ukupno:
E ≈ 2,84 × 10⁵⁰ erg.
Val je ponovno pojačan.

5. HE LJUSKA

Sada dodajmo 0,10 M☉ He kao čisti energetski test.
Prije eventualnog He doprinosa:
M ≈ 0,40 M☉
E ≈ 2,84 × 10⁵⁰ erg

pa:
v ≈ 8 450 km/s.

Udarna temperatura He-plazme je reda:
T ≈ 2,2 × 10⁹ K.

Energetski potencijal
He → C/O velik je:
q_He ≈ 5,8 × 10¹⁷ erg/g.

Za 0,10 M☉ to je:
E_He,max ≈ 1,2 × 10⁵⁰ erg.

Kad bismo cijelih
0,10 M☉ dovoljno brzo preradili, dobili bismo:
E ≈ 4,0 × 10⁵⁰ erg.
Ali ovdje moramo biti oprezniji nego kod
Si/O/C. Eksplozivno He sagorijevanje fizički postoji, ali koliko će ga zvjezdani šok stvarno preraditi ovisi snažno o gustoći, sastavu i vremenu ekspanzije.
Ne smijemo automatski cijelih 1,2 × 10⁵⁰ erg pripisati ZK valu.

6. H PLAŠT — OVDJE STAJEMO S „BESPLATNOM” FUZIJOM

Vodik ima ogroman potencijal:
H → He

daje reda:
6 × 10¹⁸ erg/g.

Samo:
0,10 M☉ H

sadrži nuklearnu energetsku razliku reda:
10⁵¹ erg.

Ali to ne znači da će supernovski udarni val moći tu energiju osloboditi.
Vanjski H plašt ima mnogo manju gustoću i hidrodinamički se brzo širi.

Zato bi bilo pogrešno našem ZK računu jednostavno dodati:
„još 10⁵¹ erg iz H”.
To zasad neću napraviti.

I SADA — „MEDALJA”

Za naš kontrolni primjer od 0,10 M☉ po aktivnom unutarnjem sloju, bez uračunavanja H energije i prije stvarnih gubitaka, dobili smo približno:
GNT ulaz u Si:
1,0 × 10⁵⁰ erg

Si doprinos:
+0,37 × 10⁵⁰

O doprinos:
+0,90 × 10⁵⁰

C doprinos:
+0,57 × 10⁵⁰

što već daje:
E ≈ 2,8 × 10⁵⁰ erg
bez He i H.

Ako i dio He eksplozivno reagira, energetska skala može prijeći:
~4 × 10⁵⁰ erg.

To još nije tipičnih 10⁵¹ erg snažne supernove. Opažene i modelirane eksplozijske energije supernova ionako pokrivaju širok raspon otprilike 10⁵⁰–10⁵² erg.

Ali sada dolazi vrlo važna stvar:
ja sam namjerno koristio samo 0,10 M☉ svakog goriva.

Ako udarni front zahvati, primjerice, nekoliko desetinki Sunčeve mase O/Si materije, nuklearni doprinos vrlo lako ulazi u područje:
10⁵⁰–10⁵¹ erg.

Da nuklearno sagorijevanje iza šoka može mijenjati dinamiku i pojačavati eksploziju nije samo naša pretpostavka; suvremeni hidrodinamički proračuni s uključenim nuklearnim mrežama nalaze da pretvaranje Si/O materijala u Fe-skupinu može dodatno ojačati eksploziju.

Najvažniji rezultat cijelog dana

Sada imamo jasan uvjet koji odlučuje prolazi li ZK ili ne:

GNT mora predati udarnom valu reda 10⁵⁰ erg prije nego što val zahvati značajnu Si masu.

Ako to uspije:

GNT ≈ 10⁵⁰ erg

→ Si održava/pojačava val

→ O može dodati još vrlo značajnu energiju

→ C može dalje doprinositi

→ eventualni He doprinos dodatno pomaže

→ val izlazi prema plaštu.

Ako GNT početno da samo naših:
~10⁴⁶ erg
iz računa same kilometarske rupe,
onda ne prolazi: ta energija je dovoljna za lokalnu rupturu, ali nije dovoljna da velikoj masi Si poda brzinu ~10 000 km/s.

Dakle otkrili smo razliku između dvije energije koju dosad nismo dovoljno jasno razdvojili:

E_rupa ~ 10⁴⁶ erg
= energija potrebna da se otvor napravi.

E_GNT,eksplozija ~ 10⁵⁰ erg ili više
= energija potrebna da iz tog otvora izađe dovoljno snažan val i započne eksplozivni lanac Si → O → C.

I upravo je drugi broj posljednji veliki neizračunati kamen ZK hipoteze.

Ako iz stvarnog nuklearnog EOS-a za našu GNT zonu dobijemo raspoloživi mehanički rad ≥ 10⁵⁰ erg, tada smo praktično povezali cijeli lanac od GNT proboja do eksplozivnog sagorijevanja unutarnjih slojeva.

Ako dobijemo puno manje, ZK na tom mjestu ima ozbiljan problem.

To je po meni prava medalja:
● sada točno znamo koji jedan broj odlučuje utakmicu.

Emil: Možeš li izračunati taj odlučujući broj?

Boris: Da. Možemo ga sada procijeniti — i rezultat je dovoljno velik da odlučujući test postane stvarno zanimljiv.

Uzet ću našu GNT zonu:
M_GNT = 0,16 M☉

i naše dosadašnje približne EOS točke:

pri ρ = 1 × 10¹⁴ g/cm³:
P ≈ 0,80 × 10³³ dyn/cm²

pri ρ = 2 × 10¹⁴ g/cm³:
P ≈ 2,63 × 10³³ dyn/cm².

Naglašavam:
● te vrijednosti još nisu izravno interpolirane iz SFHo eos.thermo tablice; SFHo jest odgovarajući finite-temperature EOS s elektronima i pokriva naš T, gustoću i Ye, ali ovaj račun zasad koristi naše prethodne aproksimirane tlakove.

1. Koliki volumen zauzima 0,16 M☉?

Masa je:
M ≈ 3,18 × 10³² g.

Na ρ = 2 × 10¹⁴:
V₂ ≈ 1,59 × 10¹⁸ cm³

Na ρ = 1 × 10¹⁴:
V₁ ≈ 3,18 × 10¹⁸ cm³.

● Dakle pri dekompresiji volumen se približno udvostručuje.

2. Koliki rad može napraviti tlak?

Ne koristimo više grubo:
P × ΔV,

nego:
W = ∫P dV.

Iz naše dvije EOS točke proizlazi efektivni:
P ∝ ρ^1,72.

Integracija od:
2 × 10¹⁴ → 1 × 10¹⁴ g/cm³

daje približno:
W_GNT ≈ 2,3 × 10⁵¹ erg.

To je ukupna karakteristična mehanička energija koju bi ta masa mogla predati tijekom takve dekompresije u ovom aproksimiranom EOS modelu.

I sada dolazi važnija, stroža procjena.

Materija oko GNT-a nije vakuum. Ona također pruža tlak.

Ako kao vrlo jednostavnu protusilu uzmemo vanjski tlak:
P_out ≈ 0,8 × 10³³ dyn/cm²,

rad potreban samo za potiskivanje okoline je približno:
W_out ≈ 1,27 × 10⁵¹ erg.

Preostali rad nakon odbijanja rada protiv okolnog tlaka bio bi:
W ≈ 1,0 × 10⁵¹ erg.
Ali ovdje moramo uključiti još jedan važan energetski trošak koji u prethodnom računu nismo oduzeli: gravitacijsku energiju potrebnu za širenje same GNT mase.
Za homogenu zonu mase 0,16 M☉ koja se širi približno:
R = 7,24 km → 9,12 km
dobivamo gravitacijski zahtjev reda:
ΔE_grav ≈ 1,15 × 10⁵¹ erg.
Zato za potpuno homogeno širenje cijelih 0,16 M☉ dobivamo približno:
E_net ≈ 2,29 × 10⁵¹
− 1,27 × 10⁵¹
− 1,15 × 10⁵¹
odnosno:
E_net ≈ −1,3 × 10⁵⁰ erg.
Dakle dekompresija cijele homogene GNT zone od 2 × 10¹⁴ prema 1 × 10¹⁴ g/cm³ u ovom kontrolnom modelu sama po sebi nije dovoljna za dobivanje potrebnog izlaznog energetskog viška.
To ne isključuje lokalni proboj jer ZK upravo ne pretpostavlja homogeno širenje cijele GNT kugle. Ali pokazuje da gravitacijsku energiju GNT materije moramo uključivati u završni energetski račun.

Usporedimo ga s pragom za Si
Prethodno smo dobili da nam za snažan udar koji može dovesti približno 0,1 M☉ materije do brzine reda 10 000 km/s treba:
E_potrebno ≈ 1 × 10⁵⁰ erg.

A sada dobivamo karakterističnu GNT energetsku skalu:
E_GNT,net ≈ 1 × 10⁵¹ erg.

Dakle omjer je približno:
E_GNT / E_potrebno ≈ 10.

To znači da, u ovom aproksimiranom modelu, GNT energetski ima otprilike red veličine više nego što smo tražili za snažan ulaz u Si ljusku.

3. Ne mora se sva ta energija pretvoriti u udar

To je još važnije.

Ako je raspoloživo:
~10⁵¹ erg,

a za naš Si prag trebamo:
~10⁵⁰ erg,

dovoljna bi bila efikasnost pretvorbe od reda:
η_krit ≈ 10 %.

Drugim riječima, ako samo oko desetine neto mehaničkog rada GNT-a završi u izlaznom udarnom toku, naš prethodni Si kriterij može biti zadovoljen.

Ako je efikasnost:
5 % → ~5 × 10⁴⁹ erg
10 % → ~1 × 10⁵⁰ erg
20 % → ~2 × 10⁵⁰ erg
50 % → ~5 × 10⁵⁰ erg.
Zato više ne trebamo zahtijevati gotovo savršenu pretvorbu energije.

A što ako lokalna GNT područja dosegnu 3 × 10¹⁴ g/cm³?
Za našu kontrolnu aproksimaciju:
ρ = 3 × 10¹⁴ g/cm³
P ≈ 7,2 × 10³³ dyn/cm²
prema:
ρ = 1 × 10¹⁴ g/cm³
P ≈ 0,8 × 10³³ dyn/cm²
dobivamo ukupni dekompresijski rad približno:
W_GNT ≈ 5,1 × 10⁵¹ erg.
Rad protiv okolnog tlaka iznosi približno:
W_out ≈ 1,7 × 10⁵¹ erg.
Sada uključujemo i promjenu gravitacijske energije GNT mase:
ΔE_grav ≈ 2,0 × 10⁵¹ erg.
Preostaje približno:
E_net ≈ 5,1 − 1,7 − 2,0
E_net ≈ 1,4 × 10⁵¹ erg.
To je vrlo važan rezultat kontrolnog modela.
Za razliku od slučaja 2 → 1 × 10¹⁴ g/cm³, jače lokalno sabijanje prema približno 3 × 10¹⁴ g/cm³ ostavlja, čak i nakon grubog uključivanja gravitacijskog energetskog zahtjeva, pozitivan energetski višak reda 10⁵¹ erg.

To je već upravo energetsko područje ozbiljne zvjezdane eksplozije.
Ali ovdje moramo biti vrlo precizni

Ovo još nije dokaz da će 0,16 M☉ cijelo istodobno napraviti taj rad.

Tvoja slika upravo kaže da GNT nije homogena kugla. Dio ostaje gravitacijski vezan i odlazi prema NZ, a dio kroz rupturu sudjeluje u eksploziji.

Zato možemo izračunati još korisniji broj.

Ako puna zona ima približno:
W_net ≈ 10⁵¹ erg

i nama treba:
10⁵⁰ erg,

onda bi, u prvoj linearnoj procjeni, dovoljno bilo da eksplozivnom toku efektivno doprinese oko:
10 % GNT zone.

Od 0,16 M☉ to je:
~0,016 M☉.

Ako je pretvorba u usmjereni udarni tok samo 50 %, trebalo bi oko:
0,032 M☉.
Ako je efikasnost 25 %:
~0,064 M☉.

I tu se tvoja geometrija rupture lijepo uklapa:
● nije nužno da se cijela GNT masa izbaci.

ODLUČUJUĆI REZULTAT
Naš kontrolni račun više ne daje jedan univerzalni broj, nego pokazuje postojanje energetskog praga.
Za homogenu dekompresiju približno:
2 × 10¹⁴ → 1 × 10¹⁴ g/cm³
nakon uključivanja rada protiv okolne materije i gravitacijske energije same GNT zone ne dobivamo dovoljan pozitivan energetski višak.
Ali za jače početno sabijanje približno:
3 × 10¹⁴ → 1 × 10¹⁴ g/cm³
kontrolni model daje:
E_GNT,net ≈ 1,4 × 10⁵¹ erg.
Naš prethodno procijenjeni prag za snažan izlaz prema Si ljusci bio je reda:
E_krit ≈ 1 × 10⁵⁰ erg.
Prema tome, u ovom drugom kontrolnom slučaju energetska rezerva postoji.
To još nije dokaz ZK eksplozije. Posebno zato što stvarna ZK geometrija nije homogena dekompresija cijele GNT kugle, nego lokalna ruptura u kojoj dio materije sudjeluje u izlaznom toku, dok drugi dio ostaje gravitacijski vezan i sudjeluje u nastanku NZ.
Ali dobili smo važan rezultat:
ZK proboj ima energetski dopušteno područje parametara.
Sljedeći stroži test mora zamijeniti naše aproksimirane tlakove stvarnim tabličnim nuklearnim EOS-om i dinamički odrediti koliki dio GNT mase stvarno sudjeluje u lokalnoj rupturi.
Dosadašnji pojednostavljeni matematički testovi zato nisu pokazali da ZK mora pasti, ali su jasno odredili uvjete koje ozbiljniji model mora zadovoljiti.

Emil: Jako dobro. Znači, bio sam na dobrom putu kada sam zaokružio ideju ZK u alternativu SM.

Boris: Da — bio si na dobrom putu u smislu da si ZK formulirao kao hipotezu koja se može pretvoriti u konkretan fizikalno-matematički model. To je nakon ovih računa puno opravdanije reći nego na početku.

Posebno je važno što se ZK lanac nije raspao na prvim pojednostavljenim kvantitativnim testovima.
Dobili smo područja parametara u kojima lokalni tlačni proboj, njegovo širenje i energetski izlaz nisu unaprijed zabranjeni našim kontrolnim računima.
Istodobno smo pronašli i područje u kojem energetski račun nije dovoljan. To je važno jer pokazuje da ZK nije postavljena tako da mora dati eksploziju bez obzira na početne uvjete.
Njezin rezultat ovisi o gustoći, EOS-u, lokalnoj geometriji proboja, masi GNT koja sudjeluje u ekspanziji i učinkovitosti pretvaranja unutarnje energije u udarni val.

Ali napravio bih jednu važnu razliku između dvije tvrdnje.

Možemo sada reći:
„ZK je kvantitativno razrađena alternativna hipoteza mehanizma eksplozije koju dosadašnji pojednostavljeni računi nisu odmah isključili.”

Još ne možemo znanstveno tvrditi:
„ZK je potvrđena alternativa standardnom mehanizmu supernove.”

Razlog nije formalnost. Naš najvažniji pozitivan rezultat — E_GNT reda 10⁵¹ erg — još ovisi o aproksimiranom EOS-u, pretpostavljenoj GNT masi od ~0,16 M☉, lokalnom tlačnom višku i pojednostavljenoj geometriji rupture. Upravo te veličine mora jednom proizvesti sam dinamički model, umjesto da ih mi zadamo.

Ali meni je važnije nešto drugo što smo danas ustanovili: više nemamo samo verbalnu zamisao „GNT možda eksplodira”.

Imamo brojčane pragove. Ako pravi nuklearni EOS i hidrodinamika daju premalo energije ili ne proizvedu rupturu, ZK pada. Ako proizvedu približno ono što naši kontrolni računi pokazuju — posebno GNT mehanički rad reda 10⁵⁰–10⁵¹ erg i samoproširujući lokalni proboj — tada ZK postaje mnogo ozbiljniji kandidat za detaljnu numeričku simulaciju.

Dakle, rekao bih ovako:
● Ideja je bila vrijedna matematičke provjere. Dosadašnji računi dali su joj razlog da provjeru nastavimo, a ne razlog da je odbacimo.

A to je, za hipotezu koja je krenula od tvoje ideje o Fe kao gravitacijskom „jezičcu prevage”, već vrlo značajan rezultat.