Optical scalars explained

In general relativity, optical scalars refer to a set of three scalar functions

\{\hat\theta

(expansion),

\hat\sigma

(shear) and

\hat\omega

(twist/rotation/vorticity)

\}

describing the propagation of a geodesic null congruence.[1] [2] [3] [4] [5]


In fact, these three scalars

\{\hat\theta,\hat\sigma,\hat\omega\}

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

\{\hat\theta\hathab,\hat\sigmaab,\hat\omegaab\}

that are adopted in tensorial equations, while the scalars

\{\hat\theta,\hat\sigma,\hat\omega\}

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

Za

, and then one could construct induced "spatial metrics" that


(1)hab=gab+ZaZb,hab=gab+ZaZb,

a
h
  b
a
=\delta
  b

+ZaZb,


where

a
h
  b
works as a spatially projecting operator. Use
a
h
  b
to project the coordinate covariant derivative

\nablabZa

and one obtains the "spatial" auxiliary tensor

Bab

,


(2)Bab

c
=h
  a

d
h
  b

\nabladZc=\nablabZa+AaZb,


where

Aa

represents the four-acceleration, and

Bab

is purely spatial in the sense that

Bab

a=B
Z
ab

Zb=0

. Specifically for an observer with a geodesic timelike worldline, we have


(3)Aa=0 ,   ⇒   Bab=\nablabZa.


Now decompose

Bab

into its symmetric and antisymmetric parts

\thetaab

and

\omegaab

,


(4)\thetaab=B(ab),\omegaab=B[ab].


\omegaab=B[ab]

is trace-free (

gab\omegaab=0

) while

\thetaab

has nonzero trace,

gab\thetaab=\theta

. Thus, the symmetric part

\thetaab

can be further rewritten into its trace and trace-free part,


(5)\thetaab=

1
3

\thetahab+\sigmaab.


Hence, all in all we have


(6)Bab=

1
3

\thetahab+\sigmaab+\omegaab,\theta=gab\thetaab=gabB(ab),\sigmaab=\thetaab-

1
3

\thetahab,\omegaab=B[ab].

For geodesic null congruences

Now, consider a geodesic null congruence with tangent vector field

ka

. Similar to the timelike situation, we also define


(7)\hat{B}ab:=\nablabka,


which can be decomposed into


(8)\hatBab=\hat\thetaab+\hat\omegaab=

1
2

\hat\theta\hathab+\hat\sigmaab+\hat\omegaab,


where


(9)\hat\thetaab=\hatB(ab),\hat\theta=\hathab\hatBab,\hat\sigmaab=\hatB(ab)-

1
2

\hat\theta\hathab,\hat\omegaab=\hatB[ab].


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

\{\hat\theta,\hat\sigma,\hat\omega\}

[1] [2] [3] [4] [5] come straightforwardly from "scalarization" of the tensors

\{\hat\theta,\hat\sigmaab,\hat\omegaab\}

in Eq(9).


The expansion of a geodesic null congruence is defined by (where for clearance we will adopt another standard symbol "

;

" to denote the covariant derivative

\nablaa

)


(10)\hat\theta=

1
2

a{}
k
;a

.

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

\theta(\ell)

and

\theta(n)

respectively, are defined by


(A.1)\theta(\ell):=hab\nablaalb,


(A.2)\theta(n):=hab\nablaanb,


where

hab=gab+lanb+nalb

represents the induced metric. Also,

\theta(\ell)

and

\theta(n)

can be calculated via


(A.3)\theta(\ell)=gab\nablaalb-\kappa(\ell),


(A.4)\theta(n)=gab\nablaanb-\kappa(n),


where

\kappa(\ell)

and

\kappa(n)

are respectively the outgoing and ingoing non-affinity coefficients defined by


(A.5)

a\nabla
l
a

lb=\kappa(\ell)lb,


(A.6)

a\nabla
n
a

nb=\kappa(n)nb.


Moreover, in the language of Newman–Penrose formalism with the convention

\{(-,+,+,+);la

a
n
a=-1,m

\bar{m}a=1\}

, we have


(A.7)\theta(l)=-(\rho+\bar\rho)=-2Re(\rho),\theta(n)=\mu+\bar\mu=2Re(\mu),


As we can see, for a geodesic null congruence, the optical scalar

\theta

plays the same role with the expansion rates

\theta(\ell)

and

\theta(n)

. Hence, for a geodesic null congruence,

\theta

will be equal to either

\theta(\ell)

or

\theta(n)

.


The shear of a geodesic null congruence is defined by


(11){\hat\sigma}

2=\hat\sigma
ab

\hat{\bar\sigma}ab=

1
2

gcagdbk(a;b)kc;d-(

1
2

a{}
k
;a

)2=gcagdb

1
2

k(a;b)kc;d-{\hat\theta}2 .


The twist of a geodesic null congruence is defined by


(12){\hat\omega}2=

1
2

k[a;b]ka;b=gcagdbk[a;b]kc;d.


In practice, a geodesic null congruence is usually defined by either its outgoing (

ka=la

) or ingoing (

ka=na

) tangent vector field (which are also its null normals). Thus, we obtain two sets of optical scalars

\{\hat\theta(\ell),\hat\sigma(\ell),\hat\omega(\ell)\}

and

\{\hat\theta(n),\hat\sigma(n),\hat\omega(n)\}

, which are defined with respect to

la

and

na

, respectively.

Applications in decomposing the propagation equations

For a geodesic timelike congruence

The propagation (or evolution) of

Bab

for a geodesic timelike congruence along

Zc

respects the following equation,


(13)

c\nabla
Z
c

Bab

c
=-B
  b

Bac+RcbadZcZd.


Take the trace of Eq(13) by contracting it with

gab

, and Eq(13) becomes


(14)

c\nabla
Z
c

\theta=\theta,\tau=-

1
3

\theta2-\sigmaab\sigmaab+\omegaab\omegaab-RabZaZb


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


(15)

c\nabla
Z
c

\sigmaab=-

2
3

\theta\sigmaab-\sigmaac

c
\sigma
b

-\omegaac

c
\omega+
b
1
3

hab(\sigmacd\sigmacd-\omegacd\omegacd)+CcbadZc

d+1
2
Z

\tilde{R}ab.


Finally, the antisymmetric component of Eq(13) yields


(16)

c\nabla
Z
c

\omegaab=-

2
3

\theta\omegaab

c
-2\sigma
[b

\omegaa]c.

For a geodesic null congruence

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


(16)

c\nabla
k
c

\hatBab=-\hat

c
B
  b

\hatBac+\widehat{Rcbadkckd}.


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


(17)

c\nabla
k
c

\hat\theta=\hat\theta,λ=-

1
2

\hat\theta2-\hat\sigmaab\hat\sigmaab+\hat\omegaab\hat\omegaab-\widehat{Rcdkckd},


(18)

c\nabla
k
c

\hat\sigmaab=-\hat\theta\hat\sigmaab+\widehat{Ccbadkckd},


(19)

c\nabla
k
c

\hat\omegaab=-\hat\theta\hat\omegaab.

For a restricted geodesic null congruence

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


(20)

c\nabla
k
c

\theta=\hat\theta,λ=-

1
2

\hat\theta2-\hat\sigmaab\hat\sigmaab-\widehat{Rcdkc

d}+\kappa
k
(\ell)

\hat\theta,


(21)

c\nabla
k
c

\hat\sigmaab=-\hat\theta\hat\sigmaab+\widehat{Ccbadkc

d}+\kappa
k
(\ell)

\hat\sigmaab,


(22)

c\nabla
k
c

\hat\omegaab=0 .

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


(23)l{L}\ell\theta(\ell)=-

1
2
2+\tilde{\kappa}
\theta
(\ell)

\theta(\ell)-\sigmaab\sigmaab+\tilde{\omega}ab\tilde{\omega}ab-Rablalb,


where

\tilde{\kappa}(\ell)

is defined such that

\tilde{\kappa}(\ell)lb:=la\nablaalb

. The quantities in Raychaudhuri's equation are related with the spin coefficients via


(24)\theta(\ell)=-(\rho+\bar\rho)=-2Re(\rho),\theta(n)=\mu+\bar\mu=2Re(\mu),


(25)\sigmaab=-\sigma\barma\barmb-\bar\sigmamamb,


(26)\tilde{\omega}ab=

1
2

(\rho-\bar\rho)(ma\barmb-\barmamb)=Im(\rho)(ma\barmb-\barmamb),


where Eq(24) follows directly from

\hat{h}ab=\hat{h}ba=mb\barma+\barmbma

and


(27)\theta(\ell)=\hat{h}ba\nablaa

b\bar
l
b=m
a\nabla
m
a

lb+\barmb

a\nabla
m
a

lb=mb\bar\deltalb+\barmb\deltalb=-(\rho+\bar\rho),


(28)\theta(n)=\hat{h}ba\nablaanb=\barmb

a\nabla
m
a
b\bar
n
b+m
a\nabla
m
a

nb=\barmb\delta

b\bar
n
b+m

\deltanb=\mu+\bar\mu.

See also

Notes and References

  1. Eric Poisson. A Relativist's Toolkit: The Mathematics of Black-Hole Mechanics. Cambridge: Cambridge University Press, 2004. Chapter 2.
  2. 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. Subrahmanyan Chandrasekhar. The Mathematical Theory of Black Holes. Oxford: Oxford University Press, 1998. Section 9.(a).
  4. Jeremy Bransom Griffiths, Jiri Podolsky. Exact Space-Times in Einstein's General Relativity. Cambridge: Cambridge University Press, 2009. Section 2.1.3.
  5. 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]