format: text

text:
# Exact Cross Sections for the Neutralino WIMP Pair-Annihilation

## 1 Introduction

The quest for identifying the nature of the dark matter (DM) in the Universe continues . It is generally believed that most of the DM is made of some hypothetical weakly-interacting massive particles (WIMPs). From the particle theory point of view, a commonly considered candidate for the WIMP is the lightest neutralino under the assumption that it is the lightest supersymmetric particle (LSP). In most approaches the LSP is stable due to an additional R-parity . The neutralino, being massive, often provides a sizeable contribution to the relic density. In addition, the requirement that the neutralino, or some other stable particle relic, does not “overclose” the Universe, often provides a strong constraints on a supersymmetric model. The two other robust candidates for the LSP and cold dark matter (CDM) are the axino (superpartner of the axion)  and the gravitino .

Continuing improvements in determining the abundance of cold dark matter, and other components of the Universe, have now reached an unprecedented precision of a few per cent . In light of this, one needs to be able to perform an accurate enough computation of the WIMP relic abundance, which would allow for a reliable comparison between theory and observation.

The literature on the relic abundance of the neutralino is vast and still growing. (For a comprehensive review, see Ref. .) A brief, and by no means complete, account of major developments, can be summarized as follows. The original paper by Goldberg  considered the neutralino in the photino limit and pointed out the strong constraints from its relic abundance. This was soon followed the first analysis by Ellis, et al.  and Krauss , of the general neutralino case. Several other early papers subsequently appeared with more detailed and elaborate analyses. In particular, Griest  was first to compute in detail the annihilation into the ordinary fermion-pair ($`f\overline{f}`$) final states through the $`Z`$–exchange, and later Griest, Kamionkowski and Turner  conducted the first more complete analysis of the general neutralino case into $`WW`$, $`ZZ`$ and Higgs-pair final states. The Higgs contribution to $`f\overline{f}`$ was first computed in Refs. . Olive and Srednicki  considered all the annhilation channels but only in the limit of the pure gaugino and higgsino cases where several important resonances and final states are absent. Drees and Nojiri  computed a first complete set of expressions for the product of the cross section times velocity using the helicity amplitude technique. When expanded in the nonrelativistic limit, these give expressions for the first two coefficients of the partial wave expansion.

In the early papers the partial wave expansion of the thermally-averaged product of the neutralino pair-annihilation cross section and their relative velocity, $`\sigma va+bx`$ was used in most cases. The method is normally expected to give an accurate enough approximation (about few per cent) but only far enough from $`s`$–channel resonances and thresholds for new final states, as was first pointed out by Griest and Seckel  and further emphasized in Refs. . In particular, it was shown in Ref.  that, because of the very narrow width of the lightest supersymmetric Higgs $`h`$, in the vicinity of its $`s`$–channel exchange the error can be as large as a few orders of magnitude. Partial remedies were suggested, for example in , by numerically integrating the full cross section near resonances only and by matching this with the expansion-based calculations further away from poles. Such methods are not fully satisfactory since they do not include interference terms. A recent detailed analysis  showed that in the case of the often wide $`s`$–channel resonance caused by the exchange of the pseudoscalar $`A`$, the expansion produces a significant error over the range of neutralino mass which can be as big as several tens of GeV. Furthermore subdominant channels and often neglected interference terms can also sometimes play a sizeable role.

A formalism for computing the relic abundance also became refined. In particular, the effect of replacing the usually assumed common heat bath for both annihilating particles  by a more accurate treatment of involving two separate thermal distributions was considered by Gondolo and Gelmini  and by Srednicki, et al. . In practical terms the numerical difference is usually negligible when one does the usual partial wave expansion. The first coefficient $`a`$ is universal while the second ones $`b`$ differ by $`3/2a`$, where usually $`ab`$. Gondolo and Gelmini  further derived a very useful compact expression for the thermally-averaged product of the neutralino pair-annihilation cross section and their relative velocity as a single integral over the cross section, as we will see below. Gondolo and collaborators next developed a Fortran code DarkSusy  where the relic density of neutralinos is numerically computed without using the partial wave expansion approximation.

An additional effect of reducing the relic abundance of WIMPs through co-annihilation was first pointed out by Griest and Seckel . In some cases there may exist some other states which are not much heavier than the stable WIMP and may therefore be still present in the thermal plasma around the WIMP decoupling. In the framework of minimal supersymmetry with the lightest neutralino LSP with a significant higgsino component, the co-annihilation with the next-to-lightest neutralino and the lightest chargino is often important . Other cases of interest involve neutralino annihilation with the lighter superpartner of the $`\tau `$–lepton  and with lighter stop .

In this paper, we will present a full set of exact analytic expressions for the cross section of the neutralino pair-annihilation in the general Minimal Supersymmetric Standard Model (MSSM) for a general neutralino case. From the point of view of low-energy supersymmetry, the most natural choice for the LSP and CDM is a nearly pure gaugino (bino) as was first shown by Roszkowski . Remarkably, just such a case of the \[electrically neutral\] LSP naturally emerges in most case in the Constrained MSSM . Nevertheless, in our analysis we will make no simplifying assumptions about the neutralino, nor will we assume the degeneracy of the left– and right–sfermion masses. We will include all tree-level final states and all intermediate states. We will also keep finite widths in $`s`$–channel resonances. We will only neglect possible CP-violating phases in the SUSY sector. We will also not consider the effect of co-annihilation here but will address it in a subsequent publication. A complete set of expressions presented here does not rely on the partial wave expansion, includes all the terms and is valid both near and further away from resonances and thresholds for new final states.

Some of the results presented here are not new but we include them here nevertheless in order to provide a complete and self-consistent reference containing the full set of exact expressions. In particular, the cross sections for the neutralino pair-annihilation into the SM fermion-pair ($`f\overline{f}`$) final states were first computed in Ref  and for the $`WW`$, $`ZZ`$ and Higgs boson-pair final states in Ref. .

Given the complexity of analytic expressions presented here and in other papers and due to often different conventions used, it is not doable to perform a real comparison with the literature. Instead, we have performed a numerical check of our cross section with the results obtained by using DarkSusy . We have found, for the same values of input parameters, an impressive agreement, at the level of a few per cent, for all the annihilation channels, which we find reassuring. (Recently another numerical code has been derived  and also numerically agreed with DarkSusy.)

While the exact analytic expressions presented here are applicable both near and away from special cases where the method of partial wave expansion fails, sometimes one may find it less CPU-time consuming to use the latter one. Starting from our exact cross sections we have therefore derived a complete set of expressions for the usual first two terms of expansion for all the dominant channels. We will present them here as well.

The annihilation into $`f\overline{f}`$ is often dominant. However, other final states can also play important role, depending on the case. In a previous paper  we performed a detailed numerical comparison of the relic abundance computed using the exact formulae with the one obtained using the expansion formulae, for all the channels, including subdominant ones. Our analysis confirmed that the expansion gives highly inaccurate results near resonances and new thresholds. We also showed that very far from such cases the error is typically rather small, of the order of a few per cent. However, we found that, because of the existence of several resonances ($`Z`$ and the Higgs bosons), the expansion produces large errors, compared to an exact treatment, over a sizeable range of the neutralino mass $`m_\chi `$, even of a several tens of GeV. In other words, the widely used method of expansion may lead to significant errors in a sizeable fraction of the neutralino mass.

The plan of the paper is as follows. In Sect. 2 we review the formalism for computing the relic density that we employ. In Sect. 3 we introduce the relevant ingredients of the MSSM and list all the neutralino pair-annihilation channels. Explicit expressions for the annihilation cross secion are given in Sect. 4. In Sect. 5 we discuss expansion and provide a list of formulae for the first two coefficients in the case of equal-mass final states. In Sect. 6 we summarize our work. Appendix A contains a complete list of Lagrangian terms and couplings which are used in the paper while in Appendix B we provide expressions for several auxiliary functions used in the text.

## 2 Calculation of the Relic Density

The relic abundance of some stable species $`\chi `$ is defined as $`\mathrm{\Omega }_\chi \rho _\chi /\rho _{crit}`$, where $`\rho _\chi =m_\chi n_\chi `$ is the relic’s mass density, $`n_\chi `$ is its number density, $`\rho _{crit}3H_0^2/8\pi G_N=1.9\times 10^{29}(h^2)g/cm^3`$ is the critical density and $`G_N`$ is the gravitational constant.(For a review of relic density calculations, see, e.g., Refs. .) The time evolution and subsequent freeze-out of $`n_\chi `$ in an expanding Universe are described by the Boltzmann equation

$`{\displaystyle \frac{dn_\chi }{dt}}`$ $`=`$ $`3Hn_\chi \sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}\left[n_\chi ^2(n_\chi ^{\mathrm{eq}})^2\right],`$ (1)
where $`n_\chi ^{\mathrm{eq}}`$ is the number density that the species would have in thermal equilibrium, $`H(T)`$ is the Hubble expansion rate, $`\sigma (\chi \chi \mathrm{all})`$ denotes the cross section of the species annihilation into ordinary particles, $`v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$ is a so-called Møller velocity  which is the relative velocity of the annihilating particles, and $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$ represents the thermal average of $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$ which will be given below. In the early Universe, the species $`\chi `$ were initially in thermal equilibrium, $`n_\chi `$ $`=n_\chi ^{\mathrm{eq}}`$. When their typical interaction rate $`\mathrm{\Gamma }_\chi `$ became less than the Hubble parameter, $`\mathrm{\Gamma }_\chi \text{ }<H`$, the annihilation process froze out. Since then their number density in a co-moving volume has remained basically constant.

The thermally-averaged product of the neutralino pair-annihilation cross section and their relative velocity $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$ is most properly defined in terms of separate thermal baths for both annihilating particles

$$\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}(T)=\frac{d^3p_1d^3p_2\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}e^{E_1/T}e^{E_2/T}}{d^3p_1d^3p_2e^{E_1/T}e^{E_2/T}}$$
(2)
where $`p_1=(E_1,𝐩_1)`$ and $`p_2=(E_2,𝐩_2)`$ are the 4-momenta of the two colliding particles, and $`T`$ is the temperature of the bath.(Note that one often uses another definition of $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$ which involves a single thermal bath for both neutralinos. Compare, e.g., Refs. . The numerical difference between the two formulae is usually rather small.) The above expression can be reduced to a one-dimensional integral which can be written in a Lorentz-invariant form as

$`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}(T)`$ $`=`$ $`{\displaystyle \frac{1}{8m_\chi ^4TK_2^2(m_\chi /T)}}{\displaystyle _{4m_\chi ^2}^{\mathrm{}}}𝑑s\sigma (s)(s4m_\chi ^2)\sqrt{s}K_1\left({\displaystyle \frac{\sqrt{s}}{T}}\right),`$ (3)
where $`s=(p_1+p_2)^2`$ is a usual Mandelstam variable and $`K_i`$ denotes the modified Bessel function of order $`i`$. In computing the relic abundance one first evaluates eq. (3) and then uses this to solve the Boltzmann eq. (1).

There are a number of methods of solving eq. (1). One often used, approximate, although in general quite accurate (for a recent discussion see Ref. ), solution to the Boltzmann equation is based on solving iteratively the equation

$`x_f^1`$ $`=`$ $`\mathrm{ln}\left({\displaystyle \frac{m_\chi }{2\pi ^3}}\sqrt{{\displaystyle \frac{45}{2g_{}G_N}}}\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}(x_f)x_f^{1/2}\right),`$ (4)
where $`g_{}`$ represents the effective number of degrees of freedom at freeze-out ($`\sqrt{g_{}}9`$). Typically one finds that the freeze-out point $`x_fT_f/m_\chi `$ is roughly given by $`1/25`$$`1/20`$. One usually introduces $`J(x_f)`$ defined as

$$J(x_f)_0^{x_f}𝑑x\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}(x),$$
(5)
where $`x=T/m_\chi `$.

The relic density at present is given by

$`\rho _\chi `$ $`=`$ $`{\displaystyle \frac{1.66}{M_{\mathrm{Pl}}}}\left({\displaystyle \frac{T_\chi }{T_\gamma }}\right)^3T_\gamma ^3\sqrt{g_{}}{\displaystyle \frac{1}{J(x_f)}},`$ (6)
where $`M_{\mathrm{Pl}}=1/\sqrt{G_N}`$ denotes the Planck mass, $`T_\chi `$ and $`T_\gamma `$ are the present temperatures of the neutralino and the photon, respectively. The suppression factor $`(T_\chi /T_\gamma )^3`$ $``$ $`1/20`$ follows from entropy conservation in a comoving volume .

## 3 WIMP Annihilation in the MSSM

In this Section we introduce the relevant parameters and definitions. We will be working in the framework of the general MSSM.(For a review, see, e.g., Ref. . We follow the conventions of Ref. .) The lightest neutralino is a mass eigenstate given by a linear combination of the bino $`\stackrel{~}{B}`$, the neutral wino $`\stackrel{~}{W}_3^0`$ and the two neutral higgsinos $`\stackrel{~}{H}_b^0`$ and $`\stackrel{~}{H}_t^0`$

$`\chi `$ $`\chi _1^0=`$ $`N_{11}\stackrel{~}{B}+N_{12}\stackrel{~}{W}_3^0+N_{13}\stackrel{~}{H}_b^0+N_{14}\stackrel{~}{H}_t^0.`$ (7)
The neutralino mass matrix is determined by the $`U(1)_Y`$ and $`SU(2)_L`$ gaugino mass parameters $`M_1`$ and $`M_2`$, respectively (and we impose the usual GUT relation $`M_1=\frac{5}{3}\mathrm{tan}^2\theta _WM_2`$), the Higgs/higgsino mass parameter $`\mu `$, the usual weak angle $`\theta _W`$ and $`\mathrm{tan}\beta =v_t/v_b`$ – the ratio of the vacuum expectation values of the two neutral Higgs fields.

The neutralino mass matrix is given by

$`_{\chi ^0}`$ $`=`$ $`\left(\begin{array}{cccc}M_1& 0& m_Z\mathrm{sin}\theta _W\mathrm{cos}\beta & m_Z\mathrm{sin}\theta _W\mathrm{sin}\beta \\ 0& M_2& m_Z\mathrm{cos}\theta _W\mathrm{cos}\beta & m_Z\mathrm{cos}\theta _W\mathrm{sin}\beta \\ m_Z\mathrm{sin}\theta _W\mathrm{cos}\beta & m_Z\mathrm{cos}\theta _W\mathrm{cos}\beta & 0& \mu \\ m_Z\mathrm{sin}\theta _W\mathrm{sin}\beta & m_Z\mathrm{cos}\theta _W\mathrm{sin}\beta & \mu & 0\end{array}\right).`$ (12)
The neutralino mass matrix is diagonalized by a unitary matrix $`N`$

$`N^{}_{\chi ^0}N^1`$ $`=`$ $`\mathrm{diag}(m_{\chi _1^0},m_{\chi _2^0},m_{\chi _3^0},m_{\chi _4^0}).`$ (14)
In the absence of possible CP violating phases, one can choose a basis such that the mixing matrix $`N`$ is real, in which case some of the neutralino masses will in general be negative.

The chargino mass matrix is given by

$`_{\chi ^\pm }`$ $`=`$ $`\left(\begin{array}{cc}M_2& \sqrt{2}m_W\mathrm{sin}\beta \\ \sqrt{2}m_W\mathrm{cos}\beta & \mu \end{array}\right).`$ (17)
The chargino mass matrix is diagonalized by two unitary matrices $`U`$ and $`V`$

$`U^{}_{\chi ^\pm }V^1`$ $`=`$ $`\mathrm{diag}(m_{\chi _1^\pm },m_{\chi _2^\pm }).`$ (18)
There are two neutral scalar Higgs bosons $`h`$ and $`H`$, a pseudoscalar $`A`$ plus a pair of charged Higgs $`H^\pm `$. (We will typically suppress the Higgs charge assignment except where this may lead to ambiguities.)

Other relevant parameters which determine the masses of scalars and various couplings are the squark soft mass parameters $`m_Q`$, $`m_U`$ and $`m_D`$, the slepton soft mass parameters $`m_L`$ and $`m_E`$, and the pseudoscalar mass $`m_A`$. We also include the trilinear terms $`A_i`$ ($`i=t,b,\tau `$) of the third generation which are important in determining the masses and couplings for the stop $`\stackrel{~}{t}_{1,2}`$, sbottom $`\stackrel{~}{b}_{1,2}`$ and stau $`\stackrel{~}{\tau }_{1,2}`$ states, respectively.

In general the flavor-violating sfermion mass–squared ($`6\times 6`$) matrices are given by

$`_{\stackrel{~}{f}}^2`$ $`=`$ $`\left(\begin{array}{cc}m_{LL}^2(\stackrel{~}{f})& m_{LR}^2(\stackrel{~}{f})\\ m_{RL}^2(\stackrel{~}{f})& m_{RR}^2(\stackrel{~}{f})\end{array}\right),`$ (21)
with

$`m_{LL}^2(\stackrel{~}{f})`$ $`=`$ $`M_{\stackrel{~}{f}L}^2+M_f^{}M_f+m_Z^2\mathrm{cos}2\beta \left(T_{3f}Q_f\mathrm{sin}^2\theta _W\right)\mathrm{𝟏},`$ (22)
$`m_{RR}^2(\stackrel{~}{f})`$ $`=`$ $`M_{\stackrel{~}{f}R}^2+M_fM_f^{}+m_Z^2\mathrm{cos}2\beta Q_f\mathrm{sin}^2\theta _W\mathrm{𝟏},`$ (23)
$`m_{LR}^2(\stackrel{~}{f})`$ $`=`$ $`\{\begin{array}{c}M_f^{}\left(A_f^{}\mu \mathrm{cot}\beta \right)\mathrm{for}T_{3f}=+1/2\\ M_f^{}\left(A_f^{}\mu \mathrm{tan}\beta \right)\mathrm{for}T_{3f}=1/2\end{array},`$ (26)
$`m_{RL}^2(\stackrel{~}{f})`$ $`=`$ $`m_{LR}^2(\stackrel{~}{f}),`$ (27)
where $`\stackrel{~}{f}`$ and $`f`$ denote here different types of (s)fermions (up– and down–type (s)quarks, charged and neutral (s)leptons). For (s)neutrinos only the $`LL`$ part of eq. (21) should be taken. $`M_{\stackrel{~}{f}L}^2`$ and $`M_{\stackrel{~}{f}R}^2`$ denote $`3\times 3`$ soft SUSY–breaking sfermion mass matrices, and $`M_f`$ denotes here a $`3\times 3`$ fermion mass matrix. Finally, $`A_f`$ is a scalar trilinear coupling matrix of the same dimension while $`Q_f`$ and $`T_{3f}`$ are the respective electric and isospin charges. All the interaction terms and couplings that we will need below are summarized in Appendix A.

In the MSSM, the neutralino LSP’s can pair-annihilate into a number of final states, if kinematically allowed. A complete list of all tree-level two-body final states is given in Table 1. We only neglect two-body loop processes into final state photon pairs and gluon pairs because they are always subdominant in computing the relic density . We also neglect three-body final states since they are unlikely to be competitive with two-body final states. They were shown to dominate in the higgsino case just below the $`WW`$ and $`t\overline{t}`$ final states  but in such regions neutralino co-annihilation with the lightest chargino and next-to-lightest neutralino reduce the relic density to very small values anyway.

The channels $`WW`$, $`ZZ`$, $`\overline{t}t`$, $`W^\pm H^{}`$, $`Zh`$, $`ZH`$, $`Ah`$ and $`AH`$, are not $`s`$-wave suppressed, and, once kinematically allowed, can give dominant contributions. But even the $`s`$-wave suppressed channels $`\overline{f}f`$ ($`ft`$), $`hh`$, $`Hh`$, $`HH`$, $`AA`$, $`H^+H^{}`$ and $`ZA`$ can play some role, especially if the other channels are not yet kinematically allowed. This in particular is the case with the light fermion-pair final states for which the cross section is suppressed by the square of the corresponding fermion mass but which are always kinematically allowed and often dominant.

## 4 Exact Expressions

We now proceed to present a full set of exact, analytic expressions for the total cross section $`\sigma (\chi \chi \mathrm{all})`$ for the neutralino pair-annihilation processes into all allowed (tree-level) two-body final states in the general MSSM. We have included all contributing diagrams as well as all interference terms and kept finite widths of all $`s`$–channel resonances. We have made no simplifying assumptions about sfermion masses although we assumed that there are no CP violating phases in SUSY parameters.

To start with, it is convenient to introduce a Lorentz-invariant function $`w(s)`$

$$w(s)=\frac{1}{4}𝑑\mathrm{LIPS}|𝒜(\chi \chi \mathrm{all})|^2$$
(28)
where $`|𝒜(\chi \chi \mathrm{all})|^2`$ denotes the absolute square of the reduced matrix element for the annihilation of two $`\chi `$ particles, averaged over initial spins and summed over final spins. The function $`w(s)`$ is related to the annihilation cross section $`\sigma (s)`$ in eq. (3) via

$$w(s)=\frac{1}{2}\sqrt{s(s4m_\chi ^2)}\sigma (s).$$
(29)
Since $`w(s)`$ receives contributions from all the kinematically allowed annihilation process $`\chi \chi f_1f_2`$, it can be written as

$`w(s)`$ $`=`$ $`{\displaystyle \frac{1}{32\pi }}{\displaystyle \underset{f_1f_2}{}}\left[c\theta \left(s(m_{f_1}+m_{f_2})^2\right)\beta _f(s,m_{f_1},m_{f_2})\stackrel{~}{w}_{f_1f_2}(s)\right],`$ (30)
where the summation extends over all possible two-body final states $`f_1f_2`$, $`m_{f_1}`$ and $`m_{f_2}`$ denote their respective masses, and

$`c`$ $`=`$ $`\{\begin{array}{cc}c_f\hfill & \text{if }f_{1(2)}=f(\overline{f})\hfill \\ 1\hfill & \text{otherwise,}\hfill \end{array}`$ (33)
where $`c_f`$ is the color factor of SM fermions ($`c_f=3`$ for quarks and $`c_f=1`$ for leptons). The kinematic factor $`\beta _f`$ is defined as

$`\beta _f(s,m_{f_1},m_{f_2})\left[1{\displaystyle \frac{(m_{f_1}+m_{f_2})^2}{s}}\right]^{1/2}\left[1{\displaystyle \frac{(m_{f_1}m_{f_2})^2}{s}}\right]^{1/2}.`$ (34)
In the CM frame, which we choose for convenience, the function $`\stackrel{~}{w}_{f_1f_2}(s)`$ can be expressed as

$`\stackrel{~}{w}_{f_1f_2}(s){\displaystyle \frac{1}{2}}{\displaystyle _1^{+1}}d\mathrm{cos}\theta _{CM}|𝒜(\chi \chi f_1f_2)|^2,`$ (35)
where $`\theta _{CM}`$ denotes the scattering angle in the CM frame. In other words, we write $`|𝒜(\chi \chi f_1f_2)|^2`$ as a function of $`s`$ and $`\mathrm{cos}\theta _{CM}`$, which greatly simplifies the computation.

We will follow Table 1 in presenting explicit expressions for $`\stackrel{~}{w}_{f_1f_2}(s)`$ for all the two-body final states. All the couplings are defined in Appendix A. All other auxiliary functions, are listed in Appendix B. Feynman diagrams corresponding to all the annihilation channels are given in Chapter 6 of Ref. .

$`\mathrm{𝟏}\mathbf{.}𝝌𝝌\mathbf{}𝒉𝑯`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{hH}=\stackrel{~}{w}_{hH}^{(h,H)}+\stackrel{~}{w}_{hH}^{(\chi ^0)}+\stackrel{~}{w}_{hH}^{(h,H\chi ^0)}:`$ (36)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{w}_{hH}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{hHr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2);`$ (37)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{w}_{hH}^{(\chi ^0)}={\displaystyle \underset{i,j=1}{\overset{4}{}}}C_S^{\chi _i^0\chi h}(C_S^{\chi _j^0\chi h})^{}C_S^{\chi _i^0\chi H}(C_S^{\chi _j^0\chi H})^{}\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{hH}+m_\chi m_{\chi _i^0}J_{ij}^{hH}+K_{ij}^{hH}\right],`$ (38)
where

$`I_{ij}^{hH}`$ $`=`$ $`(s4m_\chi ^2)\left[𝒯_0𝒴_0\right](s,m_\chi ^2,m_h^2,m_H^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$ (39)
$`J_{ij}^{hH}`$ $`=`$ $`[4𝒯_1+2(2m_\chi ^2+m_h^2+m_H^2)𝒯_02𝒴_1`$ (40)
$`2(s4m_\chi ^2)𝒴_0](s,m_\chi ^2,m_h^2,m_H^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`K_{ij}^{hH}`$ $`=`$ $`\{T_2(s+2m_\chi ^2m_h^2m_H^2)𝒯_1(m_\chi ^2m_h^2)(m_\chi ^2m_H^2)𝒯_0`$
$`𝒴_2+[(m_\chi ^2+m_h^2)(m_\chi ^2+m_H^2)+4m_\chi ^4]𝒴_0\}(s,m_\chi ^2,m_h^2,m_H^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2);`$
$``$ Higgs ($`h,H`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{hH}^{(h,H\chi ^0)}`$ $`=`$ $`2{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{hHr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}C_S^{\chi _i^0\chi h}C_S^{\chi _i^0\chi H}\right]`$ (42)
$`\times \{2m_\chi +[m_\chi (m_H^2+m_h^22m_\chi ^22m_{\chi _i^0}^2)`$
$`+m_{\chi _i^0}(s4m_\chi ^2)](s,m_\chi ^2,m_H^2,m_h^2,m_{\chi _i^0}^2)\}.`$
As mentioned above, the couplings $`C_S^{\chi _i^0\chi _j^0r}`$ and $`C^{hHr}`$ ($`r`$ $`=`$ $`h,H`$), as well as all the other couplings appearing in this Section, are defined in Appendix A. The functions $``$, $`𝒯_k`$ and $`𝒴_k`$ ($`k`$ $`=`$ $`0,\mathrm{},4`$), and all other auxiliary functions, are listed in Appendix B. By $`\mathrm{\Gamma }_h`$ and $`\mathrm{\Gamma }_H`$ we denote the widths of $`h`$ and $`H`$, respectively.

The expressions for $`hh`$ final state are obtained from the above by replacing $`m_H`$, $`C^{hHr}`$, $`C_S^{\chi _i^0\chi H}`$, $`C_S^{\chi _j^0\chi H}`$ with $`m_h`$, $`C^{hhr}`$, $`C_S^{\chi _i^0\chi h}`$, $`C_S^{\chi _j^0\chi h}`$, respectively, and multiplying $`\omega `$ by a factor of 1/2 for identical particles in the final state. The contributions for $`HH`$ final state are obtained in an analogous way.

$`\mathrm{𝟐}\mathbf{.}𝝌𝝌\mathbf{}𝑨𝑨`$

Similarly to the final state $`hH`$, this process proceeds via the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{AA}=\stackrel{~}{w}_{AA}^{(h,H)}+\stackrel{~}{w}_{AA}^{(\chi ^0)}+\stackrel{~}{w}_{AA}^{(h,H\chi ^0)}:`$ (43)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{w}_{AA}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{1}{4}}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{AAr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2);`$ (44)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{w}_{AA}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}\left(C_P^{\chi _i^0\chi A}\right)^2\left(C_P^{\chi _j^0\chi A}\right)^2\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{AA}+m_\chi m_{\chi _i^0}J_{ij}^{AA}+K_{ij}^{AA}\right],`$ (45)
where

$`I_{ij}^{AA}`$ $`=`$ $`(s4m_\chi ^2)\left[𝒯_0Y_0\right](s,m_\chi ^2,m_A^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$ (46)
$`J_{ij}^{AA}`$ $`=`$ $`\left[4𝒯_1+4(m_\chi ^2m_A^2)𝒯_0+2𝒴_1+2(s4m_\chi ^2)𝒴_0\right](s,m_\chi ^2,m_A^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`K_{ij}^{AA}`$ $`=`$ $`\{T_2[s+2(m_\chi ^2m_A^2)]𝒯_1(m_\chi ^2m_A^2)^2𝒯_0`$ (48)
$`𝒴_2+[4m_\chi ^4(m_\chi ^2+m_A^2)^2]𝒴_0\}(s,m_\chi ^2,m_A^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2);`$
$``$ Higgs ($`h,H`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{AA}^{(h,H\chi ^0)}`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{AAr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\left(C_P^{\chi _i^0\chi A}\right)^2\right]`$ (49)
$`\times \{2m_\chi +[2m_\chi (m_A^2m_\chi ^2m_{\chi _i^0}^2)`$
$`m_{\chi _i^0}(s4m_\chi ^2)](s,m_\chi ^2,m_A^2,m_A^2,m_{\chi _i^0}^2)\}.`$
$`\mathrm{𝟑}\mathbf{.}𝝌𝝌\mathbf{}𝒉𝑨`$

This process proceeds via the $`s`$-channel $`Z`$ and CP–odd Higgs boson ($`A`$) exchange as well as the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{hA}`$ $`=`$ $`\stackrel{~}{w}_{hA}^{(A)}+\stackrel{~}{w}_{hA}^{(Z)}+\stackrel{~}{w}_{hA}^{(\chi ^0)}+\stackrel{~}{w}_{hA}^{(AZ)}+\stackrel{~}{w}_{hA}^{(A\chi ^0)}+\stackrel{~}{w}_{hA}^{(Z\chi ^0)}:`$ (50)
$``$ CP–odd Higgs-boson ($`A`$) exchange:

$`\stackrel{~}{w}_{hA}^{(A)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \frac{C^{hAA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2s;`$ (51)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{w}_{hA}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{3}}\left|{\displaystyle \frac{C^{hAZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2{\displaystyle \frac{1}{m_Z^4s}}`$ (52)
$`\times \{m_Z^4s^3+s^2[6m_\chi ^2(m_A^2m_h^2)^24m_\chi ^2m_Z^42m_Z^4(m_A^2+m_h^2)]`$
$`+sm_Z^2[12m_\chi ^2(m_A^2m_h^2)^2+8m_Z^2m_\chi ^2(m_A^2+m_h^2)`$
$`+m_Z^2(m_A^2m_h^2)^2]+2m_Z^4m_\chi ^2(m_A^2m_h^2)^2\};`$
$``$ neutralino ($`\chi _i^0)`$ exchange:

$`\stackrel{~}{w}_{hA}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \underset{i,j=1}{\overset{4}{}}}C_S^{\chi _i^0\chi h}C_P^{\chi _i^0\chi A}\left(C_S^{\chi _j^0\chi h}\right)^{}\left(C_P^{\chi _j^0\chi A}\right)^{}\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{hA}+m_\chi m_{\chi _i^0}J_{ij}^{hA}+K_{ij}^{hA}\right],`$
where

$`I_{ij}^{hA}`$ $`=`$ $`s\left[𝒯_0𝒴_0\right](s,m_\chi ^2,m_h^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$ (54)
$`J_{ij}^{hA}`$ $`=`$ $`2(m_A^2m_h^2)\left[𝒯_0𝒴_0\right](s,m_\chi ^2,m_h^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$ (55)
$`K_{ij}^{hA}`$ $`=`$ $`\{T_2+(2m_\chi ^2+m_A^2+m_h^2s)𝒯_1+𝒴_2`$
$`+\{m_\chi ^2(m_A^2+m_h^2m_\chi ^2)m_A^2m_h^2\}[𝒯_0𝒴_0]\}(s,m_\chi ^2,m_h^2,m_A^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2);`$
$``$ Higgs ($`A`$)–$`Z`$ interference term:

$`\stackrel{~}{w}_{hA}^{(AZ)}`$ $`=`$ $`2Re\left[\left({\displaystyle \frac{C^{hAA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}{\displaystyle \frac{C^{hAZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right]`$ (57)
$`\times {\displaystyle \frac{m_\chi (sm_Z^2)(m_A^2m_h^2)}{m_Z^2}};`$
$``$ Higgs ($`A`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{hA}^{(A\chi ^0)}`$ $`=`$ $`2{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \frac{C^{hAA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}C_S^{\chi _i^0\chi h}C_P^{\chi _i^0\chi A}\right]`$ (58)
$`\times \left[sm_{\chi _i^0}+m_\chi (m_A^2m_h^2)\right](s,m_\chi ^2,m_h^2,m_A^2,m_{\chi _i^0}^2);`$
$``$ $`Z`$–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{hA}^{(Z\chi ^0)}`$ $`=`$ $`2{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \frac{C^{hAZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}C_S^{\chi _i^0\chi h}C_P^{\chi _i^0\chi A}\right]{\displaystyle \frac{1}{m_Z^2}}`$ (59)
$`\times [(s+2m_{\chi _i^0}^22m_\chi ^2m_A^2m_h^2)m_Z^2`$
$`+2\{m_{\chi _i^0}s[m_Z^2m_{\chi _i^0}+m_\chi (m_A^2m_h^2)]m_Z^2m_\chi m_{\chi _i^0}(m_A^2m_h^2)`$
$`+m_Z^2m_A^2m_h^2m_Z^2(m_A^2+m_h^2)(m_{\chi _i^0}^2+m_\chi ^2)+m_Z^2(m_{\chi _i^0}^2m_\chi ^2)^2`$
$`+m_\chi ^2(m_A^2m_h^2)^2\}(s,m_\chi ^2,m_h^2,m_A^2,m_{\chi _i^0}^2)].`$
The expressions for $`HA`$ final state are obtained by replacing $`m_h`$, $`C_S^{\chi _i^0\chi h}`$ and $`C^{hAZ}`$ in the above with $`m_H`$, $`C_S^{\chi _i^0\chi H}`$ and $`C^{HAZ}`$, respectively.

$`\mathrm{𝟒}\mathbf{.}𝝌𝝌\mathbf{}𝑯^\mathbf{+}𝑯^{\mathbf{}}`$

This process proceeds via the $`s`$-channel $`Z`$ and CP–even Higgs boson ($`h,H`$) exchange as well as the $`t`$\- and $`u`$-channel chargino ($`\chi _k^\pm `$, $`k=1,2`$) exchange

$`\stackrel{~}{w}_{H^+H^{}}`$ $`=`$ $`\stackrel{~}{w}_{H^+H^{}}^{(h,H)}+\stackrel{~}{w}_{H^+H^{}}^{(Z)}+\stackrel{~}{w}_{H^+H^{}}^{(\chi ^\pm )}+\stackrel{~}{w}_{H^+H^{}}^{(h,H\chi ^\pm )}+\stackrel{~}{w}_{H^+H^{}}^{(Z\chi ^\pm )}:`$ (60)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{w}_{H^+H^{}}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{H^+H^{}r}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2);`$ (61)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{w}_{H^+H^{}}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{3}}\left|{\displaystyle \frac{C^{H^+H^{}Z}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2(s4m_\chi ^2)(s4m_{H^\pm }^2);`$ (62)
$``$ chargino ($`\chi _k^\pm )`$ exchange:

$`\stackrel{~}{w}_{H^+H^{}}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \underset{k,l=1}{\overset{2}{}}}\left[m_{\chi _k^\pm }m_{\chi _l^\pm }I_{kl}^{H^+H^{}}+m_\chi m_{\chi _k^\pm }J_{kl}^{H^+H^{}}+K_{kl}^{H^+H^{}}\right],`$ (63)
where

$`I_{kl}^{H^+H^{}}`$ $`=`$ $`\{(C_k^HC_l^H+D_k^HD_l^H)(2m_\chi ^2)[𝒯_0𝒴_0]`$
$`+(C_k^HC_l^HD_k^HD_l^H)(s2m_\chi ^2)[𝒯_0𝒴_0]\}(s,m_\chi ^2,m_{H^\pm }^2,m_{H^\pm }^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`J_{kl}^{H^+H^{}}`$ $`=`$ $`C_k^HC_{+l}^H[4𝒯_14(m_\chi ^2m_{H^\pm }^2)𝒯_0`$ (65)
$`2𝒴_12(s4m_\chi ^2)𝒴_0](s,m_\chi ^2,m_{H^\pm }^2,m_{H^\pm }^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`K_{kl}^{H^+H^{}}`$ $`=`$ $`\{(C_{+k}^HC_{+l}^H+D_{+k}^HD_{+l}^H)[𝒯_2(s2m_{H^\pm }^2)𝒯_1(m_\chi ^2m_{H^\pm }^2)^2𝒯_0`$
$`+\mathrm{\hspace{0.17em}2}m_\chi ^2(m_\chi ^2m_{H^\pm }^2)𝒴_0]`$
$`+(C_{+k}^HC_{+l}^HD_{+k}^HD_{+l}^H)[2m_\chi ^2𝒯_1𝒴_2`$
$`+(m_\chi ^4m_{H^\pm }^4)𝒴_0]\}(s,m_\chi ^2,m_{H^\pm }^2,m_{H^\pm }^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
and

$`C_{\pm k}^H`$ $``$ $`|C_S^{\chi _k^+\chi H^{}}|^2\pm |C_P^{\chi _k^+\chi H^{}}|^2,`$ (67)
$`D_{\pm k}^H`$ $``$ $`C_S^{\chi _k^+\chi H^{}}\left(C_P^{\chi _k^+\chi H^{}}\right)^{}\pm C_P^{\chi _k^+\chi H^{}}\left(C_S^{\chi _k^+\chi H^{}}\right)^{};`$ (68)
$``$ Higgs ($`h,H`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{H^+H^{}}^{(h,H\chi ^\pm )}`$ $`=`$ $`2{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{H^+H^{}r}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]`$ (69)
$`\times \{C_{+k}^H(2m_\chi )+[C_{+k}^H\mathrm{\hspace{0.17em}2}m_\chi (m_{H^\pm }^2m_\chi ^2m_{\chi _k^\pm }^2)`$
$`+C_k^Hm_{\chi _k^\pm }(s4m_\chi ^2)](s,m_\chi ^2,m_{H^\pm }^2,m_{H^\pm }^2,m_{\chi _k^\pm }^2)\};`$
$``$ $`Z`$–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{H^+H^{}}^{(Z\chi ^\pm )}`$ $`=`$ $`2{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \frac{C^{H^+H^{}Z}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}\right]D_{+k}^H`$ (70)
$`\times \{s2m_{H^\pm }^22m_\chi ^2+2m_{\chi _k^\pm }^2`$
$`+2[sm_{\chi _k^\pm }^2+m_{H^\pm }^42(m_\chi ^2+m_{\chi _k^\pm }^2)m_{H^\pm }^2`$
$`+(m_\chi ^2m_{\chi _k^\pm }^2)^2](s,m_\chi ^2,m_{H^\pm }^2,m_{H^\pm }^2,m_{\chi _k^\pm }^2)\}.`$
$`\mathrm{𝟓}\mathbf{.}𝝌𝝌\mathbf{}𝑾^\mathbf{+}𝑯^{\mathbf{}}`$

This process involves the $`s`$-channel CP–even ($`h,H`$) and odd ($`A`$) Higgs boson exchange as well as the $`t`$\- and $`u`$-channel chargino ($`\chi _k^\pm `$, $`k=1,2`$) exchange

$`\stackrel{~}{w}_{WH}`$ $`=`$ $`\stackrel{~}{w}_{WH}^{(h,H)}+\stackrel{~}{w}_{WH}^{(A)}+\stackrel{~}{w}_{WH}^{(\chi ^\pm )}+\stackrel{~}{w}_{WH}^{(h,H\chi ^\pm )}+\stackrel{~}{w}_{WH}^{(A\chi ^\pm )}:`$ (71)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{w}_{WH}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{W^{}H^+r}C_S^{\chi \chi i}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2)`$ (72)
$`\times {\displaystyle \frac{s^22(m_{H^\pm }^2+m_W^2)s+(m_{H^\pm }^2m_W^2)^2}{m_W^2}};`$
$``$ CP–odd Higgs–boson ($`A`$) exchange:

$`\stackrel{~}{w}_{WH}^{(A)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \frac{C^{W^{}H^+A}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2s`$ (73)
$`\times {\displaystyle \frac{s^22(m_{H^\pm }^2+m_W^2)s+(m_{H^\pm }^2m_W^2)^2}{m_W^2}};`$
$``$ chargino ($`\chi _k^\pm )`$ exchange:

$`\stackrel{~}{w}_{WH}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{m_W^2}}{\displaystyle \underset{k,l=1}{\overset{2}{}}}\left[m_{\chi _k^\pm }m_{\chi _l^\pm }I_{kl}^{WH}+m_\chi m_{\chi _k^\pm }J_{kl}^{WH}+K_{kl}^{WH}\right],`$ (74)
where

$`I_{kl}^{WH}`$ $`=`$ $`\{(C_{+k}^{HW}C_{+l}^{HW}+D_{+k}^{HW}D_{+l}^{HW})`$ (75)
$`\times \left[𝒯_2(sm_W^2m_{H^\pm }^22m_\chi ^2)𝒯_1+G_{WH}^{T(1)}𝒯_0+6m_W^2m_\chi ^2𝒴_0\right]`$
$`+(C_{+k}^{HW}C_{+l}^{HW}D_{+k}^{HW}D_{+l}^{HW})[\mathrm{\hspace{0.17em}6}m_W^2m_\chi ^2𝒯_0𝒴_2`$
$`+G_{WH}^{Y(1)}𝒴_0]\}(s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`J_{kl}^{WH}`$ $`=`$ $`\{Re(C_{+k}^{HW}C_l^{HW}+D_{+k}^{HW}D_l^{HW})[2𝒯_2+2(2m_\chi ^2m_W^2)𝒯_1`$ (76)
$`2(m_\chi ^2m_W^2)(m_\chi ^2+2m_W^2)𝒯_03m_W^2𝒴_1G_{WH}^{Y(2)}𝒴_0]`$
$`+Re(C_{+k}^{HW}C_l^{HW}D_{+k}^{HW}D_l^{HW})[\mathrm{\hspace{0.17em}6}m_W^2𝒯_1+6m_W^2(m_\chi ^2m_{H^\pm }^2)𝒯_0`$
$`2𝒴_2(sm_{H^\pm }^22m_W^2)𝒴_1`$
$`G_{WH}^{Y(3)}𝒴_0]\}(s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`K_{kl}^{WH}`$ $`=`$ $`\{(C_k^{HW}C_l^{HW}+D_k^{HW}D_l^{HW})[(sm_\chi ^22m_W^2)𝒯_2G_{WH}^{T(2)}𝒯_1`$ (77)
$`G_{WH}^{T(3)}𝒯_0+(s2m_W^2)𝒴_2+G_{WH}^{Y(4)}𝒴_0]`$
$`+(C_k^{HW}C_l^{HW}D_k^{HW}D_l^{HW})[\mathrm{\hspace{0.17em}6}m_W^2m_\chi ^2𝒯_1m_\chi ^2𝒴_2`$
$`+G_{WH}^{Y(5)}𝒴_0]\}(s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
and

$`G_{WH}^{T(1)}`$ $`=`$ $`s(m_\chi ^2+2m_W^2)m_\chi ^4m_\chi ^2(m_{H^\pm }^2+3m_W^2)m_W^2m_{H^\pm }^2,`$
$`G_{WH}^{T(2)}`$ $`=`$ $`s(m_\chi ^2+2m_W^2)2m_\chi ^4+m_\chi ^2(m_{H^\pm }^2m_W^2)2m_W^2(m_{H^\pm }^2+m_W^2),`$
$`G_{WH}^{T(3)}`$ $`=`$ $`(m_\chi ^2m_{H^\pm }^2)(m_\chi ^2m_W^2)(m_\chi ^2+2m_W^2),`$
$`G_{WH}^{Y(1)}`$ $`=`$ $`G_{WH}^{T(1)},`$
$`G_{WH}^{Y(2)}`$ $`=`$ $`3sm_W^212m_W^2m_\chi ^2+3m_W^2m_{H^\pm }^23m_W^4,`$
$`G_{WH}^{Y(3)}`$ $`=`$ $`s^2s(2m_\chi ^2+2m_{H^\pm }^2+3m_W^2)+2m_\chi ^4+(2m_\chi ^2+m_{H^\pm }^2)(m_{H^\pm }^2+3m_W^2)2m_W^4,`$
$`G_{WH}^{Y(4)}`$ $`=`$ $`s(m_\chi ^22m_W^2)(m_\chi ^2m_{H^\pm }^2)+4m_W^2m_\chi ^4`$
$`+m_\chi ^2(m_{H^\pm }^45m_W^2m_{H^\pm }^2+2m_W^4)2m_W^4m_{H^\pm }^2,`$
$`G_{WH}^{Y(5)}`$ $`=`$ $`sm_\chi ^2(m_\chi ^2m_W^2)m_\chi ^6m_\chi ^4(m_{H^\pm }^2+3m_W^2)+m_\chi ^2m_W^2(2m_{H^\pm }^2+3m_W^2),`$
and

$`C_{\pm k}^{HW}`$ $``$ $`C_S^{\chi _k^+\chi H^{}}\left(C_V^{\chi _k^+\chi W^{}}\right)^{}\pm C_P^{\chi _k^+\chi H^{}}\left(C_A^{\chi _k^+\chi W^{}}\right)^{},`$ (78)
$`D_{\pm k}^{HW}`$ $``$ $`C_S^{\chi _k^+\chi H^{}}\left(C_A^{\chi _k^+\chi W^{}}\right)^{}\pm C_P^{\chi _k^+\chi H^{}}\left(C_V^{\chi _k^+\chi W^{}}\right)^{};`$ (79)
$``$ Higgs ($`h,H`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{WH}^{(h,H\chi ^\pm )}`$ $`=`$ $`2Re{\displaystyle \underset{k=1}{\overset{2}{}}}{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{W^{}H^+r}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}{\displaystyle \frac{1}{m_W^2}}`$ (80)
$`\times \left[C_{+k}^{HW}H_k^{(+)WH}+C_k^{HW}H_k^{()WH}\right],`$
where

$`H_k^{(+)WH}`$ $`=`$ $`2m_\chi m_{\chi _k^\pm }(s+m_W^2m_{H^\pm }^2)+m_\chi m_{\chi _k^\pm }[s^2+2s(m_{\chi _k^\pm }m_\chi ^2m_{H^\pm }^2)`$
$`+m_{H^\pm }^4m_W^4+2(m_{\chi _k^\pm }^2m_\chi ^2)(m_W^2m_{H^\pm }^2)](s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2),`$
$`H_k^{()WH}`$ $`=`$ $`s^2+s(m_W^2+m_{H^\pm }^2+2m_\chi ^2)+2m_\chi ^2(m_W^2m_{H^\pm }^2)`$ (82)
$`+\{s^2m_{\chi _k^\pm }^2+s[2m_\chi ^2(m_{\chi _k^\pm }^2m_\chi ^2)+m_\chi ^2(m_{H^\pm }^2+3m_W^2)`$
$`+m_{\chi _k^\pm }^2(m_{H^\pm }^2+m_W^2)2m_W^2m_{H^\pm }^2]`$
$`+m_\chi ^2(m_W^2m_{H^\pm }^2)(2m_{\chi _k^\pm }^22m_\chi ^2+m_{H^\pm }^23m_W^2)\}`$
$`\times (s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2);`$
$``$ Higgs ($`A`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{WH}^{(A\chi ^\pm )}`$ $`=`$ $`2Re{\displaystyle \underset{k=1}{\overset{2}{}}}\left({\displaystyle \frac{C^{W^{}H^+A}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}{\displaystyle \frac{1}{m_W^2}}`$ (83)
$`\times \left[D_{+k}^{HW}A_k^{(+)WH}+D_k^{HW}A_k^{()WH}\right],`$
where

$`A_k^{(+)WH}`$ $`=`$ $`m_\chi m_{\chi _k^\pm }[s^22s(m_{H^\pm }^2+m_W^2)`$ (84)
$`+(m_{H^\pm }^2m_W^2)^2](s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2),`$
$`A_k^{()WH}`$ $`=`$ $`s(sm_W^2m_{H^\pm }^2)+\{s^2m_{\chi _k^\pm }^2s[(m_\chi ^2+m_{\chi _k^\pm }^2)(m_{H^\pm }^2+m_W^2)`$ (85)
$`2m_W^2m_{H^\pm }^2]+m_\chi ^2(m_{H^\pm }^2m_W^2)^2\}(s,m_\chi ^2,m_{H^\pm }^2,m_W^2,m_{\chi _k^\pm }^2).`$
Note that $`\stackrel{~}{w}_{WH}`$ in eq.(71) does not include the contribution for $`W^{}H^+`$ final state. The contributions to the $`W^{}H^+`$ and $`W^+H^{}`$ final states are obviously identical and so the total contribution for $`W^\pm H^{}`$ is twice that of $`W^+H^{}`$.

$`\mathrm{𝟔}\mathbf{.}𝝌𝝌\mathbf{}𝒁𝒉`$

This process proceeds via the $`s`$-channel $`Z`$ and CP–odd Higgs boson ($`A`$) exchange as well as the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{Zh}`$ $`=`$ $`\stackrel{~}{w}_{Zh}^{(A)}+\stackrel{~}{w}_{Zh}^{(Z)}+\stackrel{~}{w}_{Zh}^{(\chi ^0)}+\stackrel{~}{w}_{Zh}^{(AZ)}+\stackrel{~}{w}_{Zh}^{(A\chi ^0)}+\stackrel{~}{w}_{Zh}^{(Z\chi ^0)}:`$ (86)
$``$ CP–odd Higgs–boson ($`A`$) exchange:

$`\stackrel{~}{w}_{Zh}^{(A)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \frac{C^{ZhA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2s{\displaystyle \frac{s^22(m_h^2+m_Z^2)s+(m_h^2m_Z^2)^2}{m_Z^2}};`$ (87)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{w}_{Zh}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{12}}\left|{\displaystyle \frac{C^{ZZh}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2{\displaystyle \frac{1}{m_Z^6s}}`$ (88)
$`\times \{m_\chi ^2[6s^412(m_h^2+2m_Z^2)s^3+(32m_Z^4+12m_Z^2m_h^2+6m_h^4)s^2`$
$`(64m_Z^420m_Z^2m_h^2+12m_h^4)sm_Z^2+2(m_h^2m_Z^2)^2m_Z^4]`$
$`+m_Z^4[s^3+(10m_Z^22m_h^2)s^2+(m_h^2m_Z^2)^2s]\};`$
$``$ neutralino ($`\chi _i^0)`$ exchange:

$`\stackrel{~}{w}_{Zh}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{m_Z^2}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}C_S^{\chi _i^0\chi h}C_S^{\chi _j^0\chi h}\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{Zh}+m_\chi m_{\chi _i^0}J_{ij}^{Zh}+K_{ij}^{Zh}\right],`$ (89)
where

$`I_{ij}^{Zh}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )`$ (90)
$`\times \left[𝒯_2(sm_Z^2m_h^22m_\chi ^2)𝒯_1+G_{Zh}^{T(1)}𝒯_0+6m_Z^2m_\chi ^2𝒴_0\right]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[\mathrm{\hspace{0.17em}6}m_Z^2m_\chi ^2𝒯_0𝒴_2`$
$`+G_{Zh}^{Y(1)}𝒴_0]\}(s,m_\chi ^2,m_h^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`J_{ij}^{Zh}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )[2𝒯_2+2(2m_\chi ^2m_Z^2)𝒯_1`$ (91)
$`2(m_\chi ^2m_Z^2)(m_\chi ^2+2m_Z^2)𝒯_03m_Z^2𝒴_1G_{Zh}^{Y(2)}𝒴_0]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[\mathrm{\hspace{0.17em}6}m_Z^2𝒯_1+6m_Z^2(m_\chi ^2m_h^2)𝒯_0`$
$`2𝒴_2(sm_h^22m_Z^2)𝒴_1`$
$`G_{Zh}^{Y(3)}𝒴_0]\}(s,m_\chi ^2,m_h^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`K_{ij}^{Zh}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )[(sm_\chi ^22m_Z^2)𝒯_2G_{Zh}^{T(2)}𝒯_1`$ (92)
$`G_{Zh}^{T(3)}𝒯_0+(s2m_Z^2)𝒴_2+G_{Zh}^{Y(4)}𝒴_0]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[\mathrm{\hspace{0.17em}6}m_Z^2m_\chi ^2𝒯_1m_\chi ^2𝒴_2`$
$`+G_{ZH}^{Y(5)}𝒴_0]\}(s,m_\chi ^2,m_h^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
and

$`G_{Zh}^{T(1)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)m_\chi ^4m_\chi ^2(m_h^2+3m_Z^2)m_Z^2m_h^2,`$
$`G_{Zh}^{T(2)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)2m_\chi ^4+m_\chi ^2(m_h^2m_Z^2)2m_Z^2(m_h^2+m_Z^2),`$
$`G_{Zh}^{T(3)}`$ $`=`$ $`(m_\chi ^2m_h^2)(m_\chi ^2m_Z^2)(m_\chi ^2+2m_Z^2),`$
$`G_{Zh}^{Y(1)}`$ $`=`$ $`G_{Zh}^{T(1)},`$
$`G_{Zh}^{Y(2)}`$ $`=`$ $`3sm_Z^212m_Z^2m_\chi ^2+3m_Z^2m_h^23m_Z^4,`$
$`G_{Zh}^{Y(3)}`$ $`=`$ $`s^2s(2m_\chi ^2+2m_h^2+3m_Z^2)+2m_\chi ^4+(2m_\chi ^2+m_h^2)(m_h^2+3m_Z^2)2m_Z^4,`$
$`G_{Zh}^{Y(4)}`$ $`=`$ $`s(m_\chi ^22m_Z^2)(m_\chi ^2m_h^2)+4m_Z^2m_\chi ^4`$
$`+m_\chi ^2(m_h^45m_Z^2m_h^2+2m_Z^4)2m_Z^4m_h^2,`$
$`G_{Zh}^{Y(5)}`$ $`=`$ $`sm_\chi ^2(m_\chi ^2m_Z^2)m_\chi ^6m_\chi ^4(m_h^2+3m_Z^2)+m_\chi ^2m_Z^2(2m_h^2+3m_Z^2),`$
and

$`C_{Vi}^\chi `$ $``$ $`C_V^{\chi _i^0\chi Z},`$ (93)
$`C_{Ai}^\chi `$ $``$ $`C_A^{\chi _i^0\chi Z};`$ (94)
$``$ Higgs ($`A`$)–$`Z`$ interference term:

$`\stackrel{~}{w}_{Zh}^{(AZ)}`$ $`=`$ $`Re\left[\left({\displaystyle \frac{C^{ZhA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}{\displaystyle \frac{C^{ZZh}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right]`$ (95)
$`\times {\displaystyle \frac{m_\chi (m_Z^2s)}{m_Z^4}}[s^22(m_h^2+m_Z^2)s+(m_h^2m_Z^2)^2];`$
$``$ Higgs ($`A`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{Zh}^{(A\chi ^0)}`$ $`=`$ $`2{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \frac{C^{ZhA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}C_S^{\chi _i^0\chi h}C_A^{\chi _i^0\chi Z}\right]{\displaystyle \frac{1}{m_Z^2}}`$ (96)
$`\times [s^2+(m_h^2+m_Z^2)s\{s^2m_{\chi _i^0}(m_\chi +m_{\chi _i^0})`$
$`+s[2m_h^2m_Z^2(m_h^2+m_Z^2)(m_\chi +m_{\chi _i^0})^2]`$
$`+m_\chi (m_\chi +m_{\chi _i^0})(m_h^2m_Z^2)^2\}(s,m_\chi ^2,m_h^2,m_Z^2,m_{\chi _i^0}^2)];`$
$``$ $`Z`$–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{Zh}^{(Z\chi ^0)}`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \frac{C^{ZZh}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}C_S^{\chi _i^0\chi h}C_A^{\chi _i^0\chi Z}\right]{\displaystyle \frac{1}{m_Z^4}}`$
$`\times \{m_\chi [2s^2+s(2m_h^2+3m_Z^2)+m_Z^2(2m_\chi ^29m_Z^2m_h^22m_{\chi _i^0}^2)]`$
$`+m_{\chi _i^0}m_Z^2[s+(2m_\chi ^2+m_Z^2+m_h^22m_{\chi _i^0}^2)]`$
$`2[m_\chi \{s^2m_{\chi _i^0}^2+s[(m_\chi ^2+m_{\chi _i^0}^2)(m_Z^2+m_h^2)+2m_Z^2m_h^2]`$
$`2m_Z^6+5m_\chi ^2m_Z^4+4m_{\chi _i^0}^2m_Z^4+m_Z^2(m_\chi ^2m_{\chi _i^0}^2)^2`$
$`m_h^2m_Z^2(3m_Z^2+2m_\chi ^2)+m_\chi ^2m_h^4\}`$
$`+m_{\chi _i^0}[s^2m_\chi ^2+s(3m_\chi ^2m_Z^2+m_Z^2m_{\chi _i^0}^22m_\chi ^2m_h^22m_Z^4)`$
$`+10m_\chi ^2m_Z^4m_{\chi _i^0}^2m_Z^4+m_Z^2(m_\chi ^2m_{\chi _i^0}^2)^2`$
$`+m_h^2m_Z^2(m_Z^2m_\chi ^2m_{\chi _i^0}^2)+m_\chi ^2m_h^4]](s,m_\chi ^2,m_h^2,m_Z^2,m_{\chi _i^0}^2)\}.`$
The expressions for $`ZH`$ final state are found by replacing $`h`$ with $`H`$ in the above.

$`\mathrm{𝟕}\mathbf{.}𝝌𝝌\mathbf{}𝒁𝑨`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{ZA}`$ $`=`$ $`\stackrel{~}{w}_{ZA}^{(h,H)}+\stackrel{~}{w}_{ZA}^{(\chi ^0)}+\stackrel{~}{w}_{ZA}^{(h,H\chi ^0)}:`$ (98)
$``$ CP-even Higgs-boson $`(h,H)`$ exchange:

$`\stackrel{~}{w}_{ZA}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{ZAr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2)`$ (99)
$`\times {\displaystyle \frac{s^22(m_A^2+m_Z^2)s+(m_A^2m_Z^2)^2}{m_Z^2}};`$
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{w}_{ZA}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{m_Z^2}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}C_P^{\chi _i^0\chi A}C_P^{\chi _j^0\chi A}\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{ZA}+m_\chi m_{\chi _i^0}J_{ij}^{ZA}+K_{ij}^{ZA}\right],`$ (100)
where

$`I_{ij}^{ZA}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )`$ (101)
$`\times \left[𝒯_2(sm_Z^2m_A^22m_\chi ^2)𝒯_1+G_{ZA}^{T(1)}𝒯_0+6m_Z^2m_\chi ^2𝒴_0\right]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[\mathrm{\hspace{0.17em}6}m_Z^2m_\chi ^2𝒯_0+𝒴_2`$
$`G_{ZA}^{Y(1)}𝒴_0]\}(s,m_\chi ^2,m_A^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`J_{ij}^{ZA}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )[2𝒯_22(2m_\chi ^2m_Z^2)𝒯_1`$ (102)
$`+2(m_\chi ^2m_Z^2)(m_\chi ^2+2m_Z^2)𝒯_0+3m_Z^2𝒴_1+G_{ZA}^{Y(2)}𝒴_0]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[\mathrm{\hspace{0.17em}6}m_Z^2𝒯_1+6m_Z^2(m_\chi ^2m_A^2)𝒯_0`$
$`2𝒴_2(sm_A^22m_Z^2)𝒴_1`$
$`G_{ZA}^{Y(3)}𝒴_0]\}(s,m_\chi ^2,m_A^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`K_{ij}^{ZA}`$ $`=`$ $`\{(C_{Vi}^\chi C_{Vj}^\chi +C_{Ai}^\chi C_{Aj}^\chi )[(sm_\chi ^22m_Z^2)𝒯_2G_{ZA}^{T(2)}𝒯_1`$ (103)
$`G_{ZA}^{T(3)}𝒯_0+(s2m_Z^2)𝒴_2+G_{ZA}^{Y(4)}𝒴_0]`$
$`+(C_{Vi}^\chi C_{Vj}^\chi C_{Ai}^\chi C_{Aj}^\chi )[6m_Z^2m_\chi ^2𝒯_1+m_\chi ^2𝒴_2`$
$`G_{ZA}^{Y(5)}𝒴_0]\}(s,m_\chi ^2,m_A^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
and

$`G_{ZA}^{T(1)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)m_\chi ^4m_\chi ^2(m_A^2+3m_Z^2)m_Z^2m_A^2,`$
$`G_{ZA}^{T(2)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)2m_\chi ^4+m_\chi ^2(m_A^2m_Z^2)2m_Z^2(m_A^2+m_Z^2),`$
$`G_{ZA}^{T(3)}`$ $`=`$ $`(m_\chi ^2m_A^2)(m_\chi ^2m_Z^2)(m_\chi ^2+2m_Z^2),`$
$`G_{ZA}^{Y(1)}`$ $`=`$ $`G_{ZA}^{T(1)},`$
$`G_{ZA}^{Y(2)}`$ $`=`$ $`3sm_Z^212m_Z^2m_\chi ^2+3m_Z^2m_A^23m_Z^4,`$
$`G_{ZA}^{Y(3)}`$ $`=`$ $`s^2s(2m_\chi ^2+2m_A^2+3m_Z^2)+2m_\chi ^4+(2m_\chi ^2+m_A^2)(m_A^2+3m_Z^2)2m_Z^4,`$
$`G_{ZA}^{Y(4)}`$ $`=`$ $`s(m_\chi ^22m_Z^2)(m_\chi ^2m_A^2)+4m_Z^2m_\chi ^4`$
$`+m_\chi ^2(m_A^45m_Z^2m_A^2+2m_Z^4)2m_Z^4m_A^2,`$
$`G_{ZA}^{Y(5)}`$ $`=`$ $`sm_\chi ^2(m_\chi ^2m_Z^2)m_\chi ^6m_\chi ^4(m_A^2+3m_Z^2)+m_\chi ^2m_Z^2(2m_A^2+3m_Z^2).`$
$`G_{ZA}^{T(13)}`$ and $`G_{ZA}^{Y(15)}`$ are obtained from $`G_{Zh}^{T(13)}`$ and $`G_{Zh}^{Y(15)}`$ by replacing $`m_h`$ with $`m_A`$. $`C_{Vi}^\chi `$ and $`C_{Ai}^\chi `$ have already been defined for the $`Zh`$ final state in eqs. (93) and (94).
$``$ Higgs $`(h,H)`$–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{ZA}^{(h,H\chi ^0)}`$ $`=`$ $`2{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{ZAr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}C_A^{\chi _i^0\chi Z}C_P^{\chi _i^0\chi A}\right]{\displaystyle \frac{1}{m_Z^2}}`$ (104)
$`\times [s^2+(m_A^2+m_Z^2)s+2m_\chi (m_{\chi _i^0}m_\chi )(sm_Z^2+m_A^2)`$
$`+\{(m_{\chi _i^0}+m_\chi )(sm_A^2)[m_{\chi _i^0}s+2m_\chi (m_{\chi _i^0}m_\chi )^2m_\chi m_A^2]`$
$`+m_Z^2m_\chi [2(m_{\chi _i^0}m_\chi )^2(m_{\chi _i^0}+m_\chi )+4m_\chi m_A^2`$
$`+m_Z^2(m_{\chi _i^0}3m_\chi )+3sm_\chi ]+sm_Z^2(m_{\chi _i^0}^22m_A^2)\}`$
$`\times (s,m_\chi ^2,m_A^2,m_Z^2,m_{\chi _i^0}^2)].`$
$`\mathrm{𝟖}\mathbf{.}𝝌𝝌\mathbf{}𝑾𝑾`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) and $`Z`$ exchange, and the $`t`$\- and $`u`$-channel chargino ($`\chi _k^\pm `$, $`k=1,2`$) exchange

$`\stackrel{~}{w}_{WW}`$ $`=`$ $`\stackrel{~}{w}_{WW}^{(h,H)}+\stackrel{~}{w}_{WW}^{(Z)}+\stackrel{~}{w}_{WW}^{(\chi ^\pm )}+\stackrel{~}{w}_{WW}^{(h,H\chi ^\pm )}+\stackrel{~}{w}_{WW}^{(Z\chi ^\pm )}:`$ (105)
$``$ CP–even Higgs–boson $`(h,H)`$ exchange:

$`\stackrel{~}{w}_{WW}^{(h,H)}`$ $`=`$ $`\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{WWr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2){\displaystyle \frac{s^24m_W^2s+12m_W^4}{8m_W^4}};`$ (106)
$``$ $`Z`$-boson exchange:

$`\stackrel{~}{w}_{WW}^{(Z)}`$ $`=`$ $`\left|{\displaystyle \frac{C^{WWZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2(s4m_\chi ^2)`$ (107)
$`\times {\displaystyle \frac{s^3+16m_W^2s^268m_W^4s48m_W^6}{12m_W^4}};`$
$``$ chargino ($`\chi _k^\pm `$) exchange:

$`\stackrel{~}{w}_{WW}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{m_W^4}}{\displaystyle \underset{k,l=1}{\overset{2}{}}}\left[m_{\chi _k^\pm }m_{\chi _l^\pm }I_{kl}^{WW}+m_\chi m_{\chi _k^\pm }J_{kl}^{WW}+K_{kl}^{WW}\right],`$ (108)
where

$`I_{kl}^{WW}`$ $`=`$ $`\{(C_k^WC_l^W+D_k^WD_l^W)`$ (109)
$`\times \left[(s4m_W^2)𝒯_2G_{WW}^{T(1)}𝒯_1+G_{WW}^{T(2)}𝒯_0+(s4m_W^2)𝒴_2G_{WW}^{Y(1)}𝒴_0\right]`$
$`+(C_k^WC_l^WD_k^WD_l^W)[18m_W^4m_\chi ^2𝒯_0`$
$`m_\chi ^2(s^2+4m_W^2s+6m_W^4)𝒴_0]\}(s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`J_{kl}^{WW}`$ $`=`$ $`\{C_k^WC_{+l}^W[12m_W^2𝒯_2+12m_W^2(m_W^22m_\chi ^2)𝒯_1`$ (110)
$`+12m_W^2(m_\chi ^4+m_W^2m_\chi ^22m_W^4)𝒯_04(sm_W^2)𝒴_2`$
$`+2m_W^2(2s3m_W^2)𝒴_1+G_{WW}^{Y(2)}𝒴_0]\}(s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
$`K_{kl}^{WW}`$ $`=`$ $`\{(C_{+k}^WC_{+l}^W+D_{+k}^WD_{+l}^W)[𝒯_4(s2m_W^24m_\chi ^2)𝒯_3G_{WW}^{T(3)}𝒯_2`$ (111)
$`G_{WW}^{T(4)}𝒯_1(m_W^2m_\chi ^2)^2(2m_W^2+m_\chi ^2)^2𝒯_0+m_\chi ^2(s4m_W^2)𝒴_2G_{WW}^{Y(3)}𝒴_0]`$
$`+(C_{+k}^WC_{+l}^WD_{+k}^WD_{+l}^W)[\mathrm{\hspace{0.17em}18}m_W^4m_\chi ^2𝒯_1+𝒴_4+(2m_\chi ^44sm_W^2+m_W^4)𝒴_2`$
$`G_{WW}^{Y(4)}𝒴_0]\}(s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2,m_{\chi _l^\pm }^2),`$
and

$`G_{WW}^{T(1)}`$ $`=`$ $`2[s(m_\chi ^2+2m_W^2)4m_W^2(m_\chi ^2+m_W^2)],`$
$`G_{WW}^{T(2)}`$ $`=`$ $`s(m_\chi ^2+2m_W^2)^2+m_W^2(4m_\chi ^410m_\chi ^2m_W^24m_W^4),`$
$`G_{WW}^{T(3)}`$ $`=`$ $`2s(m_\chi ^2+2m_W^2)+6m_\chi ^4+6m_\chi ^2m_W^2+5m_W^4,`$
$`G_{WW}^{T(4)}`$ $`=`$ $`s(m_\chi ^2+2m_W^2)^28m_W^66m_\chi ^4m_W^24m_\chi ^6,`$
$`G_{WW}^{Y(1)}`$ $`=`$ $`sm_\chi ^2(m_\chi ^2+2m_W^2)+6m_\chi ^2m_W^4+4m_W^2(m_\chi ^2m_W^2)^2,`$
$`G_{WW}^{Y(2)}`$ $`=`$ $`2[2m_W^2s^2+s(2m_\chi ^2m_W^22m_\chi ^411m_W^4)`$
$`+2m_W^2(m_\chi ^4+4m_\chi ^2m_W^2+m_W^4)],`$
$`G_{WW}^{Y(3)}`$ $`=`$ $`m_\chi ^2(m_W^2m_\chi ^2)[s(m_\chi ^2+3m_W^2)2m_W^2(2m_\chi ^2+m_W^2)],`$
$`G_{WW}^{Y(4)}`$ $`=`$ $`(m_W^2m_\chi ^2)[s(8m_W^44m_\chi ^2m_W^2)m_\chi ^2(8m_W^4m_\chi ^4m_\chi ^2m_W^2)],`$
and

$`C_{\pm k}^W`$ $``$ $`|C_V^{\chi _k^+\chi W^{}}|^2\pm |C_A^{\chi _k^+\chi W^{}}|^2,`$ (112)
$`D_{\pm k}^W`$ $``$ $`C_V^{\chi _k^+\chi W^{}}\left(C_A^{\chi _k^+\chi W^{}}\right)^{}\pm C_A^{\chi _k^+\chi W^{}}\left(C_V^{\chi _k^+\chi W^{}}\right)^{};`$ (113)
$``$ Higgs ($`h,H`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{WW}^{(h,H\chi ^\pm )}`$ $`=`$ $`{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{WWr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]{\displaystyle \frac{1}{m_W^4}}`$ (114)
$`\times \left[C_{+k}^WH_k^{(+)WW}+C_k^WH_k^{()WW}\right],`$
where

$`H_k^{(+)WW}`$ $`=`$ $`m_\chi [s^2+2s(m_{\chi _k^\pm }^2m_W^2m_\chi ^2)+8m_W^4+4m_W^2m_{\chi _k^\pm }^24m_W^2m_\chi ^2]`$ (115)
$`+2m_\chi \{s^2m_{\chi _k^\pm }^2+s[2m_W^43m_W^2m_\chi ^2m_{\chi _k^\pm }^2m_W^2+(m_\chi ^2m_{\chi _k^\pm }^2)^2]`$
$`4m_W^6+2m_W^4(m_\chi ^2+m_{\chi _k^\pm }^2)`$
$`+2m_W^2(m_\chi ^2m_{\chi _k^\pm }^2)^2\}(s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2),`$
$`H_k^{()WW}`$ $`=`$ $`m_{\chi _k^\pm }(s^22m_W^2s)+m_{\chi _k^\pm }[(m_\chi ^2+m_{\chi _k^\pm }^2)s^2+2s(2m_W^4+m_W^2m_{\chi _k^\pm }^2`$ (116)
$`+3m_W^2m_\chi ^2)24m_\chi ^2m_W^4](s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2);`$
$``$ $`Z`$–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{w}_{WW}^{(Z\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{3m_W^4}}{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \frac{C^{WWZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}\right]D_{+k}^W`$
$`\times \{s^3+s^2(18m_W^2m_\chi ^23m_{\chi _k^\pm }^2)`$
$`+2s[14m_W^418m_\chi ^2m_W^2+6m_{\chi _k^\pm }^2m_W^23(m_\chi ^2m_{\chi _k^\pm }^2)^2]`$
$`+4m_W^2[3(m_\chi ^2m_{\chi _k^\pm }^2)^2m_W^2(11m_\chi ^23m_{\chi _k^\pm }^2)6m_W^4]`$
$`+6[s^2m_{\chi _k^\pm }^2(m_\chi ^2m_{\chi _k^\pm }^2+4m_W^2)+s[8m_W^65m_W^4(m_{\chi _k^\pm }^2+3m_\chi ^2)`$
$`+2m_W^2(3m_\chi ^45m_\chi ^2m_{\chi _k^\pm }^2+2m_{\chi _k^\pm }^4)+(m_\chi ^2m_{\chi _k^\pm }^2)^3]`$
$`+m_W^2[4m_W^62m_W^4(3m_{\chi _k^\pm }^2+5m_\chi ^2)`$
$`+8m_W^2m_\chi ^2(m_\chi ^2m_{\chi _k^\pm }^2)+2(m_{\chi _k^\pm }^2m_\chi ^2)^3]](s,m_\chi ^2,m_W^2,m_W^2,m_{\chi _k^\pm }^2)\}.`$
$`\mathrm{𝟗}\mathbf{.}𝝌𝝌\mathbf{}𝒁𝒁`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$, $`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{w}_{ZZ}`$ $`=`$ $`\stackrel{~}{w}_{ZZ}^{(h,H)}+\stackrel{~}{w}_{ZZ}^{(\chi ^0)}+\stackrel{~}{w}_{ZZ}^{(h,H\chi ^0)}:`$ (118)
$``$ CP–even Higgs–boson $`(h,H)`$ exchange:

$`\stackrel{~}{w}_{ZZ}^{(h,H)}`$ $`=`$ $`\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{ZZr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2){\displaystyle \frac{s^24m_Z^2s+12m_Z^4}{16m_Z^4}};`$ (119)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{w}_{ZZ}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{2m_Z^4}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}\left[m_{\chi _i^0}m_{\chi _j^0}I_{ij}^{ZZ}+m_\chi m_{\chi _i^0}J_{ij}^{ZZ}+K_{ij}^{ZZ}\right],`$ (120)
where

$`I_{ij}^{ZZ}`$ $`=`$ $`\{(C_i^ZC_j^Z+D_i^ZD_j^Z)`$ (121)
$`\times \left[(s4m_Z^2)𝒯_2G_{ZZ}^{T(1)}𝒯_1+G_{ZZ}^{T(2)}𝒯_0+(s4m_Z^2)𝒴_2G_{ZZ}^{Y(1)}𝒴_0\right]`$
$`+(C_i^ZC_j^ZD_i^ZD_j^Z)[18m_Z^4m_\chi ^2𝒯_0`$
$`m_\chi ^2(s^2+4m_Z^2s+6m_Z^4)𝒴_0]\}(s,m_\chi ^2,m_Z^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`J_{ij}^{ZZ}`$ $`=`$ $`C_i^ZC_{+j}^Z[12m_Z^2𝒯_2+12m_Z^2(m_Z^22m_\chi ^2)𝒯_1`$
$`+12m_Z^2(m_\chi ^4+m_\chi ^2m_Z^22m_Z^4)𝒯_04(sm_Z^2)𝒴_2`$
$`+2m_Z^2(2s3m_Z^2)𝒴_1+G_{ZZ}^{Y(2)}𝒴_0](s,m_\chi ^2,m_Z^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
$`K_{ij}^{ZZ}`$ $`=`$ $`C_{+i}^ZC_{+j}^Z[𝒯_4(s2m_Z^24m_\chi ^2)𝒯_3G_{ZZ}^{T(3)}𝒯_2G_{ZZ}^{T(4)}𝒯_1`$ (123)
$`(m_Z^2m_\chi ^2)^2(2m_Z^2+m_\chi ^2)^2𝒯_0+𝒴_4`$
$`+G_{ZZ}^{Y(3)}𝒴_2G_{ZZ}^{Y(4)}𝒴_0](s,m_\chi ^2,m_Z^2,m_Z^2,m_{\chi _i^0}^2,m_{\chi _j^0}^2),`$
and

$`G_{ZZ}^{T(1)}`$ $`=`$ $`2[s(m_\chi ^2+2m_Z^2)4m_Z^2(m_\chi ^2+m_Z^2)],`$
$`G_{ZZ}^{T(2)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)^2+m_Z^2(4m_\chi ^410m_\chi ^2m_Z^24m_Z^4),`$
$`G_{ZZ}^{T(3)}`$ $`=`$ $`2s(m_\chi ^2+2m_Z^2)+6m_\chi ^4+6m_\chi ^2m_Z^2+5m_Z^4,`$
$`G_{ZZ}^{T(4)}`$ $`=`$ $`s(m_\chi ^2+2m_Z^2)^28m_Z^6+18m_Z^4m_\chi ^26m_\chi ^4m_Z^24m_\chi ^6,`$
$`G_{ZZ}^{Y(1)}`$ $`=`$ $`sm_\chi ^2(m_\chi ^2+2m_Z^2)+6m_\chi ^2m_Z^4+4m_Z^2(m_\chi ^2m_Z^2)^2,`$
$`G_{ZZ}^{Y(2)}`$ $`=`$ $`2[2m_Z^2s^2+s(2m_\chi ^2m_Z^22m_\chi ^411m_Z^4)`$
$`+2m_Z^2(m_\chi ^4+4m_\chi ^2m_Z^2+m_Z^4)],`$
$`G_{ZZ}^{Y(3)}`$ $`=`$ $`s(m_\chi ^24m_Z^2)+2m_\chi ^44m_\chi ^2m_Z^2+m_Z^4,`$
$`G_{ZZ}^{Y(4)}`$ $`=`$ $`(m_Z^2m_\chi ^2)[s(8m_Z^4+m_\chi ^4m_\chi ^2m_Z^2)+m_\chi ^2(10m_Z^4+m_\chi ^43m_\chi ^2m_Z^2)],`$
and

$`C_{\pm i}^Z`$ $``$ $`|C_V^{\chi _i^0\chi Z}|^2\pm |C_A^{\chi _i^0\chi Z}|^2,`$ (124)
$`D_{\pm i}^Z`$ $``$ $`C_V^{\chi _i^0\chi Z}\left(C_A^{\chi _i^0\chi Z}\right)^{}\pm C_A^{\chi _i^0\chi Z}\left(C_V^{\chi _i^0\chi Z}\right)^{};`$ (125)
$``$ Higgs $`(h,H)`$–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{w}_{ZZ}^{(h,H\chi ^0)}`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{ZZr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]{\displaystyle \frac{1}{2m_Z^4}}`$ (126)
$`\times \left[C_{+i}^ZH_i^{(+)ZZ}+C_i^ZH_i^{()ZZ}\right],`$
where

$`H_i^{(+)ZZ}`$ $`=`$ $`m_\chi \{s^22s(m_\chi ^2+m_Z^2m_{\chi _i^0}^2)4m_Z^2(m_\chi ^22m_Z^2m_{\chi _i^0}^2)`$ (127)
$`+2[s^2m_{\chi _i^0}^2+s\{(m_\chi ^2m_{\chi _i^0}^2)^2m_Z^2(3m_\chi ^22m_Z^2+m_{\chi _i^0}^2)\}`$
$`+2m_Z^2\{(m_\chi ^2m_{\chi _i^0}^2)^2+m_Z^2(m_\chi ^22m_Z^2+m_{\chi _i^0}^2)\}]`$
$`\times (s,m_\chi ^2,m_Z^2,m_Z^2,m_{\chi _i^0}^2)\},`$
$`H_i^{()ZZ}`$ $`=`$ $`m_{\chi _i^0}\{s^22m_Z^2s+[s^2(m_\chi ^2+m_{\chi _i^0}^2)2m_Z^2(2m_Z^2+3m_\chi ^2+m_{\chi _i^0}^2)s`$ (128)
$`+24m_Z^4m_\chi ^2](s,m_\chi ^2,m_Z^2,m_Z^2,m_{\chi _i^0}^2)\}.`$
Note that $`D_{+i}^Z`$ $`=`$ 0, because $`C_V^{\chi _i^0\chi _j^0Z}`$ is pure imaginary and $`C_A^{\chi _i^0\chi _j^0Z}`$ is real.

$`\mathrm{𝟏𝟎}\mathbf{.}𝝌𝝌\mathbf{}\overline{𝒇}𝒇`$

This process involves the $`s`$-channel Higgs boson ($`h`$, $`H`$ and $`A`$) and $`Z`$ boson exchange and the $`t`$\- and $`u`$-channel sfermion ($`\stackrel{~}{f}_a`$) exchange

$`\stackrel{~}{w}_{\overline{f}f}`$ $`=`$ $`\stackrel{~}{w}_{\overline{f}f}^{(h,H)}+\stackrel{~}{w}_{\overline{f}f}^{(A)}+\stackrel{~}{w}_{\overline{f}f}^{(Z)}+\stackrel{~}{w}_{\overline{f}f}^{(\stackrel{~}{f})}+\stackrel{~}{w}_{\overline{f}f}^{(h,H\stackrel{~}{f})}+\stackrel{~}{w}_{\overline{f}f}^{(AZ)}+\stackrel{~}{w}_{\overline{f}f}^{(A\stackrel{~}{f})}+\stackrel{~}{w}_{\overline{f}f}^{(Z\stackrel{~}{f})}:`$
$``$ CP–even Higgs–boson $`(h,H)`$ exchange:

$`\stackrel{~}{w}_{\overline{f}f}^{(h,H)}`$ $`=`$ $`\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C_S^{ffr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(s4m_\chi ^2)(s4m_f^2);`$ (130)
$``$ CP–odd Higgs–boson ($`A`$) exchange:

$`\stackrel{~}{w}_{\overline{f}f}^{(A)}`$ $`=`$ $`\left|{\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2s^2;`$ (131)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{w}_{\overline{f}f}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{4}{3}}\left|{\displaystyle \frac{C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2`$ (132)
$`\times \{12|C_A^{ffZ}|^2{\displaystyle \frac{m_\chi ^2m_f^2}{m_Z^4}}(sm_Z^2)^2`$
$`+[|C_V^{ffZ}|^2(s+2m_f^2)+|C_A^{ffZ}|^2(s4m_f^2)](s4m_\chi ^2)\};`$
$``$ sfermion ($`\stackrel{~}{f}_a`$) exchange:

$`\stackrel{~}{w}_{\overline{f}f}^{(\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{4}}{\displaystyle \underset{a,b}{}}[(C_+^aC_+^b+D_+^aD_+^b)[𝒯_22(m_\chi ^2+m_f^2)𝒯_1`$ (133)
$`+\{(m_\chi ^2+m_f^2)^2+4m_\chi ^2m_f^2\}𝒯_0]`$
$`+(C_{}^aC_{}^bD_{}^aD_{}^b)\left[𝒯_22(m_\chi ^2+m_f^2)𝒯_1+(m_\chi ^2m_f^2)^2𝒯_0\right]`$
$`(C_+^aC_+^bD_{}^aD_{}^b)\mathrm{\hspace{0.17em}4}m_\chi ^2m_f^2𝒴_0`$
$`+(C_+^aC_+^b+C_{}^aC_{}^b)sm_\chi ^2𝒴_0+(C_+^aC_+^bC_{}^aC_{}^b)sm_f^2𝒴_0`$
$`+(D_+^aD_+^bD_{}^aD_{}^b)\left[𝒴_2+(m_\chi ^2+m_f^2)^2𝒴_0\right]`$
$`+(C_+^aD_+^b+D_+^aC_+^b)m_\chi m_f\left[4(m_\chi ^2+m_f^2)𝒯_04𝒯_1+s𝒴_0\right]`$
$`(C_+^aD_+^bD_+^aC_+^b)m_\chi m_f𝒴_1](s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2,m_{\stackrel{~}{f}_b}^2),`$
where

$`C_\pm ^a`$ $`=`$ $`|\mathrm{\Lambda }_{fL}^a|^2\pm |\mathrm{\Lambda }_{fR}^a|^2,`$ (134)
$`D_\pm ^a`$ $`=`$ $`\mathrm{\Lambda }_{fL}^a(\mathrm{\Lambda }_{fR}^a)^{}\pm (\mathrm{\Lambda }_{fL}^a)^{}\mathrm{\Lambda }_{fR}^a.`$ (135)
In the above, $`a`$ is the index for sfermion mass eigenstates so that $`a`$ $`=`$ 1, $`\mathrm{}`$, 6 for squarks and charged sleptons, and $`a`$ $`=`$ 1, 2, 3 for sneutrinos.

The symbol $`f`$ represents each fermion: $`f`$ $`=`$ $`u`$, $`c`$, $`t`$, $`\mathrm{}`$. The symbol $`\stackrel{~}{f}_a`$ should be understood as follows. For up–type quarks, down–type quarks and charged leptons, the corresponding symbol $`\stackrel{~}{f}_a`$ represents $`\stackrel{~}{u}_a`$, $`\stackrel{~}{d}_a`$ and $`\stackrel{~}{e}_a`$ ($`a`$ $`=`$ 1, $`\mathrm{}`$, 6), respectively. For neutrinos, $`\stackrel{~}{f}_a`$ represents $`\stackrel{~}{\nu }_a`$ ($`a`$ $`=`$ 1, 2, 3). For eg. for $`\stackrel{~}{w}_{\overline{c}c}^{(\stackrel{~}{u})}`$, the last argument at the end of eq. (133) should read $`(s,m_\chi ^2,m_c^2,m_c^2,m_{\stackrel{~}{u}_a}^2,m_{\stackrel{~}{u}_b}^2)`$ when $`a`$ $`=`$ 1, $`\mathrm{}`$, 6.

The coupling $`\mathrm{\Lambda }_{fL}^a`$ for each fermion is defined by $`\mathrm{\Lambda }_{nai}^{(f)L}`$ in Appendix A as $`\mathrm{\Lambda }_{uL}^a`$ $`=`$ $`\mathrm{\Lambda }_{1a1}^{(u)L}`$, $`\mathrm{\Lambda }_{cL}^a`$ $`=`$ $`\mathrm{\Lambda }_{2a1}^{(u)L}`$, $`\mathrm{\Lambda }_{tL}^a`$ $`=`$ $`\mathrm{\Lambda }_{3a1}^{(u)L}`$, $`\mathrm{\Lambda }_{dL}^a`$ $`=`$ $`\mathrm{\Lambda }_{1a1}^{(d)L}`$, $`\mathrm{\Lambda }_{sL}^a`$ $`=`$ $`\mathrm{\Lambda }_{2a1}^{(d)L}`$, $`\mathrm{\Lambda }_{bL}^a`$ $`=`$ $`\mathrm{\Lambda }_{3a1}^{(d)L}`$, etc. Similarly, the coupling $`\mathrm{\Lambda }_{fR}^a`$ is defined by $`\mathrm{\Lambda }_{nai}^{(f)R}`$.
$``$ Higgs $`(h,H)`$–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{w}_{\overline{f}f}^{(h,H\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \underset{a}{}}Re\left[{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C_S^{ffr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right]`$ (136)
$`\times [C_+^a2m_\chi m_f\{2+[s2(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)]`$
$`\times (s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\}`$
$`+D_+^a\{s+[s(m_\chi ^2+m_f^2+m_{\stackrel{~}{f}_a}^2)8m_\chi ^2m_f^2]`$
$`\times (s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\}];`$
$``$ Higgs ($`A`$)–$`Z`$ interference term:

$`\stackrel{~}{w}_{\overline{f}f}^{(AZ)}`$ $`=`$ $`Re\left[\left({\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}\left({\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)\right]`$ (137)
$`\times {\displaystyle \frac{8m_\chi m_f}{m_Z^2}}s(m_Z^2s);`$
$``$ Higgs ($`A`$)–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{w}_{\overline{f}f}^{(A\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \underset{a}{}}Re\left[{\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{sm_A^2+i\mathrm{\Gamma }_Am_A}}\right]`$ (138)
$`\times [C_+^a2m_\chi m_fs(s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)`$
$`+D_+^a\{s+s(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)`$
$`\times (s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\}];`$
$``$ $`Z`$–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{w}_{\overline{f}f}^{(Z\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \underset{a}{}}\left\{Re\left[{\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right]Z_a^{(+)\overline{f}f}+Re\left[{\displaystyle \frac{C_V^{ffZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right]Z_a^{()\overline{f}f}\right\},`$
where

$`Z_a^{(+)\overline{f}f}`$ $`=`$ $`C_+^a\{s+2(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)[2s(m_\chi ^2+m_f^2)+2(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)^2`$
$`+8m_\chi ^2m_f^2(2{\displaystyle \frac{s}{m_Z^2}})](s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\}`$
$`D_+^a4m_\chi m_f\{{\displaystyle \frac{s}{m_Z^2}}3`$
$`[s+({\displaystyle \frac{s}{m_Z^2}}3)(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)](s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\},`$
$`Z_a^{()\overline{f}f}`$ $`=`$ $`C_{}^a\{s+2(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)`$ (141)
$`+[2s(m_\chi ^2m_f^2)2(m_\chi ^2+m_f^2m_{\stackrel{~}{f}_a}^2)^2+8m_\chi ^2m_f^2]`$
$`\times (s,m_\chi ^2,m_f^2,m_f^2,m_{\stackrel{~}{f}_a}^2)\}.`$
This completes the list of all the tree-level two-body neutralino pair-annihilation channels in the MSSM.

## 5 Partial Wave Expansion

In the literature one still often uses the usual approximation in terms of the expansion in powers of $`x`$ (or, equivalently, WIMP velocity-square), $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}a+bx`$. This is because in the early days it was often much easier to compute the coefficients $`a`$ and $`b`$, rather than the cross section itself . Furthermore, the partial wave expansion gives a rather good approximation to the exact result but only far enough from $`s`$–channel resonances and thresholds for new final states. (For a recent detailed study, see Ref. .)

The expansion of $`\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}`$, as defined in eq. (2), is given by

$$\sigma v_{\mathrm{M}\mathrm{ø}\mathrm{l}}=\frac{1}{m_\chi ^2}\left[w\frac{3}{2}\left(2ww^{}\right)x+𝒪(x^2)\right]_{s=4m_\chi ^2}a+bx+𝒪(x^2)$$
(142)
where $`w^{}(s)`$ denotes $`dw(s)/d(s/4m_\chi ^2)`$.

Analogously to the function $`w(s)`$ (eq. (30)), the coefficients $`a`$ and $`b`$ need to be summed over all possible final states $`f_1f_2`$. One can write them as

$`a`$ $`=`$ $`{\displaystyle \underset{f_1f_2}{}}c\theta \left(4m_\chi ^2(m_{f_1}+m_{f_2})^2\right)v_{f_1f_2}\stackrel{~}{a}_{f_1f_2},`$ (143)
$`b`$ $`=`$ $`{\displaystyle \underset{f_1f_2}{}}c\theta (4m_\chi ^2(m_{f_1}+m_{f_2})^2)v_{f_1f_2}\{\stackrel{~}{b}_{f_1f_2}+`$
$`\stackrel{~}{a}_{f_1f_2}[3+{\displaystyle \frac{3}{4}}v_{f_1f_2}^2({\displaystyle \frac{m_{f_1}^2+m_{f_2}^2}{2m_\chi ^2}}+{\displaystyle \frac{(m_{f_1}^2m_{f_2}^2)^2}{8m_\chi ^4}})]\},`$
where the summation extends over all possible two-body final states $`f_1f_2`$, the coefficient $`c`$ is defined in eq. (33), and

$`v_{f_1f_2}`$ $`=`$ $`\beta _f(4m_\chi ^2,m_{f_1},m_{f_2}),`$ (145)
and the velocity $`\beta _f(s,m_{f_1},m_{f_2})`$ was defined in eq. (34). The “reduced” coefficients $`\stackrel{~}{a}_{f_1f_2}`$ and $`\stackrel{~}{b}_{f_1f_2}`$ are given by

$`\stackrel{~}{a}_{f_1f_2}`$ $`=`$ $`{\displaystyle \frac{1}{32\pi m_\chi ^2}}\stackrel{~}{w}_{f_1f_2}(4m_\chi ^2),`$ (146)
$`\stackrel{~}{b}_{f_1f_2}`$ $`=`$ $`{\displaystyle \frac{3}{64\pi m_\chi ^2}}\stackrel{~}{w}_{f_1f_2}^{}(4m_\chi ^2),`$ (147)
where $`\stackrel{~}{w}(s)`$ was defined in eq. (35) and $`\stackrel{~}{w}^{}(s)d\stackrel{~}{w}(s)/d(s/4m_\chi ^2)`$.

In this Section, we provide a set of expressions for the coefficients $`a`$ and $`b`$ in the case of equal-mass final states. Using eq. (142), we have derived the coefficients for all the final states by using the analytic expressions for $`w(s)`$ presented in the previous Section. In the case of unequal masses of the final state particles the resulting formulae are exceedingly lengthly and we will not include them here.

In the literature one can find several analytic formulae for the expansion coefficients, including , but, due to different conventions and complexity of expressions, comparison is not always doable. We have checked our results for the $`a`$–coefficients in appropriate limits against published results and agreed in all cases.

All the couplings and auxiliary functions appearing below are listed in the Appendices.

$`\mathrm{𝟏}\mathbf{.}𝝌𝝌\mathbf{}𝒉𝒉`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$,$`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{a}_{hh}`$ $`=`$ $`\stackrel{~}{a}_{hh}^{(h,H)}+\stackrel{~}{a}_{hh}^{(\chi ^0)}+\stackrel{~}{a}_{hh}^{(h,H\chi ^0)},`$ (148)
$`\stackrel{~}{b}_{hh}`$ $`=`$ $`\stackrel{~}{b}_{hh}^{(h,H)}+\stackrel{~}{b}_{hh}^{(\chi ^0)}+\stackrel{~}{b}_{hh}^{(h,H\chi ^0)}:`$ (149)
$``$ CP-even Higgs-boson ($`h,H`$) exchange:

$`\stackrel{~}{a}_{hh}^{(h,H)}`$ $`=`$ $`0,`$ (150)
$`\stackrel{~}{b}_{hh}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{64\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{hhr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2;`$ (151)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{a}_{hh}^{(\chi ^0)}`$ $`=`$ $`0,`$ (152)
$`\stackrel{~}{b}_{hh}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i,j=1}{\overset{4}{}}}(C_S^{\chi _i^0\chi h})^2(C_S^{\chi _j^0\chi h})^2{\displaystyle \frac{1}{\mathrm{\Delta }_{hi}^2\mathrm{\Delta }_{hj}^2}}`$ (153)
$`\times [4m_\chi ^2(m_\chi ^2m_h^2)^2+4m_\chi (m_\chi ^2m_h^2)(m_\chi +m_{\chi _i^0})\mathrm{\Delta }_{hi}`$
$`+3(m_\chi +m_{\chi _i^0})(m_\chi +m_{\chi _j^0})\mathrm{\Delta }_{hi}\mathrm{\Delta }_{hj}],`$
where $`\mathrm{\Delta }_{hi}m_h^2m_\chi ^2m_{\chi _i^0}^2`$.
$``$ Higgs ($`h,H`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{a}_{hh}^{(h,H\chi ^0)}`$ $`=`$ $`0,`$ (154)
$`\stackrel{~}{b}_{hh}^{(h,H\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{hhr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}C_S^{\chi _i^0\chi h}C_S^{\chi _i^0\chi h}\right]`$ (155)
$`\times {\displaystyle \frac{[2m_\chi (m_\chi ^2m_h^2)+3(m_\chi +m_{\chi _i^0})\mathrm{\Delta }_{hi}]}{\mathrm{\Delta }_{hi}^2}}.`$
The expressions for $`HH`$ final state are obtained from the above by replacing $`m_h`$, $`C^{hhr}`$, $`C_S^{\chi _i^0\chi h}`$, $`C_S^{\chi _j^0\chi h}`$ with $`m_H`$, $`C^{HHr}`$, $`C_S^{\chi _i^0\chi H}`$, $`C_S^{\chi _j^0\chi H}`$, respectively.

$`\mathrm{𝟐}\mathbf{.}𝝌𝝌\mathbf{}𝑨𝑨`$

Similarly to the final state $`hh`$, this process proceeds via the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$,$`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{a}_{AA}`$ $`=`$ $`\stackrel{~}{a}_{AA}^{(h,H)}+\stackrel{~}{a}_{AA}^{(\chi ^0)}+\stackrel{~}{a}_{AA}^{(h,H\chi ^0)},`$ (156)
$`\stackrel{~}{b}_{AA}`$ $`=`$ $`\stackrel{~}{b}_{AA}^{(h,H)}+\stackrel{~}{b}_{AA}^{(\chi ^0)}+\stackrel{~}{b}_{AA}^{(h,H\chi ^0)}:`$ (157)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{a}_{AA}^{(h,H)}`$ $`=`$ $`0,`$ (158)
$`\stackrel{~}{b}_{AA}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{64\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{AAr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2;`$ (159)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{a}_{AA}^{(\chi ^0)}`$ $`=`$ $`0,`$ (160)
$`\stackrel{~}{b}_{AA}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i,j=1}{\overset{4}{}}}(C_P^{\chi _i^0\chi A})^2(C_P^{\chi _j^0\chi A})^2{\displaystyle \frac{1}{\mathrm{\Delta }_{Ai}^2\mathrm{\Delta }_{Aj}^2}}`$ (161)
$`\times [4m_\chi ^2(m_\chi ^2m_A^2)^2+4m_\chi (m_\chi ^2m_A^2)(m_\chi m_{\chi _i^0})\mathrm{\Delta }_{Ai}`$
$`+3(m_\chi m_{\chi _i^0})(m_\chi m_{\chi _j^0})\mathrm{\Delta }_{Ai}\mathrm{\Delta }_{Aj}],`$
where $`\mathrm{\Delta }_{Ai}m_A^2m_\chi ^2m_{\chi _i^0}^2`$.
$``$ Higgs ($`h,H`$)–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{a}_{AA}^{(h,H\chi ^0)}`$ $`=`$ $`0,`$ (162)
$`\stackrel{~}{b}_{AA}^{(h,H\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{AAr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}(C_P^{\chi _i^0\chi A})^2\right]`$ (163)
$`\times {\displaystyle \frac{[2m_\chi (m_\chi ^2m_A^2)3(m_\chi m_{\chi _i^0})\mathrm{\Delta }_{Ai}]}{\mathrm{\Delta }_{Ai}^2}}.`$
$`\mathrm{𝟑}\mathbf{.}𝝌𝝌\mathbf{}𝑯^\mathbf{+}𝑯^{\mathbf{}}`$

This process proceeds via the $`s`$-channel $`Z`$ and CP–even Higgs boson ($`h,H`$) exchange as well as the $`t`$\- and $`u`$-channel chargino ($`\chi _k^\pm `$, $`k=1,2`$) exchange

$`\stackrel{~}{a}_{H^+H^{}}`$ $`=`$ $`\stackrel{~}{a}_{H^+H^{}}^{(h,H)}+\stackrel{~}{a}_{H^+H^{}}^{(Z)}+\stackrel{~}{a}_{H^+H^{}}^{(\chi ^\pm )}+\stackrel{~}{a}_{H^+H^{}}^{(h,H\chi ^\pm )}+\stackrel{~}{a}_{H^+H^{}}^{(Z\chi ^\pm )},`$ (164)
$`\stackrel{~}{b}_{H^+H^{}}`$ $`=`$ $`\stackrel{~}{b}_{H^+H^{}}^{(h,H)}+\stackrel{~}{b}_{H^+H^{}}^{(Z)}+\stackrel{~}{b}_{H^+H^{}}^{(\chi ^\pm )}+\stackrel{~}{b}_{H^+H^{}}^{(h,H\chi ^\pm )}+\stackrel{~}{b}_{H^+H^{}}^{(Z\chi ^\pm )}:`$ (165)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{a}_{H^+H^{}}^{(h,H)}`$ $`=`$ $`0,`$ (166)
$`\stackrel{~}{b}_{H^+H^{}}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{32\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{H^+H^{}r}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2;`$ (167)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{a}_{H^+H^{}}^{(Z)}`$ $`=`$ $`0,`$ (168)
$`\stackrel{~}{b}_{H^+H^{}}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}\left|{\displaystyle \frac{C^{H^+H^{}Z}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2(m_\chi ^2m_{H^\pm }^2);`$ (169)
$``$ chargino ($`\chi _k^\pm )`$ exchange:

$`\stackrel{~}{a}_{H^+H^{}}^{(\chi ^\pm )}`$ $`=`$ $`0,`$ (170)
$`\stackrel{~}{b}_{H^+H^{}}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{k,l=1}{\overset{2}{}}}{\displaystyle \frac{1}{\mathrm{\Delta }_{H^\pm k}^2\mathrm{\Delta }_{H^\pm l}^2}}\{C_{+k}^HC_{+l}^Hm_\chi ^2[4(m_\chi ^2m_{H^\pm }^2)^2`$ (171)
$`+4(m_\chi ^2m_{H^\pm }^2)\mathrm{\Delta }_{H^\pm k}+3\mathrm{\Delta }_{H^\pm k}\mathrm{\Delta }_{H^\pm l}]`$
$`+D_k^HD_l^H\left[4m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(m_\chi ^2m_{H^\pm }^2)\right]`$
$`+C_{+k}^HC_l^H\mathrm{\hspace{0.17em}2}m_\chi m_{\chi _k^\pm }\left[2(m_\chi ^2m_{H^\pm }^2)\mathrm{\Delta }_{H^\pm k}+3\mathrm{\Delta }_{H^\pm k}^2\right]`$
$`+C_k^HC_l^H\left[\mathrm{\hspace{0.17em}3}m_{\chi _k^\pm }m_{\chi _l^\pm }\mathrm{\Delta }_{H^\pm k}\mathrm{\Delta }_{H^\pm l}\right]`$
$`+D_{+k}^HD_{+l}^H\left[2(m_\chi ^2m_{H^\pm }^2)\mathrm{\Delta }_{H^\pm k}\mathrm{\Delta }_{H^\pm l}\right]\},`$
where $`\mathrm{\Delta }_{H^\pm k}m_{H^\pm }^2m_\chi ^2m_{\chi _k^\pm }^2`$. The couplings $`C_{\pm k}^H`$ and $`D_{\pm k}^H`$ are given in eqs. (67) and (68).
$``$ Higgs ($`h,H`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{a}_{H^+H^{}}^{(h,H\chi ^\pm )}`$ $`=`$ $`0,`$ (172)
$`\stackrel{~}{b}_{H^+H^{}}^{(h,H\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[{\displaystyle \underset{r=h,H}{}}\left({\displaystyle \frac{C^{H^+H^{}r}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]{\displaystyle \frac{1}{\mathrm{\Delta }_{H^\pm }^2}}`$ (173)
$`\times \left\{C_{+k}^Hm_\chi \left[2(m_\chi ^2m_{H^\pm }^2)+3\mathrm{\Delta }_{H_k^\pm }\right]+C_k^H3m_{\chi _k^\pm }\mathrm{\Delta }_{H_k^\pm }\right\};`$
$``$ $`Z`$–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{a}_{H^+H^{}}^{(Z\chi ^\pm )}`$ $`=`$ $`0,`$ (174)
$`\stackrel{~}{b}_{H^+H^{}}^{(Z\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \frac{C^{H^+H^{}Z}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}D_{+k}^H\right]{\displaystyle \frac{(m_\chi ^2m_{H^\pm }^2)}{\mathrm{\Delta }_{H^\pm k}}}.`$ (175)
$`\mathrm{𝟒}\mathbf{.}𝝌𝝌\mathbf{}𝑾^\mathbf{+}𝑾^{\mathbf{}}`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) and $`Z`$ exchange, and the $`t`$\- and $`u`$-channel chargino ($`\chi _k^\pm `$, $`k=1,2`$) exchange

$`\stackrel{~}{a}_{WW}`$ $`=`$ $`\stackrel{~}{a}_{WW}^{(h,H)}+\stackrel{~}{a}_{WW}^{(Z)}+\stackrel{~}{a}_{WW}^{(\chi ^\pm )}+\stackrel{~}{a}_{WW}^{(h,H\chi ^\pm )}+\stackrel{~}{a}_{WW}^{(Z\chi ^\pm )},`$ (176)
$`\stackrel{~}{b}_{WW}`$ $`=`$ $`\stackrel{~}{b}_{WW}^{(h,H)}+\stackrel{~}{b}_{WW}^{(Z)}+\stackrel{~}{b}_{WW}^{(\chi ^\pm )}+\stackrel{~}{b}_{WW}^{(h,H\chi ^\pm )}+\stackrel{~}{b}_{WW}^{(Z\chi ^\pm )};`$ (177)
$``$ CP–even Higgs–boson ($`h,H`$) exchange:

$`\stackrel{~}{a}_{WW}^{(h,H)}`$ $`=`$ $`0,`$ (178)
$`\stackrel{~}{b}_{WW}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{32\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{WWr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2{\displaystyle \frac{3m_W^44m_W^2m_\chi ^2+4m_\chi ^4}{m_W^4}};`$ (179)
$``$ $`Z`$-boson exchange:

$`\stackrel{~}{a}_{WW}^{(Z)}`$ $`=`$ $`0,`$ (180)
$`\stackrel{~}{b}_{WW}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}\left|{\displaystyle \frac{C^{WWZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2{\displaystyle \frac{(3m_W^4+20m_W^2m_\chi ^2+4m_\chi ^4)(m_\chi ^2m_W^2)}{m_W^4}};`$ (181)
$``$ chargino ($`\chi _k^\pm `$) exchange:

$`\stackrel{~}{a}_{WW}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \underset{k,l=1}{\overset{2}{}}}{\displaystyle \frac{1}{m_W^4\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}}}\{C_{+k}^WC_{+l}^W\left[2m_W^4(m_\chi ^2m_W^2)\right]`$
$`+D_{+k}^WD_{+l}^W\left[m_{\chi _k^\pm }m_{\chi _l^\pm }(4m_\chi ^4+3m_W^44m_W^2m_\chi ^2)\right]\},`$
$`\stackrel{~}{b}_{WW}^{(\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{k,l=1}{\overset{2}{}}}{\displaystyle \frac{1}{m_W^4\mathrm{\Delta }_{Wk}^3\mathrm{\Delta }_{Wl}^3}}\{C_{+k}^WC_{+l}^W[16m_\chi ^2m_W^4(m_\chi ^2m_W^2)^2\mathrm{\Delta }_{Wk}^2`$ (183)
$`+4m_\chi ^2m_W^2(3m_W^2+4m_\chi ^2)(m_\chi ^2m_W^2)^2\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}`$
$`4m_\chi ^2m_W^2(7m_W^4+m_\chi ^2m_W^28m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}`$
$`+(8m_W^6+3m_W^4m_\chi ^2+8m_W^2m_\chi ^4+8m_\chi ^6)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}^2]`$
$`+C_{+k}^WC_l^W[48m_W^2m_\chi ^3m_{\chi _k^\pm }(m_\chi ^2m_W^2)^2\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}`$
$`+4m_W^2m_\chi m_{\chi _k^\pm }(3m_W^45m_W^2m_\chi ^2+2m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}`$
$`+8m_W^2m_\chi m_{\chi _k^\pm }(3m_W^48m_W^2m_\chi ^2+5m_\chi ^4)\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}^2`$
$`+2m_\chi m_{\chi _k^\pm }(5m_W^4+10m_W^2m_\chi ^28m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}^2]`$
$`+C_k^WC_l^W[16m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(m_\chi ^2m_W^2)^2(m_W^2+m_\chi ^2)\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}`$
$`+8m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(m_W^43m_W^2m_\chi ^2+2m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}`$
$`+3m_{\chi _k^\pm }m_{\chi _l^\pm }(3m_W^44m_W^2m_\chi ^2+4m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}^2]`$
$`+D_{+k}^WD_{+l}^W[16m_\chi ^2m_W^2(m_\chi ^2m_W^2)^3\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}`$
$`+32m_\chi ^2m_W^2(m_\chi ^2m_W^2)^2\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}`$
$`+2(3m_W^617m_W^4m_\chi ^2+16m_W^2m_\chi ^4+4m_\chi ^6)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}^2]`$
$`+D_k^WD_l^W[8m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(3m_W^6+7m_W^4m_\chi ^28m_W^2m_\chi ^4+4m_\chi ^6)\mathrm{\Delta }_{Wk}^2`$
$`8m_W^2m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(m_W^4+m_W^2m_\chi ^22m_\chi ^4)\mathrm{\Delta }_{Wk}\mathrm{\Delta }_{Wl}`$
$`+4m_\chi ^2m_{\chi _k^\pm }m_{\chi _l^\pm }(11m_W^418m_W^2m_\chi ^2+16m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}`$
$`+3m_{\chi _k^\pm }m_{\chi _l^\pm }(3m_W^48m_W^2m_\chi ^2+12m_\chi ^4)\mathrm{\Delta }_{Wk}^2\mathrm{\Delta }_{Wl}^2]\},`$
where $`\mathrm{\Delta }_{Wk}m_W^2m_\chi ^2m_{\chi _k^\pm }^2`$. The couplings $`C_{\pm k}^W`$ and $`D_{\pm k}^W`$ are given in eqs. (112) and (113).
$``$ Higgs ($`h,H`$)–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{a}_{WW}^{(h,H\chi ^\pm )}`$ $`=`$ $`0,`$ (184)
$`\stackrel{~}{b}_{WW}^{(h,H\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{WWr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]{\displaystyle \frac{1}{m_W^4\mathrm{\Delta }_{Wk}^2}}`$ (185)
$`\times \{C_{+k}^W[2m_\chi m_W^2(3m_W^45m_W^2m_\chi ^2+2m_\chi ^4)`$
$`+m_\chi (5m_W^4+10m_W^2m_\chi ^28m_\chi ^4)\mathrm{\Delta }_{Wk}]`$
$`+C_k^W[4m_\chi ^2m_{\chi _k^\pm }(m_W^43m_W^2m_\chi ^2+2m_\chi ^4)`$
$`+3m_{\chi _k^\pm }(3m_W^44m_W^2m_\chi ^2+4m_\chi ^4)\mathrm{\Delta }_{Wk}]\};`$
$``$ $`Z`$–chargino ($`\chi _k^\pm `$) interference term:

$`\stackrel{~}{a}_{WW}^{(Z\chi ^\pm )}`$ $`=`$ $`0,`$ (186)
$`\stackrel{~}{b}_{WW}^{(Z\chi ^\pm )}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle \underset{k=1}{\overset{2}{}}}Re\left[\left({\displaystyle \frac{C^{WWZ}C_A^{\chi \chi Z}}{sm_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)^{}\right]{\displaystyle \frac{1}{m_W^4\mathrm{\Delta }_{Wk}^2}}D_{+k}^W`$ (187)
$`\times [8m_W^2m_\chi ^2(m_W^2m_\chi ^2)^2`$
$`+(m_W^2m_\chi ^2)(3m_W^4+20m_W^2m_\chi ^2+4m_\chi ^4)\mathrm{\Delta }_{Wk}].`$
$`\mathrm{𝟓}\mathbf{.}𝝌𝝌\mathbf{}𝒁𝒁`$

This process involves the $`s`$-channel CP–even Higgs boson ($`h`$ and $`H`$) exchange and the $`t`$\- and $`u`$-channel neutralino ($`\chi _i^0`$,$`i=1,\mathrm{},4`$) exchange

$`\stackrel{~}{a}_{ZZ}`$ $`=`$ $`\stackrel{~}{a}_{ZZ}^{(h,H)}+\stackrel{~}{a}_{ZZ}^{(\chi ^0)}+\stackrel{~}{a}_{ZZ}^{(h,H\chi ^0)},`$ (188)
$`\stackrel{~}{b}_{ZZ}`$ $`=`$ $`\stackrel{~}{b}_{ZZ}^{(h,H)}+\stackrel{~}{b}_{ZZ}^{(\chi ^0)}+\stackrel{~}{b}_{ZZ}^{(h,H\chi ^0)};`$ (189)
$``$ CP–even Higgs–boson $`(h,H)`$ exchange:

$`\stackrel{~}{a}_{ZZ}^{(h,H)}`$ $`=`$ $`0,`$ (190)
$`\stackrel{~}{b}_{ZZ}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{64\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{ZZr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2{\displaystyle \frac{3m_Z^44m_Z^2m_\chi ^2+4m_\chi ^4}{m_Z^4}};`$ (191)
$``$ neutralino ($`\chi _i^0`$) exchange:

$`\stackrel{~}{a}_{ZZ}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \underset{i,j=1}{\overset{4}{}}}C_{+i}^ZC_{+j}^Z{\displaystyle \frac{(m_\chi ^2m_Z^2)}{\mathrm{\Delta }_{Zi}\mathrm{\Delta }_{Zj}}},`$ (192)
$`\stackrel{~}{b}_{ZZ}^{(\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i,j=1}{\overset{4}{}}}{\displaystyle \frac{1}{m_Z^4\mathrm{\Delta }_{Zi}^3\mathrm{\Delta }_{Zj}^3}}`$
$`\times C_{+i}^ZC_{+j}^Z\left[D_{ij}^{(1)}\mathrm{\Delta }_{Zi}^2+D_{ij}^{(2)}\mathrm{\Delta }_{Zi}\mathrm{\Delta }_{Zj}+D_{ij}^{(3)}\mathrm{\Delta }_{Zi}^2\mathrm{\Delta }_{Zj}+D_{ij}^{(4)}\mathrm{\Delta }_{Zi}^2\mathrm{\Delta }_{Zj}^2\right],`$
where

$`D_{ij}^{(1)}`$ $`=`$ $`16m_\chi ^2m_Z^4(m_Z^2m_\chi ^2)^2,`$ (194)
$`D_{ij}^{(2)}`$ $`=`$ $`4m_\chi ^2(m_Z^2m_\chi ^2)^2\{3m_Z^4+4m_\chi ^2m_{\chi _i^0}m_{\chi _j^0}`$ (195)
$`+m_Z^2[4m_\chi ^2+4m_{\chi _i^0}m_{\chi _j^0}6m_\chi (m_{\chi _i^0}+m_{\chi _j^0})]\},`$
$`D_{ij}^{(3)}`$ $`=`$ $`4m_\chi (m_Z^2m_\chi ^2)\{4m_\chi ^3m_{\chi _i^0}m_{\chi _j^0}+m_Z^4(7m_\chi +3m_{\chi _i^0}+6m_{\chi _j^0})`$ (196)
$`+2m_Z^2m_\chi [4m_\chi ^2m_{\chi _i^0}m_{\chi _j^0}m_\chi (m_{\chi _i^0}+5m_{\chi _j^0})]\},`$
$`D_{ij}^{(4)}`$ $`=`$ $`8m_Z^6+2m_Z^2m_\chi ^2[4m_\chi ^26m_{\chi _i^0}m_{\chi _j^0}5m_\chi (m_{\chi _i^0}+m_{\chi _j^0})]`$ (197)
$`+4m_\chi ^4[2m_\chi ^2+3m_{\chi _i^0}m_{\chi _j^0}+2m_\chi (m_{\chi _i^0}+m_{\chi _j^0})]`$
$`+m_Z^4[3m_\chi ^2+9m_{\chi _i^0}m_{\chi _j^0}+5m_\chi (m_{\chi _i^0}+m_{\chi _j^0})],`$
where $`\mathrm{\Delta }_{Zi}m_Z^2m_\chi ^2m_{\chi _i^0}^2`$. The couplings $`C_{\pm i}^Z`$ and $`D_{\pm i}^Z`$ are given in eqs. (124) and (125).
$``$ Higgs $`(h,H)`$–neutralino ($`\chi _i^0`$) interference term:

$`\stackrel{~}{a}_{ZZ}^{(h,H\chi ^0)}`$ $`=`$ $`0,`$ (198)
$`\stackrel{~}{b}_{ZZ}^{(h,H\chi ^0)}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \underset{i=1}{\overset{4}{}}}Re\left[\left({\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C^{ZZr}C_S^{\chi \chi r}}{sm_r^2+i\mathrm{\Gamma }_rm_r}}\right)^{}\right]{\displaystyle \frac{1}{m_Z^4\mathrm{\Delta }_{Zi}^2}}C_S^{\chi \chi r}`$ (199)
$`\times \{2m_\chi (m_\chi ^2m_Z^2)[3m_Z^44m_\chi ^3m_{\chi _i^0}+2m_Z^2m_\chi (m_\chi +m_{\chi _i^0})].`$
$`+\mathrm{\Delta }_{Zi}[4m_\chi ^4(2m_\chi +3m_\chi )+2m_Z^2m_\chi ^2(5m_\chi +6m_{\chi _i^0})`$
$`m_Z^4(5m_\chi +9m_{\chi _i^0})]\}.`$
$`\mathrm{𝟔}\mathbf{.}𝝌𝝌\mathbf{}\overline{𝒇}𝒇`$

This process involves the $`s`$-channel Higgs boson ($`h`$, $`H`$ and $`A`$) and $`Z`$ boson exchange and the $`t`$\- and $`u`$-channel sfermion ($`\stackrel{~}{f}_a`$) exchange

$`\stackrel{~}{a}_{\overline{f}f}`$ $`=`$ $`\stackrel{~}{a}_{\overline{f}f}^{(h,H)}+\stackrel{~}{a}_{\overline{f}f}^{(A)}+\stackrel{~}{a}_{\overline{f}f}^{(Z)}+\stackrel{~}{a}_{\overline{f}f}^{(\stackrel{~}{f})}+\stackrel{~}{a}_{\overline{f}f}^{(h,H\stackrel{~}{f})}+\stackrel{~}{a}_{\overline{f}f}^{(AZ)}+\stackrel{~}{a}_{\overline{f}f}^{(A\stackrel{~}{f})}+\stackrel{~}{a}_{\overline{f}f}^{(Z\stackrel{~}{f})},`$ (200)
$`\stackrel{~}{b}_{\overline{f}f}`$ $`=`$ $`\stackrel{~}{b}_{\overline{f}f}^{(h,H)}+\stackrel{~}{b}_{\overline{f}f}^{(A)}+\stackrel{~}{b}_{\overline{f}f}^{(Z)}+\stackrel{~}{b}_{\overline{f}f}^{(\stackrel{~}{f})}+\stackrel{~}{b}_{\overline{f}f}^{(h,H\stackrel{~}{f})}+\stackrel{~}{b}_{\overline{f}f}^{(AZ)}+\stackrel{~}{b}_{\overline{f}f}^{(A\stackrel{~}{f})}+\stackrel{~}{b}_{\overline{f}f}^{(Z\stackrel{~}{f})}:`$ (201)
$``$ CP–even Higgs–boson $`(h,H)`$ exchange:

$`\stackrel{~}{a}_{\overline{f}f}^{(h,H)}`$ $`=`$ $`0,`$ (202)
$`\stackrel{~}{b}_{\overline{f}f}^{(h,H)}`$ $`=`$ $`{\displaystyle \frac{3}{4\pi }}\left|{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C_S^{ffr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right|^2(m_\chi ^2m_f^2);`$ (203)
$``$ CP–odd Higgs–boson ($`A`$) exchange:

$`\stackrel{~}{a}_{\overline{f}f}^{(A)}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}\left|{\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2m_\chi ^2,`$ (204)
$`\stackrel{~}{b}_{\overline{f}f}^{(A)}`$ $`=`$ $`{\displaystyle \frac{3}{2\pi }}\left|{\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A}}\right|^2{\displaystyle \frac{m_\chi ^2m_A^2(m_A^24m_\chi ^2+\mathrm{\Gamma }_A^2)}{(4m_\chi ^2m_A^2)^2+(\mathrm{\Gamma }_Am_A)^2}};`$ (205)
$``$ $`Z`$–boson exchange:

$`\stackrel{~}{a}_{\overline{f}f}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}\left|{\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2{\displaystyle \frac{m_f^2(m_Z^24m_\chi ^2)^2}{m_Z^4}},`$ (206)
$`\stackrel{~}{b}_{\overline{f}f}^{(Z)}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}\left|{\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right|^2{\displaystyle \frac{1}{m_Z^2((4m_\chi ^2m_Z^2)^2+(\mathrm{\Gamma }_Zm_Z)^2)}}`$ (207)
$`\times [2|C_A^{ffZ}|^2\{m_Z^2(m_\chi ^2m_f^2)(m_Z^24m_\chi ^2)^2`$
$`+\mathrm{\Gamma }_Z^2[m_\chi ^2m_Z^4+m_f^2(24m_\chi ^46m_Z^2m_\chi ^2m_Z^4)]\}`$
$`+m_Z^2|C_V^{ffZ}|^2\left\{(2m_\chi ^2+m_f^2)[(4m_\chi ^2m_Z^2)^2+m_Z^2\mathrm{\Gamma }_Z^2]\right\}];`$
$``$ sfermion ($`\stackrel{~}{f}_a`$) exchange:

$`\stackrel{~}{a}_{\overline{f}f}^{(\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{32\pi }}{\displaystyle \underset{a,b}{}}{\displaystyle \frac{(m_fC_+^a+m_\chi D_+^a)(m_fC_+^b+m_\chi D_+^b)}{\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}}},`$ (208)
$`\stackrel{~}{b}_{\overline{f}f}^{(\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{64\pi }}{\displaystyle \underset{a,b}{}}{\displaystyle \frac{1}{\mathrm{\Delta }_{\stackrel{~}{f}_a}^3\mathrm{\Delta }_{\stackrel{~}{f}_b}^3}}\{C_+^aC_+^b[8m_f^2m_\chi ^2(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2`$
$`4m_\chi ^2(m_f^4+m_\chi ^2m_f^22m_\chi ^4)\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}`$
$`+4m_\chi ^2(m_f^2+2m_\chi ^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}+4(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}^2]`$
$`+D_+^aD_+^b[8m_\chi ^4(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2+4m_\chi ^2(m_\chi ^4+m_\chi ^2m_f^22m_f^4)\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}`$
$`+4m_\chi ^2(5m_\chi ^22m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}3(m_f^23m_\chi ^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}^2]`$
$`+C_{}^aC_{}^b[8m_\chi ^2(m_\chi ^2m_f^2)^2\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}+8m_\chi ^2(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}`$
$`+2(m_f^2+2m_\chi ^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}^2]`$
$`+C_+^aD_+^b\mathrm{\hspace{0.17em}2}m_fm_\chi [8m_\chi ^2(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2+12m_\chi ^2(m_\chi ^2m_f^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}`$
$`+3\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}^22(m_f^24m_\chi ^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}\mathrm{\Delta }_{\stackrel{~}{f}_b}^26(m_f^22m_\chi ^2)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2\mathrm{\Delta }_{\stackrel{~}{f}_b}]\},`$
where $`\mathrm{\Delta }_{\stackrel{~}{f}_a}m_f^2m_\chi ^2m_{\stackrel{~}{f}_a}^2`$. The index $`a`$ counts sfermions so that $`a`$ $`=`$ 1, $`\mathrm{}`$, 6 for squarks and charged sleptons, and $`a`$ $`=`$ 1, 2, 3 for sneutrinos. The couplings $`C_\pm ^a`$ and $`D_\pm ^a`$ are given in eqs. (134) and (135).
$``$ Higgs $`(h,H)`$–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{a}_{\overline{f}f}^{(h,H\stackrel{~}{f})}`$ $`=`$ $`0,`$ (210)
$`\stackrel{~}{b}_{\overline{f}f}^{(h,H\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{a}{}}Re\left[{\displaystyle \underset{r=h,H}{}}{\displaystyle \frac{C_S^{ffr}C_S^{\chi \chi r}}{4m_\chi ^2m_r^2+i\mathrm{\Gamma }_rm_r}}\right]{\displaystyle \frac{(m_\chi ^2m_f^2)}{\mathrm{\Delta }_{\stackrel{~}{f}_a}^2}}`$ (211)
$`\times \left[C_+^a\mathrm{\hspace{0.17em}2}m_fm_\chi +D_+^a(2m_\chi ^2+3\mathrm{\Delta }_{\stackrel{~}{f}_a})\right];`$
$``$ Higgs ($`A`$)–$`Z`$ interference term:

$`\stackrel{~}{a}_{\overline{f}f}^{(AZ)}`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}Re\left[\left({\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}\left({\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)\right]`$ (212)
$`\times {\displaystyle \frac{m_\chi m_f(m_Z^24m_\chi ^2)}{m_Z^2}},`$
$`\stackrel{~}{b}_{\overline{f}f}^{(AZ)}`$ $`=`$ $`{\displaystyle \frac{3}{2\pi }}Re[\left({\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A}}\right)^{}\left({\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right)`$
$`\times {\displaystyle \frac{m_\chi m_f}{m_Z^2}}[m_A(m_Z^24m_\chi ^2)^2(m_A+i\mathrm{\Gamma }_A)`$
$`+m_Z\mathrm{\Gamma }_Z\{m_A\mathrm{\Gamma }_A(m_Z^28m_\chi ^2)i[16m_\chi ^4+m_A^2(m_Z^28m_\chi ^2)]\}]];`$
$``$ Higgs ($`A`$)–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{a}_{\overline{f}f}^{(A\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \underset{a}{}}Re\left[{\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A}}\right]{\displaystyle \frac{(m_fC_+^a+m_\chi D_+^a)}{\mathrm{\Delta }_{\stackrel{~}{f}_a}}},`$ (214)
$`\stackrel{~}{b}_{\overline{f}f}^{(A\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{a}{}}Re[\left({\displaystyle \frac{C_P^{ffA}C_P^{\chi \chi A}}{(4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A)^2}}\right){\displaystyle \frac{m_\chi }{\mathrm{\Delta }_{\stackrel{~}{f}_a}^3}}`$ (215)
$`\times \{C_+^a[4m_fm_\chi ^2(m_\chi ^2m_f^2)P_A+6m_fm_\chi ^2P_A\mathrm{\Delta }_{\stackrel{~}{f}_a}`$
$`3m_fm_A(m_Ai\mathrm{\Gamma }_A)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2]`$
$`+D_+^a[4m_\chi ^3(m_\chi ^2m_f^2)P_A2m_\chi (m_f^24m_\chi ^2)P_A\mathrm{\Delta }_{\stackrel{~}{f}_a}`$
$`+6m_\chi (2m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A)\mathrm{\Delta }_{\stackrel{~}{f}_a}^2]\}];`$
where $`P_A4m_\chi ^2m_A^2+i\mathrm{\Gamma }_Am_A.`$
$``$ $`Z`$–sfermion ($`\stackrel{~}{f}_a`$) interference term:

$`\stackrel{~}{a}_{\overline{f}f}^{(Z\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \underset{a}{}}Re\left[{\displaystyle \frac{C_A^{ffZ}C_A^{\chi \chi Z}}{4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z}}\right]{\displaystyle \frac{m_f(m_Z^24m_\chi ^2)}{m_Z^2}}`$ (216)
$`\times {\displaystyle \frac{(m_fC_+^a+m_\chi D_+^a)}{\mathrm{\Delta }_{\stackrel{~}{f}_a}}},`$
$`\stackrel{~}{b}_{\overline{f}f}^{(Z\stackrel{~}{f})}`$ $`=`$ $`{\displaystyle \frac{1}{8\pi }}{\displaystyle \underset{a}{}}Re[\left({\displaystyle \frac{C_A^{\chi \chi Z}}{(4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z)^2}}\right){\displaystyle \frac{1}{m_Z^2\mathrm{\Delta }_{\stackrel{~}{f}_a}^3}}`$
$`\times [C_V^{ffZ}C_{}^a\left\{2m_Z^2P_Z\mathrm{\Delta }_{\stackrel{~}{f}_a}[2m_\chi ^2(m_\chi ^2+\mathrm{\Delta }_{\stackrel{~}{f}_a})+m_f^2(2m_\chi ^2+\mathrm{\Delta }_{\stackrel{~}{f}_a})]\right\}`$
$`+C_A^{ffZ}\{C_+^a[2m_f^2m_\chi ^2(m_\chi ^2m_f^2)(m_Z^24m_\chi ^2)P_Z`$
$`+m_\chi ^2[m_f^2m_Z^2+\mathrm{\hspace{0.17em}2}m_\chi ^2(m_Z^26m_f^2)]P_Z\mathrm{\Delta }_{\stackrel{~}{f}_a}`$
$`+\mathrm{\hspace{0.17em}2}m_Z\{m_Z(m_\chi ^2m_f^2)(m_Z^24m_\chi ^2)`$
$`+i\mathrm{\Gamma }_Z[m_Z^2m_\chi ^2m_f^2(m_Z^2+3m_\chi ^2)]\}\mathrm{\Delta }_{\stackrel{~}{f}_a}^2]`$
$`+m_fm_\chi D_+^a[4m_\chi ^2(m_\chi ^2m_f^2)(m_Z^24m_\chi ^2)P_Z`$
$`+2[6m_Z^2m_\chi ^216m_\chi ^4m_f^2(3m_Z^24m_\chi ^2)]P_Z\mathrm{\Delta }_{\stackrel{~}{f}_a}`$
$`3[(m_Z^24m_\chi ^2)^2im_Z\mathrm{\Gamma }_Z(m_Z^28m_\chi ^2)]\mathrm{\Delta }_{\stackrel{~}{f}_a}^2]\}]];`$
where $`P_Z4m_\chi ^2m_Z^2+i\mathrm{\Gamma }_Zm_Z`$.

## 6 Summary

The neutralino is undoubtedly the most popular candidate for a WIMP dark matter in the Universe. Supersymmetry remains arguably the most promising extension of the Standard Model. The next several years will witness extensive searches for supersymmetry in colliders as well as for WIMPs in underground detectors. Measurements of the cosmological parameters, and in particular of the relic abundance of the dark matter, have already reached the accuracy of a few per cent and more progress is expected.

In light of this, theoretical computations of the neutralino relic abundance need to be now performed with at least the same, if not better, level of precision, if one wants to reliably compare theoretical predictions with observations. Motivated by this goal, we have derived a full set of exact, analytic expressions for the neutralino pair-annihilation cross sections into all tree-level two-body final states in the framework of the MSSM.

###### Acknowledgments.

T.N. is grateful to Lancaster University for kind hospitality extended during his visit. The authors appreciate Dr. A. Birkedal-Hansen to point out several minor typographical errors in the first version of this paper.

## Appendix A The MSSM Lagrangian

In this Appendix, the MSSM Lagrangian is given explicitly in the mass eigenstates, which makes it easier to read out relevant Feynman rules.

Z-boson-Fermion-Fermion

$``$ $`=`$ $`{\displaystyle \underset{f}{}}\overline{f}\gamma ^\mu \left[C_V^{ffZ}C_A^{ffZ}\gamma _5\right]fZ_\mu ,`$ (218)
where

$`C_V^{ffZ}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}(T_{3f}2\mathrm{sin}^2\theta _WQ_f),`$ (219)
$`C_A^{ffZ}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}T_{3f},`$ (220)
and $`g`$ is the usual gauge coupling of $`SU(2)_L`$ whereas $`Q_f`$ and $`T_{3f}`$ were defined below eq. (21).

W-Chargino-Neutralino

$``$ $`=`$ $`{\displaystyle \underset{k=1}{\overset{2}{}}}{\displaystyle \underset{i=1}{\overset{4}{}}}\overline{\chi _k^{}}\gamma ^\mu \left[C_V^{\chi _k^+\chi _i^0W^{}}C_A^{\chi _k^+\chi _i^0W^{}}\gamma _5\right]\chi _i^0W_\mu ^{}+h.c.,`$ (221)
where

$`C_V^{\chi _k^+\chi _i^0W^{}}`$ $`=`$ $`\frac{1}{2}g(O_{ik}^L+O_{ik}^R),`$ (222)
$`C_A^{\chi _k^+\chi _i^0W^{}}`$ $`=`$ $`\frac{1}{2}g(O_{ik}^LO_{ik}^R),`$ (223)
and

$`O_{ik}^L`$ $`=`$ $`\sqrt{\frac{1}{2}}N_{i4}V_{k2}^{}+N_{i2}V_{k1}^{},`$ (224)
$`O_{ik}^R`$ $`=`$ $`\sqrt{\frac{1}{2}}N_{i3}^{}U_{k2}+N_{i2}^{}U_{k1}^{},`$ (225)
where the matrices $`N,U,V`$ were defined via eqs. (14) and (18).

Higgs-Chargino-Neutralino

$``$ $`=`$ $`{\displaystyle \underset{k=1}{\overset{2}{}}}{\displaystyle \underset{i=1}{\overset{4}{}}}\overline{\chi _k^{}}\left[C_S^{\chi _k^+\chi _i^0H^{}}C_P^{\chi _k^+\chi _i^0H^{}}\gamma _5\right]\chi _i^0H^{}+h.c.,`$ (226)
where

$`C_S^{\chi _k^+\chi _i^0H^{}}`$ $`=`$ $`\frac{1}{2}(Q_{ik}^{}_{}{}^{}L+Q_{ik}^{}_{}{}^{}R),`$ (227)
$`C_P^{\chi _k^+\chi _i^0H^{}}`$ $`=`$ $`\frac{1}{2}(Q_{ik}^{}_{}{}^{}LQ_{ik}^{}_{}{}^{}R),`$ (228)
and

$`Q_{ik}^{}_{}{}^{}L`$ $`=`$ $`g\mathrm{cos}\beta \left[N_{i4}^{}V_{k1}^{}+\sqrt{\frac{1}{2}}(N_{i2}^{}+N_{i1}^{}\mathrm{tan}\theta _W)V_{k2}^{}\right],`$ (229)
$`Q_{ik}^{}_{}{}^{}R`$ $`=`$ $`g\mathrm{sin}\beta \left[N_{i3}U_{k1}\sqrt{\frac{1}{2}}(N_{i2}+N_{i1}\mathrm{tan}\theta _W)U_{k2}\right].`$ (230)
Z-Neutralino-Neutralino

$``$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}\overline{\chi _i^0}\gamma ^\mu \left[C_V^{\chi _i^0\chi _j^0Z}C_A^{\chi _i^0\chi _j^0Z}\gamma _5\right]\chi _j^0Z_\mu ,`$ (231)
where

$`C_V^{\chi _i^0\chi _j^0Z}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}(O_{ij}^{{}_{}{}^{\prime \prime }L}O_{ij}^{{}_{}{}^{\prime \prime }L}),`$ (232)
$`C_A^{\chi _i^0\chi _j^0Z}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}(O_{ij}^{{}_{}{}^{\prime \prime }L}+O_{ij}^{{}_{}{}^{\prime \prime }L}),`$ (233)
and

$`O_{ij}^{{}_{}{}^{\prime \prime }L}`$ $`=`$ $`\frac{1}{2}(N_{i3}N_{j3}^{}+N_{i4}N_{j4}^{}).`$ (234)
We neglect CP violation hence $`C_V^{\chi _i^0\chi _j^0Z}`$ $`=`$ 0.

Gauge-Higgs-Higgs

$``$ $`=`$ $`i\left[C^{hAZ}(A\stackrel{}{^\mu }h)+C^{HAZ}(A\stackrel{}{^\mu }H)+C^{H^+H^{}Z}(H^+\stackrel{}{^\mu }H^{})\right]Z_\mu `$ (235)
$`+\{iW_\mu ^{}[C^{W^{}H^+h}(h\stackrel{}{^\mu }H^+)+C^{W^{}H^+H}(H\stackrel{}{^\mu }H^+)`$
$`+C^{W^{}H^+A}(A\stackrel{}{^\mu }H^+)]+h.c.\},`$
where

$`C^{hAZ}`$ $`=`$ $`{\displaystyle \frac{ig}{2\mathrm{cos}\theta _W}}\mathrm{cos}(\alpha \beta ),`$ (236)
$`C^{HAZ}`$ $`=`$ $`{\displaystyle \frac{ig}{2\mathrm{cos}\theta _W}}\mathrm{sin}(\alpha \beta ),`$ (237)
$`C^{H^+H^{}Z}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}(12\mathrm{sin}^2\theta _W),`$ (238)
$`C^{W^{}H^+h}`$ $`=`$ $`{\displaystyle \frac{g}{2}}\mathrm{cos}(\alpha \beta ),`$ (239)
$`C^{W^{}H^+H}`$ $`=`$ $`{\displaystyle \frac{g}{2}}\mathrm{sin}(\alpha \beta ),`$ (240)
$`C^{W^{}H^+A}`$ $`=`$ $`{\displaystyle \frac{ig}{2}}.`$ (241)
Higgs-Fermion-Fermion

$``$ $`=`$ $`{\displaystyle \underset{f}{}}\left[C_S^{ffh}\overline{f}fh+C_S^{ffH}\overline{f}fH+C_P^{ffA}\overline{f}\gamma _5fA\right]`$
$``$ $`{\displaystyle \underset{n=1}{\overset{3}{}}}[{\displaystyle \frac{gm_{e_n}\mathrm{sin}\alpha }{2m_W\mathrm{cos}\beta }}\overline{e}_ne_nh{\displaystyle \frac{gm_{e_n}\mathrm{cos}\alpha }{2m_W\mathrm{cos}\beta }}\overline{e}_ne_nH+i{\displaystyle \frac{gm_{e_n}}{2m_W}}\mathrm{tan}\beta \overline{e}_n\gamma _5e_nA`$
$`{\displaystyle \frac{gm_{u_n}\mathrm{cos}\alpha }{2m_W\mathrm{sin}\beta }}\overline{u}_nu_nh{\displaystyle \frac{gm_{u_n}\mathrm{sin}\alpha }{2m_W\mathrm{sin}\beta }}\overline{u}_nu_nH+i{\displaystyle \frac{gm_{u_n}}{2m_W}}\mathrm{cot}\beta \overline{u}_n\gamma _5u_nA`$
$`+{\displaystyle \frac{gm_{d_n}\mathrm{sin}\alpha }{2m_W\mathrm{cos}\beta }}\overline{d}_nd_nh{\displaystyle \frac{gm_{d_n}\mathrm{cos}\alpha }{2m_W\mathrm{cos}\beta }}\overline{d}_nd_nH+i{\displaystyle \frac{gm_{d_n}}{2m_W}}\mathrm{tan}\beta \overline{d}_n\gamma _5d_nA],`$
where $`n`$ $`=`$ $`1,2,3`$ is an index for generation, so $`e_n=e,\mu ,\tau `$, etc.

Higgs-Gauge-Gauge

$``$ $`=`$ $`\left[C^{W^+W^{}h}h+C^{W^+W^{}H}H\right]W_\mu ^+W^\mu +{\displaystyle \frac{1}{2}}\left[C^{ZZh}h+C^{ZZH}H\right]Z_\mu Z^\mu ,`$
where

$`C^{W^+W^{}h}`$ $`=`$ $`gm_W\mathrm{sin}(\alpha \beta ),`$ (245)
$`C^{W^+W^{}H}`$ $`=`$ $`gm_W\mathrm{cos}(\alpha \beta ),`$ (246)
$`C^{ZZh}`$ $`=`$ $`{\displaystyle \frac{gm_Z}{\mathrm{cos}\theta _W}}\mathrm{sin}(\alpha \beta ),`$ (247)
$`C^{ZZH}`$ $`=`$ $`{\displaystyle \frac{gm_Z}{\mathrm{cos}\theta _W}}\mathrm{cos}(\alpha \beta ).`$ (248)
Higgs-Neutralino-Neutralino

$``$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{i,j=1}{\overset{4}{}}}\overline{\chi _i^0}\left[C_S^{\chi _i^0\chi _j^0H}H+C_S^{\chi _i^0\chi _j^0h}h+C_P^{\chi _i^0\chi _j^0A}\gamma _5A\right]\chi _j^0,`$ (249)
where

$`C_S^{\chi _i^0\chi _j^0H}`$ $`=`$ $`{\displaystyle \frac{g}{2}}\{[N_{i2}N_{i1}\mathrm{tan}\theta _W][\mathrm{cos}\alpha N_{j3}\mathrm{sin}\alpha N_{j4}]+(ij)\},`$ (250)
$`C_S^{\chi _i^0\chi _j^0h}`$ $`=`$ $`{\displaystyle \frac{g}{2}}\{[N_{i2}N_{i1}\mathrm{tan}\theta _W][\mathrm{sin}\alpha N_{j3}\mathrm{cos}\alpha N_{j4}]+(ij)\},`$ (251)
$`C_P^{\chi _i^0\chi _j^0A}`$ $`=`$ $`{\displaystyle \frac{ig}{2}}\{[N_{i2}N_{i1}\mathrm{tan}\theta _W][\mathrm{sin}\beta N_{j3}\mathrm{cos}\beta N_{j4}]+(ij)\}.`$ (252)
Note that we assume no CP violating phases in the neutralino mass matrix and we use the convention that $`N_{ij}`$ is real.
Gauge-Gauge-Gauge

$``$ $`=`$ $`iC^{WWZ}[(_\mu W_\nu ^{}_\nu W_\mu ^{})W^{+\mu }Z^\nu (_\mu W_\nu ^+_\nu W_\mu ^+)W^\mu Z^\nu `$ (253)
$`(_\mu Z_\nu _\nu Z_\mu )W^{+\mu }W^\nu ],`$
where

$`C^{WWZ}`$ $`=`$ $`g\mathrm{cos}\theta _W.`$ (254)
Neutral Higgs-Higgs-Higgs

$``$ $`=`$ $`{\displaystyle \frac{1}{6}}C^{hhh}h^3+{\displaystyle \frac{1}{6}}C^{HHH}H^3+{\displaystyle \frac{1}{2}}C^{hhH}h^2H+{\displaystyle \frac{1}{2}}C^{HHh}H^2h`$ (255)
$`+{\displaystyle \frac{1}{2}}(C^{AAh}h+C^{AAH}H)A^2,`$
where

$`C^{hhh}`$ $`=`$ $`{\displaystyle \frac{3g}{2\mathrm{cos}\theta _W}}m_Z\mathrm{cos}2\alpha \mathrm{sin}(\alpha +\beta ),`$ (256)
$`C^{HHH}`$ $`=`$ $`{\displaystyle \frac{3g}{2\mathrm{cos}\theta _W}}m_Z\mathrm{cos}2\alpha \mathrm{cos}(\alpha +\beta ),`$ (257)
$`C^{AAh}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}m_Z\mathrm{cos}2\beta \mathrm{sin}(\alpha +\beta ),`$ (258)
$`C^{AAH}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}m_Z\mathrm{cos}2\beta \mathrm{cos}(\alpha +\beta ),`$ (259)
$`C^{hhH}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}m_Z[2\mathrm{sin}2\alpha \mathrm{sin}(\alpha +\beta )\mathrm{cos}2\alpha \mathrm{cos}(\alpha +\beta )],`$ (260)
$`C^{HHh}`$ $`=`$ $`{\displaystyle \frac{g}{2\mathrm{cos}\theta _W}}m_Z[2\mathrm{sin}2\alpha \mathrm{cos}(\alpha +\beta )+\mathrm{cos}2\alpha \mathrm{sin}(\alpha +\beta )].`$ (261)
H<sup>+</sup>-H<sup>-</sup>-Higgs

$``$ $`=`$ $`(C^{H^+H^{}h}h+C^{H^+H^{}H}H)H^+H^{},`$ (262)
where

$`C^{H^+H^{}h}`$ $`=`$ $`gm_Z\left[\mathrm{cos}\theta _W\mathrm{sin}(\alpha \beta ){\displaystyle \frac{1}{2\mathrm{cos}\theta _W}}\mathrm{sin}(\alpha +\beta )\mathrm{cos}2\beta \right],`$ (263)
$`C^{H^+H^{}H}`$ $`=`$ $`gm_Z\left[\mathrm{cos}\theta _W\mathrm{cos}(\alpha \beta ){\displaystyle \frac{1}{2\mathrm{cos}\theta _W}}\mathrm{cos}(\alpha +\beta )\mathrm{cos}2\beta \right].`$ (264)
Neutralino-Fermion-Sfermion

$``$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{4}{}}}{\displaystyle \underset{n=1}{\overset{3}{}}}\{{\displaystyle \underset{f=u,d,e}{}}{\displaystyle \underset{a=1}{\overset{6}{}}}\overline{\chi _i^0}(\mathrm{\Lambda }_{nai}^{(f)L}{\displaystyle \frac{1\gamma _5}{2}}+\mathrm{\Lambda }_{nai}^{(f)R}{\displaystyle \frac{1+\gamma _5}{2}})f_n\stackrel{~}{f}_a^{}`$ (265)
$`+{\displaystyle \underset{a=1}{\overset{3}{}}}\overline{\chi _i^0}\mathrm{\Lambda }_{nai}^{(\nu )L}{\displaystyle \frac{1\gamma _5}{2}}\nu _n\stackrel{~}{\nu }_a^{}+h.c.\},`$
where

$`\mathrm{\Lambda }_{nai}^{(f)L}`$ $`=`$ $`{\displaystyle \frac{g}{\sqrt{2}}}\left(\mathrm{\Gamma }_{(f)L}^{an}\mathrm{}^{\chi _i^0f_nL}+\mathrm{\Gamma }_{(f)R}^{an}\mathrm{}^{\chi _i^0f_nR}\right),`$ (266)
$`\mathrm{\Lambda }_{nai}^{(f)R}`$ $`=`$ $`{\displaystyle \frac{g}{\sqrt{2}}}\left(\mathrm{\Gamma }_{(f)L}^{an}r^{\chi _i^0f_nL}+\mathrm{\Gamma }_{(f)R}^{an}r^{\chi _i^0f_nR}\right).`$ (267)
Note that in eqs. (266) and (267), $`f`$ $`=`$ $`u`$, $`d`$, $`e`$, $`\nu `$ represents the type of the fermion, while $`f`$ in the main text denotes each fermion $`f`$ $`=`$ $`u`$, $`c`$, $`t`$, etc. The slepton and squark mass eigenstates $`\stackrel{~}{f}_a`$ ($`\stackrel{~}{\nu }_a`$ with $`a=1,2,3`$ and $`\stackrel{~}{e}_a`$,$`\stackrel{~}{u}_a`$ and $`\stackrel{~}{d}_a`$ with $`a=1,\mathrm{},6`$) are related to the sfermion gauge eigenstates $`\stackrel{~}{f}_{nL}`$ and $`\stackrel{~}{f}_{nR}`$ ($`n=1,2,3`$) via

$`\stackrel{~}{f}_{nL}={\displaystyle \underset{a}{}}\stackrel{~}{f}_a\mathrm{\Gamma }_{(f)L}^{an},`$ (268)
$`\stackrel{~}{f}_{nR}={\displaystyle \underset{a}{}}\stackrel{~}{f}_a\mathrm{\Gamma }_{(f)R}^{an},`$ (269)
where $`\stackrel{~}{V}_f`$ $`=`$ $`(\mathrm{\Gamma }_{(f)L},\mathrm{\Gamma }_{(f)R})`$ denotes a $`6\times 6`$ matrix which diagonalizes the sfermions mass matrix given by eqs. (22)–(27): $`\stackrel{~}{V}_f`$ $`_{\stackrel{~}{f}}^2`$ $`\stackrel{~}{V}_f^{}`$ $`=`$ $`\mathrm{diag}(m_{\stackrel{~}{f}_1}^2,\mathrm{},m_{\stackrel{~}{f}_6}^2)`$. Note that in the case of squarks and charged sleptons the mixing matrices $`\mathrm{\Gamma }_{(u)L,R},\mathrm{\Gamma }_{(d)L,R}`$ and $`\mathrm{\Gamma }_{(e)L,R}`$ have dimension $`6\times 3`$, while for sneutrinos $`\mathrm{\Gamma }_{(\nu )L}`$ is a mixing matrix of order $`3\times 3`$. (We neglect here the CKM matrix for simplicity.) Finally $`\mathrm{}^{\chi _i^0f_nL,R}`$ and $`r^{\chi _i^0f_nL,R}`$ are given as follows:

$`\mathrm{}^{\chi _i^0\nu _nL}`$ $`=`$ $`N_{i2}\mathrm{tan}\theta _WN_{i1},`$ (270)
$`\mathrm{}^{\chi _i^0e_nL}`$ $`=`$ $`N_{i2}\mathrm{tan}\theta _WN_{i1},`$ (271)
$`\mathrm{}^{\chi _i^0u_nL}`$ $`=`$ $`N_{i2}+{\displaystyle \frac{1}{3}}\mathrm{tan}\theta _WN_{i1},`$ (272)
$`\mathrm{}^{\chi _i^0d_nL}`$ $`=`$ $`N_{i2}+{\displaystyle \frac{1}{3}}\mathrm{tan}\theta _WN_{i1},`$ (273)
$`r^{\chi _i^0\nu _nL}`$ $`=`$ $`0,`$ (274)
$`r^{\chi _i^0e_nL}`$ $`=`$ $`{\displaystyle \frac{m_{e_n}}{m_W\mathrm{cos}\beta }}N_{i3},`$ (275)
$`r^{\chi _i^0u_nL}`$ $`=`$ $`{\displaystyle \frac{m_{u_n}}{m_W\mathrm{sin}\beta }}N_{i4},`$ (276)
$`r^{\chi _i^0d_nL}`$ $`=`$ $`{\displaystyle \frac{m_{d_n}}{m_W\mathrm{cos}\beta }}N_{i3},`$ (277)
$`\mathrm{}^{\chi _i^0\nu _nR}`$ $`=`$ $`0,`$ (278)
$`\mathrm{}^{\chi _i^0e_nR}`$ $`=`$ $`{\displaystyle \frac{m_{e_n}}{m_W\mathrm{cos}\beta }}N_{i3},`$ (279)
$`\mathrm{}^{\chi _i^0u_nR}`$ $`=`$ $`{\displaystyle \frac{m_{u_n}}{m_W\mathrm{sin}\beta }}N_{i4},`$ (280)
$`\mathrm{}^{\chi _i^0d_nR}`$ $`=`$ $`{\displaystyle \frac{m_{d_n}}{m_W\mathrm{cos}\beta }}N_{i3},`$ (281)
$`r^{\chi _i^0\nu _nR}`$ $`=`$ $`0,`$ (282)
$`r^{\chi _i^0e_nR}`$ $`=`$ $`2\mathrm{tan}\theta _WN_{i1},`$ (283)
$`r^{\chi _i^0u_nR}`$ $`=`$ $`{\displaystyle \frac{4}{3}}\mathrm{tan}\theta _WN_{i1},`$ (284)
$`r^{\chi _i^0d_nR}`$ $`=`$ $`{\displaystyle \frac{2}{3}}\mathrm{tan}\theta _WN_{i1}.`$ (285)
## Appendix B Auxiliary Functions

Here we give expressions for the auxiliary functions used in the text. First, we define

$`D(s,x,y_1,y_2)`$ $``$ $`x+{\displaystyle \frac{y_1+y_2}{2}}{\displaystyle \frac{s}{2}},`$ (286)
$`F(s,x,y_1,y_2)`$ $``$ $`{\displaystyle \frac{1}{2}}\sqrt{s4x}\sqrt{s(\sqrt{y_1}+\sqrt{y_2})^2}\sqrt{1{\displaystyle \frac{(\sqrt{y_1}\sqrt{y_2})^2}{s}}}.`$ (287)
If we define $`DD(s,x,y_1,y_2)`$, $`FF(s,x,y_1,y_2)`$, $`t_\pm (s,x,y_1,y_2)D\pm F`$ and $`(𝒯_i,𝒴_i)`$ $``$ $`(𝒯_i,𝒴_i)(s,x,y_1,y_2,z_1,z_2)`$ ($`i`$ $`=`$ $`0,\mathrm{},4`$), then

$`(s,x,y_1,y_2,z)`$ $`=`$ $`{\displaystyle \frac{1}{2F}}\mathrm{ln}\left|{\displaystyle \frac{t_+(s,x,y_1,y_2)z}{t_{}(s,x,y_1,y_2)z}}\right|,`$ (288)
and

$`𝒯_0`$ $`=`$ $`{\displaystyle \frac{1}{z_1z_2}}[(s,x,y_1,y_2,z_1)(s,x,y_1,y_2,z_2)],`$ (289)
$`𝒯_1`$ $`=`$ $`{\displaystyle \frac{1}{z_1z_2}}[z_1(s,x,y_1,y_2,z_1)z_2(s,x,y_1,y_2,z_2)],`$ (290)
$`𝒯_2`$ $`=`$ $`1+{\displaystyle \frac{1}{z_1z_2}}[z_1^2(s,x,y_1,y_2,z_1)z_2^2(s,x,y_1,y_2,z_2)],`$ (291)
$`𝒯_3`$ $`=`$ $`D+z_1+z_2+{\displaystyle \frac{1}{z_1z_2}}[z_1^3(s,x,y_1,y_2,z_1)z_2^3(s,x,y_1,y_2,z_2)],`$ (292)
$`𝒯_4`$ $`=`$ $`D^2+D(z_1+z_2)+(z_1^2+z_2^2+z_1z_2)+{\displaystyle \frac{1}{3}}F^2`$ (293)
$`+{\displaystyle \frac{1}{z_1z_2}}[z_1^4(s,x,y_1,y_2,z_1)z_2^4(s,x,y_1,y_2,z_2)],`$
$`𝒴_0`$ $`=`$ $`{\displaystyle \frac{1}{z_1+z_22D}}[(s,x,y_1,y_2,z_1)+(s,x,y_1,y_2,z_2)],`$ (294)
$`𝒴_1`$ $`=`$ $`{\displaystyle \frac{1}{z_1+z_22D}}[2(z_1D)(s,x,y_1,y_2,z_1)2(z_2D)(s,x,y_1,y_2,z_2)],`$
$`𝒴_2`$ $`=`$ $`1+{\displaystyle \frac{1}{z_1+z_22D}}[z_1(z_12D)(s,x,y_1,y_2,z_1)`$ (296)
$`+z_2(z_22D)(s,x,y_1,y_2,z_2)],`$
$`𝒴_3`$ $`=`$ $`2(z_1z_2)+{\displaystyle \frac{1}{z_1+z_22D}}[z_1(2z_1^26z_1D+4D^2)(s,x,y_1,y_2,z_1)`$ (297)
$`z_2(2z_2^26z_2D+4D^2)(s,x,y_1,y_2,z_2)],`$
$`𝒴_4`$ $`=`$ $`D^2D(z_1+z_2)+(z_1^2+z_2^2z_1z_2)+{\displaystyle \frac{1}{3}}F^2`$ (298)
$`+{\displaystyle \frac{1}{z_1+z_22D}}[z_1^2(z_1^24z_1D+4D^2)(s,x,y_1,y_2,z_1)`$
$`+z_2^2(z_2^24z_2D+4D^2)(s,x,y_1,y_2,z_2)].`$