Light-by-Light Scattering Constraint on Born-Infeld Theory

TL;DR

The paper uses ATLAS light-by-light scattering data to set new lower bounds on the mass scale of Born-Infeld extensions of QED and the Standard Model, significantly tightening previous constraints.

Contribution

It provides the first direct experimental constraints on the Born-Infeld theory parameters using LHC data, improving previous bounds by a large margin.

Findings

Mass scale of Born-Infeld QED > 100 GeV

Mass scale of nonlinear U(1)_Y > 90 GeV

Lower limit on magnetic monopole mass > 11 TeV

Abstract

The recent measurement by ATLAS of light-by-light scattering in LHC Pb-Pb collisions is the first direct evidence for this basic process. We find that it requires the mass scale of a nonlinear Born-Infeld extension of QED to be $\gtrsim 100$~GeV, a much stronger constraint than those derived previously. In the case of a Born-Infeld extension of the Standard Model in which the U(1)$_{\rm Y}$ hypercharge gauge symmetry is realized nonlinearly, the limit on the corresponding mass scale is $\gtrsim 90$~GeV, which in turn imposes a lower limit of $\gtrsim 11$~TeV on the magnetic monopole mass in such a U(1)$_{\rm Y}$ Born-Infeld theory.

Figures

Equations

$${\mathcal{L}}_{\mathrm{QED}}=-\frac{1}{4}{F}_{\mu \nu }{F}^{\mu \nu }\to$$

$$\displaystyle{\cal L}_{\rm BI}\;=\;\beta^{2}\Big{(}1-\sqrt{1+\frac{1}{2\beta^{2}}F_{\mu\nu}F^{\mu\nu}-\frac{1}{16\beta^{4}}(F_{\mu\nu}\tilde{F}^{\mu\nu})^{2}}\;\Big{)},$$

$$M\gtrsim 100\mathrm{GeV},$$

$$M_Y=\cos {\theta }_{W}M\gtrsim 90\mathrm{GeV},$$

$${\sigma }_{\text{BI}}(\gamma \gamma \to \gamma \gamma )=\frac{1}{2}\int d\Omega \frac{d{\sigma }_{\text{BI}}}{d\Omega }=\frac{7}{1280\pi }\frac{m_{\gamma \gamma }^6}{{\beta }^4},$$

$$\frac{d{\sigma }_{\text{BI}}}{d\Omega }=\frac{1}{4096{\pi }^2}\frac{m_{\gamma \gamma }^6}{{\beta }^4}{\left(3+\cos \theta \right)}^2.$$

$${\sigma }_{\text{excl.}}={\int }_{{\tau }_0}^{1}d\tau \frac{dL}{d\tau }{\sigma }_{\gamma \gamma \to \gamma \gamma }(\tau ).$$

$$\frac{dL}{d\tau }={\int }_{\tau }^{1}d{x}_{1}d{x}_{2}f(x_1)f(x_2)\delta (\tau -x_1{x}_2),$$

Full text

Over 80 years ago, soon after Dirac proposed his relativistic theory of the electron Dirac28 and his interpretation of ‘hole’ states as positrons Dirac31 , Halpern Halpern33 in 1933 and Heisenberg Heisenberg34 in 1934 realized that quantum effects would induce light-by-light scattering, which was first calculated in the low-frequency limit by Euler and Kockel EulerKockel35 in 1935. Subsequently, Heisenberg and Euler HeisenbergEuler36 derived in 1936 a more general expression for the quantum nonlinearities in the Lagrangian of Quantum Electrodynamics (QED), and a complete calculation of light-by-light scattering in QED was published by Karplus and Neuman KarplusNeuman51 in 1951. However, measurement of light-by-light scattering has remained elusive until very recently. In 2013 d’Enterria and Silveira dES13 proposed looking for light-by-light scattering in ultraperipheral heavy-ion collisions at the LHC, and evidence for this process was recently presented by the ATLAS Collaboration ATLAS17 , at a level consistent with the QED predictions in dES13 and K-GLS16 .

In parallel with the early work on light-by-light scattering in QED, and motivated by a ‘unitarian’ idea that there should be an upper limit on the strength of the electromagnetic field, Born and Infeld BornInfeld34 proposed in 1934 a conceptually distinct nonlinear modification of the Lagrangian of QED:

[TABLE]

where $\beta$ is an a priori unknown parameter with the dimension of [Mass]2 that we write as $\beta\equiv M^{2}$, and $\tilde{F}_{\mu\nu}$ is the dual of the field strength tensor ${F_{\mu\nu}}$. Interest in Born-Infeld theory was revived in 1985 when Fradkin and Tseytlin FT85 discovered that it appears when an Abelian vector field in four dimensions is coupled to an open string, as occurs in models inspired by M theory in which particles are localized on lower-dimensional ‘branes’ separated by a distance $\simeq 1/\sqrt{\beta}=1/M$ in some extra dimension 111Remarkably, the maximum field strength is related to the fact that the brane velocity is limited by the velocity of light Bachas95 , confirming the insight of Born and Infeld BornInfeld34 .. Depending on the specific brane scenario considered, $M$ might have any value between a few hundred GeV and the Planck scale $\sim 10^{19}$ GeV. For the purposes of this paper, we consider only the relevant terms of fourth order in the gauge field strengths in (1).

Until now, there has been no strong lower limit on the Born-Infeld scale $\beta$ or, equivalently, the brane mass scale $M$ and the brane separation $1/M$. A constraint corresponding to $M\gtrsim 100$ MeV was derived in RSG from electronic and muonic atom spectra, though the derivation has been questioned in CK . Measurements of photon splitting in atomic fields split were considered in Davilaetal , where it was concluded that they provided no limit on the Born-Infeld scale and it was suggested that measurements of the surface magnetic field of neutron stars neutronstar might be sensitive to $M=\sqrt{\beta}\sim 1.4\times 10^{-5}$ GeV. More recently, measurements of nonlinearities in light by the PVLAS Collaboration PVLAS14 are somewhat more sensitive to the individual nonlinear terms in (1), but are insensitive to the particular combination appearing in the Born-Infeld theory, as discussed in Fouche16 where more references can be found.

Here we show that the agreement of the recent ATLAS measurement of light-by-light scattering with the standard QED prediction provides the first limit on $M$ in the multi-GeV range, excluding a significant range extending to

[TABLE]

entering the range of interest to brane theories. This limit is obtained under quite conservative assumptions, and plausible stronger assumptions would strengthen our lower bound to $M\gtrsim 200$ GeV.

One may also consider a Born-Infeld extension of the Standard Model in which the hypercharge U(1)Y gauge symmetry is realised non-linearly, in which case the limit (2) is relaxed to

[TABLE]

where we have used $B_{Y}^{\mu}={\rm cos}\theta_{W}\,A_{\rm EM}^{\mu}-{\rm sin}\theta_{W}\,Z^{\mu}$ and $\sin^{2}{\theta_{W}}\simeq 0.23$, with $\theta_{W}$ the weak mixing angle. As a corollary of this lower limit on the U(1)Y brane scale, we recall that Arunasalam and Kobakhidze recently pointed out ArunasalamKobakhidze17 that the Standard Model modified by a Born-Infeld U(1)Y theory has a finite-energy electroweak monopole cho-maison ; ChoKimYoon13 solution ${\cal M}$, whose mass they estimated as $M_{\cal M}\simeq 4~{}{\rm TeV}+72.8\,M_{Y}$. Such a monopole is less constrained by Higgs measurements than electroweak monopoles in other extensions of the Standard Model EMY , and hence of interest for potential detection by the ATLAS atlasmono , CMS and MoEDAL experiments at the LHC MoEDAL . However, our lower limit $M_{Y}\gtrsim 90$ GeV (2) corresponds to a 95% CL lower limit on the mass of this monopole $M_{\cal M}\gtrsim 11$ TeV, excluding its production at the LHC.

Following the suggestion of dES13  222These authors also considered the observability of light-by-light scattering in pp collisions at the LHC, but concluded that they were less promising than Pb-Pb collisions., we consider ultra-peripheral heavy-ion collisions in which the nuclei scatter quasi-elastically via photon exchange: Pb + Pb ($\gamma\gamma$)$\to$Pb*(∗)* + Pb*(∗)*+ X, as depicted in Fig. 1, effectively acting via the equivalent photon approximation (EPA) EPA as a photon-photon collider. The EPA allows the electromagnetic field surrounding a highly-relativistic charged particle to be treated as equivalent to a flux of on-shell photons. Since the photon flux is proportional to $Z^{2}$ for each nucleus, the coherent enhancement in the exclusive $\gamma\gamma$ cross-section scales as $Z^{4}$, where $Z=82$ for the lead (Pb) ions used at the LHC. This is why heavy-ion collisions have an advantage over proton-proton or proton-lead collisions for probing physics in electromagnetic processes dES13 . Photon fusion in ultra-peripheral heavy-ion collisions has been suggested as a way of detecting the Higgs boson Higgsion ; DEZ and, more recently, the possibility of constraining new physics beyond the Standard Model (BSM) in this process was studied in Fichet ; Knapenetal .

As already mentioned, the possibility of directly observing light-by-light scattering at the LHC was proposed in dES13 , and this long-standing prediction of QED was finally measured earlier this year with $4.4\sigma$ significance by the ATLAS Collaboration ATLAS17 at a level in good agreement with calculations in dES13 ; K-GLS16 . The compatibility with the Standard Model constrains any possible contributions from BSM physics. Born-Infeld theory is particularly interesting in this regard, as constraints from low-energy optical and atomic experiments have yet to reach the sensitivity of interest for measuring light-by-light scattering PVLAS14 ; Fouche16 .

The leading-order cross-section for unpolarised light-by-light scattering in Born-Infeld theory in the $\gamma\gamma$ centre-of-mass frame is given by Davilaetal ; TurkRebhan :

[TABLE]

where $m_{\gamma\gamma}$ is the diphoton invariant mass and the differential cross-section is

[TABLE]

We recall that the parameter $\beta=M^{2}$ enters as a dimensionful parameter in the Born-Infeld theory of non-linear QED defined by the Lagrangian (1). If this originates from a Born-Infeld theory of hypercharge then the corresponding mass scale is $M_{Y}=\cos{\theta_{W}}M$.

We plot in Fig. 2 the angular distributions as functions of $\cos\theta$ in the centre-of-mass frame (where $\theta$ is the polar angle) for the leading-order differential cross-sections in both Born-Infeld theory and QED (with arbitrary normalisations), as solid blue and dashed red lines, respectively. We see that the Born-Infeld distribution is less forward peaked than that for QED. For the latter, we used the leading-order amplitudes for the quark and lepton box loops in the ultra-relativistic limit from Bernetal , omitting the percent-level effects of higher-order QCD and QED corrections, as well the $W^{\pm}$ contribution that is negligible for typical diphoton centre-of-mass masses at the LHC.

The total exclusive diphoton cross-section from Pb+Pb collisions is obtained by convoluting the $\gamma\gamma\to\gamma\gamma$ cross-section with a luminosity function $dL/d\tau$ CahnJackson ,

[TABLE]

We have introduced here a dimensionless measure of the diphoton invariant mass, $\tau\equiv m_{\gamma\gamma}^{2}/s_{NN}$, where $\sqrt{s_{NN}}=5.02$ TeV is the centre-of-mass energy per nucleon pair in the ATLAS measurement. The luminosity function, derived for example in CahnJackson , can be written as an integral over the number distribution of photons carrying a fraction $x$ of the total Pb momentum:

[TABLE]

where the distribution function $f(x)$ depends on a nuclear form factor. We follow CahnJackson in adopting the form factor proposed in DEZ , while noting that variations in the choice leads to $\sim 20\%$ uncertainties in the final cross-sections dES13 . A contribution with a non-factorisable distribution function should also be subtracted to account for the exclusion of nuclear overlaps, but this is not a significant effect for the relevant kinematic range, causing a difference within the $20\%$ uncertainty Knapenetal from the photon luminosity evaluated numerically using the STARlight code starlight . For $\sqrt{s_{NN}}=5.5$ TeV and $m_{\gamma\gamma}>5$ GeV we obtain a QED cross-section of $\sigma^{\text{QED}}_{\text{excl.}}=385\pm 77$ nb, in good agreement with dES13 . The ATLAS measurement is performed at $\sqrt{s_{NN}}=5.02$ TeV and for $m_{\gamma\gamma}>6$ GeV, for which we find $\sigma^{\text{QED}}_{\text{excl.}}=220\pm 44$ nb.

This total $\gamma\gamma\to\gamma\gamma$ cross-section is reduced by the fiducial cuts of the ATLAS analysis, which restrict the phase space to a photon pseudorapidity region $|\eta|<2.4$, and require photon transverse energies $E_{T}>3$ GeV and the diphoton system to have an invariant mass $m_{\gamma\gamma}>6$ GeV with a transverse momentum $p_{T}^{\gamma\gamma}<2$ GeV and an acoplanarity $A\text{co}=1-\Delta\phi/\pi<0.01$. We simulate the event selection using Monte-Carlo sampling, implementing the cuts with a 15% Gaussian smearing in the photon transverse energy resolution at low energies and 0.7% at higher energies ATLAS17 ; ATLASphoton above 100 GeV. Since the differential cross-section does not depend on $\phi$ we implement the acoplanarity cut as a fixed 85% efficiency in the number of signal events after the $p_{T}^{\gamma\gamma}$ selection, following the ATLAS analysis ATLAS17 . The total reduction in yield for the QED case is a factor $\epsilon\sim 0.30$, which results in a fiducial cross-section $\sigma^{\text{QED}}_{\text{fid.}}=53\pm 11$ nb for $\sqrt{s}_{NN}=5.02$ TeV, in good agreement with the two predictions of 45 and 49 nb quoted by ATLAS ATLAS17 .

Following this validation for the QED case, we repeat the procedure for the Born-Infeld cross-section. Since the Born-Infeld $\gamma\gamma\to\gamma\gamma$ cross-section grows with energy, the dominant contribution to the cross-section comes from the $\tau\lesssim 0.2$ part of the integral, compared with $\tau\lesssim 10^{-4}$ for the QED case. We show in Fig. 3 the distributions of the $\sigma(\gamma\gamma\to\gamma\gamma)$ cross-section multiplied by the photon flux luminosity factor – normalised by the total exclusive cross-section – as functions of the invariant diphoton mass distribution, for the QED case in the left panel and in Born-Infeld theory with $M=\sqrt{\beta}=200$ GeV in the right panel.

We see that the invariant-mass distribution in the Born-Infeld case extends to $m_{\gamma\gamma}>M$, where the validity of the tree-level Born-Infeld Lagrangian may be questioned because the Taylor expansion of the square root in the non-polynomial Born-Infeld Lagrangian (1) could break down. With this in mind, we use two approaches to place plausible limits on $M=\sqrt{\beta}$. In the first and most conservative method we consider $\gamma\gamma$ scattering only for $m_{\gamma\gamma}\leq M$, while in the second approach we integrate the $\gamma\gamma$ cross-section (4) up to the diphoton invariant mass where the unitarity limit $\sigma_{\text{BI}}\sim 1/m_{\gamma\gamma}^{2}$ is saturated, beyond which we assume that the cross-section saturates the unitarity limit and falls as $\sim 1/m_{\gamma\gamma}^{2}$.

We find fiducial efficiencies for the cut-off and unitarization approaches to be $\epsilon\sim 0.39$ and $0.14$, respectively. Whilst the $E_{T}$ and $\eta$ cuts have much less effect than for QED, as expected from the difference in the angular distributions visible in Fig. 2, the larger invariant masses appearing in the Born-Infeld case are much more affected by the $p_{T}^{\gamma\gamma}$ requirement.

Our calculations of the corresponding U(1)EM Born-Infeld fiducial cross-sections are plotted in the left panel of Fig. 4 as a function of $M=\sqrt{\beta}$: the green curve is for the more conservative cut-off approach, and the blue curve assumes that unitarity is saturated. These calculations are confronted with the ATLAS measurement of $\sigma_{\text{fid.}}=70\pm 24\text{ (stat.)}\pm 17\text{ (sys.)}$ nb ATLAS17 , assuming that these errors are Gaussian and adding them in quadrature with a theory uncertainty of $\pm 10$ nb. We perform a $\chi^{2}$ fit to obtain the 95% CL upper limit on a Born-Infeld signal additional to the 49 nb Standard Model prediction 333We neglect possible interference effects that are expected to be small due to the different invariant mass and angular distributions involved.. This corresponds to the excluded range shaded in pink above $\sigma_{\text{fid.}}^{95\%\text{CL}}\sim 65$ nb in the left panel of Fig. 4, which translates into the limit $M=\sqrt{\beta}\gtrsim 100\,(190)$ GeV in the cut-off (unitarized) approach, as indicated by the green (blue) vertical dashed line in Fig. 4, respectively.

These limits could be strengthened further by considering the $m_{\gamma\gamma}$ distribution shown in Fig. 3(b) of ATLAS17 , where we see that all the observed events had $m_{\gamma\gamma}<25$ GeV, in line with expectations in QED, whereas in the Born-Infeld theory most events would have $m_{\gamma\gamma}>25$ GeV. Calculating a ratio of the total exclusive cross-section of QED for $m_{\gamma\gamma}>6$ GeV and $>25$ GeV as $\sigma_{\text{excl.}}^{m_{\gamma\gamma}>25\,\text{GeV}}/\sigma_{\text{excl.}}^{m_{\gamma\gamma}>6\,\text{GeV}}\sim 0.02$, we estimate a 95% CL upper limit of $\sim 2$ nb for $m_{\gamma\gamma}>25$ GeV. The corresponding exclusion plot is shown in the right panel of Fig. 4, where we see a stronger limit $M=\sqrt{\beta}\gtrsim 210\,(330)$ GeV in the cut-off (unitarized) approach with the same colour coding as previously.

Our lower limit on the QED Born-Infeld scale $M=\sqrt{\beta}\gtrsim 100$ GeV is at least 3 orders of magnitude stronger than the lower limits on $M=\sqrt{\beta}$ obtained from previous measurements of nonlinearities in light RSG ; CK ; split ; Davilaetal ; PVLAS14 ; Fouche16 . Because of the kinematic cuts made in the ATLAS analysis, our limit does not apply to a range of values of $M\lesssim 10$ GeV for which the nonlinearities in (1) should be taken into account. However, our limit is the first to approach the range of potential interest for string/M theory constructions, since models with (stacks of) branes separated by distances $1/M:M={\cal O}(1)$ TeV have been proposed in that context Tseytlin . Our analysis could clearly be refined with more sophisticated detector simulations and the uncertainties reduced. However, in view of the strong power-law dependence of the Born-Infeld cross-section on $M=\sqrt{\beta}$ visible in (4), the scope for significant improvement in our constraint is limited unless experiments can probe substantially larger $m_{\gamma\gamma}$ ranges. In this regard, it would be interesting to explore the sensitivities of high-energy $e^{+}e^{-}$ machines considered as $\gamma\gamma$ colliders.

As mentioned in the Introduction, Arunasalam and Kobakhidze have recently pointed out ArunasalamKobakhidze17 that the Standard Model modified by a Born-Infeld theory of the hypercharge U(1)Y contains a finite-energy monopole solution with mass $M_{\cal M}=E_{0}+E_{1}$, where $E_{0}$ is the contribution associated with the Born-Infeld U(1)Y hypercharge, and $E_{1}$ is sssociated with the remainder of the Lagrangian. Arunasalam and Kobakhidze have estimated ArunasalamKobakhidze17 that $E_{0}\simeq 72.8\,M_{Y}$, where $M_{Y}=\cos{\theta_{W}}M$, and Cho, Kim and Yoon had previously estimated ChoKimYoon13 that $E_{1}\simeq 4$ TeV 444Both these estimates are at the classical level, and quantum corrections have yet to be explored.. Combining these calculations and using our lower limit $M\gtrsim 100$ GeV (2), we obtain a lower limit $M_{\cal M}\gtrsim 11$ TeV on the U(1)Y Born-Infeld monopole mass 555For completeness, we recall that it was argued in ArunasalamKobakhidze17 that nucleosynthesis constraints on the abundance of relic monopoles require $M_{\cal M}\lesssim 23$ TeV.. Unfortunately, this is beyond the reach of MoEDAL MoEDAL or any other experiment at the LHC atlasmono , but could lie within reach of a similar experiment at any future 100-TeV $pp$ collider FCC-hh , or of a cosmic ray experiment.

In this paper we have restricted our attention to possible Born-Infeld modifications of $U(1)$ gauge factors and their constraints in light-by-light scattering only. We plan to examine in the future the experimental constraints from measurements at the LHC on possible Born-Infeld extensions of the $SU(3)_{C}$ and $SU(2)_{L}$ gauge symmetries of the Standard Model.

Acknowledgements

The work of JE and NEM was supported partly by the STFC Grant ST/L000326/1. The work of TY was supported by a Junior Research Fellowship from Gonville and Caius College, Cambridge. We thank Vasiliki Mitsou for drawing our attention to ArunasalamKobakhidze17 , and her and Albert De Roeck, Igor Ostrovskiy and Jim Pinfold of the MoEDAL Collaboration for their interest and relevant discussions. TY is grateful for the hospitality of King’s College London where part of this work was completed.