## Towards a piecewise-homogeneous metamaterial model of the collision of two linearly polarized gravitational plane waves

Tom G. Mackay ∗

School of Mathematics and Maxwell Institute for Mathematical Sciences University of Edinburgh, Edinburgh EH9 3JZ, UK

and

NanoMM - Nanoengineered Metamaterials Group Department of Engineering Science and Mechanics Pennsylvania State University, University Park, PA 16802-6812, USA

## Akhlesh Lakhtakia †

NanoMM - Nanoengineered Metamaterials Group Department of Engineering Science and Mechanics Pennsylvania State University, University Park, PA 16802-6812, USA

## Abstract

We considered the experimental realization of a Tamm medium that is optically equivalent to the collision of two linearly polarized gravitational plane waves as a piecewise-homogeneous metamaterial. Our formulation was based on the homogenization of remarkably simple arrangements of oriented ellipsoidal nanoparticles of isotropic dielectric-magnetic mediums. The inverse Bruggeman homogenization was used to estimate the constitutive parameters, volume fractions, and shape parameters for the component mediums. The presented formulation is appropriate for the regions of spacetime where the two gravitational plane waves interact, excluding the immediate vicinity of the nonsingular Killing-Cauchy horizon at the focusing point of the two plane waves.

∗ E-mail: T.Mackay@ed.ac.uk.

† E-mail: akhlesh@psu.edu

Keywords: colliding gravitational plane waves; homogenization; Killing-Cauchy horizon; metamaterial; inverse Bruggeman formalism

## 1 Introduction

Metamaterials are engineered materials which, through judicious design, can facilitate the realization of exotic phenomenons such as negative refraction and cloaking [1]. Metamaterials also present opportunities to study general-relativistic scenarios. This is because a formal analogy exists between light propagation in empty curved spacetime and light propagation in a certain nonhomogeneous anisotropic or bianisotropic medium, called a Tamm medium [2, 3, 4]. For examples, theoretical descriptions of metamaterial analogs of black holes [5], de Sitter spacetime [6, 7] including Schwarzschild-(anti-)de Sitter spacetime [8], strings [9] including cosmic strings [10], and wormholes [11], have been investigated. Such analogs may offer insights into spacetime scenarios which cannot be practicably studied by direct methods.

In theory, analogs of curved spacetime may be constructed from metamaterials. However, achieving this in practice presents a major challenge to experimentalists. Concrete experimental proposals are conspicuously absent from the literature, but there are some exceptions. Lu et al. [12] presented a detailed metamaterial representation of a two-dimensional black hole. Their metamaterial took the form of a homogenized composite material (HCM) comprising relatively simple component mediums. We established a metamaterial formulation for the Tamm medium representing Schwarzschild-(anti-)de Sitter spacetime [8]. In our approach the metamaterial was also a HCM, arising from the homogenization of isotropic dielectric and isotropic magnetic component mediums which are distributed randomly as oriented spheroidal particles. In this paper, we apply the same approach to the case of the Tamm medium corresponding to the collision of two linearly polarized gravitational plane waves [13].

The collision process is a complex one from the perspective of both theoretical and numerical analyses: its intrinsic nonlinearity gives rise to either a spacetime singularity or a nonsingular Killing-Cauchy horizon at the focusing point of the two plane waves [13]. As such, a convenient experimental analog could be particularly useful, as was succinctly pointed by Bini et al. [14]. Colliding gravitational plane waves may be used to study Cauchy horizons in classical general relativity [15, 16], as well as string behaviour in gravitational fields [17, 18]. These plane waves may also assist in the study of black hole collisions [19] and travelling waves on cosmic strings [20].

Motivation for the present paper was chiefly provided by a recent study [14] in which the Tamm medium analogy was exploited to to highlight some interesting optical properties of the spacetime of two colliding gravitational plane waves.

The plan of this paper is as follows: Section 2 describes the Tamm medium for the region of spacetime in which the two plane waves interact to form the nonsingular Killing-Cauchy horizon. Section 3 succinctly describes the inverse Bruggeman homogenization formalism used for the piecewise-homogeneous implementation of the Tamm medium, followed by illustrative numerical results in Sec. 4. The paper concludes with a discussion in Sec. 5. In the units adopted, the Newtonian constant and the speed of light in vacuum are both equal to unity.

## 2 Tamm medium for colliding-gravitational-wave spacetime

The collision of two oppositely-directed gravitational planes waves may be represented as an exact solution of the Einstein field equations [13]. Following the collision, focussing effects give rise to either a nonsingular Killing-Cauchy horizon or a spacetime singularity. Conventionally, the corresponding spacetime, prior to the creation of the nonsingular Killing-Cauchy horizon or the spacetime singularity, is partitioned into four regions: Region I, where the two plane waves interact; Regions II and III, each corresponding to a single plane wave before the interaction; and Region IV, which is a flat spacetime region representing the initial state before the passage of the two plane waves. We focus on Region I bounded by

$$\begin{array} {$$

 for propagation along the z axis. Here the instant of collision occurs at t = 0 while at t = π/ 2 either a nonsingular Killing-Cauchy horizon or a spacetime singularity is created.

Provided that the plane waves are linearly polarized and they propagate in opposite directions along the z axis, the corresponding line element may be expressed as [21, 22, 23]

where the scalar function

$$d s ^ { 2 } = F _ { + } ^ { 2 } ( t ) \left ($$

$$F _ { \pm } ( t ) = 1 \pm \sigma \sin \left ( t$$

Herein σ = ± 1, with σ = +1 corresponding to the nonsingular Killing-Cauchy horizon solution at t = π/ 2 whereas σ = -1 corresponds to the spacetime singularity solution at t = π/ 2. In the following we consider the σ = +1 option; the prescription of a metamaterial analog for the σ = -1 case follows in a similar vein to the σ = +1 case but the constitutive-parameter regimes for σ = -1 would be more challenging to realize in practice than those for σ = +1.

The optical response of vacuum in the curved spacetime represented by the line element (2) is formally equivalent to the optical response of a spatially and temporally nonhomogeneous medium in flat spacetime

[2, 3, 4], known as the Tamm medium. As recently shown by Bini et al. [14], the Tamm medium here is an orthorhombic dielectric-magnetic medium [25], which may be characterized by the electromagnetic constitutive relations

 wherein SI units are implemented with /epsilon1 0 = 8 854 . × 10 -12 F m -1 and µ 0 = 4 π × 10 -12 H m -1 , and

$$\underline { D } ( x, y, z, t ) & = \epsilon _$$

$$\underline { \gamma } ( z, t ) = \text{diag} \left$$

with

In Fig. 1, the constitutive scalars γ x,y,z are plotted versus z ∈ ( -t, t ) for t = π/ 6 (green dashed curves), π/ 3 (red solid curves), and 4 π/ 9 (blue broken dashed curves). The region (4 π/ 9) &lt; t ≤ ( π/ 2), which includes the nonsingular Killing-Cauchy horizon solution at t = π/ 2, is avoided in Fig. 1 (and later in Fig. 2). We have that γ x &gt; 1 while γ y,z &lt; 1 for t &gt; 0. The constitutive scalars increasingly deviate from unity as t increases from zero. In particular, γ x very rapidly increases in magnitude as t increases, while the magnitudes of γ y and γ z decrease very rapidly as t increases. Notice that γ y = γ z in the plane z = 0; consequently, the Tamm medium is uniaxial dielectric-magnetic in that plane.

$$\gamma ^ { x } ( z, t ) & = \cos ( z ) \sqrt { \frac { F _ { 3 } ^ { 3 } ( t ) } { F _ { - } ( t ) } } \\ \gamma ^ { y } ( z, t ) & = \frac { 1 } { \cos ( z ) } \sqrt { \frac { F _ { - } ( t ) } { F _ { + } ^ { 3 } ( t ) } } \ \Big. \\ \gamma ^ { z } ( z, t ) & = \cos ( z ) \sqrt { \frac { F _ { - } ( t ) } { F _ { + } ^ { 3 } ( t ) } } \ \Big. \\ \intertext { e scalars } \gamma ^ { x, y, z } \, \text{are plotted versus } z \in ( - t, t ) \text{ for } t = \pi / 6 \text{ (green dashed curves),} \\ \i \pi / 9 \text{ (blue broken dashed curves).  The region } ( 4 \pi / 9 ) < t \leq ( \pi / 2 ), \text{ which includes}$$

A different perspective on the dependency of the constitutive parameters upon z and t is provided in Fig. 2, where the constitutive scalars γ x,y,z are plotted versus t ∈ ( z, 4 π/ 9) for z = 0 (green dashed curves), π/ 6 (red solid curves), and π/ 3 (blue broken dashed curves). We see that γ x becomes unbounded as t approaches π/ 2, whereas both γ y and γ z vanish as t approaches π/ 2, for all values of z .

## 3 Inverse Bruggeman formalism

The Tamm medium characterized by the constitutive parameters (6) is an orthorhombic dielectric-magnetic medium with identical relative permittivity and relative permeability dyadics. Such a medium may be conceptualized as a homogenized composite material (HCM), using the well-established Bruggeman homogenization formalism [26].

Figure 1: The parameters γ x,y,z plotted versus z ∈ ( -t, t ) for t = π/ 6 (green dashed curves), π/ 3 (red solid curves), and 4 π/ 9 (blue broken dashed curves).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0601v1-figure-077f7ec6-d3aa-4d23-8f09-4d33cf0a3fab.jpg)

We consider the homogenization of two component mediums, identified by the labels a and b . Both component mediums are isotropic dielectric-magnetic mediums with relative permittivities /epsilon1 a and /epsilon1 b , and

Figure 2: The parameters γ x,y,z plotted versus t ∈ ( z, 4 π/ 9) for z = 0 (green dashed curves), π/ 6 (red solid curves), and π/ 3 (blue broken dashed curves).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0601v1-figure-7cfb265c-4250-4428-90c6-ccdd2d9e923c.jpg)

relative permeabilities µ a and µ b . In consonance with the relative permittivity and relative permeability dyadics of the Tamm medium being identical, we consider the component mediums characterized by /epsilon1 a =

µ a = γ a and /epsilon1 b = µ b = γ b . Component mediums with such characteristics do not generally occur in nature but the sought after component mediums could be envisaged as HCMs themselves, arising from the homogenization of isotropic dielectric and isotropic magnetic mediums [8].

Both component mediums are distributed randomly, and their respective volume fractions are f a = 1 -f b and f b ∈ (0 , 1). Each component medium is present as an assembly of electrically small ellipsoidal particles. For example, at optical wavelengths, component particles with linear dimensions less than approximately 40-70 nanometers would be required. In keeping with the orthorhombic symmetry of the Tamm medium, the major and minor axes of the component ellipsoidal particles are aligned with the coordinate axes. All ellipsoidal particles are assumed to have the same shape. Thus, the vector r s = ρ /lscript U r · ˆ, ( /lscript = a, b ), prescribes the surface of each ellipsoidal particle relative to its centroid, with ˆ = (sin r θ cos φ, sin θ sin φ, cos φ ) as the radial unit vector, and θ ∈ [0 , π ] and φ ∈ [0 , 2 π ) being the polar and azimuthal angles, respectively. A linear measure of the size of the component ellipsoidal particles is provided by ρ /lscript . The shape of the component particles is captured by the positive-definite dyadic U = diag ( U ,U ,U x y z ).

The Bruggeman formalism requires the polarizability density dyadics [26]

$$\underline { a } _ { \ell } = \left ( \gamma _ { \ell } \underline { I } - \underline { \gamma } \right ) \cdot \left [ \underline { I } + \underline { D } _ { \ell } \cdot \left ( \gamma _ { \ell } \underline { I } - \underline { \gamma } \right ) \right ] ^ { - 1 }, \quad ( \ell = a, b ) \,. \quad \, \intertext { s u d e d i n d } \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \, \ s u d e d i n d \end{cases}$$

The depolarization dyadics D /lscript herein are given by the double integrals [27, 28]

$$\underline { D } _ { \ell } = \frac { 1 } { 4 \pi } \int _ { \phi } ^ { 2 \pi } d \phi \int _ { \theta } ^ { \pi } d \theta \, \sin \theta \, \frac { \left ( \underline { U } ^ { - 1 } \cdot \hat { \underline { r } } \right ) \left ( \underline { U } ^ { - 1 } \cdot \hat { \underline { r } } \right ) } { \left ( \underline { U } ^ { - 1 } \cdot \hat { \underline { r } } \right ) \cdot \underline { \gamma } \cdot \left ( \underline { U } ^ { - 1 } \cdot \hat { \underline { r } } \right ) }, \\ \left ( \ell = a, b \right ) \,,$$

which may be conveniently represented in terms of incomplete elliptic integrals [29]. The constitutive parameters of the HCM are related to the constitutive and shape parameters, and volume fractions, of the component mediums by the Bruggeman equation [25, 26]

$$f _ { a } \, \underline { a } _ { a } + f _ { b } \, \underline { a } _ { b } = \underline { 0 },$$

which is equivalent to three independent scalar equations, which are coupled via the constitutive parameters for the HCM and the shape parameters for the component mediums.

Without loss of generality, we fix U x = 1. Then we recast Eq. (9) to solve the following problem: With γ z, t ( ) known, determine γ a ( z, t ), f a ( z, t ) and U y ( z, t ) in terms of γ b ( z, t ) and U z ( z, t ) for any combination of z and t in Region 1. The inverse Bruggeman formalism is then conveniently implemented for any specified z and t by exploiting direct numerical methods [30].

## 4 Numerical results

Let us now provide representative numerical illustrations of component mediums which could be homogenized to realize the Tamm medium specified by Eqs. (5) and (6), based on the inverse Bruggeman formalism described in Sec. 3. The Tamm medium for the Region I spacetime was considered for two cases: the first concerns fixed t while the second concerns fixed z . For both cases the immediate vicinity of the KillingCauchy horizon at t = π/ 2 was excluded, in order to avoid extreme constitutive parameter regimes for the Tamm medium.

We begin by considering the Tamm medium at fixed t . Estimates of γ a , f a and U y are plotted against z ∈ (0 , t ) for t = π/ 6, π/ 4 and π/ 3 in Fig. 3. The shape parameter U z = 0 004 for these calculations, while . γ b = 0 3 for . t = π/ 6, γ b = 0 14 for . t = π/ 4, and γ b = 0 03 for . t = π/ 3. The values of U y rise sharply as z increases, for all values of t considered. In the limit z → 0, for all values of t the values of the shape parameter U y converge on 0 004, which is the value of . U z . This is a reflection of the fact that the Tamm medium is a uniaxial dielectric-magnetic medium in the z = 0 plane. The variation in γ a as z varies depends upon the value of t : at larger values of t , γ a is more sensitive to changes in z . The volume fraction f a varies relatively little as z varies.

Next, let us consider the Tamm medium at fixed z . In Fig. 4 estimates of γ a , f a and U y are plotted against t ∈ ( z, π/ 3) for z = 0, π/ 6 and π/ 4. Here γ b = 0 05 for . z = 0, γ b = 0 08 for . z = π/ 6, and γ b = 0 05 . for t = π/ 4. As in Fig. 3, the shape parameter U z = 0 004. . When z = 0, the shape parameter U y again remains constant at U z . Aside from this exception, the parameters γ a , f a and U y all vary markedly as t increases, for all values of z considered.

## 5 Discussion

Graded dispersals of electrically small ellipsoids of two isotropic dielectric-magnetic component mediums may be used to realize the Tamm medium representing the collision of two linearly polarized gravitational plane waves. By varying the volume fraction and the relative permittivity (or permeability) of each component mediums, along with the common ellipsoidal shape, Region I in spacetime can be accessed except the immediate vicinity of the Killing-Cauchy horizon at t = π/ 2. Thus, the inverse homogenization formulation paves the way for an experimental analog for the collision of two linearly polarized gravitational plane waves. For example, the experimental analog could be used to study either the temporal evolution of the colliding plane waves at a fixed location, or the spatial portrait of the colliding plane waves at a fixed time. The nonhomogeneous nature of the underlying spacetime may be accommodated by stitching together

Figure 3: The parameters γ a , f a , and U y plotted versus z ∈ ( -t, t ) for t = π/ 6 (green dashed curves), π/ 4 (red solid curves), and π/ 3 (blue broken dashed curves).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0601v1-figure-c11ebd2b-542f-4c24-8408-d7b9a10c56ad.jpg)

adjacent spacetime regions that are sufficiently small that they could be approximated by homogeneous

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0601v1-figure-53438274-b02a-4433-92c8-3458258100e3.jpg)

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0

0.009

0.008

0.007

0.006

0.005

0.004

a

f

U y

0

10

10

20

t

30

/LParen1

20

180

/Slash1

Π

30

40

40

50

50

60

60

t

/LParen1

180

/Slash1

Π

/RParen1

Figure 4: The parameters γ a , f a , and U y plotted versus t ∈ ( z, π/ 3) for z = 0 (green dashed curves), π/ 6 (red solid curves), and π/ 4 (blue broken dashed curves).

Tamm mediums. Experimental realization for optical experimentation would then involve nanocomposite

/RParen1

## materials [31].

The presented metamaterial model is essentially wavelength independent. Frequencies higher or lower than optical frequencies could be used in order to access regimes where the values of relative permittivity and/or relative permeability required for the component mediums may be more readily attained. However, if higher frequencies were used, then the component medium particles would have to be smaller in order for the concept of homogenization to remain valid [25].

The particular inverse homogenization formulation chosen for presentation in Sec. 3 was based on two isotropic dielectric-magnetic component mediums, each being an assembly of oriented ellipsoidal particles. However, alternative formulations may be envisaged. For example, the homogenization of four component mediums could also yield an HCM representing the Tamm medium. In this case, two isotropic dielectric mediums and two isotropic magnetic mediums, with each component medium comprising oriented spheroidal particles, could be used; or instead two uniaxial dielectric mediums and two uniaxial magnetic mediums, with each component medium comprising spherical particles, could be used [32].

Let us comment on the restriction of the inverse homogenization approach to Region I in spacetime away from the immediate vicinity of the Killing-Cauchy horizon at t = π/ 2. As the point t = π/ 2 is approached, the Tamm medium becomes extremely anisotropic (as may be appreciated from Figs. 1 and 2). Accordingly, component mediums specified by extreme values of relative permittivity (and permeability)i.e., both extremely high values and values exceedingly close to zero-would be needed to implement the corresponding HCM. At first glance, the realization of such parameter values would pose a major challenge to experimentalists. However, the relative permittivities (and permeabilities) represented in Figs. 3-4 are not beyond the grasp of contemporary materials researchers, and these correspond to a wide area of the Region I spacetime. Indeed, considerable progress been made within recent years as regards metamaterials with extremely high relative permittivities and permeabilities [33] and with relative permittivities and relative permeabilities very close to zero [34]. Furthermore, extremely high degrees of electric and magnetic anisotropy can be achieved via the homogenization of highly elongated component particles [35, 36].

Lastly, in principle, the metamaterial model described here could be extended, using the same techniques as presented in Secs. 2 and 3, to rather more complicated scenarios; i.e., ones even less well suited to analytical and numerical study than the collision of gravitational plane waves described herein. For example, by combining with results presented previously [8], the collision of gravitational plane waves in a Schwarzschild(anti-)de Sitter spacetime background could be contemplated.

## References

- [1] T. J. Cui, D. R. Smith, and R. Liu (Eds.), Metamaterials: Theory, Design, and Applications , Springer, New York, 2010.
- [2] G. V. Skrotskii, 'The influence of gravitation on the propagation of light,' Soviet Phys.-Dokl. , vol. 2, pp. 226-229, 1957.
- [3] J. Pl´banski, 'Electromagnetic waves in gravitational fields,' e Phys. Rev. , vol. 118, pp. 1396-1408, 1960.
- [4] W. Schleich and M. O. Scully, 'General relativity and modern optics,' New Trends in Atomic Physics , G. Grynberg and R. Stora (Eds.), Elsevier, Amsterdam, pp. 995-1124, 1984.
- [5] I. I. Smolyaninov, 'Surface plasmon toy model of a rotating black hole,' New J. Phys. , vol. 5, 147, 2003.
- [6] M. Li, R.-X. Miao, and Y. Pang, 'Casimir energy, holographic dark energy and electromagnetic metamaterial mimicking de Sitter,' Phys. Lett. B , vol. 689, pp. 55-59, 2010.
- [7] M. Li, R.-X. Miao, and Y. Pang, 'More studies on metamaterials mimicking de Sitter space,' Opt. Express , vol. 18, pp. 9026-9033, 2010.
- [8] T. G. Mackay and A. Lakhtakia, 'Towards a realization of Schwarzschild-(anti-)de Sitter spacetime as a particulate metamaterial,' Phys. Rev. B , vol. 83, 195424, 2011.
- [9] R.-X. Miao, R. Zheng, and M. Li, 'Metamaterials mimicking dynamic spacetime, D-brane and noncommutativity in string theory,' Phys. Lett. B , vol. 696, pp. 550-555, 2011.
- [10] T. G. Mackay and A. Lakhtakia, 'Towards a metamaterial simulation of a spinning cosmic string,' Phys. Lett. A , vol. 374, pp. 2305-2308, 2010.
- [11] A. Greenleaf, Y. Kurylev, M. Lassas, and G. Uhlmann, 'Cloaking devices, electromagnetic wormholes, and transformation optics,' SIAM Rev. , vol. 51, pp. 3-33, 2009.
- [12] W. Lu, J. Jin, Z. Lin, and H. Chen, 'A simple design of an artificial electromagnetic black hole,' J. Appl. Phys. , vol. 108, 064517, 2010.
- [13] J. B. Griffiths, Colliding Plane Waves in General Relativity , Clarendon, Oxford, 1991.
- [14] D. Bini, A. Geralico, and M. Haney, 'Refraction index analysis of light propagation in a colliding gravitational wave spacetime,' Gen. Relativ. Gravit. , vol. 46, 1644, 2014.
- [15] A. Ori, 'Structure of the singularity inside a realistic rotating black hole,' Phys. Rev. Lett. , vol. 68, pp. 2117-2120, 1992.

- [16] U. Yurtsever, 'Comments on the instability of black-hole inner horizons,' Class. Quantum Grav. , vol. 10, L17, 1993.
- [17] H. J. de Vega and N. S´nchez, 'Strings falling into spacetime singularities,' a Phys. Rev. D , vol. 45, 2783-2793, 1992.
- [18] O. Jofre and C. N´˜ez, 'Strings in plane wave backgrounds reexamined,' un Phys. Rev. D , vol. 50, pp.52325240, 1992.
- [19] V. Ferrari, P. Pendenza, and G. Veneziano, 'Beam-like gravitational waves and their geodesies,' Gen. Relativ. Gravit. , vol. 20, pp. 1185-1191, 1988.
- [20] D. Garfinkle and T. Vachaspati, 'Cosmic-string traveling waves,' Phys. Rev. D , vol. 42, pp. 1960-1963, 1990.
- [21] V. Ferrari and J. Iba˜ez, 'A new exact solution for colliding gravitational plane waves,' n Gen. Relativ. Gravit. , vol. 19, pp. 383-404, 1987.
- [22] V. Ferrari and J. Iba˜ez, 'On the collision of gravitational plane waves: n a class of soliton solutions,' Gen. Relativ. Gravit. , vol. 19, pp. 405-425, 1987.
- [23] V. Ferrari and J. Iba˜ez, 'Type-D solutions describing the collision of plane-fronted gravitational waves,' n Proc. R. Soc. Lond. A , vol. 417, pp. 417-431, 1988.
- [24] T. G. Mackay, A. Lakhtakia, and S. Setiawan, 'Gravitation and electromagnetic wave propagation with negative phase velocity,' New J. Phys. , vol. 7, 75, 2005.
- [25] T. G. Mackay and A. Lakhtakia, Electromagnetic Anisotropy and Bianisotropy: A Field Guide , Word Scientific, Singapore, 2010.
- [26] W. S. Weiglhofer, A. Lakhtakia, and B. Michel, 'Maxwell Garnett and Bruggeman formalisms for a particulate composite with bianisotropic host medium,' Microw. Opt. Technol. Lett. , vol. 15, pp. 263266, 1997; corrections: vol. 22, p. 221, 1999.
- [27] B. Michel, 'A Fourier space approach to the pointwise singularity of an anisotropic dielectric medium,' Int. J. Appl. Electromagn. Mech. , vol. 8, pp. 219-227, 1997.
- [28] B. Michel and W. S. Weiglhofer, 'Pointwise singularity of dyadic Green function in a general bianisotropic medium,' Arch. Elektr. ¨ bertrag. U , vol. 51, pp. 219-223, 1997; corrections: vol. 52, p. 310, 1998.

- [29] W. S. Weiglhofer, 'Electromagnetic depolarization dyadics and elliptic integrals,' J. Phys. A: Math. Gen. , vol. 31, pp. 7191-7196, 1998.
- [30] T. G. Mackay and A. Lakhtakia, 'Determination of constitutive and morphological parameters of columnar thin films by inverse homogenization,' J. Nanophotonics , vol. 4, 041535, 2010.
- [31] T. G. Mackay, 'On extended homogenization formalisms for nanocomposites,' J. Nanophotonics , vol. 2, 021850, 2008.
- [32] T. G. Mackay and W. S. Weiglhofer, 'Homogenization of biaxial composite materials: nondissipative dielectric properties,' Electromagnetics , vol. 21, pp. 15-26, 2001.
- [33] J. Shin, J.-T. Shen, and S. Fan, 'Three-dimensional metamaterials with an ultrahigh effective refractive index over a broad bandwidth,' Phys. Rev. Lett. , vol. 102, 093903, 2009.
- [34] M. Navarro-Cia, M. Beruete, I. Campillo, and M. Sorolla, 'Enhanced lens by /epsilon1 and µ near-zero metamaterial boosted by extraordinary optical transmission,' Phys Rev. B , vol. 83, 115112, 2011.
- [35] A. H. Gevorgyan, 'Magneto-optics of a thin film layer with helical structure and enormous anisotropy,' Mol. Cryst. Liq. Cryst. , vol. 382, pp. 1-19, 2002.
- [36] T. G. Mackay, 'Giant dielectric anisotropy via homogenization,' http : // arxiv org . / abs / 1405 1198 .
Tom G. Mackay Akhlesh Lakhtakia
Towards a piecewise-homogeneous metamaterial model of the collision of two linearly polarized gravitational plane waves
We considered the experimental realization of a Tamm medium that is optically
equivalent to the collision of two linearly polarized gravitational plane waves
as a piecewise homogeneous metamaterial. Our formulation was based on the
homogenization of remarkably simple arrangements of oriented ellipsoidal
nanoparticles of isotropic dielectric-magnetic mediums. The inverse Bruggeman
homogenization was used to estimate the constitutive parameters, volume
fractions, and shape parameters for the component mediums. The presented
formulation is appropriate for the regions of spacetime where the two
gravitational plane waves interact, excluding the immediate vicinity of the
nonsingular Killing-Cauchy horizon at the focusing point of the two plane
waves.