Astronomy:Non-expanding horizon

From HandWiki
Revision as of 00:58, 24 October 2022 by MainAI (talk | contribs) (link)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

A non-expanding horizon (NEH) is an enclosed null surface whose intrinsic structure is preserved. An NEH is the geometric prototype of an isolated horizon which describes a black hole in equilibrium with its exterior from the quasilocal perspective. It is based on the concept and geometry of NEHs that the two quasilocal definitions of black holes, weakly isolated horizons and isolated horizons, are developed.

Definition of NEHs

A three-dimensional submanifold ∆ is defined as a generic (rotating and distorted) NEH if it respects the following conditions:[1][2][3]


(i) ∆ is null and topologically [math]\displaystyle{ S^2\times\mathbb{R} }[/math];
(ii) Along any null normal field [math]\displaystyle{ l }[/math] tangent to ∆, the outgoing expansion rate [math]\displaystyle{ \displaystyle \theta_{(l)}:=\hat{h}^{ab}\hat{\nabla}_a l_b }[/math] vanishes;
(iii) All field equations hold on ∆, and the stress–energy tensor [math]\displaystyle{ T_{ab} }[/math] on ∆ is such that [math]\displaystyle{ V^a:=-T^a_b l^b }[/math] is a future-directed causal vector ([math]\displaystyle{ V^a V_a\leq 0 }[/math]) for any future-directed null normal [math]\displaystyle{ l^a }[/math].


Condition (i) is fairly trivial and just states the general fact that from a 3+1 perspective[4] an NEH ∆ is foliated by spacelike 2-spheres ∆'=S2, where S2 emphasizes that ∆' is topologically compact with genus zero ([math]\displaystyle{ g=0 }[/math]). The signature of ∆ is (0,+,+) with a degenerate temporal coordinate, and the intrinsic geometry of a foliation leaf ∆'=S2 is nonevolutional. The property [math]\displaystyle{ \theta_{(l)}=0 }[/math] in condition (ii) plays a pivotal role in defining NEHs and the rich implications encoded therein will be extensively discussed below. Condition (iii) makes one feel free to apply the Newman–Penrose (NP) formalism[5][6] of Einstein-Maxwell field equations to the horizon and its near-horizon vicinity; furthermore, the very energy inequality is motivated from the dominant energy condition[7] and is a sufficient condition for deriving many boundary conditions of NEHs.


Note: In this article, following the convention set up in refs.,[1][2][3] "hat" over the equality symbol [math]\displaystyle{ \hat= }[/math] means equality on the black-hole horizons (NEHs), and "hat" over quantities and operators ([math]\displaystyle{ \hat h^{ab} }[/math], [math]\displaystyle{ \hat\nabla }[/math], etc.) denotes those on a foliation leaf of the horizon. Also, ∆ is the standard symbol for both an NEH and the directional derivative ∆[math]\displaystyle{ :=n^a\nabla_a }[/math] in NP formalism, and we believe this won't cause an ambiguity.

Boundary conditions implied by the definition

Now let's work out the implications of the definition of NEHs, and these results will be expressed in the language of NP formalism with the convention[5][6] [math]\displaystyle{ \{(-,+,+,+); l^an_a=-1, m^a\bar{m}_a=1\} }[/math] (Note: unlike the original convention[8][9] [math]\displaystyle{ \{(+,-,-,-); l^an_a=1, m^a\bar{m}_a=-1\} }[/math], this is the usual one employed in studying trapped null surfaces and quasilocal definitions of black holes[10]). Being a null normal to ∆, [math]\displaystyle{ l^a }[/math] is automatically geodesic, [math]\displaystyle{ \kappa:= -m^al^b\nabla_b l_a\,\hat{=}\,0 }[/math], and twist free, [math]\displaystyle{ \text{Im}(\rho)=\text{Im}(-m^a \bar{m}^b \nabla_b l_a )\,\hat{=}\,0 }[/math]. For an NEH, the outgoing expansion rate [math]\displaystyle{ \theta_{(l)} }[/math] along [math]\displaystyle{ l^a }[/math] is vanishing, [math]\displaystyle{ \theta_{(l)}\,\hat{=}\,0 }[/math], and consequently [math]\displaystyle{ \text{Re}(\rho)=\text{Re}(-m^a \bar{m}^b \nabla_b l_a )=- \frac{1}{2}\theta_{(l)}\,\hat{=}\,0 }[/math]. Moreover, according to the Raychaudhuri-NP expansion-twist equation,[11]

[math]\displaystyle{ (1)\qquad D\rho=\rho^2+\sigma\bar{\sigma}+\frac{1}{2}R_{ab}l^a l^b\,\hat{=}\,0\,, }[/math]

it follows that on ∆

[math]\displaystyle{ (2)\qquad \sigma\bar{\sigma}+\frac{1}{2}R_{ab}l^a l^b\,\hat{=}\,0\,, }[/math]

where [math]\displaystyle{ \sigma:= -m^b m^a \nabla_a l_b }[/math] is the NP-shear coefficient. Due to the assumed energy condition (iii), we have [math]\displaystyle{ R_{ab}l^a l^b=R_{ab}l^a l^b-\frac{1}{2}Rg_{ab}l^al^b=8\pi \, T_{ab}l^a l^b }[/math] ([math]\displaystyle{ c=G=1 }[/math]), and therefore [math]\displaystyle{ R_{ab}l^a l^b }[/math] is nonnegative on ∆. The product [math]\displaystyle{ \sigma\bar{\sigma} }[/math] is of course nonnegative, too. Consequently, [math]\displaystyle{ \sigma\bar{\sigma} }[/math] and [math]\displaystyle{ R_{ab}l^a l^b }[/math] must be simultaneously zero on ∆, i.e. [math]\displaystyle{ \sigma\,\hat{=}\,0 }[/math] and [math]\displaystyle{ R_{ab}l^a l^b\,\hat{=}\,0 }[/math]. As a summary,

[math]\displaystyle{ (3)\qquad \kappa\,\hat{=}\,0\,,\quad\text{Im}(\rho)\,\hat{=}\,0\,,\quad \text{Re}(\rho)\,\hat{=}\,0\,,\quad\sigma\,\hat{=}\,0\,,\quad R_{ab}l^a l^b\,\hat{=}\,0. }[/math]

Thus, the isolated horizon ∆ is nonevolutional and all foliation leaves ∆'=S2 look identical with one another. The relation [math]\displaystyle{ R_{ab}l^a l^b=8\pi \cdot T_{ab}l^a l^b=8\pi \cdot T^a_{b} l^b\cdot l_a\,\hat{=}\,0 }[/math] implies that the causal vector [math]\displaystyle{ -T^a_b l^b }[/math] in condition (iii) is proportional to [math]\displaystyle{ l^a }[/math] and [math]\displaystyle{ R_{ab}l^b }[/math] is proportional to [math]\displaystyle{ l_a }[/math] on the horizon ∆; that is, [math]\displaystyle{ -T^a_b l^b\,\hat{=}\,cl^a }[/math] and [math]\displaystyle{ R_{ab}l^b\,\hat{=}\,cl_a }[/math], [math]\displaystyle{ c\in\mathbb{R} }[/math]. Applying this result to the related Ricci-NP scalars, we get [math]\displaystyle{ \Phi_{00}:= \frac{1}{2}R_{ab}l^a l^b\,\hat{=}\,\frac{c}{2}\,l_b l^b\,\hat{=}\,0 }[/math], and [math]\displaystyle{ \Phi_{01}=\overline{\Phi_{10}}:= \frac{1}{2}R_{ab}l^a m^b\,\hat{=}\,\frac{c}{2}\,l_b m^b\,\hat{=}\,0 }[/math], thus

[math]\displaystyle{ (4)\qquad R_{ab}l^b\,\hat{=}\,cl_a\,,\quad \Phi_{00}\,\hat{=}\,0\,,\quad \Phi_{10}=\overline{\Phi_{01}}\,\hat{=}\,0\,. }[/math]

The vanishing of Ricci-NP scalars [math]\displaystyle{ \{\Phi_{00}\,,\Phi_{01}\,, \Phi_{10}\} }[/math] signifies that, there is no energy–momentum flux of any kind of charge across the horizon, such as electromagnetic waves, Yang–Mills flux or dilaton flux. Also, there should be no gravitational waves crossing the horizon; however, gravitational waves are propagation of perturbations of the spacetime continuum rather than flows of charges, and therefore depicted by four Weyl-NP scalars [math]\displaystyle{ \Psi_i\; (i=0,1,3,4) }[/math] (excluding [math]\displaystyle{ \Psi_2 }[/math]) rather than Ricci-NP quantities [math]\displaystyle{ \Phi_{ij} }[/math].[5] According to the Raychaudhuri-NP shear equation[11]

[math]\displaystyle{ (5)\qquad D\sigma=\sigma(\rho+\bar\rho)+\Psi_0=-2\sigma\theta_{(l)}+\Psi_0\,, }[/math]

or the NP field equation on the horizon

[math]\displaystyle{ (6)\qquad D\sigma-\delta\kappa=(\rho+\bar{\rho})\sigma+(3\varepsilon-\bar{\varepsilon})\sigma-(\tau-\bar{\pi}+\bar{\alpha}+3\beta)\kappa+\Psi_0\,\hat{=}\,0\,, }[/math]

it follows that [math]\displaystyle{ \Psi_0:= C_{abcd} l^a m^b l^c m^d\,\hat{=}\,0 }[/math]. Moreover, the NP equation

[math]\displaystyle{ (7)\qquad \delta\rho-\bar{\delta}\sigma=\rho(\bar{\alpha}+\beta)-\sigma(3\alpha-\bar{\beta})+(\rho-\bar{\rho})\tau+(\mu-\bar{\mu})\kappa-\Psi_1+\Phi_{01}\,\hat{=}\,0 }[/math]

implies that [math]\displaystyle{ \Psi_1:= C_{abcd} l^a n^b l^c m^d\,\hat{=}\,0 }[/math]. To sum up, we have

[math]\displaystyle{ (8)\qquad \Psi_0\,\hat{=}\,0\,,\quad \Psi_1\,\hat{=}\,0\,, }[/math]

which means that,[5] geometrically, a principal null direction of Weyl's tensor is repeated twice and [math]\displaystyle{ l^a }[/math] is aligned with the principal direction; physically, no gravitational waves (transverse component [math]\displaystyle{ \Psi_0 }[/math] and longitudinal component [math]\displaystyle{ \Psi_1 }[/math]) enter the black hole. This result is consistent with the physical scenario defining NEHs.

Remarks: Spin coefficients related to Raychaudhuri's equation

For a better understanding of the previous section, we will briefly review the meanings of relevant NP spin coefficients in depicting null congruences.[7] The tensor form of Raychaudhuri's equation[12] governing null flows reads

[math]\displaystyle{ (9)\qquad \mathcal{L}_{\ell}\theta_{(l)}=-\frac{1}{2}\theta_{(l)}^2+\tilde{\kappa}_{(l)}\theta_{(l)}-\sigma_{ab}\sigma^{ab}+\tilde{\omega}_{ab}\tilde{\omega}^{ab}-R_{ab}l^a l^b\,, }[/math]

where [math]\displaystyle{ \tilde{\kappa}_{(l)} }[/math] is defined such that [math]\displaystyle{ \tilde{\kappa}_{(l)}l^b:= l^a \nabla_a l^b }[/math]. The quantities in Raychaudhuri's equation are related with the spin coefficients via[5][13][14]

[math]\displaystyle{ (10)\qquad \theta_{(l)}=-(\rho+\bar\rho)=-2\text{Re}(\rho)\,,\quad \theta_{(n)}=\mu+\bar\mu=2\text{Re}(\mu)\,, }[/math]
[math]\displaystyle{ (11)\qquad \sigma_{ab}=-\sigma \bar m_a \bar m_b-\bar\sigma m_a m_b\,, }[/math]
[math]\displaystyle{ (12)\qquad \tilde{\omega}_{ab}=\frac{1}{2}\,\Big(\rho-\bar\rho \Big)\,\Big(m_a \bar m_b-\bar m_a m_b \Big)=\text{Im}(\rho)\cdot\Big(m_a \bar m_b-\bar m_a m_b \Big)\,, }[/math]
where Eq(10) follows directly from [math]\displaystyle{ \hat{h}^{ab}=\hat{h}^{ba}=m^b\bar m^a+\bar m^b m^a }[/math] and
[math]\displaystyle{ (13)\qquad \theta_{(l)}=\hat{h}^{ba}\nabla_a l_b=m^b\bar m^a\nabla_a l_b+\bar m^b m^a\nabla_a l_b =m^b\bar \delta l_b+\bar m^b \delta l_b=-(\rho+\bar\rho)\,, }[/math]
[math]\displaystyle{ (14)\qquad \theta_{(n)}=\hat{h}^{ba}\nabla_a n_b=\bar m^b m^a\nabla_a n_b+m^b\bar m^a\nabla_a n_b=\bar m^b \delta n_b+m^b\bar \delta n_b=\mu+\bar\mu\,. }[/math]

Moreover, a null congruence is hypersurface orthogonal if [math]\displaystyle{ \text{Im}(\rho)=0 }[/math].[5]

Constraints from electromagnetic fields

Vacuum NEHs on which [math]\displaystyle{ \{\Phi_{ij}\hat{=}0\,,\Lambda_{}\hat{=}0\} }[/math] are the simplest types of NEHs, but in general there can be various physically meaningful fields surrounding an NEH, among which we are mostly interested in electrovacuum fields with [math]\displaystyle{ \Lambda\hat{=}0 }[/math]. This is the simplest extension of vacuum NEHs, and the nonvanishing energy-stress tensor for electromagnetic fields reads

[math]\displaystyle{ (15)\qquad T_{ab}=\frac{1}{4\pi}\,\Big(\, F_{ac}F_b^c -\frac{1}{4}g_{ab}F_{cd}F^{cd} \Big)\,, }[/math]


where [math]\displaystyle{ F_{ab} }[/math] refers to the antisymmetric ([math]\displaystyle{ F_{ab}=-F_{ba} }[/math], [math]\displaystyle{ F^a_a=0 }[/math]) electromagnetic field strength, and [math]\displaystyle{ T_{ab} }[/math] is trace-free ([math]\displaystyle{ T^a_a=0 }[/math]) by definition and respects the dominant energy condition. (One should be careful with the antisymmetry of [math]\displaystyle{ F_{ab} }[/math] in defining Maxwell-NP scalars [math]\displaystyle{ \phi_{i} }[/math]).

The boundary conditions derived in the previous section are applicable to generic NEHs. In the electromagnetic case, [math]\displaystyle{ \Phi_{ij} }[/math] can be specified in a more particular way. By the NP formalism of Einstein-Maxwell equations, one has[5]

[math]\displaystyle{ (16)\qquad \Phi_{ij}=\, 2\, \phi_i\, \overline{\phi_j}\,,\quad i,j\in\{0,1,2\}\,, }[/math]

where [math]\displaystyle{ \phi_i }[/math] denote the three Maxwell-NP scalars. As an alternative to Eq(), we can see that the condition [math]\displaystyle{ \Phi_{00}=0 }[/math] also results from the NP equation

[math]\displaystyle{ (17)\qquad D\rho -\bar{\delta}\kappa=(\rho^2+\sigma\bar{\sigma})+(\varepsilon+\bar{\varepsilon})\rho-\bar{\kappa}\tau-(3\alpha+\bar{\beta}-\pi)\,\kappa+\Phi_{00}\,\hat{=}\,0\, }[/math]
as [math]\displaystyle{ \kappa_{}\,\hat{=}\,\rho\,\hat{=}\,\sigma=0 }[/math], so
[math]\displaystyle{ (18)\qquad \Phi_{00}\,\hat{=}\,0\;\;\Leftrightarrow\;\;2\,\phi_{0}\,\overline{\phi_0}\,\hat{=}\,0\;\;\Rightarrow \;\;\phi_0=\overline{\phi_0}\,\hat{=}\,0\,. }[/math]

It follows straightforwardly that

[math]\displaystyle{ (19)\qquad \Phi_{01}=\overline{\Phi_{10}}=\,2\,\phi_0\,\overline{\phi_1}\,\hat{=}\,0\,,\quad \Phi_{02}=\overline{\Phi_{20}}=\,2\,\phi_0\,\overline{\phi_2}\,\hat{=}\,0\,. }[/math]


These results demonstrate that, there are no electromagnetic waves across ([math]\displaystyle{ \Phi_{00} }[/math], [math]\displaystyle{ \Phi_{01} }[/math]) or along (\Phi_{02}) the NEH except the null geodesics generating the horizon. It is also worthwhile to point out that, the supplementary equation [math]\displaystyle{ \Phi_{ij}=2\,\phi_i\, \overline{\phi_j} }[/math] in Eq() is only valid for electromagnetic fields; for example, in the case of Yang–Mills fields there will be [math]\displaystyle{ \Phi_{ij}=\,\text{Tr}\,\big(\,\digamma_i \,\bar{\digamma}_j \,\big) }[/math] where [math]\displaystyle{ \digamma_i(i\in\{0,1,2 \} }[/math] are Yang–Mills-NP scalars.[15]

Adapted tetrad on NEHs and further properties

Usually, null tetrads adapted to spacetime properties are employed to achieve the most succinct NP descriptions. For example, a null tetrad can be adapted to principal null directions once the Petrov type is known; also, at some typical boundary regions such as null infinity, timelike infinity, spacelike infinity, black hole horizons and cosmological horizons, tetrads can be adapted to boundary structures. Similarly, a preferred tetrad[1][2][3] adapted to on-horizon geometric behaviors is employed in the literature to further investigate NEHs.

As indicated from the 3+1 perspective from condition (i) in the definition, an NEH ∆ is foliated by spacelike hypersurfaces ∆'=S2 transverse to its null normal along an ingoing null coordinate [math]\displaystyle{ v }[/math], where we follow the standard notation of ingoing Eddington–Finkelstein null coordinates and use [math]\displaystyle{ v }[/math] to label the 2-dimensional leaves [math]\displaystyle{ S^2_v }[/math] at [math]\displaystyle{ v=\text{constant} }[/math]; that is, [math]\displaystyle{ {}\Delta={}\Delta'x[v_o,v_1]=S^2[v_o,v_1] }[/math]. [math]\displaystyle{ v }[/math] is set to be future-directed and choose the first tetrad covector [math]\displaystyle{ n_a }[/math] as [math]\displaystyle{ n_a=-dv }[/math],[2][3] and then there will be a unique vector field [math]\displaystyle{ l^a }[/math] as null normals to [math]\displaystyle{ S^2_v }[/math] satisfying the cross-normalization [math]\displaystyle{ l^an_a=-1 }[/math] and affine parametrization [math]\displaystyle{ Dv=1 }[/math]; such choice of [math]\displaystyle{ \{l^a,n^a\} }[/math] would actually yields a preferred foliation of ∆. While [math]\displaystyle{ \{l^a\,,n^a\} }[/math] are related to the extrinsic properties and null generators (i.e. null flows/geodesic congruence on ∆), the remaining two complex null vectors [math]\displaystyle{ \{m^a,\bar{m}^a\} }[/math] are to span the intrinsic geometry of a foliation leaf [math]\displaystyle{ S^2_v }[/math], tangent to ∆ and transverse to [math]\displaystyle{ \{l^a\,,n^a\} }[/math]; that is, [math]\displaystyle{ \mathcal{L}_{\ell}m\,\hat{=}\,\mathcal{L}_{\ell}\bar{m}\hat{=}0 }[/math].

Now let's check the consequences of this kind of adapted tetrad. Since

[math]\displaystyle{ (20)\qquad \mathcal{L}_{\ell}m=[\ell,m]\,\hat{=}\,0\;\Rightarrow\; \delta D-D\delta=(\bar{\alpha}+\beta-\bar{\pi})D+\kappa_{}\Delta-(\bar{\rho}+\varepsilon-\bar{\varepsilon})\delta-\sigma\bar{\delta} \,\hat{=}\,0 \,, }[/math]

with [math]\displaystyle{ \kappa_{}\,\hat{=}\,\rho\,\hat{=}\,\sigma\,\hat{=}\,0 }[/math], we have

[math]\displaystyle{ (21)\qquad \pi\,\hat{=}\,\alpha+\bar{\beta} \,,\quad \varepsilon\,\hat{=}\,\bar{\varepsilon} \,. }[/math]

Also, in such an adapted frame, the derivative [math]\displaystyle{ \mathcal{L}_{\bar{m}}m }[/math] on [math]\displaystyle{ {}\Delta={}\Delta'x[v_o,v_1]=S^2[v_o,v_1] }[/math] should be purely intrinsic; thus in the commutator

[math]\displaystyle{ (22)\qquad \mathcal{L}_{\bar{m}}m=[\bar{m}, m]=\bar{\delta}\delta-\delta\bar{\delta}= (\bar{\mu}-\mu)D+(\bar{\rho}-\rho)\Delta-(\bar{\beta}-\alpha)\delta-(\bar{\alpha}-\beta)\bar{\delta}\,, }[/math]

the coefficients for the directional derivatives [math]\displaystyle{ D }[/math] and ∆ must be zero, that is

[math]\displaystyle{ (23)\qquad \bar{\mu}\,\hat{=}\,\mu \,,\quad \mathcal{L}_{\bar{m}}m \,\hat{=}\,(\alpha-\bar{\beta})\delta-(\bar{\alpha}-\beta)\bar{\delta}\,, }[/math]

so the ingoing null normal field [math]\displaystyle{ n^a }[/math] is twist-free by [math]\displaystyle{ \text{Im}(\mu)=\text{Im}(\bar{m}^a m^b\nabla_b n_a)=0 }[/math], and [math]\displaystyle{ 2\mu=2\text{Re}(\mu) }[/math] equals the ingoing expansion rate [math]\displaystyle{ \theta_{(n)} }[/math].

Discussion

So far, the definition and boundary conditions of NEHs have been introduced. The boundary conditions include those for an arbitrary NEH, specific characteristics for Einstein-Maxwell (electromagnetic) NEHs, as well as further properties in an adapted tetrad. Based on NEHs, WIHs which have valid surface gravity can be defined to generalize the black hole mechanics. WIHs are sufficient in studying the physics on the horizon, but for geometric purposes,[2] stronger restrictions can be imposed to WIHs so as to introduce IHs, where the equivalence class of null normals [math]\displaystyle{ [\ell] }[/math] fully preserves the induced connection [math]\displaystyle{ \mathcal {D} }[/math] on the horizon.

References

  1. 1.0 1.1 1.2 Abhay Ashtekar, Christopher Beetle, Olaf Dreyer, et al. "Generic isolated horizons and their applications". Physical Review Letters, 2000, 85(17): 3564-3567. arXiv:gr-qc/0006006v2
  2. 2.0 2.1 2.2 2.3 2.4 Abhay Ashtekar, Christopher Beetle, Jerzy Lewandowski. "Geometry of generic isolated horizons". Classical and Quantum Gravity, 2002, 19(6): 1195-1225. arXiv:gr-qc/0111067v2
  3. 3.0 3.1 3.2 3.3 Abhay Ashtekar, Stephen Fairhurst, Badri Krishnan. "Isolated horizons: Hamiltonian evolution and the first law". Physical Review D, 2000, 62(10): 104025. gr-qc/0005083
  4. Thomas W Baumgarte, Stuart L Shapiro. Numerical Relativity: Solving Einstein's Equations on the Computer. Cambridge: Cambridge University Press, 2010. Chapter 2: The 3+1 decomposition of Einstein's equations, page 23.
  5. 5.0 5.1 5.2 5.3 5.4 5.5 5.6 Jeremy Bransom Griffiths, Jiri Podolsky. Exact Space-Times in Einstein's General Relativity. Cambridge: Cambridge University Press, 2009. Chapter 2.
  6. 6.0 6.1 Valeri P Frolov, Igor D Novikov. Black Hole Physics: Basic Concepts and New Developments. Berlin: Springer, 1998. Appendix E.
  7. 7.0 7.1 Eric Poisson. A Relativist's Toolkit: The Mathematics of Black-Hole Mechanics. Cambridge: Cambridge University Press, 2004. Chapters 2 and 3.
  8. Ezra T Newman, Roger Penrose. "An Approach to Gravitational Radiation by a Method of Spin Coefficients". Journal of Mathematical Physics, 1962, 3(3): 566-768.
  9. Ezra T Newman, Roger Penrose. "Errata: An Approach to Gravitational Radiation by a Method of Spin Coefficients". Journal of Mathematical Physics, 1963, 4(7): 998.
  10. Ivan Booth. "Black hole boundaries". Canadian Journal of Physics, 2005, 83(11): 1073-1099. [arxiv.org/abs/gr-qc/0508107 arXiv:gr-qc/0508107v2]
  11. 11.0 11.1 Subrahmanyan Chandrasekhar. The Mathematical Theory of Black Holes. Chicago: University of Chicago Press, 1983. Section 9(a), page 56.
  12. Sayan Kar, Soumitra SenGupta. The Raychaudhuri equations: a brief review. Pramana, 2007, 69(1): 49-76. [arxiv.org/abs/gr-qc/0611123v1 gr-qc/0611123]
  13. David McMahon. Relativity Demystified - A Self-Teaching Guide. New York: McGraw-Hill, 2006. Chapter 9.
  14. Alex Nielsen. PhD thesis: Black Hole Horizons and Black Hole Thermodynamics. University of Canterbury, 2007. Section 2.3. Available online: http://ir.canterbury.ac.nz/handle/10092/1363.
  15. E T Newman, K P Tod. Asymptotically Flat Spacetimes. page 27, Appendix A.2. In A Held (Editor): General relativity and gravitation: one hundred years after the birth of Albert Einstein. Vol(2). New York and London: Plenum Press, 1980.