Monopole production via photon fusion and Drell–Yan processes: MadGraph implementation and perturbativity via velocity-dependent coupling and magnetic moment as novel features

In this work we consider point-like monopole production via photon-fusion and Drell–Yan processes in the framework of an effective U(1) gauge field theory obtained from conventional models describing the interaction of spin magnetically-charged fields with ordinary photons, upon electric-magnetic dualisation. We present arguments based on such dualities which support the conjecture of an effective monopole-velocity-dependent magnetic charge. For the cases of spin- and spin-1 monopoles, we also include a magnetic-moment term \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ, which is treated as a new phenomenological parameter and, together with the velocity-dependent coupling, allows for a perturbative treatment of the cross-section calculation. We discuss unitarity issues within these effective field theories, in particular we point out that in the spin-1 monopole case only the value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa =1$$\end{document}κ=1 may restore unitarity. However from an effective-field-theory point of view, this lack of unitarity should not be viewed as an impediment for the phenomenological studies and experimental searches of generic spin-1 monopoles, given that the potential appearance of new degrees of freedom in the ultraviolet completion of such models might restore it. The second part of the paper deals with an appropriate implementation of photon-fusion and Drell–Yan processes based on the above theoretical scenarios into MadGraph UFO models, aimed to serve as a useful tool in interpretations of monopole searches at colliders such as LHC, especially for photon fusion, given that it has not been considered by experimental collaborations so far. Moreover, the experimental implications of such perturbatively reliable monopole searches have been laid out.

Abstract In this work we consider point-like monopole production via photon-fusion and Drell-Yan processes in the framework of an effective U (1) gauge field theory obtained from conventional models describing the interaction of spin 0, 1 /2, 1 magnetically-charged fields with ordinary photons, upon electric-magnetic dualisation. We present arguments based on such dualities which support the conjecture of an effective monopole-velocity-dependent magnetic charge. For the cases of spin-1 /2 and spin-1 monopoles, we also include a magnetic-moment term κ, which is treated as a new phenomenological parameter and, together with the velocity-dependent coupling, allows for a perturbative treatment of the cross-section calculation. We discuss unitarity issues within these effective field theories, in particular we point out that in the spin-1 monopole case only the value κ = 1 may restore unitarity. However from an effectivefield-theory point of view, this lack of unitarity should not be viewed as an impediment for the phenomenological studies and experimental searches of generic spin-1 monopoles, given that the potential appearance of new degrees of freedom in the ultraviolet completion of such models might restore it. The second part of the paper deals with an appropriate implementation of photon-fusion and Drell-Yan processes based on the above theoretical scenarios into MadGraph UFO models, aimed to serve as a useful tool in interpretations of monopole searches at colliders such as LHC, especially for photon fusion, given that it has not been considered

Introduction
Eighty seven years since its concrete formulation by Dirac [1,2] as a quantum mechanical source of magnetic poles, the magnetic monopole remains a hypothetical particle. Although there are concrete field-theoretical models beyond the Standard Model (SM) of particle physics which contain concrete monopole solutions [3][4][5][6][7][8][9][10]13], these are extended objects with complicated substructure, and their production at collider is either impossible, as their mass range is beyond the capabilities of the latter [3,4,13], or extremely suppressed, due to their underlying composite nature [14]. On the other hand, point-like monopoles, originally envisaged by Dirac, are sources of singular magnetic fields for which the underlying theory, if any, is completely unknown, even though in principle (due to their point-like nature) they could avoid suppression in production.
In this respect, the spin of the monopole remains a free parameter. One may attempt to obtain an, admittedly heuristic, understanding of their production by considering effective field theoretic models for such production mechanisms based on electric-magnetic duality. That is deriving the corresponding cross sections from perturbative field-theoretical models describing the interaction of fields of various spins, S = 0, 1 /2 and 1 with photons, upon the replacement of the electric charge q e by the magnetic charge g, the latter obeying Dirac's quantisation rule where c is the speed of light in vacuo,h is the reduced Planck constant, 0 is the vacuum permittivity and Z is the set of integers (with n = 0 denoting the absence of magnetic charge). The quantity ξ depends on the system of units used, with ξ = 0 representing the CGS Gaussian system of units, and ξ = 1 the SI system of units. In natural SI units (h = c = 0 = 1), which we adopt here, the fine-structure constant at zero energy scales is given by α = e 2 4π = 1 137 with e > 0 the positron charge, from which (1) yields g = 1 2α n e q e e = 68.5e e q e n ≡ n e q e g D , n ∈ Z, with g D = 68.5e the fundamental Dirac charge. (We note, for completeness, that in composite monopole models [3][4][5], the magnetic charge can be viewed as a collective coupling of 1/α quanta of constituent (non-Abelian W ± -boson and Higgs) fields to a soft photon [14], which is consistent with the charge quantisation condition (2).) The electric-magnetic duality replacement, which obeys the quantisation rule (1), may be used as a basis for the evaluation of monopole-production cross sections from collisions of SM particles (quarks and leptons). Unfortunately, due to the large value of the magnetic charge (2), such a replacement renders the corresponding production process non-perturbative, consequently the strong-magneticcoupling limit dual theory is not well defined. Nevertheless, one may attempt to set benchmark scenarios for the cross sections by using tree-level Feynman-like graphs from such dual theories. This is standard practice in all point-like monopole searches at colliders so far [15].
Depending on the spin of the monopole field M, typical graphs participating in monopole-antimonopole pair production at LHC from proton-proton ( pp) collisions are given in Fig. 1. There are two kinds of such processes: the Drell-Yan (DY) (see Fig. 1b) and the photon-fusion (PF) induced production (see Fig. 1c, d). We also mention, to be complete, that in photon-photon production we have not only elastic but also semi-elastic and inelastic photon-fusion processes. For spin-1 /2 monopoles relevant in this discussion, a comparison between the respective perturbative cross sections has been provided first in Ref. [16,17] for pp collisions at a centre-of-mass energy √ s = 1.96 TeV, and subsequently for pp collisions of √ s up to 14 TeV in Ref. [18]. The conclusion from such analyses was that for √ s = 1.96 TeV the two cross sections are of comparable magnitude [18], whilst for √ s = 14 TeV PF dominates DY by a factor > 50, thus stressing the need to utilise the latter in monopole-search interpretations.
In this work we extend such combined (DY plus PF) studies to also incorporate spin-0 and spin-1 monopoles. In Ref. [16,17], the authors have calculated (in appropriate dualised models) the total cross sections for monopole pair production by photon fusion for three different spin models, spins 0, 1 /2 and 1. Theoretically, the cross sections increase with increasing spin, for a fixed (common) value of the monopole mass. It should be stressed that the expressions for the cross sections are specific to particular definitions of the monopole interactions. Specifically, the spin-1 /2 model fixes the particle in a minimally coupled theory, dual to standard QED, thus mirroring the observed behaviour of the electron with a gyromagnetic ratio g Re = 2. The corresponding magnetic moment κ is assumed zero. The magnetically charged monopole with spin 1 considered in [16,17], on the other hand, is characterised by a non zero magnetic moment term κ = 1, which is the value that characterises the charged W ± bosons in the SM. In fact the model mirrors the interactions of such bosons with a photon, but in the dual theory, where the electric charge is replaced by the magnetic charge. The value κ = 1 is the only one that respects unitarity [19]. For spin-0 monopoles, the dual theory of which resembles the Scalar Quantum Electrodynamics (SQED), no magnetic moment is allowed.
In this work, we generalise the discussion to include an arbitrary value for the magnetic dipole moment for monopoles with spin. The detailed reasoning for this is given in the next section. Hence, we shall treat κ as a new phenomenological parameter. 1 We note that, setting the issue of unitarity aside, phenomenological models of charged Wbosons interacting with photons with κ = 1 have been considered in the past [20], where it was demonstrated that the behaviour of the total cross section for the W -boson pair production for the κ = 1 case is quite distinct from the κ = 1 cases. In the current work, we shall dualise such models to use them as effective theories for the κ = 1 spin-1 monopole case, generalising the work of [16,17]. As we shall see, one may allow for some formal large-κ limit where, despite the strong magnetic coupling, the associated monopole-pair production cross sections can be made finite. In fact one may give meaning -under some circumstances to be specified below (specifically, the production of slow monopoles) -to the perturbative tree-level Feynman-like graphs of the effective theory. The relevant formalism for arbitrary κ and various monopole spins will then be used as a guide for the construction of appropriate MadGraph [21] algorithms that can be used as tools for data analyses in monopole searches at colliders. We remark for completeness that, for fast relativistic monopoles (characterised by a relative velocity β 1) Fig. 1 Feynman-like tree-level graphs for production processes of monopoles with generic spin S. a Typical SM Drell-Yan process describing charged lepton production from quark-antiquark annihilation; b DY monopole-antimonopole pair production from quark annihilation; c monopole-antimonopole pair production via photon fusion (for monopole spins 0, 1 /2 and 1); d additional (contact) diagram for monopole-antimonopole pair production via photon-fusion (for spins 0, and 1). The quantities q e and g denote the electric and magnetic charge, respectively passing through materials, the (large) number of electronpositron pairs produced can be used as a signal for the presence of the monopole [22]. The perturbativity conditions discussed in the present article, which pertain to slowly moving monopoles (with β 1), of relevance to MoEDAL-LHC searches [23,24], do not apply to such cases, the study of which requires non-perturbative treatments.
The structure of the article is the following: Sect. 2 is a review of the formal procedure to construct perturbative cross sections from the scattering amplitudes relevant to monopole production. Particular attention is given to a discussion of a rather unresolved current issue, regarding pointlike monopole production through scattering of SM particles, namely the use of an effective magnetic charge coupling that depends on the relative velocity β between the monopole and the centre-of-mass of the producing particles (quarks in the case of interest here, see Fig. 1). We also motivate the introduction of the magnetic dipole moment for monopoles with spin. In Sect. 3, the differential and total cross sections for monopole-antimonopole pair production via PF and DY processes are derived for various monopole spins. Combined limits of large magnetic-moment parameter κ and small β can lead to finite perturbatively valid results, providing some support for the effective formalism. The pertinent Feynman rules are implemented in a dedicated MadGraph model, which is described in detail in Sect. 4. In Sect. 5, the monopole phenomenology at the LHC is discussed, utilising the Mad-Graph UFO models developed in this work, in the context of the above theoretical considerations. Conclusions and outlook are given in Sect. 6, while details on the calculations are provided in appendices A and B.

From amplitudes to kinematic distributions
We are interested in the electromagnetic interactions of a monopole of spin S = 0, 1 /2, 1 with ordinary photons. The corresponding theory is an effective U (1) gauge theory which is obtained after appropriate dualisation of the pertinent field theories describing the interactions of charged fields of spin-S with photons. However, there is a subtle issue here, which we now proceed to discuss, and which will be relevant to our subsequent studies. It concerns a potential dependence of the effective magnetic charge on the relative velocity of the monopole pair and the centre-of-mass of the producing particles, as noted in [25,26].
To this end, we need first to recall some basic facts of the theory of point-like monopole-matter scattering. Classical (tree-level) scattering of charged particles off massive magnetic monopoles has been studied extensively, using modern quantum field theory (even relativistic) treatments [25][26][27][28]. As the monopoles of interest are relatively heavy compared to electrons, with masses of order at least TeV, the use of non-relativistic scattering suffices for our purposes in this section. The differential cross section for the classical (nonrelativistic) scattering of electrons, representing matter with electric charge e and mass m, off a magnetic monopole, with magnetic charge g and mass M, reads [25,27]: in a frame where the monopole is initially at rest and the electron has an initial velocity v 0 . In the above formula, c denotes the speed of light in vacuo and θ is the scattering angle, given in [25]: m+M is the reduced mass of the two body problem at hand, with μ m in the cases of monopoles where M m which is of interest here. The angle χ = 2 arccot (μv 0 b/|κ|) defines the (Poincaré) cone on which the (classical) trajectory of the electron in the background of the monopole is confined, with b the impact parameter. We note for completeness, that the cross section diverges in two occasions: (i) sin θ = 0 and sin χ = 0, θ = π , which occurs for a certain discrete set of cone angles [25,26], and (ii) when dθ dχ = 0, which occurs for a certain discrete set of scattering angles θ [25,26].
For small scattering angles, θ 1, the differential cross section (3) can be approximated by: The reader should notice that this exhibits a scaling with the inverse square power of the velocity v 0 , which is different from the standard Rutherford formula of the scattering of electrically-charged particles entailing a v −4 0 scaling.

Velocity-dependent magnetic charge
As can be readily seen from (4), this expression reduces to the standard Rutherford formula for electron-electron scattering upon the replacement This prompted some authors, including Milton, Schwinger and collaborators [25,26], to conjecture, upon invoking electric-magnetic duality, that when discussing the interaction of monopole with matter (electrons or quarks), e.g.
when discussing propagation of monopoles in matter media used for detection and capture of monopoles, or considering monopole-antimonopole pair production through DY or PF processes, a monopole-velocity dependent magnetic charge has to be considered in the corresponding cross section formulae: We stress again that the above substitution is based on the assumption of electric-magnetic duality, which would lead to the equivalence of the electron-monopole scattering crosssection (4) with the corresponding Rutherford formula upon applying (6). It is not known at present how to derive such an effective magnetic charge in the context of effective field theories (one might think of applying Schwinger-Dyson techniques which could resum the (large) magnetic charge couplings, but such theories are not available, and one cannot simply extend a strongly coupled QED model to the monopole case, upon replacing the electric charge with a magnetic one). Hence the substitution (6) should only be viewed at present as a conjecture, motivated by electricmagnetic duality symmetry. The replacement (6) was then used to interpret the experimental data in collider searches for magnetic monopoles [16][17][18]25,[29][30][31][32]. Due to the lack of a concrete theory for magnetic sources, the results of the pertinent experimental searches can be interpreted in terms of both a β-independent and a β-dependent magnetic charges, and then one may compare the corresponding bounds, as done in the recent searches by the MoEDAL Collaboration [32]. The monopole velocity β used in this work is given by the Lorentz invariant expression in terms of the monopole mass M and the Mandelstam variable s (A8), where s representing the square of the centreof-mass energy of the fusing incoming particles (photons or (anti)quarks) ( The Lagrangian for each effective model describes the propagation and interactions of a massive monopole field and a photon field. The field theory is chosen according to the spin of the monopole, S = 0, 1 /2 or 1. As a result of (6), the coupling of the monopole to the photon g(β) is linearly dependant on the particle boost, β = | p|/E p , where | p| and E p are the monopole's three-momentum and energy, respectively.
To incorporate in a unified way both the β-dependent and β-independent magnetic couplings of the monopole to photons, we define the magnetic fine structure constant as α g = g 2 β 2δ /(4π), where δ = 1(0) for β-(in)dependent cou-plings (6). The monopole Lagrangian of the effective theory can then be recognised as the Lagrangian in an electromagnetic field theory describing the interaction of a spin-S field with photons, with the following substitutions: From such a Lagrangian, Feynman rules can be extracted and observables are computed as in standard perturbative treatments. But, as stressed above, this is only a formal procedure, since, given the large value of the magnetic charge, as a consequence of the quantisation rule (2), the interaction coupling is in the non-perturbative regime. Hence, truncated processes, like DY or PF, have no meaning, unless, as we shall discuss below, certain formal limits are considered in some special cases. An important remark is in order at this point. Since the 'velocity' β (7) is expressed in terms of Lorentz-invariant Mandelstam variables, the Lorentz invariance of the effective field theory action of the monopole is not affected by the introduction of the effective β-dependent magnetic charge (8). However, there is a well known paradox, due to Weinberg [33], who pointed out that the amplitude for a single photon exchange between and electric and a magnetic current (of relevance to DY monopole-antimonopole production processes) is neither Lorentz nor gauge invariant due to the rôle played by the Dirac string, which contradicts the fact that monopoles appear as consistent soliton solutions in Lorentz and gauge invariant field theories [3,4]. It is only recently [34], that this paradox was arguably resolved, albeit within toy models of monopoles, coupled to photons, with perturbative electric and magnetic couplings. The resolution of Weinberg's paradox in such models is provided by a resummation of soft photons, which was possible due to the pertubatively small magnetic charges involved. Such a resummation resulted in the exponentiation of the Lorentz (and gauge) non-invariant terms pointed out in [33] to a (Bohm-Aharonov type) phase factor in the respective amplitude. In this sense, the modulus of the amplitude (and hence the associated physical observables, such as cross sections) are Lorentz (and gauge) invariant, with the important result that, upon Dirac quantisation (1), the (resummed over soft photons) amplitude itself is Lorentz invariant. It is in such Lorentz (and gauge) invariant frameworks that, as we conjecture, the considerations of the effective, velocity-dependent magnetic charge (8), may apply, which by the way also implies perturbative magnetic couplings for sufficiently small production velocities of the monopoles in the laboratory frame.
2.2 The magnetic dipole moment as a novel free parameter for monopoles with spin As already mentioned in the introductory Sect. 1, in previous effective-field theory treatments of spin-1 Dirac monopoles [16,17], a magnetic dipole moment with the value κ = 1 has been introduced mimicking the unitary SM case of W ± bosons (representing the monopoles), interacting with photons. Lacking an underlying microscopic model for point-like monopoles, the above restrictions in the value of the magnetic moment may not necessarily be applied to the monopole field. Indeed, for monopoles with spin, such a parameter may arise, e.g. by quantum corrections, in similar spirit to the electron case. The difference of course is that in the latter case it is the electric charge of the electron that would play a rôle, while in the monopole case it is the magnetic charge which, in view of its large value following the quantisation rule (2), cannot be treated perturbatively. Nonetheless, a non trivial (possibly large) magnetic moment might be induced in such a case, which might also be responsible for the restoration of unitarity of the effective theory. For example, one might hope that the apparent unitarity issues for generic κ = 1 values in the case of spin-1 monopoles can be remedied by embedding the corresponding theory in microscopic ultraviolet complete models beyond the SM, in much the same way as unitarity is restored in the case of the W ± gauge bosons interacting with photons in the SM case.
There is an additional reason as to why a non-trivial value for the parameter κ for the monopoles might be feasible. Although the full microscopic rôle of κ still needs to be determined, something which is obscured at present due to the lack of a microcsopic theory of point-like Dirac monopoles, nonetheless its presence is arguably consistent with the charge quantisation rule (1). This is due to the fact that, as we shall discuss, a magnetic dipole moment does not contribute to the singular part of the magnetic field of the monopole, which is responsible for the charge quantisation [27]. This can be argued as follows: as it is well known [27], Dirac's quantisation of charge in the presence of a monopole can be discussed by considering a particle of electric charge q e (e.g. electron) moving along a loop far away from the monopole centre, whose area is pierced by the Dirac string. Placing a magnetic dipole moment vector μ D at the origin of a coordinate system, where the monopole is assumed at rest, will induce, according to standard (classical) electromagnetism, a magnetic field where μ 0 is the magnetic permittivity of the vacuum ('free space'), the symbol denotes a unit vector, and (r, θ, φ) are the usual spherical polar coordinates. The formula is valid for large distances compared to the dipole longitudinal dimension. Writing B D = ∇ × A D , we can determine the corresponding vector potential as A D = μ 0 4π μ D × r r 3 , at large distances r . One can readily confirm that ∇ · B D = 0.
On the other hand, due to the monopole's magnetic charge, there is a magnetic field contribution, which however, due to the (singular) Dirac string, requires proper regularisation [27]. Upon doing so, one obtains for the regularised monopole magnetic field for a Dirac string along the z-axis, in which case the unit vector n = (0, 0, 1) also lies along that axis. The regularised form of the monopole's magnetic field intensity yields the correct formula ∇ · B reg monopole = 4π g δ (3) (r), implying that the magnetic monopole is the source of a field.
If one considers a charged particle looping the Dirac string far away from the position of the monopole, one would then observe that it is the singular part of the magnetic field (10), B sing , which contributes to the phase of the electron wavefunction [27], q e L d x · A = (L) dσ · B sing = 4πq e g. The magnetic dipole moment does not contribute to the singular part of the magnetic field, and thus the charge quantisation (1) is not affected. We note that, as a result of the r 3 suppression, the contributions of (9) would be subdominant, at large distances r from the monopole centre, as compared to those of (10).
It is worth remarking at this stage that one can also view [35,36] the quantisation rule (1) itself as a consequence of representing the (non-physical) singular string solenoid, assumed in the original Dirac's construction [1,2], as a collection of (small) fictitious current loops (with an area perpendicular to the solenoid's axis). Each one of these loops will induce a magnetic moment I A, with I the current and A the area of the loop (assumed vanishing in this case). Assuming a uniform magnetic moment per unit length M for the Dirac string, then, and taking into account that the solenoid may be viewed as the limiting case of a magnetic dipole of infinite length, one may apply the aforementioned formula (9) in this case to derive the singular magnetic field of the monopole itself, in which case the magnetic charge is obtained as g ∝ M [35,36]. However, we stress, that, the contribution of the induced magnetic moment of the quantum effective theory of a monopole with spin, whose strength depends on the parameter κ, which we shall discuss in this work, is independent of that due to the magnetic charge g, as explained above. Lacking though an underlying fundamental theory for the point-like monopole the determination of κ is at present not possible.
Before closing this section, we would like to present an equivalent, yet less elaborate, way 2 to see the irrelevance of the magnetic dipole moment for the quantisation rule (1), which avoids the use of Dirac strings. To this end, one covers the three-space surrounding the monopole by two hemispheres, with appropriate gauge potentials defined in each of them, whose curl yields the corresponding magnetic field strengths. For the magnetic monopole gauge potential one has the expressions [27]: for the south hemisphere, which is singular at the south pole θ = π , and for the north hemisphere, which is singular at the north pole θ = 0. These two patches overlap π 2 −δ < θ < π 2 +δ, δ → 0 + , and, as is well known, the difference of yields a singular gauge transformation at θ = 0, π, which contributes to the phase q e L d x · A of the charged particle wavefunction, the requirement of single-valuedness of which yields the rule (1).
On the other hand, as already mentioned, the vector potential corresponding to the magnetic moment, for large distances r from the centre of the sphere where the monopole is located, is of the form with η the unit vector perpendicular to the plane of r and μ D (assumed parallel to the z-axis); this is not singular at the poles θ = 0, π (in fact it vanishes there). The total potential in each hemisphere is then given by the corresponding sum A i + A D , i = S, N . Hence, the magnetic moment does not contribute to the difference, and thus it does not affect the wavefunction phase, which is associated only with the monopole part (13).

Cross sections for spin-S monopole production
In this section we derive the pertinent Feynman rules and then proceed to give expressions for the associated differential and total production cross sections for monopole fields of various spins. The pertinent expressions are evaluated using the package FeynCalc [37,38] in Mathematica. We commence the discussion with the well-studied cases of scalar (spin-0) and fermion (spin-1 /2) monopole cases, but extend the fermion-monopole case to include an arbitrary magnetic moment term κ = 0. Then we proceed to discuss the less studied case of a spin-1 monopole including an arbitrary magnetic moment term κ. We consider both β-dependent and β-independent magnetic couplings. In an attempt to make some sense of the perturbative estimates, we discuss, where appropriate, various formal limits of weak β and large κ for which the pertinent cross sections remain finite. In each spin case we present both PF and DY cross sections, which will help us present a comparison at the end of the section.

Scalar monopole
The first model studied in this work, is the one for massive spin-0 monopole interacting with a massless U (1) gauge field representing the photon. Searches at the LHC [31,32,[39][40][41] and other colliders have set upper cross-section limits is this scenario assuming Drell-Yan production. The Lagrangian describing the electromagnetic interactions of the monopole is given simply by a dualisation of the SQED Lagrangian where D μ = ∂ μ − ig(β)A μ , A μ is the photon field, whose field strength (Maxwell) tensor is F μν = ∂ μ A ν − ∂ ν A μ and φ is the scalar monopole field. There are two interaction vertices associated with theory. The three-and four-point vertices are illustrated in Fig. 29 in "Appendix B", where the cross-section calculations are detailed. These interactions are the only couplings generated between the spin-0 monopole and the U (1) gauge field.

Pair production of spin-0 monopoles via photon fusion
There are three possible graphs contributing to scalar monopole production by PF, a t-channel, u-channel and seagull graph shown in Fig. 2. Their respective matrix amplitudes are given by Eq. (B2) in "Appendix B". M is the spin-0 boson mass, ε λ (q 1 ) and ε λ (q 2 ) are the photon polarisations, p 1 and p 2 are the monopoles four-momenta such that p 2 i μ = M 2 , and q 1 and q 2 are the photons four-momenta such that q 2 i μ = 0, as defined in Fig. 2.
After some calculation, detailed in "Appendix B", the differential cross section for the spin-0 monopole-antimonopole production is reduced to: which in terms of the pseudorapidity η, defined in "Appendix A", becomes: The differential distributions as a function of the scattering angle and the pseudorapidity are shown in Fig. 3 for a hypothetical monopole of mass M = 1.5 TeV at The pair production is mostly central as is the case for various beyond-SM scenarios. Such distributions are used for the validation of the simulation package discussed in Sect. 4.
After integrating over the solid angle, as discussed in Sect. 2, the total cross section becomes [16,17] and is shown in Fig. 4 for the selected energy of √ s γ γ = 4 TeV. The cross section drops rapidly with increasing mass and disappears sharply at the kinematically forbidden limit of M > √ s γ γ /2. Again this result is compared against the MadGraph implementation prediction in Sect. 4.

Pair production of spin-0 monopoles via Drell-Yan
The Feynman-like diagram in Fig. 5 shows the DY process in the case of a scalar monopole. The quarks qq annihilate to a photon A π , which decays to a M M pair in the s-channel.
The quark lines are supplemented by momentum 4-vectors q 1μ and q 2μ , where q 2 1,2 = m 2 and the scalar monopole lines have momentum 4-vectors p 1μ and p 2μ , where p 2 1,2 = M 2 on shell. The centre-of-mass energy of the colliding quarks is k π k π = s qq . The three-point vertex in this model is illustrated again in Fig. 30a, and the vertex for the qqA μ coupling in Fig. 30b.  After some calculations detailed in "Appendix B" and in particular in Eqs. (B9)-(B15), the differential cross section for DY scalar monopole production yields where β is defined in Eq. (7). The distribution is shown in Fig. 6 for a monopole with a mass M = 1.5 TeV at a centreof-mass energy √ s qq = 2E q where E q = 6M. In terms of the pseudorapidity η, the differential cross section reads: and is also plotted in Fig. 6. As expected for the production of scalar particles, the distribution is almost flat with respect to the scattering angle θ . When compared to the corresponding kinematic distributions from the PF case in Fig. 3, the DY case exhibits a more central distribution. Finally, the total cross section is evaluated by integrating Eq. (19) over the solid angle d = dφ dcos θ : and is shown graphically at a centre-of-mass energy √ s qq = 4 TeV in Fig. 7. A note of caution here: in an experimental setup for high energy collisions, the result (21) is valid for opposite-sign hadrons, e.g. pp at the Tevatron; it will be doubled in the case of same-hadron colliding bunches, such as in proton-proton collisions at the LHC. A comparison of the DY versus the PF cases, as far as the differential/total cross sections is concerned, is reserved for Sect. 3.4.

Spin-1 /2 point-like monopole with arbitrary magnetic moment term
The phenomenology of a monopole with spin 1 /2, examined in this section, is the most thoroughly studied case [16- , however so far only the Drell-Yan production process has been explored in collider searches [31,32,[39][40][41]. This type of monopole resembles a magnetic dual to the electron. The electromagnetic interactions of the monopole with photons, are described by model U (1) gauge theory for a spinor field ψ representing the monopole interacting with the massless U (1) gauge field A μ , representing the photon. In the cases discussed in the existing literature, the effective Lagrangian describing the interactions of the spinor monopole with photons is taken from standard QED upon imposing electric-magnetic duality, in which there is no bare magnetic-moment term (at tree level). However for our analysis here, we shall insert in the Lagrangian a magnetic moment generating term, to keep the treatment general. The origin of the magnetic moment of the monopole is not known, so it would be naïve to assume it is generated only through anomalous quantumlevel spin interactions as for the electron in QED. 3 In the event κ = 0, the Dirac Lagrangian is recovered. Thus, the effective Lagrangian for the spinor-monopole-photon interactions takes the form where F μν is the electromagnetic field strength tensor, the total derivative is / is a commutator of γ matrices. The magnetic coupling g(β) is given in (8), and is at most (depending on the case considered) linearly dependent on the monopole boost, β = | p|/E p , where | p| and E p are the monopole momentum and energy, respectively. The effect of the magnetic-moment term is observable through its influence on the magnetic moment at tree level which is where M is the spinor-monopole mass,Ŝ is the spin expectation value and the corresponding "gyromagnetic ratio" Noticeably, the parameter κ in (23) is not dimensionless, but has units of inverse mass which breaks the renormalisability of the theory. This may not be a serious obstacle for considering the case κ = 0, if the pertinent model is considered in the context of an effective field theory embedded in some yet unknown microscopic theory, in which the renormalisation can be recovered. 4 3 For instance, it is known that such terms have a geometrical (gravitational) origin in 4-dimensional effective field theories obtained from (Kaluza-Klein) compactification of higher-dimensional theories, such as brane/string universes [42]. Moreover, as mentioned in the introduction, one could also include in the Lagrangian a CP-violating electric dipole moment (EDM) term for the spinor monopole, parametrised by a parameter η, In this work, we assume such EDM terms suppressed, compared to the magneticmoment-κ terms, which can be arranged by assuming appropriate limits of parameters in the underlying microscopic theory, e.g. [42]. Nonetheless, our analysis can be extended appropriately to include both κ and η parameters. 4 In case one adds an EDM η-term for the spinor monopole, the corresponding dimensionless parameterη can also be defined in analogy with (25), i.e. η =η M .

Spinor pair production via photon fusion
The monopole couples to the photon at a three-particle vertex through a Feynman rule that is shown in Fig. 31 in "Appendix B" together with details of the amplitudes that lead to the cross-section calculation. The parameter κ influences the amplitude in the second term of (B17), ensuring the κ-dependence of the observables. PF occurs through tchannel and u-channel processes as shown in Fig. 8, with the κ-dependent matrix amplitudes stated in Eq. (B16) of "Appendix B". The photon momenta are q 1 , q 2 , the monopole momenta are p 1 , p 2 , whilst k,k are the t-and u-channel exchange momenta, respectively.
The differential cross section distributions are computed as described in Sect. 2, with more detail given in "Appendix B". The κ-dependent forms are given below: with the standard (dual QED) case κ = 0 being given in Refs. [16][17][18]: The angular distributions for monopole-antimonopole pair production by PF are shown in Fig. 9 for spin 1 /2 and for various values of the parameterκ. The monopole mass is set to M = 1.5 TeV and the photon energy is E γ = 6M. This is the expectation if the monopole is truly dual to the electron in that it is in every way just the magnetic counterpart of well-known electric sources. The distribution shape is quite distinct from that of the scalar-monopole case (see Fig. 3) exhibiting a depression at η 0 for allκ values. The casẽ κ = 0 represents the SM expectation for electron-positron pair production if the coupling is substituted for the electric charge e, i.e. restoring the Lagrangian to simple Dirac QED. This case is clearly distinctive as the only unitary and renormalisable case. It is observed that the vertical heights of the curves for the differential cross section change with κ. Furthermore, the κ = 1 and κ = −1 cases are totally equivalent, so a degeneration between positive and negative κ values is evident. The total cross section for arbitrary κ, also shown graphically in Fig. 10, is given by: Setting κ = 0, the expression (28) reduces to that given in Refs. [16][17][18]30], used in standard monopole searches at colliders for data interpretation [31,32,[39][40][41]. The high-energy limit for (26), expressed by s γ γ → ∞, leads to the approximations β n → 1 for n > 2 and β 1 − 2M 2 s γ γ . In this limit, we observe that the differential cross section (26) diverges, except, unsurprisingly, when recovering the Dirac model with κ = 0: × κ 4 s γ γ cos 4 θ + 6κ 4 s γ γ cos 2 θ + κ 4 s γ γ +8κ 2 cos 2 θ + 28κ 2 + 1 s γ γ 4 cos 2 θ + 1 The total cross section (28) carries the same high energy behaviour. It is finite for κ = 0, and diverges for all other values of κ: For a constant value of s γ γ , as assumed in Fig. 10, the highenergy limit may be approximated by M → 0, where the cross section becomes finite only for the value κ = 0 and diverges in other κ values, as expected from (30). Hence, a unitarity requirement may isolate the κ = 0 model from the others as the only viable theory for the spin-1 /2 monopole, unless the model is viewed as an effective field theory. In that case, the value of κ can be used as a window to extrapolate some characteristics of the extended model in which unitarity is restored. Also, as already mentioned, κ is not dimensionless, hence, a non-zero κ clearly makes this (effective) theory non-renormalisable.

Pair production of spin-1 /2 monopoles via Drell-Yan
Monopole pair production through the s-channel is also possible for fermionic monopoles through the annihilation of quarks into a photon, which decays to the monopoleantimonopole pair. The relevant Feynman rules are displayed in Fig. 32 with the κ-dependent matrix amplitude given in Eq. (B20), both in "Appendix B". The complete DY process  Fig. 11. The exchange energy in the centre-ofmass frame is k π k π = s qq .
The kinematic distributions are computed with details given in "Appendix B" and are shown in Fig. 12. The monopole is treated as the magnetic dual to the electron. Analytically, the differential cross section becomes: where β q = 1 on the right as the quarks are assumed to be massless (compared to the heavy monopoles), and the index has been dropped on β p → β. Still, the magnetic coupling depends on β as α g (β) ∝ β 2 in the velocity-dependent magnetic charge case. As observed in the previous section, the value of κ affects the unitarity and renormalisability of the model. In the spin-1 /2 model, this parameter has dimensions, destroying the renormalisability except when κ = 0, in which case the magnetic moment is only a by-product of anomalous spin interactions within a minimally coupling field theory. But unitarity for this process is maintained for all values of κ. This becomes apparent in the high energy limit s qq → ∞. Taking the expression (30) to first order in β in this limit, we observe that it is finite for all κ, The total cross section follows the same trend. Its full form is which is displayed graphically in Fig. 13. In the high-energy limit, it reduces to which is finite for all values of κ. In Fig. 13, where s γ γ is constant, the high-energy limit is approximated at M → 0, where the cross section becomes finite for the value κ = 0 and diverges as ∼ M −2 in other κ values, bearing in mind that κ is related to the monopole mass as in (25). It should be noted however that DY alone is not the only mechanism for producing ψψ pairs and on its own cannot define the unitarity of the theory as a whole. In the previous section, it was shown how this model violates unitarity when considering pair production by PF. It is worth noting that monopole production in a high energy collider sees twice the production cross section for collisions with same-sign incoming beams, such as protonproton collisions at the LHC, in contrast to opposite-sign hadron colliders, such as the Tevatron, which maintain the exact cross section as given in (32).

Vector monopole with arbitrary magnetic moment term
Monopoles of spin-1 have been addressed for the first time in colliders recently by the MoEDAL experiment for the Drell-Yan production [32]. A monopole with a spin S = 1 is postulated as a massive vector meson W μ interacting only with a massless gauge field A μ in the context of a gauge invariant Proca field theory. As mentioned previously, lacking a fundamental theory for point-like magnetic poles, we keep the treatment general by including a magnetic moment term in the effective Lagrangian, proportional to κ, which is a free phenomenological parameter. Unlike the spin-1 /2 monopole case, however, for the vector monopole the magnetic moment parameter κ is dimensionless. The case κ = 0 corresponds to a pure Proca Lagrangian, and κ = 1 to that of the SM W μ boson in a Yang-Mills theory with spontaneous symmetry breaking. In this respect, our approach resembles early phenomenological studies of charged W ± -boson production in the SM through PF, where the magnetic moment of the W -boson was kept free [20], different from the value κ = 1 dictated by unitarity. The aim of such analyses was to determine measurable (physical) quantities in purely electromagnetic SM processes, that were sensitive to the value of κ, and more or less independent of the Higgs field and the neutral gauge boson Z 0 . These quantities were the angular distributions at sufficiently high energies, whose behaviour for the unitarity-imposed value κ = 1 was found to be quite distinct from the case κ = 1. As we shall see in our case, for certain formal limits of large κ and slowly moving monopoles, one may also attempt to make sense of the perturbative DY or PF processes of monopole-antimonopole pair production, when velocity-dependent magnetic charges are employed. The pertinent effective Lagrangian, obtained by imposing electric-magnetic duality on the respective Lagrangian terms for the interaction of W ± gauge bosons with photons in the generalised SM framework, as described above [20], takes the form: where the symbol † denotes the hermitian conjugate, and 1) covariant derivative, which provides the coupling of the (magnetically charged) vector field W μ to the gauge field A μ , playing the role of the ordinary photon. The parameter ξ is a gauge-fixing parameter. The magnetic coupling is considered in the general form of (8), so as to cover both the β-dependent and β-independent cases in a unified formalism. The tensor F μν represents the Abelian electromagnetic field strength (Maxwell). The magnetic and quadrupole moments are given respectively by whereŜ is the monopole spin expectation value and the corresponding "gyromagnetic ratio" g R = 1 + κ. The phenomenological moment term −ig(β)κ F μν W † μ W ν is highly divergent and contributes correction terms to the magnetic and quadruple moments in Eqs. (35a) and (35b), respectively. As mentioned already, following Ref. [20], we treat κ as a free phenomenological parameter of the theory. 5 For κ = 1, the theory is known to be non unitary, and is plagued by ultraviolet divergences in the self-energy loop graphs in this model, making the quadrupole moment infinite in a non-renormalisable way. To tackle such divergences, in the pre-SM era, Lee and Yang [19] proposed the effective Lagrangian (34), and demonstrated that such divergences are removed through the inclusion of the gauge fixing term with the gauge fixing parameter ξ = 0, but at a cost of introducing a negative metric (and thus ghosts, reflecting the unitarity issue for κ = 1). In this "ξ -limiting formalism", as it is called, the observables are evaluated from the ξ -dependent Lagrangian before taking the limit ξ → 0. The quadruple moment becomes finite at one-loop level [20]. Unitarity in this formalism is held only at energy scales E 2 ≤ M 2 /ξ , the rest mass energy of a single ghost state. For our purposes, this could be an acceptable assumption for an effective field theory considered valid up to a cut-off scale 2 = M 2 /ξ . In general, lacking a concrete fundamental theory on magnetic poles, we shall ignore the unitarity issue when we consider the incorporation of an arbitrary magnetic moment κ in our construction.
At this point, we should also mention that in general, there is another unitarity issue that arises within the context of the PF or DY cross sections in the dual effective theories we consider here. Even in the unitary κ = 1 case inspired by SM physics, the associated cross section for large values of β violates the unitarity bound for a cross section dominated by a single partial wave of angular momentum J : For sufficiently small values of β on the other hand, such unitarity bounds are respected. In fact this is a feature that characterises the cross sections of all three cases of monopole spin S = 0, 1 /2, 1 and not only the vector case [16,17]. To tackle such an issue, within a phenomenological effective model for the monopoles, one may assume the existence of appropriate form factors that depend on the energy of the photon-monopole interaction [16,17]. We shall not pursue such issues further in this section, but we mention that such 5 Corrections to the magnetic moment of the monopole could also arise through anomalous spin interactions at a quantum level. For βindependent magnetic charges, these are uncontrollable, as perturbation theory fails. However, if one accepts velocity-dependent couplings (5), then for slowly-moving monopoles such loop corrections can be made subleading. Moreover, as for the case of spinor monopoles, one could also add an EDM term for the vector monopole [43], μνρσ F ρσ , is the dual Maxwell tensor, with μνρσ the totally antisymmetric Levi-Civita tensor in four space-time dimensions. Such terms are assumed suppressed in our analysis, although the latter can be straightforwardly extended to include them. form factors will be included in the MadGraph generator simulating monopole production at colliders, as we shall discuss in Sect. 5.
After this parenthesis, we come back to the spin-1 monopole-production process via PF within the context of the model (34). Restricting our attention to tree level, the gauge fixing parameter ξ is redundant and the ξ -independent interaction vertices are given in Fig. 33. These are used to evaluate the κ-dependent PF Born amplitudes, As already mentioned, we introduce a velocity-dependent magnetic coupling g = g(β) corresponding to a magnetic charge linearly dependent on the monopole boost β = | p|/E p where p and E p are the monopole momentum and energy, respectively.

Pair production of spin-1 monopoles via photon fusion
Monopole-antimonopole pairs are generated at tree level by photons fusing in the t-channel, u-channel and at a 4-point vertex depicted by the Feynman-like graphs in Fig. 14. The matrix amplitude for each process is given in Eq. (B24) in Fig. 33 in "Appendix B". k andk are the exchange momenta of the t-and u-channel processes, respectively. The monopoles have polarisation vectors ϒ( p 1 ) κ , ϒ(p 2 ) κ and momentum 4-vectors p 1μ and p 2μ , where (on massshell): p 2 1,2 = M 2 . The photons polarisation vectors are (q 1 ) λ , (q 2 ) λ and their momentum 4-vectors are q 1μ and q 2μ , q 2 1,2 = 0. Details of the calculations of the analytic expressions for the kinematic distributions and cross section are given in "Appendix B".
The phenomenological parameter κ enters in the expressions of the differential cross sections given in (37). These kinematic distributions are plotted as functions of the kinematic variables θ and η in Fig. 15 assuming a monopole mass M = 1.5 TeV and centre-of-mass energy √ s γ γ = 2E γ with E γ = 6M. The monopole boost β is defined in (7), where the centre-of-mass energy is understood to be the energy of the fusing photons, i.e.
− 2β 2 11κ 4 + 60κ 3 + 58κ 2 + 12κ + 35 The kinematic distributions change as the parameter κ is varied. The reader can readily see from Fig. 15 the distinct behaviour of the kinematic distribution in the unitarityrespecting case κ = 1, that shows a depression around η = 0, as compared to the peaks in the cases where κ = 1. This is in agreement with the situation characterising W + W − production in the SM case [20]. However as we shall see in Sect. 5.2, when the PDF of the photon in the proton is taken into account, these shape differences are smoothed out in pp collisions including the (discernible here) κ = 1 case. For κ = 1 the expression (37) becomes: with β given in (7) and α g (β) in (8).
As mentioned previously, the parameter κ influences the unitarity and renormalisability of the effective theory. It is important to notice that the differential cross section in the κ = 1 case (38) is the only finite solution in the ultraviolet limit s γ γ → ∞. Indeed, in the high energy limit, s γ γ → ∞, one may approximate β 4 1 and β 1 − 2M 2 s γ γ , implying that the angular distribution (38) falls off as s −1 γ γ : For all other κ values one obtains a differential cross section proportional to s γ γ , which diverges linearly with s γ γ → ∞: The total cross section for general κ is given by which for the κ = 1 case reduces to [16,17] The total cross section is shown graphically in Fig. 16 for various κ values. In the high energy limit s γ γ → ∞, only the total cross section for κ = 1 is finite, in similar spirit to the differential cross section behaviour: We stress once again that the κ = 1 case is the SM result for vector boson scattering with photons in which the coupling term −igF μν W † μ W ν naturally arises in SU (2) ×U (1) gauge theory with spontaneous symmetry breaking. This value for κ also restores renormalisability and unitarity in the absence  The quarks annihilate to a photon through an electromagnetic process, which subsequently decays to a monopoleantimonopole pair. The variable definitions are given in the text of the negative metric and the ξ -gauge fixing. The meson adopts a (tree-level) gyromagnetic ratio g R = 2, which is the value associated with the W ± gauge boson of the SM.

Pair production of spin-1 monopoles via Drell-Yan
As discussed previously, another mechanism contributing to the production of monopole-antimonopole pairs is quarkantiquark annihilation through the s-channel, also known as Drell-Yan, drawn in Fig. 17, for which the relevant Feynman rules are given in Fig. 34. The quarks each have a mass m considered small compared to the monopole mass, M, and are characterised by momentum 4-vectors q 1μ and q 2μ , where on mass shell one has q 2 1,2 = m 2 . Similarly, the mesons have mass M each and are characterised by momentum 4-vectors p 1μ and p 2μ , where on mass-shell one has p 2 1,2 = M 2 . The centre-of-mass energy of the quark-antiquark pair is k ν k ν = s qq .
The differential cross section distributions are computed as defined in Sect. 2. The expression for β (7) is formally valid here as well, with the understanding that s = s qq now represents the Mandelstam variable for the initial quark-antiquark pair (see Fig. 17). In the approximation of negligible quark masses when compared to the monopole mass M, the differential cross section reads: Neglecting quark masses set β q = 1 so the index on the monopole boost β p has been dropped. Also, the β dependance of the magnetic coupling is made apparent.
The kinematic distributions are plotted in Fig. 18 for various values of the parameter κ with β q = 1 as expected for massless quarks. As in previous spin models, a monopole mass of M = 1.5 TeV and √ s qq = 2E q with E q = 6M is assumed. This parameter influences the shape of the distributions and the convergence of the cross section as s qq → ∞. Indeed, where the limit was taken such that β 4 1, β 3 1 and keeping β 2 = 1 − 4M 2 s qq . The total cross section for arbitrary values of κ is given by: and it is plotted as a function of the monopole mass in Fig. 19.
In the high energy limit, it becomes which converges only for κ = 0. In Fig. 13, we observe the κ-dependent behaviour at M → 0 as expected from (45): the and It should be noted that the unitarity of the κ = 1 case becomes apparent only if other SM processes are included in the total amplitude, hence the divergence seen for κ = 1 in the high energy limit (cf. (44), (45)) is not surprising. We also remark that the DY cross section increases by a factor of two when considering processes generated by high energy collisions of identical hadrons, e.g. protons to protons.

Comparison of various spin models and production processes
This section is brought to a close by briefly comparing the cross section distributions discussed in this work. Firstly, it is important to point out that the dominant production process is PF by a large margin at √ s qq/γ γ = 4 TeV. This is seen in Fig. 20 for all spin models. In particular, graphs 20a and 20b show the SM-like cases for which κ = 1 andκ = 0 represent the S = 1 and S = 1 /2 SM-like cases, respectively. Graph 20(c) shows the spin-0 monopole case, the only one for which there is no magnetic moment. This corroborates the assertion in [18]. Graphs 20d and 20e demonstrate that this behaviour is maintained for all values of κ andκ by choosing the non-distinctive value equal to two.
Secondly, it is apparent that the cross section for monopole production increases with the spin of the monopole most of the mass range, as observed in Fig. 21, if the SM-like cases for the magnetic-moment parameters are chosen. This cases, the SM valuesκ = 0 and κ = 1, respectively, are drawn, while there is no magnetic moment in the spin-0 case observation supports the findings of Ref. [18]. As shown in Figs. 10 and 13 for a fermionic monopole, the trend is maintained forκ > 0 for all masses. For a vector monopole, on the other hand, the cross-section ordering is not consistent across the monopole mass for varying κ values, as evident from Figs. 16 and 19. More discussion on the phenomenological implications of the magnetic-moment parameter will follow in Sect. 5.2, this time in the context of proton-proton collisions.

Perturbatively consistent limiting case of large κ and small β
As discussed in Sect. 1, the non-perturbative nature of the large magnetic Dirac charge of the monopole invalidate any perturbative treatment based on Drell-Yan calculations of the pertinent cross sections and hence any result based on the latter is only indicative, due to the lack of any other concrete theoretical treatment. This situation may be resolved if thermal production in heavy-ion collisions -that does not rely on perturbation theory -is considered [44][45][46]. Another approach is discussed here involving a specific limit of the parameters κ and β of the effective models of vector and spinor monopoles, used above, in the case of a velocity-dependent magnetic charge (8). In this limit, the perturbative truncation of the monopole pair production processes, described by the Feynman-like graphs of Figs. 8, 11, 14 and 17, becomes meaningful provided the monopoles are slowly moving, that is β 1. In terms of the centre-of-mass energy √ s γ γ /qq , such a condition on β implies, on account of Eq. (7), that the monopole mass is around 2M √ s γ γ /qq + O(β 2 ). It should be noted at this point that, in collider production of monopole-antimonopole pairs considered in this work, s γ γ /qq is not definite but follows a distribution, according to the parton (or photon) distribution function (PDF) for the DY or PF processes.
In the absence of a magnetic-moment parameter, κ, or for the unitary value κ = 1 in the case of spin-1 monopoles studied in [16,17], the condition β 1 would lead to strong suppression of the pertinent cross sections beyond the current experimental sensitivities, thereby rendering the limit β → 0 experimentally irrelevant for placing bounds on monopole masses or magnetic charges. Indeed, the various total cross sections discussed so far behave as follows, when β → 0 (using the definition (8) of the magnetic fine structure constant): However, in the case of non-trivial and large (dimensionless) magnetic-moment-related parameters κ,κ, relevant for the cases of vector and spinor monopoles, the situation changes drastically, as we shall now argue. To this end, we consider the limits withκ defined in (25), but in such a way that the strength of the derivative magnetic-moment couplings, given in Eq. (B17) (see Fig. 31) for spin-1 /2 monopoles and on the left side of (B25) (see Fig. 33) for spin-1 monopoles is perturbatively small. Since the magnitude of the monopole momentum is proportional to Mβ, one expects that the condition of a perturbatively small derivative coupling for magnetic moment is guaranteed if, by order of magnitude, one has: where κ =κ for spin-1 /2 monopole and κ = κ for spin-1 monopole.
For the spin-1 /2 monopole, from (28) and (32), one observes that in the limit (49) and respecting (50), upon requiring the absence of infrared divergences in the cross sections as β → 0, and postulating that: so that (50) is trivially satisfied, sinceκ gβ 2 = |c 1 | and respectively, where we used (8) and (25). Hence for slowlymoving spinor-monopoles, with velocity-dependent magnetic charge, and large magnetic moment parameters, it is the PF cross section which is the dominant one of relevance to collider experiments. Similar results characterise the spin-1 monopole. Indeed, we observe that in the limit (49), (50), the dominant contributions to the total cross sections for the PF (see (41)) and DY (see (44)) processes are such that and respectively, where we used (8).
We can see that, by requiring the absence of infrared (β → 0) divergences in the total cross sections, one may consistently arrange that the PF cross section (54) acquires a non-zero (finite) value as β → 0, whilst the DY cross section (55) vanishes in this limit: In such a limit, the quantity κgβ 2 = |c 1 |  (50) is trivially satisfied, and thus the perturbative nature of the magnetic moment coupling is guaranteed. Hence in this limiting case of velocity-dependent magnetic charge, large magnetic moment couplings and slowly moving vector monopoles, again the PF cross section is the dominant one relevant to searches in current colliders and can be relatively large (depending on the value of the phenomenological parameter c 1 ). This argument is successfully tested with simulated events in Sect. 5.2.

MADGRAPH implementation
The MadGraph generator [21] is used to simulate the generation of monopoles. In this section, we briefly present the development of the MadGraph Universal FeynRules Output (UFO) model [47] used to simulate different production mechanisms of monopole. This includes both the monopole velocity (β) dependent and independent photon-monopolemonopole coupling. Three different spin cases have been included: spins 0, 1 /2 and 1.

Monopole couplings
In Dirac's model, the relation between the elementary electric charge q e and the basic magnetic charge g is given in Eq. (2) in Gaussian units. However in MadGraph, Heaviside-Lorentz units are used, where (2) becomes q e g = 2π n, n ∈ Z. (57) Hence, the unit of the magnetic charge is The electromagnetic vertex in Heaviside-Lorentz units simply becomes Similarly, the monopole-photon vertex becomes In Heaviside-Lorentz units, electric charge q e is given by √ 4πα where α is the fine-structure constant. Hence, Eq. (58) turns out to be In (61), the monopole velocity β is not used. But if we want to consider the monopole velocity dependent coupling, the value of g D simply becomes β √ π/α. The velocity β is expressed in (7) as a function of monopole mass M.

Implementation of the monopole Lagrangians in MadGraph
In this section, the practical details on the use of Mad-Graph [21] to simulate the photon-fusion production mechanism of monopoles is described. The MadGraph was downloaded and installed following the instructions given in [48]. The general procedure to simulate a model with the help of MadGraph is the following: 1. Create a model 6 with all the user defined fields, parameters and interactions. Lately, the use of the UFO format [47] is strongly encouraged for such models. 2. In the MadGraph command prompt, import that model. 3. Generate the process which will be simulated using the generate command and create an output folder. 4. Fix the centre-of-mass energy, colliding particles, parton distribution functions in the run card. 5. Fix the parameters (electric and magnetic charges, masses, etc) of the colliding and generated particles in the parameter card. 6. Launch the output folder in order to compile the model and create the Les Houches Event (LHE) files [49]. 7. These LHE files will be used to produce simulated results.

Generating and validating the UFO models
To generate a UFO model from the Lagrangian, Feyn-Rules [50], an interface to Mathematica, was utilised.
Here the parameters of a model (mass of a particle, spin, electric charge, magnetic charge, coupling constant, fermionic or bosonic field, etc) and the corresponding Lagrangians are written in a text file. From the Lagrangian, the UFO model is generated with the help of FeynRules. The velocity β, defined in Eq. (7), is defined as a form factor in the generated UFO models. The instructions given in [51], specially the 'Method 2: Fortran Way' were followed to get a proper form factor. To get the value ofŝ inside the Fortran function, we used this formula (for elastic collision): where P 1 and P 2 are the 4-momenta of the two colliding particles. For scalar monopoles, the inclusion in the simulation of the four-particle vertex shown in Fig. 2c in addition to the u-and t-channel, shown in Figs. 2a, b respectively, led to the necessary use of UFO model written as a Python object and abandon the rather older method in Fortran code. The implementation of the four-vertex diagram proved to be nontrivial due to the g 2 coupling. The Lagrangian, which takes the form given in (15), is rewritten in a Mathematica format so that FeynRules can understand the variables. We created the text file containing all the information related to the field, mass, spin and charge, etc, following the instructions given in Ref. [52].
To validate the MadGraph UFO model for monopoles, we compare the cross sections for the photon-fusion process from the theoretical calculation derived in Sect. 3 to those obtained from simulation. Since the theoretical calculations consider bare photon-to-photon scattering, we chose in MadGraph the no-PDF option, i.e. we assume direct γ γ collisions at √ s γ γ = 13 TeV. Also, the coupling used here is assumed to depend on β.
The cross-section values for spin-0 monopoles are shown in the first columns of Table 1. The UFO-model-over-theory ratio values, also shown in the fourth column of the table, are very close to unity. This clearly shows the validity of spin-0 monopole UFO model.
In a similar fashion, spin-1 /2 monopole Lagrangians (23) are also rewritten in a Mathematica format. The magneticmoment parameterκ is also implemented in the model. No additional diagram was added to the u/t-channels already described in the UFO model. Again, the cross sections from theoretical calculations and MadGraph UFO models (for no PDF) for spin-1 /2 monopoles were compared, are shown in Table 1. The comparison clearly shows the validity of the MadGraph UFO model for spin-1 /2 monopoles.
Finally, the Lagrangian (34) for spin-1 monopoles is also written in Mathematica code. The possibility for choosing the value of the κ parameter in (34) exists, yet for validation purposes, the value of ξ is taken to be zero and the value of κ is taken to be one. The cross sections for spin-1 monopoles from the theoretical calculations and MadGraph UFO models (for no PDF), shown in Table 1, match which satisfactorily proved the validity of the MadGraph UFO model for the spin-1 monopole.
Apart from the photon-fusion production mechanism, the UFO models were also tested for the DY production mechanism as well. The DY process for monopoles was already implemented in MadGraph both for β-dependent and β-independent coupling using a Fortran code setup for spin-0 and spin-1 /2 monopoles. This setup was utilised by ATLAS [31,39] and MoEDAL [32,40,41] to interpret their search results in terms of monopoles under these assumptions. These Fortran setups were rewritten in the context of this work as UFO models, following their PF counterparts, and they were extended to include the spin-1 case. The latter was used in the latest MoEDAL monopole-search analysis [32]. After validating the UFO models against their Fortran implementations for scalar and fermionic monopoles, all spin cases were confronted by the theoretical predictions given in Sect. 3. Here again the MadGraph UFO models produced satisfactorily close cross-section values with that predicted by the theory, as shown in Table 2.

Comparison between the photon fusion and the Drell-Yan production mechanisms
Apart from the total cross sections, it is important to study the angular distributions of the generated monopoles. This is of great interest to the interpretation of the monopole searches in collider experiments, given that the geometrical acceptance and efficiency of the detectors is not uniform as a function of the solid angle around the interaction point. The kinematic distributions for the direct γ γ and qq scattering obtained with the UFO models were also compared against the calculated differential cross sections of Sect. 3 and showed good agreement with respect to the pseudorapidity, η, and the transverse momentum, p T , of the monopole. It is worth noting that the differential cross sections in Sect. 3 are plotted for a specific value of β 0.986, while in this section we consider a range of monopole velocities connected to the ratio of the selected monopole mass over the proton-proton collision energy, thus some differences in the PF-vs-DY comparison are expected.
Here the kinematic distributions are compared between the photon-fusion (γ γ ) and the Drell-Yan mechanisms. For this purpose, the β-dependent UFO monopole model was used in MadGraph. Monopole-antimonopole pair events have been generated for proton-proton collisions at √ s = 13 TeV, i.e. for the LHC Run-2 operating energy. The PDF was set to NNPDF23 [53] at LO for the Drell-Yan and LUXqed [54] for the photon-fusion mechanism. The latter choice is made due to the relatively small uncertainty in the photon distribution function in the proton provided by LUXqed [55]. The monopole magnetic charge is set to 1 g D , yet the kinematic spectra are insensitive to this parameter. The distributions are normalised to the same number of events, in order to facilitate the shape comparison.
The distributions of the monopole velocity, an important parameter for the detection of monopoles in experi-ments such as MoEDAL [23], are depicted in Fig. 22. The velocity β, which is calculated in the centre-of-mass frame of the colliding protons, largely depends on the PDF of the photon (quark/antiquark) in the proton for the photonfusion (Drell-Yan) process. For scalar monopoles, Fig. 22 (left) shows that slower-moving monopoles are expected for PF than DY, an observation favourable for the discovery potential of MoEDAL NTDs, the latter being sensitive to low-β monopoles. The comparison is reversed for fermionic monopoles, where PF yields faster monopoles than DY (Fig. 22 (centre)). Last, as deduced from Fig. 22 (right), the β distributions for PF and DY are very similar. It should be remarked that the plotted β expresses the monopole boost in the laboratory frame and, in the context of proton-proton collisions, it is distinct from the parameter β defined in (7) that enters into the monopole-photon coupling.
The kinetic energy 7 spectra are shown in Fig. 23. We choose to show distributions of the kinetic energy because it is relevant for the monopole energy loss in the detector material, hence important for the detection efficiency. The kineticenergy spectrum is slightly softer for PF than DY for scalar (left panel) and vector (right panel) monopoles, whereas it is significantly harder for fermions (central panel). This difference may be also due to the four-vertex diagram included in the bosonic monopole case. This observation is in agreement with the one made for β previously. We have also compared MadGraph predictions for kinetic-energy and p T distributions between with-and without-PDF cases, the latter also against analytical calculations (cf. Sect. 3) across different spins and production mechanisms. As expected, some features seen in the direct γ γ or qq production are attenuated in the pp production due to the sampling of different β values in the latter as opposed to the fixed value in the former.
As far as the pseudorapidity is concerned, its distributions are shown in Fig. 24. The spin-0 (left panel) and spin-1 (right  15, 18 (vector), respectively, as far as the production "centrality" is concerned. In the PF process, in particular, we observe that the depressions at η 0 for SM values of the κ parameters for both spin-1 /2 (see Fig. 9) and spin-1 (see Fig. 15) have been converted to flat tops when photon PDFs are also considered. This is normal taking into account the boost of the γ γ centreof-mass frame with respect to the pp (laboratory) frame and the event-by-event variation of the monopole velocity β that yields different event weight. The total cross sections for the various spin cases, assuming SM magnetic-moment values for spin 1 /2 and spin 1 and β-dependent coupling are drawn in Fig. 25 for photon fusion and Drell-Yan processes, as well as their sum. At tree level there is no interference between the PF and DY diagrams, so the sum of cross sections expresses the sum of the corresponding amplitudes. The PF mechanism is the dominant at the LHC energy of 13 TeV throughout the whole mass range of interest of 1 ÷ 6 TeV for the bosonic-monopole case. However if the monopole has spin 1 /2 the PF dominates for masses up to ∼ 5 TeV, while DY takes over for M 5 TeV. This results underlines the importance that the photon-fusion production mechanism has for LHC without, however, overlooking the DY process.  (56) cases. In this section, we focus on this aspect of the photon-fusion production mechanism, since it dominates at LHC energies. We first put to test this theoretical claim utilising the MadGraph implementation and later we discuss the kinematic distributions and comment on experimental aspects of a potential perturbatively-consistent search in colliders to follow in this context.

Spin-1 /2 case
For a spin-1 /2 monopole, the dimensionless parameterκ = κ M, with M the mass of the monopole, is varied from zero (the SM scenario) to 10,000 for γ γ collisions at √ s γ γ = 13 TeV. The cross-section of the photon fusion process for κ = 0 is going to zero very fast as β → 0, as can be seen in the third column of Table 3. However for non-zeroκ, the cross-section values remain finite even if β goes to zero, as expected, as becomes evident from the last row of Table 3. The same conclusion is drawn from Fig. 26 (left), where the cross sections are plotted for pp collisions at √ s = 13 TeV, i.e. with PDF. For masses M 6 TeV, the monopole production, although still rare, remains at detectable limits for the LHC experiments.
The central and right-hand-side plots of Fig. 26 depict a comparison of the p T and η distributions, respectively, between the SM valueκ = 0 and much higher values up toκ = 10 4 . The SM-like case is characterised by a distinguishably "softer" p T spectrum and a less central angular distribution than the large-κ case. 8 The latter case, on the   other hand, seems to converge to a common shape for the kinematic variables asκ increases to very large values. This is not the case forκ values distinct, yet near, the SM value, where the angular distributions are not considerably different, as shown in Fig. 9. The common kinematics among largẽ κ values would greatly facilitate an experimental analysis targeting perturbatively reliable results. We note here that the DY process, which dominates the cross section for heavy monopoles for the SM magnetic-moment value (see Fig. 25), vanishes asκ acquires large values as shown in Sect. 3.5, rendering the study of photon-fusion process sufficient at the perturbative-coupling limit.

Spin-1 case
Repeating the same procedure this time for spin-1 monopoles, we vary the dimensionless parameter κ from unity (the SM scenario) to 10,000 and we list the γ γ → MM cross sections in Table 4 for the photon-fusion process. Similar to the spin-1 /2 monopole case, the cross section for κ = 1, i.e. the SM scenario, is going to zero very fast as β → 0. However, for κ > 1, the cross section becomes finite even if β goes to 0, as seen in the last row of Table 4. This observation remains valid when the cross section for pp collisions, instead of γ γ scattering, is considered. Indeed, Fig. 27 shows that for large values of κ and M √ s/2, which is equivalent to β 1, the cross section, although very small, remains finite.
In Fig. 26, a comparison of the p T (centre) and η (right) distributions, between the SM value κ = 1 and much higher  values up to κ = 10 4 is given for spin 1. As for spin 1 /2, the large-κ curves converge to a single shape independent of the actual κ value, however the SM-case here also yield similar distributions as the for κ 1. 9 Therefore the ηdistribution features shown in Fig. 15 without PDF are completely smoothed out by folding with the photon distribution function in the proton. As discussed in Sect. 5.2.1 for fermionic monopoles, such an experimental analysis can be concentrated on the κ-dependence of the total cross section and the acceptance for very slow monopoles to provide perturbatively valid mass limits in case of non-observation of a monopole signal, since the kinematic distributions are κinvariant in pp collisions. The MoEDAL experiment [23,24], in particular, being sensitive to slow monopoles can make 9 As in the spin-1 /2-monopole case, we also observe here that, in the absence of PDF, in the non-unitary cases κ = 1, there is a fast-increasing distribution of events at high p T -values up to a cutoff of p T = √ s γ γ /2. the best out of this new approach in the interpretation of monopole-search results.

Conclusions
The work described in this article consists of two parts. In the first part, we dealt with the computation of differential and total cross sections for pair production of monopoles of spin S = 0, 1 /2, 1, through either photon-fusion or Drell-Yan processes. We have employed duality arguments to justify an effective monopole-velocity-dependent magnetic charge in monopole-matter scattering processes. Based on this, we conjecture that such β-dependent magnetic charges might also characterise monopole production. A magnetic-moment term proportional to a new phenomenological parameter κ is added to the effective Lagrangians describing the interactions of these monopoles with photons for spins 1 /2 and 1. The lack of unitarity and/or renormalisability is restored when the monopole effective theory adopts a SM form, that is when the bare magneticmoment parameter takes on the values κ = 0 for spin-1 /2 monopoles, and κ = 1 for spin-1 monopoles. However we remark that the lack of unitarity and renormalisability is not necessarily an issue, from an effective-field-theory point of view. Indeed, given that the microscopic high-energy (ultraviolet) completion of the monopole models considered above is unknown, one might not exclude the possibility of restoration of unitarity in extended theoretical frameworks, where new degrees of freedom at a high-energy scale might play a role. In this sense, we consider the spin-1 monopole as a potentially viable phenomenological case worthy of further exploration.
The motivation behind the magnetic-moment introduction is to enrich the monopole phenomenology with the (undefined) correction terms to the monopole magnetic moment to be treated as free parameters potentially departing from the ones prescribed for the electron or W ± bosons in the SM. Lacking a fundamental microscopic theory of magnetic poles, such an addition appears reasonable. This creates a dependence of the scattering amplitudes of processes on this parameter, which is passed on to the total cross sections and, in some cases, to kinematic distributions. Therefore the parameter κ is proposed as a tool for monopole searches which can be tuned to explore different models.
Moreover, even more intriguing is the possibility to use the parameter κ in conjunction with the monopole velocity β to achieve a perturbative treatment of the monopolephoton coupling. Indeed, in general the large value of the magnetic charge prevents any perturbative treatment of the monopole interactions limiting us to a necessary truncation of the Feynman-like diagrams at the tree level. By limiting the discussion to very slow (β 1) monopoles, the perturbativity is guaranteed, however, at the expense of a vanishing cross section. Nonetheless it turns out that the photon-fusion cross section remains finite and the coupling is perturbative at the formal limits κ → ∞ and β → 0. This ascertainment opens up the possibility to interpret the cross-section bounds set in collider experiments, such as MoEDAL, in a proper way, thus yielding sensible monopole-mass limits.
In the second part of this article, a complete implementation in MadGraph of the monopole production is performed both for the photon-fusion and the Drell-Yan processes, also including the magnetic-moment terms. The UFO models were successfully validated by comparing cross-section values obtained by the theoretical calculations and the Mad-Graph UFO models. Kinematic distributions, relevant for experimental analyses, were compared between the photonfusion and the Drell-Yan production mechanism of spins 0, 1 /2 and 1 monopoles. This tool will allow to probe for the first time the -dominant at LHC energies -photonfusion monopole production. Furthermore, the experimental aspects of a perturbatively valid monopole search for large values of the magnetic-moment parameters and slow-moving monopoles have also been outlined, based on these kinematic distributions.
The sum over massive (of mass M) vector meson polarisation states reads: The spinor sum rules and properties of Dirac's γ matrices are: Tr[odd number of γ matrices] = 0, from which we obtain for the trace over spinor polarisations: 1. The angular differential cross section distribution dσ d is defined as usual where |M| 2 is the squared matrix amplitude, averaged over initial spin or polarisation states (depending on the process) and summed over final ones, with the matrix element M denoting the sum of the matrix amplitudes of all the relevant processes contributing to the interaction. This expression is evaluated in terms of the Mandelstam variables defined as where m i are the masses of initial and final state particles i for which q 2 = m 2 q and p 2 = m 2 p as defined in Fig. 28. 2. The kinematics of the scattering of two initial-state particles is shown in Fig. 28 in the centre-of-mass frame. Fig. 28 Definition of coordinates pertaining to the particle-antiparticle pair production from the interaction of two initial-state particles in the centre-of-mass frame. The process has cylindrical symmetry in the longitudinal direction (q i ), perpendicular to which the azimuth φ is defined. The scattering angle θ and pseudorapidity η of the (final-state) M and M particles is also defined These states with momenta q 1 , q 2 interact with the same energy and with three momentum vectors that have opposite directions and equal magnitudes. In this frame, the energies E q 1 = E q 2 and E p 1 = E p 2 . Finally, it is worth noting that the s-channel exchange energy s = ( p 1 + p 2 ) 2 = (q 1 + q 2 ) 2 demands that E q = E p . Using the geometry of the system, |M| 2 is reduced to its final form as a function of θ , defined as the angle between the axes of propagation of the initial state particles and the scattered monopoles (see Fig. 28). The Mandelstam variables in this frame take the form t = m 2 q +m 2 p − s 2 (1−β p β q cos θ) and u = m 2 q + m 2 p − s 2 (1 + β p β q cos θ). 3. At this point, |M| 2 is inserted into the expression for the angular cross section distribution in (A7). The behaviour of the cross section can be studied as a function of kinematic variables other than θ by a change of variable. Specifically, experimental monopole searches are usually interested in expressing the pertinent cross sections in terms of the pseudorapidity η, defined as The latter represents an angular coordinate relative to the beam axis as identified in Fig. 28, or alternatively the significance of the longitudinal boost relative to the total momentum. It is also the high energy limit of the rapidity y, which is the relativistic, i.e. Lorentz-invariant, realisation of velocity. It becomes the new kinematic variable through the coordinate transformation 1 cosh 2 η dη = 1 2π d .
In the centre-of-mass frame, s γ γ = 4E q 1 E q 2 and the photon three-momenta q 1 = q 2 are defined as | q| = E q for any massless photon. It follows from this that the s-channel exchange energy s γ γ = ( p 1 + p 2 ) 2 = (q 1 + q 2 ) 2 requires that E q = E p . As discussed in Sect. 2, this squared matrix amplitude is used to evaluated the kinematic distributions in terms of the angle θ and the pseudorapidity η. The total cross section in the centre-of-mass frame is obtained by integrating (B7) over the solid angle d = dφ dcos θ : (B8) Scalar monopole pairs produced by the DY digram shown in Fig. 5 have kinematic distributions the analytic forms of which are derived analogously to the PF case. The total matrix amplitude, derived from the Feynman rules in Fig. 30, is given by where u αi , v β j , u αi , v β j are the quark fermionic spinors with δ i j the colour index factor and g μν is the Minkowski spacetime metric tensor. The quarks have a mass m each and are characterised by momentum 4-vectors q 1μ and q 2μ , where q 2 1,2 = m 2 . Similarly, the scalar monopoles have mass M and momentum 4-vectors p 1μ and p 2μ , where, on shell, p 2 1,2 = M 2 . The quark masses are negligible compared to a monopole mass. The quantity e represents the positron charge and Q is the charge fraction relevant for quarks, Q = 1 3 , 2 3 . The squared matrix amplitude, averaged over quark spins and colours, |M 2 | DY is then computed: The 2!, 1 4 1 3 factors account for the symmetry factor of the final states, and the averaging over spins and colour states. The factor of 3 includes the contribution from 3 quark flavours. Finally, the sum over quark charges adds a factor of 5 9 to the result. The rest of the expression was found using the spin sum rules and properties of γ -matrices listed in "Appendix A" (see Eqs. (A5), (A6)).
Using the relationship between the Mandelstam variables (A8) and the scattering angle between the trajectories of the incoming quarks and outgoing scalar bosons as depicted in Fig. 28, with β q = |q|/E q and β p = |p|/E p , (B11) simplifies to: which is used to compute the differential cross section.
The quark masses are assumed negligible at this point, hence β q = 1 and the subscript is dropped on β p . The total cross section is obtained by integrating over the solid angle as in the photon fusion case.

Spin-1 /2 monopole cross section
Starting with pair production by PF, the kinematic distributions for spin-1 /2 monopoles are derived. The total matrix amplitude for the process is the sum of a t-and a u-channel process, shown in Fig. 8, M P H = M t + M u . Each amplitude is derived from the Feynman rule (B17), shown in Fig. 31. The fermion in the latter may represent a spin-1 /2 monopole with magnetic moment κ = 0, upon dualising the theory by replacing the electric charge by a (in general β-dependent) magnetic charge g(β). The amplitude in the latter case is given in Eq. (B17) where the influence of the added magnetic moment term is visible. The term [γ μ , γ ν ] is a commutator of γ matrices and q μ is the photon fourmomentum.
(B17) where the monopole boost is β = | p|/E p and the photon momentum is | q| = E q , recognising also that E q = E p in the centre-of-mass frame. It is expressed in terms of the scattering angle θ and the pseudorapidity η as defined in Fig. 28 and explained in Sect. 2.
Moving on now to monopole pair production by DY, the matrix amplitude is derived from the Feynman rule (B21) shown in Fig. 32.