Breit-Wigner approximation for propagators of mixed unstable states

For systems of unstable particles that mix with each other, an approximation of the fully momentum- dependent propagator matrix is presented in terms of a sum of simple Breit-Wigner propagators that are multiplied with finite on-shell wave function normalisation factors. The latter are evaluated at the complex poles of the propagators. The pole structure of general propagator matrices is carefully analysed, and it is demonstrated that in the proposed approximation imaginary parts arising from absorptive parts of loop integrals are properly taken into account. Applying the formalism to the neutral MSSM Higgs sector with complex parameters, very good numerical agreement is found between cross sections based on the full propagators and the corresponding cross sections based on the described approximation. The proposed approach does not only technically simplify the treatment of propagators with non-vanishing off-diagonal contributions, it is shown that it can also facilitate an improved theoretical prediction of the considered observables via a more precise implementation of the total widths of the involved particles. It is also well-suited for the incorporation of interference effects arising from overlapping resonances.


JHEP09(2017)079
As a consequence of electroweak symmetry breaking, particles with the same quantum numbers of the electric charge and colour can mix with each other. This is reflected in their propagators contributing to production or decay processes in particle physics. Beyond lowest order those propagators receive contributions from all possible mixings with other particles, such that the propagators of the system of particles that can mix with each other are of matrix-type. Since the involved particles are in general unstable, a proper treatment of imaginary parts is necessary. The mass of an unstable particle is determined from the real part of the (gauge-invariant) complex pole [1][2][3][4], while the imaginary part yields the total decay width of the particle. The eigenstates associated with the masses obtained from the complex poles of the propagator matrix in general differ from the interaction eigenstates that contribute at lowest order.
As an example, in the Standard Model (SM) the neutral interaction eigenstates arising from the SU(2) I and U(1) Y gauge groups mix with each other to form the mass eigenstates photon, γ, and Z. Higher-order contributions lead to a mixing between these states, giving rise to a γ − Z propagator matrix. In a renormalisable gauge, the (would-be) Goldstone boson as an unphysical degree of freedom does not vanish, and mixing contributions with this unphysical scalar particle need to be taken into account as well. In the quark sector of the SM the relation between the interaction and the mass eigenstates is encoded in terms of the CKM matrix [5,6]. The additional states in models of physics beyond the SM (BSM) usually follow the same pattern, i.e. interaction eigenstates mix with each other to form mass eigenstates, and the propagators of the latter are of matrix-type as a consequence of mixing contributions at higher orders (in some cases it can also be convenient to derive loop-corrected mass eigenstates directly from lowest-order interaction eigenstates, avoiding the introduction of lowest-order mass eigenstates). In particular, extended Higgs sectors usually contain several neutral Higgs states that can mix with each other. If CP-conservation is assumed, CP-even and CP-odd states can only mix among themselves. In the general case where the possibility of CP-violating interactions affecting the Higgs sector is taken into account, non-vanishing mixing contributions can occur between all neutral Higgs bosons of the model. Examples are general two-Higgs-doublet models, the Higgs sector of the minimal supersymmetric extension of the SM (MSSM), as well as the Higgs sectors of non-minimal supersymmetric models. Similarly, propagator mixing between vector resonances can occur for instance for the Kaluza-Klein excitations of the Z-boson and the photon (Z and A , respectively), which can have similar masses and have a large mixing with each other [7]. New, possibly nearly degenerate vector bosons can also arise in extended gauge symmetries and theories with a strongly coupled sector, e.g. composite Higgs models. For the case of the mixing of two Z -like vectors, see e.g. ref. [8]. An example involving mixing with a new fermion is the mixing between the top quark t and a top partner t [9].
In order to properly treat external particles in physical processes, the correct on-shell properties of in-and outgoing particles have to be ensured. While the S-matrix theory is formulated only for stable in-and outgoing states [10], in practical applications one often JHEP09(2017)079 needs to deal with the case of unstable external particles. The proper normalisation of external particles requires in particular that all external particles are on their mass shell, that their residues are equal to one, and that the mixing contributions with all other particles vanish on-shell at the considered order. This requirement holds independently from the specific renormalisation scheme that has been adopted for a certain calculation. In case the renormalisation conditions of the chosen scheme do not impose the proper normalisation of external particles, the correct on-shell properties need to be established via UV-finite wave function normalisation factors ("Z-factors"). In the case of unstable particles, finite-width effects have to be taken into account, and more generally imaginary parts arising from the absorptive parts of loop integrals need to be properly treated. A further complication arises if CP-violating effects are taken into account, since in this case the presence of complex parameters yields an additional source of imaginary parts. In the SM the different sources of imaginary parts affect for instance the renormalisation of the CKM matrix from the two-loop level onwards, see e.g. the discussion in ref. [11] and references therein. In the MSSM with complex parameters the corresponding effects arising from absorptive parts of loop integrals and complex model parameters enter predictions for physical observables in the chargino and neutralino sector already at the one-loop level, see refs. [12][13][14].
In the present paper we demonstrate that on-shell wave function normalisation factors, evaluated at the complex poles, are well suited for approximating the full mixing propagators also off-shell. In order to achieve this, the elements of the propagator matrix are expressed in terms of simple Breit-Wigner propagators that are multiplied with wave function normalisation factors. We show for various examples that this approximation works very well in practice. In this context we discuss in detail the uncertainties that are associated with the neglected terms in the approximation and provide estimates of their impact, which we find to be in good agreement with the observed numerical deviations between the full result and the approximation. We furthermore demonstrate that besides very significant technical simplifications our approach also has conceptual advantages.
Concerning technical simplifications, the evaluation of the propagator matrix for a given value of the external momentum squared requires the full momentum dependence of all contributing diagonal and off-diagonal self-energies. In our approximation based on wave function normalisation factors that are evaluated at the complex poles, on the other hand, the only momentum dependence that needs to be taken into account is the one of the simple Breit-Wigner propagators. We emphasize that our approach goes significantly beyond a simple "effective coupling"-type approximation where all external momenta at the loop level are neglected. While the latter approximation would disregard the important absorptive parts of the loop functions, those contributions are taken account in our approach through the evaluation of the Z-factors at the complex poles. We will demonstrate in particular that expanding the full propagators around all of their complex poles indeed results in the sum of Breit-Wigner propagators in combination with on-shell Z-factors. Approximating the full propagator matrix in this simple and convenient way has several conceptual advantages. In particular, the formulation in terms of Breit-Wigner propagators and Z-factors facilitates the implementation of a more precise total width than the JHEP09(2017)079 width obtained from the imaginary part of the complex pole with self-energies evaluated at the same order. Furthermore, for states with masses that are nearly degenerate (i.e. the difference between the masses of the two states is of the same order as the sum of the total widths of the two states) the mixing becomes resonantly enhanced and the resonances overlap, see e.g. refs. [7][8][9][15][16][17]. In such a case a single-pole approximation is not applicable [7,8,16]. Instead, the contributions from the relevant complex poles of the propagators have to be taken into account. The approximation of the propagator matrix in terms of on-shell wave function normalisation factors evaluated at the complex poles and Breit-Wigner propagators allows the implementation of interference contributions using the prescription of ref. [16].
We perform in this paper a careful treatment of imaginary parts and analyse the pole structure of matrix-type propagators in detail. While in the no-mixing case every element of an (n × n) propagator matrix has a single pole, we will demonstrate that in general every single element of an (n × n) propagator matrix (diagonal and off-diagonal) by itself has n poles. This fact makes the association of the (loop-corrected) mass eigenstates with the propagator poles non-trivial. We will show that in principle all permutations of the n propagator poles with the n mass eigenstates are physically equivalent. For practical purposes, however, it is important to make a choice that leads to numerically stable and well-behaved results.
To be specific, in the following we use the language of the Higgs sector of the MSSM with complex parameters. In this case CP-violating effects entering at the loop level give rise to a mixing between the lowest-order mass eigenstates h, H, A, which are CP-eigenstates (h, H are CP-even states, while A is a CP-odd state). The (3×3) mixing between the neutral Higgs states is the smallest type of propagator mixing that exhibits the qualitative features of general (n × n) propagator matrices. Furthermore, the MSSM example is well-suited to emphasise the distinction between interaction eigenstates (in this case the component field of the two Higgs doublets, e.g. φ 1 and φ 2 ), lowest-order mass eigenstates (h, H, A) and loop-corrected mass eigenstates, which in the MSSM case are called h 1 , h 2 , h 3 . The MSSM case also illustrates the possibility of mixing with other spin states (γ, Z) and with unphysical degrees of freedom (G). An important feature of models like the MSSM is the fact that the masses are predicted from the input parameters of the model, which makes the need for an appropriate assignment of mass eigenstates with propagator poles particularly apparent. The generalisation of our results to other models, cases with (n × n) mixing for n > 3 and other spin states should be straightforward.
The paper is organised as follows. After reviewing aspects of higher-order contributions and mixing effects in a system of three scalar states in section 2, we derive the pole structure of the propagator matrix and resulting relations between the wave function normalisation factors in section 3. An analytical derivation of the multi-Breit-Wigner approximation of the full propagators including an uncertainty estimate is given in section 4. We perform a detailed numerical comparison of the two approaches in section 5 before we conclude in section 6.

JHEP09(2017)079
2 Loop-level mixing of scalar propagators In a system of n particles which have the same conserved quantum numbers in general a non-trivial mixing between the different states will occur. While at lowest order the mass matrix given in terms of the interaction eigenstates can be diagonalised in order to obtain the mass eigenstates, at the loop level momentum-dependent self-energy corrections enter. These contributions give rise to a momentum-dependent (n × n) propagator matrix which is in general non-diagonal. The physical masses that are associated with the loop-corrected propagator matrix need to be determined from the (complex) poles of the propagator matrix.
While the case of (2 × 2) mixing provides the minimal setup for studying such a propagator structure, the issue of how to associate the mass eigenstates with the n poles of the propagator matrix becomes fully apparent only from the (3 × 3) case onwards. The case of (3 × 3) mixing captures all relevant aspects of the pole structure and permutation of the states so that the conclusions can be generalised for larger mixing systems.
For illustration and definiteness we use in this paper the language and the formalism of the neutral Higgs sector of the MSSM with CP-violating mixing, i.e. for the three lowestorder mass eigenstates i, j, k = h, H, A and their mixing into the loop-corrected mass eigenstates h a , a = 1, 2, 3. Since the analytic discussions in sections 2-4 do not rely on any model-dependent relations of the MSSM, the results can be transferred also to scalar sectors of other models by replacing h, H, A by a different set of i, j, k, and analogously h a by a general X a , where the generalisation to the case of mixing among more than three particles should be straightforward. Upon incorporation of the appropriate spin structure, the results can also be generalised to the case of mixing propagators of vector bosons or fermions.
The Higgs propagators in the MSSM in fact receive contributions from the mixing with the longitudinal components of the neutral gauge bosons Z and γ. While the Zboson propagator contributes to the diagonal elements of the Higgs propagator matrix at the two-loop level, the photon enters in those elements only at the four-loop level. Besides the mixing with the longitudinal component of the Z boson, also the mixing with the unphysical (would-be) Goldstone boson G needs to be taken into account. Instead of considering a (6 × 6) propagator matrix of the states {h, H, A, G, Z, γ}, we focus on the full mixing contributions of the physical Higgs fields, i.e. we treat the case of a (3 × 3) propagator matrix (which would generalise to (n × n) for the case of n physical scalar fields). For the (3 × 3) propagator matrix of physical Higgs fields, higher-order effects from the inversion of the matrix of the irreducible 2-point vertex functions are taken into account. In contrast, we treat the mixing contributions with the gauge bosons and the (would-be) Goldstone boson in a strict perturbative expansion up to the desired order. It should be noted that Higgs-G/Z mixing contributions already appear at the one-loop level in processes with external Higgs bosons, see e.g. refs. [12,14,18,19], and therefore need to be included in order to obtain a complete one-loop result for processes of this kind. For contributions to the propagator matrix with an incoming and outgoing Higgs boson the mixing contributions with G and Z enter from the order of (Σ i,G/Z ) 2 onwards, which is of JHEP09(2017)079 subleading two-loop order. For our numerical analysis below, see section 5, we will use the program FeynHiggs [20][21][22][23], where dominant two-loop corrections and leading higher-order contributions are incorporated in the irreducible self-energies. Since sub-leading two-loop contributions that are of the same order as the mixing contributions with G and Z [15] are neglected in this analysis, we include no further corrections beyond the (3 × 3) propagator matrix of the physical Higgs fields in this case.

Propagator matrix and the effective self-energy
As explained above, we use the case of the MSSM with complex parameters in order to illustrate propagator mixing between three physical scalar fields. If CP were assumed to be conserved in the MSSM, the CP-violating self-energies would vanish,Σ hA =Σ HA = 0, so that only the two CP-even states h and H would mix with each other. Our treatment corresponds to the case where non-zero phases from complex parameters are taken into account. Hence all renormalised self-energiesΣ ij (p 2 ) of the Higgs bosons i, j = h, H, A are in general non-vanishing, so that the matrix M of mass squares consists of the tree-level masses m 2 i on the diagonal and renormalised self-energies on the diagonal and off-diagonal entries. Expressed in terms of the lowest-order mass eigenstates h, H, A, which are also CP eigenstates, the matrix takes the form (2.1) The renormalised irreducible 2-point vertex functionŝ can be collected in the 3 × 3 matrixΓ hHA in terms of M aŝ Finally, the propagator matrix ∆ hHA equals, up to the sign, the inverse ofΓ hHA , Accordingly, the matrix inversion yields the individual propagators ∆ ij (p 2 ) as the the (ij) elements of the 3 × 3 matrix ∆ hHA (p 2 ), (2.5)

JHEP09(2017)079
The off-diagonal entries (for i = j) result in: (2.6) All 2-point vertex functionsΓ(p 2 ) depend on p 2 via eq. (2.2). Here we do not write the p 2 -dependence explicitly for the purpose of a simpler notation, but the full p 2 -dependence is implied also below. The solutions of the diagonal propagators, ∆ ii , can be expressed in the following compact way: , where the effective self-energy is introduced, . (2.9) It contains the diagonal self-energy,Σ ii (which exists already at 1-loop order), and the mixing 2-point functions (whose products only contribute toΣ eff ii from 2-loop order onwards). Hence, replacing the pure self-energyΣ ii by the effective one,Σ eff ii , includes also the 3 × 3 mixing contributions to the diagonal propagator in eq. (2.8) while preserving formally the structure of the propagator as in the unmixed case. In the limit of no mixing, the second term in eq. (2.9) vanishes.
The composition ofΣ eff ii in terms of the unmixed and the mixing contributions in eq. (2.9) can be written in the alternative way [13,30,31] for i, j, k all different. eq. (2.10) represents the sum of the diagonal and off-diagonal selfenergies involving a Higgs boson i where the off-diagonal contributions are weighted by the ratios of the respective off-diagonal and the diagonal propagators. This expression follows from eq. (2.9) with the help of the equality where the shorthand has been used, analogously for j ↔ k, andΣ ij = −iΓ ij from eq. (2.2) with i = j (for the special case of (2 × 2) mixing, this relation can be directly read off from eqs. (3.9)-(3.11) given below).

Complex poles
The complex poles M 2 a , where a = 1, 2, 3, are determined as solutions of the equations with i = h, H, A, see section 3 below. While the propagators have poles for p 2 = M 2 a , it is interesting to note that the ratios of off-diagonal and diagonal propagators (2.14) stay finite at the complex poles M 2 a , a = 1, 2, 3.

2.2.2Ẑ-matrix for on-shell properties of external particles
Higgs bosons appearing as external particles in a process need the appropriate on-shell properties for a correct normalisation of the S-matrix. In the hybrid on-shell/DR renormalisation scheme [15], the masses are renormalised on-shell, but the DR renormalisation conditions for the fields and tan β [15], where α is the mixing angle between h and H, see eq. (5.7) below, do not ensure proper onshell properties of the Higgs bosons. In fact, the loop-corrected mass eigenstates h 1 , h 2 , h 3 , which occur as external, on-shell, particles e.g. in decay processes, are a mixture of the lowest order states h, H, A. Thus, finite wave function normalisation factors need to be introduced in order to guarantee that the mixing vanishes on-shell and that the propagators of the external particles have unit residue. These so-called Z-factors for a neutral Higgs boson i = h, H, A on an external line are obtained from the residue of the propagators at the complex pole M 2 a , a = 1, 2, 3, [25,27] ExpandingΣ eff ii (p 2 ) around the complex pole p 2 = M 2 a , one obtains the diagonal propagator . (2.19) Performing the limit p 2 → M 2 a yields the residue of eq. (2.18), . (2.20)

JHEP09(2017)079
Considering a diagram with the Higgs boson i on an external line, whose propagator has three poles (see section 3), there are three possibilities which residue to compute. If the amputated Green's function is evaluated at M 2 a , the external i-line has to be multiplied by Ẑa i for the correct S-matrix normalisation. 1 So the resulting mass eigenstate as an outgoing particle is h a . Alternatively, if the Green's function is evaluated at M 2 b , it has to be normalised by to achieve the correct S-matrix element. For this choice, the external mass eigenstate is h b . Moreover, the wave function normalisation factor for i − j mixing on an external onshell line at M 2 a is composed of the overall normalisation factor Ẑa i times the on-shell transition ratioẐ (2.22) Since ∆ ii and ∆ ij have -in the case of 3 × 3 mixing -3 complex poles (see section 3), any of them can be chosen for the evaluation ofẐ i ,Ẑ ij andẐ ik ; in the considered example this is M 2 a . Correspondingly,Ẑ j ,Ẑ ji andẐ jk will be evaluated at M 2 b andẐ k ,Ẑ ki andẐ kj at M 2 c where a, b, c are a permutation of 1,2,3, and i, j, k a permutation of h, H, A [18,31]. For 2 × 2 mixing, of course only the two indices involved in the mixing can be permuted.
All choices allowed by the mixing structure are generally possible because of the pole structure of each propagator. However, they might not be equally numerically stable. If the loop-corrected mass eigenstate h a contains only a small admixture of the lowest-order state i, the propagator still has a pole at M 2 a , but the contribution of i to h a is suppressed at p 2 = M 2 a . In order to be definite, it is in any case necessary to define at which pole to evaluate which normalisation and mixingẐ-factor. This choice corresponds to fixing an assignment between the indices i, j, k of the lowest-order states and the indices a, b, c of the higherorder mixed states and then using it consistently. The assignment (i, a), (j, b), (k, c), which we label as scheme I, prescribes to evaluateẐ i ,Ẑ ij andẐ ik at M 2 a . Once the indices have been assigned we can clear up the notation by writinĝ and accordingly for the other indices such that the first index always refers to a loopcorrected mass eigenstate (a, b, c ∈ {1, 2, 3}) 2 and the second index to a lowest-order state (i, j, k ∈ {h, H, A}). Note thatẐ ai =Ẑ bj =Ẑ ck ≡ 1 in the index scheme I defined above. Once the index scheme has been specified, one can leave out the subscript I. For the scheme-independence of physical results, see section 3.2.

JHEP09(2017)079
Furthermore, it is convenient [15] to arrange the products of the normalisation factors Ẑ a and transition ratiosẐ aj asẐ aj = Ẑ aẐaj (2.24) (note the difference betweenẐ aj andẐ aj ) into a non-unitary matrix: TheẐ-matrix defined in eq. (2.25) fulfils the on-shell conditions (unit residue and vanishing mixing), which can be written in the following compact form [13,18,19,30]: with vanishing off-diagonal entries. It is equally possible to begin with these equations (2.26)-(2.28), i.e. to require unit residues and to demand vanishing mixing on-shell, to derive the elements of theẐ-matrix whose solutions are given in eqs. (2.20) and (2.22). In ref. [18] the evaluation at the full complex poles and the inclusion of imaginary parts were introduced (see eq. (2.35) below for the treatment of loop functions of complex arguments). The evaluation at the complex poles leads to numerically more stable results than the evaluation at real p 2 = M 2 ha , as well as to the physically equivalent choices of index assignments. The calculation of theẐ-factors of MSSM Higgs bosons can be performed with the program FeynHiggs.

Use ofẐ-factors for external states
TheẐ-factors have been introduced for the correct normalisation of matrix elements with external (on-shell) Higgs bosons h a , p 2 = M 2 a . It should be noted thatẐ does not provide a unitary transformation between the basis of lowest-order states and the basis of loopcorrected mass eigenstates. The fact thatẐ is a non-unitary matrix is related to the imaginary parts appearing in the propagators of unstable particles. Using theẐ-matrix, one can express the one-particle irreducible (1PI) vertex functionsΓ ha involving a loopcorrected mass eigenstate h 1 , h 2 , h 3 as an external particle as a linear combination of the 1PI vertex functions of the lowest-order states,Γ i : = Ẑ a Ẑ ahΓh +Ẑ aHΓH +Ẑ aAΓA + . . . , (2.30) JHEP09(2017)079 where the ellipsis refers to additional terms arising from the mixing with Goldstone and vector bosons, which are not described by theẐ-matrix. Thus, the overall normalisation factor Ẑ a accounts for the particle h a appearing at an external line. In addition, the factorsẐ ai given in eqs. (2.22) and (2.23) as ratios of propagators at p 2 = M 2 a describe the transition between the states h a and i. The transition factorẐ ai occurs in a diagram where h a is the external particle, but i directly couples to the vertex. All possibilities for i = h, H, A need to be included for each h a , hence the sum arises. This is depicted in figure 1 (cf. also refs. [30,32]).
Conveniently, eq. (2.30) can be written in matrix form for all h 1 , h 2 , h 3 as In this way, propagator corrections at external legs are effectively absorbed into the vertices of neutral Higgs bosons. IfẐ-factors are applied to supplement the Born result, only non-Higgs propagator type corrections (such as mixing with the Goldstone and Z-bosons) as well as vertex, box and real corrections need to be calculated individually.
As it is the case for the usual (UV-divergent) on-shell field renormalisation constants, the finite wave function normalisation factors are in general gauge-dependent quantities. While the gauge-parameter dependence drops out at a fixed order in perturbation theory, the inclusion of higher-order contributions into theẐ-factors described above can give rise to a residual gauge-parameter dependence, which is formally of sub-leading higher order. This is caused by the fact that theẐ-factors contain not only gauge-parameter independent leading higher-order contributions but also sub-leading higher-order contributions that are in general gauge-parameter dependent. Those sub-leading higher-order contributions would cancel with corresponding sub-leading higher-order contributions from other Green functions such as vertex and box contributions.

Effective couplings
Since theẐ-matrix is not unitary, it does not represent a unitary transformation between the {h, H, A} and the {h 1 , h 2 , h 3 } basis. However, it is not necessary to diagonalise the mass matrix for the determination of the poles of the propagators. Hence there is a priori no need to introduce a unitary transformation. Though, if a unitary matrix U is desired for the definition of effective couplings, an approximation of the momentum dependence ofẐ is required. There is no unique prescription of how to achieve a unitary mixing matrix as an JHEP09(2017)079 approximation of theẐ-matrix, but a possible choice is the p 2 = 0 approximation [15,32]. As in the effective potential approach, the external momentum p 2 is set to zero in the renormalised self-energiesΣ ij (p 2 ) →Σ ij (0) so that they become real. Then U diagonalises the real matrix M(0). U can be chosen real and it transforms the CP-eigenstates into the mass eigenstates, so that U 2 aA quantifies the admixture of a CP-odd component inside h a [15]. The elements of U can then be used to introduce effective couplings of the loop-corrected states h a to any other particles X in terms of the couplings of the unmixed states i by the relation absorbing some higher-order corrections, but neglecting imaginary parts and the full momentum dependence of the self-energies. Hence, the application of U resembles the use of Z-factors in eq. (2.29) for the purpose of implementing partial higher-order effects into an improved Born result. Yet, the rotation matrix U introduced for effective couplings as a unitary approximation is conceptually different from theẐ-matrix arising from propagator corrections and introduced for the correct normalisation of the S-matrix.

Technical treatment of imaginary parts
In order to account for complex momenta and imaginary parts of self-energies, we employ an expansion of the self-energies around the real part of the complex momentum, This expansion enables the calculation of self-energies at complex momentum in terms of self-energies evaluated at real momentum (e.g. by FeynHiggs). For the inclusion of all products of real and imaginary parts, we do not expand the effective self-energy from eq. (2.9) directly according to eq. (2.34). Instead, in the same way as in refs. [13,30,31], we expand allΓ ij (p 2 ) in eq. (2.9) individually before combining them intoΣ eff ii .

Pole structure of propagator matrices
In this section, we will first discuss in detail the solutions of the poles of the propagator matrix depending on the mixing. Subsequently we will prove the equivalence of different assignments between the lowest-order states and the loop-corrected mass eigenstates for the physicalẐ-matrix.

Determination of the poles
Imaginary parts of the self-energies lead to propagator poles in the complex momentum plane. The Higgs masses are determined from the real part of the complex poles M 2 of the diagonal propagators. Equivalently, the complex poles are obtained as the zeros of the inverse diagonal propagators. Due to because of the relation between ∆ hHA andΓ hHA from eq. (2.4). 3 Then the loop-corrected masses M are obtained from the real parts of the complex poles and the total widths Γ from the imaginary parts via In the following, the impact of higher-order and mixing contributions on the pole structure of the propagators will be discussed.
Lowest order. At lowest order, the self-energy contributions in eq. (2.2) are absent and the matrixΓ simply readŝ where the shorthand of eq. (2.12) has been used. The solutions of eq. (3.1) are the three tree-level masses m 2 i .
Higher order without mixing. Beyond tree-level, the self-energies are added at the available order. Restricting them to the unmixed case,Σ ij = 0 for i = j, leads tô is achieved if p 2 fulfils the following on-shell relation for any i = h, H, A. Thus, the full propagator matrix ∆ has three poles and each propagator ∆ ii (p 2 ) has exactly one pole p 2 = M 2 i that solves eq. (3.6) in this case.

JHEP09(2017)079
Higher order with 2 × 2 mixing. If now the mixing between h and H is taken into account, corresponding to the CP-conserving case, the matrixΓ becomes block-diagonal with the 2 × 2 matrixΓ hH and the 2-point vertex function of A, which does not mix with the other states:Γ For a closer look at the relation between the roots of the determinant and the roots of the inverse propagator, we write down the propagators and the effective self-energy of the {h, H} system explicitly. They follow from eqs. (2.6), (2.7) and (2.9) by settingΣ hA = Σ HA = 0 or equivalently from the inversion of the 2 × 2 submatrixΓ hH : .
Comparing the inverse diagonal propagators with the determinant of the submatrixΓ hH , we find for i, j ∈ {h, H} , i = j, Eq. (3.12) reveals that both inverse diagonal propagators, 1/∆ hh and 1/∆ HH , are proportional to the determinant ofΓ hH , which has two zeros. As opposed to the unmixed case, both zeros of det Γ hH (p 2 ) are poles of each of the diagonal propagators ∆ hh , ∆ HH . The lowest-order states h and H are mixed into the loop-corrected mass eigenstates h 1 and h 2 with the loop-corrected masses M h 1 , M h 2 . The corresponding poles for p 2 = M 2 ha in any combination of i = h, H and a = 1, 2, where the effective self-energy is given in eq. (3.11). In the 2 × 2 mixing system, it is convenient to choose M h 1 ≤ M h 2 . As for the nomenclature in the 2 × 2 case, the lighter mass eigenstate h 1 is often denoted as h and the heavier one as H because both are CP-even states. It should be noted that both roots ofΓ hH , M 2 h 1 and M 2 h 2 , are also complex poles of the off-diagonal propagators ∆ hH (p 2 ) ≡ ∆ Hh (p 2 ) due to (3.14)

JHEP09(2017)079
Since in this case A does not mix with h and H, the third pole M 2 A solely solves (3.15) but no other combination of A and h a satisfies the on-shell condition. M A is the loopcorrected mass of the (mass and interaction) eigenstate A.
Higher order with 3 × 3 mixing. Now we turn to the case where complex MSSM parameters lead to CP-violating self-energiesΣ hA ,Σ HA . Thus, all three neutral Higgs lowestorder and CP-eigenstates h, H, A mix into the loop-corrected mass eigenstates h 1 , h 2 , h 3 , which have no longer well-defined CP quantum numbers, but are admixtures of CP-even and CP-odd components. In this framework,Γ hHA is a full 3 × 3 matrix with the determinant where we dropped the explicit p 2 -dependence of each term for an ease of notation. Comparing eq. (3.16) with the diagonal and off-diagonal propagators from eqs. (2.7) and (2.6), respectively, we see that their inverse is proportional to the determinant ofΓ hHA : where i = j = k = i and without summation over the indices. From eq. (3.17) we conclude that all three roots p 2 = M 2 ha , a = 1, 2, 3, of det Γ hHA (p 2 ) are complex poles of each of the three diagonal propagators ∆ ii , i = h, H, A. This means that holds for any combination of i and a in the presence of 3 × 3 mixing. Moreover, eq. (3.18) implies that also the off-diagonal propagators have as many poles as the determinant has zeros, namely three in the case of CP-violating mixing. In the unmixed case,Σ eff ii =Σ ii and each propagator has exactly one pole so that there is a unique mapping between i and a, see eq. (3.6). On the other hand, for the general mixing case it is not unique how to relate the mass eigenstates to the interaction eigenstates. An assignment will be needed for the definition of on-shell wave function normalisation factors in section 2.2.

Index scheme independence of theẐ-matrix
As discussed above, the pole structure of the full propagators provides the freedom at which pole to evaluate whichẐ-factor, or equivalently which loop-corrected mass eigenstate index JHEP09(2017)079 (a, b, c) to assign to a lowest-order mass eigenstate index (i, j, k). We denote two choices for such index schemes by This initial ambiguity, however, results in physically equivalent results. Using properties of ratios of diagonal and off-diagonal propagators and exploiting relations at a complex pole (for details, see ref. [31]), we are able to show that This equality provides a transformation between scheme I (where i and a are associated, henceẐ ai ≡ 1) and scheme II (where j and a are matched): While the values ofẐ a andẐ ai depend on the choice of the index mapping, eq. (3.23) ensures that the elements of theẐ-matrix, which appear in physical processes, are invariant under the choice of a, b, c as a permutation of 1, 2, 3. We also tested this relation numerically for various parameter points and always found agreement within the numerical precision.

Breit-Wigner approximation of the full propagators
The propagator matrix depends on the momentum p 2 in a twofold way. On the one hand, the tree-level propagator factors D i (p 2 ) = p 2 − m 2 i give rise to an explicit p 2 -dependence. On the other hand, the self-energiesΣ ij (p 2 ) depend on the momentum as well, but away from thresholds their p 2 dependence is not particularly pronounced. In this section, we will develop an approximation of the full mixing propagators with the aim to maintain the leading momentum dependence, but to greatly simplify the mixing contributions by making use of the on-shellẐ-factors.

Unstable particles and the total decay width
In the context of determining complex poles of propagators, we now briefly discuss resonances and unstable particles, see e.g. refs. [33][34][35][36][37]. Stable particles are associated with a real pole of the S-matrix, whereas the self-energies of unstable particles develop an imaginary part, so that the pole of the propagator is located within the complex momentum plane off the real axis. For a single pole, the scattering amplitude A can be schematically written near the complex pole M 2 a in a gauge-invariant way as

JHEP09(2017)079
where s is the squared centre-of-mass energy, R denotes the residue, while F represents non-resonant contributions. The mass M ha of the unstable particle h a is obtained from the real part of the complex pole M 2 a = M 2 ha − iM ha Γ ha , while the imaginary part gives rise to the total width. Accordingly, the expansion around the complex pole M 2 a leads to a Breit-Wigner propagator with a constant width, In the following, we will use a Breit-Wigner propagator of this form to describe the contribution of the unstable scalar h a with mass M ha and total width Γ ha in the resonance region.

Expansion of the full propagators around the complex poles
Eqs. (3.17) and (3.18) imply for 3 × 3 mixing that each propagator ∆ ii , ∆ ij has a pole at M 2 1 , M 2 2 and M 2 3 . Because of this structure, an expansion of the full propagators near one single pole is not expected to yield a sufficient approximation. Instead, we will perform an expansion of the full propagators around all of their complex poles. The final expression obtained from combining the contributions from the different poles will constitute a main result of the present paper.

Expansion of the diagonal propagators
We begin with an expansion of ∆ ii (p 2 ) in the vicinity of M 2 a as in eq. (2.19), where the first factor equals the definition of the Breit-Wigner propagator of the state h a , and the second factor corresponds toẐ a in scheme I where i and a are associated indices. On top of that, Z a I =Ẑ 2 ai as defined in eq. (2.24), and the elements of theẐ-matrix are independent of the index scheme (see eq. (3.23)). Thus, the following scheme-independent approximation holds for p 2 M 2 a : In this approach, the mixing contributions are summarised in the on-shell Z-factor evaluated at M 2 a . In contrast, the leading momentum dependence is contained in the Breit-Wigner propagator parametrised by the loop-corrected mass M ha and the total width Γ ha from the complex pole. In addition, ∆ ii (p 2 ) has a second pole at M 2 b because p 2 − m 2 i +Σ eff ii (p 2 ) = 0 holds also at p 2 = M 2 b . Analogously, we can expandΣ eff ii around M 2 b and obtain for the diagonal propagator . (4.6)

JHEP09(2017)079
Formally, has the structure of the definition of aẐ-factor from eq. (2.20), but in the index scheme II where b is assigned to i, whereas eq. (4.3) has been obtained in scheme I with the (i, a) assignment. Using the relation (3.23), we can rewrite eq. (4.6) as where theẐ-factors in eq. (4.9) are expressed in the same scheme as in eq. (4.3). Hence, in the vicinity of p 2 M 2 b , the diagonal propagator ∆ ii can be approximated by the Breit-Wigner propagator of h b weighted by the square ofẐ bi that ensures the coupling to the incoming fields as Higgs boson i, propagation as the mass eigenstate h b and the coupling to the outgoing fields again as Higgs boson i. In the same manner, ∆ ii can be expanded around the third complex pole, M 2 c , yielding Thus, close to one of the complex poles (e.g. M 2 a ), the dominant contribution to the full propagator ∆ ii can be approximated by the corresponding Breit-Wigner propagator (∆ BW a ) multiplied by the square of the respectiveẐ-factor, (Ẑ 2 ai ). However, close-by poles may cause overlapping resonances. In order to include this possibility and to extend the range of validity of the Breit-Wigner approximation to a more general case, we take the sum of all three Breit-Wigner contributions into account: (4.14)

Expansion of the off-diagonal propagators
We proceed similarly for the off-diagonal propagators, which also have three complex poles so that we can expand the propagators around them. Note thatẐ ai = Ẑ a andẐ aj = Ẑ aẐaj as defined in eq. (2.24). Starting at p 2 M 2 a , we express theẐ-factors in scheme I,  Next, we approximate ∆ ij near p 2 = M 2 b : For p 2 M 2 c , we switch to a scheme where the indices i and c belong together. Thereby we can write which is expressed in terms of scheme-invariantẐ-factors. Finally, we take the sum of eqs. (4.15)-(4.17) to obtain This sum is illustrated diagrammatically in figure 2.
Eq. (4.18) represents the central result of this section, covering also the diagonal propagators in the special case of i = j. It shows how the full propagator can be approximated by the contributions of the three resonance regions, expressed by the Breit-Wigner propagators ∆ a (p 2 ), a = 1, 2, 3, reflecting the main momentum dependence. The mixing among the Higgs bosons is comprised in theẐ-factors which are evaluated on-shell. Nonetheless, even a part of the momentum dependence of the self-energies is accounted for because the derivation of eq. (4.14) is based on a first-order expansion of the momentum-dependent effective self-energies. Furthermore, theẐ-factors serve as transition factors between the loop-corrected mass eigenstates h a and the lowest-order states i (althoughẐ is not a unitary matrix transforming the states into each other). Pictorially, ∆ ij is a propagator that begins on the state i and ends on j while mixing occurs in between. Also on the r.h.s. of figure 2 and eq. (4.18) the propagator in the h a -basis begins with i and ends on j. Thus, the coupling to the rest of any diagram connected to the propagator is well-defined. In between, each of the h a can propagate, and the correct transition is ensured byẐ ai and Z aj . All three combinations are visualised in figure 2.
If CP is conserved and only h and H mix, or if two states are nearly degenerate and their resonances widely separated from the remaining complex pole, the full 3 × 3 mixing is (exactly or approximately) reduced to the 2 × 2 mixing case. Then the off-diagonal Z-factors involving the unmixed state vanish or become negligible so that some terms in eq. (4.18) approach zero.

JHEP09(2017)079
Beyond that, if no mixing occurs among the neutral Higgs bosons, all off-diagonal full propagators as well as the off-diagonalẐ-factors vanish and each diagonal full propagator consists of only a single Breit-Wigner term where theẐ-factor is based on the diagonal self-energy instead of the effective self-energy. Thus, eq. (4.18) covers all special cases of the a priori 3 × 3 mixing among the neutral Higgs bosons.

Uncertainty estimate
As for the usual Breit-Wigner approximation in the case of a single pole, the approximation of eq. (4.18) contains all resonant contributions for the case with mixing, while the nonresonant contributions that are neglected in this approximation are suppressed by powers of Γ ha /M ha . While a narrow width enhances the resonance contribution near its pole, it also causes it to fall rapidly for momenta away from the pole and thereby suppresses the nonresonant, far off-shell propagators, as well as their possible interference with a resonant term. The approximation in the off-shell regime is therefore of a comparable quality, determined by the ratios Γ ha /M ha for the different resonances, as the standard narrowwidth approximation without mixing. Interference effects between different resonances are accounted for by employing the approximation of eq. (4.18). As we will discuss in more detail below, interference effects can be large for resonances that are close in mass, while on the other hand interferences between widely separated resonances are suppressed. The effects of the interference of resonances with a non-resonant background continuum can be inferred from the case without mixing.
While the on-shell approximation of eq. (4.18) is based on the expansion around the complex pole, the description of the momentum dependence can be improved by incorporating higher powers of (p 2 −M 2 a ). It should be noted in this context that apart from thresholds the self-energies generally depend only rather mildly on the incoming momentum.
In the following, we will assess the main sources of uncertainties that apply in the vicinity of the complex poles. The approximation given in eq. (4.18) of the fully momentumdependent propagators in terms of Breit-Wigner propagators and on-shellẐ-factors introduces an uncertainty by neglecting higher orders in the expansion of the effective selfenergies and in the evaluation of the self-energies depending on complex momenta in eq. (2.35). Besides, the pole condition is numerically not exactly fulfilled, which may represent an additional source of uncertainty.
Expansion of self-energies around the real part of a complex momentum. As defined in eq. (2.35), we expand all self-energies as functions of a complex squared momentum p 2 around the real part of p 2 up to the first order in the imaginary part p 2 i ≡ Im p 2 . Furthermore, higher powers of p 2 i can enter, In order to expand the propagators in terms of the second-order contribution we define , (4.20)

JHEP09(2017)079
and calculate ∆ ii to first order in Im ij . One should keep in mind, though, that we also evaluate the full propagators using the above-mentioned expansion of self-energies around the real part of the momentum. Thus, the comparison between the full and the approximated propagators is not directly affected by omitting terms of O( Im ij ), but it represents a source of uncertainty in both methods.
Expansion of effective self-energies around a complex pole. Regarding the effective self-energies, the expansion ofΣ eff ii (p 2 ) around M 2 a in eq. (2.19) includes the term of O p 2 − M 2 a . Additional terms in the expansion up to O (p 2 − M 2 a ) 2 , contribute if the term involving the second derivative of the effective self-energy at the complex pole,Σ eff ii (M 2 a ), is non-negligible. We define where the approximation holds for ReΣ eff ii (M 2 a ), ImΣ eff ii (M 2 a ) 1 and the variation of p 2 = xM 2 a around M 2 a . The impact of the neglected terms on the propagator ∆ ii reads at first order in p a,i : The higher powers of (p 2 −M 2 a ) directly translate to an expansion in terms of Γ ha /M ha via Im M 2 a = −iM ha Γ ha , analogously to the unmixed case. Expanding the diagonal propagator ∆ ii (p 2 ) around the complex pole p 2 ∼ M 2 a yields Beyond that, the correction scales with powers of Γ ha /M ha multiplied by derivatives ofΣ eff ii . While the term of the first derivative giving rise to the diagonal Z-factor Z 2 ai is included in the approximation in eq. (4.14), the term linear in Γ ha /M ha times the second derivative, and higher-order terms, are omitted and therefore represent a source of uncertainty. These higher terms exist at the real part of the resonance, p 2 = M 2 ha , but vanish at the complex pole, p 2 = M 2 a . In any case, they are not only suppressed by a  Figure 3. Contributions from the full mixing propagators ∆ ij (p 2 ) for i, j = h, H, A to a generic amplitude (cf. ref. [13]). If theẐ-factor approach is applied, each of the 9 full propagators needs to be approximated by the sum of the three corresponding Breit-Wigner diagrams as shown in figure 2.
narrow width compared to the mass, but also by the smallness of the higher derivatives at the pole. We confirmed numerically that already the second derivatives are suppressed by several orders compared to the first derivatives. The expansion in eq. (4.24) indeed corresponds exactly to the well-known unmixed case for the simplification of M ha → M, Γ ha → Γ andΣ eff ii → Σ. In contrast, the presence of mixing necessitates the expansion around all complex poles. The mixing is accounted for by the effective self-energies, but the structure of the dependence on Γ ha /M ha remains in analogy to the unmixed case.
Even in the case of no mixing, the inclusion of the diagonalẐ-factor does not only represent a conceptual, but also a numerical improvement of the agreement between the full and the approximate propagator compared to the pure Breit-Wigner propagator.

Amplitude with mixing based on full or Breit-Wigner propagators
In a physical process where neutral Higgs bosons can appear as intermediate particles, all of them need to be included in the prediction, see figure 3 and ref. [13]. The Higgs part of the amplitude then contains a sum over the irreducible vertex functionsΓ X i (for a coupling of Higgs i at the first vertex X) andΓ Y j (for a coupling of Higgs j at the second vertex Y ) times the fully momentum-dependent mixing propagators, Applying eq. (4.18), the amplitude in eq. (4.26) can be approximated by the sum over Breit-Wigner propagators multiplied by on-shell Z-factors, in agreement with ref. [13], The first bracket in eq. (4.28) representsΓ X ha , i.e., the vertex X connected to the mass eigenstate h a as for an external Higgs boson in eq. (2.30). Subsequently, the second bracket is equal to the coupling of h a at vertex Y ,Γ Y ha . As opposed to section 2.3, the h a is not on-shell here, but a propagator with momentum p 2 between the vertices X and Y , represented by the Breit-Wigner propagator ∆ BW a (p 2 ). So theẐ-factors are not only useful for the on-shell properties of external Higgs bosons, but they can also be used as an onshell approximation of the mixing between Higgs propagators. This will be investigated numerically in section 5.
Concerning the propagators of unstable particles, the long-standing problems related to the treatment of unstable particles in quantum field theory should be kept in mind. While the introduction of a non-zero width into the propagator of an unstable particle is indispensable for treating the resonance region where the unstable particle is close to its mass shell, such a modification of the propagator in general mixes orders of perturbation theory. This can give rise to a violation of unitarity and / or to gauge-dependent results, see e.g. ref. [38] for a recent discussion of this issue and references therein. The approximation of the full off-shell propagators, where the effective self-energy appears in the denominator of the diagonal propagator according to eq. (2.8), in terms of on-shellẐ-factors and Breit-Wigner propagators is not meant to provide a solution to this long-standing problem. It should be noted, however, that the expression of the off-shell propagator in terms of on-shell Z-factors and simple Breit-Wigner factors significantly remedies the gauge dependence of the off-shell propagator. This is due to the fact that the approximation of eq. (4.3) and eq. (4.18) is on the one hand numerically close to the full propagator (in the considered kinematic region and applying the Feynman gauge) but on the other hand does not contain gauge-dependent contributions that depend on the squared momentum p 2 of the propagator. While those p 2 -dependent contributions of the gauge part of a Dyson-resummed self-energy in the denominator of a propagator often lead to violations of unitarity for high values of p 2 (see e.g. ref. [39]), the Breit-Wigner factors multiplied withẐ-factors that are evaluated on-shell are much better behaved in the limit where p 2 gets large. We defer a more detailed discussion of this issue to future work.

Calculation of interference terms in the Breit-Wigner formulation
In eq. (4.18), the Breit-Wigner propagators are combined such that they approximate a given full propagator. Conversely, we will now separate the h a part from the contribution JHEP09(2017)079 of the other mass eigenstates in the amplitude with Higgs exchange between the vertices X and Y : i.e. the exchange of the state h a coupling with the mixed verticesΓ ha from eq. (2.30) as for an external Higgs. In order to calculate the squared amplitude as a coherent sum, all contributions of h 1 , h 2 , h 3 are summed up first before taking the absolute square, (4.31) On the contrary, the incoherent sum is the sum of the squared individual amplitudes, which misses the interference contribution, (4.32) Thus, an advantage of the Breit-Wigner propagators is also the possibility to conveniently discern the interference of several resonances from their individual contributions in a squared amplitude

JHEP09(2017)079
In contrast, the squared amplitude based on the full propagators does not allow for a straightforward determination of the pure interference term.

Numerical comparison in the MSSM with complex parameters
After the analytical considerations so far, we will now compare the full propagators with their approximation as a combination of Breit-Wigner propagators andẐ-factors numerically. For definiteness, we will evaluate the propagators in the MSSM with complex parameters. We use FeynHiggs-2.10.2 4 [20][21][22][23] for the numerical evaluation of Higgs masses, widths and mixing properties. In order to investigate the applicability of the expansion of the full propagators in one or all three resonance regions, we will first use as input a complex squared momentum around the three complex poles. For the later application to physical processes where the squared momentum equals the centre-of-mass energy s, we will also evaluate the propagators at p 2 = s near the real parts of the complex poles. In section 5.2 we choose a scenario where all three Higgs bosons are relatively light so that we can study their mutual overlap. As a test of theẐ-factor approximation, we work in a scenario with large mixing between H and A in section 5.3.

The MSSM with complex parameters at tree level
In this section we fix the notation for the different particle sectors of the MSSM with complex parameters, following ref. [15].
Sfermion sector. The mixing of sfermionsf L ,f R within one generation into mass eigenstatesf 1 ,f 2 is parametrised by the matrix cot β, f = up-type tan β, f = down-type.
The trilinear couplings A f = |A f |e iφ A f , as well as µ = |µ|e iφµ , can be complex. These phases enter the Higgs sector via sfermion loops starting at one-loop order. Diagonalising M 2 f for allf separately, one obtains the sfermion masses mf 1 ≤ mf 2 . Gluino sector. The gluinog a , a = 1, 2, 3, has a mass of mg = |M 3 |, where M 3 = |M 3 | e iφ M 3 is the possibly complex gluino mass parameter. Since the gluino does not directly couple to Higgs bosons, the phase φ M 3 enters the Higgs sector only at the two-loop level, but has an impact for example on the bottom Yukawa coupling already at one-loop order.

JHEP09(2017)079
Neutralino and chargino sector. At tree-level, mixing in the chargino sector is governed by the higgsino and wino mass parameters µ and M 2 , respectively. In the neutralino sector it additionally depends on the bino mass parameter M 1 . The charginos χ ± i , i = 1, 2, as mass eigenstates are superpositions of the charged winos W ± and higgsinos H ± , with the mass matrix, Likewise in the neutralino sector, the neutral electroweak gauginos B, W 3 and the neutral Higgsinosh 0 d ,h 0 u mix into the mass eigenstates χ 0 i , i = 1, . . . , 4. The mixing is encoded in the gaugino mass matrix Y , The gaugino mass parameters M 1 and M 2 as well as the higgsino mass parameter can in principle be complex. However, only two of the phases are independent, and a frequently used convention is to set φ M 2 = 0.

(5.6)
A relative phase between both Higgs doublets vanishes at the minimum of the Higgs potential, and the possible phase of the coefficient of the bilinear term in the Higgs potential can be rotated away. Hence, the Higgs sector conserves CP at lowest order, and the 4 × 4 mass matrix of the neutral states M φφχχ becomes block-diagonal. The neutral tree level mass eigenstates are obtained by a diagonalisation of M φφχχ : where we introduced the short-hand notation s x ≡ sin x, c x ≡ cos x. For later use we define t β ≡ tan β. The mixing angle α is applied for the CP-even Higgs bosons h, H; β n for the neutral CP-odd Higgs A and Goldstone boson G, and β c for the charged Higgs H ± and the charged Goldstone boson G ± . The minimum conditions for V H lead to β = β n = β c at tree level. At higher orders, however, tan β is renormalised whereas the mixing angles α, β n

JHEP09(2017)079
and β c are not renormalised. The masses of the CP-odd and the charged Higgs bosons are at tree level related by At lowest order, the Higgs sector is fully determined by the two SUSY input parameters (in addition to SM masses and gauge couplings) tan β and m H ± (or, for conserved CP, equivalently m A ).
Higher-order corrections. Higher-order corrections in the MSSM Higgs sector are very relevant. Particles from other sectors contribute via loop diagrams to Higgs observables, in particular the trilinear couplings A f , the stop and sbottom masses, and in the sub-leading terms the higgsino mass parameter µ. Thus, beyond the lowest order, the Higgs sector is influenced by more parameters than only M H ± (or M A ) and tan β. We adopt the hybrid on-shell and DR-renormalisation scheme defined in ref. [15].

Scenario with three light Higgs bosons
For the numerical evaluation of the propagators, we work in the M mod+ h scenario [40]. In this example, we fix the variable parameters µ = 200 GeV, M H ± = 160 GeV, tan β = 50, (5.9) and choose for the complex phase φ At = π/4 to allow for CP-violating mixing. This parameter choice is not meant to be the experimentally viable; in fact it has been experimentally ruled out by searches for H, A → τ τ . The purpose is rather to provide a setting for illustration with nearby, but resolvable resonances so that the test of the Breit-Wigner approximation is not limited to well separated poles. These parameter values result in the following complex poles: is only significant if the productẐ aiẐaj is not suppressed. In this case, we can already estimate that, for example, h 3 hardly contributes to ∆ AA . We have analysed all propagators ∆ ij around M 2 1 , M 2 2 and M 2 3 as well as for real momenta. In the following, we show and discuss a selection of these cases.

Propagators depending on complex momenta
The analytical derivation of eq. (4.18) builds on the expansion of the full propagators around the complex poles, and the on-shell condition in eq. (3.19) holds exactly only at complex momentum. Therefore, we evaluate the self-energies and propagators around the complex poles. such that x 1 = 1 corresponds to the complex pole p 2 = M 2 1 . The black line (labelled by ∆ full) represents the fully momentum-dependent mixing propagator from eq. (2.8). Since the three poles do not have the same ratio between the real and imaginary parts, scaling x 1 does not run into M 2 2 and M 2 3 . ∆ hh has a pole at x = 1 and a second peak at x 1.1 which is close to the real part of M 2 3 . This structure is very precisely reproduced (for a discussion of the uncertainties, see below) by the sum 3 a=1Ẑ 2 ah ∆ BW a (p 2 ) according to eq. (4.14) -as can be seen from the red dotted line (labelled as BW · Z), which lies directly on top of the black solid line.
In order to understand which of the Breit-Wigner propagators andẐ-factors dominate at which momentum, the individual contributions are shown by the dashed curves. The blue line (labelled by h 1 ) represents the contribution of h 1 to ∆ hh , i.e.,Ẑ 2 1h ∆ BW 1 (p 2 ). It clearly displays the pole at x = 1, but strongly deviates from the full propagator at momenta away from M 2 1 . The orange line (labelled by h 2 ) representsẐ 2 2h ∆ BW 2 (p 2 ). Sincê Z 2h is small in this scenario, the contribution of h 2 to ∆ hh is numerically suppressed, but a tiny share is visible. The green line (labelled by h 3 ) stands forẐ 2 3h ∆ BW 3 (p 2 ) and it contributes significantly to ∆ hh near M 2 3 becauseẐ 3h = −0.35 − 0.09i is sizeable in this scenario. So we notice that none of the individual Breit-Wigner propagators multiplied by the appropriateẐ-factors suffices to approximate the full propagator, which has three complex poles. On the other hand, the sum of all three Breit-Wigner propagators timeŝ Z-factors yields an accurate approach to the full mixing. This holds for the real and the imaginary part.
Having discussed the example of a diagonal propagator, we will now assess whether thê Z-factor approximation succeeds also for off-diagonal propagators. For instance, figure 6 depicts ∆ HA versus x 2 = p 2 /M 2 2 such that x = 1 matches p 2 = M 2 2 where the propagator diverges. As above, owing to the different ratio between the real and imaginary part of each complex pole, scaling x 2 does not run into M 2 1 and M 2 3 , but ∆ HA peaks close to their real parts. As in figure 5 representing the sum of Breit-Wigner propagators according to eq. (4.18) agree very well. Additionally, one can see the individual Breit-Wigner shapes. Because the products of the relevantẐ-factors, hereẐ aHẐaA , are non-negligible for all a = 1, 2, 3, each Breit-Wigner propagator is important in the approximation of both the real part (figure 6(a)) and the imaginary part ( figure 6(b)) of ∆ HA . The other diagonal and off-diagonal propagators in this scenario which are not displayed here have an equally good agreement between the full calculation and the approximation.
In order to quantify this agreement, we compute the relative deviation of the approximated propagators from the full ones, for each diagonal and off-diagonal element of the propagator matrix, depending on the momentum around each of the complex poles. For p 2 = x · M 2 1 , 0.5 ≤ x ≤ 1.5 as in figure 5, the relative deviation of the approximated propagators from the full ones is below 4 · 10 −3 for ∆ hh , ∆ HH and below 9 · 10 −3 for ∆ AA (only ∆ hh is displayed here), apart from the regions where the propagator itself goes through zero and changes sign. A similar agreement holds both for the diagonal and off-diagonal propagators such as in figure 6, and for their real and imaginary parts. As expected, the relative deviation is minimal at the complex pole (x = 1). The largest deviation between the full and the approximated propagators occurs (apart from the sign flips) at momenta far away from the pole where the full propagators approach zero and in some cases the individual propagators of the states h a contribute with opposite signs to sum.
The relative deviation for the displayed momenta should be compared to the expected uncertainty based on the estimates in section 4.2.3. The higher powers of the imaginary part of the momentum that are neglected in the approximate evaluation of the self-energies at complex momenta give rise to only a small effect ( Im ij ∼ O(10 −7 ). This induces an effect on the propagator of O(10 −9 − 10 −6 )), which does not account for the actual difference (as explained above, both the full and the approximated propagators are subject to this kind of uncertainties arising from the evaluation at complex momenta). In contrast, the relative effect of the expansion of the momentum-dependent effective self-energy around the complex pole has a more significant effect of up to p a,i 6 · 10 −3 for the analysed momentum range, where the largest deviations occur furthest away from the poles. As a conclusion, the uncertainty arising from the approximate evaluation of the propagators at complex momenta is found to be negligible, while the expected uncertainty from the expansion around a complex pole is very close to the observed deviation.

Propagators depending on the real momentum p 2 = s
The calculation of the propagators at and around the complex poles together with the evaluation of the self-energies at complex momenta according to eq. (2.35) was needed to fulfil the assumptions of the approximation. However, in collider processes, the Higgs propagator might appear for example in the s-channel of a 2 → 2 scattering process where the squared momentum equals the square of the centre-of-mass energy s. Therefore here we will check the Breit-Wigner approximation around the real parts of the complex poles. Re [∆ AA ] with a broad peak at M h 2 . The black curve of the full propagator is again directly underneath the red, dotted curve of the Breit-Wigner approximation, which in this case stems nearly entirely from h 2 because ofẐ 2A 1. Within ∆ AA , the contribution of h 1 only has a minor impact, which can be seen as a small kink in figure 7(b). Although ∆ BW 1 (p 2 ) gets close to its pole, the resonance of h 1 is strongly suppressed by the smallẐ-factor Z 1A = 0.05(1 + i). As we anticipated above from the structure ofẐ in eq. (5.12), ∆ BW 3 is a negligible component of ∆ AA for this parameter point.
So far we have seen that the Breit-Wigner formulation combined with on-shellẐfactors accurately reproduces (up to a relative deviation of O(10 −3 −10 −2 )) the momentum dependence of the full diagonal and off-diagonal propagators by adding the contributions from all resonance regions. If all resonances are sufficiently separated (not shown in this example) or if all but one contributing products ofẐ-factors are negligible, a single Breit-Wigner term is enough to approximate the full propagator in one of the resonance regions. In the general case, however, all Breit-Wigner terms need to be included. Even if the peaks are not located very close to each other compared to their widths, the tail of one resonance supported by a substantial product ofẐ-factors can leak into another resonance region.

Scenario with large mixing
While the scenario in the previous section was characterised by three relatively similar masses, we now choose a setting with quasi degenerate heavy states h 2 and h 3 . In section 5.2 we considered the M mod+ h -scenario with the standard value of µ = 200 GeV in combination with the phase φ At = π/4, leading to a moderate mixing predominantly between h and A. In ref. [40] it was suggested to choose also different values, µ = ±200, ±500, ±1000 GeV. So in addition to the choice above, we now apply the following modification of the parameters in eq. (5.9): µ = 1000 GeV, M H ± = 650 GeV, tan β = 20. (5.14)

JHEP09(2017)079
This results in the complex poles Hence it is useful to examine whether the Breit-Wigner propagators withẐ-factors still yield a viable approximation of the fully momentum-dependent propagators in this scenario with large mixing.

Propagators depending on complex momenta
Due to the large difference between M h 1 and M h 2 M h 3 , the propagators ∆ hh , ∆ hH and ∆ hA are strongly suppressed around M 2 2 . Figure 8(a) shows Re [∆ HH ] for complex momentum around M 2 2 . The full calculation (black) and the Breit-Wigner approximation (red, dotted) are in relatively good agreement with each other although the curves do not lie directly on top of one another as in the scenario of section 5.2. The two heavy states h 2 and h 3 with a mass difference of less than 0.2 GeV and total widths of Γ h 2 = 2.6 GeV and Γ h 3 = 2.9 GeV are too close to be resolved.Ẑ 2 2H ∆ h 2 (orange) andẐ 2 3H ∆ h 3 (green) contribute with similar magnitude, but opposite signs so that the result differs strongly from the single terms. Indeed, only their sum provides a reliable approximation. A comparable situation is shown in figure 8 In this scenario of highly admixed, almost mass-degenerate states h 2 , h 3 , the relative deviation of the approximated propagators from the full ones is larger, in particular in the JHEP09(2017)079 region where the propagator itself is close to zero. The uncertainty of the evaluation of complex momenta is negligible here, as in the scenario of section 5.2. Instead, the reasons for the resulting deviation are the following. The neglected terms from the expansion around the complex pole (see eq. (4.22)) amount to p a,i ∼ O(10 −2 − 10 −1 ) for the states i = H, A involved in the mixing around the nearby complex poles M 2 a , a = 2, 3. The relative effect of the deviation is enhanced by the fact that the full propagator is close to zero. The approximation consists of the sum of two resonant contributions (h 2 , h 3 ), which are in this case several orders of magnitude larger than the full propagator and occur with opposite signs. Hence, despite a cancellation of several orders of magnitude the approximation yields the correct sign of the propagator and a good description of its behaviour and size. In contrast, a single-resonance approach would result in a prediction that would be off by several orders of magnitude. As a result, although slight deviations between the approximated propagator and the full one are visible in this scenario that is characterised by large mixing effects, the approximation based on eq. (4.18) still represents a very significant improvement compared to a single-resonance treatment. a strong resonance-enhanced mixing. We note, however, that the agreement between the full propagators and those with on-shell mixing factors is less accurate here than in the scenario with moderate, nearly unitary mixing.

Comparison ofẐ-factors with effective couplings
The effective coupling approach mentioned in section 2.4 makes use of the unitary, real Umatrix instead of theẐ-matrix. However, U is not evaluated at the complex pole, but at p 2 = 0 (see section 2.4) and it does not comprise the imaginary parts of the self-energies. We compare the effective coupling approach based on the matrix U with theẐ-factor approach which does take the imaginary parts into account, but cannot be directly interpreted as a unitary transformation between the states of the different bases. Based onẐ-and U-factors from FeynHiggs, figure 10 displays in the scenario of eq. (5.14) real parts of ∆ HH and ∆ HA at real momentum around M h 2 M h 3 , calculated as a fully momentum-dependent mixing propagator (black), using theẐ-matrix approach (red, dotted) defined in eq. (4.18), i.e. ∆ Z ij (p 2 ) 3 a=1Ẑ ai ∆ BW a (p 2 )Ẑ aj , and the U-matrix variant (grey, dashed) in the p 2 = 0 approximation as The black curves in figure 10 representing the full propagators are identical to those shown in figure 9. While in figure 10(a) the curve representing the approximation in terms ofẐ-factors is almost identical to the curve of the full ∆ HH (with a relative deviation at the peak of 0.8%), the approximation in terms of U-factors differs from the full result by up to 14%. In figure 10(b), not even the shape of ∆ HA is correctly approximated by the approach using the U-matrix whereas the approach using theẐ-matrix comes close to the full ∆ HA (up to the small deviations that can be seen in the plots). This analysis indicates that theẐ-factors combined with Breit-Wigner propagators are well-suited to describe the Higgs propagators including their mixing also in scenarios with close-by resonances and strong mixing. This approach captures the leading momentum dependence and adequately accounts for the imaginary parts. In contrast, the combination of U-factors and Breit-Wigner propagators is -despite its unitary nature -incomplete with respect to the mixing effects in the resonance region and regarding the significance of imaginary parts.

Breit-Wigner and full propagators in cross sections
As an application of the derivations above, we calculate a cross section with Higgs exchange. We study the example process bb → h, H, A → τ + τ − , where the intermediate Higgs bosons are once represented by the full mixing propagators ∆ ij and once by Breit-Wigner propagators multiplied byẐ-factors. In order to disentangle this investigation from other higher-order effects, we restrict the Higgs-fermion-fermion vertices to the tree-level and do not include the emission of real particles in the initial or final state, but focus on the propagator corrections.
For the implementation of the full propagator method, we extended a FeynArts model file. New scalars ij are introduced that correspond to the full propagator ∆ ij (p 2 ) and couple to the first vertex as the interaction eigenstate i and to the second vertex as j. Those propagators are used in the FormCalc calculation supplemented by self-energies from FeynHiggs incorporating corrections up to the two-loop level and the full momentum dependence at the one-loop level.
Considering only mixing between h and H at this point, we choose as a CP-conserving scenario the M max h -scenario [41,42] with tan β = 50, M H ± = 153 GeV, but we modify it by setting A f 3 = 2504 GeV. As before, this scenario has been selected for illustration purposes and it is not meant to be phenomenologically viable.   . The resonances of h 1 and h 2 partly overlap as the mass difference is of the order of the total widths, but the two peaks can still be distinguished. The h 3 contribution peaks at a lower mass in this scenario, but for completeness it is also shown (purple, dotted). The incoherent sum |h 1 | 2 +|h 2 | 2 +|h 3 | 2 (grey, dash-dotted) from eq. (4.32) clearly overestimates JHEP09(2017)079 the full cross section as a consequence of the missing interference term that turns out to be destructive in this case. It is taken into account in the full calculation and in the coherent sum of Breit-Wigner propagators withẐ-factors. For the efficient calculation of interference terms of quasi degenerate resonances, see e.g. refs. [16,43].
While the comparison in figure 11 is restricted to the case of 2 × 2-mixing due to a scenario with real parameters, the agreement of the cross sectionσ(bb → h, H, A → τ + τ − ) calculated with the full or Breit-Wigner propagators can be seen also in a scenario with complex parameters in figure 12. For the numerical evaluation, we choose a scenario with large mixing defined in eq. (5.14), i.e. the M mod+ h scenario with µ = 1000 GeV and the phase φ At = π/4 where we vary M H ± and t β .
We investigate the relative deviation between the cross section σ full based on the full propagators and the cross section σ BWẐ coh based on the coherent sum of Breit-Wigner propagators withẐ-factors, where the total widths are obtained from the imaginary parts of the complex poles, Γ Im , defined in eq. (5.23), (5.21) Figure 12 reveals that both methods agree very well, with a maximum deviation of ±2% around M H ± = 500 GeV, tan β = 28 and of about 0.8% along the green band from this parameter point to larger values of M H ± and lower values of t β . Otherwise the two calculations lead to the same results within 0.1%. Hence the use of Breit-Wigner propagators is suitable also for phenomenological applications such as the calculation of cross sections even in a scenario with complex parameters and large, CP-violating mixing. A phenomenological investigation of such a scenario will be addressed in a forthcoming publication [44][45][46].

Impact of the total width
This section addresses the impact of the precise value of the total width. So far, we have obtained the Higgs widths from the imaginary part of the complex poles as in eq. (3.2) in order to consistently compare with the full propagator mixing. If the self-energies inΣ eff ii are calculated at the one-loop level, the total width extracted from a complex pole of ∆ ii is then a tree-level width. Correspondingly, partial two-loop contributions to the imaginary parts of the self-energies give rise to partial one-loop corrections of the decay width. However, two-loop self-energies evaluated at p 2 = 0, as they are approximated in FeynHiggs, do not contribute to the imaginary part of the pole so that the width determined from the imaginary part of the complex pole remains at its tree-level value. Corrections to Higgs boson decays in the MSSM at and beyond the one-loop level are known and have been found to be important, see e.g. refs. [19,25,47,48]. Thus, the sum of the partial decay widths into any final state X of a Higgs boson h a , even if M 2 a is based on self-energies at the same order as used for the calculation of Γ tot ha as in eq. (5.22). FeynHiggs contains the partial Higgs decay widths and their sum at the leading two-loop order. Having checked the compelling agreement between the full propagators and the Breit-Wigner propagators with the width from the imaginary part of the complex pole in the previous sections, now we implement the total width from FeynHiggs into the Breit-Wigner propagators in order to obtain the most precise phenomenological prediction.
In the modified M max h scenario, the higher-order corrections have a significant impact on the Higgs decay widths so that Γ tot h 1 = 2.55 GeV and Γ tot h 2 = 3.24 GeV are much larger than the widths obtained from the imaginary part of the complex pole. This affects the order of magnitude of the cross-sectionσ bb → τ + τ − and the structure of the resonances, as can be seen in figure 13. The coherent sum of Breit-Wigner propagators including the interference term (red, dashed) and the incoherent sum without the interference term (grey, dash-dotted) using Γ Im from the imaginary parts of the complex poles are the same as in figure 11. In contrast, the total widths Γ tot FH obtained from FeynHiggs as the sum of higher-order partial widths are implemented into the Breit-Wigner propagators in the cross section based on the coherent sum of all h a -contributions (black, solid) and the incoherent sum (black, dotted).
The large increase in the widths from Γ Im ha to Γ tot ha has a very significant effect in figure 13, since the resonant behaviour is smeared out by the larger widths. As a consequence, the separate resonances are less pronounced, and the cross section is suppressed. Here, the incoherent sum without the interference term again overestimates the cross-sections. In addition, it lacks the two-peak structure. This observation emphasizes the importance of  figure 11 with tan β = 50 and M H + = 153 GeV. The Breit-Wigner propagators with the total widths from the imaginary part of the complex pole including (red, dashed) and excluding (grey, dash-dotted) the interference term (as in figure 11) are compared with the Breit-Wigner propagators where the total widths are obtained from FeynHiggs. The corresponding results are shown including (black, solid) and excluding (black, dotted) the interference term.
including the total width at the highest available precision and to take the interference term into account. One can also see that two resonances might be too close to be resolved if they are smeared by large widths.

Conclusions
We have shown that the momentum-dependent propagator matrices of systems of unstable particles that mix with each other can be accurately approximated by a combination of Breit-Wigner propagators and wave function normalisation factors, where the latter are evaluated at the complex poles of the propagator matrix. For illustration, we have applied this approach to the example of the neutral Higgs bosons of the MSSM with complex parameters, where the CP-violating interactions give rise to a 3 × 3 mixing between the lowest-order mass eigenstates h, H, A and the loop-corrected mass eigenstates h 1 , h 2 , h 3 (further mixing contributions with unphysical Goldstone bosons and vector bosons can be incorporated separately up to the considered order in perturbation theory). In the special case of CP-conservation the mixing between the physical Higgs bosons is reduced to a 2 × 2 mixing between the states h and H.
Analysing the pole structure of propagator matrices, we have shown that for the case of 3 × 3 mixing each entry of the propagator matrix has three complex poles, while for the 2×2 mixing case each entry of the propagator matrix has two complex poles. Consequently, JHEP09(2017)079 a single-pole approximation is not sufficient to approximate the momentum dependence of the full propagators. In particular for close-by states with sizeable widths, the different resonance regions may overlap.
We have derived in a process-independent way how the full propagators can be expanded around all of their complex poles. This approximation results in the sum of Breit-Wigner propagators of the corresponding resonances, weighted by wave function normalisation factors which encompass the relation between the lowest-order mass eigenstates and the loop-corrected mass eigenstates. As a key result, we have demonstrated that wave function normalisation factors that have been derived to ensure the correct on-shell properties of external particles in physical processes at their (in general complex) poles are a useful tool also for describing off-shell propagators, i.e. including momentum dependence. It has been shown that the momentum dependence of the full propagators can be accurately approximated by simple Breit-Wigner propagators. The complex wave function normalisation factors properly incorporate the imaginary parts of the self-energies arising from absorptive parts of loop integrals, which in contrast are neglected in effective coupling approaches where the contributing self-energies are evaluated at vanishing external momentum.
In our numerical analysis we have found very good agreement between the approximation presented in this paper and the full propagators. Besides the numerical comparison of the full and approximated propagators around the complex poles or depending on real momenta, we have provided a detailed discussion of the uncertainties involved in the approximation. In particular, the omitted higher-order terms from the expansion of self-energies around the real part of complex momenta, the expansion of the momentum-dependent effective self-energies around the complex poles and numerical effects in the on-shell condition have been quantified. The uncertainty estimates of these sources are in general well reflected in our numerical results. We have also discussed the validity of the approximation far away from the poles.
The formalism of Breit-Wigner propagators and on-shell wave function normalisation factors has several appealing advantages in describing the propagators of unstable particles that mix with each other. It avoids the momentum-dependent evaluation of self-energies and thereby significantly simplifies and accelerates the calculation of higher-order contributions. Besides, it enables the separation of the individual resonant contributions of the mass eigenstates h a within a given process and the straightforward calculation of their interference term. Moreover, the Breit-Wigner propagator turns into a δ-distribution in the limit of a vanishing width, thus facilitating the separate calculation of the production and decay of an intermediate unstable particle by means of the narrow-width approximation and its generalisation to the case of overlapping and interfering resonances. In addition, the Breit-Wigner formulation allows the implementation of a more precise total width by incorporating important higher-order effects from the partial widths that are not included in the imaginary part of the complex pole with self-energies evaluated at the same perturbative order as for the partial widths. This feature is very useful for phenomenological predictions of processes involving the exchange of unstable particles, benefitting from the use of quantities computed at the highest available order in perturbation theory.
The explicit calculations presented in this paper have been performed in the MSSM,

JHEP09(2017)079
but the introduced method itself can be easily extended to the cases of particles with nonzero spin and mixing among more than three particles, such as, for example, in different models with a non-minimal scalar sector or new vector resonances.