Slepton Production at e+e- Colliders in the Complex MSSM: A Full One-Loop Analysis

For the search for scalar leptons in the Minimal Supersymmetric Standard Model (MSSM) as well as for future precision analyses of these particles an accurate knowledge of their production and decay properties is mandatory. We evaluate the cross sections for the slepton production at e+e- colliders in the MSSM with complex parameters (cMSSM). The evaluation is based on a full one-loop calculation of the production mechanisms e+e- ->slepton_s slepton_s' including soft and hard photon radiation. The dependence of the slepton production cross sections on the relevant cMSSM parameters is analyzed numerically. We find sizable contributions to many production cross sections. They amount to roughly 15% of the tree-level results but can go up to 40% or higher in extreme cases. Also the complex phase dependence of the one-loop corrections for charged slepton production was found non-negligible. The full one-loop contributions are thus crucial for physics analyses at a future linear e+e- collider such as the ILC or CLIC.


Introduction
One of the important tasks at the LHC is to search for physics beyond the Standard Model (SM), where the Minimal Supersymmetric Standard Model (MSSM) [1][2][3][4] is one of the leading candidates. Two related important tasks are the investigation of the mechanism of electroweak symmetry breaking, including the identification of the underlying physics of the Higgs boson discovered at ∼ 125 GeV [5,6], as well as the production and measurement of the properties of Cold Dark Matter (CDM). Here the MSSM offers a natural candidate for CDM, the Lightest Supersymmetric Particle (LSP), the lightest neutralino,χ 0 1 [7,8] (see below). These three (related) tasks will be the top priority in the future program of particle physics.
Supersymmetry (SUSY) predicts two scalar partners for all SM fermions as well as fermionic partners to all SM bosons. Contrary to the case of the SM, in the MSSM two Higgs doublets are required. This results in five physical Higgs bosons instead of the single Higgs boson in the SM. These are the light and heavy CP-even Higgs bosons, h and H, the CP-odd Higgs boson, A, and the charged Higgs bosons, H ± . In the MSSM with complex parameters (cMSSM) the three neutral Higgs bosons mix [9][10][11][12][13], giving rise to the CP-mixed states h 1 , h 2 , h 3 . The neutral SUSY partners of the (neutral) Higgs and electroweak gauge bosons are the four neutralinos,χ 0 1,2,3,4 . The corresponding charged SUSY partners are the charginos,χ ± 1,2 . If SUSY is realized in nature and the scalar quarks and/or the gluino are in the kinematic reach of the (HL-)LHC, it is expected that these strongly interacting particles are eventually produced and studied. On the other hand, SUSY particles that interact only via the electroweak force, e.g., the scalar leptons, have a much smaller production cross section at the LHC. Correspondingly, the LHC discovery potential as well as the current experimental bounds are substantially weaker [14,15].
At a (future) e + e − collider sleptons, depending on their masses and the available centerof-mass energy, could be produced and analyzed in detail. Corresponding studies can be found for the ILC in Refs. [16][17][18][19][20][21] and for CLIC in Refs. [21][22][23]. (Results on the combination of LHC and ILC results can be found in Refs. [24][25][26].) Such precision studies will be crucial to determine their nature and the underlying (SUSY) parameters.
In order to yield a sufficient accuracy, one-loop corrections to the various slepton production and decay modes have to be considered. Full one-loop calculations in the cMSSM to (heavy) scalar tau decays was evaluated in Ref. [27], where the calculation can easily be taken over to other slepton decays. Sleptons can also be produced in SUSY cascade decays. Full one-loop evaluations in the cMSSM exist for the corresponding decays of Higgs bosons [28] as well as from charginos and neutralinos [29][30][31][32]. In this paper we take the next step and concentrate on the slepton production at e + e − colliders, i.e. we calculate σ(e + e − →ẽ ± gsẽ ∓ gs ) s, s = 1, 2 , σ(e + e − →ν gν * g ) g = 1, 2, 3 , In this paper we present for the first time a full and consistent one-loop calculation in the cMSSM for scalar lepton production at e + e − colliders. We take into account soft, hard and collinear QED radiation and the treatment of collinear divergences. Again, here it is crucial to stress that the same renormalization scheme as for the decay of sleptons [27] (and for slepton production from Higgs boson decays [28] as well as from chargino and neutralino decays [29][30][31][32]) has been used. Consequently, the predictions for the production and decay can be used together in a consistent manner (e.g., in a global phenomenological analysis of the slepton sector at the one-loop level). We analyze all processes w.r.t. the most relevant parameters, including the relevant complex phases. In this way we go substantially beyond the existing analyses (see above). In Sect. 2 we briefly review the renormalization of the relevant sectors of the cMSSM and give details as regards the calculation. In Sect. 3 various comparisons with results from other groups are given. The numerical results for the production channels (1) and (2) are presented in Sect. 4. The conclusions can be found in Sect. 5.
They will be further explained in the text below.

Calculation of diagrams
In this section we give some details regarding the renormalization procedure and the calculation of the tree-level and higher-order corrections to the production of sleptons in e + e − collisions. The diagrams and corresponding amplitudes have been obtained with FeynArts (version 3.9) [45][46][47], using our MSSM model file (including the MSSM counterterms) of Ref. [48]. The further evaluation has been performed with FormCalc (version 9.5) and LoopTools (version 2.14) [49,50].

The complex MSSM
The cross sections (1) and (2) are calculated at the one-loop level, including soft, hard and collinear QED radiation; see the next section. This requires the simultaneous renormalization of the gauge-boson sector, the lepton sector as well as the slepton sector of the cMSSM. We give a few relevant details as regards these sectors and their renormalization. More details and the application to Higgs-boson and SUSY particle decays can be found in Refs. [27-32, 48, 51-55]. Similarly, the application to Higgs-boson and chargino/neutralino production cross sections at e + e − colliders are given in Refs. [56][57][58].
The renormalization of the fermion and gauge-boson sector follows strictly Ref. [48] and the references therein (see especially Ref. [59]). This defines in particular the counterterm δt β , as well as the counterterms for the Z boson mass, δM 2 Z , and for the sine of the weak mixing angle, δs w (with s w = 1 − c 2 w = 1 − M 2 W /M 2 Z , where M W and M Z denote the W and Z boson masses, respectively). For the fermion sector we use the default values as given in Ref. [48].
The "m b , A b DR" scheme of Refs. [52][53][54], extended to all generations is herein after referred to as mixed scheme, of course not to be confused with the mixed scheme of Ref. [48].
The schemes affectingẽ g andd g are chosen with the variable $SfScheme[t, g]: The non-diagonal entries of the up-type mass matrix are determined by [52,53] dMSfsq1 [1,2,3, g] ≡ δY ug = 1 2 Re Σũ g (m 2 u g1 ) + Σũ g (m 2 u g2 ) 12 , (4e) For clarity of notation we furthermore define the auxiliary constants For the bottom quark we choose: δm b = δm DR b . In the mixed scheme the down-type off-diagonal mass counterterms are related as Uẽ g 11 Uẽ g * 21 δm 2 e g1 − δm 2 e g2 + (6a) The trilinear couplings A ftg ≡ A ft gg are renormalized by , where the subscripted [div] means to take the divergent part in the mixed scheme only, to effect DR renormalization of A eg and A dg [52,53]. The squark and slepton Z-factors are derived in the OS scheme and can be found in Section 3.6.1 and 3.6.2 of Ref. [48].
As now all the sfermion masses are renormalized as on-shell an explicit restoration of the SU (2) L relation is needed. Requiring the SU (2) L relation to be valid at the one-loop level induces the following shifts in the soft SUSY-breaking parameters: with Now M 2 L (ẽ g ) and M 2 Q (d g ) are used in the scalar electron-and down-type mass matrix instead of the parameters M 2 L,Q in the sfermion mass matrix when calculating the values of mẽ gs and md gs . However, with this procedure, both (s = 1, 2) sfermion masses are shifted, which contradicts our choice of independent parameters. To keep this choice, also the right-handed soft SUSY-breaking mass parameters MẼ ,D receive a shift: with our choice of mass ordering, mf 1 < mf 2 , we have Taking into account the shift Eq. (11) in MẼ and Eq. (12) in MD, up to one-loop order 2 , the new resulting mass parameters mẽ gs and md gs are the same as the on-shell masses: where δm dep. egs and δm dep. dgs are the dependent mass counterterms. A slightly different slepton sector renormalization is also described in detail in Ref. [48] and references therein.

Contributing diagrams
Sample diagrams for the process e + e − →ẽ ± gsẽ ∓ gs and e + e − →ν gν * g are shown in Fig. 1. Not shown are the diagrams for real (hard and soft) photon radiation. They are obtained from the corresponding tree-level diagrams by attaching a photon to the (incoming/outgoing) electron or slepton. The internal particles in the generically depicted diagrams in Fig. 1 are labeled as follows: F can be a SM fermion f , charginoχ ± c or neutralinoχ 0 n ; S can be a sfermionf s or a Higgs (Goldstone) boson h 0 , H 0 , A 0 , H ± (G, G ± ); U denotes the ghosts u V ; V can be a photon γ or a massive SM gauge boson, Z or W ± . We have neglected all electron-Higgs couplings and terms proportional to the electron mass whenever this is safe, i.e. except when the electron mass appears in negative powers or in loop integrals. We have verified numerically that these contributions are indeed totally negligible. For internally appearing Higgs bosons no higher-order corrections to their masses or couplings are taken into account; these corrections would correspond to effects beyond one-loop order. 3 Moreover, in general, in Fig. 1 we have omitted diagrams with self-energy type corrections of external (on-shell) particles. While the contributions from the real parts of the loop functions are taken into account via the renormalization constants defined by OS renormalization conditions, the contributions coming from the imaginary part of the loop functions can result in an additional (real) correction if multiplied by complex parameters. In the analytical and numerical evaluation, these diagrams have been taken into account via the prescription described in Ref. [48].
Within our one-loop calculation we neglect finite width effects that can help to cure threshold singularities. Consequently, in the close vicinity of those thresholds our calculation does not give a reliable result. Switching to a complex mass scheme [60] would be another possibility to cure this problem, but its application is beyond the scope of our paper.
The tree-level formulas σ tree (e + e − →ẽ ± gsẽ ∓ gs ) and σ tree (e + e − →ν gν * g ) are given in Refs. [34,36] and Ref. [40], respectively. Concerning our evaluation of σ(e + e − →ẽ ± gsẽ ∓ gs ) we define: if not indicated otherwise. Differences between the two charge conjugated processes can appear at the loop level when complex parameters are taken into account, as will be discussed in Sect. 4.2.

Ultraviolet, infrared and collinear divergences
As regularization scheme for the UV divergences we have used constrained differential renormalization [61], which has been shown to be equivalent to dimensional reduction [62,63] at the one-loop level [49,50]. Thus the employed regularization scheme preserves SUSY [64,65] and guarantees that the SUSY relations are kept intact, e.g., that the gauge couplings of the SM vertices and the Yukawa couplings of the corresponding SUSY vertices also coincide to one-loop order in the SUSY limit. Therefore no additional shifts, which might occur when using a different regularization scheme, arise. All UV divergences cancel in the final result. Soft photon emission implies numerical problems in the phase space integration of radiative processes. The phase space integral diverges in the soft energy region where the photon  Figure 1: Generic tree, self-energy, vertex, box, and counterterm diagrams for the process e + e − →l gslgs (l gs = {ẽ gs ,ν g }; g = 1, 2, 3; s, s = 1, 2). The additional diagrams, which occur only in the case of first generation slepton production, are denoted withl 1s . F can be a SM fermion, chargino or neutralino; S can be a sfermion or a Higgs/Goldstone boson; V can be a γ, Z or W ± . It should be noted that electron-Higgs couplings are neglected. 3.359 ± 0.015 10 −5 3.361 ± 0.020 10 −6 3.342 ± 0.023 10 −7 3.299 ± 0.028 10 −8 3.231 ± 0.028 momentum becomes very small, leading to infrared (IR) singularities. Therefore the IR divergences from diagrams with an internal photon have to cancel with the ones from the corresponding real soft radiation. We have included the soft photon contribution via the code already implemented in FormCalc following the description given in Ref. [66]. The IR divergences arising from the diagrams involving a photon are regularized by introducing a photon mass parameter, λ. All IR divergences, i.e. all divergences in the limit λ → 0, cancel once virtual and real diagrams for one process are added. We have numerically checked that our results do not depend on λ or on ∆E = δ s E = δ s √ s/2 defining the energy cut that separates the soft from the hard radiation. As one can see from the example in Fig. 2 this holds for several orders of magnitude. Our numerical results below have been obtained for fixed δ s = 10 −3 . Numerical problems in the phase space integration of the radiative process arise also through collinear photon emission. Mass singularities emerge as a consequence of the collinear photon emission off massless particles. But already very light particles (such as electrons) can produce numerical instabilities. For the treatment of collinear singularities in the photon radiation off initial state electrons and positrons we used the phase space slicing method [67][68][69][70], which is not (yet) implemented in FormCalc and therefore we have developed and implemented the code necessary for the evaluation of collinear contributions; see also Refs. [56,57].
In the phase space slicing method, the phase space is divided into regions where the integrand is finite (numerically stable) and regions where it is divergent (or numerically unstable). In the stable regions the integration is performed numerically, whereas in the unstable regions it is carried out (semi-) analytically using approximations for the collinear photon emission.  The collinear part is constrained by the angular cut-off parameter ∆θ, imposed on the angle between the photon and the (in our case initial state) electron/positron.
The differential cross section for the collinear photon radiation off the initial state e + e − pair corresponds to a convolution dσ coll (s) = α π with P ee (z) = (1 + z 2 )/(1 − z) denoting the splitting function of a photon from the initial e + e − pair. The electron momentum is reduced (because of the radiated photon) by the fraction z such that the center-of-mass frame of the hard process receives a boost. The integration over all possible factors z is constrained by the soft cut-off δ s = ∆E/E, to prevent over-counting in the soft energy region.
We have numerically checked that our results do not depend on the angular cut-off parameter ∆θ over several orders of magnitude; see the example in Fig. 3. Our numerical results below have been obtained for fixed ∆θ/rad = 10 −2 .
The one-loop corrections of the differential cross section are decomposed into the virtual, soft, hard, and collinear parts as follows: The hard and collinear parts have been calculated via Monte Carlo integration algorithms of the CUBA library [71,72] as implemented in FormCalc [49,50].

Comparisons
In this section we present the comparisons with results from other groups in the literature for slepton production in e + e − collisions. These comparisons were restricted to the MSSM with real parameters, with one exception for tree-level tau slepton pair production. The level of agreement of such comparisons (at one-loop order) depends on the correct transformation of the input parameters from our renormalization scheme into the schemes used in the respective literature, as well as on the differences in the employed renormalization schemes as such. In view of the non-trivial conversions and the large number of comparisons such transformations and/or change of our renormalization prescription are beyond the scope of our paper. In the following we list all relevant papers in the literature and explain either our comparison, or why no (meaningful) comparison could be performed.
• In Ref. [33] the production of slepton and squark pairs in e + e − annihilation and Z decay have been calculated in the rMSSM at tree level (including arbitrarily polarized beams). Unfortunately, in Ref. [33] not sufficient information as regards their input parameters where given, rendering a comparison impossible.
• Selectron pair production at e + e − colliders in SUGRA models were presented in Ref. [34] at tree level. We omitted a comparison with Ref. [34], since implementing the SUGRA spectrum is beyond the scope of our paper.
• Tree-level tau slepton pair production at muon colliders in the cMSSM were analyzed in Ref. [35] including CP violation. The center-of-mass energy was assumed to be around the resonances of the heavy neutral Higgs bosons, at √ s 500 GeV. We used their input parameters as far as possible, but we differ from their results, especially in the case of complex input parameters. This can most likely be attributed to the differences in the Higgs-boson mass calculations employed in Ref. [35] and the current version 2.13.0 of FeynHiggs [59,[73][74][75][76][77][78][79], which we use. The cross section (close to the heavy Higgs boson thresholds) depends strongly on tiny mass differences of the two heavy neutral Higgs bosons, which deviate clearly between the two employed calculations, rendering this comparison not significant.
• In Refs. [36,37] pair production of smuons and selectrons near threshold in e + e − and e − e − collisions have been computed in the rMSSM, including Coulomb rescattering effects. 4 Because these near production threshold effects are beyond the scope of our paper we have omitted a comparison with Refs. [36,37].
• Ref. [38] deals with selectron pair production at e − e − and e + e − colliders with polarized beams. The calculations have been computed in the rMSSM near the production threshold. Again, these effects are beyond the scope of our paper and we have omitted a comparison with Ref. [38].
• Also Ref. [39] is (mainly) based on Ref. [37]. The general production of scalar leptons at linear e + e − colliders has been computed in Ref. [39], including multi-photon initial  state radiation and polarized beams. The authors used (older versions of) FeynArts, FeynCalc [80] and LoopTools for their calculations. As input parameters they used the mSUGRA parameter point SPS1a [81], translated from the DR to on-shell values; see the appendix of Ref. [39]. We also used this parameter point (as far as possible) and reproduced successfully their tree-level results in their Figs. 2 and 3a (see Ref. [39]) in the upper row of our Fig. 4. Our (relative) one-loop results are in qualitative agreement with the ones in their Figs. 17 and 18a of Ref. [39]; see the lower row in our Fig. 4. The quantitative numerical differences can be explained with the different renormalization schemes, slightly different input parameters, and the different treatment of the photon bremsstrahlung, where they have included multi-photon emission while we kept our calculation at O(α). It should also be kept in mind that the relative one-loop corrections are sensitive to every kind of difference.
• Ref. [40] is an addendum to Ref. [39] dealing with e + e − →ν gν * g . As input parameters e + e − →ν gν * g they used the SUSY parameter point SPS1a translated from the DR to on-shell values; see Ref. [81]. We reproduced successfully e + e − →ν µν * µ and e + e − →ν τν * τ of Ref. [40] in our Fig. 5. The (quantitative) difference in the relative loop corrections can be explained with the different renormalization schemes, slightly different input parameters, and the different treatment of the photon radiation, where they have included multi-photon emission while we kept our calculation at O(α). The process e + e − →ν eν * e , on the other hand, while in rather good qualitative agreement, differ quantitative significantly already at the tree level. Unfortunately, we were not able to trace back the source of the difference. However, since (in our automated approach) we agree with other tree-level calculations, we are confident that our results are correct.
• We performed a comparison with Ref. [41] using their input parameters (as far as possible). They have calculated (third generation) scalar fermion production in the rMSSM within an on-shell scheme close to ours at the "full" one-loop level (but without explicit QED radiation). They also used (older versions of) FeynTools for their calculations. We found very good agreement with their Fig. 13, as can be seen in our Fig. 6. The tiny differences can easily be explained with the slightly different SM input parameters and the slightly different renormalization scheme.
• In Refs. [42,43] the "complete" one-loop corrections to e + e − →f sf * s (third generation) in the rMSSM were analyzed numerically including photon corrections. 5 However, their numerical results have been presented taking into account only weak corrections. They used (an older version of) LoopTools for their calculations and their own on-shell renormalization procedure together with a α(M Z )| MS = 1/127.934 scheme. We used their "Scenario 1-gaugino" (as far as possible) for our comparison. We are in rather good agreement with their Figs. 4c, 4f, 7a, and 7d which can be seen in our Fig. 7. The minor differences can be explained (as usual) with the slightly different input parameters and the different renormalization scheme. 6 • Finally, in Refs. [43,44] full O(α) corrections to e + e − →f sf * s have been calculated in the MSSM with real parameters. Polarized electrons and multi-photon bremsstrahlung were included in the phenomenological analysis. As input parameters they used the mSUGRA parameter point SPS1a [81], translated from the DR to on-shell values. We also used this parameter point (as far as possible) and reproduced parts of their results in their Figs. 7 and 8 (see Ref. [44]) in our Fig. 8. We are in good agreement for the processes e + e − →τ + 2τ − 2 and e + e − →τ ± 1τ ∓ 2 , while we disagree in e + e − →τ + 1τ − 1 and e + e − →ν τν * τ already at tree-level. Unfortunately, we were not able to track down the 5 It should be noted that Ref. [43] is the "source" of Ref. [42] and Ref. [44]. 6 It should be noted that the sum, but not the individual contributions, of vertex (vert) and propagator (prop) contributions (in our upper right plot of Fig. 7) are (nearly) the same as the corresponding sum in Fig. 4f of Ref. [42]. This is because of the different renormalization of the (charged) slepton sector. source of the differences. However, since we agree (in our automated approach) with other tree-level calculations, we are confident that our results are correct.
To conclude, in most cases where a meaningful comparison could be performed, we found good agreement with the literature where expected, and the encountered differences can be traced back to different renormalization schemes, corresponding mismatches in the input parameters and small differences in the SM parameters. Nevertheless, in some cases we disagree already significantly at tree-level but we were not able to track down the source of these differences. This does not disprove the reliability of our calculation because our computational method/code has already been successfully tested and compared with quite a few other programs; see Refs. [27-29, 31, 48, 51, 52, 54, 55]. After comparing with the existing literature we would like to stress again that here we present for the first time a full one-loop calculation of σ(e + e − →ẽ ± gsẽ ∓ gs ) and σ(e + e − →ν gν * g ) in the cMSSM, using the scheme that was employed successfully already for the full one-loop decays of the (produced) sleptons. The various calculations can readily be used together for the full production and decay chain.

Numerical analysis
In this section we present our numerical analysis of slepton production at e + e − colliders in the cMSSM. In the figures below we show the cross sections at the tree level ("tree") and at the full one-loop level ("full"), which is the cross section including all one-loop corrections as described in Sect. 2. The CCN [1] renormalization scheme (i.e. OS conditions for the two charginos and the lightest neutralino) has been used for most evaluations. In cases where the CCN[1] scheme is divergent (i.e. M 2 = µ) and/or unreliable, also some CNN[c, n, n ] schemes (OS conditions for one chargino and two neutralinos) have been used, as indicated below; see Ref. [48] for further informations to these renormalization schemes. It should be noted that within the CCN[1] scheme at µ = M 2 a divergence already in the tree-level result can be induced by δµ through the shifted scalar lepton masses 7 : δµ in Eq. (5a) enters via Eq. (6a) into Eq. (10) from which the slepton shifts are calculated. The divergence is suppressed with m eg (see Eq. (5a)). In several analyses in Sect. 4 in order to overcome the problem with a divergent tree-level result, we switched to the tree-level result of the CNN[1,2,3] scheme, which is free of such a divergence. When several schemes are shown in one plot, a full comparison would require the transition of the relevant input parameters (which are varied). However, we do not intend to perform an analysis of the advantages and disadvantages of the various renormalization schemes. We want to demonstrate, however, that it is always possible to choose a "good" renormalization scheme, i.e. a scheme that leads to stable and not excessively large higher-order corrections. Consequently, the above mentioned parameter conversion is not (yet) included in our calculation. We first define the numerical scenario for the cross section evaluation. Then we start the numerical analysis with the cross sections of e + e − →ẽ ± gsẽ ∓ gs (g = 1, 2, 3; s, s = 1, 2) in Sect. 4.2, evaluated as a function of √ s, MẼ, µ, |M 1 | and/or M 2 , ϕ Ae g or ϕ M 1 , the phase of M 1 . In some cases also the t β dependence is shown. Then we turn to the processes 7 We use the same (shifted) sfermion masses for the tree-level and the full one-loop corrected cross section.

Parameter settings
The renormalization scale µ R has been set to the center-of-mass energy, √ s. The SM parameters are chosen as follows; see also [82]: • Fermion masses (on-shell masses, if not indicated differently): According to Ref. [82], m s is an estimate of a so-called "current quark mass" in the MS scheme at the scale µ ≈ 2 GeV. m c ≡ m c (m c ) is the "running" mass in the MS scheme and is the DR bottom quark mass as calculated in Ref. [54]. m u and m d are effective parameters, calculated through the hadronic contributions to ∆α (5) • Gauge-boson masses: • Coupling constant: The SUSY parameters are chosen according to the scenario S, shown in Tab. 1. This scenario is viable for the various cMSSM slepton production modes, i.e. not picking specific parameters for each cross section. They are in particular in agreement with the relevant SUSY searches of ATLAS and CMS: Our electroweak spectrum is not covered by the latest ATLAS/CMS exclusion bounds. Two limits have to be distinguished. The limits not taking into account a possible intermediate slepton exclude a lightest neutralino only well below 300 GeV [83, 84], whereas in S we have mχ0 1 ≈ 323 GeV. Limits with intermediary sleptons often assume a chargino decay to lepton and sneutrino, while in our scenario mχ± 1 < mν e,µ,τ . Furthermore, the exclusion bounds given in the mχ0 1 -mχ0 2 mass plane (with mχ0 2 ≈ mχ± 1 assumed) above mχ0 2 ∼ 300 GeV do not cover a compressed spectrum [83, 85] forχ 0 1 ,χ 0 2 , and χ ± 1 . In particular our scenario S assumed masses of mχ0 1 ≈ 323 GeV and mχ0 2 ≈ 354 GeV, which are not excluded. It should be noted that higher-order corrected Higgs-boson masses do not enter our calculation. 8 However, we ensured that over larger parts of the parameter space the lightest Higgs-boson mass is around ∼ 125 ± 3 GeV to indicate the phenomenological validity of our scenarios. In our numerical evaluation we will show the variation with √ s (up to 3 TeV, shown in the upper left plot of the respective figures), MẼ (from 100 to 500 GeV, upper right plot), µ (starting at µ = 100 GeV up to µ = 1000 GeV, shown in the middle/lower left plots), |M 1 | or M 2 (from 100 to 1000 GeV, middle/lower right plots), and ϕ Ae g or ϕ M 1 (between 0 • and 360 • , lower right plots). The dependence of t β turned out to be rather small, therefore we show it only in a few cases, where it is of special interest. Concerning the complex parameters, some more comments are in order. Potentially complex parameters that enter the selectron and electron sneutrino production cross sections at tree level (via the t-channel exchange of a neutralino or chargino) are the soft SUSYbreaking parameters M 1 and M 2 as well as the Higgs mixing parameter µ. Also trilinear slepton couplings A eg enter the tree-level production cross sections. However, when performing an analysis involving complex parameters it should be noted that the results for physical observables are affected only by certain combinations of the complex phases of the parameters µ, the trilinear couplings A f and the gaugino mass parameters M 1,2,3 [94,95]. It is possible, for instance, to rotate the phase ϕ M 2 away. Experimental constraints on the (combinations of) complex phases arise, in particular, from their contributions to electric dipole moments of the electron and the neutron (see Refs. [96][97][98] and the references therein), of the deuteron [99] and of heavy quarks [100,101]. While SM contributions enter only at the three-loop level, due to its complex phases the MSSM can contribute already at one-loop order. Large phases in the first two generations of sfermions can only be accommodated if these generations are assumed to be very heavy [102,103] or large cancellations occur [104][105][106]; see, however, the discussion in Ref. [107]. A review can be found in Ref. [108]. Recently additional constraints at the two-loop level on some CP phases of SUSY models have been investigated in Ref. [109]. Accordingly (using the convention that ϕ M 2 = 0, as done in this paper), in particular, the phase ϕ µ is tightly constrained [110], and we set it to zero. On the other hand, the bounds on the phases of the third-generation trilinear couplings are much weaker. Consequently, the largest effects on the slepton production cross sections at the tree level are expected from the complex gaugino mass parameter M 1 , i.e. from ϕ M 1 . As mentioned above, the only other phase entering at the tree level, is ϕ Ae g . This motivates our choice of ϕ M 1 and ϕ Ae g as parameters to be varied.
Since now complex parameters can appear in the couplings, contributions from absorptive parts of self-energy type corrections on external legs can arise. The corresponding formulas for an inclusion of these absorptive contributions via finite wave function correction factors can be found in Refs. [48,54].
The numerical results shown in the next subsections are of course dependent on the choice of the SUSY parameters. Nevertheless, they give an idea of the relevance of the full one-loop corrections.
We start with the process e + e − →ẽ + 1ẽ − 1 shown in Fig. 9. Away from the production threshold, loop corrections of ∼ +6 % at √ s = 1000 GeV are found in scenario S (see Tab.  The cross sections are decreasing with increasing MẼ due to kinematics, and the full one-loop result has its maximum of ∼ 28 fb at MẼ = 100 GeV. Analogously the relative corrections are decreasing from ∼ +13 % at MẼ = 100 GeV to ∼ −30% at MẼ = 490 GeV. The tree crossing takes place at MẼ ≈ 375 GeV. For higher MẼ values the loop corrections are negative, where the relative size becomes large due to the (relative) smallness of the tree-level results, which goes to zero for MẼ ≈ 500 GeV due to kinematics. For the other parameter variations one can conclude that a cross section roughly twice as large can be possible for very low MẼ (which, however, are challenged by the current ATLAS/CMS exclusion bounds).
With increasing µ in S (middle left plot) we find a decrease of the loop corrected production cross section within CCN [1] (green dashed line). The relative loop corrections reach ∼ +18 % at µ = 100 GeV (at the border of the experimental exclusion bounds) and ∼ +7 % at µ = 600-1000 GeV. One can see the expected breakdown of the CCN[1] scheme for µ = M 2 , i.e. in our case at µ ≈ M 2 = 600 GeV (see also Refs. [31,32]). Therefore, in the middle left plot of Fig. 9   dash-dotted line) and CNN[2,1,3] (black dash-dotted line) schemes, which are smooth at µ = M 2 . Outside the region of µ ∼ M 2 the scheme CCN [1] is expected to be reliable, since each of the three OS conditions is strongly connected to one of the three input parameters, |M 1 |, M 2 and µ. Similarly, CNN[1,2,3] is expected to be reliable for µ smaller than M 2 , as in this case, again each of the three OS renormalization conditions is strongly connected to the three input parameters. This behavior can be observed in the plot: for µ < ∼ M 2 = 600 GeV CNN[1,2,3] is reliable, while for µ > ∼ M 2 = 600 GeV it becomes unreliable. While δµ in the CNN[2,1,3] scheme has a strong minimum at µ ≈ 421 GeV, dominating the loop corrections, it approximates CCN[1] very good for µ > ∼ 500 GeV. A rising deviation between the schemes can be observed for µ < 250 GeV, where the schemes CNN[1,2,3] and CNN[2,1,3] are nearly constant, i.e. independent of µ. Overall, it is possible to find for every value of µ a renormalization scheme that behaves stable and "flat" w.r.t. the tree-level cross section.
The dependence M 2 of the cross section in S is shown in the lower left plot. One can clearly see the expected breakdown of the CCN[1] scheme for M 2 = µ, i.e. in our case at M 2 ≈ µ = 350 GeV (see also Refs. [31,32]) and the smooth behavior of CNN[2,2,3] (red dashed line) around M 2 ∼ µ = 350 GeV. Within CCN[1] the cross section even turns out to be negative for M 2 < 154 GeV due to a maximum of δM 1 at M 2 ≈ 20 GeV, dominating the loop corrections. This renders the CCN[1] scheme to be unreliable for M 2 < ∼ µ. For M 2 > ∼ µ, on the other hand, the scheme CCN [1] is expected to be reliable, since each of the three OS conditions is strongly connected to one of the three input parameters, |M 1 |, M 2 and µ. The loop corrections at the level of ∼ +6 % are found to be nearly independent of M 2 for M 2 > ∼ µ within the CCN[1] scheme. And corrections of ∼ +7 % are found to be (nearly) independent of M 2 within the CNN[2,2,3] scheme. Now we turn to the complex phase dependence. We find that the phase dependence ϕ M 1 of the cross section in S is large for CCN [1] (lower right plot, green dashed line). Loop corrections at the level of ∼ +6 % at ϕ M 1 = 0 • , 360 • and ∼ −48 % at ϕ M 1 = 180 • are found. It should be noted here, that there is no divergency or threshold at ϕ M 1 = 180 • , the cross section is smooth/finite; see the inlay in the lower right plot. This large structure is caused by a strong (local) minimum of the renormalization constant δM 1 at ϕ M 1 = 180 • , dominating the loop corrections. Using another renormalization scheme (e.g. CNN [2,2,3]) the dip disappears, as can be seen in the plot. The one-loop corrections then reach ∼ +7 % and are nearly independent from ϕ M 1 . The loop effects of ϕ Ae g , on the other hand, are tiny and therefore not shown explicitely.
The relative corrections for the process e + e − →ẽ ± 1ẽ ∓ 2 , as shown in Fig. 10, are rather large for the parameter set chosen; see Tab. 1. In the upper left plot of Fig. 10    The dependence on MẼ is shown in the upper right plot of Fig. 10 and follows the same pattern as for e + e − →ẽ + 1ẽ − 1 , i.e. a strong decrease with increasing MẼ as obvious from kinematics. The loop corrections decrease from ∼ +28% at MẼ = 100 GeV to ∼ −31% at MẼ = 470 GeV (the latter is due to the smallness of the tree-level cross section), with a tree crossing at MẼ = 425 GeV.
The Higgs mixing parameter µ dependence is shown in the middle left plot. It is rather linear and decreasing from ∼ 48.6 fb at small µ down to ∼ 46.9 fb at µ = 1000 GeV. The relative corrections are ∼ +17.5 % at µ = 100 GeV and ∼ 16.6 % at µ = 1000 GeV. In this case the (expected) breakdown of the CCN[1] scheme for µ ≈ M 2 = 600 GeV is rather weakly pronounced. The corresponding results are shown for the CNN[1,1,2] and CNN [2,1,3] schemes, which are smooth at µ = M 2 . For small and large values of µ CNN[1,1,2] differs from CCN [1], whereas the other scheme, CNN [2,1,3], is very close to CCN [1]. While the CNN[2,1,3] scheme has a strong (local) minimum of δµ at µ ≈ 421 GeV, dominating the loop corrections, it approximates CCN [1] rather good for all other values of µ.
The dependence on |M 1 | is shown in the middle right plot of Fig. 10. A strong dependence of the tree-level cross section can be observed (from the dominant neutralino t-channel exchange), which is amplified at the one-loop level. The size of the loop corrections varies from ∼ +16.6% at |M 1 | = 350 GeV (i.e. S) to ∼ +15.2% at |M 1 | = 750 GeV and then increase again to ∼ +16.1% at |M 1 | = 1000 GeV.
The The phase dependence ϕ M 1 of the cross section in S is shown in the lower right plot of Fig. 10. In this case it turns out to be substantial, already changing the tree-level cross section by up to 3.7 %. The relative loop corrections (σ loop /σ tree ) vary with ϕ M 1 between ∼ +14.6 % at ϕ M 1 = 180 • and ∼ +17.3 % at ϕ M 1 = 158 • , 203 • for CCN [1]. Again, the cross section does not diverge at ϕ M 1 = 180 • . This structure is caused by a minimum of δM 1 at ϕ M 1 = 180 • , dominating the loop corrections. Using the CNN[2,2,3] scheme (red dashed line) the dip disappears, as can be seen in the plot. The relative loop corrections reach ∼ +17.7 % at ϕ M 1 = 180 • . For our parameter set S, with the complex phase ϕ M 1 , the CP asymmetry turn out to be (numerically) zero, which can be seen from the identical green dashed (ẽ + 1ẽ − 2 ) and yellow dash-dotted (ẽ − 1ẽ + 2 ) lines. Finally the variation with ϕ Ae g is again negligible and not shown here.
We now turn to the process e + e − →ẽ + 2ẽ − 2 shown in Fig. 11, which is found to be of O(10 fb). As a function of √ s (upper left plot) the loop corrections range from ∼ −3 % at √ s = 1000 GeV (i.e. S) to ∼ −10 % at √ s = 3000 GeV.
The cross sections are decreasing with increasing MẼ due to kinematics, and the full one-loop result has its maximum of ∼ 11.6 fb at MẼ = 100 GeV. Analogously the relative corrections are decreasing from ∼ +4 % at MẼ = 100 GeV to ∼ −23% at MẼ = 440 GeV.    parameter variations one can conclude that a cross section twice as large can be possible for very low MẼ (which however are challenged by the current ATLAS/CMS exclusion bounds).
The dependence on µ (middle left plot) is rather small. The one-loop corrections within the CCN [1] scheme are ∼ +2 % at µ = 100 GeV, ∼ −2 % at µ = 1000 GeV and have a tree crossing at µ ≈ 150 GeV. The corresponding smooth 9 results CNN [1,1,2] and CNN [1,2,3] are shown, but both differ for small values (and CNN [1,2,3] also for large values) of µ significantly from CCN [1]. Therefore, in addition we also show the CNN[2,1,3] scheme (black dash-dotted line) which approximates CCN [1] very good for all values of µ except for µ ≈ 421 GeV where δµ has a strong (local) minimum, that dominates the loop corrections in this part of the parameter space.
With increasing |M 1 | in S (middle right plot) we find again a strong decrease of the production cross section, due to the change in the interference of the (dominant) neutralinos in the t-channel. It should be noted that there is no tree crossing in this plot. The loop corrections are rather small and decrease from ∼ −0.6% at |M 1 | = 100 GeV to ∼ −3% at We now turn to the second generation charged slepton production. The process e + e − → µ + 1μ − 1 is shown in Fig. 12, which is found in S at the level of 10 fb, but can be substantially larger by roughly a factor of two for small MẼ; see below. Away from the production threshold, loop corrections of ∼ +14 % at √ s = 1000 GeV (i.e. S) are found. They reach their maximum of ∼ +35 % at √ s = 3000 GeV. The tree crossing takes place at √ s ≈ 725 GeV. The cross section depends strongly on MẼ, as can be seen in the right plot. It is decreasing with increasing MẼ and the full correction has its maximum of ∼ 28 fb at MẼ = 100 GeV. The variation of the relative corrections are rather large, ∼ +33 % at MẼ = 100 GeV, with a tree crossing at MẼ = 440 GeV, and ∼ −25 % at MẼ = 490 GeV where the cross section goes to zero due to kinematics.
The dependence on the remaining parameters is (rather) negligible and therefore we have omitted showing the corresponding plots here.
The process e + e − →μ ± 1μ ∓ 2 is shown in Fig. 13. It should be noted that the smuon (and stau) tree-level process consist of only one Z exchange diagram (see Fig. 1) and is 10 and consequently For vanishing m the cross section can be relatively large (but remains finite). In our scenario S we have chosen a (more realistic) setting with ML = MẼ (with an off-set of m = 50 GeV). Therefore within this scenario the cross section for e + e − →μ ± below the reach of a linear collider. For this reason we refrain from a more detailed discussion here. The only expected exception is the variation with t β (see footnote 10) which we show in the right plot of Fig. 13. The loop corrected cross section increases from O(0.01 ab) at small t β to ∼ 23 ab at t β = 50, as expected. The relative corrections for the t β dependence are increasing from ∼ −12 % at t β = 10 to ∼ −7 % at t β = 50. Now we turn to the process e + e − →μ + 2μ − 2 shown in Fig. 14. As a function of √ s (left plot) we find relative corrections of ∼ −5 % at √ s = 1000 GeV (i.e. S), and ∼ +10 % at √ s = 3000 GeV with a tree crossing at √ s ≈ 1200 GeV.
In the analysis as a function of MẼ (right plot) the cross section is decreasing with increasing MẼ, but can vary roughly by a factor of two w.r.t. S. The full (relative) one-loop correction has its maximum of ∼ 26 fb (∼ +12 %) at MẼ = 100 GeV, decreasing to ∼ 0.1 fb (∼ −40 %) at ML = 440 GeV with a tree crossing at MẼ ≈ 245 GeV.
The dependence on the other parameters is again (rather) negligible and therefore not shown here.
Turning to the third generation charged slepton production, the process e + e − →τ + 1τ − 1 is shown in Fig. 15. As a function of √ s we find loop corrections of ∼ +14 % at √ s = 1000 GeV (i.e. S), a tree crossing at √ s ≈ 725 GeV (where the one-loop corrections are between ±10 % for √ s < ∼ 900 GeV) and ∼ +35 % at √ s = 3000 GeV, very similar to e + e − →μ + 1μ − 1 . In the analysis as a function of MẼ (upper right plot) the cross sections are decreasing with increasing MẼ as obvious from kinematics and the full corrections have their maximum of ∼ 28 fb at MẼ = 100 GeV, more than two times larger than in S. The relative corrections are changing from ∼ +33 % at MẼ = 100 GeV to ∼ −25 % at MẼ = 490 GeV with a tree crossing at MẼ = 415 GeV.
Here we show in the lower left plot of Fig. 15 the dependence on t β . Contrary to other slepton production cross sections analyzed before 11 , σ full (e + e − →τ + 1τ − 1 ) increases with t β . 11 We have omitted showing these plots because the dependence on t β was indeed negligible (with exception The relative corrections for the t β dependence vary between ∼ +14.2 % at t β = 5 and ∼ +13.4 % at t β = 50.
The phase dependence ϕ Ae g of the cross section in S is shown in the lower right plot of Fig. 15. The loop correction increases the tree-level result by ∼ +14 %. The phase dependence of the relative loop correction is very small and found to be below 0.2 %. The variation with ϕ M 1 is negligible and therefore not shown here.
The process e + e − →τ ± 1τ ∓ 2 is shown in Fig. 16. The overall size of this cross section turns out to be rather small, including all analyzed parameter variations, but enhanced w.r.t. e + e − →μ ± 1μ ∓ 2 by about a factor of m 2 τ /m 2 µ ; see Eq. (25). The loop corrections have a noticeable impact, as can be seen in all six panels of Fig. 16, but never lift the cross section above 2 fb.
As a function of √ s (upper left plot) we find relative corrections of ∼ −11 % at √ s = 1000 GeV (i.e. S), and ∼ +19 % at √ s = 3000 GeV with a tree crossing at √ s ≈ 1300 GeV.
As for other slepton production cross sections, σ(e + e − →τ ± 1τ ∓ 2 ) depends strongly on MẼ, where values one order of magnitude larger than in S (with MẼ = 300 GeV) are possible for small MẼ. One can see that the full corrections have their maximum of ∼ 0.8 fb at MẼ = 100 GeV. The relative corrections are decreasing from ∼ +6 % at MẼ = 100 GeV to ∼ −51 % at MẼ = 470 GeV with a tree crossing at MẼ ≈ 170 GeV.
With increasing µ in S (middle left plot) we find a nearly linear increase of the production cross section. While the CCN[1] scheme has its expected singularity 12 at µ = 600 GeV, the of e + e − →μ ± 1μ ∓ 2 ; see above). 12 It should be noted that, as discussed in the beginning of Sect. 4, within the CCN[1] scheme a small divergence in the tree-level result (at µ = M 2 = 600 GeV) is induced by δµ through the shifted scalar tau masses. δµ in Eq. (5a) enters via Eq. (6a) into Eq. (10) from which the slepton shifts are calculated. The divergence is suppressed with m eg (see Eq. (5a)) and therefore not visible in the tree-level results for selectron production; see above. In order to overcome the problem with a divergent tree-level result, we used here the tree-level result of the CNN[1,2,3] scheme, which is free of such a divergence.     at t β = 50.
The phase dependence ϕ Ae g of the cross section in S is shown in the lower right plot of Fig. 16. It is very pronounced and can vary σ full (e + e − →τ ± 1τ ∓ 2 ) from 0.03 fb to 0.3 fb. The (relative) loop corrections are at the level of ∼ 10 % and vary with ϕ Ae g below ±4.5 % w.r.t. the tree cross section. For our parameter set S, with the complex phase ϕ Ae g , the CP asymmetry turn out to be very small, as can be seen in the inlay. Again, here we show in the lower left plot of Fig. 17 the dependence on t β . Contrary to other slepton production cross sections analyzed before, σ full (e + e − →τ + 2τ − 2 ) decreases with t β . The relative corrections for the t β dependence vary between ∼ −5.4 % at t β = 5 and ∼ −4.4 % at t β = 34. A dip (not visible) in the dotted line at t β ≈ 46 is the threshold The phase dependence ϕ Ae g of the cross section in S is shown in the lower right plot of Fig. 17. The loop correction decreases the tree-level result by ∼ −5 %. The phase dependence of the relative loop correction is small and found to be below 0.5 %. The variation with ϕ M 1 is (again) negligible and therefore not shown here.
Overall, for the charged slepton pair production we observed a decreasing cross section ∝ ln (s)/s for the first and ∝ 1/s for the second and third slepton generations for s → ∞; see Ref. [39]. The (loop corrected) cross sections for the slepton pair production can reach a level of O(10 fb), depending on the SUSY parameters, but is very small for the production of two different smuonsμ 1μ2 at the O(10 ab). This renders these processes difficult to observe at an e + e − collider. 13 The full one-loop corrections are very roughly 15 % of the tree-level results, but vary strongly on the size of MẼ and in the case of selectrons also on the size of M 1 and M 2 . Depending on the size of in particular these parameters the loop corrections can be either positive or negative. This shows that the loop corrections, while being large, have to be included point-by-point in any precision analysis. The dependence on ϕ M 1 (ϕ Ae g ) was found at the level of ∼ 15 % (∼ 8 %), but can go up to ∼ 18 % (∼ 14 %) for the extreme cases. The relative loop corrections varied by up to 2 % (4 %) with ϕ M 1 (ϕ Ae g ). Consequently, the complex phase dependence must be taken into account as well. Finally, for our parameter set S the CP asymmetries turn out to be very small, well below ±1 % (hardly measurable in future e + e − collider experiments).
For all parameter choices it was possible to identify at least one renormalization scheme that exhibited a "smooth" behavior. It appears to be possible for any parameter variation to find a combination of schemes that yield a numerically stable (and nearly constant) one-loop level contribution. A detailed analysis of which scheme yields this desired behavior as a function of the underlying SUSY parameters, however, is beyond the scope of this paper.

4.3
The process e + e − →ν gν * g In Figs. 18 -(20) we present the results for scalar neutrino production at e + e − colliders. It should be noted that for s → ∞ decreasing cross sections ∝ ln (s)/s for the first and ∝ 1/s for the second and third slepton generations are expected; see Ref. [40]. We start with the process e + e − →ν eν * e that is shown in Fig. 18. In the analysis of the production cross section as a function of √ s (upper left plot) we find the expected behavior: a strong rise close to the production threshold, followed by a decrease with increasing √ s, where the s-channel dominates. We find a very small shift w.r.t.  Fig. 18. The masses of the electron sneutrinos are governed by ML = MẼ + 50 GeV. Consequently, a strong, nearly linear decrease can be observed from ∼ 140 fb down to zero for MẼ = 450 GeV (i.e. the sneutrino production threshold), as can be expected from kinematics. In scenario S we find a non-negligible decrease of the cross sections from the loop corrections. They start at MẼ = 100 GeV with ∼ −3 % and reach ∼ −53 % at MẼ ≈ 450 GeV. In the latter case these large loop corrections are due to the (relative) smallness of the tree-level results, which goes to zero at the sneutrino production threshold.
The dependence M 2 of the cross section in S is shown in the lower right plot of Fig. 18, where again the chargino t-channel exchange plays a dominant role (χ ± 1 is dominant for M 2 < ∼ 380 GeV andχ ± 2 dominant for M 2 > ∼ 380 GeV). Again, one can see the (expected) breakdown of the CCN[1] scheme for M 2 = µ = 350 GeV. Outside the region M 2 ≈ µ the scheme CCN [1] is expected to be reliable, since each of the three OS conditions is strongly connected to one of the three input parameters, |M 1 |, M 2 and µ. The loop corrections are ∼ −12 % at M 2 = 600 GeV (i.e. S) and ∼ −14 % at M 2 = 1000 GeV, with a tree crossing at M 2 ≈ 130 GeV. The CNN[2,2,3] scheme (red dashed line) is smooth for all values of M 2 shown and a perfect approximation for CCN [1].
Due to the absence of ϕ M 1 in the tree-level production cross section the effect of this complex phase is expected to be small. Correspondingly we find that the phase dependence ϕ M 1 of the cross section in our scenario is tiny. The same holds for the variation with ϕ Ae g , it also remains tiny and unobservable. Therefore we omit showing these two complex phases explicitely.
In Fig. 19 we present the cross sections σ(e + e − →ν µν * µ ). In the analysis as a function of √ s (left plot) we find as before a tiny shift w.r.t. The dependence on the other parameters is (again) negligible and therefore not shown here.
We finish the e + e − →ν gν * g analysis with σ(e + e − →ν τν * τ ) in Fig. 20. The results are very similar as for σ(e + e − →ν µν * µ ) (see above), since the sneutrino masses are the same for all three generations, as are the contributing higher-order diagrams.
Overall, for the sneutrino pair production we observed a decreasing cross section ∝ ln (s)/s for the first and ∝ 1/s for the second and third slepton generations for s → ∞; see Ref. [40]. The full one-loop corrections are very roughly 10 % of the tree-level results, but depend strongly on the size of MẼ, where larger values result even in negative loop corrections. The cross sections are largest for e + e − →ν eν * e and roughly smaller by one order of magnitude for e + e − →ν µν * µ and e + e − →ν τν * τ . This is caused by the absence of the chargino t-channel diagram (i.e. aχ ∓ c e ±ν g coupling), which only contributes to e + e − →ν eν * e ; see Fig. 1. The variation of the cross sections with ϕ M 1 and ϕ Ae g is found extremely small and have not been shown explicitely.
As for charged slepton production, for all parameter choices it was possible to identify at least one renormalization scheme that exhibited a "smooth" behavior. Also for sneutrino production it appears to be possible for any parameter variation to find a combination of schemes that yield numerically stable (and nearly constant) one-loop level contributions.

Conclusions
We have evaluated all slepton production modes at e + e − colliders with a two-particle final state, i.e. e + e − →ẽ ± gsẽ ∓ gs and e + e − →ν gν * g allowing for complex parameters. In the case of discovery of sleptons a subsequent precision measurement of their properties will be crucial to determine their nature and the underlying (SUSY) parameters. In order to yield sufficient accuracy, one-loop corrections to the various slepton production modes have to be considered. This is particularly the case for the anticipated high accuracy of the slepton property determination at e + e − colliders [21].
The evaluation of the processes (1) and (2) is based on a full one-loop calculation, also including hard, soft and collinear QED radiation. The renormalization is chosen to be identical as for the slepton decay calculations [27], or slepton production from heavy Higgsboson decay [28], as well as from chargino and neutralino decays [29][30][31][32]. Consequently, the predictions for the production and decay can be used together in a consistent manner (e.g., in a global phenomenological analysis of the slepton sector at the one-loop level).
We first briefly reviewed the relevant sectors including some details of the one-loop renormalization procedure of the cMSSM, which are relevant for our calculation. In most cases we follow Ref. [48]. We have discussed the calculation of the one-loop diagrams, the treatment of UV, IR, and collinear divergences that are canceled by the inclusion of (hard, soft, and collinear) QED radiation. As far as possible we have checked our result against the literature, and in most cases where a meaningful comparison could be performed we found good agreement; parts of the differences can be attributed to problems with input parameters and/or different renormalization schemes (conversions).
For the analysis we have chosen a standard parameter set (see Tab. 1), which allows for the production of all combinations of sleptons at an e + e − collider with a center-of-mass energy up to √ s = 1000 GeV. In the analysis we investigated the variation of the various production cross sections with the center-of-mass energy √ s, the Higgs mixing parameter µ, the gaugino mass parameters |M 1 | and/or M 2 , the slepton soft SUSY-breaking parameter MẼ (ML = MẼ + 50 GeV) and the complex phases ϕ Ae g (of the trilinear Higgs-slepton coupling, A eg ; g = 1, 2, 3) and ϕ M 1 (of the gaugino mass parameter M 1 ), respectively. Where relevant we also showed the variation with t β . In our numerical scenarios we compared the tree-level production cross sections with the full one-loop corrected cross sections. The numerical results we have shown are, of course, dependent on the choice of the SUSY parameters. Nevertheless, they give an idea of the relevance of the full one-loop corrections. For the slepton pair production, e + e − →ẽ ± gsẽ ∓ gs and e + e − →ν gν * g we observed for s → ∞ a decreasing cross section ∝ ln (s)/s for the first and ∝ 1/s for the second and third slepton generations. The (loop corrected) cross sections for the charged slepton pair production can reach a level of 10 fb, depending on the SUSY parameters, but is very small for the production of two different smuonsμ 1μ2 at the O(10 ab) (and with the cross section ofτ 1τ2 enhanced by a factor of m 2 τ /m 2 µ ). This renders these processes difficult to observe at an e + e − collider. 14 The full one-loop corrections are very roughly 15 % of the tree-level results, but vary strongly on the size of MẼ and in the case of selectrons also on the size of M 1 and M 2 . Depending on the size of in particular these parameters the loop corrections can be either positive or negative. The dependence on ϕ M 1 (ϕ Ae g ) was found at the level of ∼ 15 % (∼ 8 %), but can go up to ∼ 18 % (∼ 14 %) for the extreme cases. The relative loop corrections varied by up to 2 % (4 %) with ϕ M 1 (ϕ Ae g ). This shows that the loop corrections, including the complex phase dependence, have to be included point-by-point in any precision analysis, or any precise determination of SUSY parameters from the production of cMSSM sleptons at e + e − linear colliders.
Concerning scalar neutrino production, the full one-loop corrections are very roughly 10 % of the tree-level results, but depend strongly on the size of MẼ, where larger values result even in negative loop corrections. The cross sections are largest for e + e − →ν eν * e and roughly smaller by one order of magnitude for e + e − →ν µν * µ and e + e − →ν τν * τ . The variation of the cross sections with ϕ M 1 and ϕ Ae g is found extremely small. Also for scalar neutrino production the loop corrections have to be included point-by-point in any precision analysis at e + e − linear colliders.
For all cross section calculations and for all parameter choices it was possible to identify at least one renormalization scheme that exhibited a "smooth" behavior. It appears to be possible for any parameter variation to find a combination of schemes that yield numerically stable (and nearly constant) one-loop level contributions. A detailed analysis of which scheme yields this desired behavior as a function of the underlying SUSY parameters, however, is beyond the scope of this paper.