Next Article in Journal
Measuring the Lense–Thirring Orbital Precession and the Neutron Star Moment of Inertia with Pulsars
Previous Article in Journal
A Needle in a Cosmic Haystack: A Review of FRB Search Techniques
Previous Article in Special Issue
Towards an Acoustic Geometry in Slightly Viscous Fluids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boyer–Lindquist Space-Times and Beyond: Metamaterial Analogues for Arbitrary Space-Times

by
Sebastian Schuster
1,2,† and
Matt Visser
3,*,†
1
Ústav Teoretické Fyziky, Matematicko-Fyzikální Fakulta, Univerzita Karlova, V Holešovičkách 2, 180 00 Prague, Czech Republic
2
School of Mathematics and Statistics, University of Sheffield, Hounsfield Road, Sheffield S3 7RH, UK
3
School of Mathematics and Statistics, Victoria University of Wellington, P.O. Box 600, Wellington 6140, New Zealand
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2024, 10(4), 159; https://doi.org/10.3390/universe10040159
Submission received: 31 October 2022 / Revised: 27 February 2024 / Accepted: 19 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Analogue Gravity)

Abstract

:
Analogue space-times (and in particular metamaterial analogue space-times) have a long varied and rather complex history. Much of the previous related work to this field has focused on spherically symmetric models; however, axial symmetry is much more relevant for mimicking astrophysically interesting systems that are typically subject to rotation. Now it is well known that physically reasonable stationary axisymmetric space-times can, under very mild technical conditions, be put into Boyer–Lindquist form. Unfortunately, a metric presented in Boyer–Lindquist form is not well adapted to the “quasi-Cartesian” metamaterial analysis that we developed in our previous articles on “bespoke analogue space-times”. In the current article, we shall first focus specifically on various space-time metrics presented in Boyer–Lindquist form, and subsequently determine a suitable set of equivalent metamaterial susceptibility tensors in a laboratory setting. We shall then turn to analyzing generic space-times, not even necessarily stationary, again determining a suitable set of equivalent metamaterial susceptibility tensors. Perhaps surprisingly, we find that the well-known ADM formalism proves to be not particularly useful, and that it is instead the dual “threaded” (Kaluza–Klein–inspired) formalism that provides much more tractable results. While the background laboratory metric is (for mathematical simplicity and physical plausibility) always taken to be Riemann flat, we will allow for arbitrary curvilinear coordinate systems on the flat background space-time. Finally, for completeness, we shall reconsider spherically symmetric space-times, but now in general spherical polar coordinates rather than quasi-Cartesian coordinates. In summary, this article provides a set of general-purpose calculational tools that can readily be adapted for mimicking various interesting (curved) space-times by using nontrivial susceptibility tensors in general (background-flat) laboratory settings.

1. Introduction

In three earlier articles [1,2,3], the present authors have carefully reanalyzed and re-explored the notion of electromagnetic analogue space-times. In all of these metamaterial analogues, we were very careful to distinguish the background (flat) laboratory metric from the effective (curved space-time) analogue metric.
In the first article [1], we readdressed the one-century-old question of just when a (possibly moving) electromagnetic medium is equivalent (in a purely algebraic sense [3]) to an effective Lorentzian metric [4]. In a simplified setting (static with appropriate coordinate choices), Landau and Lifschitz [5], in their influential book “The classical theory of fields”, write (§90, p. 257 of the 3rd edition, 1971):
The reader should note the analogy (purely formal, of course) of [the Maxwell equations in a gravitational field] to the Maxwell equations for the electromagnetic field in material media. In particular, in a static gravitational field [the constitutive relations] reduce to D = E / h , B = H / h . We may say that with respect to its effect on the electromagnetic field a static gravitational field plays the role of a medium with electric and magnetic permeabilities ϵ = μ = 1 / h .
In the discussion below, we shall go far beyond this simplified setting, allowing for generic gravitational fields. That is, the electromagnetic medium is taken to be characterized by effectively three-dimensional permittivity ϵ , permeability μ , magnetoelectric ζ tensors, and 4-velocity V a . We used these four quantities to construct an effective Lorentzian metric—an analogue space-time. Specifically, we explicitly constructed (up to an undetermined conformal factor) the effective space-time metric in terms of the optical tensors and vice versa.
In the second article [2], we explicitly constructed the flat-space laboratory (Cartesian coordinate) susceptibility tensors appropriate for mimicking the Schwarzschild geometry in various quasi-Cartesian forms. (Specifically, we analyzed curvature, isotropic, Kerr–Schild, Painlevé–Gullstrand, and Gordon forms of the Schwarzschild geometry). We also analyzed general static spherically symmetric space-times, and the Kerr geometry [6,7,8,9], in both Kerr–Schild and Doran forms [10,11] (both of which are easily put into quasi-Cartesian form [2]).
In the third article [3], we contrasted the essentially algebraic nature of the present electromagnetic analogue space-time model with other potential electromagnetic analogues, which would be more in line with the traditional, analytic approach: Instead of building on tensorial similarities in the action (see the next section), one would traditionally look for a specific differential equation (usually the Laplace–Beltrami equation on a curved background) for making comparisons between electromagnetism and gravitation. This also allowed us to demonstrate simple cartographic effects on the analogue Hawking temperature that are to be expected from simple spatial rescaling. This furthermore links back to existing work on cartography in electromagnetic space-time analogues [12,13].
Unfortunately, working with the Kerr geometry (or more generally, arbitrary axially symmetric space-times) in Boyer–Lindquist coordinates is technically much messier than working in Kerr–Schild and Doran forms, or in spherical symmetry. The difficulty arises because there is no longer any natural unit-determinant quasi-Cartesian form for the metric, and we shall turn to this topic in the current article.
One particular reason that it is desirable to work with space-times in Boyer–Lindquist form is that there is then only one off-diagonal metric component. In contrast, dealing with Kerr–Schild or Doran forms, there is an explicit trade-off: While Kerr–Schild or Doran forms can easily be put in quasi-Cartesian form (so that the background (flat) laboratory metric can be put in Cartesian form), they have multiple off-diagonal elements in the metric.
Therefore, the Boyer–Lindquist form minimizes the number of off-diagonal components in the effective analogue metric, at the cost of making the background laboratory metric trickier to deal with. The background laboratory metric is still Riemann flat, just presented in unusual coordinates. (Some parts of the calculation to be presented below are much simpler than the general discussion in [1], or the quasi-Cartesian discussion in [2]; other parts of the discussion are considerably more subtle).
After first dealing with generic Boyer–Lindquist space-times, we turn to arbitrary space-times, not even necessarily stationary, and perform a similar analysis. While the background laboratory metric is always taken to be Riemann flat, we will now focus on arbitrary curvilinear coordinate systems.
Finally, we shall reanalyze static spherically symmetric space-times eschewing quasi-Cartesian coordinates. Whereas quasi-Cartesian coordinates are particularly useful for discussing laboratory physics, certain purely theoretical computations can be more clearly carried out using spherical polar coordinates. This means that one has to keep track of some scalar densities arising from metric determinants (both effective metric and background metric determinants) and be much more careful when raising and lowering indices.
For further background on these topics, see our three earlier articles [1,2,3], the very early 1923 article by Gordon [4], and the early discussion by Landau and Lifshitz [5]. Other relevant, relatively early articles come from the general relativity community [14,15,16,17,18,19,20,21]. More recently, the optics community has developed the closely related notion of “transformation optics”, leading to an extensive body of work [12,13,22,23,24,25,26,27]. The electromagnetic analogue space-times are a natural complement to the acoustic analogue space-times of [28,29,30,31,32,33,34,35,36,37]. For a much more general background and history, see [38,39,40,41,42,43].

2. General Strategy

Consider electromagnetism in curved space-time. Let the coordinates be x a and take the effective metric to be g a b . We shall explicitly write the contravariant inverse effective metric as [ g 1 ] a b . Then the curved space-time electromagnetic Lagrangian is [1,3]:
S space - time = 1 8 det ( g a b ) [ g 1 ] a c [ g 1 ] b d [ g 1 ] a d [ g 1 ] b c F a b F c d d 4 x .
We wish to mimic this in some laboratory setting, using the same coordinates but now a laboratory (background) metric ( g 0 ) a b —which will typically be flat space-time, but possibly in unusual coordinates. For convenience, we shall assume that the laboratory (background) metric takes the block-diagonal time+space form:
[ g 0 ] a b = 1 0 0 [ g 0 ] i j ( ξ ) ; det ( [ g 0 ] a b ) = det ( [ g 0 ] i j ) .
Then for a general (linear) metamaterial,
S laboratory = 1 2 det ( [ g 0 ] a b ) ϵ i j E i E j + 2 ζ i j B i E j [ μ 1 ] i j B i B j d 4 x .
Equivalently, to be very explicit,
S laboratory = 1 2 + det ( [ g 0 ] i j ) ϵ i j E i E j + 2 ζ i j B i E j [ μ 1 ] i j B i B j d t d 3 ζ .
We wish to match up these two actions, S space - time and S laboratory , to develop an equivalent laboratory analogue to the curved space-time—at least insofar as it comes to the behavior of electromagnetic fields. Specifically, we wish to mimic the effective metric g a b (up to some undetermined conformal factor) by using specified susceptibility tensors { ϵ i j , [ μ 1 ] i j , ζ i j } and a suitable background metric [ g 0 ] a b for laboratory physics. (Conformal invariance of the susceptibility tensors { ϵ i j , [ μ 1 ] i j , ζ i j } is an automatic consequence of the conformal invariance of electromagnetism in (3 + 1) dimensions. In contrast, the conformal invariance of the Hawking temperature is a somewhat deeper result, applicable not just to electromagnetism [44]). Indices will always be raised and lowered using the laboratory metric g 0 (which is why we need to use the explicit notation [ g 1 ] a b for the inverse of the effective metric we want to mimic).
We now consider the constitutive tensor [1]:
Z a b c d = 1 2 det ( g a b ) det ( [ g 0 ] a b ) [ g 1 ] a c [ g 1 ] b d [ g 1 ] a d [ g 1 ] b c ,
which mimics the electromagnetic properties of the effective metric g a b . (Note that Z a b c d is, by construction, a true tensor, not a tensor density). When no confusion can arise, we simplify det ( g a b ) det ( g ) and det ( [ g 0 ] a b ) det ( g 0 ) . The laboratory permittivity, permeability, and magnetoelectric tensors { ϵ i j , [ μ 1 ] i j , ζ i j } are then [1]:
ϵ i j = 2 Z i 0 j 0 ; [ μ 1 ] i j = 1 2 ε i k l ε j m n Z k l m n ; ζ i j = ε i k l Z k l j 0 .
The ε i k l appearing here have to be interpreted as three-dimensional Levi–Civita pseudo-tensors defined in terms of the background spatial 3-metric [ g 0 ] i j and its metric density det ( [ g 0 ] i j ) , with indices raised and lowered using the background spatial 3-metric. This guarantees that ϵ i j and [ μ 1 ] i j are true T 0 2 tensors under spatial coordinate transformations, while ζ i j is a pseudo-tensor under parity reversal. (The computations are much simpler if the laboratory metric is Cartesian and the space-time metric quasi-Cartesian; see reference [2]). The net result of this procedure is that
Z a b c d = Z a b c d ( ε i j , [ μ 1 ] i j , ζ i j , [ g 0 ] i j )
can be written as an explicit function of the optical tensors, while conversely, the optical tensors
{ ϵ i j , [ μ 1 ] i j , ζ i j } = f ( g a b , [ g 0 ] a b )
can be written as explicit functions of the effective metric g a b and background metric [ g 0 ] a b . This implies a stringent compatibility condition (linking the permittivity, permeability, and magnetoelectric tensors) that must be satisfied in order for the analogy to work in the sense of Gordon [4]. For scalar permittivity and permeability, this gives the compatibility condition ϵ = μ of the above quote from Landau and Lifshitz [5]. The tensorial generalization to ϵ i j = μ i j was then rederived several times in the literature with different inspirations [16,45,46]. Even the most general cases can always (at least locally) be reduced to this case by appropriate frame changes [1]. This is the general algorithm—some of the intermediate stages in the specific computations below are somewhat tedious—but the final results are clean, useful, and easy to interpret.

3. Boyer–Lindquist Space-Times

Essentially, all physically interesting stationary axisymmetric space-times (not just Kerr or Kerr–Newman space-times) can be put into Boyer–Lindquist form. Adopting quasi-spherical-polar coordinates, with the coordinates ordered as ( t , r , θ , ϕ ) , with θ and ϕ having the usual interpretation as polar and azimuthal angles, respectively, it is sufficient to note that without significant loss of generality, any stationary axisymmetric geometry can be put in the form [47,48]:
g a b = g t t 0 0 g t ϕ 0 g r r 0 0 0 0 g θ θ 0 g t ϕ 0 0 g ϕ ϕ .
Here, the metric components depend only on ( r , θ ) and are independent of ( t , ϕ ) . The inverse metric is easily computed:
[ g 1 ] a b = g ϕ ϕ / g 2 0 0 g t ϕ / g 2 0 1 / g r r 0 0 0 0 1 / g θ θ 0 g t ϕ / g 2 0 0 g t t / g 2 = [ g 1 ] 00 [ g 1 ] 0 j [ g 1 ] 0 i [ g 1 ] i j .
Here, g 2 = g t t g ϕ ϕ g t ϕ 2 , and det ( g a b ) = g 2 g r r g θ θ . Note that g t t = 0 at any ergo-surface, while it is easy to convince oneself that g 2 = 0 at any horizon; equivalently, g r r = at any horizon. (See, for instance, references [6,9]).
The “laboratory” will be assumed to be flat Minkowski space, in coordinates that we will label ( t , r , θ , ϕ ) = ( t , ξ i ) . We assume that the background metric takes the following block-diagonal form:
[ g 0 ] a b = 1 0 0 [ g 0 ] i j ( ξ ) , det ( [ g 0 ] a b ) = det ( [ g 0 ] i j ) .
Here, [ g 0 ] i j ( ξ ) is the 3-metric corresponding to some (essentially arbitrary) coordinate representation of flat 3-space. For instance, these coordinates might be spherical polar, oblate spheroidal, prolate spheroidal, or cylindrical coordinates (in which case, one would more likely relabel θ z ), or something even more exotic (for example, parabolic cylindrical, paraboloidal, elliptic cylindrical, ellipsoidal, bipolar, toroidal, conical, or general orthogonal coordinates). Let us now explicitly calculate the laboratory permittivity, permeability, and magnetoelectric tensors { ϵ i j , [ μ 1 ] i j , ζ i j } .

3.1. Permittivity Tensor

We start by noting the following:
ϵ i j = det ( g ) det ( g 0 ) [ g 1 ] i j [ g 1 ] 00 [ g 1 ] 0 i [ g 1 ] 0 j = det ( g ) det ( g 0 ) 1 g 2 g ϕ ϕ / g r r 0 0 0 g ϕ ϕ / g θ θ 0 0 0 1 i j .
Now in the physically interesting region (in the domain of outer communication, outside the outermost noncosmological horizon), both det ( g ) and g 2 are negative, while det ( g 0 ) is negative everywhere. This allows us to rewrite ϵ i j as:
ϵ i j = g r r g θ θ g 2 det ( g 0 ) g ϕ ϕ / g r r 0 0 0 g ϕ ϕ / g θ θ 0 0 0 1 i j .
Note that this is manifestly conformally invariant (under conformal rescaling of the effective space-time metric g a b one is mimicking) as it must be (due to the conformal invariance of electromagnetism in (3 + 1) dimensions); see, for instance, references [1,2].

3.2. Permeability Tensor

We start by noting the following:
[ μ 1 ] i j = 1 2 ε i k l ε j m n Z k l m n = 1 2 ε i k l ε j m n det ( g ) det ( g 0 ) [ g 1 ] k m [ g 1 ] l n .
Now in terms of the Levi–Civita tensor density ε ¯ i j k = signum ( i j k ) , we can write the true tensor identities:
ε i j k = det ( g 0 ) ε ¯ i j k ,
and
ε i j k = [ g 0 ] i p det ( g 0 ) ε ¯ p j k ,
whence
[ μ 1 ] i j = 1 2 [ g 0 ] i p [ g 0 ] j q [ det ( g 0 ) ] det ( g ) det ( g 0 ) ε ¯ p k l ε ¯ q m n [ g 1 ] k m [ g 1 ] l n .
Now we know that [ g 1 ] l n is diagonal, so ε ¯ p k l ε ¯ q m n [ g 1 ] k m [ g 1 ] l n is also guaranteed to be diagonal. Then,
[ μ 1 ] i j = [ g 0 ] i p [ g 0 ] j q det ( g ) det ( g 0 ) g t t / ( g θ θ g 2 ) 0 0 0 g t t / ( g r r g 2 ) 0 0 0 1 / ( g r r g θ θ ) p q .
(Note in passing the conformal invariance under rescaling of the effective metric g a b ).
Matrix inversion (and raising indices using the background metric g 0 ) now yields:
μ i j = g r r g θ θ g 2 det ( g 0 ) g 2 / ( g t t g r r ) 0 0 0 g 2 / ( g t t g θ θ ) 0 0 0 1 i j .
Note that this is conformally invariant (under conformal rescaling of the effective metric g a b one is mimicking) as it must be (due to the conformal invariance of electromagnetism in (3 + 1) dimensions); see, for instance, references [1,2].

3.3. Magnetoelectric Tensor

We start by noting
ζ i j = 1 2 det ( g ) det ( g 0 ) ε i k l [ g 1 ] 0 l [ g 1 ] k j .
Inserting the explicit expressions for [ g 1 ] 0 l and ε i k ϕ and simplifying gives
ζ i j = 1 2 g t ϕ g r r g θ θ g 2 0 1 / g θ θ 0 1 / g r r 0 0 0 0 0 i j .
Note again the conformal invariance under rescaling of the effective metric g a b that one is trying to mimic. Note further that the magnetoelectric tensor is proportional to the only off-diagonal metric component g t ϕ . Consequently, magnetoelectric effects vanish in the absence of rotation.

3.4. Summary for Metrics in Boyer–Lindquist Coordinates

Collecting results, for Boyer–Lindquist space-times, we have the following:
ϵ i j = g r r g θ θ g 2 det ( g 0 ) g ϕ ϕ / g r r 0 0 0 g ϕ ϕ / g θ θ 0 0 0 1 i j .
μ i j = g r r g θ θ g 2 det ( g 0 ) g 2 / ( g t t g r r ) 0 0 0 g 2 / ( g t t g θ θ ) 0 0 0 1 i j .
ζ i j = 1 2 g t ϕ g r r g θ θ g 2 0 1 / g θ θ 0 1 / g r r 0 0 0 0 0 i j .
While the abstract computation required to get to this stage has been slightly tedious, the final results are fully explicit, and quite general. A number of interesting implications can immediately be read off.

3.5. Implications for Metrics in Boyer–Lindquist Coordinates

  • First, note that ϵ ϕ ϕ = μ ϕ ϕ . (This is related to what we saw happened for spherical symmetry in reference [2]. The general point is that electromagnetic properties in the direction of the 3-vector [ g 1 ] 0 i are degenerate). Indeed, all the components of the permittivity tensor ϵ i j are well defined down to the outermost noncosmological horizon ( g 2 = 0 , g r r = ).
  • Second, note that while μ ϕ ϕ is well defined all the way to the outermost noncosmological horizon, μ r r and μ θ θ are only well defined down to the outermost noncosmological ergo-surface (where g t t = 0 ).
  • Third, note that the magnetoelectric tensor ζ i j is well defined down to the outermost noncosmological horizon ( g 2 = 0 , g r r = ).
  • Fourth, note that the difference ϵ i j μ i j is relatively simple:
    ϵ i j μ i j = g r r g θ θ g 2 det ( g 0 ) g t ϕ 2 g t t 1 / g r r 0 0 0 1 / g θ θ 0 0 0 0 i j .
    Observe that as g t ϕ 0 (that is, as the rotation is switched off), we find ϵ i j = μ i j and ζ i j = 0 , the standard compatibility condition for static space-times [1,2].
  • Observe that the magnetoelectric tensor always has a zero-eigenvalue eigenvector, currently in the ϕ direction, and so det ( ζ i j ) = 0 for all Boyer–Lindquist mimics.
  • Observe that:
    ( ζ 2 ) i k = ζ i j ζ j k = 1 4 g t ϕ 2 g 2 1 0 0 0 1 0 0 0 0 i k .
    This is actually proportional to a projection operator onto the directions perpendicular to the 3-vector [ g 1 ] 0 i . (We saw similar things happen in the quasi-Cartesian analysis of reference [2]).
  • Observe that:
    tr ( ζ 2 ) = ζ i j ζ j i = 1 2 g t ϕ 2 g 2 .
    This is a nice scalar invariant describing the strength of the magnetoelectric effect, well defined down to the outermost noncosmological horizon (where g 2 = 0 ).
  • Finally, note that, for small g t ϕ , that is, | g t ϕ |   | g t t g ϕ ϕ | , we can explicitly write the following:
    ϵ i j = det ( g ) det ( g 0 ) 1 g t t 1 / g r r 0 0 0 1 / g θ θ 0 0 0 1 / g ϕ ϕ i j 1 + O ( g t ϕ 2 ) = μ i j .
    ζ i j = 1 2 g t ϕ g r r g θ θ g t t g ϕ ϕ 0 1 / g θ θ 0 1 / g r r 0 0 0 0 0 i j 1 + O ( g t ϕ 2 ) .

4. Specific Concrete Examples

Let us now consider some specific concrete examples.

4.1. Kerr

The Kerr space-time represents a rotating vacuum black hole in standard general relativity [6,7,8,9]. We have already investigated the Cartesian Kerr–Schild and Cartesian Doran forms of the Kerr metric in reference [2], but now turn our attention to the Boyer–Lindquist form of the Kerr metric. Defining as usual
Δ = r 2 2 m r + a 2 ; Σ = r 2 + a 2 cos 2 θ ;
the Kerr line element can be written as [6,9]:
d s 2 = 1 2 m r Σ d t 2 4 m a r sin 2 θ Σ d ϕ d t + r 2 + a 2 + 2 m a 2 r sin 2 θ Σ sin 2 θ d ϕ 2 + Σ Δ d r 2 + Σ d θ 2 .
For convenience, we can explicitly take the background metric to be the m = 0 limit of Kerr, which is known to be Riemann-flat Minkowski space in oblate spheroidal coordinates [6,9]:
( d s 2 ) 0 = d t 2 + r 2 + a 2 sin 2 θ d ϕ 2 + Σ r 2 + a 2 d r 2 + Σ d θ 2 .
Here, the free parameter a specifies the distance of the oblate spheroid’s two foci from the origin, and for a 0 , one regains spherical coordinates with a single focus at the origin. Then,
g 2 Δ sin 2 θ ; det ( g a b ) Σ 2 sin 2 θ ; det ( [ g 0 ] a b ) Σ 2 sin 2 θ .
Furthermore,
g r r g θ θ g 2 det ( g 0 ) 1 Δ sin 2 θ ; g r r g θ θ g 2 Σ Δ sin θ .
Collecting results, for Kerr space-time, we have the following:
ϵ i j = 1 Δ Σ Δ r 2 + a 2 + 2 m a 2 r sin 2 θ Σ 0 0 0 r 2 + a 2 + 2 m a 2 r sin 2 θ Σ 0 0 0 Σ sin 2 θ i j ,
and
μ i j = Δ Σ 2 m r 0 0 0 1 Σ 2 m r 0 0 0 1 Δ sin 2 θ i j .
Finally, for the magnetoelectric tensor, we have
ζ i j = m a r sin θ Σ 0 1 / Δ 0 1 0 0 0 0 0 i j .

4.2. Lense–Thirring

The Lense–Thirring slow rotation approximation to Kerr has a long and complex history in its own right [49,50,51]. However, for current purposes, it is sufficient to take a < m r and Taylor expand. Explicitly keeping the zeroth-order and first-order terms in a, one has the following:
d s 2 = ( 1 2 m / r ) d t 2 + d r 2 1 2 m / r + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) 4 m a sin 2 θ r d ϕ d t + O ( a 2 ) ,
with the electromagnetic tensors simplifying to the following:
ϵ i j = 1 0 0 0 1 r 2 ( 1 2 m / r ) 0 0 0 1 r 2 ( 1 2 m / r ) sin 2 θ i j + O ( a 2 ) = g i j 1 2 m / r + O ( a 2 ) .
μ i j = 1 0 0 0 1 r 2 ( 1 2 m / r ) 0 0 0 1 r 2 ( 1 2 m / r ) sin 2 θ i j + O ( a 2 ) = g i j 1 2 m / r + O ( a 2 ) .
ζ i j = m a sin θ r 0 1 r 2 ( 1 2 m / r ) 0 1 0 0 0 0 0 i j + O ( a 2 ) .
That is,
ζ i j = m a sin θ r 0 1 0 1 0 0 0 0 0 i j + O ( a 2 ) .
Note that this is just what one would expect for Schwarzschild (in curvature coordinates), plus a first-order magnetoelectric effect to mimic slow rotation. Additionally, note that the naming conventions of electromagnetism and its constitutive tensors differ here compared with the conventions commonly adopted in gravitomagnetism; the closeness of the Lense–Thirring metric to the latter warrants an extra warning concerning this point.

5. Going beyond Boyer–Lindquist: Arbitrary Space-Times

A surprise for arbitrary space-times is that the usual ADM formalism proves to be less than useful, whereas the “threaded” (Kaluza–Klein–inspired) form of the metric is more amenable to developing a metamaterial analogue. To now proceed beyond Boyer–Lindquist space-times, we first write the metric to be mimicked in threaded form [52,53,54] (also known as Kaluza–Klein form [1,2]):
[ g 1 ] a b = α 2 + γ k l 1 β k β l β j β i γ i j ; det ( g ) = α 2 det ( γ 1 ) .
Equivalently, for the covariant metric, one has:
[ g ] a b = α 2 α 2 β j α 2 β i [ γ 1 ] i j α 2 β i β j ; β i = [ γ 1 ] i k β k .
Note that the entries of the covariant and contravariant metrics are exactly swapped compared with the ADM formalism, hence the “duality” of Kaluza–Klein and ADM.
This is not at all a restriction on the metric, merely a convenient way of writing it. (We have not enforced the unimodular condition of reference [2] since we are now not using quasi-Cartesian coordinates, and the unimodular condition is now neither necessary nor useful). Note that because the analysis below is purely algebraic, there is no need to demand stationarity (time independence)—such an assumption may be convenient, but it is not necessary.

5.1. Permittivity Tensor

We start by noting the following:
ϵ i j = det ( g ) det ( g 0 ) [ g 1 ] i j [ g 1 ] 00 [ g 1 ] 0 i [ g 1 ] 0 j
= + det ( g ) det ( g 0 ) ( α 2 γ k l 1 β k β l ) γ i j + β i β j .
We can write this as:
ϵ i j = 1 α | det ( g 0 ) | det ( γ p q ) ( 1 α 2 [ γ 1 ] k l β k β l ) γ i j + α 2 β i β j
Note that this is a true tensor equation (under arbitrary spatial coordinate transformations), and that it is conformally invariant (under conformal rescaling of g a b , the effective metric to be mimicked).

5.2. Permeability Tensor

We start by noting the following:
[ μ 1 ] i j = 1 2 ε i k l ε j m n Z k l m n
= 1 2 [ g 0 ] i p [ g 0 ] j q [ det ( g 0 ) ] det ( g ) det ( g 0 ) ε ¯ p k l ε ¯ q m n γ k m γ l n .
Lowering indices (using the background metric),
[ μ 1 ] i j = 1 2 [ det ( g 0 ) ] det ( g ) det ( g 0 ) ε i k l ε j m n γ k m γ l n .
Now purely as a matter of algebra, for 3 × 3 matrices, we have
ε ¯ i k l ε ¯ j m n γ k m γ l n = 2 det ( γ p q ) [ γ 1 ] i j .
Using this, together with a matrix inversion, now yields
μ i j = 1 α | det ( g 0 ) | det ( γ p q ) γ i j .
Note that this is a true tensor equation (under arbitrary spatial coordinate transformations), and that it is conformally invariant (under conformal rescaling of g a b , the effective metric to be mimicked).

5.3. Magnetoelectric Tensor

Lowering the first index, the best we can do in full generality is (after a very short calculation) as follows:
ζ i j = 1 2 det ( g ) det ( g 0 ) ε i k l β l γ k j .
Note again the conformal invariance under rescaling of g a b . Furthermore, this is a true pseudo-tensor equation (under spatial coordinate transformations).

5.4. Summary for Arbitrary Space-Times

Collecting results, for a generic space-time, we have the following:
ϵ i j = 1 α | det ( g 0 ) | det ( γ p q ) ( 1 α 2 γ k l 1 β k β l ) γ i j + α 2 β i β j
μ i j = 1 α | det ( g 0 ) | det ( γ p q ) γ i j .
ζ i j = 1 2 det ( g ) det ( g 0 ) ε i k l β l γ k j .
Again, getting to this stage has been slightly tedious, but the final results are fully explicit and relatively simple. Several interesting implications immediately follow.

5.5. Implications for Arbitrary Space-Times

  • First, note that the difference ϵ i j μ i j is quite simple:
    ϵ i j μ i j = α | det ( g 0 ) | det ( γ p q ) ( γ k l 1 β k β l ) γ i j β i β j .
    Observe that as β i 0 , we again find ϵ i j = μ i j (and ζ i j = 0 ), the standard compatibility condition for static space-times [1,2].
  • Second, note that if we define | | β | | 2 = ( γ k l 1 β k β l ) and β ^ i = β i / | | β | | , then
    ϵ i j μ i j = α | | β | | 2 | det ( g 0 ) | det ( γ p q ) γ i j β ^ i β ^ j .
    Then in the direction of the 3-vector [ γ 1 ] j k β k , we have
    ( ϵ i j μ i j ) [ γ 1 ] j k β k = 0 .
    (This is related to what we saw happen for Boyer–Lindquist above (where we found ϵ ϕ ϕ = μ ϕ ϕ ), and also for spherical symmetry in [2]. The general point is that electromagnetic properties in the direction of the 3-vector β i = [ g 1 ] 0 i are degenerate).
  • For the magnetoelectric tensor, observe that β i ζ i j = 0 , so the direction β i is again special. This implies that the magnetoelectric tensor always has a left eigenvector with eigenvalue zero, and det ( ζ i j ) = 0 . (This is a general property of electromagnetic media with a single light cone. See [19], p. 282). The corresponding right eigenvector with eigenvalue zero is [ γ 1 ] j k β k . In fact, the magnetoelectric tensor can always be written in the form ζ i k = A i j S j k with A i j antisymmetric and S j k symmetric 3 × 3 matrices, which makes the existence of these eigenvectors obvious.
  • Observe that
    ( ζ 2 ) i j = ζ i k ζ k j = 1 4 det ( g ) det ( g 0 ) ε i p l β l γ p k ε k m n β n γ m j .
    Hence,
    tr ( ζ 2 ) = 1 2 α 2 [ γ 1 ] i j β i β j .
    This is a simple scalar invariant (under spatial coordinate transformations) describing the strength of the magnetoelectric effect.
  • If we choose to rotate our coordinate system so that β i = ( 0 , 0 , β ) , then things simplify somewhat. The eigenvalues of ζ i j are then proportional to ( 0 , + β , β ) , and consequently, the eigenvalues of ( ζ 2 ) i j are then proportional to ( 0 , + β 2 , + β 2 ) . Indeed,
    ( ζ 2 ) i j = tr ( ζ 2 ) 2 1 0 0 1 0 0 0 i j ,
    which is almost a projection operator.

6. Revisiting Spherically Symmetric Space-Times in Spherical Polar Coordinates

Finally, it is worth revisiting the analysis of reference [2] for static spherically symmetric space-times, but now eschewing the use of quasi-Cartesian coordinates, and allowing for time dependence. While the quasi-Cartesian coordinates of reference [2] correspond to Cartesian coordinates for the background metric describing the laboratory, and so are particularly useful for presentational purposes when phrasing laboratory-based questions, sometimes explicit computations are more cleanly carried out in spherical polar coordinates.
There is, however, a price to be paid: one has to keep track of some scalar densities arising from metric determinants (both effective and background) and be much more careful raising and lowering indices.
First, let us adopt ( t , r , θ , ϕ ) coordinates, and write the flat background metric in the following form:
( d s 0 ) 2 = d t 2 + [ R ( r ) d r ] 2 + R ( r ) 2 { d θ 2 + sin 2 θ d ϕ 2 } .
Note det ( [ g 0 ] a b ) = R ( r ) 2 R ( r ) 4 sin 2 θ 1 ; the background metric is not unimodular. This is the most general form of a flat space-time metric compatible with explicit spherical symmetry in the sense of being based on spherical polar coordinates. Then for the metric to be mimicked, g a b we can write without any loss of generality:
( d s ) 2 = g t t d t 2 + 2 g t r d t d t + g r r d r 2 + R ( r ) 2 { d θ 2 + sin 2 θ d ϕ 2 } .
Note that we are using the same coordinates in both the laboratory and the metric to be mimicked, and we are making extensive use of the assumed spherical symmetry. We are also keeping the metric to be mimicked in as general a form as possible—so that we can simultaneously deal with curvature coordinates (where R ( r ) = r ), isotropic coordinates (where g r r = R ( r ) 2 ), or various off-diagonal coordinates (such as Kerr–Schild coordinates, Painlevé–Gullstrand coordinates, or Gordon coordinates [42,43,55]).
To proceed, we again rewrite the metric to be mimicked in threaded form:
[ g 1 ] a b = α 2 + β 2 / γ β 0 0 β γ 0 0 0 0 R 2 0 0 0 0 R 2 ( sin 2 θ ) 1 .
As before, this is a mere rewriting of the metric. Note that
det ( g a b ) = α 2 γ 1 R ( r ) 4 sin 2 θ 1 ,
and
det ( g a b ) / det ( [ g 0 ] a b ) = α 2 γ 1 R ( r ) 2 .
(We again emphasize that we have not enforced the unimodular condition of [2] since we are now not using quasi-Cartesian coordinates, and the unimodular condition is now neither necessary nor useful).

6.1. Permittivity Tensor

We start by noting
ϵ i j = + γ α | R ( r ) | 1 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 ( sin 2 θ ) 1 i j .
Note that this is conformally invariant (under conformal rescaling of g a b ). The somewhat ugly ( sin 2 θ ) 1 factor can be simplified away by going to a local orthonormal basis in the angular coordinates (adopting an orthonormal dyad, or zweibein) and writing
ϵ i ^ j ^ = + γ α | R ( r ) | 1 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 i ^ j ^ .

6.2. Permeability Tensor

It is most efficient to recall the result we obtained for general space-times:
μ i j = 1 α | det ( g 0 ) | det ( γ p q ) γ i j ,
and to simply unpack the various contributions to obtain the following:
μ i j = + 1 α γ | R ( r ) | γ 0 0 0 R 2 0 0 0 R 2 ( sin 2 θ ) 1 i j .
Alternatively, we can rewrite this as
μ i j = + γ α | R ( r ) | 1 0 0 0 γ 1 R 2 0 0 0 γ 1 R 2 ( sin 2 θ ) 1 i j .
Note that this is conformally invariant (under conformal rescaling of g a b ). Going to a local orthonormal basis for the angular coordinates, adopting an orthonormal (dyad/zweibein), this simplifies to the following:
μ i ^ j ^ = + γ α | R ( r ) | 1 0 0 0 γ 1 R 2 0 0 0 γ 1 R 2 i ^ j ^ .

6.3. Magnetoelectric Tensor

We start by noting
ζ i j = β 2 det ( g ) det ( g 0 ) ε i k r γ k j = β 2 det ( g ) ε ¯ i k r γ k j .
That is,
ζ i j = β 2 α | R | γ 0 0 0 0 0 ( sin θ ) 1 0 sin θ 0 i j .
Note again the conformal invariance under rescaling of g a b . Adopting an orthonormal dyad in the angular directions, we have the relatively simple form:
ζ i ^ j ^ = β 2 α | R | γ 0 0 0 0 0 1 0 1 0 i ^ j ^ .

6.4. Summary for Spherically Symmetric Space-Times

Collecting results, for a generic spherically symmetric space-time we have the following:
ϵ i ^ j ^ = + γ α | R ( r ) | 1 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 0 0 0 ( 1 α 2 β 2 / γ ) γ 1 R 2 i ^ j ^ .
μ i ^ j ^ = + γ α | R ( r ) | 1 0 0 0 γ 1 R 2 0 0 0 γ 1 R 2 i ^ j ^ .
ζ i ^ j ^ = β 2 α | R | γ 0 0 0 0 0 1 0 1 0 i ^ j ^ .
Yet again, getting to this stage has been slightly tedious, but the final results are fully explicit and relatively easy to work with. Several interesting implications immediately follow.

6.5. Implications for Spherically Symmetric Space-Times

  • First, note that ϵ r r = μ r r . This is similar to something that we have seen several times before (the general point being that electromagnetic properties in the direction of the 3-vector [ g 1 ] 0 i are degenerate).
  • Second, note the following:
    ϵ i ^ j ^ μ i ^ j ^ = + γ α | R ( r ) | α 2 β 2 γ 2 R 2 0 0 0 0 1 0 0 0 1 i ^ j ^ .
    Observe that as β i 0 , we find ϵ i j = μ i j and ζ i j = 0 , the standard compatibility condition for static space-times [1,2].
  • Observe that the magnetoelectric tensor always has a zero-eigenvalue eigenvector, now the radial direction, and so det ( ζ i j ) = 0 , as required for materials with a single light cone.
  • Observe that ( ζ 2 ) i j is again proportional to a projection operator:
    ( ζ 2 ) i j = ζ i j ζ j i = 1 4 α 2 β 2 | R | 2 γ 0 0 0 0 1 0 0 0 1 .
    Then,
    tr ( ζ 2 ) = 1 2 α 2 β 2 γ ( R ) 2 .
    This is again a simple scalar invariant describing the strength of the magnetoelectric effect.

7. Discussion

When using metamaterial models to mimic interesting general relativistic space-times, there is often a trade-off between simplicity of presentation and simplicity of calculation. Certainly, easily accessible experimental laboratories are, for all practical purposes, living in flat space-time, and most typically for presentational purposes, one might like to deal with simple Cartesian coordinates, which underpinned the quasi-Cartesian analysis we carried out in a previous article on bespoke analogue space-times [2], wherein the metrics to be mimicked were all cast into unimodular form det ( g a b ) = 1 and raising and lowering indices with the Cartesian background metric was trivial.
In contrast, in the current article, we avoid any quasi-Cartesian assumptions, at the cost of having to deal with and carefully keep track of metric tensor densities (for both the effective metric to be mimicked and the background laboratory metric)—one also has to be careful raising and lowering indices using the background laboratory metric. While intermediate stages of the computations were algebraically tedious, the final answers were both tractable and physically interesting.
The reason for going to this extra effort is essentially a theoretical one—some calculations eventually are more tractable in symmetry adapted coordinate systems, though the initial barrier to setting up the formalism is higher. We hope to return to these issues in future work.

8. Conclusions

This article provides numerous general-purpose calculational tools for electromagnetically mimicking the metrics (line elements) of curved space-times in terms of 3 × 3 electromagnetic permeability, permittivity, and magnetoelectric tensors. Specifically, given any explicit space-time metric (implicitly making a coordinate choice) and choosing any explicit Riemann-flat laboratory metric (background metric), the permeability, permittivity, and magnetoelectric tensors can be uniquely defined and extracted from the formalism. Several special cases are carefully worked out.
Indeed, the current article deals with three explicit themes: Boyer–Lindquist space-times (suitable for dealing with axisymmetric stationary space-times), generic space-times, and finally (as a consistency check) spherically symmetric space-times represented in spherical polar coordinates. In all of these situations, we have been able to present quite specific and easily applicable formulae specifying the constitutive tensors (the permittivity, permeability, and magnetoelectric tensors) required to mimic the given space-time geometry. The long-term goal is to apply these ideas to analogue Hawking radiation [56].

Author Contributions

Conceptualization, S.S. and M.V.; methodology, S.S. and M.V.; software, S.S. and M.V.; validation, S.S. and M.V.; formal analysis, S.S. and M.V.; investigation, S.S. and M.V.; resources, M.V.; writing—original draft preparation, S.S. and M.V.; writing—review and editing, S.S. and M.V.; supervision, M.V.; project administration, M.V.; funding acquisition, M.V. All authors have read and agreed to the published version of the manuscript.

Funding

S.S. acknowledges financial support via OP RDE Project No. CZ.02.2.69/0.0/0.0/18_053/0016976 (international mobility of research) and the technical and administrative staff at the Charles University. M.V. was supported by the Marsden Fund via a grant administered by the Royal Society of New Zealand.

Data Availability Statement

All relevant data are reported in the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Schuster, S.; Visser, M. Effective metrics and a fully covariant description of constitutive tensors in electrodynamics. Phys. Rev. D 2017, 96, 124019. [Google Scholar] [CrossRef]
  2. Schuster, S.; Visser, V. Bespoke analogue space-times: Meta-material mimics. Gen. Relativ. Gravit. 2018, 50, 55. [Google Scholar] [CrossRef]
  3. Schuster, S.; Visser, M. Electromagnetic analogue space-times, analytically and algebraically. Class. Quantum Gravity 2018, 36, 134004. [Google Scholar] [CrossRef]
  4. Gordon, W. Zur Lichtfortpflanzung nach der Relativitätstheorie, (The propagation of light according to the theory of relativity). Ann. Phys. 1923, 377, 421–456. [Google Scholar] [CrossRef]
  5. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields, 3rd ed.; Pergamon Press: Oxford, UK, 1971. [Google Scholar]
  6. Wiltshire, D.L.; Visser, M.; Scott, S.M. (Eds.) The Kerr Space-Time: Rotating Black Holes in General Relativity; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  7. Kerr, R.P. Gravitational field of a spinning mass as an example of algebraically special metrics. Phys. Rev. Lett. 1963, 11, 237–238, Reprinted in ibid. [Google Scholar] [CrossRef]
  8. Kerr, R.P. Gravitational collapse and rotation. 1965; Reprinted in ibid. [Google Scholar]
  9. Visser, M. The Kerr space-time: A brief introduction. [arXiv:0706.0622 [gr-qc]]. Published in ibid.
  10. Doran, C. A New form of the Kerr solution. Phys. Rev. D 2000, 61, 067503. [Google Scholar] [CrossRef]
  11. Hamilton, A.J.S.; Lisle, J.P. The river model of black holes. Am. J. Phys. 2008, 76, 519. [Google Scholar] [CrossRef]
  12. Thompson, R.T.; Fathi, M. Shrinking cloaks in expanding space-times: The role of coordinates and the meaning of transformations in transformation optics. Phys. Rev. A 2015, 92, 013834. [Google Scholar] [CrossRef]
  13. Fathi, M.; Thompson, R.T. Cartographic distortions make dielectric space-time analog models imperfect mimickers. Phys. Rev. D 2016, 93, 124026. [Google Scholar] [CrossRef]
  14. Plebański, J. Electromagnetic waves in gravitational fields. Phys. Rev. 1960, 118, 1396–1408. [Google Scholar] [CrossRef]
  15. Plebanśki, J. Lectures on Nonlinear Electrodynamics; Nordita: Copenhagen, Denmark, 1970. [Google Scholar]
  16. De Felice, F. On the gravitational field acting as an optical medium. Gen. Relativ. Gravit. 1971, 2, 347–357. [Google Scholar] [CrossRef]
  17. Balazs, N.L. Effect of a gravitational field, due to a rotating body, on the plane of polarization of an electromagnetic wave. Phys. Rev. 1958, 110, 236–239. [Google Scholar] [CrossRef]
  18. Anderson, J.L.; Spiegel, E.A. Radiative transfer through a flowing refractive medium. Astrophys. J. 1975, 202, 454–464. [Google Scholar] [CrossRef]
  19. Hehl, F.W.; Obukhov, Y.N. Foundations of Classical Electrodynamics: Charge, Flux and Metric; Birkhäuser: Boston, MA, USA, 2003. [Google Scholar]
  20. Pham, Q.M. Sur les équations de l’electromagné dans la materie. Comptes Rendus Hebd. Seances Acad. Sci. 1956, 242, 465–467. [Google Scholar]
  21. Skrotskii, G.V. The influence of gravitation on the propagation of light. Sov. Phys. Dokl. 1957, 2, 226–229. [Google Scholar]
  22. Thompson, R.T.; Frauendiener, J. Dielectric analog space-times. Phys. Rev. 2010, D82, 124021. [Google Scholar] [CrossRef]
  23. Thompson, R.T.; Cummer, S.A.; Frauendiener, J. Generalized transformation optics of linear materials. J. Opt. 2011, 13, 055105. [Google Scholar] [CrossRef]
  24. Thompson, R.T.; Cummer, S.A.; Frauendiener, J. A completely covariant approach to transformation optics. J. Opt. 2011, 13, 024008. [Google Scholar] [CrossRef]
  25. Thompson, R.T. Covariant electrodynamics in linear media: Optical metric. Phys. Rev. D 2018, 97, 065001. [Google Scholar] [CrossRef]
  26. Xiao, W.; Chen, H. Analogy of the interior Schwarzschild metric from transformation optics. Opt. Express 2023, 31, 11490–11498. [Google Scholar] [CrossRef] [PubMed]
  27. Leonhardt, U.; Piwnicki, P. Optics of nonuniformly moving media. Phys. Rev. A 1999, 60, 4301. [Google Scholar] [CrossRef]
  28. Unruh, W.G. Experimental black hole evaporation. Phys. Rev. Lett. 1981, 46, 1351. [Google Scholar] [CrossRef]
  29. Visser, M. Acoustic propagation in fluids: An unexpected example of Lorentzian geometry. arXiv 1993, arXiv:gr-qc/9311028. [Google Scholar]
  30. Visser, M. Acoustic black holes: Horizons, ergospheres, and Hawking radiation. Class. Quantum Gravity 1998, 15, 1767. [Google Scholar] [CrossRef]
  31. Barceló, C.; Liberati, S.; Visser, M. Analog gravity from Bose-Einstein condensates. Class. Quantum Gravity 2001, 18, 1137. [Google Scholar] [CrossRef]
  32. Barceló, C.; Liberati, S.; Visser, M. Towards the observation of Hawking radiation in Bose-Einstein condensates. Int. J. Mod. Phys. A 2003, 18, 3735. [Google Scholar] [CrossRef]
  33. Barceló, C.; Liberati, S.; Visser, M. Probing semiclassical analog gravity in Bose-Einstein condensates with widely tunable interactions. Phys. Rev. A 2003, 68, 053613. [Google Scholar] [CrossRef]
  34. Visser, M.; Molina-París, C. Acoustic geometry for general relativistic barotropic irrotational fluid flow. New J. Phys. 2010, 12, 095014. [Google Scholar] [CrossRef]
  35. Fischer, U.R.; Visser, M. Riemannian geometry of irrotational vortex acoustics. Phys. Rev. Lett. 2002, 88, 110201. [Google Scholar] [CrossRef] [PubMed]
  36. Visser, M.; Weinfurtner, S. Massive Klein-Gordon equation from a BEC-based analogue spacetime. Phys. Rev. D 2005, 72, 044020. [Google Scholar] [CrossRef]
  37. Barceló, C.; Liberati, S.; Sonego, S.; Visser, M. Quasi-particle creation by analogue black holes. Class. Quantum Gravity 2006, 23, 5341–5366. [Google Scholar] [CrossRef]
  38. Barceló, C.; Liberati, S.; Visser, M. Analogue gravity. Living Rev. Relativ. 2011, 14, 3. [Google Scholar] [CrossRef] [PubMed]
  39. Visser, M.; Barceló, C.; Liberati, S. Analog models of and for gravity. Gen. Relativ. Gravit. 2002, 34, 1719. [Google Scholar] [CrossRef]
  40. Visser, M. Survey of analogue space-times. Lect. Notes Phys. 2013, 870, 31. [Google Scholar] [CrossRef]
  41. Cropp, B.; Liberati, S.; Turcati, R. Vorticity in analog gravity. Class. Quantum Gravity 2016, 33, 125009. [Google Scholar] [CrossRef]
  42. Giacomelli, L.; Liberati, S. Rotating black hole solutions in relativistic analogue gravity. Phys. Rev. D 2017, 96, 064014. [Google Scholar] [CrossRef]
  43. Liberati, S.; Schuster, S.; Tricella, G.; Visser, M. Vorticity in analogue space-times. Phys. Rev. D 2019, 99, 044025. [Google Scholar] [CrossRef]
  44. Jacobson, T.; Kang, G. Conformal invariance of black hole temperature. Class. Quantum Gravity 1993, 10, L201. [Google Scholar] [CrossRef]
  45. Rubilar, G.F. Linear pre-metric electrodynamics and deduction of the light cone. Ann. Phys. 2002, 11, 717–782. [Google Scholar] [CrossRef]
  46. Balakin, A.B.; Zimdahl, W. Optical metrics and birefringence of anisotropic media. Gen. Relativ. Gravit. 2005, 37, 1731–1751. [Google Scholar] [CrossRef]
  47. Wald, R.M. General Relativity; University of Chicago Press: Chicago, IL, USA, 1984. [Google Scholar] [CrossRef]
  48. Misner, C.W.; Thorne, K.S.; Wheeler, J.A. Gravitation; Freeman: San Francisco, CA, USA, 1973. [Google Scholar]
  49. Thirring, H.; Lense, J. Über den Einfluss der Eigenrotation der Zentralkörper auf die Bewegung der Planeten und Monde nach der Einsteinschen Gravitationstheorie. Phys. Z. Leipz. Jg 1918, 19, 156–163. [Google Scholar]
  50. Mashoon, B.; Hehl, F.W.; Theiss, D.S. (Translators) On the influence of the proper rotations of central bodies on the motions of planets and moons in Einstein’s theory of gravity. (Translation of Lense–Thirring 1918, ibid.). Gen. Relativ. Gravit. 1984, 16, 727–741. [Google Scholar]
  51. Pfister, H. On the history of the so-called Lense–Thirring effect. Gen. Relativ. Gravit. 2007, 39, 1735–1748. [Google Scholar] [CrossRef]
  52. Boersma, S.; Dray, T. Slicing, threading & parametric manifolds. Gen. Relativ. Gravit. 1995, 27, 319. [Google Scholar] [CrossRef]
  53. Bejancu, A.; Calin, C. On the (1+3) threading of space-time with respect to an arbitrary timelike vector field. Eur. Phys. J. C 2015, 75, 159. [Google Scholar] [CrossRef]
  54. Gharechahi, R.; Nouri-Zonoz, M.; Tavanfar, A. A tale of two velocities: Threading vs slicing. Int. J. Geom. Meth. Mod. Phys. 2018, 15, 1850047. [Google Scholar] [CrossRef]
  55. Rosquist, K. A Moving medium simulation of Schwarzschild black hole optics. Gen. Relativ. Gravit. 2004, 36, 1977. [Google Scholar] [CrossRef]
  56. Reznik, B. Origin of the thermal radiation in a solid state analog of a black hole. Phys. Rev. D 2000, 62, 044044. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Schuster, S.; Visser, M. Boyer–Lindquist Space-Times and Beyond: Metamaterial Analogues for Arbitrary Space-Times. Universe 2024, 10, 159. https://doi.org/10.3390/universe10040159

AMA Style

Schuster S, Visser M. Boyer–Lindquist Space-Times and Beyond: Metamaterial Analogues for Arbitrary Space-Times. Universe. 2024; 10(4):159. https://doi.org/10.3390/universe10040159

Chicago/Turabian Style

Schuster, Sebastian, and Matt Visser. 2024. "Boyer–Lindquist Space-Times and Beyond: Metamaterial Analogues for Arbitrary Space-Times" Universe 10, no. 4: 159. https://doi.org/10.3390/universe10040159

APA Style

Schuster, S., & Visser, M. (2024). Boyer–Lindquist Space-Times and Beyond: Metamaterial Analogues for Arbitrary Space-Times. Universe, 10(4), 159. https://doi.org/10.3390/universe10040159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop