The role of meson exchanges in light-by-light scattering

TL;DR

This paper investigates how meson exchange mechanisms, including s-, t-, and u-channel contributions with off-shell corrections, influence light-by-light scattering cross sections, highlighting potential interference effects relevant for collider experiments.

Contribution

It introduces the first comprehensive analysis of t- and u-channel meson exchanges with off-shell effects in light-by-light scattering, expanding beyond previous s-channel focused studies.

Findings

Meson exchanges can significantly alter cross sections depending on vertex form factors.

Interference effects among meson groups and with fermion-box amplitudes are notable.

Potential relevance for interpreting recent ultraperipheral heavy-ion collision measurements.

Abstract

We discuss the role of meson exchange mechanisms in $\gamma \gamma \to \gamma \gamma$ scattering. Several pseudoscalar ($\pi^0$, $\eta$, $\eta'(958)$, $\eta_c(1S)$, $\eta_c(2S)$), scalar ($f_{0}(500)$, $f_0(980)$, $a_0(980)$, $f_0(1370)$, $\chi_{c0}(1P)$) and tensor ($f_2(1270)$, $a_2(1320)$, $f_2'(1525)$, $f_2(1565)$, $a_2(1700)$) mesons are taken into account. We consider not only $s$-channel but also for the first time $t$- and $u$-channel meson exchange amplitudes corrected for off-shell effects including vertex form factors. We find that, depending on not well known vertex form factors, the meson exchange amplitudes interfere among themselves and could interfere with fermion-box amplitudes and modify the resulting cross sections. The meson contributions are shown as a function of collision energy as well as angular distributions are presented. Interesting interference effects separately for light pseudoscalar, scalar and tensor meson groups are discussed. The meson exchange contributions may be potentially important in the context of a measurement performed recently in ultraperipheral collisions of heavy ions by the ATLAS collaboration. The light-by-light interactions could be studied in future in electron-positron collisions by the Belle II at SuperKEKB accelerator.

Tables

Table 1. Table 1: A list of mesons that contribute to the γ ​ γ → γ ​ γ → 𝛾 𝛾 𝛾 𝛾 \gamma\gamma\to\gamma\gamma scattering. The meson masses, their total widths Γ M subscript Γ 𝑀 \Gamma_{M} and branching fractions, ℬ ​ ( M → γ ​ γ ) = Γ ​ ( M → γ ​ γ ) / Γ M ℬ → 𝑀 𝛾 𝛾 Γ → 𝑀 𝛾 𝛾 subscript Γ 𝑀 {\cal B}(M\to\gamma\gamma)=\Gamma(M\to\gamma\gamma)/\Gamma_{M} , are taken from PDG  PDG .

Meson	$m_M$ (MeV)	${\Gamma }_M$ (MeV)	$ℬ\left(M\to \gamma \gamma \right)$	$\Gamma \left(M\to \gamma \gamma \right)$ (keV)
${\pi }^0$	$134.9766\pm 0.0006$		$\left(98.823\pm 0.034\right)\times {10}^{-2}$	$7.8\times {10}^{-3}$
$\eta$	$547.862\pm 0.017$	$\left(1.31\pm 0.05\right)\times {10}^{-3}$	$\left(39.41\pm 0.20\right)\times {10}^{-2}$	0.516
${\eta }^{\prime }\left(958\right)$	$957.78\pm 0.06$	$0.197\pm 0.009$	$\left(2.21\pm 0.08\right)\times {10}^{-2}$	4.354
${\eta }_c\left(1S\right)$	$2983.4\pm 0.5$	$31.8\pm 0.8$	$\left(1.59\pm 0.13\right)\times {10}^{-4}$	5.056
${\eta }_c\left(2S\right)$	$3639.2\pm 1.2$	${11.3}_{-2.9}^{+3.2}$	$\left(1.9\pm 1.3\right)\times {10}^{-4}$	2.147
$f_0\left(500\right)$	$400-550$	$400-700$	$\mathrm{seen}$	$2.05\pm 0.21$ Dai_Pennington:2014
$f_0\left(980\right)$	$990\pm 20$	$10-100$	$\mathrm{seen}$	$0.32\pm 0.05$ Dai_Pennington:2014
$a_0\left(980\right)$	$980\pm 20$	$50-100$	$\mathrm{seen}$	$0.30\pm 0.10$ Amsler98
$f_0\left(1370\right)$	$1200-1500$	$200-500$	$\mathrm{seen}$	$4.0\pm 1.9$ Dai_Pennington:2014
${\chi }_{c0}\left(1P\right)$	$3414.75\pm 0.31$	$10.5\pm 0.6$	$\left(2.23\pm 0.13\right)\times {10}^{-4}$	2.342
$f_2\left(1270\right)$	$1275.5\pm 0.8$	${186.7}_{-2.5}^{+2.2}$	$\left(1.42\pm 0.24\right)\times {10}^{-5}$	2.651; $3.14\pm 0.20$ Pennington:2008xd
$a_2\left(1320\right)$	${1318.3}_{-0.6}^{+0.5}$	$107\pm 5$	$\left(9.4\pm 0.7\right)\times {10}^{-6}$	$1.00\pm 0.06$
$f_2^{\prime }\left(1525\right)$	$1525\pm 5$	${73}_{-5}^{+6}$	$\left(1.10\pm 0.14\right)\times {10}^{-6}$	$0.081\pm 0.009$; $0.13\pm 0.03$ Shchegelsky
$f_2\left(1565\right)$	$1562\pm 13$	$134\pm 8$		$0.70\pm 0.14$ Shchegelsky
$a_2\left(1700\right)$	$1732\pm 16$	$194\pm 40$		$0.30\pm 0.05$ Shchegelsky
$f_4\left(2050\right)$	$2018\pm 11$	$237\pm 18$		0.7 KS2013

Equations

$$\gamma (p_1,{\lambda }_1)+\gamma (p_2,{\lambda }_2)\to \gamma (p_3,{\lambda }_3)+\gamma (p_4,{\lambda }_4)$$

$$\frac{d\sigma }{d\cos \theta }=\frac{1}{32\pi s}\frac{1}{2}\frac{1}{4}\sum _{{\lambda }_1,{\lambda }_2,{\lambda }_3,{\lambda }_4}|{\mathcal{M}}_{{\lambda }_1{\lambda }_2\to {\lambda }_3{\lambda }_4}{|}^2,$$

$$s=(p_1+p_2{)}^2=(p_3+p_4{)}^2,$$

$$t=(p_1-p_3{)}^2=(p_2-p_4{)}^2,$$

$$u=(p_2-p_3{)}^2=(p_1-p_4{)}^2,$$

$$p_s=p_1+p_2=p_3+p_4,$$

$$p_t=p_2-p_4=p_3-p_1,$$

$$p_u=p_1-p_4=p_3-p_2,$$

$$p_s^2=s,p_t^2=t,p_u^2=u.$$

$${\mathcal{M}}_{{\lambda }_1{\lambda }_2\to {\lambda }_3{\lambda }_4}=$$

$$i{\mathcal{M}}_{{\lambda }_1{\lambda }_2\to {\lambda }_3{\lambda }_4}^{(M_{PS})}$$

$$i{\Delta }^{(M_{PS})}(k)=\frac{i}{k^2-m_{M_{PS}}^2+i{m}_{M_{PS}}{\Gamma }_{M_{PS}}},$$

$$i{\Gamma }_{\mu \nu }^{(M_{PS}\gamma \gamma )}(k_1,k_2)=-i{e}^2{\epsilon }_{\mu \nu \kappa \lambda }{k}_1^{\kappa }{k}_2^{\lambda }{F}_{{\gamma }^{\ast }{\gamma }^{\ast }\to M_{PS}}(0,0)F^{(M_{PS}\gamma \gamma )}(k^2),$$

$$F_{{\gamma }^{\ast }{\gamma }^{\ast }\to M_{PS}}(0,0)=\frac{1}{4{\pi }^2{f}_{M_{PS}}}=\frac{1}{4\pi {\alpha }_{em}}\sqrt{\frac{64\pi \Gamma (M_{PS}\to \gamma \gamma )}{m_{M_{PS}}^3}}.$$

$$F^{(M\gamma \gamma )}(t/u)=\exp \left(\frac{t/u-m_M^2}{{\Lambda }_{exp}^2}\right),$$

$$F^{(M\gamma \gamma )}(s)=\exp \left(-\frac{(s-m_M^2{)}^2}{{\Lambda }_{exp}^4}\right),$$

$$i{\Gamma }_{\mu \nu }^{(M_S\gamma \gamma )}(k_1,k_2)=\frac{e^2}{2m_{M_S}}\left(k^2{g}_{\mu \nu }-2k_{1\nu }{k}_{2\mu }\right)F_{{\gamma }^{\ast }{\gamma }^{\ast }\to M_S}(0,0)F^{(M_S\gamma \gamma )}(k^2),$$

$$F_{{\gamma }^{\ast }{\gamma }^{\ast }\to M_S}(0,0)=\frac{1}{4\pi {\alpha }_{em}}\sqrt{\frac{64\pi \Gamma (M_S\to \gamma \gamma )}{m_{M_S}}}.$$

$$k_{1\mu }{\Gamma }^{\mu \nu }(k_1,k_2)=0,k_{2\nu }{\Gamma }^{\mu \nu }(k_1,k_2)=0.$$

$$i{\mathcal{M}}_{{\lambda }_1{\lambda }_2\to {\lambda }_3{\lambda }_4}^{(M_T)}$$

$$\begin{array}{r}i{\Gamma }_{\mu \nu \kappa \lambda }^{(f_2\gamma \gamma )}(k_1,k_2)=i\left[2a_{f_2\gamma \gamma }{\Gamma }_{\mu \nu \kappa \lambda }^{(0)}(k_1,k_2)-b_{f_2\gamma \gamma }{\Gamma }_{\mu \nu \kappa \lambda }^{(2)}(k_1,k_2)\right]F^{(f_2\gamma \gamma )}(k^2).\end{array}$$

$$i{\Delta }_{\kappa \lambda ,\alpha \beta }^{(f_2)}(k)$$

$$a_{f_2^{\prime }(1525)\gamma \gamma }=\pm \frac{e^2}{4\pi }0.12{\mathrm{GeV}}^{-3},b_{f_2^{\prime }(1525)\gamma \gamma }=\pm \frac{e^2}{4\pi }0.30{\mathrm{GeV}}^{-1}.$$

$$b_{f_2(1565)\gamma \gamma }=\pm \frac{e^2}{4\pi }0.93{\mathrm{GeV}}^{-1},b_{a_2(1700)\gamma \gamma }=\pm \frac{e^2}{4\pi }0.52{\mathrm{GeV}}^{-1}.$$

$${\sigma }_{\gamma \gamma \to R\to \gamma \gamma }(s)$$

Full text

I Introduction

In the Standard Model the light-by-light scattering is usually assumed to proceed through lepton, quark and $W$ gauge boson loops Bohm:1994sf ; Jikia:1993tc ; Bern:2001dg ; Bardin:2009gq . The QCD and QED corrections to the fermion-loop contributions were discussed in Bern:2001dg and their contribution was found to be very small. In KLS2017 our group considered also double fluctuations of photons into light vector mesons and their subsequent soft interactions treated in the Regge theory. The correction turned out to be important at higher center-of-mass energies ($\sqrt{s}>2$ GeV) and very small scattering angles. In KSS2017 we considered in addition two-gluon exchange contribution. This mechanism, in contrast to the VDM-Regge one survives to larger scattering angles.

In the present paper we shall consider several pseudoscalar ($J^{PC}=0^{-+}$), scalar ($J^{PC}=0^{++}$), and tensor ($J^{PC}=2^{++}$) meson exchange contributions for the $\gamma\gamma\to\gamma\gamma$ elastic scattering. The meson-$\gamma$-$\gamma$ vertex functions depend on the quantum numbers of objects involved and form factors that are in principle a function of four-momenta of both photons and meson. So far the $\gamma$-$\gamma$-meson processes were investigated at $e^{+}e^{-}$ colliders for one or both virtual photons and on-shell $\pi^{0}$ BaBar_TFF_2009 ; Belle_TFF , $\eta$, $\eta^{\prime}$ BaBar_TFF_2011 , $\eta_{c}$ BaBar_TFF_2010 mesons and also from recent measurements by the Belle collaboration Masuda:2015yoh for $f_{0}(980)$ and $f_{2}(1270)$ mesons. Recently, in Danilkin the authors formulated light-by-light scattering sum rules, that lead to relations between the $\gamma^{*}\gamma$ transition form factors for $C$-even scalar, pseudoscalar, axial-vector and tensor mesons when assuming sum rule saturation.

The situation for $\gamma\gamma\to\gamma\gamma$ elastic scattering is different as here photons are on mass shell and meson is off-shell. We do not know very precisely the dependence of the form factor for off-shell mesons on their virtuality. For the $s$-channel we have time-like meson while for $t$- and $u$-channels space-like one (see Fig. 1). We shall discuss here how the results depend on the form factors. In the following we shall parametrise it in either exponential or monopole type. We expect that for the $t$ and $u$ diagrams the corresponding form factors can be rather hard (larger cut-off parameter) compared to purely hadronic processes as they describe the coupling of space-like meson to semi point-like objects (photons).

II Meson exchange contributions

Below we will discuss the expressions relevant to pseudoscalar ($M_{PS}$), scalar ($M_{S}$), and tensor ($M_{T}$) as well as spin-4 $f_{4}(2050)$ meson exchange contributions for the reaction

[TABLE]

with the four-momenta $p_{i}$ $(i=1,...,4)$ and the photon helicities, $\lambda_{i}\in\{1,-1\}$, indicated in brackets. Schematic diagrams for the the process (1) are given in Fig. 1.

The differential cross section for the angular distribution for the reaction (1) is given by

[TABLE]

where the explicit factor $1/2$ takes into account identity of photons. The kinematical variables used in the present paper are

[TABLE]

The amplitude for the reaction (1) with the meson exchanges is written as

[TABLE]

In Table 1 we have collected possible potential resonances that may contribute to the process (1). The contribution of axial-vector mesons vanishes for on-shell photons due to the Landau-Yang theorem LY . The two-photon branching fractions for the resonances are relatively well known and were measured in recent years by the Belle and BaBar collaborations.

II.1 Pseudoscalar meson exchanges

The amplitude for the pseudoscalar meson exchange is written as

[TABLE]

where $\epsilon_{i}^{\mu_{i}}$ are the polarisation vectors of the photons with the helicities $\lambda_{i}$.

The propagator and $M_{PS}\gamma\gamma$ vertex function for the pseudoscalar meson with mass $m_{M_{PS}}$ are:

[TABLE]

where $k^{2}=p_{s}^{2}$, $p_{t}^{2}$, $p_{u}^{2}$ for the $s$, $t$, $u$ diagrams, respectively. The transition form factor $F_{\gamma^{*}\gamma^{*}\to M_{PS}}(0,0)$ is related to the two-photon decay width as

[TABLE]

Above $\alpha_{em}=e^{2}/(4\pi)$ and $f_{M_{PS}}$ is a decay constant of psudoscalar meson. We have, from Eq. (8), $f_{M_{PS}}$ = 0.092, 0.093, 0.074, 0.375, 0.776 GeV for $\pi^{0}$, $\eta$, $\eta^{\prime}$, $\eta_{c}(1S)$, $\eta_{c}(2S)$ mesons, respectively. The $f_{M_{PS}}$ can be tried to be calculated within quark model. For $\eta$ and $\eta^{\prime}$ this requires inclusion of $\eta$-$\eta^{\prime}$ mixing (see e.g. DGH1992 ). The two-photon decay of heavy quarkonium $\eta_{c}$ states was discussed, e.g., in LP_etac .

In our vertices (7) the photons are on mass shell and the meson is off-shell. Therefore in calculations we use off-shell meson form factors normalized to $F^{(M\gamma\gamma)}(m_{M}^{2})=1$. These form factors can be parametrised, e.g., by the exponential form for the $t$- and $u$-channel meson exchanges

[TABLE]

while for the $s$-channel meson exchange we assume

[TABLE]

where $\Lambda_{exp}$ is, in principle, a free parameter. In our calculation for pseudoscalar mesons we take $\Lambda_{exp}=2$ GeV.

II.2 Scalar meson exchanges

In our calculations here we consider four light scalar resonances: $f_{0}(500)$, $f_{0}(980)$, $a_{0}(980)$ and $f_{0}(1370)$. The first two states have been seen in $\gamma\gamma\to\pi\pi$ reactions Belle_pipi , while the $a_{0}(980)$ resonance was observed in the $\gamma\gamma\to\pi^{0}\eta$ reaction Belle_pieta . The situation for light scalar mesons was discussed, e.g., in AS_first . While the direct coupling is expected to be rather small, the rescattering of mesons ($\pi\pi$, $KK$ or $\pi\eta$) leads to two-photon widths of the order of a fraction of keV or even larger. The partial decay widths were estimated, e.g., in AABN1999 ; Dai_Pennington:2014 and are collected in Table 1. Note that the masses and widths of low-mass resonances are known with a large uncertainties PDG ; see also analysis in ALEPH for the $f_{0}(1500)$ and $f_{0}(1710)$ states. The two-photon decay of heavy quarkonium states ($\chi_{c0,2}$ and $\chi_{b0,2}$) was discussed, e.g., in LP_chic ; Danilkin_chic .

For the scalar meson $M_{S}$ exchange the amplitude is similar as (LABEL:amplitude_2to2_PS) with $\Delta^{(M_{PS})}$ and $\Gamma^{(M_{PS}\gamma\gamma)}_{\mu\nu}$ replaced by $\Delta^{(M_{S})}$ and $\Gamma^{(M_{S}\gamma\gamma)}_{\mu\nu}$, respectively. The corresponding vertex can be written as:

[TABLE]

where

[TABLE]

It is easy to check that

[TABLE]

Thus our vertices and as a consequence our $s$-, $t$- and $u$-channel amplitudes are gauge invariant.

For light scalar $f_{0}(500)$, $f_{0}(980)$, $a_{0}(980)$ and $f_{0}(1370)$ meson exchanges we take a smaller value of $\Lambda_{exp}=1$ GeV. 222This is related to the fact that the $\gamma\gamma$ coupling to the scalar mesons is dominantly not to quark-antiquark system but rather to mesonic loops. Their contribution will be shown in final summary plots.

II.3 Tensor meson exchanges

The tensor meson exchanges are much more complicated. The $f_{2}(1270)$ and $a_{2}(1320)$ couplings to two photons were studied in detail, e.g., in EMN14 (see section 5). Using the formalism from EMN14 the amplitude for the tensor meson $M_{T}$ exchange can be written as

[TABLE]

The $f_{2}(1270)\gamma\gamma$ vertex is given as (see formulae (3.39), (3.40) and section 5.3 of EMN14 )

[TABLE]

Here the so-called helicity-0 and helicity-2 $f_{2}\to\gamma\gamma$ amplitudes are parametrised by the $a$ and $b$ constants, respectively. Two rank-four tensor functions, $\Gamma_{\mu\nu\kappa\lambda}^{(0)}$ and $\Gamma_{\mu\nu\kappa\lambda}^{(2)}$, are defined by (3.18) and (3.19) in EMN14 . In our calculation for the tensor mesons ($q\bar{q}$ states) we take $\Lambda_{exp}=2$ GeV in form factors given by (9) and (10). We adopt also the numerical values of the $a_{2}(1320)\gamma\gamma$ coupling constants discussed in sections 5.4 and 7.2 of EMN14 .

For a better analysis we shall use model for the $f_{2}$ propagator considered in EMN14 ; see Appendix A there. The Breit-Wigner formula is used here

[TABLE]

where $\hat{g}_{\mu\nu}=-g_{\mu\nu}+k_{\mu}k_{\nu}/k^{2}$.

In our work we consider also the $f_{2}^{\prime}(1525)$ meson contribution. We assume the same ratio of the two-photon decay widths of the specific helicity to the total two-photon decay width as for the $f_{2}(1270)$ meson. We have

[TABLE]

A contribution of heavier tensor-meson states are also possible, see Table 1. Recently in Danilkin a role of $f_{2}(1565)$ and $a_{2}(1700)$ mesons was discussed in the context of transition form factors. Here we also take into account these two states. In order to estimate the strength of couplings we assume only helicity-2 contributions. From Eq. (5.28) of EMN14 and with the parameters of Table 1 we get

[TABLE]

In calculation we choose the positive values of above coupling constants.

The $t/u$-channel exchanges of spin-2 particles leads to a growing of the cross section as a function of $\sqrt{s}$. In practice the rapid growth appears quickly above the $s$-channel resonance. No clear unitarization procedure is known to as. From phenomenology we know that the exchange of the meson must be replaced by the exchange of the $f_{2}$ reggeon EMN14 . In the following we shall therefore omit the $t$- and $u$-channel meson exchange contributions.

II.4 $f_{4}(2050)$ meson exchange

The $f_{4}(2050)$ resonance was identified in $\gamma\gamma\to\pi\pi$ processes KS2013 . However, the detailed tensorial coupling for $f_{4}(2050)\to\gamma\gamma$ is not known, so as a consequence we do not know corresponding angular distributions. We shall therefore neglect $t$- and $u$-channel contributions. In the following we use a simple Breit-Wigner resonance form for a spin-$J$ resonance $R$: 333In the general case the initial particles have spins $s_{1}$ and $s_{2}$. In our case, for a real photons, a factors $(2s_{i}+1)$ are replaced by 2.

[TABLE]

which represents only $s$-channel contribution to the total cross section. The resonance form reproduces result of our diagrammatic calculation for lower-spin resonances at $\sqrt{s}\simeq m_{R}$. In (LABEL:sigma_f4) the resonance width $\Gamma_{R}$ is a constant. An improved description of resonance shapes can be obtained when the width is made $s$-dependent.

III Results

In this section we present first numerical results for meson exchange contributions to the $\gamma\gamma\to\gamma\gamma$ scattering. We shall include not only $s$-channel resonant contributions but also $t$- and $u$-channel exchanges, not discussed so far in the literature.

We start with contributions of the pseudoscalar mesons that were observed in two-photon processes in the $e^{+}e^{-}$ collisions (Belle, BaBar). In the left panel of Fig. 2 we show, for example, separate contributions of $s$, $t$ and $u$ channels and their coherent sum for $\eta^{\prime}(958)$ meson. The $t$ and $u$ channels play a role only for $\sqrt{s}>2$ GeV. In the right panel we show separate contributions of $\pi^{0}$, $\eta$ and $\eta^{\prime}$ exchanges as well as their sum. A large interference effects can be seen.

In Fig. 3 we present an example of angular distribution at $\sqrt{s}=1.3$ GeV for the light pseudoscalar mesons (similar distribution for the same energy will be discussed for tensor meson contributions). The resulting distribution is relatively flat which is caused by the dominance of the $s$-channel exchange at the selected energy (see the previous figure).

In the left panel of Fig. 4 we show separate contributions of tensor meson exchanges ($s$-channel only) as well as their coherent sum. Similarly as for the pseudoscalar meson exchange contributions relatively large interference effects can be seen. In the right panel we show angular distributions for $\sqrt{s}=1.3$ GeV. The resulting tensor meson distribution shows an enhancement at $\cos\theta\approx\pm 1$ (compare Fig. 3 for pseudoscalar meson contributions). The subdominant (helicity-0) $f_{2}(1270)\gamma\gamma$ coupling ($\Gamma^{(0)}$ in (15)) plays important role for $\cos\theta\approx 0$.

Finally, in Fig. 5 we show energy dependence of all meson exchange contributions integrated over full $z=\cos\theta$ range (left panel) and for $|z|<0.6$ (right panel). For reference we show also the fermion-box (lepton and quark) contributions separately and added coherently together. The $W$-boson loops contribution plays a role only for $\sqrt{s}>200$ GeV LPS14 . The mesonic contribution exceed the quark-box one for different intervals of $\sqrt{s}$. Only a few narrow resonances clearly exceed the total fermion continuum. The reader is asked to note a sizeable interference effect between pseudoscalar and tensor meson contributions.

The $\pi^{0}$, $\eta$, $\eta^{\prime}(958)$ states were already observed in the $e^{+}e^{-}$ collisions. It seems rather difficult to observe the heavy quarkonia $\eta_{c}(1S)$, $\eta_{c}(2S)$ and $\chi_{c0}$. In heavy ion ultraperipheral collisions at CERN this is because of poor two-photon invariant mass resolution of the order of 0.5 GeV while in low-energy $e^{+}e^{-}$ collisions because of limited phase space and the presence of two-photon bremsstrahlung background. The region of $f_{2}(1270)$ seems quite interesting as here some enhancement could be potentially identified by the Belle II at SuperKEKB for instance. Imposing a cut $|z|<0.6$ (see the right panel of Fig. 5) improves the signal (meson exchanges) to background (boxes) ratio.

The meson exchange contributions are limited only to $\sqrt{s}<4$ GeV, and should not influence the recent ATLAS experimental result Aaboud:2017bwk .

IV Conclusions

In the present work we have discussed for the first time the role of $s$-, $t$-, $u$-channel meson exchanges in elastic $\gamma\gamma\to\gamma\gamma$ scattering. We have included a few mesonic states with large two-photon branching fractions, see Table 1. Large interference effects for light pseudoscalar and tensor mesons have been found. The interference of $\pi^{0}$, $\eta$ and $\eta^{\prime}$ amplitudes increases the cross section for $\sqrt{s}>m_{\eta^{\prime}}$ by almost an order of magnitude. The interference of $f_{2}(1270)$, $a_{2}(1320)$ and the other tensor mesons leads to interesting spectral shape for collision energies in the window between 1 – 2 GeV. Could this be measured at SuperKEKB in a future? Not excluded, provided two-photon bremsstrahlung eliminated experimentally or is taken into account in the calculations.

We have included not only the $s$-channel diagrams (leading to peaks at $\sqrt{s}\simeq m_{M}$) but also the meson exchanges in $t$- and $u$-channels (leading to broad continua). In general, the $t$ and $u$ diagrams contribute above the resonance peaks associated with the $s$-channel exchanges which is caused by kinematics of the process and gauge invariance explicitly fulfilled (imposed) in our calculation. This observation is true for the pseudoscalar, scalar and tensor meson exchanges. Only the $s$-channel contributions play a role for heavy pseudoscalar ($\eta_{c}(1S)$, $\eta_{c}(2S)$) and scalar ($\chi_{c0}$) mesons and the $t/u$-channel contributions may be safely neglected.

The results related to meson exchanges have been compared with the standard, in the context of $\gamma\gamma\to\gamma\gamma$ scattering, fermion-box continuum. The mesonic contributions concentrate in the region of $\sqrt{s}\in(0.1,4.0)$ GeV and are in general smaller than box contributions, except of some specific regions of the phase space. For instance, at the resonance positions the meson exchange contributions sometimes even strongly exceed the standard box ones but experimental observation may depend on diphoton invariant mass ($M_{\gamma\gamma}$) resolution of a given experiment.

The exchanges of light pseudoscalar ($\pi^{0}$, $\eta$ and $\eta^{\prime}$) and tensor mesons is also important in the context of possible interference with the fermion-box contributions. This goes, however, beyond the scope of the present paper and will be studied elsewhere.

The meson exchange contributions discussed in the present paper could potentially influence the ${}^{208}Pb+^{208}Pb\to^{208}Pb+^{208}Pb+\gamma\gamma$ ultraperipheral collisions measured recently by the ATLAS collaboration Aaboud:2017bwk . Our present analysis suggest, however, that the meson exchange contributions do not play important role for the recent ATLAS measurement where $M_{\gamma\gamma}>6$ GeV.

Acknowledgements.

We are indebted to Sadaharu Uehara for a discussion on the Belle measurements of mesons in the $\gamma\gamma$ channel. This research was partially supported by the Polish Ministry of Science and Higher Education Grant No. IP2014 025173 (Iuventus Plus), the Polish National Science Centre Grant No. DEC-2014/15/B/ST2/02528 (OPUS) and by the Center for Innovation and Transfer of Natural Sciences and Engineering Knowledge in Rzeszów.