EXCLUSIVE PHOTOPRODUCTION OF A DI-MESON PAIR WITH LARGE INVARIANT MASS

Abstract

We consider the exclusive photoproduction of a di-meson pair with large invariant mass, $\gamma N\rightarrow N^{\prime}M_{1}M_{2}$, in the framework of collinear factorisation. The mesons considered $M_{1}$ and $M_{2}$ are either pions or rho mesons, charged or neutral. We consider the kinematic regime characterised by a large invariant mass of the two-meson system, and a small deflection of the nucleon in the centre-of-mass frame. In this kinematic domain, the amplitude factorises into a perturbative hard part and non-perturbative parts described by Generalised Parton Distributions (GPDs) and Distribution Amplitudes (DAs). We automate the calculation of the fully differential cross section at leading twist and leading order, and we present some numerical results at JLab 12 GeV kinematics. This class of processes provides yet more exclusive $2\to 3$ channels that can be used to extract GPDs.

1 Introduction

Understanding the internal dynamics of nucleons has been one of the most active research fields for the last few decades. Beyond the Parton Distribution Functions (PDFs) which describe the probability of finding a parton carrying longitudinal momentum fraction $x$ inside the nucleon, it has now become possible to access the transverse spatial distribution through Generalised Parton Distributions (GPDs) [1, 2]. These functions appear as the non-perturbative part of the amplitude for exclusive processes where collinear factorisation is expected to hold.

Here, we focus our attention on the family of $2\to 3$ exclusive processes of the form $\gamma N\rightarrow N^{\prime}M_{1}M_{2}$, namely the exclusive photoproduction of a di-meson pair [3]. We work in the kinematics where the di-meson pair has a large invariant mass (provided by the large relative transverse momentum of the final-state mesons), and small squared momentum transfer to the nucleon. In such a case, it can be shown that the process admits a collinear factorisation approach, and thus provides access to GPDs [4, 5]. Depending on the spin and polarisation of the two mesons, chiral-even (helicity non-flip) and/or chiral-odd (helicity-flip) GPDs contribute to the amplitude.

Our work aims to provide theoretical predictions for the considered class of exclusive processes, which will be indispensable for future experimental studies, by assessing their sensitivity to GPD modeling, especially with respect to chiral-odd GPD which are poorly constrained experimentally. In this report, we restrict ourselves to Leading-Order (LO) in $\alpha_{s}$ at the leading twist, and to processes that do not involve 2-gluon exchanges in the $t$-channel, in order to avoid the factorisation breaking effects described in [6, 7].

2 Kinematics

We label the momenta and polarisation vectors in the process as follows:

$$\gamma(q,\epsilon_{q})+N(p_{1})\rightarrow N^{\prime}(p_{2})+M_{1}(p_{M_{1}},\epsilon_{M_{1}})+M_{2}(p_{M_{2}},\epsilon_{M_{2}})\,,$$ (1)

where $\epsilon_{1}$ and $\epsilon_{2}$ are the polarisations of the two outgoing mesons if they are vector mesons.

We define the momenta

$$q^{\mu}=\frac{\sqrt{s}}{2}(1,0,0,-1)\,,\qquad p^{\mu}=\frac{\sqrt{s}}{2}(1,0,0,1)\,.$$ (2)

Any vector $v^{\mu}$ can be decomposed in the Sudakov basis as

$$v^{\mu}=v^{+}p^{\mu}+v^{-}q^{\mu}+v_{\perp}^{\mu}\,.$$ (3)

We further define $\Delta=p_{2}-p_{1}$, and introduce the usual skewness parameter $\xi=-\frac{\Delta^{+}}{2p^{+}}$.

In the limit where the deflection of the nucleon is negligible ($\Delta_{\perp}$ small compared to the hard scale in the process), and if we neglect all hadron masses, the momenta of the particles in the process are given by

$$p_{1}^{\mu}=(1+\xi)\,p^{\mu},\quad\;p_{2}^{\mu}=(1-\xi)\,p^{\mu}\,,\quad\;p_{M_{1}}^{\mu}=(1-\alpha)\,q^{\mu}+2\alpha\xi\,p^{\mu}-p^{\mu}_{\perp}\,,\quad\;p_{M_{2}}^{\mu}=\alpha\,q^{\mu}+2(1-\alpha)\xi\,p^{\mu}+p^{\mu}_{\perp}\,.$$ (4)

One can show (see [8] for details) that in this limit, we have:

$$S_{\gamma N}\equiv(q+p_{1})^{2}=(1+\xi)s\,,$$ (5)

$$M_{12}^{2}\equiv(p_{M_{1}}+p_{M_{2}})^{2}=2\xi s\,.$$ (6)

From these equations, it follows that

$$-t^{\prime} \equiv-(q-p_{M_{2}})^{2}=(1-\alpha)\,M_{12}^{2}\,,$$ (7)

$$-u^{\prime} \equiv-(q-p_{M_{1}})^{2}=\alpha\,M_{12}^{2}\,,$$ (8)

$$s^{\prime} \equiv(p_{M_{1}}+p_{M_{2}})^{2}=M_{12}^{2}=\frac{\vec{p_{t}}^{2}}{\alpha(1-\alpha)}\,.$$ (9)

Since $M_{12}^{2}$ is positive, the last equality implies that $0<\alpha<1$, so that $-t^{\prime}$ and $-u^{\prime}$ are also positive.

The polarisation vectors of the outgoing mesons can be written as [8]

$$\epsilon^{\mu}_{M_{i}}(p_{M_{i}},L) =\frac{1}{m_{M_{i}}}p^{\mu}_{M_{i}}-\frac{m_{M_{i}}}{p\cdot p_{M_{i}}}p^{\mu}\,,$$ (10)

$$\epsilon^{\mu}_{M_{i}}(p_{M_{i}},T) =\epsilon_{M_{i}\perp}^{\mu}-\frac{\epsilon_{M_{i}\perp}\cdot p_{M_{i}}}{p\cdot p_{M_{i}}}p^{\mu}\,,$$ (11)

where $L$ and $T$ stand for longitudinal and transverse respectively, and $m_{M_{i}}$ represents the mass of the light meson $M_{i}$.

We choose the axial gauge, $\epsilon_{q}\cdot p=0$ for the photon polarisation vector, which leads to

$$\epsilon_{q}^{\mu}=\epsilon_{q\perp}^{\mu}\,.$$ (12)

In the chosen kinematics, the amplitudes for our class of processes are functions of the three parameters $\alpha$, $\xi$, $s$ and of the possible tensor structures (scalar products or Levi-Civita tensors) that can be generated with $p_{\perp}$, $\epsilon_{1,2\perp}$, $\epsilon_{q\perp}$, $p$ and $q$.

3 Non-perturbative ingredients

In order to make predictions, we need to model the non-perturbative parts of the amplitude. Within the collinear factorisation framework, the amplitude factorises into a perturbative hard part, denoted by $T_{H}(x,v,z)$, and a non-perturbative part, which consists of a Generalised Parton Distribution (GPD), $H(x)$, and two Distribution Amplitudes (DAs), $\phi_{1}(v)$ and $\phi_{2}(z)$. The GPD encodes the internal dynamics of the nucleon, while the DAs describe those of each meson. In the GPD, $x$ represents the average longitudinal momentum of the parton probed from the nucleon, while in the DAs, $v$ and $z$ represent the momentum fraction of each meson carried by the quark.

The factorised amplitude is then given by

$$i\mathcal{M}=\int_{-1}^{1}dx\int_{0}^{1}dz\int_{0}^{1}dv\;T_{H}(x,v,z)H(x)\phi_{1}(v)\phi_{2}(z)\,.$$ (13)

There are two kinds of quark GPDs, namely chiral-even and chiral-odd GPDs. The first ones describe processes in which the total transfer of helicity between the initial and final nucleon is zero. They can be further distinguished into either vector, denoted by $H^{q},E^{q}$, or axial $\tilde{H}^{q},\tilde{E^{q}}$, with $q$ being the flavour of the active parton. In our study, only $H^{q}$ and $\tilde{H}^{q}$ are kept, the other being suppressed due to multiplication by powers of the skewness $\xi$ in the cross section [9]. The chiral odd GPDs, denoted by $H_{T}^{q},\,\tilde{H}_{T}^{q},\,E_{T}^{q}$ and $\tilde{E}_{T}^{q}$ describe a transfer of one unit of helicity in the $t$-channel. The dominant contribution comes from the transversity GPD $H_{T}^{q}$, the other ones being suppressed by powers of $\xi$.

The GPD and the DAs appear in the amplitude as matrix elements of bilocal operators at light-like separation. For example, for the vector GPD, we have

$$\langle N^{\prime}(p_{2})|\bar{q}\left(-\frac{z}{2}\right)\gamma^{+}q\left(\frac{z}{2}\right)|N(p_{1})\rangle\Big|_{\begin{subarray}{c}z_{\perp}=0\\ z^{+}=0\end{subarray}}=\int_{-1}^{1}dx\,e^{-ixP^{+}z^{-}}\bar{u}(p_{2})\gamma^{+}u(p_{1})\,H^{q}(x,\xi,t)\,,$$ (14)

while, for an axial DA, as the one of a $\pi^{+}$,

$$\langle\pi^{+}(p_{\pi})|\bar{u}(y)\gamma^{5}\gamma^{\mu}d(-y)|0\rangle=if_{\pi}p_{\pi}^{\mu}\int_{0}^{1}dz\,e^{-i(2z-1)p_{\pi}\cdot y}\phi_{\pi}(z)\,,$$ (15)

where $f_{\pi}$ is the pion decay constant.

Various models for GPDs and DAs exist. For concreteness, we adopt Radyushkin’s Double Distribution Ansatz for the GPDs [10],

$${H}^{q}(x,\xi,t)=\int_{\{|\beta|+|\alpha|\leq 1\}}d\beta\,d\alpha\,\delta(\beta+\xi\alpha-x)F^{q}(\beta,\alpha,t)\,,$$ (16)

where $F^{q}$ is constructed using a profile function $\Pi(\beta,\alpha)=\frac{3}{4}\frac{(1-\beta)^{2}-\alpha^{2}}{(1-\beta)^{3}}$ and PDFs through

$$F^{q}(\beta,\alpha,t)=\left(\Pi(\beta,\alpha){q}(\beta)\Theta(\beta)-\Pi(-\beta,\alpha){\bar{q}}(-\beta)\Theta(-\beta)\right)\frac{C^{2}}{(t-C)^{2}}\,,$$ (17)

$$\tilde{F}^{q}(\beta,\alpha,t)=\left(\Pi(\beta,\alpha){\Delta q}(\beta)\Theta(\beta)-\Pi(-\beta,\alpha){\Delta\bar{q}}(-\beta)\Theta(-\beta)\right)\frac{C^{2}}{(t-C)^{2}}\,,$$ (18)

$$F^{q}_{T}(\beta,\alpha,t)=\left(\Pi(\beta,\alpha){\delta q}(\beta)\Theta(\beta)-\Pi(-\beta,\alpha){\delta\bar{q}}(-\beta)\Theta(-\beta)\right)\frac{C^{2}}{(t-C)^{2}}\,.$$ (19)

In the above, $q$, $\Delta q$ and $\delta q$ denote respectively the unpolarised, polarised and transversity PDF of flavour $q$. The constant $C$ is fixed to $C=$0.71\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$. The parametrisation adopted here is the same as the one used in the study of the photoproduction of a photon-meson pair reported in [11, 12, 9, 8]. The manipulation of PDFs is facilitated by the use of the Mathematica package ManeParse [13].

For the DAs, we use the asymptotic form, given by

$$\phi_{i}(v)=6\,v\,(1-v)\,,\qquad i=1,2\,.$$ (20)

4 Automation of amplitude computation

4.1 Generation and storage of the amplitudes

In order to generate all possible diagrams without fixing the particle species, we generate the diagrams of the process $\gamma(q)\rightarrow u(k_{1})\,\bar{u}(k_{2})\,u(k_{3})\,\bar{u}(k_{4})\,u(k_{5})\,\bar{u}(k_{6})$ using the FeynArts and FeynCalc packages [14, 15, 16, 17, 18] in Mathematica. The momenta of the outgoing quarks and anti-quarks are represented by $k_{i}$, $i=1,...,6$. These diagrams, such as the one shown on the left-hand side of Fig. 1, are then mapped onto the desired form by extracting the chains of Dirac matrices between spinors, and assembling them into a Dirac trace. For example, for the specific diagram in Fig. 1,

Figure 1: The projection of a diagram obtained with FeynArts and FeynCalc onto the GPD and DAs.

$$\frac{{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}\bar{u}(k_{1})}\boldsymbol{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}{\gamma}_{\mu}}{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}{v(k_{4})\bar{u}(k_{5})}}\boldsymbol{\color[rgb]{1,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{1,0,0}\gamma_{\nu}}{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}v(k_{6})\bar{u}(k_{3})}{\boldsymbol{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}{\gamma}^{\mu}(\not{k}_{1}+\not{k}_{3}+\not{k}_{4})\not{\varepsilon}(q1)(-\not{k}_{2}-\not{k}_{5}-\not{k}_{6})\gamma^{\nu}}}{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}v(k_{2})}}{(k_{1}+k_{4})^{2}(k_{5}+k_{6})^{2}(k_{2}+k_{5}+k_{6})^{2}(k_{1}+k_{3}+k_{4})^{2}}$$ (21)

                        $\Downarrow$

$$\frac{\mathrm{Tr}\left(\boldsymbol{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\gamma^{\mu}(\not{k}_{1}+\not{k}_{3}+\not{k}_{4})\not{\varepsilon}(q_{1})(-\not{k}_{2}-\not{k}_{5}-\not{k}_{6})\gamma^{\nu}}{\color[rgb]{.75,.5,.25}\definecolor[named]{pgfstrokecolor}{rgb}{.75,.5,.25}\boxed{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}GPD}}\right)\mathrm{Tr}\left(\boldsymbol{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\gamma_{\mu}}{\color[rgb]{.75,.5,.25}\definecolor[named]{pgfstrokecolor}{rgb}{.75,.5,.25}\boxed{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}DA_{1}}}\boldsymbol{\color[rgb]{1,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{1,0,0}\gamma_{\nu}}{\color[rgb]{.75,.5,.25}\definecolor[named]{pgfstrokecolor}{rgb}{.75,.5,.25}\boxed{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0}DA_{2}}}\right)}{(k_{1}+k_{4})^{2}(k_{5}+k_{6})^{2}(k_{2}+k_{5}+k_{6})^{2}(k_{1}+k_{3}+k_{4})^{2}}$$ (22)

Each GPD and DA enters the scattering amplitude multiplied by a Dirac matrix, which is absorbed into the hard part $T_{H}$ in Eq. (13). These Dirac matrices are symbolically written as brown boxes in Eq. (22). Depending on the nature of the GPD or DA, those matrices can be $\gamma^{\mu}$ (vector), $\gamma^{5}\gamma^{\mu}$ (axial), or $\frac{i}{2}[\gamma^{\mu},\gamma^{\nu}]$ (tensor, which corresponds to transversity). There are $3^{3}=27$ possible combinations of Dirac structures. However, the tensor Dirac structures must appear in pairs, otherwise the amplitude vanishes due to the presence of an odd number of gamma matrices in the trace. This reduces the number of combinations to 14. Therefore, each diagram on the right-hand side of Fig. 1 can be projected onto 14 different possible combinations of Dirac structures. This gives a total of 1480 diagrams which are organised according to their topology. The five possible topologies are represented in Fig. 2.

A flag containing all the structural information is associated with each diagram. As an illustration, a generic amplitude, corresponding to the right panel of Fig. 1, is given by

$$\frac{eg^{4}C_{F}\left(\epsilon_{q\perp}\cdot p_{\perp}\right)\left(\alpha\xi\left(-4v^{2}+4v-4z^{2}+4z-2\right)+\xi(1-2v)^{2}+2(\alpha-1)vx-2\alpha xz+x\right)}{64\xi^{2}{N_{c}}^{2}s^{2}(v-1)v(z-1)z(\alpha(v+z-1)-v)(\alpha(v+z-1)-v+1)\left(i\epsilon+\frac{\xi(-\alpha(v+z-1)-2vz+v+2z-1)}{\alpha(v+z-1)-v+1}+x\right)}$$

$$\hskip 170.71652pt\times\frac{f(\{244,1,1,\{3,2\}\},N,N^{\prime},\{M_{1},P_{1}\},\{M_{2},P_{2}\})}{\left(-i\epsilon+\frac{\xi(\alpha(v+z-1)-2vz+v)}{\alpha(v+z-1)-v}+x\right)}\,.$$ (23)

The flag $f(\{244,1,1,\{3,2\}\},N,N^{\prime},\{M_{1},P_{1}\},\{M_{2},P_{2}\})$ has a list of four arguments as its first entry, which respectively indicates that this amplitude

• originates from diagram number 244,

• was calculated using only vector GPD and DAs (combination number 1, out of 14),

• belongs to the first topology (out of 5) shown in Fig. 2, and

• the photon is attached to the fermion line connecting $k_{2}$ to $k_{3}$.

The variables $N$, $N^{\prime}$, $M_{1}$, $M_{2}$, $P_{1}$, $P_{2}$ refer respectively to the incoming and outgoing nucleon species, the first and second type of meson, and their polarisations $P_{1}$ and $P_{2}$. Once these are specified, the flag becomes the prefactor of the diagram for the given process. This prefactor includes the electric charge (of the quark line which is connected to the photon), the decay constants of the mesons, and the GPD and DAs.

Figure 2: Five ways of connecting the fermion lines, giving five topologies. The unmarked blobs correspond to places where the internal gluons and/or incoming photon could attach to.

For example, for two longitudinal $\rho^{0}$ mesons produced from a proton, this prefactor reads

$$f(\{244,1,1,\{3,2\}\},p,p,\{\rho^{0},L\},\{\rho^{0},L\})=f_{\rho}^{2}\,\phi_{1}(v)\phi_{2}(z)({Q_{d}}{H_{d}}(x)+Q_{u}H_{u}(x))\,.$$ (24)

One recognises the DAs $\phi_{1}(v)$ and $\phi_{2}(z)$, the square of the rho meson decay constant $f_{\rho}$, the quark electric charges $Q_{u}$ and $Q_{d}$, and the vector GPDs $H_{u}(x)$ and $H_{d}(x)$ corresponding to the quark flavours $u$ and $d$. The factor $1/2$ that is supposed to appear due to the square root of 2 in the $\rho^{0}$ wave function is in fact compensated by the equal contribution from the $u$ and $d$ parts of this wave function. Indeed, this diagram belongs to first topology in Fig. 2, where the photon is attached to the lower blob. Thus, the charge seen by the photon, as well as the flavour of the GPD, are completely independent of the flavours taken by the fermionic lines connecting the DAs. Moreover, this specific topology constrains the flavours of both mesons to be the same, both being either $u$ or $d$. In this latter case, the minus sign coming from the two $\rho^{0}$ wave functions compensate, and therefore, each flavour of the $\rho^{0}$ meson has the same contribution to this amplitude, including the sign. It is this factor of 2 that compensates the above-mentioned factor of 1/2.

Therefore, if one chooses a specific process, all prefactors are replaced in the list of diagrams. Often, most of them actually end up being zero, since a specific process could be completely incompatible with one or more topologies in Fig. 2. For example, the first topology in Fig. 2 is incompatible with a transversity GPD (tensor Dirac structure) due to the presence of an odd number of gamma matrices in the trace.

The next step is to integrate the non-vanishing diagrams over the variables $x$, $v$, and $z$, which have support $x\in[-1,1]$, $z,v\in[0,1]$.

4.2 Numerical integration

Unlike the previous specific study of exclusive $\pi^{+}\rho^{0}_{T}$ photoproduction in [19], we perform the integration over all variables numerically. One important advantage of doing this is that one is able to investigate different models for the DAs. First, a partial fraction decomposition with respect to $x$ is performed. Second, each resulting term is further split using the Sokhotski-Plemelj formula, which results in

$$\int_{-1}^{1}dx\frac{g(x,v,z)}{x-r(v,z)\pm i\epsilon}=\int_{-1}^{1}dx\frac{g(x,v,z)-g(r,v,z)}{x-r(v,z)}+\left[\ln\left(\frac{1-r(v,z)}{1+r(v,z)}\right)\mp i\pi\right]g(r,v,z)\,,$$ (25)

where $g(x,v,z)$ denotes the numerator of the term, and $r(v,z)$ denotes the real part of the pole in $x$. Note that the only $x$ dependence in the function $g$ corresponds to that of the relevant GPD. The first term on the right-hand side has to be integrated over $x$, $v$, and $z$, whereas the second term on the right-hand side has to be integrated over only $v$ and $z$.

One has to consider that, during the partial fraction decomposition, the $i\epsilon$ terms are distributed across both the numerator and denominator of the function $g(x,v,z)$. The $i\epsilon$ dependence in the numerator can be safely eliminated. However, some terms become non-integrable in $v$ or $z$ when $i\epsilon$ is eliminated in the denominator. A careful inspection shows that such terms are always of the form

$$\frac{z-v}{(z-v)^{2}+i\epsilon}\;\;\text{or}\;\;\frac{z+v-1}{(z+v-1)^{2}+i\epsilon}\,,$$ (26)

that is, the problematic part lies on the diagonals of the $(v,z)$ plane. We stress that the divergences that have appeared in $v$ and $z$ as a result of the partial fraction are spurious - Indeed, they were absent before performing the partial fractioning procedure. One convenient way to deal with them is to perform a series of foldings of the $(v,z)$ integration region, along the lines $z=v$, $z=1-v$, and finally $z=1/2$. This essentially corresponds to

$$\int_{0}^{1}dz\int_{0}^{1}dv\,\mathcal{A}(z,v)=\int_{0}^{\frac{1}{2}}\,dz\,\int_{0}^{z}\,dv \bigg(\mathcal{A}(z,v)+\mathcal{A}(v,z)+\mathcal{A}(1-v,1-z)+\mathcal{A}(1-z,1-v)$$

$$+\mathcal{A}(1-v,z)+\mathcal{A}(z,1-v)+\mathcal{A}(1-z,v)+\mathcal{A}(v,1-z)\bigg)\,,$$ (27)

where $\mathcal{A}$ denotes the term that contains spurious divergences in $v=z$ and/or $v=1-z$. In this way, all the $i\epsilon$ factors can be safely eliminated. For example, if we have a generic term of the form $\frac{h(v,z)(z-v)}{(z-v)^{2}+i\epsilon}$, then after performing the series of folding, one gets

$$\frac{h(v,z)(z-v)}{(z-v)^{2}+i\epsilon}+\frac{h(z,v)(v-z)}{(v-z)^{2}+i\epsilon}+(z\rightarrow 1-z\;\,\text{and}\;\,v\rightarrow 1-v)\underset{v\rightarrow z}{\sim}\frac{(v-z)^{2}\;\left(\frac{\partial{h}(z,v)}{\partial{v}}-\frac{\partial{h(z,v)}}{\partial{z}}\right)}{(v-z)^{2}+i\epsilon}+...\,,$$ (28)

where $i\epsilon$ can be safely taken to 0 before integration on the right-hand side.

5 Numerical results: Unpolarised cross section

Figure 3: Fully differential cross section as a function of $-u^{\prime}$ for the photoproduction of $\pi^{-}\rho^{0}_{L}$ (top-left), $\pi^{-}\rho^{0}_{T}$ (top-right), $\pi^{+}\rho^{0}_{L}$ (bottom-left) and $\pi^{+}\rho^{0}_{T}$ (bottom-right) at $S_{\gamma N}=$20\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$ and $M^{2}_{12}=$3\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$.

Figure 4: Fully differential cross sections as a function of $-u^{\prime}$ for the photoproduction of $\pi^{0}\rho^{+}_{L}$ (top-left), $\pi^{0}\rho^{+}_{T}$ (top-right), $\pi^{0}\rho^{-}_{L}$ (bottom-left) and $\pi^{0}\rho^{-}_{T}$ (bottom-right) at $S_{\gamma N}=$20\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$ and $M^{2}_{12}=$3\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$.

Let $|\overline{\mathcal{M}}|^{2}$ denote the amplitude squared averaged over the polarisation of the incoming photon and summed over the polarisations (if any) of the outgoing mesons. The fully differential cross section, expressed in terms of the kinematical variables $-u^{\prime}$, $-t$ and $M^{2}_{12}$, is

$$\frac{d\sigma}{d(-t)d(-u^{\prime})dM_{12}^{2}}=\frac{|\overline{\mathcal{M}}|^{2}}{32(2\pi)^{3}S_{\gamma N}M_{12}^{2}}\,.$$ (29)

In this proceeding, the values of $S_{\gamma N}$ and $M_{12}^{2}$ are taken to match the set-up of JLab in Hall B, in which the incoming electron beam has an energy of ${12}\,\mathrm{GeV}$. Constraints on the phase space (see e.g. [11]) require $M_{12}^{2}$ to lie between 2 and $9\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$, so we take $M_{12}^{2}=$3\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$ in the following to illustrate the results.

Some plots of the differential cross section as functions of the Mandelstam variable $u^{\prime}$ are shown in Fig. 3 and Fig. 4. The range of $(-u^{\prime})$ is fixed by the two constraints $-u^{\prime}>1\;\hbox{GeV}^{2}$ and $-t^{\prime}>1\;\hbox{GeV}^{2}$ to ensure that the kinematics corresponds to the regime where collinear factorisation is applicable [9]. While the cross sections of the exclusive photoproduction of a photon-meson pair [11, 12, 9, 8] did not exceed a few tens of $\hbox{pb}.\hbox{GeV}^{-6}$, here, the values can reach the order of a thousand, for processes with longitudinal $\rho^{0}_{L}$ in the final state.

Meson exchange symmetry acts as a good check of our results. The cross section obtained for the process $\gamma N\rightarrow N^{\prime}M_{1}M_{2}$ is related to that of the process $\gamma N\rightarrow N^{\prime}M_{2}M_{1}$ by the simple transformation of $-u^{\prime}\rightarrow M_{12}^{2}-(-u^{\prime})$, see Eqs. (7) and (8). More precisely, this transformation maps the differential cross section onto its reflection with respect to the line $u^{\prime}=M_{12}^{2}/2$. This has been checked for all of our processes. In the special case where $M_{1}$ and $M_{2}$ are identical mesons, with the same polarisation, the cross section becomes symmetric with respect to the $-u^{\prime}=M_{12}^{2}/2$ axis. This is indeed the case for the processes $\pi^{0}\pi^{0}$, $\rho^{0}_{L}\rho^{0}_{L}$ and $\rho^{0}_{T}\rho^{0}_{T}$ photoproduction as shown in Fig. 5 for a proton target. Finally, the same symmetry implies a perfect cancellation of the cross section at $u^{\prime}=M_{12}^{2}/2$ for the processes $\pi^{0}\pi^{0}$ and $\rho^{0}_{L}\rho^{0}_{L}$, which is clearly visible from the plots.

Figure 5: Fully differential cross sections for the photoproduction of $\rho^{0}_{L}\rho^{0}_{L}$ (top-left), $\rho^{0}_{T}\rho^{0}_{T}$ (top-right) and $\pi^{0}\pi^{0}$ (bottom) on a proton target at $S_{\gamma N}=$20\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$ and $M^{2}_{12}=$3\text{\,}\mathrm{G}\mathrm{e}\mathrm{V}^{2}$$. For the two first processes, the range of $-u^{\prime}$ has been extended beyond the kinematically allowed range (given the cuts $-u^{\prime}>1\;\hbox{GeV}^{2}$ and $-t^{\prime}>1\;\hbox{GeV}^{2}$, see main text) in order to observe the reflection symmetry of the cross section about the line $-u^{\prime}=M_{12}^{2}/2$.

6 Conclusion

By fully automating the differential cross section computation, we can make predictions for the whole class of $\gamma\,N\rightarrow N^{\prime}M_{1}M_{2}$ processes, where $M_{1}$ and $M_{2}$ can be rho mesons or pions of any charge and polarisation, excluding channels which allow the exchange of 2 gluons in the $t$-channel. This paves the way to a broad range of processes for the extraction of both chiral-even and chiral-odd GPDs. The cross sections that have been computed in JLab kinematics are one or two orders of magnitude above the ones for the $\gamma\,N\rightarrow N^{\prime}\gamma M$ processes, which is quite promising. We intend to also compute the linear polarisation asymmetry with respect to the incoming photon and integrated cross sections for JLab, LHC (in ultraperipheral collisions) and future EIC kinematics.

References