Physics:Optical scalars

From HandWiki

In general relativity, optical scalars refer to a set of three scalar functions [math]\displaystyle{ \{ \hat\theta }[/math] (expansion), [math]\displaystyle{ \hat\sigma }[/math] (shear) and [math]\displaystyle{ \hat\omega }[/math] (twist/rotation/vorticity)[math]\displaystyle{ \} }[/math] describing the propagation of a geodesic null congruence.[1][2][3][4][5]
In fact, these three scalars [math]\displaystyle{ \{ \hat\theta\,,\hat\sigma\,, \hat\omega \} }[/math] can be defined for both timelike and null geodesic congruences in an identical spirit, but they are called "optical scalars" only for the null case. Also, it is their tensorial predecessors [math]\displaystyle{ \{ \hat\theta \hat h_{ab}\,,\hat\sigma_{ab}\,, \hat\omega_{ab} \} }[/math] that are adopted in tensorial equations, while the scalars [math]\displaystyle{ \{ \hat\theta\,,\hat\sigma\,, \hat\omega \} }[/math] mainly show up in equations written in the language of Newman–Penrose formalism.

Definitions: expansion, shear and twist

For geodesic timelike congruences

Denote the tangent vector field of an observer's worldline (in a timelike congruence) as [math]\displaystyle{ Z^a }[/math], and then one could construct induced "spatial metrics" that


[math]\displaystyle{ (1)\quad h^{ab}=g^{ab}+Z^a Z^b\;,\quad h_{ab}=g_{ab}+Z_a Z_b\;,\quad h^a_{\;\;b}=\delta^a_{\;\;b}+Z^a Z_b\;, }[/math]


where [math]\displaystyle{ h^a_{\;\;b} }[/math] works as a spatially projecting operator. Use [math]\displaystyle{ h^a_{\;\;b} }[/math] to project the coordinate covariant derivative [math]\displaystyle{ \nabla_b Z_a }[/math] and one obtains the "spatial" auxiliary tensor [math]\displaystyle{ B_{ab} }[/math],


[math]\displaystyle{ (2)\quad B_{ab}=h^c_{\;\;a}\, h^d_{\;\;b}\, \nabla_d Z_c = \nabla_b Z_a +A_a Z_b\;, }[/math]


where [math]\displaystyle{ A_a }[/math] represents the four-acceleration, and [math]\displaystyle{ B_{ab} }[/math] is purely spatial in the sense that [math]\displaystyle{ B_{ab}Z^a=B_{ab}Z^b=0 }[/math]. Specifically for an observer with a geodesic timelike worldline, we have


[math]\displaystyle{ (3)\quad A_a=0\;,\quad\Rightarrow\quad B_{ab}= \nabla_b Z_a\;. }[/math]


Now decompose [math]\displaystyle{ B_{ab} }[/math] into its symmetric and antisymmetric parts [math]\displaystyle{ \theta_{ab} }[/math] and [math]\displaystyle{ \omega_{ab} }[/math],


[math]\displaystyle{ (4)\quad \theta_{ab}=B_{(ab)}\;,\quad \omega_{ab}=B_{[ab]}\;. }[/math]


[math]\displaystyle{ \omega_{ab}=B_{[ab]} }[/math] is trace-free ([math]\displaystyle{ g^{ab}\omega_{ab}=0 }[/math]) while [math]\displaystyle{ \theta_{ab} }[/math] has nonzero trace, [math]\displaystyle{ g^{ab}\theta_{ab}=\theta }[/math]. Thus, the symmetric part [math]\displaystyle{ \theta_{ab} }[/math] can be further rewritten into its trace and trace-free part,


[math]\displaystyle{ (5)\quad \theta_{ab}=\frac{1}{3}\theta h_{ab} +\sigma_{ab}\;. }[/math]


Hence, all in all we have


[math]\displaystyle{ (6)\quad B_{ab}=\frac{1}{3}\theta h_{ab} +\sigma_{ab}+\omega_{ab}\;,\quad \theta=g^{ab}\theta_{ab}=g^{ab}B_{(ab)}\;,\quad \sigma_{ab}=\theta_{ab}-\frac{1}{3}\theta h_{ab}\;,\quad \omega_{ab}=B_{[ab]}\;. }[/math]

For geodesic null congruences

Now, consider a geodesic null congruence with tangent vector field [math]\displaystyle{ k^a }[/math]. Similar to the timelike situation, we also define


[math]\displaystyle{ (7)\quad \hat{B}_{ab}:= \nabla_b k_a\;, }[/math]


which can be decomposed into


[math]\displaystyle{ (8)\quad \hat B_{ab}=\hat\theta_{ab}+\hat\omega_{ab}=\frac{1}{2}\hat\theta \hat h_{ab}+\hat\sigma_{ab}+\hat\omega_{ab}\;, }[/math]


where


[math]\displaystyle{ (9)\quad \hat\theta_{ab}=\hat B_{(ab)}\;,\quad \hat\theta=\hat h^{ab} \hat B_{ab}\;,\quad \hat\sigma_{ab}=\hat B_{(ab)}-\frac{1}{2}\hat\theta \hat h_{ab}\;,\quad \hat\omega_{ab}=\hat B_{[ab]}\;. }[/math]


Here, "hatted" quantities are utilized to stress that these quantities for null congruences are two-dimensional as opposed to the three-dimensional timelike case. However, if we only discuss null congruences in a paper, the hats can be omitted for simplicity.

Definitions: optical scalars for null congruences

The optical scalars [math]\displaystyle{ \{ \hat\theta\,,\hat\sigma\,, \hat\omega \} }[/math][1][2][3][4][5] come straightforwardly from "scalarization" of the tensors [math]\displaystyle{ \{ \hat\theta\,,\hat\sigma_{ab}\,, \hat\omega_{ab} \} }[/math] in Eq(9).


The expansion of a geodesic null congruence is defined by (where for clearance we will adopt another standard symbol "[math]\displaystyle{ ; }[/math]" to denote the covariant derivative [math]\displaystyle{ \nabla_a }[/math])


[math]\displaystyle{ (10)\quad \hat\theta = \frac{1}{2}\, k^a{}_{;\,a} \;. }[/math]

Comparison with the "expansion rates of a null congruence": As shown in the article "Expansion rate of a null congruence", the outgoing and ingoing expansion rates, denoted by [math]\displaystyle{ \theta_{(\ell)} }[/math] and [math]\displaystyle{ \theta_{(n)} }[/math] respectively, are defined by


[math]\displaystyle{ (A.1)\quad \theta_{(\ell)}:=h^{ab}\nabla_a l_b\;, }[/math]


[math]\displaystyle{ (A.2)\quad \theta_{(n)}:=h^{ab}\nabla_a n_b\;, }[/math]


where [math]\displaystyle{ h^{ab}=g^{ab}+l^a n^b+n^a l^b }[/math] represents the induced metric. Also, [math]\displaystyle{ \theta_{(\ell)} }[/math] and [math]\displaystyle{ \theta_{(n)} }[/math] can be calculated via


[math]\displaystyle{ (A.3)\quad \theta_{(\ell)}=g^{ab}\nabla_a l_b -\kappa_{(\ell)}\;, }[/math]


[math]\displaystyle{ (A.4)\quad \theta_{(n)}=g^{ab}\nabla_a n_b -\kappa_{(n)}\;, }[/math]


where [math]\displaystyle{ \kappa_{(\ell)} }[/math] and [math]\displaystyle{ \kappa_{(n)} }[/math] are respectively the outgoing and ingoing non-affinity coefficients defined by


[math]\displaystyle{ (A.5)\quad l^a\nabla_a l_b=\kappa_{(\ell)}l_b\;, }[/math]


[math]\displaystyle{ (A.6)\quad n^a\nabla_a n_b=\kappa_{(n)}n_b\;. }[/math]


Moreover, in the language of Newman–Penrose formalism with the convention [math]\displaystyle{ \{(-,+,+,+); l^a n_a=-1\,,m^a \bar{m}_a=1\} }[/math], we have


[math]\displaystyle{ (A.7)\quad \theta_{(l)}=-(\rho+\bar\rho)=-2\text{Re}(\rho)\,,\quad \theta_{(n)}=\mu+\bar\mu=2\text{Re}(\mu)\,, }[/math]


As we can see, for a geodesic null congruence, the optical scalar [math]\displaystyle{ \theta }[/math] plays the same role with the expansion rates [math]\displaystyle{ \theta_{(\ell)} }[/math] and [math]\displaystyle{ \theta_{(n)} }[/math]. Hence, for a geodesic null congruence, [math]\displaystyle{ \theta }[/math] will be equal to either [math]\displaystyle{ \theta_{(\ell)} }[/math] or [math]\displaystyle{ \theta_{(n)} }[/math].


The shear of a geodesic null congruence is defined by


[math]\displaystyle{ (11)\quad {\hat\sigma} ^2=\hat\sigma_{ab}\hat{\bar\sigma}^{ab} =\frac{1}{2}\,g^{ca}\,g^{db}\,k_{(a\,;\,b)}\,k_{c\,;\,d} - \Big(\frac{1}{2}\, k^a{}_{;\,a} \Big)^2 = \,g^{ca}\,g^{db}\frac{1}{2}\,k_{(a\,;\,b)}\,k_{c\,;\,d} - {\hat\theta}^2\;. }[/math]


The twist of a geodesic null congruence is defined by


[math]\displaystyle{ (12)\quad {\hat\omega}^2 =\frac{1}{2}\,k_{[a\,;\,b]}\,k^{a\,;\,b} =g^{ca}\,g^{db}\,k_{[a\,;\,b]}\,k_{c\,;\,d}\;. }[/math]


In practice, a geodesic null congruence is usually defined by either its outgoing ([math]\displaystyle{ k^a=l^a }[/math]) or ingoing ([math]\displaystyle{ k^a=n^a }[/math]) tangent vector field (which are also its null normals). Thus, we obtain two sets of optical scalars [math]\displaystyle{ \{ \hat\theta_{(\ell)}\,,\hat\sigma_{(\ell)}\,, \hat\omega_{(\ell)} \} }[/math] and [math]\displaystyle{ \{ \hat\theta_{(n)}\,,\hat\sigma_{(n)}\,, \hat\omega_{(n)} \} }[/math], which are defined with respect to [math]\displaystyle{ l^a }[/math] and [math]\displaystyle{ n^a }[/math], respectively.

Applications in decomposing the propagation equations

For a geodesic timelike congruence

The propagation (or evolution) of [math]\displaystyle{ B_{ab} }[/math] for a geodesic timelike congruence along [math]\displaystyle{ Z^c }[/math] respects the following equation,


[math]\displaystyle{ (13)\quad Z^c\nabla_c B_{ab}=-B^c_{\;\;b}B_{ac}+R_{cbad} Z^c Z^d\;. }[/math]


Take the trace of Eq(13) by contracting it with [math]\displaystyle{ g^{ab} }[/math], and Eq(13) becomes


[math]\displaystyle{ (14)\quad Z^c\nabla_c \theta=\theta_{,\,\tau}=-\frac{1}{3}\theta^2 -\sigma_{ab}\sigma^{ab}+\omega_{ab}\omega^{ab}-R_{ab}Z^a Z^b }[/math]


in terms of the quantities in Eq(6). Moreover, the trace-free, symmetric part of Eq(13) is


[math]\displaystyle{ (15)\quad Z^c\nabla_c \sigma_{ab}=-\frac{2}{3}\theta\sigma_{ab}-\sigma_{ac}\sigma^c_{\;b}-\omega_{ac}\omega^c_{\;b}+\frac{1}{3}h_{ab}\,(\sigma_{cd}\sigma^{cd}-\omega_{cd}\omega^{cd})+C_{cbad}Z^c Z^d+\frac{1}{2}\tilde{R}_{ab}\,. }[/math]


Finally, the antisymmetric component of Eq(13) yields


[math]\displaystyle{ (16)\quad Z^c\nabla_c \omega_{ab}=-\frac{2}{3}\theta\omega_{ab}-2\sigma^c_{\;[b}\omega_{a]c}\;. }[/math]

For a geodesic null congruence

A (generic) geodesic null congruence obeys the following propagation equation,


[math]\displaystyle{ (16)\quad k^c\nabla_c \hat B_{ab}=-\hat B^c_{\;\;b}\hat B_{ac}+\widehat{R_{cbad} k^c k^d}\;. }[/math]


With the definitions summarized in Eq(9), Eq(14) could be rewritten into the following componential equations,


[math]\displaystyle{ (17)\quad k^c\nabla_c \hat\theta=\hat\theta_{,\,\lambda}=-\frac{1}{2}\hat\theta^2 -\hat\sigma_{ab}\hat\sigma^{ab}+\hat\omega_{ab}\hat\omega^{ab}-\widehat{R_{cd}k^c k^d}\;, }[/math]


[math]\displaystyle{ (18)\quad k^c\nabla_c \hat\sigma_{ab}=-\hat\theta\hat\sigma_{ab}+\widehat{C_{cbad}k^c k^d}\;, }[/math]


[math]\displaystyle{ (19)\quad k^c\nabla_c \hat\omega_{ab}=-\hat\theta\hat\omega_{ab}\;. }[/math]

For a restricted geodesic null congruence

For a geodesic null congruence restricted on a null hypersurface, we have


[math]\displaystyle{ (20)\quad k^c\nabla_c \theta=\hat\theta_{,\,\lambda}=-\frac{1}{2}\hat\theta^2 -\hat\sigma_{ab}\hat\sigma^{ab}-\widehat{R_{cd}k^c k^d}+\kappa_{(\ell)}\hat\theta\;, }[/math]


[math]\displaystyle{ (21)\quad k^c\nabla_c \hat\sigma_{ab}=-\hat\theta\hat\sigma_{ab}+\widehat{C_{cbad}k^c k^d}+\kappa_{(\ell)}\hat\sigma_{ab}\;, }[/math]


[math]\displaystyle{ (22)\quad k^c\nabla_c \hat\omega_{ab}=0\;. }[/math]

Spin coefficients, Raychaudhuri's equation and optical scalars

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


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


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


[math]\displaystyle{ (24)\quad \theta_{(\ell)}=-(\rho+\bar\rho)=-2\text{Re}(\rho)\,,\quad \theta_{(n)}=\mu+\bar\mu=2\text{Re}(\mu)\,, }[/math]


[math]\displaystyle{ (25)\quad \sigma_{ab}=-\sigma \bar m_a \bar m_b-\bar\sigma m_a m_b\,, }[/math]


[math]\displaystyle{ (26)\quad \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(24) 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{ (27)\quad \theta_{(\ell)}=\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{ (28)\quad \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]

See also

References

  1. 1.0 1.1 1.2 Eric Poisson. A Relativist's Toolkit: The Mathematics of Black-Hole Mechanics. Cambridge: Cambridge University Press, 2004. Chapter 2.
  2. 2.0 2.1 Hans Stephani, Dietrich Kramer, Malcolm MacCallum, Cornelius Hoenselaers, Eduard Herlt. Exact Solutions of Einstein's Field Equations. Cambridge: Cambridge University Press, 2003. Chapter 6.
  3. 3.0 3.1 Subrahmanyan Chandrasekhar. The Mathematical Theory of Black Holes. Oxford: Oxford University Press, 1998. Section 9.(a).
  4. 4.0 4.1 Jeremy Bransom Griffiths, Jiri Podolsky. Exact Space-Times in Einstein's General Relativity. Cambridge: Cambridge University Press, 2009. Section 2.1.3.
  5. 5.0 5.1 P Schneider, J Ehlers, E E Falco. Gravitational Lenses. Berlin: Springer, 1999. Section 3.4.2.
  6. 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]