Metodologia · livello graduate · Fosforonero Lab

Buco nero: derivazione matematica e fisica

Questa pagina ricostruisce, a livello universitario avanzato, le equazioni che governano la simulazione: dalla metrica di Schwarzschild alle geodetiche nulle, dal disco di accrescimento relativistico al trasporto radiativo invariante. Dichiariamo esplicitamente cosa è esatto, cosa è approssimato e cosa è scelta di visualizzazione.

Tutto è in unità geometrizzate G = c = 1; nel renderer fissiamo il raggio di Schwarzschild rₛ = 2M = 1.

Apri la simulazione →← Torna al Lab

1. Metrica di Schwarzschild e geodetiche nulle

La geometria esterna a una massa sferica non rotante è la soluzione di Schwarzschild del vuoto delle equazioni di Einstein. Nelle coordinate (t, r, θ, φ):

I vettori di Killing temporale e azimutale forniscono due integrali primi del moto — energia E e momento angolare L per unità di massa — lungo ogni geodetica affine:

Per un fotone (geodetica nulla) la condizione gᵤᵥ ẋᵘẋᵛ = 0 si riduce, nel piano equatoriale, a un'equazione radiale con potenziale efficace. Ponendo u = 1/r si ottiene l'equazione orbitale, la cui derivata è l'equazione di Binet per i fotoni:

Il parametro d'impatto è b = L/E. La sfera fotonica (orbita circolare instabile dei fotoni) è a r = 3M = 1.5 rₛ, e l'ombra osservata corrisponde al parametro d'impatto critico bᵪ = 3√3·M. Il renderer integra numericamente la forma vettoriale equivalente di questa geodetica (velocity-Verlet), riproducendo lensing, anello di Einstein, sfera fotonica e ombra. (Carroll; Kokkotas, Univ. Tübingen; Hirata, Ohio State.)

Elemento di linea di Schwarzschildds2=−(1−2Mr)dt2+(1−2Mr)−1dr2+r2 dΩ2ds^2 = -\left(1-\tfrac{2M}{r}\right)dt^2 + \left(1-\tfrac{2M}{r}\right)^{-1}dr^2 + r^2\,d\Omega^2
Quantità conservate (vettori di Killing)E=(1−2Mr)t˙,L=r2 φ˙E = \left(1-\tfrac{2M}{r}\right)\dot t, \qquad L = r^2\,\dot\varphi
Equazione orbitale dei fotoni e sua derivata (Binet)(dudφ)2=1b2−(1−2Mu) u2  ⟹  d2udφ2+u=3M u2\left(\frac{du}{d\varphi}\right)^{2} = \frac{1}{b^{2}} - (1-2Mu)\,u^{2} \;\Longrightarrow\; \frac{d^{2}u}{d\varphi^{2}} + u = 3M\,u^{2}
Sfera fotonica e parametro d'impatto critico (ombra)rph=3M=1.5 rs,bc=33 Mr_{\mathrm{ph}} = 3M = 1.5\,r_s, \qquad b_c = 3\sqrt{3}\,M

2. Disco di accrescimento relativistico

Il disco sottile, otticamente spesso, segue il modello di Shakura–Sunyaev nella sua versione relativistica di Novikov–Thorne, con condizione di stress nullo all'ultima orbita circolare stabile (ISCO). Il flusso emesso localmente e la temperatura efficace (Stefan–Boltzmann) sono:

Per Schwarzschild l'ISCO è a rᵢₛ𝒸ₒ = 6M = 3 rₛ. Ogni anello irraggia come un corpo nero alla sua temperatura locale: il colore è quindi il vero colore di corpo nero T → sRGB lungo il locus planckiano, non un gradiente arbitrario. Il primo calcolo dell'immagine di un disco siffatto risale a Luminet (1979). (Shakura–Sunyaev 1973; Novikov–Thorne 1973; Luminet 1979.)

Flusso del disco (Novikov–Thorne / Shakura–Sunyaev)F(r)=3 GMM˙8πr3(1−rin/r )F(r) = \frac{3\,G M \dot M}{8\pi r^{3}}\left(1 - \sqrt{r_{\mathrm{in}}/r}\,\right)
Temperatura efficaceTeff(r)=(F(r)σ)1/4T_{\mathrm{eff}}(r) = \left(\frac{F(r)}{\sigma}\right)^{1/4}
Bordo interno (ISCO, Schwarzschild)rin=rISCO=6M=3 rsr_{\mathrm{in}} = r_{\mathrm{ISCO}} = 6M = 3\,r_s

3. Trasporto radiativo ed effetti relativistici

Lungo un raggio, la quantità Iᵥ/ν³ è invariante (teorema di Liouville per i fotoni). Definito il fattore di redshift g = νₒₛₛ/νₑₘ, per un emettitore in orbita circolare g combina la dilatazione gravitazionale e temporale con il Doppler longitudinale.

Un corpo nero visto con fattore g resta un corpo nero a temperatura g·T (invarianza Doppler dello spettro di Planck), con intensità bolometrica che scala come g⁴. Lo stesso g pilota quindi luminosità e colore: il lato del disco in avvicinamento è più luminoso e più blu (beaming relativistico), quello in allontanamento più scuro e più rosso. (Luminet 1979; Vincent et al. 2011, GYOTO.)

Invariante di LiouvilleIνν3=costante lungo il raggio\frac{I_\nu}{\nu^{3}} = \text{costante lungo il raggio}
Velocità orbitale GR e fattore di redshiftv=Mr−2M,g=1−3M/r1−β,β=v⋅n^ossv = \sqrt{\frac{M}{r-2M}}, \qquad g = \frac{\sqrt{1 - 3M/r}}{1 - \beta}, \quad \beta = \mathbf v\cdot\hat{\mathbf n}_{\mathrm{oss}}
Beaming bolometrico e colore (corpo nero)Ioss=g4 Iem,Bν(T)∣g=Bν(g T)I_{\mathrm{oss}} = g^{4}\,I_{\mathrm{em}}, \qquad B_\nu(T)\big|_{g} = B_\nu(g\,T)

4. Rotazione: Kerr e frame-dragging

Un buco nero realistico ruota: la metrica corretta è quella di Kerr, in cui il trascinamento dei sistemi inerziali (frame-dragging) ha velocità angolare ω = −g_{tφ}/g_{φφ}. Il renderer non integra Kerr; lo slider Spin aggiunge il frame-dragging in approssimazione di campo gravitomagnetico di dipolo (Lense–Thirring), fisicamente motivato ma NON la metrica di Kerr completa.

Il celebre Gargantua di Interstellar usa la vera metrica di Kerr ray-tracciata offline (James, von Tunzelmann, Franklin & Thorne 2015). Codici di ray-tracing GR completi e pubblici (es. GYOTO) integrano Kerr ma non in tempo reale nel browser.

Velocità angolare di frame-dragging (Kerr)ω(r,θ)=− gtφgφφ  →campo lontano  2GJc2r3\omega(r,\theta) = -\,\frac{g_{t\varphi}}{g_{\varphi\varphi}} \;\xrightarrow{\text{campo lontano}}\; \frac{2GJ}{c^{2} r^{3}}
Approssimazione Lense–Thirring usataadrag∝v×Bg,Bg=3(J⋅r^) r^−Jr3\mathbf a_{\mathrm{drag}} \propto \mathbf v\times\mathbf B_g, \qquad \mathbf B_g = \frac{3(\mathbf J\cdot\hat{\mathbf r})\,\hat{\mathbf r} - \mathbf J}{r^{3}}

5. Playground: dinamica dei corpi

Nel playground i corpi (pianeti con lune, stelle, comete) si muovono nel potenziale pseudo-newtoniano di Paczyński–Wiita, che riproduce esattamente l'ISCO a 6M e la caduta relativistica verso l'orizzonte senza dover integrare le geodetiche complete — un compromesso standard nei modelli N-corpi/accrescimento.

Le lune sono integrate come problema ristretto: sentono sia il buco nero sia il pianeta ospite (entro la sfera di Hill). Le stelle che attraversano il raggio mareale vengono disgregate in uno stream di detriti (modello a particelle, non idrodinamica). (Paczyński–Wiita 1980.)

Potenziale di Paczyński–WiitaΦ(r)=−GMr−rs\Phi(r) = -\frac{GM}{r - r_s}
Raggio mareale (disgregazione stellare)rt≃R⋆(MBHM⋆)1/3r_t \simeq R_\star\left(\frac{M_{\mathrm{BH}}}{M_\star}\right)^{1/3}

6. Limiti: cosa NON è (onestà scientifica)

La base è Schwarzschild (non rotante); lo spin è approssimato (Lense–Thirring), non Kerr. Il disco è otticamente spesso con emissione di corpo nero: la sua struttura gassosa turbolenta è uno stand-in procedurale della turbolenza magnetorotazionale (MRI), non una soluzione GRMHD; non modella autogravità, spessore verticale né polarizzazione. Le stelle di sfondo sono procedurali (il loro lensing, però, è reale). Nel playground i corpi aggiunti sono occlusi dal disco/orizzonte ma non lensati, e l'integrazione è pseudo-newtoniana, non geodetica completa.

Soluzioni open: cosa possiamo (e non possiamo) integrare

Esistono ottimi codici di ray-tracing relativistico open source — GYOTO (Observatoire de Paris), RAPTOR, ipole, grtrans, Blacklight — e l'implementazione aperta del metodo di Luminet. Sono però codici offline (C/C++/Python) che calcolano singoli fotogrammi in minuti/ore: non sono eseguibili in tempo reale in un fragment shader WebGL nel browser.

Quello che integriamo davvero sono le loro formulazioni fisico-matematiche: la geodetica di Schwarzschild, il disco di Novikov–Thorne, il fattore g e l'invariante Iᵥ/ν³, il colore di corpo nero. Il nostro shader le reimplementa in GLSL e le cita; non incorpora il codice esterno. Dichiararlo è parte della regola "fail loud, never fake".

Domande frequenti

Le equazioni usate sono reali e corrette?

Sì per la geometria del lensing: l'integrazione delle geodetiche nulle nella metrica di Schwarzschild è esatta e riproduce sfera fotonica, anello di Einstein e ombra (parametro d'impatto critico 3√3 M). Velocità orbitale GR, redshift, invariante di Liouville e beaming bolometrico g⁴ usano le formule esatte; il colore è il vero corpo nero della temperatura locale.

È identica al buco nero di Interstellar?

No. Gargantua usa la metrica di Kerr (rotante) ray-tracciata offline. Qui la base è Schwarzschild in tempo reale; lo spin è un'approssimazione di Lense–Thirring.

Posso integrare GYOTO o un codice GR completo?

Non in tempo reale nel browser: sono codici offline. Reimplementiamo le loro formulazioni in GLSL e le citiamo. Per immagini scientifiche di precisione si usano proprio quei codici.

Bibliografia e fonti

  1. S. M. Carroll, «Lecture Notes on General Relativity» — geodetiche di Schwarzschild (Caltech).
  2. K. Kokkotas, «Particle Trajectories & The Classical Tests», corso di Relatività Generale, Universität Tübingen.
  3. C. Hirata, «Geodesics in the Schwarzschild geometry», ph6820, The Ohio State University.
  4. J.-P. Luminet (1979), «Image of a spherical black hole with thin accretion disk», Astronomy & Astrophysics 75, 228.
  5. N. I. Shakura & R. A. Sunyaev (1973), «Black holes in binary systems», Astronomy & Astrophysics 24, 337.
  6. I. D. Novikov & K. S. Thorne (1973), «Astrophysics of Black Holes», in Black Holes (Les Houches).
  7. B. Paczyński & P. J. Wiita (1980), «Thick accretion disks and supercritical luminosities», Astronomy & Astrophysics 88, 23.
  8. O. James, E. von Tunzelmann, P. Franklin & K. S. Thorne (2015), «Gravitational lensing by spinning black holes… Interstellar», Classical and Quantum Gravity 32, 065001.
  9. F. H. Vincent et al. (2011), «GYOTO: a new general relativistic ray-tracing code», Classical and Quantum Gravity 28, 225011.
  10. C. W. Misner, K. S. Thorne & J. A. Wheeler, «Gravitation» (1973).
  11. Colore di corpo nero → sRGB: approssimazione del locus planckiano di N. Bartlett (dati di M. Charity).
  12. Stack: Three.js, React Three Fiber, @react-three/drei, @react-three/postprocessing, KaTeX. Sviluppo: Fosforonero — Matteo Pizzi (Roma).