Mathematical Foundations of Galactic Missing Mass (Fundamentos matemáticos da massa galáctica ausente): Disco, Esfera, Densidade, Potencial e Escalonamento Radial
TL;DR: O problema da massa ausente aparece quando a massa inferida a partir das curvas de rotação galáctica excede a massa diretamente observada em estrelas, gás e poeira. Matematicamente, isso requer a conexão da densidade da superfície em um disco, da densidade do volume em três dimensões, do potencial gravitacional, da aceleração radial e da massa contida.
1. Coordenadas radiais e geometria
Distinguimos duas geometrias:
- Geometria do disco: a matéria galáctica visível é distribuída principalmente em um fino disco rotativo.
- Geometria esférica: a massa escura ou ausente é geralmente modelada como um halo aproximadamente esférico.
Elemento da área do disco:
\[ dA = R\,dR\,d\phi \]
Elemento de volume esférico:
\[ dV = 4\pi r^2\,dr \]
O mesmo símbolo \(r\) é frequentemente usado para o raio galactocêntrico, mas o significado depende da geometria. Em um disco, \(R\) é um raio cilíndrico. Em um halo, \(r\) geralmente é um raio esférico.
2. Massa visível em um disco galáctico
O disco visível é frequentemente aproximado por uma densidade de superfície exponencial:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
A massa de um anel entre \(R\) e \(R+dR\) é:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disco}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
A massa cumulativa do disco visível é, portanto:
\[ M_{\rm disk}(R)=2\pi\int_0^R R’\Sigma_0e^{-R’/R_d}\,dR’ \]
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Em um raio grande:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
A massa do disco visível se aproxima de um valor finito.
3. Massa esférica e densidade de volume
Para uma distribuição esférica de massa, a densidade de volume \(\rho(r)\) determina a massa contida:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
A relação inversa é:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Essa relação é fundamental para o problema da massa ausente. Se a massa inferida cresce linearmente com o raio, então a densidade esférica correspondente diminui como \(1/r^2\).
4. Massa dinâmica do movimento circular
Para o movimento circular, a aceleração gravitacional é satisfatória:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Portanto:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Para uma curva de rotação plana:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Isso fornece a escala padrão:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definição de massa ausente
A massa ausente é a diferença entre a massa dinâmica e a massa visível:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Para um disco visível exponencial:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Para \(v(r)\approx v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Em um raio grande, o termo do disco satura, enquanto o termo dinâmico continua a crescer aproximadamente como \(r\).
6. Potencial gravitacional em 3D
O potencial gravitacional newtoniano gerado por uma massa pontual é:
\[ \Phi(r)=-\frac{GM}{r} \]
O campo gravitacional correspondente é a derivada radial do potencial:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Isso explica a relação entre \(1/r\) e \(1/r^2\): o potencial de uma massa localizada cai como \(1/r\), enquanto a força ou a aceleração cai como \(1/r^2\).
7. Equação de Poisson
A densidade de massa e o potencial gravitacional estão conectados pela equação de Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Na simetria esférica, isso se torna:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Essa equação relaciona três quantidades:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Potencial de uma densidade 3D estendida
Para uma distribuição geral de densidade 3D, o potencial é:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
O kernel \(1/|\mathbf{x}-\mathbf{x}’|\) é a origem matemática do potencial \(1/r\) em três dimensões.
9. Potencial de um disco fino
Para um disco fino com densidade de superfície \(\Sigma(R’)\), o potencial gravitacional no plano do disco pode ser escrito como:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’\,dR’\,d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
A distância entre um ponto de campo no raio \(R\) e um ponto de origem no raio \(R’\) é:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
A aceleração radial no disco é obtida pela diferenciação do potencial:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
A velocidade de rotação decorre disso:
\[ v^2(R)=R\,|g_R(R)| \]
10. Projeção de uma interação 3D no disco
Se uma interação se propagar em três dimensões, mas for avaliada no plano do disco, a projeção radial introduz um fator geométrico. Para dois pontos no disco separados pela distância \(d\), o fator de projeção radial é:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Assim, um kernel radial 3D genérico \(K(d)\), projetado no disco, aparece como:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Por exemplo, um núcleo semelhante a uma força newtoniana tem:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Um kernel 3D exponencial pode ser escrito como:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Núcleos exponenciais em três dimensões
Um fator radial puramente exponencial tem a forma:
\[ e^{-r/\lambda} \]
Na teoria de campo tridimensional, um potencial exponencialmente filtrado aparece com frequência na forma semelhante à de Yukawa:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
A força radial associada contém termos \(1/r^2\) e exponenciais:
\[ g_Y(r) = -\frac{d}{dr} \left( \frac{e^{-r/\lambda}}{r} \right) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \right) \]
Isso mostra por que o comportamento radial exponencial em 3D não é independente da geometria de \(1/r\). O exponencial controla a atenuação, enquanto \(1/r\) e \(1/r^2\) surgem da propagação tridimensional.
12. Leis de escala radial
O problema da massa ausente está fortemente ligado ao escalonamento radial. Várias leis radiais importantes aparecem repetidamente:
| Quantidade | Escala típica | Significado |
|---|---|---|
| Potencial da massa pontual | \(\Phi(r)\sim 1/r\) | Função verde 3D da gravidade |
| Força da massa pontual | \(g(r)\sim 1/r^2\) | Derivada de \(1/r\) |
| Velocidade de rotação plana | \(v(r)\sim constante\) | Observada em discos galácticos externos |
| Massa dinâmica | \(M(r)\sim r\) | Exigida pela rotação plana |
| Densidade do halo | \(\rho(r)\sim 1/r^2\) | Dá \(M(r)\sim r\) |
| Disco exponencial | \(\Sigma(R)\sim e^{-R/R_d}\) | O disco visível se desvanece rapidamente |
| Potencial 3D filtrado | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Atenuação exponencial mais espalhamento 3D |
13. Da densidade à curva de rotação
Para um halo esférico com densidade:
\[ \rho(r)\propto \frac{1}{r^2} \]
a massa contida é:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Então:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\aprox constante \]
\[ v(r)\aprox constante \]
Essa é a ponte matemática entre uma densidade de halo \(1/r^2\) e uma curva de rotação galáctica plana.
14. Da massa do disco à massa ausente
A massa do disco visível cresce rapidamente no início e depois satura:
\[ M_{\rm disk}(R)\rightarrow M_d \]
A massa dinâmica inferida a partir de uma curva de rotação plana continua crescendo:
\[ M_{\rm dyn}(r)\propto r \]
Portanto, a massa ausente se comporta aproximadamente como:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
Em um raio suficientemente grande:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Aviso matemático: disco e esfera não são intercambiáveis
A equação
\[ M(r)=\frac{v(r)^2r}{G} \]
é exato para a simetria esférica. Para um disco achatado, deve-se calcular o potencial integrando-o ao disco e, em seguida, derivar a aceleração radial. A expressão esférica é frequentemente usada como uma aproximação eficaz, especialmente quando se discute a massa necessária para suportar uma determinada curva de rotação.
16. Resumo das equações-chave
Densidade da superfície do disco:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Elemento de massa do disco:
\[ dM_{\rm disco}=2\pi R\Sigma(R)dR \]
Massa do disco visível:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Volume esférico de massa:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Densidade da massa fechada:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Massa dinâmica:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Massa ausente:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Potencial newtoniano:
\[ \Phi(r)=-\frac{GM}{r} \]
Equação de Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Integral de potencial 3D:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potencial de disco fino:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’dR’d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Núcleo radial projetado:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Potencial exponencial de tela 3D:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Conclusão
O problema da massa ausente é uma incompatibilidade matemática entre dois comportamentos radiais. O disco visível segue uma densidade de superfície exponencial e atinge uma massa cumulativa finita. A massa dinâmica inferida a partir de curvas de rotação aproximadamente planas cresce aproximadamente de forma linear com o raio. Se interpretada como um halo esférico, isso corresponde a uma densidade que diminui aproximadamente como \(1/r^2\). As equações de integração de disco, cascas esféricas, potencial gravitacional e projeção radial fornecem a linguagem matemática necessária para analisar essa incompatibilidade.