agujero negro

You are currently browsing articles tagged agujero negro.

Calculamos la curvatura escalar $latex R$ para la métrica de Kerr-Newman, es decir, para la solución analítica a las ecuaciones de Einstein en presencia de momento y de carga ($latex J neq 0$ y $latex Q neq 0$).

La métrica es:

$latex g = -frac{J^2+M^2 left(Q^2+r (-2 M+r)right)+J^2 text{Sin}[theta ]^2}{M^2 r^2+J^2 text{Cos}[theta ]^2} dt otimes dt + $

$latex + 2 frac{J M left(Q^2-2 M rright) text{Sin}[theta ]^2}{M^2 r^2+J^2 text{Cos}[theta ]^2} dt tilde{otimes} dvarphi + $

$latex + frac{M^2 r^2+J^2 text{Cos}[theta ]^2}{J^2+M^2 left(Q^2+r (-2 M+r)right)} dr otimes dr +$

$latex r^2+frac{J^2 text{Cos}[theta ]^2}{M^2} dtheta otimes dtheta +$

$latex + frac{left(frac{J^2}{M^2}+r^2right)^2 text{Sin}[theta ]^2}{r^2+frac{J^2 text{Cos}[theta ]^2}{M^2}}+frac{J^2 text{Sin}[theta ]^4}{M^2} dvarphi otimes dvarphi$.

Como ya hemos hecho con la métrica de Kerr, solo mostramos una componente de cada elemento calculado debido a su extrema complejidad (como para realizar los cálculos manualmente…):

$latex R^{r}_{theta varphi t}$:

tRiemann_kn_rthvpt

$latex R_{r theta}$:

tRicci_kn_rth

Utilizando nuestras funciones, obtenemos los siguientes gráficos:

$latex M=0.9, J=0.1, Q=0.5$:

R_M09_J01_Q05

$latex M=0.9, J=0.1, Q=0.25$:

R_M09_J01_Q025

$latex M=1, J=1, Q=1$:

R_M1_Q1_J1

$latex M=0.1, J=0.9, Q=0.5$:

R_M01_J09_Q05

Tags: , , , , ,

Volvemos a representar los gráficos de este post pero ahora en coordenadas cilíndricas, que tienen mas sentido:

Tags: , ,

Calcularemos el tensor de Riemann, el de Ricci y la curvatura escalar para la métrica de Kerr correspondiente a un agujero negro en rotación y sin carga eléctrica ($latex J neq 0, Q=0$), cuya métrica ya utilizamos.

A continuación, mostramos nuestras funciones para realizar los calculos automáticamente:

tensor de Riemann:

riemanNd

tensor de Ricci:

ricciNd

curvatura escalar:

rNd

y obtenemos (solo escribiremos una elemento de cada debido a su complejidad):

$latex R^{t}_{ttvarphi}$:

tRiemann_rtttvp

donde $latex x1=t, x2=r, x3=theta, x4=varphi$.

$latex R_{rtheta}$:

tRicci_rth

y dos gráficas correspondientes a su curvatura escalar $latex R$:

curvaturaEscalar3D2

curvaturaEscalar3D

En la definición de la métrica, tenemos la restricción $latex 0 leq frac{a}{M} leq 1$, que en nuestro caso, como imponemos $latex M=1$, nos queda $latex 0 leq J leq 1$. A continuación una serie de gráficos en los que hacemos el valor de

$latex J=0.1, 0.25, 0.5, 0.75, 0.9, 1$:

Tags: , , , , , , ,

Como ya comentamos, existen diferentes soluciones analíticas de las ecuaciones de Einstein correspondientes a los diferentes tipos de BH en equilíbrio. En la tesis “Evolution formalism of Einstein equations: numerical and geometrical issues” de I. Cordero podemos encontrarlas.

Para empezar, la consideración de variedades de Lorentz con simetría esférica y tensor de Ricci nulo, nos lleva a la métrica de Schwarzschild ($latex J=0, Q=0$) que podemos escribir en diferentes sistemas de coordenadas:

Coordenadas de Schwarzschild $latex (r,theta,varphi,tau)$ con $latex r > 2M$ y siendo $latex tau$ el tiempo propio:

$latex g = frac{1}{1-frac{2M}{r}}dr otimes dr + r^2 (dtheta otimes dtheta + sin^2 theta dvarphi otimes dvarphi)-(1-frac{2M}{r})dtau otimes dtau$

que en forma matricial queda:

$latex g=begin{bmatrix} frac{1}{1-frac{2M}{r}} & 0 & 0 & 0 \ 0 & r^2 & 0 & 0 \ 0 & 0 & r^2 sin^2 theta & 0 \ 0 & 0 & 0 & -(1-frac{2M}{r}) end{bmatrix}$

y en física se suele escribir:

$latex ds^2 = frac{1}{1-frac{2M}{r}}dr^2+ r^2 (dtheta^2 + sin^2 theta dvarphi^2)-(1-frac{2M}{r})dtau^2$

Además, como $latex dOmega^2 = dtheta^2 + sin^2theta dvarphi^2$ es la métrica de $latex S^2$ ($latex S^2(theta,varphi) = (sin theta cos varphi, sin theta sin varphi, cos theta)$ en $latex ]0,pi[ times ]0,2pi[$ de manera que $latex g_{11} = S^2_theta cdot S^2_theta = 1$, $latex g_{12} = g_{21} = S^2_theta cdot S^2_varphi = 0$ y $latex g_{22} = S^2_varphi cdot S^2_varphi = sin^2 theta$, con lo que $latex g = dtheta otimes dtheta + sin^2 theta dvarphi otimes dvarphi$), tenemos:

$latex ds^2 = frac{1}{1-frac{2M}{r}}dr^2+ r^2 dOmega^2 – (1-frac{2M}{r})dtau^2$

Coordenadas isotrópicas $latex (bar{r},theta,varphi,tau)$ con $latex r = bar{r} (1 + frac{M}{2bar{r}})^2$ respecto de las de Schwarzschild:

$latex ds^2 = (1+frac{M}{2bar{r}})^4(dbar{r}^2+ bar{r}^2 dOmega^2 )- big (frac{1-frac{M}{2bar{r}}}{1+frac{M}{2bar{r}}} big) dtau^2$

Las métricas inducidas por estas dos métricas en las hipersuperficies de tiempo propio constante son conformemente planas y singulares en el horizonte.

Coordenadas de Painlevé-Gullstrand-Lemaître $latex (r,theta,varphi, T)$ con $latex dT = dtau + frac{sqrt{frac{2M}{r}}}{1-frac{2M}{r}}dr$ respecto de las de Schwarzschild:

$latex ds^2 = dr^2 + r^2 dOmega^2 + 2 sqrt{frac{2M}{r}}dTdr – (1-frac{2M}{r})dT^2$

que en forma matricial queda:

$latex g=begin{bmatrix} 1 & 0 & 0 & sqrt{frac{2M}{r}} \ 0 & r^2 & 0 & 0 \ 0 & 0 & r^2 sin^2 theta & 0 \ sqrt{frac{2M}{r}} & 0 & 0 & -(1-frac{2M}{r}) end{bmatrix}$

Coordenadas de Eddington-Finkelstein $latex (t, r, theta, varphi)$ con $latex t = tau + 2M ln |frac{r}{2M} – 1|$ respecto de las de Schwarzschild:

$latex ds^2 = frac{1}{1+frac{2M}{r}} dr^2 + r^2 dOmega^2 + frac{4M}{r} dtdr – (1-frac{2M}{r})dt^2$

que en forma matricial queda:

$latex g=begin{bmatrix} 1+frac{2M}{r} & 0 & 0 & frac{2M}{r} \ 0 & r^2 & 0 & 0 \ 0 & 0 & r^2 sin^2 theta & 0 \ frac{2M}{r} & 0 & 0 & -(1-frac{2M}{r}) end{bmatrix}$

Las métricas inducidas por estas dos métricas en las hipersuperficies de tiempo constante son planas y regulares en el horizonte.

En el llibre “Geometria diferencial i relativitat” de J. Girbau tambe comenta les coordenadas de Kruskal-Szekeres $latex (u,v,theta,varphi)$:

$latex ds^2 = frac{32M^3}{r} e^{-frac{r}{2M}} (du^2 – dv^2) + r^2 dOmega^2$

donde

$latex u=sqrt{frac{r}{2M}-1} e^{frac{r}{4M}} cosh frac{tau}{4M}$

y

$latex v=sqrt{frac{r}{2M}-1} e^{frac{r}{4M}} sinh frac{tau}{4M}$

No hay singularidad física en $latex r=2M$, pero hay dos en $latex r=0$.

En segundo lugar, tenemos la métrica de Kerr ($latex J neq 0, Q = 0$) que también podemos escribir en diferentes sistemas de coordenadas:

Coordenadas de Boyer-Lindquist $latex (r,theta,varphi,t)$:

$latex ds^2 = frac{rho^2}{Delta} dr^2 + rho^2 dtheta^2 + tilde{w}^2(dvarphi – wdt)^2 – (frac{rho sqrt{Delta}}{Sigma})^2dt^2$

donde

$latex Delta = r^2 -2Mr + a^2$

$latex rho^2 = r^2 + a^2 cos^2 theta$

$latex Sigma^2 = (r^2 + a^2)^2 – a^2 Delta sin^2 theta$

$latex w = frac{2aMr}{Sigma^2}$

$latex tilde{w} = frac{Sigma sin theta}{rho}$

y siendo $latex a$ el momento angular del BH. Fijando $latex a=0$ obtenemos el BH de Schwarzchild en coordenadas de Schwarzchild.

Coordenadas de Kerr-Schild $latex (r,theta,bar{varphi},bar{t})$:

$latex ds^2 = frac{Z^{2k}-1}{Z-1} dr^2 + rho^2 dtheta^2+ sin^2 theta rho^2 [1+Y(1+Z)] dbar{varphi}^2 – (1-Z) dbar{t}^2+$

$latex +2aepsilon sin^2 theta frac{Z^{k+1}-1}{Z-1}drdbar{varphi} -2 epsilon Z^k dr dbar{t} -2 a sin^2 theta Z dbar{varphi}dbar{t}$

donde

$latex Y = frac{a^2 sin^2 theta}{rho^2}$, $latex Z = frac{2Mr}{rho^2}$

y $latex epsilon = +1(-1)$ regulariza el horizonte futuro (pasado) del BH. La relación con las anteriores coordenadas viene dada por

$latex dbar{varphi} = dvarphi – epsilon frac{a}{Delta} dr$

$latex dbar{t} = dt – epsilon [ frac{1+Y}{1+Y-Z} – frac{1-Z^k}{1-Z} ] dr$

donde $latex Delta$ es la función horizonte, que es cero en el horizonte y hace que la componente $latex g_{tt}$ de la métrica se anule.

Tags: , , , , , , , , , , , ,

Las estrellas con una masa mayor a 8 masas solares finalizan su evolución hidroestática con una supernova. El remanente de este evento puede ser tanto un agujero negro, en el caso de las estrellas mas masivas , o una estrella compacta, comunmente llamada estrella de neutrones, para el resto de estrellas.

Existe la posibilidad de que estos objetos sean en realidad estrellas de quarks.

Tags: , , , ,

FireStats icon Powered by FireStats