## Optical Kerr effect in vacuum

Scott Robertson 1

1 LAL, Univ. Paris-Sud, CNRS/IN2P3, Universit´ e Paris-Saclay, Orsay, France (Dated: December 17, 2019)

From an effective field theory of electromagnetism in vacuum including all lowest-order nonlinear terms consistent with Lorentz invariance and locality of photon/photon interactions, we derive an effective medium description of strong background fields as regards their influence on a weak probe. We mainly consider as background a pump beam with well-defined wave vector and polarization. This leads us to define a nonlinear index of vacuum which, in the Euler-Heisenberg model derived from QED, has an optimal value of 1 555 . × 10 -33 cm 2 / W for a linearly polarized pump as seen by a counter-propagating, orthogonally polarized probe. We further generalize the model to include coupling to an axion field. In the limit where the axion mass is much smaller than the typical photon energy, this yields dispersive corrections, and the axionic signature is found to be greatly enhanced for a circularly polarized pump as compared to a linearly polarized one. The formalism here presented points to a simplification of the DeLLight experiment [X. Sarazin et al. , Eur. Phys. J. D 70 , 13 (2016)] aiming to measure the deflection of a probe by a tightly focused laser pulse.

## I. INTRODUCTION

In media, the dependence of optical properties on the electric or magnetic field has been known since the time of Faraday in the mid 19th-century [1], though it has gained particular prominence in the last sixty years with the availability of high-intensity lasers and the subsequent development of nonlinear optics [2-6]. Such field-dependent behavior arises from the nontrivial response of bound charges and currents within the medium, leading to a range of nonlinear effects which have been extensively studied in the literature. Since at optical frequencies the intensity associated with a single photon is very low, reaching the nonlinear regime requires an intense field comprising so many photons that it can be treated classically. From the point of view of a weak probe, a system comprising a strong background field and a medium can be treated as a single 'dressed' medium, with the background field contributing to the total refractive index [5, 6]. In typical dielectric media whose molecules possess inversion symmetry, the refractive index change is proportional to the square of the electric field, a phenomenon usually referred to as the Kerr effect after its discoverer John Kerr [7]. When the index change is engendered by intense light (rather than static fields), this optical Kerr effect allows the assignment of a nonlinear index n 2 to the medium [5, 6] such that the total refractive index includes a term proportional to the intensity I of the wave:

$$n ( I ) = n _ { 0 } + n _ { 2 } \, I \,,$$

where n 0 is the 'bare' index of the medium in the absence of strong fields. Though n 2 depends on the precise configuration (such as the relative polarization between the probe and the background), values typically range from 10 -16 to 10 -14 cm 2 / W (see, e.g., Table 4.1.2 of [5]).

In vacuum, by contrast, classical electrodynamics is a linear theory, where n 2 is exactly zero. By analogy with the situation in media, we may ask whether this apparent linearity is only a low-field approximation, i.e., whether the field equations might become nonlinear when the fields are strong enough. Indeed, the Standard Model already answers in the affirmative; in particular, quantum electrodynamics (QED) allows photon/photon scattering mediated by virtual electron/positron pairs, which play a role analogous to that of bound charges in media. In the long-wavelength limit, this yields a nonlinear effective field theory for the electromagnetic field, the Lagrangian of which was first derived by Euler, Kockel and Heisenberg [8-10]. In principle, however, this is but one way in which nonlinearities could be generated: there may well be as-yet-unidentified particles (such as axions [11-14]) coupling to photons and thereby contributing to the effective nonlinear response; alternatively, there may be higher-order corrections to the classical electromagnetic sector of the Lagrangian (such as proposed by Born and Infeld [15]).

It is thus of interest to experimentally probe nonlinear electrodynamics (NLED) in vacuum, in order to test our current predictions and potentially rule out alternative models [16, 17]. While there have been several experiments in the high-energy photon regime which also tend to involve charged particles of some kind (a recent example is provided by the heavy ion collisions observed at the LHC [18]; see Ref. [16] for many others), the direct elastic scattering of real photons has not yet been observed, and the low-energy photon regime - where, as in media, the strong background field is classical and nonlinearities can be treated as field-dependent contributions to the total refractive index - remains relatively unexplored. The most sensitive tests of low-energy NLED to date are those of the BMV [19, 20] and PVLAS [21, 22] experiments, which aim to detect the birefringence induced by a strong magnetic field perpendicular to the direction of the probe wave. These have not yet reached the sensitivity required to test effects on the order of

those predicted by QED. An alternative, complementary approach is to directly exploit the intensity-dependence of the refractive index by using strong background fields to deflect the trajectory of a probe wave. This was attempted by Jones in 1960 [23, 24] using a static magnetic field, while the recently proposed DeLLight experiment [25] aims to observe such a deflection using the much greater intensities within a tightly focused laser pulse.

Inspired by such proposals, we develop in this paper a theoretical framework for describing the lowest-order nonlinear interactions between electromagnetic waves in vacuum. The novelty lies in the generality of the approach, which covers two aspects. First, rather than restricting ourselves to the Euler-Heisenberg model describing the nonlinearities induced by QED, we consider instead a generalized model (in the spirit of Pleba´ski [26] and Boillat [27]) which n is consistent with Lorentz invariance and locality of effective photon/photon interactions. As far as we are aware, this generalized Lagrangian was previously adopted only in [28] (for the particular case of static background fields), though the approach used in [29] is reminiscent of it. Moreover, being mindful of proposals for the detection of axions [30-34], we also allow a relaxation of the assumption of local photon/photon interactions through coupling to an axion-like field. This allows an explicit demonstration that, while the inclusion of the axion coupling would lead to a straightforward renormalization of the model parameters if the axion mass were large enough, the breakdown of locality of photon/photon interactions when this is not the case yields dispersive corrections to the behavior of a probe wave. It also, in the case of an elliptically polarized background, induces some elliptical birefringence that is completely absent in the purely local theory.

The second aspect in which this paper generalizes the standard treatment of NLED concerns its explicitization of the effective medium engendered by an arbitrary configuration of strong background fields. That is, assuming an explicit decomposition of the total fields into those of a strong background and a weak probe, the field equations can be linearized in those of the latter, in which case the 'dressed' vacuum comprising the vacuum and strong background fields together behaves as an optical medium in its own right. To this effective medium is assigned a set of well-defined susceptibility tensors. This description is more intuitive than the manifestly covariant formalism commonly used in the literature (as in, e.g., [35]). The general expression for the susceptibility tensors encompasses any background field configuration, describing both static and propagating fields; in the latter case, elliptical polarization is included in a relatively straightforward manner. The behavior of a probe wave is completely determined by the optical properties encoded in the effective susceptibilities. In particular, they allow a straightforward determination of the nonlinear index n 2 of vacuum by analogy with Eq. (1). While a few previous works have given expressions for n 2 of vacuum [36, 37], they have worked with the Euler-Heisenberg model 1 rather than the generalized one adopted here, and explicit values given [5, 38] have only been estimates that lie somewhat below the optimized value 2 .

The paper is organized as follows. In Sec. II we lay the theoretical foundations by specifying the field normalization and Lagrangian we shall use, the latter being subject to the restrictions of Lorentz invariance and local interactions mentioned above. We also indicate how particularly important models of NLED fit into this generalized framework. In Sec. III, we develop the effective medium description by explicitly separating the total fields into a strong background and a weak probe, then linearizing the wave equations in the fields of the latter. We pay particular attention to plane probe waves and their eigenstates within the effective medium, i.e., their refractive indices and polarizations, and we use this formalism to derive some known results in the case of static background fields. We consider an intense pump wave as a background in Sec. IV, showing how elliptical polarization of the pump can be taken into account and deriving the nonlinear index of vacuum by analogy with its definition in standard optical media. In Sec. V, we generalize the Lagrangian to include coupling to an axion field of arbitrary mass, yielding an effective theory of NLED which is nonlocal and thus characterized by dispersion. The analysis is carried through as before, and the key differences are emphasized. We summarize our findings and conclude in Sec. VI.

## II. PRELIMINARIES

We begin by establishing some theoretical foundations. First, we introduce a convenient normalization for the electromagnetic fields which simplifies the writing of many equations. We then state and discuss the most general form of the Lagrangian for the fields given a set of reasonable constraints. This Lagrangian has three free parameters, and we finish this section by identifying the subsets of these parameters corresponding to two particularly important models of NLED: those of Euler-Heisenberg (EH) and Born-Infeld (BI).

1 In the Euler-Heisenberg context, Eq. (2.21) of [37] gives a general result for propagating background fields that includes an arbitrary angle between the propagation directions of pump and probe, and even elliptical polarization, though how arbitrary polarization states enter is rather implicit.

2 Equation (13.8.9) of [5] gives only an estimate for n 2 as it neglects the tensorial nature of the response and is derived from a nonlinear permittivity that does not apply to weak probe waves on a strong background. An order of magnitude is extracted from numerical simulations in [38], but not in the optimized scenario with the probe and background orthogonally polarized and counter-propagating.

## A. Field normalization and equations

We work in Minkowski (flat) space, so that a 3-vector description of the electromagnetic fields may be straightforwardly applied. To avoid overuse of the fundamental constants glyph[epsilon1] 0 and µ 0 (respectively, the permittivity and permeability of free space), it is convenient to use the following rescaled definitions for the electric and magnetic fields:

$$E = \sqrt { \epsilon _ { 0 } } \, E _ { \text{SI} } \,, \quad D = \frac { D _ { \text{SI} } } { \sqrt { \epsilon _ { 0 } } } \,, \quad B = \frac { B _ { \text{SI} } } { \sqrt { \mu _ { 0 } } } \,, \quad \text{H} = \sqrt { \mu _ { 0 } } \, \text{H} _ { \text{SI} } \,,$$

where the subscript 'SI' indicates the corresponding fields expressed in SI units. With these definitions, each of the fields E D B , , and H has exactly the same units (the square root of an energy density), and the Maxwell equations in the absence of free charges and currents take the following form:

$$\nabla \cdot { \mathbf B } = 0 \,, \quad \nabla \times { \mathbf E } + \partial _ { c t } { \mathbf B } = 0 \,,$$

$$\nabla \cdot { \mathbf D } = 0 \,, \quad \nabla \times { \mathbf H } - \partial _ { c t } { \mathbf D } = 0 \,,$$

where c = 1 / √ glyph[epsilon1] 0 µ 0 is the speed of light in vacuum. Equations (3a) are automatically satisfied when E and B are defined in the standard covariant formulation as components of the antisymmetric tensor F µν , itself defined as the exterior derivative of the four-potential A µ [39]; in effect they are consistency conditions that allow such a writing to take place. By contrast, Eqs. (3b) are a convenient writing of the Euler-Lagrange equations found by extremizing the action with respect to variations of A µ , where D and H are defined as

$$D \doteq \frac { \partial \mathcal { L } } { \partial \mathbf E } \,, \quad \mathbf H \doteq - \frac { \partial \mathcal { L } } { \partial \mathbf B } \,.$$

While wave equations (3) can only be fully solved once the constitutive equations relating D and H to E and B are specified, Eqs. (4) indicate that these relations are fully determined once the Lagrangian L ( E , B ) is.

## B. Parametrized Lagrangian for NLED

There are only two scalar quantities invariant under proper orthochronous Lorentz transformations (i.e., those continuously connected to the identity, requiring no spatial reflection or time reversal) which can be constructed from the electromagnetic fields alone [39]:

$$\mathcal { F } \doteq - \frac { 1 } { 4 } F _ { \mu \nu } F ^ { \mu \nu } = \frac { 1 } { 2 } \left ( E ^ { 2 } - B ^ { 2 } \right ) \,, \quad \mathcal { G } \doteq - \frac { 1 } { 4 } F _ { \mu \nu } \widetilde { F } ^ { \mu \nu } = \mathbf E \cdot \mathbf B \,,$$

where ˜ F µν = 1 2 glyph[epsilon1] µναβ F αβ is the Hodge dual of F µν , glyph[epsilon1] µναβ being the completely antisymmetric Levi-Civita symbol with glyph[epsilon1] 0123 = 1. Therefore, a Lorentz invariant Lagrangian containing the electromagnetic fields alone must depend only on F and G , and if we further assume that the Lagrangian is spatially and temporally local (i.e., interactions are purely of the 'contact' type), then its value at a given point of spacetime is straightforwardly a function of F and G at the same point: L ( x ) = L F ( ( x , ) G ( x )). Such a Lagrangian for NLED is said to be of the Pleba´ski class [26]. n We recognize F itself as the standard Lagrangian for electrodynamics in vacuum; it yields the trivial constitutive relations D = E and H = B when plugged into Eqs. (4), whence we recover the usual Maxwell equations in vacuum when inserted into Eqs. (3). Moreover, if G is added to the Lagrangian with some constant coefficient, it is straightforward to show that it has no effect on the field equations 3 , and we are thus free to exclude the occurrence of such a term. We thereby conclude that, at lowest order, the Lagrangian is simply F .

In the weakly nonlinear regime, the Lagrangian can be expanded in powers of F and G [16, 17]. The first nonlinearities will be due to terms quadratic in F and G (i.e., quartic in the fields), and we parametrize their contribution as follows:

$$\mathcal { L } = \mathcal { F } + \delta _ { 1 } \, \frac { 1 } { 2 } \mathcal { F } ^ { 2 } + \delta _ { 2 } \, \frac { 1 } { 2 } \mathcal { G } ^ { 2 } + \delta _ { 3 } \, \mathcal { F } \mathcal { G } \, + \dots$$

The 'post-Maxwellian' parameters [29] δ 1 , δ 2 and δ 3 have units of inverse energy density. They are not completely arbitrary: we will later see that they must satisfy certain inequalities in order for causality to be respected. The

3 This is related to the fact that, when written in terms of A µ , G turns out to be a total derivative [39].

term proportional to δ 3 is often neglected as it breaks invariance under spatial reflection ( P ) or time reversal ( T ) transformations, which preserve the sign of F while causing G to flip sign. Although QED is invariant under P and T , the full Standard Model is not, P -invariance being broken by weak interactions [40, 41]. Therefore, if we wish to include possible deviations from QED in our description, there is no fundamental reason why δ 3 should vanish, and we keep it here for the sake of completeness.

Substituting into Eqs. (4), we find the nonlinear constitutive relations

$$\mathbf D = \frac { \partial \mathcal { L } } { \partial \mathcal { F } } \, \mathbf E + \frac { \partial \mathcal { L } } { \partial \mathcal { G } } \, \mathbf B \,, \quad \mathbf H = \frac { \partial \mathcal { L } } { \partial \mathcal { F } } \, \mathbf B - \frac { \partial \mathcal { L } } { \partial \mathcal { G } } \, \mathbf E \,,$$

where

$$\frac { \partial \mathcal { L } } { \partial \mathcal { F } } = 1 + \delta _ { 1 } \, \mathcal { F } + \delta _ { 3 } \, \mathcal { G } + \dots, \quad \frac { \partial \mathcal { L } } { \partial \mathcal { G } } = \delta _ { 2 } \, \mathcal { G } + \delta _ { 3 } \, \mathcal { F } + \dots$$

The first nonlinear terms in D and H are thus of third order in the fields, and can be thought of as defining the third-order susceptibilities of vacuum by analogy with media. There is a key difference, of course: nonlinearities in media tend to be dominated by a single third-order susceptibility, governing the contribution to D which is cubic in the electric field 4 . In vacuum, however, the requirement of Lorentz invariance restricts the possible nonlinearities to those described by Eqs. (7)-(8), and thus does not allow this term to occur on its own. An important consequence of this difference is the following: since plane waves satisfy F = 0 = G = 0, they behave in vacuum NLED just as in the linear theory of Maxwell, with the same dispersion relation, ω = ck . In the language of nonlinear optics, we may say that vacuum does not exhibit self-phase modulation (SPM), but only cross-phase modulation (XPM): interactions are induced only between different plane waves [38]. In media, on the other hand, SPM and XPM both occur, with different (though closely related) nonlinear indices (as defined in Eq. (1)) [6].

## C. Euler-Heisenberg and Born-Infeld models

Particular models yield particular values of and/or relations between the coefficients δ i . We here consider two specific cases.

The EH effective Lagrangian [9, 10] is derived from QED by summing over all Feynman diagrams containing a single electron-positron loop. In performing the summation, it is assumed that the electromagnetic fields themselves are constant over the loop, which leads to a local effective theory 5 . Therefore, when expanded to quartic order in the fields [8], the EH Lagrangian takes the form (6), with δ 3 = 0 due to the P -and T -invariance of QED, and with the particular values

$$\delta _ { 1 } ^ { ( \text{EH} ) } = \frac { 1 6 } { 4 5 } \alpha ^ { 2 } \frac { \lambda _ { e } ^ { 3 } } { m _ { e } c ^ { 2 } } \approx 1 3. 3 \times 1 0 ^ { - 1 2 } \mu m ^ { 3 } / J \,, \quad \delta _ { 2 } ^ { ( \text{EH} ) } = \frac { 7 } { 4 } \delta _ { 1 } ^ { ( \text{EH} ) } \approx 2 3. 3 \times 1 0 ^ { - 1 2 } \mu m ^ { 3 } / J \,.$$

Here, α ≈ 1 / 137 is the fine structure constant, m e is the mass of the electron, and λ e = glyph[planckover2pi1] /m c e is the reduced Compton wavelength of the electron.

The BI model [15] is derived from the postulate that there exists a fundamental upper limit on the field strength, thus regularizing the self-energy of charged point particles. This model is also P - and T -invariant so that δ 3 = 0, but it predicts δ 1 = δ 2 . = δ (BI) , in strict disagreement with the EH result given above. It thus contains one free parameter, usually written as the maximum absolute field strength b , where δ (BI) = 1 /b 2 . No precise value is predicted, though Born and Infeld considered that the absolute field strength should be approximately that produced by an electron at its own classical radius, and using this prescription one finds

$$\delta ^ { ( \text{BI} ) } \sim 4 \pi \, \alpha ^ { 3 } \frac { \lambda _ { e } ^ { 3 } } { m _ { e } c ^ { 2 } } \approx 3. 4 3 \times 1 0 ^ { - 1 2 } \, \mu m ^ { 3 } / J \,,$$

about a factor of 4 smaller than δ (EH) 1 (or a factor of 7 smaller than δ (EH) 2

).

4 In standard terminology, it is only this term that is referred to as the 'Kerr effect'.

5 This amounts to assuming that the typical photon wavelengths are much larger than the Compton wavelength of the electron, which gives the characteristic size of the electron-positron loop. The loop can then be considered as a point-like vertex. Since λ e ∼ 10 -12 m, this is a good approximation at optical wavelengths glyph[greaterorsimilar] 10 -7 m.

## III. EFFECTIVE MEDIUM DESCRIPTION

In this section, we develop the analogy between the 'dressed' vacuum (including strong electromagnetic fields) and an optical medium. The Lagrangian and wave equations are explicitly decomposed into a background term describing the strong fields alone, and the lowest-order correction due to the presence of the probe. A general equation for the probe wave eigenstates is derived, and some known results for the case of static background fields are reproduced.

## A. Decomposition into background and probe fields

Much like in gravity, where we consider test particles assumed light enough not to have any significant effect on the gravitational field and whose motion is thus entirely determined by the spacetime metric already present, we wish here to consider probe waves propagating in a vacuum whose optical properties have been altered by the presence of strong fields, the probe waves being too weak to contribute to this alteration themselves. To this end, we decompose the total field into a sum of two terms: a background field, much the stronger of the two, entirely responsible for the alteration of the optical properties of the vacuum; and a significantly weaker probe field whose propagation through the altered vacuum we wish to solve for. The 'dressed' vacuum (i.e., the combination of vacuum plus background fields) can be considered as a medium in its own right. The insensitivity of the properties of this effective medium to the presence of the probe implies that the wave equations for the probe fields will be linear, or equivalently that the part of the Lagrangian relevant to the probe will be quadratic in those same fields.

Explicitly, let us write the total fields as E = E 0 + e and B = B 0 + b , where E 0 and B 0 represent the background fields while e and b are the probe fields. The Lagrangian is written as a Taylor series in the latter:

$$\mathcal { L } = \mathcal { L } _ { 0 } + \frac { \partial \mathcal { L } } { \partial E _ { i } } \Big | _ { 0 } \ e _ { i } + \frac { \partial \mathcal { L } } { \partial B _ { i } } \Big | _ { 0 } \ b _ { i } + \frac { 1 } { 2 } \, \frac { \partial ^ { 2 } \mathcal { L } } { \partial E _ { i } \partial E _ { j } } \Big | _ { 0 } \ e _ { i } e _ { j } + \frac { 1 } { 2 } \, \frac { \partial ^ { 2 } \mathcal { L } } { \partial B _ { i } \partial B _ { j } } \Big | _ { 0 } \ b _ { i } b _ { j } + \frac { \partial ^ { 2 } \mathcal { L } } { \partial E _ { i } \partial B _ { j } } \Big | _ { 0 } \ e _ { i } b _ { j } + \dots \quad ( 1 1 )$$

The subscript '0' indicates that the quantity in question is to be evaluated for the background fields E 0 and B 0 , which are taken to be solutions of the full nonlinear wave equations, extremizing the action by definition. Therefore, the terms linear in e and b in Eq. (11), representing the first-order variation of L 0 , give zero contribution to the action and can be removed. The first non-trivial terms involving the probe are those quadratic in the e and b fields, and we define this quadratic part as the effective Lagrangian for the probe:

$$\mathcal { L } _ { \text{probe} } \doteq \frac { 1 } { 2 } \left ( e ^ { 2 } - b ^ { 2 } \right ) + \frac { 1 } { 2 } \text{e} ^ { T } \delta \mathcal { L } _ { \text{ee} } \text{e} + \frac { 1 } { 2 } \text{b} ^ { T } \delta \mathcal { L } _ { \text{bb} } \text{b} + \text{e} ^ { T } \delta \mathcal { L } _ { \text{eb} } \text{b} \,,$$

where the superscript 'T' indicates the transpose, and where we have introduced the matrices

$$\ e p v e \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d } { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n D } { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { [ \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d n \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d {  e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { $ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { c e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { n e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { 1 } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }, } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }$$

These matrices clearly vanish when the full Lagrangian takes the standard Maxwell form (i.e., when L = F ), in which case L probe is simply the Maxwell Lagrangian for the probe fields. Differentiating Eqs. (7) as prescribed by the definitions in Eqs. (13), we may write explicit expressions for the matrices δ L ee , δ L bb and δ L eb :

$$[ \delta \mathcal { L } _ { \text{ee} } ] _ { i j } & = ( \delta _ { 1 } \mathcal { F } _ { 0 } + \delta _ { 3 } \mathcal { G } _ { 0 } ) \, \delta _ { i j } + \delta _ { 1 } E _ { 0, i } E _ { 0, j } + \delta _ { 2 } B _ { 0, i } B _ { 0, j } + \delta _ { 3 } \left ( E _ { 0, i } B _ { 0, j } + B _ { 0, i } E _ { 0, j } \right ) \,, \\ [ \delta \mathcal { L } _ { \text{bb} } ] _ { i j } & = - \left ( \delta _ { 1 } \mathcal { F } _ { 0 } + \delta _ { 3 } \mathcal { G } _ { 0 } \right ) \delta _ { i j } + \delta _ { 1 } B _ { 0, i } B _ { 0, j } + \delta _ { 2 } E _ { 0, i } E _ { 0, j } - \delta _ { 3 } \left ( E _ { 0, i } B _ { 0, j } + B _ { 0, i } E _ { 0, j } \right ) \,, \\ [ \delta \mathcal { L } _ { \text{eb} } ] _ { i j } & = ( \delta _ { 2 } \mathcal { G } _ { 0 } + \delta _ { 3 } \mathcal { F } _ { 0 } ) \, \delta _ { i j } - \delta _ { 1 } E _ { 0, i } B _ { 0, j } + \delta _ { 2 } B _ { 0, i } E _ { 0, j } + \delta _ { 3 } \left ( E _ { 0, i } E _ { 0, j } - B _ { 0, i } B _ { 0, j } \right ) \,,$$

where δ ij is the Kronecker delta. Equivalently, using vector and matrix notation, we have

$$\delta \mathcal { L } _ { \text{ec} } & = ( \delta _ { 1 } \mathcal { F } _ { 0 } + \delta _ { 3 } \mathcal { G } _ { 0 } ) \, 1 + \delta _ { 1 } \mathbf E _ { 0 } \mathbf E _ { 0 } ^ { T } + \delta _ { 2 } \mathbf B _ { 0 } \mathbf B _ { 0 } ^ { T } + \delta _ { 3 } \left ( \mathbf E _ { 0 } \mathbf B _ { 0 } ^ { T } + \mathbf B _ { 0 } \mathbf E _ { 0 } ^ { T } \right ), \\ \delta \mathcal { L } _ { \text{bb} } & = - \left ( \delta _ { 1 } \mathcal { F } _ { 0 } + \delta _ { 3 } \mathcal { G } _ { 0 } \right ) \, 1 + \delta _ { 1 } \mathbf B _ { 0 } \mathbf B _ { 0 } ^ { T } + \delta _ { 2 } \mathbf E _ { 0 } \mathbf E _ { 0 } ^ { T } - \delta _ { 3 } \left ( \mathbf E _ { 0 } \mathbf B _ { 0 } ^ { T } + \mathbf B _ { 0 } \mathbf E _ { 0 } ^ { T } \right ), \\ \delta \mathcal { L } _ { \text{eb} } & = \left ( \delta _ { 2 } \mathcal { G } _ { 0 } + \delta _ { 3 } \mathcal { F } _ { 0 } \right ) \, 1 - \delta _ { 1 } \mathbf E _ { 0 } \mathbf B _ { 0 } ^ { T } + \delta _ { 2 } \mathbf B _ { 0 } \mathbf E _ { 0 } ^ { T } + \delta _ { 3 } \left ( \mathbf E _ { 0 } \mathbf E _ { 0 } ^ { T } - \mathbf B _ { 0 } \mathbf B _ { 0 } ^ { T } \right ),$$

where 1 is the 3 × 3 identity matrix, and where the vector bi-products of the form uv T are 'outer products' yielding matrices rather than scalars. Importantly, the δ L -matrices can be non-zero even when F 0 and G 0 vanish. Therefore, a single plane wave, though in some sense a 'linear' solution of the wave equations, will nonetheless generate an effective medium as it will affect a probe wave whose propagation direction is not equal to its own 6 .

The probe Lagrangian (12) allows us to treat e and b as the 'full' electromagnetic fields, the background fields no longer being treated dynamically but rather having been subsumed into the definition of the effective medium. That is, Eqs. (3) may be applied to the probe fields alone, and the associated constitutive relations are found by inserting L probe into Eqs. (4):

$$\mathrm d = ( 1 + \overline { \delta \mathcal { L } } _ { \text{ee} } ) \, \mathrm e + \overline { \delta \mathcal { L } } _ { \text{eb} } \mathbf b \,, \quad \mathbf h = ( 1 - \overline { \delta \mathcal { L } } _ { \text{bb} } ) \, \mathbf b - \overline { \delta \mathcal { L } } _ { \text{eb} } ^ { T } \mathbf e \,.$$

Here, we have introduced an overbar on the δ L -matrices to indicate a spacetime average over the wavelength and period of the probe; equivalently, the overbar selects their 'slowly-varying' component with respect to the oscillations of the probe. This ensures that d and h inherit the same carrier wave as e and b , differing only in the form of their slowly-varying envelope. In general (and particularly when they are provided by a propagating wave) the background fields are highly oscillatory, and the δ L -matrices defined in Eqs. (14)-(15) will inherit some of this oscillatory behavior. However, since the δ L -matrices influence the probe via an accumulated phase, the highly oscillatory terms can (to a good approximation) usually be neglected 7 .

Rearranging Eqs. (16), and neglecting products of the δ L -matrices (to be consistent with our neglect of higher-order terms in the Lagrangian (6)), we find

$$\mathrm d = ( 1 + \overline { \delta \mathcal { L } } _ { \text{ee} } ) \, \mathrm e + \overline { \delta \mathcal { L } } _ { \text{eb} } \mathbf h \,, \quad \mathbf b = ( 1 + \overline { \delta \mathcal { L } } _ { \text{bb} } ) \, \mathbf h + \overline { \delta \mathcal { L } } _ { \text{eb} } ^ { \text{T} } \mathbf e \,.$$

Equations (17) are in the standard form with respect to which the susceptibilities of an optical medium are defined. We may thus identify the effective electric, magnetic and magnetoelectric susceptibilities of the 'dressed' vacuum:

$$\chi ^ { e } = \overline { \delta \mathcal { L } } _ { \text{ee} } \,, \quad \chi ^ { m } = \overline { \delta \mathcal { L } } _ { \text{bb} } \,, \quad \alpha = \overline { \delta \mathcal { L } } _ { \text{eb} } \,.$$

## B. Plane probe waves in the effective medium

For definiteness, and without loss of generality, we take the probe to be propagating in the -z -direction. The convenience of this choice stems from the fact that, when using a right-handed coordinate system, projections onto the xy -plane (with the z -axis pointing out of the page) intuitively represent what is 'seen' by the probe during its propagation. Using the fact that the averaged δ L -matrices are (by definition) slowly-varying with respect to the wavelength and period of the probe, we may locally decompose its electric field into a slowly-varying envelope and a carrier wave:

$$\mathrm e \left ( z, t \right ) = \mathrm R e \left \{ \mathrm e ^ { ( 0 ) } e ^ { - i k z - i \omega t } \right \} = \frac { 1 } { 2 } \mathrm e ^ { ( 0 ) } e ^ { - i k z - i \omega t } + \mathrm c. c. \,,$$

where k &gt; 0 and ω &gt; 0. Analogous expressions hold for b , d and h . The vector e (0) determines the amplitude and polarization of the electric field, and likewise for d (0) , etc., while the ratio of ω to k gives the phase velocity of the wave: ω/k = c/n , where n is the refractive index. In the weakly nonlinear regime we are considering, n will remain very close to 1, in which case it is more convenient to express this relation in the form

$$1 + \delta n = \frac { c k } { \omega } \,. \quad \quad \quad$$

We wish to determine the refractive index variation δn , which satisfies 0 ≤ δn glyph[lessmuch] 1 8 and which will (to lowest order) be quadratic in the background fields E 0 and B 0 .

6 In this respect, the original DeLLight proposal [25] is over-complicated as it suggests using two counter-propagating pump beams to engender a nontrivial refractive index profile as seen by a probe. One of the key points of this paper (examined further in Sec. IV) is that a single pump beam is sufficient for this purpose.

7 This 'rotating wave approximation' is a standard procedure in nonlinear optics; see, e.g., Ref. [6]. The rapidly oscillating terms in δ L become significant only when phase matching occurs, i.e., when a certain combination of the wave vectors and frequencies involved generates another carrier wave which is itself 'on-shell', with frequency and wave vector approximately satisfying the dispersion relation ω = ck . For the quartic nonlinear Lagrangian considered here, there will be a total of four such waves in any given combination, and these processes are typically referred to as four-wave mixing . See [37, 42] for four-wave mixing in vacuum NLED.

8 The positivity of δn stems from the requirement of causality in Special Relativity, i.e., that signals cannot propagate faster than the speed of light in vacuum, c . Generally speaking, this applies not to the phase velocity ω/k but to the group velocity d ω/ k d , for which we may define a 'group index' n g = n + ω d n/ d ω such that the group velocity is c/n g . Then the causality condition is simply n g ≥ 1. In the present case, due to the local nature of the electromagnetic self-interaction encoded in the Lagrangian (6), n will be independent of ω , so the phase and group velocities are identical and this condition reduces to δn ≥ 0. In Sec. V, we shall examine a dispersive case where δn can be negative, yet the positivity of δn g is still respected.

Even before accounting for the constitutive equations (16) relating e and b to d and h , the probe must satisfy the independent set of equations (3). For a plane wave, these become

$$k \cdot b ^ { ( 0 ) } = 0 \,, \quad c \, \mathbf k \times \math$$

$$\mathbf k \cdot \mathbf d ^ { ( 0 ) } = 0 \,, \quad c \mathbf k$$

where in the present case we have k = -k z ˆ. In each line of Eqs. (21), the second equation implies the first, so that these give only two independent equations rather than four. Considering therefore only the second equation of each line, and using the definition of δn given in Eq. (20), Eqs. (21) may be written as

$$b ^ { ( 0 ) } = - \left ( 1 + \delta n \right ) \Omega _ { z } \$$

$$\mathfrak { d } ^ { ( 0 ) } & = \left ( 1 + \delta n \right ) \O$$

where we have defined the 3 × 3 matrix

$$\Omega _ { z } = \left [ \begin{array} { c c c } 0 & - 1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 0 \end{array} \right ]$$

such that, when acting on a 3-dimensional vector v , we have Ω z v = ˆ z × v . We also note that, since Ω T z = -Ω , the z ordering is faithfully represented when acting on a transposed vector: v T Ω z = -(Ω z v ) T = -(ˆ z × v ) T = ( v × ˆ) z T .

It now remains to impose the constitutive relations (16) for the probe. Using Eq. (22a) so that both d (0) and h (0) can be written directly in terms of e (0) (with no reference to b (0) ), and neglecting products of small quantities, we have:

$$\mathrm d ^ { ( 0 ) } = \left [ \mathbb { 1 } + \overline { \delta \mathcal { L } } _ { \text{ee} } - \overline { \delta \mathcal { L } } _ { \text{eb} } \Omega _ { z } \right ] \mathrm e ^ { ( 0 ) } \,,$$

$$\mathbf h ^ { ( 0 ) } & = \left [ - \left ( 1 + \delta n \right ) \Omega _ { z } + \overline { \delta \mathcal { L } } _ { \mathbf b b } \Omega _ { z } - \overline { \delta \mathcal { L } } _ { \mathbf e b } ^ { T } \right ] \mathbf e ^ { ( 0 ) } \,.$$

Equations (24) are related to each other using Eq. (22b), thus yielding a single homogeneous equation involving the vector e (0) :

$$\left [ \mathfrak { 1 } + \Omega _ { z } ^ { 2 } + 2 \, \delta n \, \Omega _ { z } ^ { 2 } + \overline { \delta \mathcal { L } } _ { \text{ec} } - \Omega _ { z } \, \overline { \delta \mathcal { L } } _ { \text{bb} } \, \Omega _ { z } - \overline { \delta \mathcal { L } } _ { \text{eb} } \, \Omega _ { z } - \left ( \overline { \delta \mathcal { L } } _ { \text{eb} } \, \Omega _ { z } \right ) ^ { T } \right ] \text{e} ^ { ( 0 ) } = 0 \,.$$

The value of δn is found by requiring the determinant of the operator in square brackets in Eq. (25) to vanish. Note that, as Ω 2 z = diag {-1 , -1 0 , } , δn will appear only at quadratic order in the determinant (as could have been expected, since the two solutions correspond to the two possible polarizations of the probe wave). Moreover, since 1 + Ω 2 z = diag { 0 0 1 , , } and all other terms are of first order in δn and the δ L -matrices, then to lowest nontrivial order only the xy -projection of Eq. (25) need be considered 9 . Using expressions (15) for the δ L -matrices and the above-mentioned identification of Ω z with the cross product operator ˆ z × , Eq. (25) reduces to the following 2 × 2 eigenvalue problem:

$$\frac { 1 } { 2 } \left [ \delta _ { 1 } \, \overline { \mathcal { E } \mathcal { E } ^ { T } } + \delta _ { 2 } \, \overline { \mathcal { B } \mathcal { B } ^ { T } } + \delta _ { 3 } \left ( \overline { \mathcal { E } \mathcal { B } ^ { T } } + \overline { \mathcal { B } \mathcal { E } ^ { T } } \right ) \right ] \mathbf e _ { \perp } ^ { ( 0 ) } = \delta n \, \mathbf e _ { \perp } ^ { ( 0 ) } \,,$$

where the subscript ' ⊥ ' indicates the projection onto the xy -plane, and where we have defined (in very similar fashion to Eq. (20) of [28])

$$\mathcal { E } \ & = \mathbf E _ { 0, \perp } - \hat { z } \times \mathbf B _ { 0, \perp } \ = - \hat { z } \times ( \hat { z } \times \mathbf E _ { 0 } ) - \hat { z } \times \mathbf B _ { 0 } \,, \\ \mathcal { B } \ & = \mathbf B _ { 0, \perp } + \hat { z } \times \mathbf E _ { 0, \perp } \ = - \hat { z } \times ( \hat { z } \times \mathbf B _ { 0 } ) + \hat { z } \times \mathbf E _ { 0 } \,.$$

$$( 2 7 )$$

It is clear that E and B lie in the xy -plane, and from their definition it immediately follows that B = ˆ z × E and E = - ×B ˆ z . Therefore, for a given propagation direction, the behavior of the probe wave is determined by a single orthogonal vector formed from E 0 and B 0 .

Equations (26) and (27) are one of the main results of this paper, giving the eigenstates of the probe in the effective medium generated by an arbitrary configuration of strong background fields. They form the basis of the analysis up to the end of Sec. IV (before nonlocal corrections induced by an axion field are considered in Sec. V).

9 To see this, write Eq. (25) in the form ( ˆˆ zz T + δM ) e (0) = 0, where δM is of lowest order in δn and the δ L -matrices. At zeroth-order, the longitudinal component e (0) z vanishes, while the transverse component e (0) ⊥ is determined by the xy -projection of δM .

## C. Constant background fields: the DC Kerr and Cotton-Mouton effects

As an illustrative example, we consider the simplest case where the background fields are constant, or at least slowly-varying with respect to the wavelength and period of the probe. This case has already been analysed in some detail in the literature [28, 43-45], though with particular emphasis on the predictions of the EH model. We give a quick run-through of the various results here, showing that they are indeed reproduced by the effective medium framework we have used, and paving the way for the analysis of a propagating wave as background in Sec. IV.

## 1. Refractive indices as eigenvalues

Since the background fields are slowly-varying, the overbars in Eq. (26) are redundant, the 2 × 2 eigenvalue equation reducing to

$$\frac { 1 } { 2 } \left [ \delta _ { 1 } \, \mathcal { E } ^ { T } + \delta _ { 2 } \, \mathcal { B B } ^ { T } + \delta _ { 3 } \left ( \mathcal { E B } ^ { T } + \mathcal { B E } ^ { T } \right ) \right ] \mathbf e _ { \perp } ^ { \left ( 0 \right ) } = \delta n \, \mathbf e _ { \perp } ^ { \left ( 0 \right ) } \,,$$

where the vectors E and B are now to be considered as constant. Since E and B have the same magnitude and are perpendicular to each other, we may use the orthonormal vectors ˆ = E E / |E| and ˆ = B B / |B| as a basis in the xy -plane. The ordered vectors { ˆ E , ˆ B , ˆ z } form a right-handed orthonormal basis, ˆ and E ˆ being analogous to the standard B

Cartesian basis vectors ˆ and ˆ, respectively. x y In the { ˆ E , ˆ B } basis, Eq. (28) can be written in matrix notation as

$$\left [ \begin{array} { c c } \delta _ { 1 } & \delta _ { 3 } \\ \delta _ { 3 } & \delta _ { 2 } \end{array} \right ] \left [ \begin{array} { c } e _ { \mathcal { E } } ^ { ( 0 ) } \\ e _ { \mathcal { B } } ^ { ( 0 ) } \end{array} \right ] = \frac { \delta n } { \frac { 1 } { 2 } \left | \mathcal { E } \right | ^ { 2 } } \left [ \begin{array} { c } e _ { \mathcal { E } } ^ { ( 0 ) } \\ e _ { \mathcal { B } } ^ { ( 0 ) } \end{array} \right ].$$

This is readily solved. There are two refractive indices (corresponding to two polarizations of the probe wave),

$$\delta n _ { \pm } = \delta _ { \pm } \cdot \frac { 1 } { 2 } \left | \mathcal { E } \right | ^ { 2 } \,,$$

where the coefficients δ ± are the eigenvalues of the matrix on the left-hand side of Eq. (29):

$$\delta _ { \pm } = \frac { 1 } { 2 } \left ( \delta _ { 1 } + \delta _ { 2 } \pm \sqrt { ( \delta _ { 1 } - \delta _ { 2 } ) ^ { 2 } + ( 2 \delta _ { 3 } ) ^ { 2 } } \right ) \,.$$

By definition, δ + ≥ δ -and δn + ≥ δn -. Equality only holds when δ 1 -δ 2 = δ 3 = 0 (as in the BI model); otherwise the presence of the background fields makes the vacuum birefringent [44], a phenomenon referred to as the DC Kerr effect when the static external field is a pure electric field and the Cotton-Mouton effect when it is a pure magnetic field. The vacuum Cotton-Mouton effect is the basis for the BMV [19, 20] and PVLAS [21, 22] experiments. For a magnetic field oriented perpendicular to the direction of the probe wave, we have |E| 2 = B 2 0 = B 2 0 SI , /µ 0 , and the difference in the two refractive indices is

$$\Delta n = \sqrt { \left ( \delta _ { 1 } - \delta _ { 2 } \right ) ^ { 2 } + \left ( 2 \delta _ { 3 } \right ) ^ { 2 } } \, \frac { B _ { 0, \text{SI} } ^ { 2 } } { 2 \mu _ { 0 } } \,.$$

For a magnetic field of 1 Tesla, this gives, for the EH and BI models,

$$\Delta n ^ { ( \text{EH} ) } = 3. 9 8 \times 1 0 ^ { - 2 4 } \,, \quad \Delta n ^ { ( \text{BI} ) } = 0 \,.$$

The EH value is in agreement with the predictions of Refs. [19-22], while the vanishing of ∆ n (BI) indicates the absence of birefringence in the BI model [26-28].

As noted in footnote 8, respecting causality requires the avoidance of a negative value of δn ± , or equivalently of δ ± . It is straightforward to show that this implies the inequalities δ 1 ≥ 0, δ 2 ≥ 0 and δ 1 δ 2 -δ 2 3 ≥ 0 (in agreement with Eqs. (25) of [28]). Interestingly, using the identification of the effective susceptibilities made in Eqs. (18), and using the simplifying assumption F 0 = G 0 = 0, we find

$$\chi ^ { e } _ { i i } \chi ^ { m } _ { j j } - \left ( \alpha _ { i j } \right ) ^ { 2 } = \left ( \delta _ { 1 } \delta _ { 2 } - \delta _ { 3 } ^ { 2 } \right ) \left ( E _ { 0, i } E _ { 0, j } + B _ { 0, i } B _ { 0, j } \right ) ^ { 2 } \,.$$

So in this case the inequality δ 1 δ 2 -δ 2 3 ≥ 0, derived here from the requirement of causality, is equivalent to χ χ e ii m jj -( α ij ) 2 ≥ 0, previously derived (for any optical medium) in [46] from the requirement of thermodynamic stability.

While the factors δ ± are fixed by the post-Maxwellian parameters entering Eq. (6), the strength of the refractive index change is also proportional to 1 2 |E| 2 . This is simply quadratic in the background fields, but because of the projection and combination required to form E and B , the dependence on relative orientation (between E 0 and B 0 , as well as between these fields and the probe wave vector k ∝ -ˆ) can be rather complicated. After a bit of algebra, z it can be shown that

$$\frac { 1 } { 2 } \left | \mathcal { E } \right | ^ { 2 } = \frac { 1 } { 2 } \left ( \left | \hat { z } \times \mathbf E _ { 0 } \right | ^ { 2 } + \left | \hat { z } \times \mathbf B _ { 0 } \right | ^ { 2 } \right ) + \hat { z } \cdot \left ( \mathbf E _ { 0 } \times \mathbf B _ { 0 } \right ) \,.$$

The first term here is rather simple, in that E 0 and B 0 contribute separately, and with the squared magnitude of their projections onto the xy -plane. The second term is more subtle, as it depends on the relative orientation of E 0 and B 0 . Moreover, it is directionally dependent: whereas the first term depends only on the line along which the wave vector k lies (defined to be the z -axis) and does not vary under the transformation k →-k (i.e., ˆ z →-ˆ), the z second term changes sign under this transformation. In this sense the effective medium behaves as if it were moving with a velocity proportional to the 'Poynting vector' E 0 × B 0 . This anisotropy (which is independent of the probe polarization as it stems only from the magnitude of the vector E ) was described in Ref. [45]. We shall see that it is also present when the background is a plane wave, with co-propagating probe waves seeing no refractive index change while counter-propagating waves experience the strongest effect.

## 3. Eigenpolarizations

The eigenvectors of Eq. (29) give the two eigenpolarizations of e (0) . Since the matrix on the left-hand side is real and symmetric, the eigenvectors are necessarily real and (when normalized) form the columns of a two-dimensional rotation matrix

$$R \left ( \varphi \right ) = \left [ \begin{array} { c c } \cos \varphi & - \sin \varphi \\ \sin \varphi & \cos \varphi \end{array} \right ] \,.$$

The parameter ϕ is simply the angle through which the eigenvectors are rotated with respect to the basis { ˆ E , ˆ B }

(see Fig. 1), and is defined up to a multiple of π since a half-rotation simply flips the signs of the eigenvectors without changing their orientation. ϕ can thus be chosen to lie in the half-open interval ( -π/ 2 , π/ 2], and the matrix on the left-hand side of Eq. (29) can be written as R ϕ DR ( ) -1 ( ϕ ), where D is a diagonal matrix whose entries are the eigenvalues of Eq. (31). We take δ + to be the first diagonal component of D , so that the left column of R ϕ ( ) gives the polarization with the larger refractive index change, δn + . A direct calculation shows that ϕ must satisfy

$$\delta _ { 1 } - \delta _ { 2 } = ( \delta _ { + } - \delta _ { - } ) \cos \left ( 2 \varphi \right ) \,, \quad 2 \delta _ { 3 } = ( \delta _ { + } - \delta _ { - } ) \sin \left ( 2 \varphi \right ) \,.$$

glyph[negationslash]

There are three special cases. First, when δ 3 = 0 and δ 1 -δ 2 = 0 (as in the EH model), we have ϕ = 0 or π/ 2 (depending on the sign of δ 1 -δ 2 ) and the polarizations are aligned with the vectors ˆ and E ˆ . B Second, when δ 1 -δ 2 = 0 and δ 3 = 0, we have ϕ = ± π/ 4, so the polarizations are at 45 ◦ to ˆ and E ˆ . B Third, when both δ 3 = 0 and δ 1 -δ 2 = 0 (as in the BI model), ϕ is undefined, but this is not a problem as there is an absence of birefringence (and hence also of well-defined eigenpolarizations) in this case.

glyph[negationslash]

## IV. OPTICAL KERR EFFECT

In this section we turn to the main focus of this paper: the refractive index change of vacuum engendered by an intense propagating wave or 'pump', which provides the strong background fields described in Sec. III. We derive the dependence on the tilt angle (between the propagation directions of pump and probe), as well as the effect of elliptical polarization of the pump. Finally, we express the results in terms of the wave intensity (rather than the field strength) in order to extract the equivalent of the nonlinear index for vacuum by analogy with Eq. (1).

## A. Fields of a monochromatic pump wave

Let us consider then the fields of a propagating beam, which we assume can be approximated as monochromatic over spacetime regions much larger than the wavelength and longer than the period of the probe. We may again use

'

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00896v3-figure-623f7996-15c7-46ef-90d3-1c4494988e91.jpg)

'

Figure 1. Eigenpolarizations for the electric field amplitude e (0) of the probe, which propagates in the -z -direction (i.e., into the page). They lie in the xy -plane, rotated with respect to ˆ and E ˆ by the angle B ϕ whose value is determined by Eqs. (37). The polarizations corresponding to refractive index changes δn + and δn -of Eq. (30) are given by the left and right columns of R ϕ ( ) of Eq. (36) and shown here in solid and dashed line, respectively. The left and right panels differ by a rotation of π/ 2, which (in effect) leaves the eigenpolarization directions invariant but switches the associated refractive indices. In the EH model, we find ϕ = π/ 2, so the δn + polarization is aligned with ˆ while the B δn -polarization is aligned with ˆ . E

Eqs. (26) and (27), though now the overbars extracting the 'slowly-varying' components of the outer products in Eq. (26) will come into play. We write the pump fields in the form

$$\mathbf E _ { 0 } & = \frac { 1 } { 2 } \, \mathbf E _ { 0 } ^ { ( 0 ) } \, e ^ { i k _ { 0 } \cdot r - i \omega _ { 0 } t } + c. c. \,, \\ \mathbf B _ { 0 } & = \frac { 1 } { 2 } \, \mathbf B _ { 0 } ^ { ( 0 ) } \, e ^ { i k _ { 0 } \cdot r - i \omega _ { 0 } t } + c. c. \,,$$

where, in order to satisfy the Maxwell equations (3) with D = E and H = B (as a single plane wave must, having F 0 = G 0 = 0), we have ω 0 = ck 0 (where k 0 = | k 0 | ), and

$$B _ { 0 } ^ { ( 0 ) } = \frac { k _ { 0 } } { k _ { 0 } } \times \mathbf E _ { 0 } ^ { ( 0 ) } \,.$$

Note that we do not specify the direction of k 0 , whose orientation with respect to the probe wave vector k = -kz ˆ is taken to be arbitrary.

The next step is to work out the vectors E and B entering the matrix in Eq. (26), before application of the overbars. These will be oscillatory just as E and B are, and we may write:

$$\mathcal { E } & = \frac { 1 } { 2 } \, \mathcal { E } ^ { ( 0 ) } \, e ^ { i \mathbf k _ { 0 } \cdot \mathbf r - i \omega _ { 0 } t } + c. c. \,, \\ \mathcal { B } & = \frac { 1 } { 2 } \, \mathcal { B } ^ { ( 0 ) } \, e ^ { i \mathbf k _ { 0 } \cdot \mathbf r - i \omega _ { 0 } t } + c. c. \,.$$

where the amplitude vectors are given by

$$\mathcal { E } ^ { ( 0 ) } & = \mathbf E _ { 0, \perp } ^ { ( 0 ) } - \hat { z } \times \mathbf B _ { 0, \perp } ^ { ( 0 ) } = - \hat { z } \times \left [ \left ( \hat { z } + \frac { k _ { 0 } } { k _ { 0 } } \right ) \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] \,, \\ \mathcal { B } ^ { ( 0 ) } & = \mathbf B _ { 0, \perp } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0, \perp } ^ { ( 0 ) } = - \hat { z } \times \left [ \left ( \hat { z } + \frac { k _ { 0 } } { k _ { 0 } } \right ) \times \mathbf B _ { 0 } ^ { ( 0 ) } \right ] \,,$$

Here, we have used Eq. (39), as well as standard identities concerning two successive applications of the cross product. These vectors evidently lie in the xy -plane, and by construction we again have B (0) = ˆ z ×E (0) and E (0) = - ×B ˆ z (0) . After some further algebra, the squared magnitude of E (0) (and hence also of B (0) ) can be shown to be

$$\dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow
\mathcal { E } ^ { ( 0 ) * } \cdot \mathcal { E } ^ { ( 0 ) } & = \left ( 1 + \frac { \hat { z } \cdot \mathbf k _ { 0 } } { k _ { 0 } } \right ) ^ { 2 } \left | \mathbf E _ { 0 } ^ { ( 0 ) } \right | ^ { 2 } \\ & = ( 1 + \cos \theta ) ^ { 2 } \left | \mathbf E _ { 0 } ^ { ( 0 ) } \right | ^ { 2 } \\ & = 4 \cos ^ { 4 } \frac { \theta } { 2 } \left | \mathbf E _ { 0 } ^ { ( 0 ) } \right | ^ { 2 } \,.$$

Here, we have introduced the tilt angle θ between the wave vectors of the pump and probe (illustrated in Fig. 2). This is defined to be zero when the pump and probe are exactly counter-propagating and ± π when they are exactly co-propagating. Note that the magnitude of E (0) vanishes in the latter case, this being consistent with the observation made after Eqs. (8) that there is no SPM in vacuum; on the other hand, it is maximal when the pump and probe are exactly counter-propagating.

## B. Accounting for elliptical polarization of the pump

Inserting Eqs. (40) into Eq. (26), and implementing the overbars by dropping all rapidly oscillating terms, we are led to the following eigenproblem:

$$\frac { 1 } { 4 } \, R e \left \{ \delta _ { 1 } \, \mathcal { E } ^ { ( 0 ) * } \mathcal { E } ^ { ( 0 ) T } + \delta _ { 2 } \, \mathcal { B } ^ { ( 0 ) * } \mathcal { B } ^ { ( 0 ) T } + \delta _ { 3 } \left ( \mathcal { E } ^ { ( 0 ) * } \mathcal { B } ^ { ( 0 ) T } + \mathcal { B } ^ { ( 0 ) * } \mathcal { E } ^ { ( 0 ) T } \right ) \right \} e ^ { ( 0 ) } = \delta n \, e ^ { ( 0 ) } \,.$$

glyph[negationslash]

Compared to the case of constant background fields studied in Sec. III C, we have here a complication in that the field amplitudes E (0) 0 and B (0) 0 , and by extension the vectors E (0) and B (0) , are generally complex. It is thus no longer convenient to use an orthonormal basis aligned with E (0) and B (0) , since the matrix on the left-hand side of Eq. (43) also depends on E (0) glyph[star] and B (0) glyph[star] , and in general E (0) glyph[star] = E (0) and B (0) glyph[star] = B (0) . The issue is not with overall phases - it is clear that the matrix in question is invariant under equal overall phase rotations of E (0) and B (0) - but with relative phases between the components of these vectors. Such relative phases are directly related to the degree of elliptical polarization of the pump wave. If we were to restrict our attention to a linearly polarized pump, there would be no such relative phase, E (0) and B (0) could be defined to be real, and the eigenproblem would be equivalent to that of Eq. (28) (except for an overall factor of 1 / 2 stemming from the average over rapidly oscillating terms). We can thus expect to recover almost the same results as those of Sec. III C when the pump is linearly polarized. However, with a proper treatment of the complex vectors entering the left-hand side of Eq. (43), the effects of a general elliptical polarization can be fully included in the analysis.

glyph[negationslash]

Up to an overall phase, the complex field amplitudes E (0) 0 and B (0) 0 can always be written as a sum of two orthogonal real vectors, one (which for definiteness we take to be the larger in magnitude of the two) with a purely real coefficient, the other with a purely imaginary coefficient. These vectors necessarily lie in the plane perpendicular to the wave vector k 0 . They describe the principle directions of the polarization, i.e., the major and minor axes of the ellipse formed by the oscillating electric or magnetic field (see Fig. 2). The field lies along the minor axis exactly a quarter of a period after it lies along the major axis, and so the complex components of the field along these two directions appear with a relative phase of π/ 2. We may thus write the field amplitudes using the real vectors E (0) abs and B (0) abs , and again utilizing the fact that B (0) 0 = ( k 0 /k 0 ) × E (0) 0 :

$$\mathbf E _ { 0 } ^ { ( 0 ) } & = e ^ { i \phi _ { 0 } } \left ( \cos \chi \, \mathbf E _ { a b s } ^ { ( 0 ) } + i \sin \chi \, \mathbf B _ { a b s } ^ { ( 0 ) } \right ) \,, \\ \mathbf B _ { 0 } ^ { ( 0 ) } & = e ^ { i \phi _ { 0 } } \left ( \cos \chi \, \mathbf B _ { a b s } ^ { ( 0 ) } - i \sin \chi \, \mathbf E _ { a b s } ^ { ( 0 ) } \right ) \,.$$

where B (0) abs = ( k 0 /k 0 ) × E (0) abs and E (0) abs = -( k 0 /k 0 ) × B (0) abs . With this writing, E (0) 0 and E (0) abs have exactly the same magnitude. The angle χ ∈ -[ π/ 4 , π/ 4] is the so-called ellipticity angle of the polarization ellipse (also illustrated in Fig. 2). When χ = 0, the wave is linearly polarized; when χ = ± π/ 4, it is circularly polarized, in which case the directions of E (0) abs and B (0) abs in the 2D plane can be chosen arbitrarily. χ is defined to be positive when the sense of rotation of the fields is from E (0) abs towards B (0) abs , corresponding to right-handed polarization of the pump.

Figure 2. Relative orientations for a monochromatic pump wave. In ( a ) are shown the vectors E ( ) and t B ( ) (in dashed and t dotted line, respectively) in the xy -plane, with ˆ pointing out of the page and the wave vector of the probe pointing into the z page. These oscillate in time so as to trace out ellipses, rotated with respect to each other by 90 ◦ . The constant real vectors E (0) abs and B (0) abs are defined to lie along the semi-major axes of their respective ellipses, with magnitude equal to the length of the hypotenuse of the right-angled triangle shown. The ellipticity angle χ is formed by the same right-angled triangle, though its sign is determined by the sense of rotation of E ( ) and t B ( ). t In ( b ) are shown the wave vectors of the pump and probe waves forming the tilt angle θ , equal to 0 for exactly counter-propagating waves and ± π for exactly co-propagating waves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00896v3-figure-3c44a74b-2e14-4128-9e88-b5cfc392a5da.jpg)

By a straightforward application of Eqs. (41), we may now write

$$\mathcal { E } ^ { ( 0 ) } & = e ^ { i \phi _ { 0 } } \, \left ( \cos \chi \, \mathcal { E } ^ { ( 0 ) } _ { \text{abs} } + i \sin \chi \, \mathcal { B } ^ { ( 0 ) } _ { \text{abs} } \right ) \,, \\ \mathcal { B } ^ { ( 0 ) } & = e ^ { i \phi _ { 0 } } \, \left ( \cos \chi \, \mathcal { B } ^ { ( 0 ) } _ { \text{abs} } - i \sin \chi \, \mathcal { E } ^ { ( 0 ) } _ { \text{abs} } \right ) \,,$$

where we have defined

$$\mathcal { E } _ { \text{abs} } ^ { ( 0 ) } & = - \hat { z } \times \left [ \left ( \hat { z } + \frac { \mathbf k _ { 0 } } { k _ { 0 } } \right ) \times \mathbf E _ { \text{abs} } ^ { ( 0 ) } \right ] \,, \\ \mathcal { B } _ { \text{abs} } ^ { ( 0 ) } & = - \hat { z } \times \left [ \left ( \hat { z } + \frac { \mathbf k _ { 0 } } { k _ { 0 } } \right ) \times \mathbf B _ { \text{abs} } ^ { ( 0 ) } \right ] \,.$$

It is again fairly straightforward to show that B (0) abs = ˆ z × E (0) abs and E (0) abs = - × B ˆ z (0) abs , so that these are orthogonal and have equal magnitude. Therefore, analogously to Eqs. (44), Eqs. (45) define an ellipse of ellipticity angle χ in the xy -plane. It is particularly interesting that the ellipticity angle remains χ when passing to E (0) and B (0) , since these lie in the xy -plane rather than the plane containing E (0) abs and B (0) abs . One might naively have expected the shape of the ellipse to depend on the orientation of k 0 with respect to ˆ (as the orthogonal projections of the ellipses traced out z by the electric and magnetic fields certainly do depend on the angle they are viewed from). Remarkably, however, it turns out that E (0) and B (0) combine the electric and magnetic fields of the pump wave in just the right way so that χ is invariant with respect to the relative orientation of pump and probe. Only the magnitude of E (0) (and hence of E (0) abs ) varies with the orientation, as shown by Eq. (42).

## C. Nonlinear index of vacuum

Finally, Eqs. (45) are inserted into Eq. (43), upon which it becomes the following eigenproblem:

$$\frac { 1 } { 4 } \left [ \delta _ { 1 } ^ { \prime } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } + \delta _ { 2 } ^ { \prime } \mathcal { B } _ { \text{abs} } ^ { ( 0 ) } \mathcal { B } _ { \text{abs} } ^ { ( 0 ) T } + \delta _ { 3 } ^ { \prime } \, \left ( \mathcal { E } _ { \text{abs} } ^ { ( 0 ) } \mathcal { B } _ { \text{abs} } ^ { ( 0 ) T } + \mathcal { B } _ { \text{abs} } ^ { ( 0 ) } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } \right ) \right ] e ^ { ( 0 ) } = \delta n \, e ^ { ( 0 ) } \,,$$

where

$$\delta _ { 1 } ^ { \prime } + \delta _ { 2 } ^ { \prime } & = \delta _ { 1 } + \delta _ { 2 } \,, \\ \delta _ { 1 } ^ { \prime } - \delta _ { 2 } ^ { \prime } & = ( \delta _ { 1 } - \delta _ { 2 } ) \cos { ( 2 \chi ) } \, \\ \delta _ { 3 } ^ { \prime } & = \delta _ { 3 } \cos { ( 2 \chi ) } \. \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon$$

This is now in a form completely analogous to Eq. (28), with the vectors E (0) abs and B (0) abs being equal in magnitude and orthogonal to each other. They can thus be used to define basis vectors in the xy -plane, ˆ = E E (0) abs / ∣ ∣ ∣ E (0) abs ∣ ∣ ∣ and ˆ = B B (0) abs / ∣ ∣ ∣ B (0) abs ∣ ∣ ∣ , so that the vectors { ˆ E , ˆ B , ˆ z } form a right-handed orthonormal basis. Restricting our attention to the xy -plane in the basis { ˆ E , ˆ B } , Eq. (47) can be written in matrix form as follows:

$$\left [ \begin{array} { c c } \delta _ { 1 } ^ { \prime } & \delta _ { 3 } ^ { \prime } \\ \delta _ { 3 } ^ { \prime } & \delta _ { 2 } ^ { \prime } \end{array} \right ] \left [ \begin{array} { c } e _ { \mathcal { E } } ^ { ( 0 ) } \\ e _ { \mathcal { B } } ^ { ( 0 ) } \end{array} \right ] = \frac { \delta n } { \frac { 1 } { 4 } \left | \mathcal { E } _ { \text{abs} } ^ { ( 0 ) } \right | ^ { 2 } } \left [ \begin{array} { c } e _ { \mathcal { E } } ^ { ( 0 ) } \\ e _ { \mathcal { B } } ^ { ( 0 ) } \end{array} \right ] \,.$$

The eigenvalues of the matrix on the left-hand side of this equation are readily found, and give the possible values of δn . Using Eq. (42) and the fact that E (0) and E (0) abs are defined to have the same magnitude, we have

$$\delta n _ { \pm } = \delta _ { \pm } ^ { \prime } \cos ^ { 4 } \frac { \theta } { 2 } \left | \mathbf E _ { 0 } ^ { ( 0 ) } \right | ^ { 2 },$$

where

$$\delta _ { \pm } ^ { \prime } & = \frac { 1 } { 2 } \left ( \delta _ { 1 } ^ { \prime } + \delta _ { 2 } ^ { \prime } \pm \sqrt { ( \delta _ { 1 } ^ { \prime } - \delta _ { 2 } ^ { \prime } ) ^ { 2 } + ( 2 \delta _ { 3 } ^ { \prime } ) ^ { 2 } } \right ) \\ & = \frac { 1 } { 2 } \left ( \delta _ { 1 } + \delta _ { 2 } \pm \cos { ( 2 \chi ) } \ \sqrt { ( \delta _ { 1 } - \delta _ { 2 } ) ^ { 2 } + ( 2 \delta _ { 3 } ) ^ { 2 } } \right ) \,. \\ \omega _ { 1 } \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \end{cases}$$

We thus see that the ellipticity angle χ directly affects the strength of the birefringence, which vanishes completely in the case of circular polarization. Moreover, Eqs. (37) still hold (being simply multiplied by an overall factor of cos (2 χ )), so that the rotation angle ϕ of the eigenpolarizations with respect to the { ˆ E , ˆ B } basis is independent of χ . Note that the probe eigenstates are linearly polarized no matter the polarization state of the pump, a direct result of the fact that the matrix on the left-hand side of Eq. (43) is real and symmetric. (In Sec. V we shall examine a model where this is no longer the case, allowing complex eigenvectors which encode states of elliptical polarization.)

We are now in a position to give explicit values for the corresponding nonlinear index n 2 , defined (by analogy with Eq. (1)) such that the refractive index change δn = n 2 I , where I is the intensity of the pump wave. For definiteness and simplicity, we consider pump and probe to be exactly counter-propagating ( θ = 0), this being the optimal geometry 10 according to Eq. (50). The energy density of the pump (after averaging over rapidly oscillating terms) is ∣ ∣ ∣ E (0) 0 ∣ ∣ ∣ 2 / 2, and its intensity is found upon multiplication by c . The corresponding values of n 2 are then simply δ ′ × 2 /c . In the EH model we have, for a linearly polarized pump beam,

±

$$n _ { 2, \| } ^ { ( \text{EH} ) } \approx 0. 8 8 8 \times 1 0 ^ { - 3 3 } \, \text{cm} ^ { 2 } / W \,, \quad n _ { 2, \perp } ^ { ( \text{EH} ) } \approx 1. 5 5 5 \times 1 0 ^ { - 3 3 } \, \text{cm} ^ { 2 } / W \,,$$

where the subscripts ' ‖ ' and ' ⊥ ' refer to the probe and pump fields being equally and orthogonally polarized, respectively. On the other hand, for a circularly polarized pump, n 2 no longer depends on the polarization of the probe and is simply the arithmetic mean of the two values given above:

$$n _ { 2, c i r c } ^ { ( E H ) } \approx 1. 2 2 2 \times 1 0 ^ { - 3 3 } \, c m ^ { 2 } / W \,.$$

10 In the context of the EH model, Ref. [36] seems to have been the first to study this particular geometry. Ref. [38] studied a suboptimal geometry in which the waves have parallel polarizations and a tilt angle of 90 ◦ , yielding an overall reduction factor of 7 in the value of n 2 with respect to the optimal value, n (EH) 2 , ⊥ in Eqs. (52).

It can thus be seen that the EH prediction for the nonlinear index of vacuum, though it depends on pump polarization and tilt angle, is on the order of 10 -33 cm 2 / W.

In the BI model instead, δ ′ ± of Eqs. (51) are both equal to δ (BI) and independent of the ellipticity angle χ . We thus have simply

$$n _ { 2 } ^ { ( \text{BI} ) } = \frac { 2 } { c \, b ^ { 2 } } \,,$$

where b 2 is the square of the critical field parametrizing the BI model (expressed in units of energy density). Using the value of δ (BI) in Eq. (10), with b being the field at the classical radius of the electron, this gives

$$n _ { 2 } ^ { ( \text{BI} ) } \approx 0. 2 2 9 \times 1 0 ^ { - 3 3 } \text{cm} ^ { 2 } / W \,.$$

## V. A DISPERSIVE MODEL: COUPLING TO AXIONS

Experiments in NLED have been considered as potentially enabling the detection of the axion [30-34], a hypothetical particle introduced as a possible explanation for strong CP invariance in quantum chromodynamics [11-14], and which has been proposed as a candidate for dark matter [47]. As far as electromagnetism is concerned, the axion field couples directly to G = E B · , and will thus contribute to the effective photon/photon interaction in NLED. However, there are compelling astrophysical [48] and cosmological [49] reasons to consider an axion mass which is significantly smaller than 1 eV, the energy scale of an optical photon. In this case, the Compton wavelength of the axion is long compared to the typical photon wavelength, and the assumption of purely local effective photon/photon interactions made in Sec. II is explicitly broken. Coupling to the axion field is thus not only of potential experimental relevance (though the question of experimental feasibility is beyond the scope of this paper), but is also of theoretical interest as it will lead to a dispersive model of NLED (as was recently illustrated in [50]).

In this section, starting from the electromagnetic Lagrangian of Eq. (6), we couple the electromagnetic field to an axion field. The analysis of previous sections is carried through in a similar manner; the post-Maxwellian parameters and the polarization of the pump are left unspecified, thereby generalizing the results of [50] (which considered the EH model and a linearly polarized pump). We shall assume that the various plane waves are infinite in extent and duration. This allows us to neglect retardation effects due to the non-instantaneous nature of the axion response, including photon-axion oscillations [34]. Instead, we here focus solely on the refractive index change induced by the axion coupling.

## A. Lagrangian and constitutive relations

The total effective Lagrangian, including the coupling to axions, may be written

$$\mathcal { L } = \mathcal { L } _ { \text{EM} } + \mathcal { L } _ { \text{ax} } + \mathcal { L } _ { \text{int} } \,.$$

Here, L EM can be considered as the 'local part' of the effective Lagrangian containing the electromagnetic fields alone, and is just that used in previous sections (and given in Eq. (6)). The next term is the Lagrangian of the free axion field:

$$\mathcal { L } _ { \text{ax} } = \frac { 1 } { 2 } \left ( \partial _ { c t } \phi \right ) ^ { 2 } - \frac { 1 } { 2 } \left ( \nabla \phi \right ) ^ { 2 } - \frac { 1 } { 2 } k _ { \text{ax} } ^ { 2 } \phi ^ { 2 } \,,$$

where k ax = m c/ ax glyph[planckover2pi1] is the wave vector associated to the Compton wavelength of the axion. Finally, the interaction between the axion and electromagnetic fields is described by

$$\mathcal { L } _ { \text{int} } = - \eta \, \phi \, \mathcal { G } \,,$$

where the inverse square of the coupling constant, η -2 , has dimensions of energy per unit length.

Employing the separation into background and probe fields described in Sec. III A, and using a similar decomposition for the axion field φ = φ 0 + δφ , we may write the total Lagrangian as

$$\mathcal { L } & = \mathcal { L } _ { \text{EM} } \left ( \mathbf E _ { 0 } + \mathbf e \,, \, \mathbf B _ { 0 } + \mathbf b \right ) + \mathcal { L } _ { \text{ax} } \left ( \phi _ { 0 } + \delta \phi \right ) + \mathcal { L } _ { \text{int} } \left ( \mathbf E _ { 0 } + \mathbf e \,, \, \mathbf B _ { 0 } + \mathbf b \,, \, \phi + \delta \phi \right ) \\ & \approx \mathcal { L } _ { 0 } + \mathcal { L } _ { \text{probe} }$$

where

$$\mathcal { L } _ { 0 } = \mathcal { L } _ { \text{EM} } \left ( \mathbf E _ { 0 } \,, \, \mathbf B _ { 0 } \right ) + \mathcal { L } _ { \text{ax} } \left ( \phi _ { 0 } \right ) + \mathcal { L } _ { \text{int} } \left ( \mathbf E _ { 0 } \,, \, \mathbf B _ { 0 } \,, \, \phi _ { 0 } \right )$$

is the Lagrangian associated with the background alone and, keeping only terms which are quadratic in the weak fields of the probe,

$$\mathcal { L } _ { \text{probe} } = \mathcal { L } _ { \text{EM,probe} } \left ( \mathbf E _ { 0 } \,, \, \mathbf B _ { 0 } \, ; \, e, \, \mathbf b \right ) + \mathcal { L } _ { \text{ax} } \left ( \delta \phi \right ) - \eta \left ( \phi _ { 0 } \, \mathbf e \cdot \mathbf b + \delta \phi \, \mathbf E _ { 0 } \cdot \mathbf b + \delta \phi \, \mathbf B _ { 0 } \cdot \mathbf e \right ) \,.$$

L EM probe , is exactly the probe Lagrangian derived in Sec. III A, while L ax (which is already purely quadratic) is again just the Lagrangian for a free axion field. The terms proportional to η describe the interplay between the electromagnetic and axion fields associated with the passage of the probe wave. Since we deal with a restricted class of background configurations, we can simplify this term further. Noting that we deal either with plane wave background fields, which satisfy G 0 = 0 and hence φ 0 = 0, or with constant background fields, in which case the term in e · b is a total derivative when expressed in terms of the vector potential (see the discussion following Eqs. (5), including footnote 3), the term φ 0 e · b can be removed from the Lagrangian with no effect on the field equations. Therefore, we consider only the terms proportional to δφ as far as the coupling to the axion field is concerned.

Applying definitions (4) to the probe fields, we have

$$d & = d _ { \text{EM} } - \eta \, \delta \phi \text{B} _ { 0 } \,, \\ \text{h} & = \text{h} _ { \text{EM} } + \eta \, \delta \phi \text{E} _ { 0 } \,,$$

where d EM and h EM are due to L EM probe , alone and are already defined in Eqs. (16). Our aim is to subject these relations to the same treatment as in Sec. III; in particular, to find the new forms of Eqs. (24) relating the amplitudes d (0) and h (0) directly to e (0) , so that a homogeneous linear equation analogous to Eq. (25) is obtained. Our first task, then, is to determine the axion field δφ generated by a probe wave of amplitude e (0) .

## B. Response of axion field to passage of probe wave

The response of δφ to the presence of electromagnetic fields is determined by the following equation of motion:

$$\left [ \partial _ { c t } ^ { 2 } - \nabla ^ { 2 } + k _ { \text{ax} } ^ { 2 } \right ] \delta \phi = - \eta \, \delta \mathcal { G } \,,$$

where we have defined δ G = E 0 · b + B 0 · e . The probe-induced δ G thus acts as a source for δφ , with a simple relationship between their Fourier components:

$$\delta \phi _ { \omega ^ { \prime }, k ^ { \prime } } = \eta \, \frac { \delta \mathcal { G } _ { \omega ^ { \prime }, k ^ { \prime } } } { \left ( \omega ^ { \prime } / c \right ) ^ { 2 } - \left ( k ^ { \prime } \right ) ^ { 2 } - k _ { \text{ax} } ^ { 2 } } \,,$$

glyph[negationslash]

assuming of course that we are not at resonance, i.e., ( ω /c ′ ) 2 -( k ′ ) 2 -k 2 ax = 0. Using Eqs. (38) for the background/pump fields (noting that they reduce to constant fields when k 0 and ω 0 vanish), we have, to lowest order,

$$\delta \mathcal { G } \approx \frac { 1 } { 4 } \left \{ \left [ \mathbf B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] e ^ { i k _ { 0 }, r - i \omega _ { 0 } t } + \left [ \mathbf B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] ^ { * } e ^ { - i k _ { 0 }, r + i \omega _ { 0 } t } \right \} \cdot e ^ { ( 0 ) } e ^ { - i k z - i \omega t } + c. c.,$$

where we have used Eq. (22a), as well as the cyclic invariance of the vector triple product, to write E (0) 0 · b (0) 0 ≈ ( ˆ z × E (0) 0 ) · e (0) . Using relation (64), we can immediately write down the generated axion field:

$$\left ( \stackrel { \lambda } { \lambda } \cdots \overline { - 0 } } { 0 } \right ) \, \stackrel { \lambda } { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \overbrace { \lambda } \cdots \over \quad \, \stackrel { \lambda } { 4 } } { \left ( \omega _ { 0 } + \omega \right ) ^ { 2 } / c ^ { 2 } - \left ( k _ { 0 } - k \acute { z } \right ) ^ { 2 } - k _ { 2 \max } ^ { 2 } } + \frac { \left [ B _ { 0 } ^ { ( 0 ) } + \acute { z } \times \acute { E } _ { 0 } ^ { ( 0 ) } \right ] ^ { * } e ^ { - i k _ { 0 } } r + i \omega _ { 0 } t } { \left ( \omega _ { 0 } - \omega \right ) ^ { 2 } / c ^ { 2 } - \left ( k _ { 0 } + k \acute { z } \right ) ^ { 2 } - k _ { 2 \max } ^ { 2 } } \right ) \cdot e ^ { ( 0 ) } ( e ^ { - i k z - i \omega t } + c. c. \\ \approx \frac { \eta } { 4 } \left ( \frac { \left [ B _ { 0 } ^ { ( 0 ) } + \acute { z } \times \acute { E } _ { 0 } ^ { ( 0 ) } \right ] e ^ { i k _ { 0 } } r + i \omega _ { 0 } t } { 4 \cos ^ { 2 } \frac { \theta } { 2 } \omega _ { 0 } \omega / c ^ { 2 } - k _ { 2 \max } ^ { 2 } } - \frac { \left [ B _ { 0 } ^ { ( 0 ) } + \acute { z } \times \acute { E } _ { 0 } ^ { ( 0 ) } \right ] ^ { * } e ^ { - i k _ { 0 } } r + i \omega _ { 0 } t } { 4 \cos ^ { 2 } \frac { \theta } { 2 } \omega _ { 0 } \omega / c ^ { 2 } + k _ { 2 \max } ^ { 2 } } \right ) \cdot e ^ { ( 0 ) } e ^ { - i k z - i \omega t } + c. c.$$

In the second line, we have expanded the squares in the denominators, neglecting the variation of the refractive index here so that k ≈ ω/c , and used the fact that k 0 · ˆ = z k 0 cos θ where θ is the tilt angle between the pump and probe waves introduced in the previous section (and illustrated in Fig. 2).

## C. Backreaction of axion field on probe fields

Substituting expression (66) for δφ back into the constitutive relations (62) removes the explicit dependence on the axion field, giving d and h in terms of the probe amplitude e (0) alone. Since d EM and h EM are already given in Eqs. (24), we need focus here only on the additional terms δ d = -η δφ B 0 and δ h = η δφ E 0 . Since δφ , E 0 and B 0 are all generally oscillatory, the substitution generates rapidly oscillating terms that are far off-shell, much like when calculating d and h in Sec. III. As there, we retain only those terms whose oscillations are synchronized with those of the probe, in which case we may write δ d = Re { δ d (0) exp( -ikz -iωt ) } and δ h = Re { δ h (0) exp( -ikz -iωt ) } , with amplitudes δ d (0) and δ h (0) that are linearly related to e (0) :

$$\delta \mathbf d ^ { ( 0 ) } = \delta D _ { \text{ax} } \, \mathbf e ^ { ( 0 ) } \,, \quad \delta \mathbf h ^ { ( 0 ) } = \delta H _ { \text{ax} } \, \mathbf e ^ { ( 0 ) } \,.$$

δD ax and δH ax are 3 × 3 matrices, and direct substitution shows that

$$\L x \, \text{auru or nav are } \L o \times \L o \text{ maurices,auru direction shows each} \\ \delta D _ { \text{ax} } & = \frac { 1 } { 8 } \, \frac { \eta ^ { 2 } } { k _ { \text{ax} } ^ { 2 } } \left \{ \frac { B _ { 0 } ^ { ( 0 ) * } \left [ B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] ^ { T } } { 1 - 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } + \frac { B _ { 0 } ^ { ( 0 ) } \left [ B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] ^ { * T } } { 1 + 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } \right \}, \\ \delta H _ { \text{ax} } & = - \frac { 1 } { 8 } \, \frac { \eta ^ { 2 } } { k _ { \text{ax} } ^ { 2 } } \left \{ \frac { \mathbf E _ { 0 } ^ { ( 0 ) * } \left [ B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] ^ { T } } { 1 - 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } + \frac { \mathbf E _ { 0 } ^ { ( 0 ) } \left [ B _ { 0 } ^ { ( 0 ) } + \hat { z } \times \mathbf E _ { 0 } ^ { ( 0 ) } \right ] ^ { * T } } { 1 + 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } \right \}, \\, & \quad \Gamma \, \L \, \text{.}$$

where ω ax = ck ax = m c / ax 2 glyph[planckover2pi1] .

## D. Axionic contribution to the refractive index

Writing Eqs. (67) as linear in e (0) is particularly convenient as it immediately allows us to combine these results with those of Eqs. (24), and then to obtain the modification to Eq. (25). Indeed, it is straightforward to show that the matrix on the left-hand side of Eq. (25) is changed simply by addition of δD ax -Ω z δH ax , where (in accordance with the observations made just after Eq. (23)) Ω z acting on δH ax replaces the E (0) 0 and E (0) glyph[star] 0 of the second of Eqs. (68) by ˆ z × E (0) 0 and ˆ z × E (0) glyph[star] 0 , respectively. Furthermore, as discussed after Eq. (25), we need only consider the xy -projection of the matrix in Eq. (25) when working to first order in δn . The relevant combination of δD ax and δH ax , once projected onto the xy -plane, depends only on the vector B (0) of Eqs. (41), and not on E (0) 0 and B (0) 0 separately. In short, to the matrices on the left-hand sides of Eqs. (43) and (47) must be added the following:

$$( \delta D _ { \text{ax} } & - \Omega _ { z } \, \delta H _ { \text{ax} } ) | _ { x _ { y } } = \frac { 1 } { 8 } \frac { \eta ^ { 2 } } { k _ { \text{ax} } ^ { 2 } } \left \{ \frac { \mathcal { B } ^ { ( 0 ) * } \mathcal { B } ^ { ( 0 ) T } } { 1 - 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } + \frac { \mathcal { B } ^ { ( 0 ) } \mathcal { B } ^ { ( 0 ) * T } } { 1 + 4 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } \right \} \\ & \equiv \frac { 1 } { 4 } \left [ \Delta _ { 2, \text{ax} } \left ( \omega \right ) \left \{ \cos ^ { 2 } \chi \, \mathcal { B } _ { \text{abs} } ^ { ( 0 ) } B _ { \text{abs} } ^ { ( 0 ) T } + \sin ^ { 2 } \chi \, \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } \right \} + i \, \Delta _ { 3, \text{ax} } \left ( \omega \right ) \, \sin \left ( 2 \chi \right ) \, \left \{ \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } B _ { \text{abs} } ^ { ( 0 ) T } - B _ { \text{abs} } ^ { ( 0 ) T } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } \right \} \right ],$$

where in the last line we have used the decomposition of B (0) given in Eqs. (45), and defined

$$\Delta _ { 2, \text{ax} } \left ( \omega \right ) = \frac { \eta ^ { 2 } } { k _ { \text{ax} } ^ { 2 } } \, \frac { 1 } { 1 - \cos ^ { 4 } \frac { \theta } { 2 } \left ( 4 \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } \right ) ^ { 2 } } \,, \quad \Delta _ { 3, \text{ax} } \left ( \omega \right ) = \frac { \eta ^ { 2 } } { k _ { \text{ax} } ^ { 2 } } \, \frac { 2 \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } } { 1 - \cos ^ { 4 } \frac { \theta } { 2 } \left ( 4 \omega _ { 0 } \omega / \omega _ { \text{ax} } ^ { 2 } \right ) ^ { 2 } } \,.$$

The eigenproblem thus becomes:

$$\frac { 1 } { 4 } \left [ \delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) \mathcal { E } _ { \text{abs} } ^ { ( 0 ) } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } + \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) \mathcal { B } _ { \text{abs} } ^ { ( 0 ) } \mathcal { B } _ { \text{abs} } ^ { ( 0 ) T } + \delta _ { 3 } ^ { \prime \prime } \left ( \omega \right ) \mathcal { B } _ { \text{abs} } ^ { ( 0 ) T } \mathcal { E } _ { \text{abs} } ^ { ( 0 ) T } \right ] \mathrm e ^ { ( 0 ) } = \delta n \left ( \omega \right ) \mathrm e ^ { ( 0 ) } \,,$$

where

$$\delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) + \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) & = \delta _ { 1 } + \delta _ { 2 } + \Delta _ { 2, \text{ax} } \left ( \omega \right ) \,, \\ \delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) - \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) & = \left ( \delta _ { 1 } - \delta _ { 2 } - \Delta _ { 2, \text{ax} } \left ( \omega \right ) \right ) \cos \left ( 2 \chi \right ) \,, \\ \delta _ { 3 } ^ { \prime \prime } \left ( \omega \right ) & = \delta _ { 3 } \cos \left ( 2 \chi \right ) + i \, \Delta _ { 3, \text{ax} } \left ( \omega \right ) \sin \left ( 2 \chi \right ) \,.$$

As a matrix equation, this is simply

$$\left [ \begin{array} { c c } \delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) & \delta _ { 3 } ^ { \prime \prime } \left ( \omega \right ) \\ \delta _ { 3 } ^ { \prime \prime * } \left ( \omega \right ) & \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) \end{array} \right ] \left [ \begin{array} { c } e _ { \xi } ^ { \left ( 0 \right ) } \\ e _ { B } ^ { \left ( 0 \right ) } \end{array} \right ] = \frac { \delta n \left ( \omega \right ) } { \frac { 1 } { 4 } \left | \xi _ { \text{abs} } ^ { \left ( 0 \right ) } \right | ^ { 2 } } \left [ \begin{array} { c } e _ { \xi } ^ { \left ( 0 \right ) } \\ e _ { B } ^ { \left ( 0 \right ) } \end{array} \right ].$$

The eigenvalues of the matrix on the left-hand side of this equation are calculated as before, and the refractive index changes δn ± are again given by Eq. (50) with the now frequency-dependent factors δ ′′ ± ( ω ) taking the form

$$\delta _ { \pm } ^ { \prime \prime } \left ( \omega \right ) & = \frac { 1 } { 2 } \left ( \delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) + \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) \pm \sqrt { \left ( \delta _ { 1 } ^ { \prime \prime } \left ( \omega \right ) - \delta _ { 2 } ^ { \prime \prime } \left ( \omega \right ) \right ) ^ { 2 } + \left | 2 \delta _ { 3 } ^ { \prime \prime } \left ( \omega \right ) \right | ^ { 2 } } \right ) \\ & = \frac { 1 } { 2 } \left ( \delta _ { 1 } + \delta _ { 2 } + \Delta _ { 2, \text{ax} } \left ( \omega \right ) \pm \sqrt { \left [ \left ( \delta _ { 1 } - \delta _ { 2 } - \Delta _ { 2, \text{ax} } \left ( \omega \right ) \right ) ^ { 2 } + \left ( 2 \delta _ { 3 } \right ) ^ { 2 } \right ] \cos ^ { 2 } ( 2 \chi ) + \left ( 2 \Delta _ { 3, \text{ax} } \left ( \omega \right ) \right ) ^ { 2 } \sin ^ { 2 } ( 2 \chi ) } \right ) \,.$$

In the EH model with a linearly polarized pump, the only change with respect to previous sections is that δ 2 → δ 2 +∆ 2 ax , ( ω ). This is in agreement with Eq. (13) of [50] (an apparent difference by a factor of 4 arising only because of different definitions of the pump amplitude).

We briefly mention here that the group index n g = n + ω d n/ d ω , though generally complicated, can be fairly easily calculated in the limits of linear and circular polarization. Although n &lt; 1 for some frequencies, we find (at least in these two limits) that n g &gt; 1 for all frequencies, so that relativistic causality is respected (as explained in footnote 8).

## E. Eigenpolarizations

Generally speaking, the matrix on the left-hand side of Eq. (73) is hermitian but not symmetric, thanks to the imaginary contribution to the off-diagonal component δ ′′ 3 ( ω ). Unlike in previous sections, this means that its eigenvectors are generally complex, which in turn means that the eigenpolarizations are elliptically polarized. Two angles are required to describe these eigenpolarizations: an ellipticity angle ψ ∈ -[ π/ 4 , π/ 4] giving the degree of elliptical polarization, and a rotation angle ϕ ∈ ( -π/ 2 , π/ 2] giving the orientation of the major axes of the ellipses traced by e ( ) and t b ( ) with respect to the basis t { ˆ E , ˆ B } . The (normalized) eigenvectors thus form the columns of a unitary matrix U ψ, ϕ ( ) which can be decomposed as follows:

$$U \left ( \psi, \varphi \right ) = E \left ( \psi \right ) R \left ( \varphi \right ) = \left [ \begin{array} { c c } \cos \psi & i \sin \psi \\ i \sin \psi & \cos \psi \end{array} \right ] \left [ \begin{array} { c c } \cos \varphi & - \sin \varphi \\ \sin \varphi & \cos \varphi \end{array} \right ] \,.$$

The matrix on the left-hand side of Eq. (73) can be written as U ψ, ϕ DU ( ) -1 ( ψ, ϕ ), where D is a diagonal matrix whose entries are the eigenvalues δ ′′ + and δ ′′ -of Eq. (74). As before, we take δ ′′ + to be the first diagonal component of D , so that the left column of U corresponds to the polarization with the larger refractive index change δn + . Explicit calculation of this form of the matrix yields the following relations, which determine ψ and ϕ :

$$( \delta _ { 1 } - \delta _ { 2 } - \Delta _ { 2, \text{ax} } ( \omega ) ) \, \cos \left ( 2 \chi \right ) & = \left ( \delta _ { + } ^ { \prime \prime } \left ( \omega \right ) - \delta _ { - } ^ { \prime \prime } \left ( \omega \right ) \right ) \cos \left ( 2 \varphi \left ( \omega \right ) \right ) \cos \left ( 2 \psi \left ( \omega \right ) \right ) \,, \\ - 2 \Delta _ { 3, \text{ax} } ( \omega ) \, \sin \left ( 2 \chi \right ) & = \left ( \delta _ { + } ^ { \prime \prime } \left ( \omega \right ) - \delta _ { - } ^ { \prime \prime } \left ( \omega \right ) \right ) \cos \left ( 2 \varphi \left ( \omega \right ) \right ) \sin \left ( 2 \psi \left ( \omega \right ) \right ) \,, \\ 2 \delta _ { 3 } \cos \left ( 2 \chi \right ) & = \left ( \delta _ { + } ^ { \prime \prime } \left ( \omega \right ) - \delta _ { - } ^ { \prime \prime } \left ( \omega \right ) \right ) \sin \left ( 2 \varphi \left ( \omega \right ) \right ) \,.$$

Whilst in general the solutions ψ ω ( ) and ϕ ω ( ) are quite complicated, they become rather simple in the two polarization limits of the pump. When the pump is linearly polarized ( χ = 0), we find ψ = 0 so that the probe polarizations are also linearly polarized; the only non-triviality is in their rotation with respect to the pump fields, the rotation angle being given by Eqs. (37) with δ 2 → δ 2 +∆ 2 ax , ( ω ). On the other hand, for a circularly polarized pump ( χ = ± π/ 4), we find ψ = ± π/ 4 too, i.e., the eigenpolarizations of the probe are themselves circularly polarized. In the high-frequency regime where ∆ 3 ax , ( ω ) &lt; 0, it is straightforward to show 11 that the larger refractive index change δn + is felt by the probe rotating in the same sense as E ( ) and t B ( ), or equivalently (since in Fig. 2 the probe propagates t into the page) by the probe with opposite handedness to that of the pump.

11 There is a slight complication as there exists a degeneracy in ϕ and ψ (as defined by Eqs. (76)) when χ = ± π/ 4: either ϕ = 0 and ψ = χ , or ϕ = π/ 2 and ψ = -χ . However, plugging these solutions into Eq. (75), we find that switching between them simply amounts to multiplication of the columns of U ψ , ϕ ( ) by ± i , so that the eigenvectors remain unchanged (up to an unimportant phase).

## F. Discussion of key axionic effects

The terms ∆ 2 ax , ( ω ) and ∆ 3 ax , ( ω ) describe two key effects induced by the coupling to axions. Moreover, each is dominant in a particular regime, for we have

$$\frac { \Delta _ { 3, \text{ax} } \left ( \omega \right ) } { \Delta _ { 2, \text{ax} } \left ( \omega \right ) } = 2 \cos ^ { 2 } \frac { \theta } { 2 } \, \frac { \omega _ { 0 } \omega } { \omega _ { \text{ax} } ^ { 2 } } \,.$$

Therefore, whether ω ω/ω 0 2 ax is large or small compared to 1 determines which of the two terms is dominant. (Both can be considered large when ω ω/ω 0 2 ax ≈ 1, but then we are close to resonance, and the analysis performed here will cease to be valid.)

## 1. Renormalization of δ 2

First, the parameter δ 2 entering the local part of the effective Lagrangian is renormalized in a frequency-dependent way to become δ 2 +∆ 2 ax , ( ω ). When ω ω/ω 0 2 ax is very large, it is clear from Eqs. (70) that ∆ 2 ax , ( ω ) → 0. We may thus focus on the opposite limit, ω ω/ω 0 2 ax → 0, where the Compton wavelength of the axion is much smaller than the typical photon wavelengths (or indeed when the background fields are static). In this case, we have ∆ 2 ax , ( ω ) → η /k 2 2 ax , with no residual dependence on frequency. This is the limit in which the effective photon/photon interaction becomes local; η /k 2 2 ax can be absorbed into the definition of δ 2 , and the analysis of Sec. IV carries through as before. We are thus provided with an explicit demonstration of the results of the local effective theory emerging in the correct limit.

For the first experimental probes of NLED, this renormalization is only expected to be significant if it is comparable to the 'bare' EH value, δ (EH) 2 (see Eqs. (9)). In units with glyph[planckover2pi1] = c = 1, this requires η/m ax glyph[greaterorsimilar] √ δ (EH) 2 ≈ 3 × 10 4 GeV -2 . This has not yet been ruled out by current experimental tests of axions but, since values of η larger than about 10 -10 GeV -1 have been essentially excluded (see Fig. 111.1 of [51]), this would require an axion mass smaller than a few times 10 -6 eV. This makes it less likely that the condition ω ω/ω 0 2 ax glyph[greatermuch] 1 will be satisfied, unless the background field is essentially static. For the DeLLight experiment [25], which uses optical wavelengths with glyph[planckover2pi1] ω ∼ 1 eV for both the background and the probe fields, this renormalization of δ 2 is unlikely to be significant.

## 2. Elliptical birefringence

The second key axionic effect is the contribution ∝ ∆ 3 ax , ( ω ) to the off-diagonal terms in Eq. (73). Interestingly, this contribution is imaginary, and leads to probe eigenstates which are elliptically polarized. It is clear that ∆ 3 ax , ( ω ) → 0 when ω ω/ω 0 2 ax → 0, as it should, since we then recover the results of the local effective theory where (as we have seen in Sec. IV) the relevant matrix is purely real. In the opposite limit, we have already seen (in Eq. (77)) that ∆ 3 ax , ( ω ) becomes the dominant signature of the axion coupling, and yet, since it enters Eqs. (72) with a factor of sin (2 χ ), it completely drops out when the pump is linearly polarized. We thus conclude that, when working in the limit ω ω/ω 0 2 ax glyph[greatermuch] 1, a circularly polarized pump is much more efficient than a linearly polarized one at inducing an axionic signature in the refractive index. This result is made all the more interesting by the fact that this axionic signature includes birefringence, since (as seen in Sec. IV C) models with local effective photon/photon interactions show no birefringence at all when the pump is circularly polarized.

Working in the appropriate limit where ω ω/ω 0 2 ax is very large,

$$\Delta _ { 3, \text{ax} } ( \omega ) \rightarrow - \frac { \eta ^ { 2 } c ^ { 2 } } { 8 } \frac { 1 } { \cos ^ { 2 } \frac { \theta } { 2 } \, \omega _ { 0 } \omega } \,,$$

and we see explicitly that this limiting case is independent of the actual value of m ax . Moreover, in the optimized scenario with θ = 0 and χ = ± π/ 2 (i.e., a counter-propagating probe and a circularly polarized pump), we have from Eq. (74) that the birefringence δ + -δ -= 2 | ∆ 3 ax , ( ω ) . | It is appropriate to compare this with δ (EH) 2 -δ (EH) 1 , the birefringence predicted by the EH model (in the absence of axions) for a linearly polarized pump. In units with glyph[planckover2pi1] = c = 1, they are comparable when η/ √ ω ω 0 ∼ 2 √ δ (EH) 2 -δ (EH) 1 ≈ 4 × 10 4 GeV -2 . At optical frequencies (as used in the DeLLight experiment [25]), and given that current exclusion plots indicate η is at most ∼ 10 -10 GeV -1 (see Fig. 111.1 of [51]), this circular birefringence turns out to be at least six orders of magnitude smaller than the linear birefringence of the EH model.

## VI. SUMMARY AND CONCLUSION

Starting from the most general Lagrangian for electromagnetic fields consistent with Lorentz invariance, locality of effective interactions, and weak nonlinearities (so that the lowest nonlinear contributions are sufficient to describe the physics), we have derived an effective medium description for the propagation of weak probe waves in the presence of strong background fields, whose effects on the probe can be incorporated through well-defined electric, magnetic and magnetoelectric susceptibility tensors. This description allows the assignment of a refractive index to the effective medium, though the index typically exhibits anisotropy and birefringence. In the case where the background fields are provided by an intense propagating wave or 'pump', there is a further dependence of the refractive index on the degree of elliptical polarization of the pump. The effects of wave vector direction and polarization turn out to be neatly separated: the misalignment of the wave vectors of pump and probe is equivalent to an overall reduction in the pump intensity by a factor of cos 4 ( θ/ 2), where θ is the tilt angle between the two wave vectors; and the ellipticity angle χ enters into the strength of the birefringence with a factor of cos (2 χ ), being maximum for a linearly polarized pump and vanishing for a circularly polarized one. Finally, factoring out the intensity of the pump allows us to extract the nonlinear index of vacuum, which in the Euler-Heisenberg model derived from QED is typically on the order of 10 -33 cm 2 / W, around 18 orders of magnitude smaller than in nonlinear optical media.

Locality of effective photon/photon interactions was relaxed as a constraint through coupling the electromagnetic field to an axion field of unspecified mass. Generation of an axion field by the interaction between pump and probe, followed by backreaction of the axion field on the probe, yields a contribution to the effective photon/photon coupling which is dispersive when the axion's Compton wavelength is larger than or of the same order as the typical photon wavelength. In the case of a linearly polarized pump, this amounts to a straightforward renormalization of one of the post-Maxwellian parameters entering the Lagrangian, but when the pump is circularly polarized, there exists a residual birefringence that (as mentioned above) would vanish if the effective photon/photon interactions were purely local. Whether the axionic contribution to the birefringence is larger for a linearly or circularly polarized pump depends on the typical photon wavelength (defined via the geometric mean of the pump and probe frequencies): if it is much smaller than the Compton wavelength of the axion, the effect is larger for a circularly polarized pump; conversely, if the typical photon wavelength is much larger than the Compton wavelength of the axion, the effect is larger for a linearly polarized pump.

Concerning the DeLLight experiment [25], whose aim is to detect the deflection of a probe wave by the index variation induced by a tightly focused laser pulse, the results presented here confirm that an effect of this kind should indeed be seen, and any effects due to axions are expected to be subdominant. However, as discussed after Eqs. (15) (particularly in footnote 6), we have shown that only a single pump pulse is required, rather than two counterpropagating pump pulses as proposed in [25]. We conclude that the proposal presented in [25] can be simplified by keeping only the pump pulse which is counter-propagating with respect to the probe. We are currently performing numerical simulations of the simplified DeLLight experiment, whose results will be published in a future work. It is worth noting that the two-pump proposal of [25] would be an interesting case in which four-wave mixing processes could turn out to be important, as the two counter-propagating pump pulses would engender a stationary oscillation in space that could act as a diffraction grating for the probe [52]; this, however, is beyond the scope of the present work.

## ACKNOWLEDGMENTS

I thank Xavier Sarazin and Fran¸ cois Couchot for many fruitful discussions on the DeLLight experiment, which motivated the present work, and for their careful reading of and constructive comments on the manuscript. I thank Renaud Parentani for many discussions on the theoretical side, and I am grateful to Laurent Schoeffel for a very useful exchange on how to include the coupling to axions. This work was mainly done during two short-stay postdocs at LAL in 2018, funded by LAL and LPT (Laboratoire de Physique Th´ eorique). It is continuing to be supported by the French National Research Agency via the grant no. ANR-18-CE31-0005-01.

[1] M. Faraday, Phil. Mag. 28 , 294 (1846).

[2] N. Bloembergen, Rev. Mod. Phys. 54 , 685 (1982).

[3] N. Bloembergen, Nonlinear Optics , 4th ed. (World Scientific, 1996).

[4] Y. R. Shen, The Principles of Nonlinear Optics , 1st ed. (Wiley, 1984).

[5] R. Boyd, Nonlinear Optics , 3rd ed. (Academic Press, 2008).

- [6] G. P. Agrawal, Nonlinear Fiber Optics , 4th ed. (Academic Press, 2007).
- [7] J. Kerr, Phil. Mag. 50 , 337 (1875).
- [8] H. Euler and B. Kockel, Naturwissenschaften 23 , 246 (1935).
- [9] H. Euler, Ann. Phys. (Leipzig) 26 , 398 (1936).
- [10] W. Heisenberg and H. Euler, Z. Phys. 98 , 714 (1936), arXiv:physics/0605038 [physics.hist-ph].
- [11] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38 , 1440 (1977).
- [12] R. D. Peccei and H. R. Quinn, Phys. Rev. D 16 , 1791 (1977).
- [13] F. Wilczek, Phys. Rev. Lett. 40 , 279 (1978).
- [14] S. Weinberg, Phys. Rev. Lett. 40 , 223 (1978).
- [15] M. Born and L. Infeld, Proc. Roy. Soc. A 144 , 425 (1934).
- [16] R. Battesti and C. Rizzo, Rep. Prog. Phys. 76 , 016401 (2013), arXiv:1211.1933 [physics.optics].
- [17] M. Fouch´ e, R. Battesti, and C. Rizzo, Phys. Rev. D 93 , 093020 (2016), arXiv:1605.04102 [physics.optics].
- [18] M. Aaboud et al. (ATLAS Collaboration), Nature Physics 13 , 852 (2017).
- [19] R. Battesti, B. Pinto Da Souza, S. Batut, C. Robilliard, G. Bailly, C. Michel, M. Nardone, L. Pinard, O. Portugall, G. Tr´ enec, J.-M. Mackowski, G. L. J. A. Rikken, J. Vigu´, e and C. Rizzo, Eur. Phys. J. D 46 , 323 (2008), arXiv:0710.1703 [physics.optics].
- [20] A. Cad` ene, P. Berceau, M. Fouch´ e, R. Battesti, and C. Rizzo, Eur. Phys. J. D 68 , 16 (2014), arXiv:1302.5389 [physics.optics].
- [21] G. Zavattini, U. Gastaldi, R. Pengo, G. Ruoso, F. Della Valle, and E. Milotti, Int. J. Mod. Phys. A 27 , 1260017 (2012), arXiv:1201.2309 [hep-ex].
- [22] F. Della Valle, A. Ejlli, U. Gastaldi, G. Messineo, E. Milotti, R. Pengo, G. Ruoso, and G. Zavattini, Eur. Phys. J. C 76 , 24 (2016), arXiv:1510.08052 [physics.optics].
- [23] R. V. Jones, Nature 186 , 706 (1960).
- [24] R. V. Jones, Proc. R. Soc. Lond. A 260 , 47 (1961).
- [25] X. Sarazin, F. Couchot, A. Djannati-Ata¨ ı, O. Guilbaud, S. Kazamias, M. Pittman, and M. Urban, Eur. Phys. J. D 70 , 13 (2016), arXiv:1507.07959 [physics.optics].
- [26] J. Pleba´ski, n Lectures on Nonlinear Electrodynamics (NORDITA, Copenhagen, 1970).
- [27] G. Boillat, J. Math. Phys. 11 , 941 (1970).
- [28] Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. D 2 , 2341 (1970).
- [29] V. I. Denisov and I. P. Denisova, Dokl. Phys. 46 , 377 (2001).
- [30] P. Sikivie, Phys. Rev. Lett. 51 , 1415 (1983).
- [31] L. Maiani, R. Petronzio, and E. Zavattini, Phys. Lett. B 175 , 359 (1986).
- [32] M. Gasperini, Phys. Rev. Lett. 59 , 396 (1987).
- [33] K. van Bibber, N. R. Dagdeviren, S. E. Koonin, A. K. Kerman, and H. N. Nelson, Phys. Rev. Lett. 59 , 759 (1987).
- [34] G. Raffelt and L. Stodolsky, Phys. Rev. D 37 , 1237 (1988).
- [35] E. Brezin and C. Itzykson, Phys. Rev. D 3 , 618 (1971).
- [36] E. B. Aleksandrov, A. A. Ansel'm, and A. N. Moskalev, Zh. Eksp. Teor. Fiz. ´ 89 , 1181 (1985).
- ´
- [37] N. N. Rozanov, Zh. Eksp. Teor. Fiz. 103 , 1996 (1993).
- [38] A. Ferrando, H. Michinel, M. Seco, and D. Tommasini, Phys. Rev. Lett. 99 , 150404 (2007), arXiv:0703215 [physics].
- [39] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields , 4th ed. (Butterworth-Heinemann, Oxford, 1980).
- [40] T.-D. Lee and C.-N. Yang, Phys. Rev. 104 , 254 (1956).
- [41] C. S. Wu, E. Ambler, R. W. Hayward, D. D. Hoppes, and R. P. Hudson, Phys. Rev. 105 , 1413 (1957).
- [42] F. Moulin and D. Bernard, Opt. Commun. 164 , 137 (1999), arXiv:0203069 [physics].
- [43] S. L. Adler, Ann. Phys. (NY) 67 , 599 (1971).
- [44] G. L. J. A. Rikken and C. Rizzo, Phys. Rev. A 63 , 012107 (2000).
- [45] G. L. J. A. Rikken and C. Rizzo, Phys. Rev. A 67 , 015801 (2003).
- [46] W. F. Brown, Jr., R. M. Hornreich, and S. Shtrikman, Phys. Rev. 168 , 574 (1968).
- [47] P. Arias, D. Cadamuro, M. Goodsell, J. Jaeckel, J. Redondo, and A. Ringwald, J. Cosmol. Astropart. Phys. 6 , 13 (2012), arXiv:1201.5902 [hep-ph].
- [48] G. G. Raffelt, 'Astrophysical axion bounds,' in Axions , edited by M. Kuster, G. Raffelt, and B. Beltr´n (Springer Berlin a Heidelberg, Berlin, Heidelberg, 2008) pp. 51-71, arXiv:hep-ph/0611350 [hep-ph].
- [49] D. J. E. Marsh, Phys. Rep. 643 , 1 (2016), arXiv:1510.07633 [astro-ph.CO].
- [50] C. Baldenegro, S. Hassani, C. Royon, and L. Schoeffel, Phys. Lett. B 795 , 339 (2019).
- [51] M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98 , 030001 (2018).
- [52] A. Di Piazza, K. Z. Hatsagortsyan, and C. H. Keitel, Phys. Rev. Lett. 97 , 083603 (2006).
Scott Robertson
Optical Kerr effect in vacuum
From an effective field theory of electromagnetism in vacuum including all
lowest-order nonlinear terms consistent with Lorentz invariance and locality of
photon/photon interactions, we derive an effective medium description of strong
background fields as regards their influence on a weak probe. We mainly
consider as background a pump beam with well-defined wave vector and
polarization. This leads us to define a nonlinear index of vacuum which, in the
Euler-Heisenberg model derived from QED, has an optimal value of $1.555 \times
10^{-33} \, {\rm cm}^{2}/{\rm W}$ for a linearly polarized pump as seen by a
counter-propagating, orthogonally polarized probe. We further generalize the
model to include coupling to an axion field. In the limit where the axion mass
is much smaller than the typical photon energy, this yields dispersive
corrections, and the axionic signature is found to be greatly enhanced for a
circularly polarized pump as compared to a linearly polarized one. The
formalism here presented points to a simplification of the DeLLight experiment
[X. Sarazin et al., Eur. Phys. J. D 70, 13 (2016)] aiming to measure the
deflection of a probe by a tightly focused laser pulse.