New Series Representations for the Two-Loop Massive Sunset Diagram

We derive new convergent series representations for the two-loop sunset diagram with three different propagator masses m1, m2 and m3 and external momentum p by techniques of analytic continuation on a well-known triple series that corresponds to the Lauricella Fc function. The convergence regions of the new series contain regions of interest to physical problems. These include some ranges of masses and squared external momentum values which make them useful from Chiral Perturbation Theory to some regions of the parameter space of the Minimal Supersymmetric Standard Model. The analytic continuation results presented on the Lauricella series could be used in other settings as well.

function. The convergence regions of the new series contain regions of interest to physical problems. These include some ranges of masses and squared external momentum values which make them useful from Chiral Perturbation Theory to some regions of the parameter space of the Minimal Supersymmetric Standard Model. The analytic continuation results presented on the Lauricella series could be used in other settings as well.

Introduction
The sunset integral is among the simplest of two-loop integrals that appear in the perturbation expansion of quantities in various quantum field theories, including the Standard Model (SM). In the convention of the classic work [1], the general massive sunset diagram of Fig. 1 is defined as T (p 2 , m 2 1 , m 2 2 , m 2 3 , α, β, γ) ≡ where the external momentum p 2 = s can take any value, n = 4 − 2 and m 1 , m 2 , m 3 are the masses of the three distinct particles in the propagators. The most general sunset integral, therefore, can have up to four independent mass scales. The sunset integral with general powers of the propagators (α, β, γ) can be reduced using integration by parts into a linear combination of a maximum of four master integrals (MI) [2]. These master integrals are the sunset integral with the following propagator powers: (α, β, γ) = (1, 1, 1), (2, 1, 1), (1, 2, 1), (1,1,2). It may also be p p q + r − p, m 3 q, m 1 r, m 2 Figure 1: The sunset diagram noted that all tensor and vector integrals, defined with momentum factors and appropriate Lorentz indices inside the integrals, as well as arbitrary derivatives with respect to the momentum and / or the masses can also be reduced to the same set of MI. Thus the study of the MI becomes of great importance in the context of the sunset integrals. Important special cases for renormalization theory include the fact that at p 2 = 0 the residues in dimensional regularization can be evaluated in closed form, see [3,4]. For configurations known as threshold and pseudo-threshold configurations the coefficient of the 0 is also known in terms of elementary functions, see [5]. Despite these early important advances, much work continues to be done on the evaluation of the sunset integrals. In particular, concerning the analytical evaluation of the sunset and more generally of the two-loop self energies, an implementation of analytical results is still far from complete, and their evaluation remains often numerical, see for example the recent phenomenological work of [6] in the context of the Minimal Supersymmetric Standard Model (MSSM).
As far as the analytic evaluation of the sunset diagrams is concerned we will not try to give a list of more recent publications dedicated to the sunset here but instead refer the reader to the extensive bibliography in [7] (see also the list given in [6]). In [7] the most general 4 mass scale sunset integral has been evaluated in terms of newly defined generalized elliptic functions. In the classic work [1], the most general 4 mass scale sunset integral has also been calculated analytically in terms of Lauricella multiple series. These can be obtained from introducing series representations for the Bessel functions in terms of which the propagators can be expressed, or by using among others Mellin-Barnes (MB) representations. It may be noted here that the analytic expressions given in [1] can be used only for restricted ranges of the values of the masses and squared external momentum values (p 2 ) derived from the convergence regions of the Lauricella series. In what follows, we will refer to these as known series representations of the general four mass scale sunset.
It may also be noted here that the sunset integrals also appear in effective field theories of the SM. One of the important effective field theories of the SM at low-energies is Chiral Perturbation Theory (ChPT), the effective theory of the pseudo-scalar octet of the light mesons associated with the spontaneous symmetry breaking of the approximate axial-vector symmetry, part of the symmetry associated with independent left-and right-chiral rotations of the QCD Lagrangian when the three lightest quark masses are sent to 0. The pseudo-scalar octet in the absence of iso-spin breaking and electromagnetic corrections has three masses, namely the pion, kaon and eta masses. At two-loop order, the masses and the decay constants of these particles involve sunset diagrams that themselves require the masses to be inserted into the two-loop integrals, see [8].
Driven by the needs of ChPT, in recent works [9,10] (see also [11,12]), analytic representations of the masses and decay constants of the light pseudoscalar mesons have been derived at the twoloop level. The relevant sunset diagrams that appear here have special values of the masses m 1 , m 2 and m 3 and squared external momentum p 2 (all in terms of m π , m K and m η ) that do not belong to the convergence regions of the analytic expressions given in [1]. One of the aims of the works [9,10] was to give the ChPT results in terms of series representations which are easy to handle and computationally efficient at the same time. Therefore, to get the required results in [9,10] we have followed the spirit of [1]. However, instead of considering the most general two-loop sunset integral as in [1], we have focused on the particular sunset configurations needed in ChPT: we have developed the Mellin-Barnes approach to evaluate the ChPT sunsets and analytical results were presented in terms of double series of Kampé de Fériet type (details of this derivation will be given in [13]).
In the present work we return to the most general massive sunset and we show that, as discussed in the conclusion and outlook section of [1], an analytic continuation procedure applied to the expressions presented in terms of Lauricella series in the latter paper allows one to provide new series representations of the sunset diagram, giving access to some ranges of values of the masses and squared external momentum which are not reachable using the results of [1]. Let us mention that these new ranges of values include, among others, the masses and squared external momentum needed in the ChPT context, which gives a cross-check of the results given in [9], [10] and [13]. The analytic continuation expressions that will be presented here are however more general than the ChPT results since p 2 is a free parameter whereas in ChPT it is fixed to the value of one of the squares of the meson masses, as we have discussed above. Therefore, we will see that our formulas may be applied in completely different contexts, and an illustration will be presented in the MSSM case [6]. Of interest also in physical applications is the expansion of the formulas as Laurent series in . We will discuss this point too, at any order in . However, to avoid very lenghty expressions, this will be detailed on a single particular example only. This paper is organized as follows: in Section 2, we recall the well-known representations of the 4-mass scale sunset diagram in terms of Lauricella functions. In Section 3, we develop the analytic continuation procedure on the latter, to obtain new series representations of the sunset diagram that are presented in Section 4 and in the Appendix. The -expansion procedure is briefly shown in Section 5. Then, in Section 6, some physical situations in ChPT and in the MSSM are considered and we show which series can be applied to them. Some conclusions follow.
2 Some known series representations of the sunset diagram Let us begin by recalling some known series representations of the sunset integral, first presented in [1] where small and large momentum expansions have been derived for the sunset MI T (p 2 , m 2 1 , m 2 2 , m 2 3 , 1, 1, 1). The small momentum expansion of the sunset diagram T (p 2 , m 2 1 , m 2 2 , m 2 3 , 1, 1, 1) is given as [1] where In Eq.
(2) appears the well-known Lauricella triple series where the superscript (3) indicates the number of its independent variables and (α) m ≡ Γ(α + m)/Γ(α) is the Pochhammer symbol. In the most general case, the Lauricella series are generalizations of hypergeometric series to multiple variables and can be of four different types: (see for instance [14,15,16] for details) but it is the case of three variables that is commonly referred to as a Lauricella series. The convergence region of the triple series in Eq.(4) is given by the inequality which is shown graphically in Fig. 2. In this figure, the convergence region is shown in R 3 0 , where namely the first (closed) octant (+ + +) of the (z 1 , z 2 , z 3 )-space. Indeed, since the convergence of an N -dimensional power series in z 1 , ..., z N is in general absolute, its convergence region D, which is a subset of C N , may be adequately described by the domain D in absolute space R N 0 , obtained by the natural projection (z 1 , ..., z N ) → (|z 1 |, ..., |z N |) [16]. The full convergence region (in the other octants) can be simply obtained by symmetry.
In the rest of this paper, however, giving the convergence regions only in the first octant will not be sufficient, because the analytic continuation formulas will be expressed as power series not in the z i but in some intricate combinations of the latter. It will therefore not be obvious to deduce regions of convergence in the full space from the first octant only. However, in order to simplify the figures as much as possible, we will show in this paper only the first octant and the fifth (the latter being defined as (+ + −)), since these are the most relevant from the physics point of view (indeed, from Eq.(3) we see that z 1 and z 2 are positive quantities while z 3 may be either positive or negative).
From Eqs. (3) and (5), it is clear that each of the series of Eq.(2) converges in the region shown in Fig. 3, which is the same as Fig. 2 with the addition of the fifth octant (z 1 ≥ 0, z 2 ≥ 0, z 3 ≤ 0). To derive the sunset expression for large-p 2 values, the authors of [1] have used the standard analytic continuation formula for the Lauricella F which is valid, due to Eq.(9), for |p 2 | > (m 1 + m 2 + m 3 ) 2 . This region is shown in Fig. 4. Note that Eq.(10) slightly differs from Eq.(24) of [1]. Indeed, in the latter, terms like (−z 3 ) −a coming from Eq.(8) were rewritten as (−1/z 3 ) a , which is not true when a is not an integer and z 3 > 0. Therefore we have preferred to leave these terms as they naturally appear when one applies Eq.(8) on Eq.(2). We have also corrected several misprints in the last two terms of Eq.(24) of [1]. Eqs. (2) and (10) are the main results of [1] for the sunset diagram (with unit powers of the propagators).

Invariance of the sunset diagram under permutation of masses
Due to the invariance of the Lauricella series under the permutation of its variables (and of the corresponding parameters) one can obtain from Eq.(8) the following alternative analytic continuation formula Therefore, from the physics point of view, one cannot learn anything new from S 3 and S 4 compared to S 1 , even if in the (z 1 , z 2 , z 3 )-space the convergence regions of the former are, obviously, not the same as the one of the latter. For illustration, these regions are depicted in Fig  It should be noted that, contrary to the cases of S 3 and S 4 considered in the present section, the invariance of the Lauricella series under the permutation of its variables will give birth to non trivial series representations of the sunset in Section 4.
It is clear from Fig. 5 that a large part of the (z 1 , z 2 , z 3 )-space cannot be reached by any of the series representations that we have discussed so far. It is argued in [1] (although not quantitatively shown) that by analyticity the convergence region of the total sum in the r.h.s. of Eq.(2) may be extended beyond the convergence region of each of its individual series, at the price of a careful numerical treatment at those points which belong to the former and not to the latter. This extended region is given by the values of the masses and external momentum that simultaneously C series in order to derive from them some series representations of the sunset diagram valid in the domain of the (z 1 , z 2 , z 3 )-space where none of the series representations presented until here can be used. We also wish to access the light green region of Fig. 6 without the careful numerical treatment mentioned in [1], and with a robust analytic method.

New analytic continuations for the Lauricella F
(3)

C series
To derive analytic continuation formulas for a multiple series, one possibility is to write some of the infinite sums of the multiple series in terms of a relatively common function for which analytic continuations are known or can be obtained, replace this function by one of these analytic continuations, and find the convergence region of the newly obtained multiple series, hoping that it will extend beyond the convergence region of the initial multiple series. This is the strategy that we adopt here for the analytic continuation of the Lauricella series and, more precisely, we will rewrite the latter as an infinite sum of Appell F 4 double hypergeometric functions. We indeed have derived new analytic continuation formulas for the Appell F 4 function in [17]. Note however that instead of making the Appell F 4 function appear in our calculations, one could have chosen to write single sums as p F q hypergeometric functions and performed a similar analysis (see [18] for detailed calculations with this alternative approach). In general, the analytic continuation of a multiple series is a sum of several multiple series and the convergence region of the resulting expression is given by the intersection of the convergence regions associated with each of the multiple series obtained at the end. It may also be necessary, for the validity of the new expression, to take into account the constraints coming from the range of validity of the analytic continuation formulas or other relations that have been used in the intermediate steps of the derivation (if any).
As a first simple illustration of this general procedure, we will now derive Eq.(8) from Eq.(4). To this end, one can either write one of the three sums of the Lauricella F C series as a Gauss 2 F 1 hypergeometric series, and then use a well-known analytic continuation formula of 2 F 1 or, equivalently, one can write two sums of the Lauricella series as an Appell F 4 series and use a known analytic continuation formula of the latter. Let us briefly demonstrate the second way.
Let us write where We now replace the Appell F 4 series in Eq.(16) using the analytic continuation formula [14] F 4 (α, β;γ, δ; x, y) = Note that F 4 = F one may then recognize Lauricella F C series in each of the two terms which appear at the end. This yields Eq. (8).
As a remark, let us mention that one can prove Eq.(18) in a similar way to what we just did, using the well-known analytic continuation of the Gauss hypergeometric function and, as said above, it is also possible to directly derive Eq.(8) by inserting Eq. (20) in the Lauricella series written as a double sum of Gauss hypergeometric series. It is remarkable that the r.h.s. of the analytic continuation formula Eq.(8) (Eq. (18)) is expressed only in terms of the same type of Lauricella function (Appell function) as its l.h.s 1 . This gives the convergence region of the r.h.s in a straightforward way, since it is sufficient to replace . Doing this, we obtain Eq.(9). Concerning Feynman integrals and other physical quantities of interest that are expressed as a sum of several different multiple series, one has to apply the analytic continuation procedure to each of the involved multiple series that need it, and the complexity of the intermediate steps can differ significantly from one multiple series to the other. However, in the case of the sunset, we have seen that the results given in [1] are written in terms of only one type of multiple series, the Lauricella F C Lauricella series are obtained, it will be possible to apply them to each of the different terms of Eqs.(2), (10), (12) and to the symmetrical equations obtained by the exchange of the variables (and corresponding coefficients) of the latter. Doing so, we will derive 1 Notice however that this nice property does not hold for non-trivial analytic continuations and, in fact, not even when one applies to above simple analysis on the other Lauricella series F B and F several new series representations of the sunset, which will greatly enlarge the reachable region of the (z 1 , z 2 , z 3 )-space shown in Fig. 5. In this obtained set of new formulas, some will be valid, among others, for ChPT and MSSM configurations.
A particularly relevant analytic continuation formula can be obtained for the Lauricella F (3) C series if one notices that the parameters of all the Lauricella series involved in the series representations of the sunset presented in Section 2 (in Eqs. (2), (10), (12) as well as in their symmetrical expressions) admit a common pattern. Indeed, owing to the general symmetry property one can see that, for appropriate choices of a and b, each of the four Lauricella series appearing in Eq.(2) can be rewritten in the form each of these appearing in Eq.(10) can be rewritten as etc.
In fact, it would also be possible to get the following patterns for all Lauricella series in Eq.(2), and similar relations for those of Eq.(10), etc. but due to the invariance properties of the sunset, one will not learn anything new from these other patterns compared to the one of Eq. (22), therefore we focus now on the latter. To obtain the analytic continuation of F C with this particular pattern of its parameters, one starts once more from Eq.(16) with e = a − b + 1 but then, instead of Eq.(18), one uses the following analytic continuation formula derived in [17] for the Appell F 4 function F 4 (α, β; γ, α − β + 1; x, y) = M (α, β, γ; x, y) where with N (α, β, γ; δ, λ; x, y) = ∞ n,p=0 x n n! y p p!
Eq.(26) is valid in the region where and Let us now proceed to the derivation of the analytic continuations of the Lauricella function in Eq.(22) from the Appell F 4 expression above.
We have A careful look at Eqs.(30), (31) and (32) shows that Eq.(26) is valid only for 0 < y < 1 [17]. Therefore, we can already guess from Eq.(33) that the analytic continuation of the Lauricella series deduced from Eq.(26) will be valid only in the 0 < z 3 < 1 region. The z 3 > 1 region case as well as negative values of z 3 will be discussed further.
3.1 0 < z 3 < 1 case After substituting Eq.(26) in Eq.(33), we get the following analytic continuation where we have defined theM function as the sum of the following three terms and We now have to determine the convergence regions of each of these three new series. To determine the convergence region of T 1 , we note that this series has the same convergence behaviour as the series This is due to the fact that the convergence properties of multiple gaussian hypergeometric series are independent of their parameters (excluding exceptional values of the latter) [16]. By expressing the sum over p as a 2 F 1 hypergeometric function after using the duplication formula on Γ(B + m + n + 2p), the above may be written as ∞ m,n,p=0 Now, using the following result from [19] 2 F 1 which is valid for |X| < 1, we can write if we set α = B 2 + m 2 + n 2 and F = B. As already mentioned, the latter choice does not affect the convergence of the underlying series in Eq.(40), which is dependent upon the structure of the summation parameters m, n and p rather than on the values of A, B, etc. The series can then be said to go as The series in Eq.(43) is an Appell F 4 double series whose convergence region is |r| + |s| < 1.
It is therefore established that the series in T 1 converges in the intersectionR 1 of this region with the region (1−z3) 2 (1+z3) 2 < 1 that comes from Eq.(42):R 1 ≡ ( |r| + |s| < 1) ∩ (1−z3) 2 (1+z3) 2 < 1. We now have to deal with the convergence properties of the series in T 2 . Its convergence region is determined following the same steps as for the case of T 1 above. Finally, one gets The region of convergence of the series in T 3 follows from a slightly different derivation, where the use of valid for |X| < 1, is needed [19]. At the end, one obtains the region of convergence of the series in T 3 asR 3 ≡ (|r 1 | + |s 1 | < 1) ∩ (|r 2 | + |s 2 | < 1) ∩ | z1 z2 | < 1 with and The intersection ofR 1 ,R 2 andR 3 is shown in Fig. 7 (red region). We see, as already guessed just before the present section, that this region is indeed restricted to positive values of z 3 . However, the final region of validity of the analytic continuation formula in Eq.(34) must also take into account the constraints coming from the intermediate steps of its derivation, namely those in Eq.(30). The R 0 constraint does in fact not provide any additional restriction, therefore the main constraint comes from R 1 and we conclude that the region of validity of Eq.(34) is the restriction to 0 < z 3 < 1 of the region plotted in red in Fig. 7. It is explicitly given as and shown in Fig. 8.

z 3 > 1 case
It is easy to derive the z 3 > 1 case from the results shown in the preceding section. To this end, we follow exactly the same steps of derivation as those presented in [17] in order to obtain a new analytic continuation formula for the Appell F 4 (α, β; γ, α − β + 1; x, y) function valid for y > 1: in the latter case, we first applied Eq.(18) on the l.h.s. of Eq.(26) and then used a slightly rewritten form of the M (α, β, γ; x, y) function, defined in [17] as M M (α, β, γ; x, y), to analytically continue both of the two F 4 functions that result.
In the Lauricella case, which generalizes the Appell ones, we obtained in the r.h.s of Eq.(34) an analytic continuationM (a, b, c, d; z 1 , z 2 , z 3 ), valid in the 0 < z 3 < 1 region, which in fact generalises the M (α, β, γ; x, y) analytic continuation, valid for 0 < y < 1, to the triple series case. Therefore, to obtain a series representation of the Lauricella F (3) C function valid in the |z 3 | > 1 region (i.e. the part of the red region of Fig. 7 that does not belong to the red region of Fig. 8), it is sufficient to first apply Eq.
This gives the following analytic continuation formula and it is straightforward to show that it is valid in the regioñ shown in Fig. 9. It is therefore just the restriction to z 3 > 1 of the region plotted in red in Fig.  7.

Other analytic continuation formulas
An obvious new analytic continuation formula for the Lauricella series with the specific pattern under study in this paper can be obtained by performing exactly the same analysis as in the preceding sections with the starting point instead of Eq.(33). This is equivalent to use the invariance of the Lauricella series under the exchange of z 1 and z 2 as well as their corresponding coefficients c and d, which preserves the pattern of the other coefficients.
One then obtains the two following analytic continuation formulas which is valid in the regionR(z 1 ↔ z 2 ) shown in blue in Fig. 10, and which is valid in the regionR (z 1 ↔ z 2 ) shown in purple in Fig. 10.

New series representations for the sunset diagram
We will now use the results presented in Section 3 to derive several new series representations of the sunset diagram. Let us begin by replacing each of the Lauricella series in Eq.(2) by their analytic continuation obtained from Eq.(34). We get a new series representation of the sunset diagram T (p 2 , m 2 1 , m 2 2 , m 2 3 , 1, 1, 1), valid in the regionR, defined in Eq.(49) and shown in Fig. 8, which can be written as S 7 cannot be obtained from any other series representations of the sunset shown in this paper by using the permutation of masses invariance property of the sunset diagram. From the phenomenological point of view, S 7 is quite efficient because one can use it at the same time for the kaon and eta cases of Chiral Perturbation Theory, as well as for all the "heavy" sunsets of MSSM considered in the next section. Moreover, S 7 is also able to reach negative values of z 3 (p 2 < 0) which is not the case for S 5 and S 6 .
Three other series representations of the sunset (S 8 , S 9 and S 10 ), derived using the invariance of the Lauricella series under the permutation of some of their variables (and which cannot be obtained from the others series S i (i = 1, ..., 7) by permutations of the masses) are given in the Appendix with their respective regions of validity, see Eqs.(81)-(83). S 7 (S 8 ) is obtained from analytically continuing S 1 after using the permutation of the variables z 2 (z 1 ) and z 3 , while S 9 (S 10 ) is derived from S 2 after permutating the variables z1 z3 ( z2 z3 ) and 1 z3 (as well as permuting the corresponding parameters). These 4 permutations of variables are the only non trivial ones in the set of the 10 possible permutations. The 6 others permutations give series representations that can be obtained from suitable permutations of the masses in the S i (i = 5, ..., 10).
Combining the regions of validity of our new series representations for the sunset with 4 mass scales, a large new area of the (z 1 , z 2 , z 3 )-space is now reachable (see Fig. 11 and compare with Fig. 5 which combines the regions of convergence of the 2 series representations given in [1] as well as those of the series obtained from them by permuting the variables).
The ellipsis in the arguments of the T i imply that they are equivalent to the argument of theM on the LHS of the equation.
T 1 (1 − ν, 1, 1 + ν, 1 + ν; z 1 , z 3 , z 2 ), after combining with the factors preceding theM of Eq.(60) and simplifying it, can be expressed as The -expansion of the above involves expanding the dependent part by means of the generalized product rule The derivatives of the Pochhammer, the inverse Pochhammer, and the quotient of two Pochhammer functions, are given in [20,21]: where and r j,m,n (α, β) = To be able to use the above formulae for the derivatives, it is necessary to re-express the Pochhammer functions of Eq.(62) in the following manner Applying Eq.(63) to the above, we get the expansion of Eq.(62) about = 0 A similar procedure yields the expansion of as and of as where The full expansion of Eq.(60) is therefore given by T 1, + T 2, + T 3, . The above exercise can be carried out for the other terms of S 7 and the other sunsets representations, as required by physical considerations.
In this section, S 7 was used as a representative sunset as it covers the "heavy" MSSM and the kaon and eta ChPT cases considered in the next section.

Some physical applications
In this section, we show two situations where the analytic continuation formulas of Section 4 and of the Appendix can be used. In each case, checks have been performed by a numerical comparison between our new series representations of the sunset diagram and the numerical integration of the corresponding Feynman parameter integral representations (choosing a finite value for ).

Chiral perturbation theory
The values m 1 = m 2 = m K , m 3 = m η (or any permutation of them) and p 2 = m 2 π define the pion sunset configuration in ChPT, whereas m 1 = m π , m 2 = m K , m 3 = m η (or any permutation of them) and p 2 = m 2 K define the kaon one, and m 1 = m π , m 2 = m 3 = m K (or any permutation of them) and p 2 = m 2 η the eta one. These sunsets diagrams, which depend on three different masses and which have been studied in [9,10,11,12], were the last missing pieces in the analytic expressions of the masses and decay constants of the light pseudoscalars at two-loop level in ChPT.
In Section 2, we recalled the well-known result for the general four mass scale sunset T 123 that converges for |z 1 | + |z 2 | + |z 3 | < 1, and which is given as a linear combination of the Lauricella F C series (see Eq.(16) of [1] or Eq.(2) in the present paper). Its analytic continuation to the region |z 3 | > |z 1 |+ |z 2 |+1 is given in Eq.(24) of [1] (see also Eq.(10) in the present paper, where some misprints have been corrected), again in terms of the F C function. The authors of [1] stated that provided that m 1 + m 2 < m 3 , one could extend the region of convergence of Eq.(2) to |p 2 | < (m 1 + m 2 + m 3 ) 2 , with careful numerical treatment in the annulus ( As already discussed in the introduction, with z 1 ≡ m 2 1 /m 2 3 , z 2 ≡ m 2 2 /m 2 3 and z 3 ≡ p 2 /m 2 3 , neither of the above cases corresponds to the physical kaon, eta or pion sunset configurations. This is true under any exchange of the masses in the propagators, which can be shown to leave the T 123 integral invariant. This is due to the fact that although the chiral perturbation theory sunsets do satisfy |p 2 | < (m 1 + m 2 + m 3 ) 2 , they always violate the m 1 + m 2 < m 3 condition, as well as m 1 + m 2 + |p 2 | < m 3 and m 1 + m 2 + m 3 < |p 2 |, for any of the permutations of the masses.
In fact, it is easy to check that the pion configuration, with the choice of masses above, does fall into the region of validity of S 8 which can then be used to evaluate it. This series is the only one that can be used for the pion case, even after permutations of the masses.
The kaon configuration does fall into the regions of validity of two of our new series which do have a non-empty intersection: S 5 and S 7 . One can therefore choose either one or the other to evaluate it.
In passing, although the particular set of masses and squared external momentum defined above for the kaon configuation m 1 = m π , m 2 = m K , m 3 = m η and p 2 = m 2 K does not belong to the region of validity of S 10 defined after Eq.(83), the numerical evaluation of the latter for this particular set is in fact in perfect agreement with the result coming from the corresponding Feynman parameter integral representation. This is due to the fact that this configuration of masses is exactly on the boundary z3 z2 = 1 of the region of validity of S 10 where the latter is in fact also valid. In Sections 3.1 and 3.2, we have performed an analysis for 0 < z 3 < 1 and for z 3 > 1, but we did not try to prove anything at z 3 = 1. The reasons are that on one hand results are in general harder to prove there and on the other hand, as we have seen in the kaon example, it is not necessary to use S 10 on its boundary because S 5 and S 7 can very well do the job. As a last remark, by appropriate permutations of the masses, it is also possible to use S 6 or S 9 , on their boundary, to evaluate the kaon configuration.
The eta configuration with m 1 = m π , m 2 = m 3 = m K and p 2 = m 2 η does fall into the regions of validity of S 6 and S 10 and on the boundary of S 7 . By an appropriate permutation of the masses, it can be shown that S 9 can also be used to evaluate the eta configuration.

Minimal Supersymmetric Standard Model
In the recent work [6], the full 2-loop QCD contributions to the lightest Higgs-boson mass in the MSSM with complex parameters have been computed. This implies sunset diagrams evaluations, which the authors of [6] computed, as well as other 2-loop self energy master integrals, using numerical packages such as SecDec [22] (we refer the reader to [6] for details). Among the MI considered in [6], those of the sunset type are the following 24 different configurations: for p 2 = m 2 H + or p 2 = m 2 W on the external legs, the possible corresponding internal masses are given by the triplets (mb 1 , m t , mg), (mb 2 , m t , mg), (mt 1 , m b , mg) and (mt 2 , m b , mg). For p 2 = m 2 h1 , p 2 = m 2 h2 , p 2 = m 2 h3 or p 2 = m 2 Z , the triplets are (mb 1 , m b , mg), (mb 2 , m b , mg), (mt 1 , m t , mg) and (mt 2 , m t , mg). Using the numerical values of this particular parameter point, it is easy to check that all the "heavy" sunsets (i.e those with p 2 = m 2 H + , p 2 = m 2 h2 and p 2 = m 2 h3 ) cannot be computed from series others than S 5 or S 7 , whereas the "light" ones do fall into the convergence region of S 1 .
It is interesting to note that S 7 is able to handle the kaon and eta cases in ChPT as well as all the "heavy" configurations of MSSM.

Conclusions and outlook
In this work, we have provided new series representations S i (i = 5, ..., 10) for the two-loop sunset diagram with four mass scales, following an analytic continuation procedure that we have applied to an old series representation of the sunset diagram given in [1] in terms of Lauricella F (3) C series. The S i (i = 5, ..., 10) are obtained, in our analytic continuation procedure, by the application of non-trivial analytic continuation formulas derived in [17] for the Appell F 4 function. Basing ourself on the example of S 7 , we then have shown how the -expansion of the S i can be performed, at any order in .
These new series representations give access to large regions of the mass scale parameters space that were unreachable with the representations presented in [1]. It may be recalled that the new regions encompass realistic physical situations, including chiral perturbation theory and related investigations on the lattice. Furthermore, in extensions of the SM such as the MSSM, one also encounters situations where the original series representations of [1] fail to converge but where we have shown that ours can be used. Although the subject is by now quite classical, the present investigation recalls that one must have a handle on rigorous analysis of convergence properties, analytic continuation, of even relatively simple configurations arising in field theory. In this context, although Lauricella functions are known for a long time, the mathematical study of their analytic continuations in terms of series representations is still, to our knowledge, far from complete (see [15] for some results). Since these objects appear in varied contexts of the physical science, from experimental physics in the description of the solid angle subtended at a disk source by a non-coaxial parallel-disk detector [23], to string theory, in the expression of exact string scattering amplitudes [24], in passing by general relativity, in the solutions of the equations of motion of test particle and photon in Kerr and Kerr-(anti) de Sitter spacetimes [25] and in many other situations, it would be interesting to try to fill this gap. Motivated by quantum field theory, the present investigation is a first step in this direction for the F C case, the generalisation to situations where d is not restricted to be equal to a − b + 1 being under way. Concerning the sunset diagram in particular, the exploration of the regions not covered by Fig.11 and Fig.5, as well as the L-loop generalisation, are left for future investigations.