Fine-Tuning in the 2HDM

The Two-Higgs Doublet Model (2HDM) is one of the most popular and natural extensions of the Higgs sector; but it has two potential fine-tuning problems, related to the electroweak (EW) breaking and the requirement of alignment with the SM Higgs boson. We have quantified the fine-tunings obtaining analytical expressions, both in terms of the initial 2HDM parameters and the physical ones (masses, mixing angles, etc.). We also provide simple approximate expressions. We have taken into account that the fine-tunings are not independent and removed the"double counting"by projecting the variations of the alignment onto the constant-${\rm v}^2$ hypersurface. The EW and the alignment fine-tunings become severe in different, even opposite, regions of the parameter space, namely in the regimes of large and small extra-Higgs masses, respectively; emerging an intermediate region, $500\ {\rm GeV} \leq \{m_H, m_A, m_{H^\pm}\} \leq 700\ {\rm GeV}$, where both are acceptably small. We also discuss a remarkable trend that is not obvious at first glance. Namely, for large $\tan\beta$ both the EW and the alignment fine-tunings become mitigated. In consequence, the 2HDM becomes quite natural for $\tan\beta\geq {\cal O} (10)$, even if $m_H, m_A, m_{H^\pm}$ are as large as 1500 GeV. We explain why this is not the case for the 2HDM stemming from supersymmetry. We have illustrated all these trends by numerically analyzing several representative scenarios.


Introduction
The Two-Higgs Doublet Model (2HDM) [1][2][3] is one of the most popular extensions of the Higgs sector. It naturally stems in prominent beyond-the-Standard-Model (BSM) scenarios, such as grand unification [4,5], supersymmetry [6] and axion models [7]. Besides, it has the potential to improve the performance of the conventional Higgs sector for important issues, e.g. the appealing possibility of Higgs-portal dark matter (see [8][9][10][11][12][13]). On top of all this, the 2HDM represents the simplest non-trivial extension of the SM Higgs sector, and thus a test bed of BSM physics, which could be probed in present and future experiments.
The scalar sector of the 2HDM consists of two scalar doublets, which leads to five physical scalar fields (h, H 0 , H ± , A), unlike the sole SM Higgs boson. On the other hand, the current experimental data require that the properties of one of these fields (typically the lightest CP-even Higgs mass eigenstate, h), such as quantum numbers, mass and couplings, are equal or very similar to those of the SM Higgs boson. This is called the alignment limit [14][15][16][17][18][19], which is essentially imposed by the observation.
Apart from its theoretical appeal, the 2HDM has two potential fine-tuning problems. The first one is related to the electroweak (EW) breaking, when the magnitude of the VEV, v 2 = (246 GeV) 2 , is much smaller than the squared-mass terms entering the theory. The second one is related to the alignment condition. The main goal of this paper is to explore these fine-tunings, looking for regions of the parameter space where both are mild or irrelevant.
In this analysis we will adopt an "agnostic" point of view, in the sense that we will not consider additional symmetries which constrain the model (beyond the usual Z 2 parity to avoid FCNC). There have been several analyses exploring this alternative direction, in particular to get an exact or approximate alignment [20][21][22][23][24][25]. As pointed out in ref. [22], an approximate alignment can arise from a softly broken global symmetry of the scalar potential, but this requires to extend the Yukawa sector, e.g. with vector-like top quark partners.
In section 2 we present the generic 2HDM, fixing the notation and providing analytical expressions for the Higgs VEVs and the alignment parameter. In section 3 we explain the method to evaluate the fine-tunings and discuss the expected regions where they become severe. In section 4 we provide analytical expression for all the potential fine-tunings, both in terms of the initial 2HDM parameters and the physical ones (masses, mixing angles, etc.) We also explain how to remove the "double counting" of tunings by projecting the variations of the alignment onto the constant-v 2 hypersurface. In addition, simple approximate expressions are also provided. The analysis shows that the electroweak and the alignment fine-tunings become severe in different regions of the parameter space. In section 5 we illustrate all these trends by numerically analyzing several representative scenarios. This allows to identify the regions where both fine-tunings are acceptably mild or even irrelevant. Finally, in section 6 we present our conclusions. The complete analytical expressions for all the fine-tunings are given in the Appendix.

The Two-Higgs Doublet Model
Let us briefly review the general formulation of the 2HDM following the notation and conventions of ref. [26]. Denoting by the two complex Y = 1/2, SU (2) L doublet scalar fields, the most general gauge-invariant renormalizable scalar potential is given by: where the m 2 11 , m 2 22 mass terms and the λ 1,2,3,4 couplings are real, while m 2 12 and λ 5,6,7 could be complex. The absence of flavour changing neutral currents (FCNCs) essentially requires each type of quark to couple to just one scalar doublet [27,28]. This is accomplished by means of a softly-broken Z 2 symmetry, which implies λ 6 = λ 7 = 0, while a non-zero m 2 12 value is still possible, see table 1 below.
Furthermore, to avoid dangerous CP-violating phenomena, we will assume throughout the paper that m 2 12 , λ 5 are real. In addition, stability bounds on the quartic coefficients imply (see e.g. [29]) We have also imposed constraints on the size of the various couplings from perturbativity and the requirement that unitarity is not violated in 2HDM scalar scattering processes, along the lines of Refs [30][31][32].
As long as the Higgs mass matrix possesses at least one negative eigenvalue, the scalar fields develop non-zero vacuum expectation values (VEVs): which must satisfy v 2 1 + v 2 2 = v 2 (246 GeV) 2 . As usual, we define the β angle such that The vacuum is CP-conserving provided |m 2 12 | ≥ λ 5 |v 1 ||v 2 | [14]. Besides, the fields can be redefined where λ 345 = λ 3 + λ 4 + λ 5 . A special instance takes place when either v 1 = 0 or v 2 = 0, i.e. t β = tan β = 0, ∞. This requires m 2 12 = 0 and corresponds to the so-called inert-doublet model [8]. This case is protected by the symmetry Φ 1 → Φ 1 , Φ 2 → −Φ 2 , and represents a somehow trivial situation where one of the two Higgs fields plays exactly the SM-Higgs role while the other one does not couple to the rest of the SM fields.
For the non-trivial cases (t β = 0, ∞), we can write the previous conditions in terms of v and t β , Eliminating v from Eqs. (8) and (9), we get a quartic equation for t β : which allows the replacement of either P 1 = 0 or P 2 = 0 by P β = 0 in the couple of minimization equations (8,9). The β angle introduced before relates the initial basis (in which the Z 2 -symmetry is apparent) to the so-called Higgs-basis which satisfies H 0 The associated quadratic (quartic) parameters of the Higgs scalar potential in this basis are denoted by Y i (Z j ), where i ∈ {1, 2, 3} (j ∈ {1, . . . , 7}) (expressions for Y i , Z j in terms of the initial parameters, m ij , λ j can be found in ref. [26]).
From the original eight scalar degrees of freedom in the two Higgs doublets, three Goldstone bosons are absorbed by the electroweak bosons, W ± and Z, and the remaining degrees of freedom correspond to five physical Higgs particles: two CP-even scalars (h and H, with the convention m h ≤ m H ), one CP-odd scalar (A) and one pair of charged Higgses (H ± ).
The A and H ± fields stem directly from the above H 2 doublet, namely A = √ 2 Im(H 0 2 ), with masses given by On the other hand, in the initial basis the CP-even neutral Higgs fields, ( which is diagonalized by (with the usual notation, c α = cos α, s α = sin α).
Obviously, for the trivial cases s β = 0 (c β = 0), i.e. t β = 0 (∞), the mass matrix (14) is diagonal from the beginning, since m 2 12 = 0. Then the SM-Higgs corresponds to Φ 1 (Φ 2 ) and α = −π/2 (0) . Otherwise, α is determined by the equation where we have used Eqs. (8,9) to eliminate m 2 11 , m 2 22 . The α angle is defined modulo π and we have chosen the convention −π/2 ≤ α ≤ π/2, thus c α ≥ 0. The mass eigenvalues (15) for the light and heavy neutral Higgses are then given by corresponding to the physical mass-eigenstates This is usually called the physical basis, which is related to the Higgs-basis (11) by a β − α rotation: The couplings of these Higgses to the gauge bosons, V = W ± , Z, are usually parametrized by the coefficients C h V = C H V = sin(β − α), which relate the former to the SM Higgs couplings where h SM is the SM Higgs. Similarly, the couplings of these Higgses to the fermions are parametrized by the C h F , C H F coefficients, which in turn depend on the initial couplings of the two Φ 1 , Φ 2 doublets to the fermions. The latter are strongly restricted by the requirement of the absence of dangerous FCNCs. As mentioned above, this is guaranteed by imposing a softly-broken Z 2 symmetry, which affects both the Higgs doublets and the fermion fields. The possible Z 2 charge assignments lead to the well-known four types of 2HDM, shown in Table 1 [33].
Then, the corresponding C V and C F coefficients are given in Table 2 for the four types of 2HDM.

The Alignment Limit
This limit occurs when the light CP-even Higgs, h, behaves as the SM Higgs, i.e. it presents the same couplings as the latter to all the SM fields. From the couplings of h to vector bosons and fermions shown in Table 2, it is clear that this is achieved for Alignment limit:  Since 0 ≤ β ≤ π/2, −π/2 ≤ α ≤ π/2, this is accomplished for α = β − π/2. In this limit the Higgs-basis and the physical basis are equivalent. Namely, from Eq. (22), h becomes aligned with ( √ 2Re H 0 1 − v) and H 1 plays the role of the ordinary SM doublet. Note also that in this limit the first term in the numerator of Eq.(18) vanishes, so m 2 h = O(λv 2 ), as in the SM. Up to now, the Higgs boson observed in the LHC looks remarkably similar to the SM one [34]. This is why the 2HDM is usually considered at the (either exact or approximate) alignment limit [16]. More precisely, ATLAS [35] and CMS [36] Run 2 Higgs data allow to determine (C h V , C h F ) within a ∼ 10% error (reduced to ∼ 7% when combining both analysis). These bounds are more restrictive when evaluated in a given model since C F 's are in general non-universal and correlated with C V . In the 2HDM, these correlations are given in Table 2. In particular, for type I, deviations from C h V = 1 are ∼ 1% and [37,38] |c β−α | 0.15 .
The range is even more reduced in type II. It is worth commenting that there are symmetries of the potential [39][40][41][42] that lead to an exact alignment, i.e. c β−α = 0. A rather trivial instance is the above-mentioned inert-doublet scenario, where t β = 0, ∞ and there is an exact Φ 1 → Φ 1 , Φ 2 → −Φ 2 symmetry, so that one of the Higgs doublets precisely corresponds to the SM-Higgs, coupled to all fermions, while the other one is completely decoupled.
Beside this special case, the existence of symmetries that lead to alignment can be explored by eliminating t β in {P α , P β } = 0 evaluated at the alignment limit, t α = −1/t β ; namely subject to the consistency condition Whenever Eq. (27) is fulfilled by the initial parameters of the potential, there exists a minimum with t β given by Eq. (26) where the alignment is exact. Notice, in particular, that Eq (27) is satisfied for some simple relations between the potential parameters, such as {λ 1 = λ 2 & m 2 11 = m 2 22 }, which can be achieved by imposing certain global symmetries in the potential [20][21][22][23]. However, for phenomenological reasons these symmetries cannot be exact. Still, an approximate alignment can arise if the symmetry is softly broken in the scalar sector. Nevertheless, as pointed out in ref. [22], this requires to extend the Yukawa sector, e.g. with vector-like top quark partners coupled to the Higgses, which goes beyond the scope of this paper.
In any case, even if the parameters are in a suitable combination to get exact or approximate alignment, if the underlying symmetry is not exact one must vary the parameters in a free way to evaluate the fine-tuning. This is exactly what we have done in the present analysis, which we expose below, so these possibilities are taken into account.
3 Fine-tuning in the 2HDM. Decoupling and non-decoupling regimes A theoretical model presents fine-tuning (or, equivalently, absence of naturalness) when some observable quantity depends critically on a fine adjustment (or "conspiracy") of the fundamental parameters. Such adjustment in the parameter space is conceptually problematic, as it is implausible unless it can be explained from the theory itself.
Regarding the 2HDM, we distinguish two potential fine-tunings. The first one is related to the EW breaking, when the magnitude of the vacuum expectation value, v 2 , is much smaller than the other squared-mass terms entering the theory. Namely, from Eq.(8) or Eq. (9): Since , then v 2 is likely to be fine-tuned. This is sometimes called the "little hierarchy problem" [43]. Given that in the 2HDM there are two expectations values, v 1 and v 2 , or equivalently v 2 , tan β; one can wonder if there are two independent fine-tunings, i.e. if tan β might also be a fine-tuned parameter. We will examine this issue throughout the paper. Let us note that the mass of the SM-like Higgs, m h , has electroweak size, however it is not a fine-tuned parameter since, as mentioned above, m 2 h = O(λv 2 ), as in the SM. Therefore, once the value of v 2 has been set, m 2 h , similarly to M 2 W , M 2 Z , is naturally of electroweak size. The second potential fine-tuning is related to the magnitude of c β−α in the alignment regime, i.e. |c β−α | 1, since, in principle, there is no reason why the initial parameters should yield such small value. To see this more closely, let us use the expression given in ref. [26] that links c β−α with m 2 h , m 2 H and the Z 1,6 parameters: where Demanding an approximate (exact) alignment limit is equivalent to demand a small (vanishing) value of |c β−α |. This can be achieved in two ways: , the denominator of Eq.(29) becomes in this case much larger than the numerator, leading to a small |c β−α |. This instance is called in the literature alignment by decoupling [14,17,19]. It is worth-noticing that this limit requires large m 2 12 v 2 (see Eq. (17)). Then, from Eqs. (12,13,17), all the extra Higgs states get similar masses: Actually, in this regime one can integrate out the heavy Higgs states and the resultant effective theory is essentially the SM (more precisely, an SMEFT), with h playing the role of the SMlike Higgs boson. This may seem quite natural, but, as discussed above, it implies a potential fine-tuning related to the EW breaking.
2. Z 6 1. Then, |c β−α | 1 without requiring large Higgs masses. This regime is commonly called alignment without decoupling [16,18,44]. However, as it is clear from Eq. (30), this typically requires a precise cancellation inside Z 6 = i O(λ i ). A possible exception, however, occurs when s 2β 0, i.e. for t β 1. Note here that the t β 1 case is excluded because it leads to a non-perturbative Yukawa coupling for the top quark (which, by definition, is coupled to the Φ 2 doublet).
Hence, the alignment generically requires fine-tuning, whether it is achieved by decoupling or not. Thus, evaluating both potential fine-tunings in detail may lead to a better understanding of the parameter space, and to find regions in which none of them is too large. This is the main goal of this paper.
In order to quantify the fine-tuning associated with a generic observable, Ω, with respect to an initial parameter, θ i , we will use the somewhat standard criterion proposed by by Ellis et al. [45] and Barbieri and Giudice [46]. Namely, we define the fine-tuning parameters as Then ∆Ω ∼ 10 (100) denotes a fine-tuning of about 10% (1%), etc [46]. As discussed above, the set of (potentially fine-tuned) observables that we will consider is while the set of initial parameters is Our aim is to find analytical expressions for ∆ θ i v 2 , ∆ θ i t β and ∆ θ i c β−α in terms of the initial parameters or (most conveniently) the physical parameters, and discuss their magnitude in different regimes.
4 Analytical expressions for the Fine-tuning

Fine-tuning in terms of the initial parameters
The most straightforward way to get analytical expressions for the fine-tuning on any observable, Ω, is to start with its explicit dependence on the initial parameters, Ω(θ i ), and evaluate the derivatives dΩ/dθ i involved in Eq. (32). However, in our case these kinds of expressions are cumbersome; e.g. t β is given by the quartic equation (10). Fortunately, we can still extract analytical expressions for the fine-tuning using the constraints P 1,2,β,α = 0 (see Eqs. (8,9,10,16)), together with the Implicit Function Theorem. Next we do this for the three potentially fine-tuned observables, v 2 , t β , c β−α .

Fine-tuning in t β
For convenience, we start with the fine-tuning in t β . Taking derivatives in Eq. (10), Hence, replacing here the explicit expression (10) for P β we get analytical formulas for all ∆ θ i t β .
The above expression is enough to see that, in general, ∆ θ i t β < ∼ O(1), and thus t β is not a fine-tuned parameter. To check this, note that the size of the denominator is typically > ∼ O(m 2 ij ), barring accidental cancellations; and this is also the expected size of the numerator for any θ i . E.g. let θ i be one of the initial mass-parameters (the ones more directly involved in the electroweak fine-tuning), say θ = m 2 12 . Then Eq.(36) reads which is typically < ∼ O(1). This is obviously the case for λ i , t β = O(1), and also for t β 1. (In the last instance |∆ m 2 12 t β | ∼ 1/4.) The fact that t β is not fine-tuned in most cases can be traced back to the origin of this observable. If m 2 11 > 0, the VEV of Φ 1 is triggered by that of Φ 2 through the linear term in the potential (2), Hence, v 1 can have a similar size as v 2 with no need of tunings or cancellations (whether or not v 2 has electroweak size); thus t β can be O(1), or actually any size, with no fine-tuning. Note that a very large (or small) t β implies a strong hierarchy between m 2 11 and m 2 12 . This may be considered as an odd fact (as it is e.g. the hierarchy of masses of the fermionic generations), but it is not a fine-tuning.

Fine-tuning in v 2
For the fine-tuning in v 2 , we proceed in a similar way. From Eqs. (8,9), we get which leads to: Thus, replacing the explicit expressions (8,9) for P 1,2 , we get two alternative but equivalent analytical forms for each ∆ θ i v 2 . Now it is easy to see that v 2 is generically fine-tuned with respect to the initial mass-parameters. E.g. for θ i = m 11 , we replace into Eq. (41), getting . Hence, ∆ θ v 2 captures the whole fine-tuning associated to the scale of the EW breaking (the little hierarchy problem).

Fine-tuning in c β−α
Let us finally consider the fine-tuning in the alignment, i.e. ∆ θ i c β−α . Using the trigonometric identity in Eq.(32) (for Ω = c β−α ) we can express the fine-tuning in c β−α as where ∆ θ i t α ≡ θ i t α dt α dθ i is the "fine-tuning in α", which can be obtained from Eq.(16) following similar steps as for the other fine-tunings: Hence, Replacing here the explicit expression (16) for P α , beside Eqs. (36,41), we get analytical formulas for all ∆ θ i t α . Finally, the alignment fine-tuning reads which, upon replacement of Eqs. (36,49), has an explicit analytical expression for all θ i . Note that, in the alignment limit, 1 + t α t β → 0, thus we expect c β−α to be generically fine-tuned. It is worth-noticing that the term proportional to ∆ θ i t α has "inside" the electroweak fine-tuning, ∆ θ i v 2 , through Eq. (49). Actually, it can be checked from the previous expressions that for m 2 ij v 2 (decoupling limit) and 1 + t α t β → 0 (alignment limit) the coefficient of ∆ θ i v 2 is 1: This remarkably simple result can be understood as follows. Let us recall from Eq.(29) that Hence, in absolute value, The first term is exactly ∆ θ i v 2 , the second one represents the fine-tuning in c β−α when there is no decoupling and the third term is Eqs. (52, 53) show in a transparent way why the electroweak fine-tuning, ∆v 2 , appears involved in the the alignment one. The definition of fine-tuning takes into account all the adjustments required to get the desired value of the observable under consideration. In our case, if the alignment is achieved through decoupling, there is a fine-tuning price, namely the one required to get v 2 much smaller than the mass parameters of the potential.
From the above discussion it is clear that the fine-tuning in v 2 and c β−α are not independent. Thus, it would be incorrect to multiply them to get the "total" fine-tuning. In particular, in the decoupling limit, once the initial parameters have been adjusted to get the right size of v 2 , there is no need of extra adjustments to get alignment. Hence, a practical way to remove the "double counting" in the evaluation of the c β−α fine-tuning is to discard variations in the θ i parameters that change the value of v 2 . In the {log θ i } space, this is equivalent to project the gradient ∆ θ i c β−α = ∂ log c β−α /∂ log θ i onto the constant-v 2 hypersurface [47]: where

Fine-tuning in terms of the physical parameters
Very often it is convenient and more meaningful to express the fine-tuning in terms of the physical parameters of the model, rather than in terms of the initial ones. To this end, we have to trade the eight initial parameters, θ i , of Eq.(34) for eight physical observables, χ i . Specifically, we have chosen Although λ 5 is not strictly an observable, it is related to the triple scalar self-couplings, e.g.
Now the fine-tunings ∆ θ i t β , ∆ θ i v 2 , ∆ θ i c β−α , given by Eqs. (36,41,50), can be written as functions of the physical parameters, χ i , by replacing the initial parameters in terms of the latter, i.e. θ i (χ j ). This can be done with the help of the relations obtained in Sect.2. Namely, from Eqs. (12,13) we derive m 2 12 (χ i ), λ 4 (χ i ). Then, from Eqs. (16,17,18) we derive λ 1,2,3 (χ i ); and finally, from Eqs. (8,9) we get m 2 11 (χ i ), m 2 22 (χ i ). Altogether, the explicit expressions for θ i (χ j ) read Incidentally, the change of variables from the initial parameters to the physical observables corresponds to a remarkably simple Jacobian, We present the final analytical expressions for all the fine-tunings in the Appendix. This is one of the main results of the present work.

Results in −expansion
Due to the smallness of c β−α in the alignment regime, we can rewrite all the fine-tuning expressions of the Appendix as power series in c β−α = , using This expansion yields simpler formulas, which capture the general behaviour of the fine-tuning and allow for a closer examination. Next we gather and discuss the expressions of the various fine-tunings at the lowest order in . We will focus the discussion on the two instances that, as argued in Sect.3, naturally lead to alignment, namely the decoupling limit and the t β 1 regime.

Fine-tuning in
In the decoupling regime, m 2 indicating that generically the EW fine-tuning diverges as the masses of the extra Higgses increase, as expected from the discussion in section 4.1. More precisely, in that regime In contrast ∆ λ i v 2 remain O(1). Let us consider now the t β 1 regime. It can be checked from Eqs. (60) that in this case all ∆ θ i v 2 fine-tunings tend to zero, except ∆ m 2 22 v 2 , ∆ λ 2 v 2 , which become O(1), i.e. irrelevant. This remarkable absence of EW fine-tuning comes from the following. In the t β → ∞ limit (this can be seen, e.g. from Eqs. (7) or (9) As already mentioned, a hierarchy of mass terms may be or may be not seen as an odd fact, but it does not imply a fine-tuning between parameters.
This situation contrasts with what happens in the supersymmetric case. In the minimal supersymetric Standard Model (MSSM), the t β 1 regime (which is desirable in order to reproduce the ordinary Higgs mass) implies a notorious fine-tuning. This occurs because m 2 22 at low energy contains several contributions proportional to the soft supersymmetric masses, e.g. a large and negative contribution proportional to the mass-squared of the gluinos, and a positive contribution equal to the µ−term (squared). Then, the smallness of |m 2 22 | can only be achieved by a cancellation between all these large contributions, and this does imply a fine-tuning. However, this does not need to be the case in a generic 2HDM.
Of course, the previous discussion has nothing to do with the famous hierarchy problem [48,49], which is related to the stability of the scalar mass-parameters under quadratic radiative corrections.
In summary, from Eq.(61) the electroweak fine-tuning in the 2HDM has the approximate size which shows that it increases with the masses of the extra Higgses (as expected), but it can be compensated by a large t β .

Fine-tuning in
It is clear from the previous expressions that t β is not a fine-tuned parameter in any instance, as anticipated in section 4.1. In particular, in the decoupling regime, m 2 . Something similar occurs in the large t β regime, where all the fine-tuning parameters are smaller than 1, except ∆ m 2 11 t β ∼ ∆ m 2 12 t β ∼ 1. This comes from the fact that if m 2 11 > 0, then t β ∼ m 2 11 /m 2 12 , as discussed around Eq. (39).
In summary, t β is not a fine-tuned parameter and typically

Fine-tuning in c β−α
In this case the lowest order in the expansion is always 1/ , reflecting the fact that we are evaluating the fine-tuning in c β−α = itself, which is logically more severe as shrinks: Let us note that the O( 0 ) term, that has been dropped above, precisely contains the ∆ θ i v 2 contribution to ∆ θ i c β−α , see the discussion around Eq.(54). Actually, Eqs.(69) are a good approximation of the independent ∆ θ i c β−α fine-tuning, i.e. once it is projected onto the constant-v 2 hypersurface, according to Eq.(54). Consistently with that, we note that at this order in the −expansion, , thus the alignment becomes natural in the decoupling limit.
showing that the alignment becomes non-fine-tuned in that regime. The origin of this remarkable fact was already mentioned in Sect.3, when discussing Eq. (29).
In summary, the alignment fine-tuning in the 2HDM is of the order which shows that it is mitigated in the decoupling limit and for large t β .

Numerical Analysis
We perform in this section a numerical exploration of the the fine-tuning(s) in the 2HDM parameterspace. We will focus on their dependence on the masses of the extra Higgses and on t β ; which, as discussed above, are the main determinants of the size of the fine-tuning. We will also consider two possible circumstances: that the alignment is sharp, i.e. quite smaller than the present experimental upper bounds, or close to them.
For these purposes we have designed five scenarios which capture all the relevant features of the problem, see Table 3. We examine them in order.  Table 3: Benchmark Scenarios used in the numerical analysis to explore the electroweak and alignment fine-tunings in the 2HDM. The fact that m A = m H ± implies λ 4 = λ 5 in all cases.

Scenario 1
This scenario is chosen to show the dependence of the fine-tuning on the extra-Higgs masses in the case of sharp alignment, = c β−α = 0.01 and low t β . As discussed above, in the alignment limit the larger the extra-Higgs masses, the closer each other become. This feature is trivially fulfilled here as we have set m H = m A = m H ± ≡ m 0 . On the other hand, the value of λ 5 has been chosen to enhance the range of m 0 that preserves perturbativity and stability, but it does not play an important role in the discussion. The electroweak fine-tuning, ∆ θ i v 2 vs. m 0 is shown in Fig. 1 for the various θ i parameters. As expected, the only important fine-tunings are those associated with mass parameters, m 2 ij , which, for large m 0 increase according to the trend of Eq.(61), becoming as large as 250 for m 0 ∼ 2 TeV. In contrast, for m 0 < ∼ 700 GeV the electroweak fine-tuning is quite mild. Concerning the alignment fine-tuning, we have projected the complete expressions for ∆ θ i c β−α given in the Appendix onto the constant-v 2 surface, according to Eq.(54), in order to show only the independent fine-tuning. As commented at the end of the previous section, this is essentially equivalent to use the lowest-order expressions for ∆ θ i c β−α , Eqs.(69). The corresponding results are shown in Fig. 2. As expected, the alignment fine-tuning increases for decreasing masses of the extra Higgses, becoming as large as ∼ 1000 at m 0 = 200 GeV. . Finally, the fine-tuning in t β is shown in Fig. 3. As expected, t β is not a fine-tuned quantity, i.e. the associated fine-tuning parameters, ∆ θ i t β are O(1) in all cases.

Scenario 2
This is similar to Scenario 1, but with a weak alignment; namely, = c β−α = 0.1, about as large as allowed by experimental constraints. An important consequence of the mild alignment is that now m 0 cannot be too large. The reason is that for very large masses of the extra-Higgses we enter the decoupling regime, which automatically leads to alignment. Under such circumstances, it is very difficult to keep a certain misalignment (typically it requires non-perturbative λ−couplings).
In consequence the available range for m 0 lies below ∼ 1 TeV. As we will see later, for higher t β , this upper bound becomes stronger.
Apart from this limitation, the trends of the fine-tuning are similar to those for Scenario 1, as shown in Figs. 4-6. The main difference has to do with the alignment fine-tuning, ∆c β−α , which is less severe than before (compare Figs. 2 and 5) because the required value of c β−α is not that small. Again, there is an intermediate region, 500 GeV < ∼ m 0 < ∼ 700 GeV, where both the electroweak and the alignment fine-tunings remain at acceptable values, < ∼ O(10). same as Fig. 1 for Scenario 2.

Scenario 3
This scenario is chosen to show the dependence of the fine-tuning on t β in the case of sharp alignment, = c β−α = 0.01 and low extra-Higgs masses, ∼ 250 GeV. The value of λ 5 has been picked to ensure perturbativity and stability in a significant range of t β , namely t β ∈ [1,15].
The electroweak fine-tuning, ∆ θ i v 2 vs. t β is shown in Fig. 7. As discussed in the previous section, the fine-tuning decreases with increasing t β . Since this scenario has a rather low value of m 0 , the electroweak fine-tuning is never large or even significant.
The fine-tuning in the alignment, ∆ θ i c β−α vs. t β , is shown in Fig. 8. As discussed above, this fine-tuning also gets less severe as t β increases, and for t β > ∼ 10 it becomes ≤ O (10). It is worth-noticing that several fine-tuning parameters, namely ∆ m 2 12 c β−α , ∆ λ 3 c β−α , ∆ λ 4 c β−α and ∆ λ 5 c β−α , abruptly fall to zero for t β → 1. This is a general fact, as can be verified from the general expressions (69) for these quantities, which have a (t 2 β − 1) factor. The origin of this feature is the following. It is easy to check from Eqs. (8,9,26) that t β 1 and a fine alignment require λ 1 λ 2 , m 2 11 m 2 22 . This corresponds to a situation where the consistency condition for perfect alignment, (27) is approximately fulfilled in an obvious way. In that case, it is clear that changes in m 2 12 , λ 3 , λ 4 , λ 5 do not spoil condition (27); in other words, the fine alignment is insensitive to these changes, thus the vanishing of the associated fine-tuning parameters at t β = 1. This rule does not apply to the other initial parameters, m 2 11 , m 2 22 , λ 1 , λ 2 . As it was commented after Eq. (27), in the λ 1 = λ 2 , m 2 11 = m 2 22 . limit (which corresponds to t β = 1 and perfect alignment) there arises a symmetry in the potential, which is however broken by the Yukawa couplings.
Finally, as usual, there is no fine-tuning in t β , as shown in Fig. 9.

Scenario 4
This scenario explores the dependence of the various fine-tunings on t β , in a scenario with large extra-Higgs masses, ∼ 1500 GeV, and very sharp alignment, = 10 −3 . The results are shown in Figs. 10-12. Such a strong alignment is actually the natural consequence of two features that push in that direction: the decoupling due to the large extra masses, and the sizeable t β . This is demonstrated by the fact that, even with such alignment, there is no relevant ∆c β−α fine-tuning (Fig. 11). For t β ≤ O(10), there is still a (strong) electroweak fine-tuning, ∆v 2 , see Fig. 10, due to the large extra Higgs masses. This fine-tuning is compensated by the competing effect of a large t β for t β ≥ O (10).
Finally, once more, t β is not fine-tuned, as shown in Fig. 12.   This scenario explores the dependence of the various fine-tunings on t β , for small extra-Higgs masses, ∼ 300 GeV, and mild alignment, = 0.1 (the opposite to scenario 4). The results are shown in Figs. 13-15. In this case, due to the low extra-Higgs masses, there is no appreciable electroweak fine-tuning, Fig. 13. Likewise, due to the (as large as possible) value of c β−α , there is no fine-tuning for the alignment, Fig. 14. Of course, t β is not fine-tuned either, Fig. 15.

Conclusions
The Two-Higgs Doublet Model is one of the most popular and natural extensions of the Higgs sector. It arises in many BSM scenarios, such as grand unification, supersymmetry and axion models. But, besides its theoretical appeal, the 2HDM has two potential fine-tuning problems. The first one is related to the electroweak breaking, when the magnitude of the VEV, v 2 = (246 GeV) 2 , is much smaller than the squared-mass terms entering the theory. This is the familiar electroweak fine-tuning or little-hierarchy problem. The second one is related to the alignment condition. Namely, current experimental data require that the properties of one of the Higgs fields are equal or very similar to those of the SM Higgs boson, i.e. the model must be close to the alignment limit.
In this paper we have explored how these fine-tunings arise in the 2HDM (with the usual Z 2 parity implemented to avoid danegerous FCNC), looking for regions of the parameter space where both are mild or irrelevant. In order to quantify the fine-tunings we have used the fairly standard Barbieri-Giudice criterion; e.g. the fine-tuning in v 2 respect to some initial parameter θ is quantified by the fine-tuning parameter ∆ θ v 2 = ∂ log v 2 /∂ log θ.
Concerning the electroweak fine-tuning, since there are two vacuum expectation values, v 2 1 + v 2 2 = v 2 , one may wonder if there are fine-tunings associated with both v 2 1 , v 2 2 , or equivalently with v 2 and t β = v 2 /v 1 . We have shown that, in all cases, t β is not a fine-tuned parameter, so the whole electroweak fine-tuning is captured by ∆v 2 . Regarding the alignment, the relevant fine-tuned observable is c β−α = cos(β − α), which vanishes for perfect alignment and is currently constrained by experimental bounds to be |c β−α | < ∼ 0.1. We have obtained analytical expressions for ∆ θ i v 2 , ∆ θ i t β , ∆ θ i c β−α both in terms of the initial parameters of the theory, {m 2 i,j , λ i } and on the physical observables, {m H , m A , m H ± , t β , etc.}. As already mentioned, there is no fine-tuning for t β , i.e. ∆ θ i t β = O(1) in all cases. On the other hand, the electroweak and alignment fine-tunings are not independent and they should not be multiplied without further ado to get the "total" fine-tuning.
In particular, in the decoupling limit, once the initial parameters have been adjusted to get the right v 2 , there is no need of extra adjustments to get alignment. To remove the "double counting" in the evaluation of the c β−α fine-tuning one has to discard variations in the θ i parameters that change the value of v 2 . This is equivalent to project the gradient ∆ θ i c β−α = ∂ log c β−α /∂ log θ i onto the constant-v 2 hypersurface. This is the procedure we have followed to present the ∆ θ i v 2 , ∆ θ i c β−α fine-tunings as independent parameters.
We have also obtained approximate expressions for the various fine-tunings performing an expansion in c β−α . The general trend for them is: We have investigated numerically the behaviour of the fine-tunings in the the 2HDM parameter space, by using representative benchmark scenarios that explore different relevant regions.
As it is clear from the previous general trends, for moderate values of t β , the electroweak and the alignment fine-tunings, ∆v 2 and ∆c β−α , become severe in different regions of the parameter space, namely in the regimes of large and small extra-Higgs masses, respectively. The first one corresponds to the decoupling limit. This means that there is an intermediate region, more precisely 500 GeV < ∼ {m H , m A , m H ± } < ∼ 700 GeV, where both tunings are acceptably small. It should be mentioned here that, if one does not take into account the electroweak fine-tuning it is not correct to say that the alignment becomes natural in the decoupling limit. In that case, one should not project ∆c β−α onto the constant-v 2 hypersurface. Then it contains a contribution of order ∆ θ i v 2 that shows that it is in fact fine-tuned. This occurs because the Barbieri-Giudice parameters incorporate all the adjustments required to get the desired size for the observable under consideration.
The analytical expressions and the numerical results show an interesting trend that is not obvious at first glance. Namely, for large t β both the electroweak and the alignment fine-tunings become mitigated. We have discussed in the paper the reason for this feature. In consequence, the 2HDM becomes quite natural for t β ≥ O(10), even if m H , m A , m H ± are as large as 1500 GeV.
This situation contrasts with the supersymmetric case, where the t β 1 regime implies a notorious fine-tuning. The reason for that was that the mass parameter m 2 22 at low energy contained several contributions proportional to the soft supersymmetric masses and the µ−term. Then, the smallness of |m 2 22 | could only be achieved by a fine cancellation of such contributions. This does not need to be the case in a generic 2HDM.
In summary, the 2HDM has potentially severe fine-tunings associated with the smallness of the electroweak scale and the sharpness of the alignment, but there are interesting regions in the parameter space where both remain at acceptable levels.