Scalar Democracy

We conjecture that there exists a scalar bound state for every pair of fundamental fermions at a UV (`composite') scale, $\Lambda\gg v_{\text{weak}}$. This implies a large number of universally coupled, sub-critical Higgs doublets. All but the Standard Model Higgs are `dormant,' with large positive squared masses and each receives a small vacuum expectation values via mixing with the Standard Model Higgs. Universal couplings, modulo renormalization group running effects, flips the flavor problem into the masses and mixings of the Higgs system. Doublets associated with heavy fermion masses, $b,c, \tau$ likely lie in the multi-TeV range, but may be observable at the current LHC, or a high-luminosity and/or an energy-upgraded LHC. In the lepton sector we are lead to a Higgs seesaw for neutrino masses, and corollary processes of observable flavor violation. The observation of the first sequential doublet coupled to $\bar{b}b$ with masses $\lesssim 3.5$ TeV would lend credence to the hypothesis.


I. INTRODUCTION
In the present paper we propose that every fermion pair binds to form a complex scalar boson, due to a universal attractive interaction at a very high scale, Λ. Amongst many new states, including lepto-quarks, colored isodoublets and singlets, etc., this hypothesis implies the existence of a large number of Higgs bosons.
We assume the lightest of the these is the Standard Model Higgs (SMH). The remaining doublets are sequentially heavier with positive M 2 's, i.e., "dormant." 1 They will have universal couplings to their constituent fermions at Λ, but renormalized couplings at the electroweak scale that are g 1 for quarks, and g 0.7 for leptons. Each Standard Model (SM) fermion acquires its mass through its coupling to the particular Higgs doublet it comprises, which in turn mixes with the SMH to acquire a perturbative "tadpole" mass. In particular our present model is "subcritical:" the negative M 2 of the SMH arises from mixing with dormant Higgses.
The associated Higgs bosons thus become lighter for the heavier fermions, and very heavy for the neutrinos. Though we have no theory of the masses and mixings of the large array of composite scalars, we can make use of phenomenology. The model then becomes predictive, essentially because the Yukawa couplings are determined. We call this system "Scalar Democracy" as it harkens back to the "Nuclear Democracy" of the late 1960's. Sequential Higgs bosons have certainly been considered previously, and we cannot review the vast literature.
Nonetheless, few theorists venture beyond a few Higgs bosons. Many Higgs bosons arise in the context of an extended gravity, such as a scheme described to us by Bjorken [1], which first inspired the idea of Scalar Democracy. A model that presages some of our present discussion is the "Private Higgs" [2,3]. A Coleman-Weinberg mechanism and Higgs portal interactions lead one of us to a second Higgs doublet coupled to bb with g = 1 at about 400 GeV, (see [4] and a long list of references therein). In this model, the SMH is a tt composite state [5][6][7], but the predictions of minimal composite tt are significantly modified by mixing. We note that recent work attempting to construct an "asymptotically safe" UV completion of the SM arrives at a similar structure in the Higgs sector to ours [8][9][10], but for different reasons. The phenomenology of many Higgs bosons is, to our knowledge, essentially unexplored.
In our particular scenario, we count 18 sequential colorless doublets in the quark sector, and 18 in the lepton sector. In the quark sector this includes the lightest doublet, H 0 , associated with the top quark, (t L , b L )H 0 t R (wherẽ H 0 is the charge conjugated Higgs), which is identified with the SM Higgs doublet. This establishes the universal quark sector Yukawa coupling to be g = 1. The quark masses and mixings are then determined by the spectrum of the 18 quark-sector Higgs doublets. We will see the Higgs doublet associated with the b-quark, g(t L , b L )H b b R , is expected to have a mass 3.5 TeV. We also have 18 doublets in the charged lepton sector. With the exception of H τ these tend to be much heavier as lepton masses are small. This framework provides three alternatives for the neutrino mass generation. Interestingly, neutrinos could be Dirac or Majorana via a type I or II seesaw mechanism.
Many of the scalars providing quark and lepton masses are well beyond present collider reach. The lighter ones, associated with the b and c quarks, as well as the τ lepton may have thresholds in the sub-10-TeV range and thus accessible to the LHC and future upgrades and higher energy machines. However, the heavier states may leave indirect imprints on flavor-changing observables in the quark and lepton sectors.
Scalar Democracy for us is new dynamics with subcritical compositeness due to a universal interaction in the far UV. We thus blend a few key ideas from compositeness and mainly emphasize that the Yukawa couplings are universal at Λ, subject only to renormalization group evolution from Λ to the weak scale. This, together with the input masses and mixings of the Higgs bosons, flips the quark and lepton flavor problems away from the issue of understanding a fundamental Yukawa-coupling matrix. The puzzle of fermion mass hierarchy becomes one of understanding and disentangling the multi-Higgs mass spectrum. This conversion of the flavor problem is an interesting exercise in its own right.
It should be emphasized that it is hard to understand the small Yukawa couplings in the SM, such as g electron ∼ 10 −6 . This cannot be generated perturbatively from zero, owing to the chiral symmetry of the electron. It is natural that g electron ∼ 1 in its coupling to a new Higgs H e , but through mass mixing, the induced coupling to the SMH becomes small. This is a key motivation for a scheme such as the one presented here.
As we have emphasized, the present scheme does not provide any explanation of the multiple Higgs masses and mixings. We will not concern ourselves with the overall naturalness, and treat symmetry-breaking effects as inherent in the Higgs mass terms. This is analogous to "soft-symmetry-breaking" in chiral Lagrangians and resonance models of the 1960's. We do impose fine-tuning constraints in the low-energy effective theory of H b , where the mixing generates the largest feedback on the SMH, H 0 (level repulsion), and the negative M 2 of the SMH can arise from this effect.
It is our conclusion, at least at a first pass with various simplifying assumptions made along the way, that such a theory can exist. New phenomena may show up at a high luminosity and/or energy-doubled LHC, and certainly at a 100 TeV collider. H b in the lightest mass limit ∼ 1 TeV is already accessible to the LHC, while it would be seen in the higher mass range ∼ 3 TeV in upgraded LHC runs. In our opinion, a robust theoretical spectroscopy of multiple Higgs bosons offers a rationale for luminosity and energy upgrades of the LHC and future ultra energetic machines. If the upgraded LHC were to fail to discover a pair of isodoublet Higgs bosons with universal coupling, such as the lightest H b or H τ states in our model, then this scheme would be disfavored.
Analyzing, at least schematically, the phenomenological consequences and constraints of the Scalar Democracy hypothesis is our main goal. We begin in Section II with a theoretical "motivation" for this perspective. Then in Section III, we summarize the dynamics of our model at low energies and arrive at a fairly simple effective Lagrangian describing the couplings to quarks and leptons of the multi-Higgs spectrum. This is followed by a more detailed discussion of phenomenology in the subsequent section, Section IV. Here we discuss the main production channels and collider prospects for the dormant Higgses. We will also discuss the implications of the dormant Higgses on flavor physics and the resulting bounds on their masses.
A reader interested in a summary of the observable features of the model, including neutrino masses, may skip directly to Section IV.

II. THEORETICAL MOTIVATION
A. Are There Many Scalars in Nature?
Gravity is a universal attractive interaction. All pairs of fundamental fermions must have attractive gravitational scattering amplitudes. Near the Planck scale this may involve exotic, new strong dynamics, new condensates, instantons, etc., and perhaps a new way to generate hierarchies.
For example, enhanced gravitational interactions may trigger condensates, dynamically generating Majorana masses for the right-handed neutrinos [11]. Similar effects may arise through gravitational instantons [12]. Alternatively, intriguing extensions of gravity, such as brane-world models, extra dimensions [13], or "bigravity" [14,15], have ingredients that may likewise produce universal attractive interactions at various scales. We can also generate a large subset of these scalars by postulating new, strong gauge dynamics, which is more concrete but will not be developed presently.
We assume that a universal attractive interaction generates bound state scalar fields at a high energy scale Λ (which may be of order M Planck , but could be lower). Our hypothesis is general and transcends a wide class of possible models. While we invoke a universal pairing force, such as gravity, this is nonetheless a schematic proposal. However, when we tie this to the SMH it becomes predictive.
We suppose the new scalar bosons are field-theoretic bound states of pairs of SM fermions. In analogy to the Nambu-Jona-Lasinio (NJL) model [16], the constituent fermions are free (unconfined) and any fermion pair couples with a common Yukawa coupling g at scale Λ to its bound-state scalar (the usual non-relativistic intuition does not apply.). If the new interaction is medium-strong (subcritical), then bound state scalars will form with masses which will be below Λ and positive. The formation of these bound states does not break any symmetries since a composite scalar inherits the quantum numbers of its constituents. Hence the scalars, in the absence of mixing effects, are presumably degenerate with mass M 2 < Λ 2 and cannot break flavor or gauge symmetries which would trigger proton decay, etc.
We do not have a theory of the scalar mass spectrum, in particular, and do not provide a fine-tuning mechanism to generate a large hierarchy. Moreover, we require symmetry-breaking and mixing effects to further empower the scenario. All of these effects will be considered to be "soft-symmetry-breaking" or "relevant" operators, and simply inserted by hand.
While the dynamics determining the spectrum is unknown, we know that it must respect the SM gauge symmetries. Rather than trying to concoct a theory of masses and mixing angles amongst a large number of Higgs doublets over a large range of scales, we will assume such a theory exists and let phenomenology guide us.
We connect the hypothesis with the SM by postulating that the SMH is the lightest scalar doublet. This would then be the (t L , b L )t R bound-state element of the scalar complex [5,6]. The top-quark Yukawa coupling, g, is then the universal coupling for all quark pairs to their particular Higgs doublets. Leptons will likewise have a universal coupling to their Higgs bosons, g , but we expect that the renormalization group (RG) running will yield a ratio g /g ∼ 0.7, due mainly to QCD (this is the analogue of the SU(5) relation for m b /m τ [17]). In fact g = 1 is close to the RG fixed point for the top quark [18,19], which is not far from the prediction of compositeness [6]. In fact, this scheme can bring the RG fixed point into concordance with m t = 173 GeV (see Eq. (16) and discussion below).

B. Sketch of the Dynamics
The dynamics of the lowest-lying sequential Higgs bosons is rather simple. Let us anticipate the dynamical implications in the case of the top and bottom quarks, and the three generations of neutrino masses.
The SMH, H 0 , in isolation has the usual potential: where −M 2 0 is negative, and phenomenologically M 0 88 GeV, λ 1/4 (the observed physical Higgs boson mass is √ 2M 0 125 GeV). New sequential Higgs doublets, H x , are "dormant," meaning they have the usual SMH electroweak quantum numbers, but owing to large, positive mass terms, M 2 x H † x H x , they do not directly undergo condensation. However, in order to generate the light quark and lepton masses, they must have small mixings to the SMH, parametrized by µ x . The vacuum expectation values (VEV) of the dormant Higgses are small, so the quartic terms should generally be negligible, and we ignore it. The mass mixing causes each H x to acquire a small VEV ("tadpole") of order where v 174 GeV is the electroweak VEV. Due to the mixing H 0 is then "level-repelled" down by an amount or order µ 4 x /M 2 x : We refer the interested reader to Ref. [20] for a similar implementation of the mechanism in supersymmetric models. Therefore, starting from a positive mass term, M 2 H , in Eq. (2), this cumulative effect can explain why the SMH boson has the tachyonic, or negative, −M 2 0 , in Eq. (1). Hence, in this scenario there is only a single condensate, associated with a conventional "Mexican hat potential," i.e., the SMH. The rest of the sequential Higgs bosons remain approximately pure doublets acquiring small tadpoles via mixing, which is generally suppressed by 1/M 2 x . Let us presently anticipate the main discussion and illustrate how this setup operates for the top-bottom sub-system. There we have the Yukawa interactions: The top mass determines the common Yukawa coupling to be g = m t /v 1. The b-quark then receives its mass from H b . By assum- while the Higgs boson mass is 2 In order to avoid fine-tuned cancellations between the two terms on the RHS of the above equation we anticipate that providing us with the estimate for an upper "naturalness" bound on M b : 3 This is not a firm upper bound, as it is based on a tuning argument. Another possible way of parametrizing the low 2 Here and throughout the paper we will use Mx to denote Higgs masses and mx to denote fermion masses. 3 Note that, if we assumed M 2 H = 0, then we could generate all of the tachyonic Higgs mass from the mixing with H b and we would predict: M b = (88 GeV)mt/m b = 3.36 TeV. scale fine-tuning would be to assume that all µ parameters are at the electroweak scale, for instance µ = 100 GeV. We will use these two tuning criteria later as benchmarks for comparison. Furthermore, we can certainly have a lighter M b with less fine-tuning. We note that a mass M b ∼ 380 GeV was previously obtained in a scheme in which the H b with g = 1 drives a Coleman-Weinberg potential for the SMH, and this has not to our knowledge been ruled out by the LHC [4].
We will find that this simple scheme is consistent with flavor-physics constraints. Interestingly, flavor dynamics is not then a fundamental consequence of the Yukawa coupling matrices as in the SM, but rather the low-energy suppression of flavor-changing neutral currents (FCNCs) is a consequence of the heaviness of these states-we also find some natural limits in which the structure of the theory is simplified and in which the phenomenological constraints are easy to understand.
The lightest doublets beyond SMH are H b associated with the b-quark, H τ associated with τ , and H c associated with charm (and possibly H ct and H tc , depending on how the CKM matrix is generated, as we will see later). Depending on mixing assumptions, H τ could in principle be lighter than H b . We thus encourage LHC experimentalists to consider searching for these objects, and we discuss collider phenomenology in Section IV A and Section IV B.
At the other end of the spectrum, we have the physics of neutrinos and charged leptons. In the Scalar Democracy, there are three alternatives for the neutrino mass mechanism. For example, where neutrinos are Dirac their masses are generated like any other fermion in the present model. The Yukawa couplings for the neutrinos are similar to those for quarks and leptons, and the appropriate neutrino mass terms are generated by mixing of the dormant and SM Higgses. We generate neutrino mixings and small neutrino masses by having ultra-large and positive M 2 for the neutrino-Higgs fields in the matrix, e.g., For neutrino masses of order 10 −2 eV we thus have M να ∼ 10 11 GeV, too heavy for any observable collider or flavor signature. Notwithstanding, the Scalar Democracy scenario may also lead to Majorana neutrinos via type-I and type-II seesaw mechanisms. In the type-I seesaw case, a light sterile neutrino is predicted, with mass in the keV−GeV range. A general discussion of all these possibilities will be provided in Section V A.

C. Scalar Democracy and Counting Scalars
If we assume one boundstate per fermion pair at some high scale Λ, then we can count the number of composite scalars in the theory.
All of the SM matter fields can be represented by 48 two-component left-handed spinors, ψ i A . This includes all the left-handed and anti-right-handed fermions. We can collect these into a large global SU(48) × U(1) multiplet, corresponding to the global symmetry assuming that the new dynamics are blind to the SM gauge interactions. We emphasize that this is a dynamical symmetry, and familiar GUT theories that contain only the SM fermions will be gauged subgroups of this SU (48). Here the indices (i, j) run over all the 48 flavor, doublet, and color degrees of freedom of the SM fermions. The most general scalar-field bilinear interaction we can construct of these fields is where Θ ij transforms as the symmetric 1176 representation of the SU(48) × U(1) (this is analogue to the sextet representation of SU (3)). The field Θ ij contains many complex scalar fields with assorted quantum numbers, including baryon and lepton number, color, and weak charges.
To make contact with the SM fields, we consider the usual 24 left-handed quarks and leptons, Ψ L,i , and the 24 right-handed counterparts, Ψ R, i . The index i now runs over the chiral SU(24) L and i over the chiral SU(24) R subgroups of SU (48). With this notation, we can construct three interactions with bilinear fermion fields, (12) where Φ i j is the (24 L , 24 R ) complex scalar field with 24 2 = 576 complex degrees of freedom. Ω and Ω are the symmetric 300 representation of SU(24) L and SU(24) R respectively 4 . Thus we have Φ (576)+Ω (300)+ Ω (300) = Θ (1176), matching the degrees of freedom of Θ ij .
Here Ω ij and Ω ij are the analogues of Majorana masses and carry fermion number, while Φ contains fermion number neutral fields, such as Higgs fields, in addition to (B − L) leptoquark multiplets. The Φ, Ω and Ω fields can be viewed as the "composite fields" arising from a NJL model effective description of the new forces. Consider just the SU (24) where the negative sign denotes an attractive interaction in the potential. It should be noted that we can equally well write current-current (and tensor-tensor) interactions, mediated by heavy spin-1 bosons (or Pauli-Fierz spin-2 4 Notation: If Ψ j L = 1−γ 5 2 Ψ j is a left-handed spinor transforming as a 24 under SU(24) L , then (Ψ j gravitons); these will generally contain scalar channels and will Fierz rearrange to effectively reduce to Eq. (13) with the attractive signs. There also exist the possibility of the following NJL models: which lead to the composite bosons Ω and Ω. Such universal master interactions may arise as sub-sectors of more general gravitational scattering amplitudes with many other effects near the Planck scale, M ∼ M Planck , including gravitational instantons or, at lower energy scales, e.g., from a strong bigravity force. Eq. (13) by itself can therefore be a starting point of a discussion of a dynamically generated extended Higgs boson spectrum. The first step to solving an NJL theory would be to factorize the interaction of Eq. (13) by introducing auxiliary scalar fields. This leads to the equation we started with, Eq. (12), where Φ, Ω and Ω are auxiliary fields. The universal interaction will bind fermion pairs into scalars that are bound states of ordinary quarks and leptons and will generate a plethora of Higgs doublets. These bound states will have a universal Yukawa coupling g at the scale M 2 . Moreover, with g taking on a near-but-subcritical value, these bound states will generally have large positive masses but can be tuned to be lighter than M .
Symmetry-breaking effects will be required to split the spectroscopy, including the SMH down to its observed negative mass term. All other doublets remain heavy, but will mix. The problem of solving an NJL model in the large-N fermion loop approximation, or equivalently by the RG, is discussed in detail in Refs. [6,7].
If g is supercritical then some or all multiplets will acquire negative renormalized masses, M 2 (µ → 0) < 0 and the theory develops a vacuum instability. For example, the field Φ ij with a supercritical coupling will generally condense into a diagonal VEV, Φ ij = V δ ij , and this would become a spontaneously broken Σ-model of SU(24) L × SU(24) R × U(1) × U(1) A . In this supercritical case, all the fermions would acquire large, diagonal constituent masses of order gV , inconsistent with observation.
The structure we have just outlined, even if subcritical, will contain many composite Higgs doublets. While we want to have a reasonably deep binding of these scalars, the system must be near-to-but-subcritical such that no large VEVs will form. We would then expect a positive, diagonal mass-squared matrix amongst the many composite scalar states and all fermions would be massless. However, the effect of mixing, i.e., off-diagonal mass terms, can arise from "extended interactions" in analogy to "extended technicolor" or latticized extra dimensions.
Exactly how the scalar mass spectrum is generated is beyond the scope of our present discussion. We will simply assume such a spectrum of masses and mixings between the bound state scalars that allows for a light sector from the SMH to multi-TeV scales exists and extends up to the highest scales. We assume that Ω ij , Ω ij and the color-carrying weak doublets have very large positive M 2 and therefore we will ignore them.
Let us examine the quantum numbers of the spectrum of states in the Φ ij system. There are 24 × 24 = 576 composite complex scalars, and these devolve into the following states upon gauging the fermions: where the brackets denote the SM quantum numbers. The first four entries in the above list are the 36 Higgs doublets in the quark and lepton sectors respectively.
If we consider the 18 scalars in the quark sector, ignoring their masses and EW charges, we will have a Yukawa interaction at the scale Λ that is SU(6) L ×SU(6) R invariant and of the form: where Ψ = (u, d, c, s, t, b) and Σ is a 6 × 6 complex matrix composed of 18 doublets. The renormalization group equation for g is then determined to be at one-loop order, where κ includes the smaller electroweak corrections (which breaks the SU(6) L × SU(6) R invariance). This describes the running of g down to scales at which the various Higgs doublets decouple. Likewise, we have an equation for the lepton sector, g → g , where the −8g 2 3 term is dropped. Assuming the Planck mass corresponds to a Landau pole in g and that all Higgs bosons are active down to the electroweak scale, and then we derive g 0.93 and g 0.71 (this also leads to some "fine-structure" as the electroweak terms in κ split degeneracy of between the Yukawa couplings for the up-and down-type quarks; the full details of this are beyond the scope of this paper).
This result implies a top quark mass of 161 GeV. This is the prediction of the modified RG fixed point (equivalent to a focus point) of [19] including additional Higgs bosons, and represents a significant improvement over the original minimal top condensation models [6]. This prediction is robust with respect to the precise values, Λ and g(Λ). If we include the masses of the heavier Higgs bosons (as discussed below) and decouple them at their thresholds, the prediction will increase and we expect it to converge on the observed top quark mass.
Grand unification is greatly complicated by such dynamics. We relegate such an investigation to future work.

III. LOW ENERGY EFFECTIVE THEORY
The observed Higgs boson must reside amongst the color singlet Q L q R doublets. We proceed under the simplifying assumption that the SM Higgs doublet, H 0 , can be identified with the doublet that couples to the fermionic combination of a top quark pair This is the unique logical choice, as it has the largest Yukawa coupling in the Standard Model and our theory dictates that all quarks will have this universal coupling. This therefore recovers in part the top-condensation models [5][6][7]. We further reduce the scope of the problem by assuming all of the dormant Higgs doublets apart from the 36 color singlet Q L q R and L L R doublets are arbitrarily heavy and therefore decoupled from the low-energy effective theory. Although the phenomenology of these other scalars could be very interesting, this framework does not provide any insight on their mass scale. Besides, we will assume for concreteness that the neutrino mass mechanism is the same as for charged fermions, and comment on alternatives in Section V A.
At this point, it is convenient to resort to a common notation for the individual Higgs bosons. Below the scale Λ, the effective theory we are considering is a modification of the SM where the Higgs sector has been replaced by The generational indices of the fermions are labeled by indices i, j and the primes indicate that we are working in the gauge eigenbasis. The Higgs doublets are denoted by H f ij in their flavor eigenbasis, where f = u, d, ν, e represents the fermion type which acquire mass from the doublet and i, j the generations they couple to. Each doublet has an upper charged component and a lower neutral component, H x = (h + x , h 0 x ) and we employ the charge conjugation conventionH = iσ 2 H * . All fermionic masses and mixings are due to the VEVs of the Higgs bosons. Universality implies that all quarks have a Yukawa coupling g 1 and leptons g 0.7.

A. The Higgs Potential
Formally, we can define the 36 Higgs doublets as a "vector," where we separate out the (mostly) SM-like Higgs: (H 0 , H a ). Thus, in addition to the SM-like Higgs, H 0 , we have 35 doublets represented as, H a = (H 1 , . . . , H 35 ). Here, we define new notation and apologize for that. Nevertheless, this index notation will only be applied in this subsection. In the remainder of the paper, we revert to the notation H f ij . We may write out the multi-Higgs boson potential as where the positive quartic interactions of the dormant Higgses are taken to be negligible due to their large masses 5 . Thus the Higgs mixing is determined by the mass matrix. The mass terms that mix H a with H 0 can be viewed as vector, µ 2 a = µ 2 1 , ..., µ 2 35 , and M 2 ab is a 35 × 35 hermitian mass matrix amongst the heavy scalars, which is taken to have positive eigenvalues.
It is important to note that there are two distinct types of mixing in this model sourced by µ 2 a and the off-diagonal part of M 2 ab respectively. The former is responsible for mixing the SMH into all the Higgs flavor states, while the latter introduces mixing between the dormant Higgses. We will regard the off-diagonal elements, µ 2 a , as perturbatively small when compared to M 2 ab . The negligible quartic couplings leads to a mass degeneracy of the charged and neutral Higgs components of the dormant doublets. Therefore, retaining the symmetrybreaking mass terms, the dormant Higgs fields will acquire tadpole VEVs via the µ 2 a H † a H 0 mixings, but their fluctuating fields will be degenerate doublets.
The Higgs mass terms are diagonalized at leading order (LO) in M 2 H , µ 2 a /M 2 ab by going to the basis defined by where H 0 , H a defines the physical Higgs doublets before EWSB. After rotating away the µ a mixing, the equations of motion for the VEVs of the Higgs fields read where is the mass term for the SMH. The Higgs tadpoles are thus given by where v 2 = −M 2 0 /λ (174 GeV) 2 is the SM VEV. From our assumptions of the Higgs potential, we ignore further mixing between the massive Higgs bosons. The dormant Higgses constitute degenerate doublets and H 0 is directly identified with the SMH, its components being the physical Higgs and the 3 EW Goldstone Bosons. Eq. (20) shows how H 0 is mixed into all the Higgs flavor states. Accordingly, the dormant Higgses acquire small VEVs which in turn produces the fermion masses in our model.

B. Fermion masses and mixing
The fermions acquire their masses from the Yukawa couplings as shown in Eq. (18). In general, the mass matrices for the fermions are determined by the VEVs of the corresponding Higgs flavor eigenstates; where g f is the Higgs coupling to quarks (g) or leptons (g ), L f , R f are unitary matrices which rotate left-and right-handed fermions respectively from the gauge to the mass basis. We will proceed to determine the physical couplings between the fermions and the new Higgses assuming the Higgs potential outlined in the previous section and furthermore take the limit where the mass matrix, M 2 ab , of the dormant Higgses is diagonal thereby eliminating dormant-Higgs mixing.
First, we focus on the quarks whose mass eigenstates are given by The couplings between the down-type quarks and the neutral components of the dormant Higgses, in the mass eigenbasis, are given by 6 6 For convenience we have identified e.g. H d 23 = H sb . We apologize for the notational inefficiencies. and likewise for the up-type quarks arises from the feedback on the SMH. Similarly in the lepton sector we set the mass eigenstates to be which gives interactions with the neutral Higgs bosons of the following form: A similar construction for the coupling of the fermions to the charged Higgses follows straightforwardly. We observe that h, the neutral component of H 0 , is completely indistinguishable from the SMH. From the previously outlined assumptions on the Higgs sector we are able to estimate the masses of the dormant Higgses. Using Eq. (24) and Eq. (20), along with diagonal dormant mass matrix, we find that the fermion mass matrix is related to the Higgs mass terms by To avoid a large fine-tuning in the SMH mass term, we expect µ 4 /M 2 (100 GeV) 2 . We thus arrive at the estimate We therefore expect the dormant Higgs masses to be inversely proportional to the corresponding entry in the fermion mass matrix. Once we have made an ansatz for the fermion rotation matrices, L f , R f , this provides us with an estimate for the Higgs masses. The bound shown in Eq. (32) should merely be viewed as a guideline, as it is based on a fine-tuning argument. Before we proceed to review two important limiting cases in the next section, we first define our notation: the Cabibbo-Kobayashi-Maskawa (CKM) matrix is given by V CKM = L u L † d , and likewise the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix is given by V PMNS = L ν L † e . In the limit of negligible dormant Higgs mixing, the required input to determine the Yukawa couplings is the set of mass rotations, L f and R f . These rotation matrices are implicitly generated by the unknown µ 2 a which gives rise to the VEV structure we have encoded into the above expressions. Examining Eqs. (26) and (30), we find that the individual matrices L f , R f are unobservable in the Standard Model. However, their combinations form the mixing matrices, V CKM and V PMNS , which enter in the charged currents of the quark and lepton sectors respectively. Given sufficient dormant Higgs data, the structure of L and R could be determined, but at present cannot be determined directly. We must therefore make an ansatz for their form.

C. R f = 13 and No-Mixing Limits
We proceed by considering the various symmetries that arise in the special limit R f → 1 3 and m q = 0. These are useful for relaxing rather stringent flavor mixing constraints.
Consider Eq. (26) with R d = 1 3 and the strange quark mass set to zero, m s = 0. Furthermore, decompose the mass matrix of the dormant Higgs fields, H a , as a sum of a diagonal M 2 diag and an off-diagonal δM 2 Hermitian matrices: In the case δM 2 = 0, we find there exists a discrete symmetry 7 , e.g., reflection of the right-handed strange quark, s R → −s R and the corresponding Higgs fields (H bs , H s , H ds ) → −(H bs , H s , H ds ). This symmetry is a generalization of the Glashow-Weinberg symmetry [21], and in our present case of many Higgs bosons this is restrictive.
First, we observe that in general H s will mediate an interaction of the form, in the gauge basis. In the case R d were not unity, this would contain mixed combinations such as (35) in the flavor basis, which results in ∆S = 2 transitions. The K L -K S mass splitting places severe limits on M s > ∼ 10 3 TeV for left-and right-handed mixings of the order of the CKM matrix, while the mass estimate of Eq. (32) suggests M s 100 TeV. With the discrete symmetry, R d = 1 3 , no such interactions are generated at tree-level, and the tension is substantially alleviated. Moreover, the discrete symmetry forbids the similarly dangerous Higgs mixing term δM 2 ds,sd H † ds H sd , which mediates interactions such as 1 directly. At the one-loop level the Higgs-Higgs box diagrams, with exchange of the full doublet H s , produces the effective interaction This operator is not restricted by R d = 1 3 , but yields a ∆S = 2 operator for nontrivial left-handed rotations, L q , Comparing this effective operator with kaon-mixing bounds on the left-handed current [22], we arrive at the limit M s 60 TeV, which is compatible with the mass estimate.
Our model can thus generate CKM mixing, with multiple flavorful Higgs doublets and yet present no large FCNC (no tree-level O(1/M 2 x ) operators) if we make R = 1 3 . Typical cases include small breaking effects, giving R = 1 3 , and flavor physics remains an important probe in a system like this.
In getting a sense of the model, it is useful to consider the limit in which L f , R f are set to unity with δM 2 ab = 0 which turns off the mixing among dormant Higgses at leading order. We still allow non-zero µ 2 a terms that mix the dormant Higgses to the SMH.
Let us focus on the up-quark system c.f. Eq. (27). The top quark has acquired its mass by direct coupling to the SMH, which defines the universal quark Yukawa coupling g = 1. In this limit we can identify the fields that develop tadpole VEVs and give the light quark masses, Eq. (20): and likewise for the d, e and ν sectors. In this limit we can compute the diagonal dormant Higgs masses. All  other fields H x have no VEVs in this limit. The full set of dormant Higgs masses are tabulated in Table I under two different case assumptions: 1. The feedback on the Higgs mass is maximal, but limited to (100 GeV) 2 for each fermion (thus choosing a larger value of µ for lighter quarks), hence M q = (m t /m q ) · 100 GeV and M = (g m t /m ) · 100 GeV.
2. µ = 100 GeV for each SMH mixing term, M q = (m t /m q ) 1/2 · 100 GeV and M = (g m t /m ) 1/2 · 100 GeV, which makes only the lightest H b , H τ , H c have any significant feedback on the SMH mass.
Hence, in the zero mixing angle limit, it is easy to see how the flavor problem is mapped into an inverted mass spectrum of Higgs doublets.

IV. QUARK-SECTOR PHENOMENOLOGY
In Section IV A, we begin by asking if relevant current LHC searches can place limits on the, possibly multi-TeV-scale, lightest dormant Higgses. We follow with a discussion of the prospects of discovery both at the LHC and at future colliders of higher luminosities and center of mass (COM) energies in Section IV B. Throughout Section IV A and Section IV B, we calculate the leading-order (LO) production cross section of the NP (new physics) processes using Universal FeynRules Output (UFO) format of FeynRules [23] implemented into MadGraph5 aMC@NLO [24] and applying the NNPDF2.3 PDF set [25]. We implement the five-flavor scheme in order to account for the intrinsic b-quark content of the proton. In addition, constraints on the dormant Higgs masses from meson mixing and dedicated searches for flavor change in the charged lepton sector are outlined in Section IV C and Section V B respectively.

A. Current Limits from LHC Searches
In this section, we investigate if current limits from the LHC can place meaningful constraints on the masses of the Higgs spectrum of the Scalar Democracy. The dormant Higgs most accessible at the LHC, is the lightest (non-SM) doublet which couples to the b-quarks, namely Let us first focus on the neutral component, h 0 b . For single h 0 b production, the dominant contribution is via b-quark fusion as shown in Fig. 1. We remind the reader, that the coupling of the b-quarks to h 0 b is g = 1. This implies that h 0 b will decay approximately to two bquarks with a branching ratio of one. As a consequence, a simple comparison of a search for new resonances decaying into jets containing b-hadrons is particularly amenable for us to investigate if the best suited current searches may be able to exclude the relevant (multi-TeV) mass regime of h 0 b . The ATLAS analysis we choose to recast corresponds to a dataset collected at √ s = 13 TeV and searches for two b-tagged jets with an invariant mass (m bb ) in the 0.57 − 5 TeV range [26]. In the "low" (0.6 ≤ m bb (TeV) < 1.25) and "high" (1.25 ≤ m bb (TeV) ≤ 5.0) invariant mass regimes, 24.3 fb −1 and 36.1 fb −1 of integrated luminosity were analyzed respectively. The primary aim of this analysis is to search for Z decays to b-quarks. Although this search is designed to place limits on a vector boson, rather than a scalar, the final states are the same in both cases so the difference in acceptance between our model and that of [26] should be small.
We note that single h 0 b production is also possible via gluon fusion mediated by a bottom loop 8 . However, due 8 The top loop will also contribute to this process because of the mixing of the SMH with  to the chiral suppression present in the loop mediated processes, these contributions are subdominant to the treelevel production by more than two orders of magnitude and therefore will be ignored presently.
The one-and two-sigma regions for the expected upper limit on cross section of di-b-jet production is shown in Fig. 2, as indicated by green and yellow respectively, while the observed upper limit, at 95% CL, is shown by the black circles. There is no significant deviation between observed and expected limits. The details of the analysis used to calculate these curves can be found in the aforementioned reference.
We have not performed a full analysis, including detector effects, to calculate signal cross section and therefore the uncorrected cross section will likely be a slight overestimate. To mitigate this issue, the LO production cross section of h 0 b was multiplied by a relatively stringent efficiency factor, | 2b−jet | = 0.2, to account for the two b-jet reconstruction efficiency. To justify this choice of efficiency, we refer to [26], which details the b-tagging efficiencies for both "low" and "high" mass regions. They found for the low-dijet-mass regions the efficiency for tagging two b-jets decreased from 0.5 to 0.2 as m bb was increased from 0.65 to 1.25 TeV. However, in the highdijet-mass region, the event tagging efficiency for two b-jets ranged from 0.4 − 0.05 for masses as m bb was increased from 1.25 to 5.0 TeV. The total inclusive cross section for charged and neutral Higgs production in heavy quark annihilation has been calculated at next-to-next-to-leading (NNLO) accuracy in QCD [28]. For masses of new Higgses of 600 GeV (which couple solely to b-quarks), the NNLO corrections are small, σ NNLO /σ LO ∼ 0.95. Moreover, it is likely this K-factor will remain small for larger Higgs masses and therefore we did not calculate higher order corrections to the cross section.
Although, our reinterpretation lacks the sophistication of a full experimental analysis, it allows for a rough comparison of our signal cross section to a currently available dataset. From Fig. 2, it seems likely that the LHC could exclude the lower mass region, M b 1.0 TeV.  Expected ± 1σ Expected ± 2σ Figure 2. The green (yellow) region shows the ±1σ (±2σ) expected limit for number of events from Z decaying to two b-quarks. We note that the observed cross section (within 95% C.L) is within 2σ of the expected cross section. The blue lines shows the predicted LO inclusive cross section of h 0 b production. The observed and expected data was analysed by an ATLAS group [26] and this data was made available via HEPData [27].

B. Prospects for Future Discovery
The next most pertinent issue is the prospect for future discovery of h 0 b and the required integrated luminosity. In the left plot of Fig. 3, we show the LO production cross section for single h 0 b production as a function of h 0 b mass for three COM energies: 13, 26 and 100 TeV. As masses below 1 TeV are disfavored by data and the mass regime of ∼ 3.5 TeV is favored theoretically, we shall discuss the latter. For masses M b = 3.5 TeV, the corresponding cross section is ∼ 2.4 × 10 −4 , 0.011 and 2.0 pb for √ s = 13, 26 and 100 TeV respectively. In the right plot of Fig. 3 we display the LO cross section as a function of COM energy for three fixed masses: M b = 1.0, 2.0, 3.5 TeV as shown in purple, green and red respectively. Naturally, the cross section increases for higher COM energies.
The most crucial plots in addressing the question of observability are the two plots shown in Fig. 4. The quantity of interest is the significance, defined as S/ √ B where S is the expected number of signal counts and B is the expected number of background (BG) counts. The left and right plots display the significance, S/ √ B, as a function of integrated luminosity for COM energies √ s = 13, 26 TeV respectively. The significance is defined to be where L is the integrated luminosity and | 2b−jet | = 0.2 is the applied reconstruction efficiency for two b-jets which we assume is the same for both signal and SM background.      In order to estimate the SM background, we applied a number of cuts on the di-b-quark production: where p T (b) is the transverse momentum of each of the b-quarks, m bb is the invariant mass of the b-quark system, Γ h 0 b = 3M b /8π is the width of h 0 b and y is the rapidity. We note that the cross sections have been calculated at the level of the hard matrix element. The above cuts were applied to the invariant mass system of the two b-jets of the BG. We did not apply the cuts to the signal process but have chosen them in such a manner as to not significantly reduce the signal strength. In order to approximate the effect of the cuts on the signal, we assumed a Breit-Wigner distribution for the differential cross section as a function of the invariant mass such that the chosen range captures approximately half of the signal events. Subsequently, we multiplied the signal by a factor of 0.5, c.f Eq. (39).
In both plots of Fig. 4, the S/ √ B value is shown for three fixed masses: 1.0, 2.0, 3.5 TeV (as indicated by purple, green and red colors respectively) and at two COM energies: 13 and 26 TeV (solid and dashed lines respectively). The integrated luminosity takes the range between 10 2 − 3.5 × 10 6 pb −1 . The value S/ √ B = 5 is indicated by the dashed black lines and provides an estimate of the required integrated luminosity for the discovery of h 0 b . For current COM energies, the masses, M b = 1.0, 2.0 TeV would be discoverable at integrated luminosities of ∼ 2 × 10 3 and 2 × 10 5 pb −1 respectively. However, for the heavier mass M b = 3.5 TeV, this would require greater than 3.5 × 10 6 pb −1 of integrated luminosity. With higher COM energies, such as 26 TeV, the masses M b = 1.0, 2.0, 3.5 TeV would all become accessible at ∼ 5 × 10 2 , 10 4 , 2 × 10 5 pb −1 respectively. Although we have not shown the projections of the significance for a COM collider with √ s = 100 TeV, we have found that h 0 b of all three masses can be discovered with less than 3 ab −1 of integrated luminosity.
In the event that h 0 b were discovered, a search and discovery of its charged counterpart, h ± b , would complete the doublet, H b . Although, h ± b is mass degenerate with h 0 b , its production cross section is smaller by at least an order of magnitude, because it couples to tb or bt and therefore cannot be produced from heavy quark fusion. The leading production channel of h ± b is in association with a bottom and a top quark, where the final state will be two b-jets and two further jets. Such a search would be more feasible than di-h ± b production, which is mainly mediated via photons or off-shell Z-boson, as this process is kinematically suppressed from the two Higgses in the final state.
One could also consider the production of the lightest dormant doublet associated to the lepton sector, H τ . Due to its large coupling to τ , g = 0.7, the neutral h 0 τ decays dominantly to τ τ . As H τ does not couple strongly to quarks, it can only be produced via electroweak interactions, making it challenging to probe at the LHC. In spite of the fact the τ reconstruction efficiency is slightly higher than that of b-quarks (| τ | ∼ 0.8 as measured at √ s = 13 TeV [29]), as the h 0 τ production is electroweak in nature, the cross section is significantly (several orders of magnitude) suppressed compared with the analogous production of h 0 b . Although this discussion has been largely schematic, and certainly an in-depth analysis would be needed to search for the lightest dormant Higgs of the Scalar Democracy, we find it an encouraging first step in demonstrating that the lightest dormant Higgs is within reach at current energies and most certainly future upgrades.

C. Meson-antimeson mixing
Having discussed the main production channels for the dormant Higgses, we now turn to the implications of the Scalar Democracy for flavor phenomenology in the simple limit discussed in Section III B, where the dormant Higgses do not mix with each other. The NP parameters are reduced from several hundreds in the generic model down to a more manageable set consisting of the dormant Higgs masses and the fermion rotation matrices, L f , R f , once we impose the requirement that SM fermion masses must be reproduced.
Potentially large FCNCs may arise when the quark fields are rotated from their gauge to their mass eigenbasis and will place stringent constraints on the mass scale of these dormant Higgses. The constraints from flavorviolating processes in both the quark and lepton sectors have been studied in depth the context of two Higgs doublet models, see e.g. Ref. [30].
Integrating out the heavy Higgs mass eigenstates gives NP contributions to the effective four-fermion operators at the low energy scale relevant for flavor physics. The neutral Higgses induce new four-fermion operators already at tree-level, contributing to the effective low-energy interactions among fermions with flavor indices (q, p, r, s), namely (40) The sum runs over all fermions except the t-quark denoted by (u,33), in order to avoid double counting, as this would be the SMH contribution. Besides, the SMH couplings are aligned with the fermion masses and thus do not lead to tree-level flavor changing transitions. The tree-level exchange is in most cases expected to be the leading NP contribution to the ∆F = 2 operators, as it is only suppressed by the mass of the dormant Higgses mediating the process.
The right-handed rotation matrices, R f , are not constrained a priori as they are unphysical in the SM. The limit R f = 1 3 is thus viable and turn off the ∆F = 2 operators in Eq. (40), as discussed in Section III C. Near this limit the leading NP contributions will not come from tree-level contributions, but will rather be induced at the one-loop level from box diagrams with neutral or charged Higgses with or without SM charged currents W bosons.
We will focus on meson-antimeson mixing which are ∆F = 2 observables that are loop and GIM suppressed in the SM. These observables are thus sensitive to NP degrees of freedom which propagate in the mixing amplitudes. As such, signals of NP may be constrained through comparison of observed meson mixing and the SM predictions.
An exact analysis of the flavor constraints are somewhat beyond the scope of this paper. We will content ourselves with a rough estimate for the allowed mass range for each of the dormant Higgses. To this end, we apply the bounds of Isidori [22] on the Wilson coefficients for the scalar exchange effective operators of Eq. (40), cf. Table  II. To assess the parameter space, we will probe the NP constraints using a number of benchmark cases for the rotation matrices. In each case, flavor bounds will place a lower limit on the masses of the dormant Higgses, and we will compare these to the expected masses obtained from the fine-tuning arguments c.f. Eq. (32).
The meson-antimeson mixing bounds on NP contribution to the ∆F = 2 operators, can be translated into bounds on the dormant Higgs masses once L q and R q are assumed. For each Higgs we determine the mass where it saturates each of the bounds in the absence of contribution from other dormant Higgses; the highest such mass is used as a lower bound for the corresponding Higgs mass. This is a somewhat simplified approach to the constraints, as it discounts possible interference between contributions from multiple Higgses, not to mention the running of the effective operators below the EW scale. Since we are only exploring a small region of parameter space of the full model in any case, the bounds we obtain using these approximations remain a good indication for the scale where the Higgses become visible in flavor physics.
The benchmark points we consider are informed by the constraint that the left-handed quark mixings must reproduce the CKM matrix. Hence, it is unavoidable to have some left-handed mixing, albeit there is still freedom to put all of the left-handed mixing into either the up or the down sector. A good middle ground would seem to be taking a half and half approach. One will quickly find that with this assumption for the left-handed rotation, the meson-antimeson mixing constraints will favor small mixing of the right-handed quarks. We consider the following benchmarks: is an approximately unitary matrix. For each benchmark we have determined the lower bound individually for each dormant Higgs mass, for all the effective operators in Table II to satisfy their bounds. The resulting lower bounds are plotted in Fig. 5. As a comparison the gray lines in the figure show the expected masses of the Higgses assuming the absence of fine-tuning, cf. Eq. (32). We proceed to investigate the benchmark points one at a time.
Benchmark 1 shows the mass spectrum of the theory in the absence of right-handed rotation. The Higgs mass estimates in Fig. 5 exhibits the inverse hierarchy of the quark masses detailed in Table I with the Higgses associated to the lighter quarks tending to be the heaviest. The non-trivial left-handed rotation will lower the mass estimates for the off-diagonal mass. In particular, H ct , must be almost as light as H b in order to generate a large off diagonal mass element in the fermion mass matrix. The corresponding plot in Fig. 5 does not show any bounds on the Higgs masses, as we have not included the loop-induced effects here. Benchmarks 2 and 3 illustrate a marginal case, for how much right-handed rotation can be allowed without tension between mass estimates and flavor bounds. The matrix V (η 1 , η 2 ) used to parametrize the right-handed matrix is applied for both R u and R d for simplicity. The (1, 2) entry governs the mixing of the first two generations of quarks. As η 1 is increased as do the flavor bounds on H u , H uc , H cu and H c , because they are governed by D 0 mixing (between the u and c quarks). Similarly the bounds on the H d , H ds , H sd and H s from K 0 are pushed up, though they are less problematic. The entry (1,3) of V (η 1 , η 2 ) governs the bound on H b from B s mixing in a similar manner. Going from the marginal case η 1 = η 2 = 1 in Benchmark 2, to the case η 1 = η 2 = 3 in Benchmark 3, we observe how all the corresponding flavor bounds are increased by a factor √ 3, giving a slight tension between mass estimates and flavor bounds. We would like to emphasize that this tension only exists in the complete absence of fine-tuning of the Higgs mass, −M 2 0 . If one were to tolerate fine-tuning to the percent level, the Higgs mass estimates would increase by a factor 10. Benchmarks 4 and 5 parametrize R q in terms of powers of CKM matrix, to illustrate how large the right-handed mixings are allowed to be in terms of a more familiar matrix. Due to the relatively large mixing, (V CKM ) 12 (V CKM ) 21 0.23, one has to go to the case R † u = R d = V 1/3000 CKM of Benchmark 4, before the flavor bounds decrease below the mass estimates of H c . On the other hand, if one is willing to accept a percent-level fine tuning of the Higgs mass we can go to the case CKM of Benchmark 5. We stress that, the mass estimates shown in Fig. 5 all are with the complete absence of fine-tuning.
Benchmark 6 puts the entire left-handed rotation into L d , and gives an alternative example of evading the flavor bounds. In this case a large R u can be tolerated without generating D 0 mixing at tree level. The main constraints in this case, is from K 0 mixing.
Though our selection of a few benchmark points does not constitute an in-depth analysis of the flavor physics of the Scalar Democracy, it demonstrates that there is room for the framework to avoid tension with flavor physics. Yet, care has to be taken to ensure that one works in an allowed regions of parameter space.

A. Neutrino sector
In our framework, the mass spectrum of SM fermions is mapped onto a scalar mass spectrum where lighter scalars develop larger VEVs and thus leads to heavier fermion masses, see Eq. (32). As the masses of charged fermions have been precisely measured, a sharper prediction of the scalar spectrum can be made. We do not know the absolute value of neutrino masses, nor the mass generation mechanism, resulting in looser possibilities for the scalar bound states associated to neutrino masses. The goal of this section is to examine these possibilities.
Although the absolute values of neutrino masses are still to be measured, neutrino oscillations give us valuable information on their masses. From the KamLAND experiment and observations of the solar neutrino spectrum [34,35], the solar mass splitting 9 has been measured to be ∆m 2 21 ≡ m 2 2 − m 2 1 +7.4 × 10 −5 eV 2 . The sign of ∆m 2 21 is known due to matter potential effects on neutrino oscillations. Besides, atmospheric and accelerator neutrino oscillation experiments [36][37][38] are compatible with |∆m 2 31 | ≡ |m 2 3 − m 2 1 | 2.5 × 10 −3 eV 2 . Therefore, the following two neutrino mass orderings are still viable experimentally, m 1 < m 2 < m 3 or m 3 < m 1 < m 2 (recent experiments exhibit a slight preference towards the former scenario [39,40]). In addition, the sum of neutrino masses is bounded from above by 0.12 eV from cosmological observations [41]. We can thus identify three extreme scenarios for the neutrino mass spectrum (m 1 , m 2 , m 3 ): normal hierarchy (0, 0.008 eV, 0.05 eV), inverted hierarchy (0.05 eV, 0.05 eV, 0) and semi-degenerate neutrinos ∼ (0.03 eV, 0.03 eV, 0.06 eV). As the heaviest neutrino needs to have a mass above 0.05 eV, that will typically translate into an upper bound for some of the scalars (see Table III).
In the Scalar Democracy there are three viable alternatives for the generation of neutrino masses. The neutrino mass generation can be similar to the charged fermion mass generation, as outlined in Eq. (30). In this case, neutrinos would be Dirac fermions. Nevertheless, due to the quantum numbers of the lepton doublet, L L , and the right-handed neutrino, N R , two other bound states can contribute to neutrino masses, namely, S N ∼ N R N C L , and ∆ ∼ L L L C R . While H ν may lead to Dirac neutrino masses, an accompanying non-zero expectation value for S N would constitute a Majorana mass term for righthanded neutrinos, realizing a type-I seesaw mechanism, whereas ∆ may be identified as the scalar triplet in type-II seesaw models (see Refs. [42][43][44][45][46][47][48][49]). Thus, the mechanism of neutrino masses depends on which dormant Higgses acquire a VEV.
In the type-I seesaw realization, the relevant terms in the Lagrangian are (generation indices have been suppressed) and hence neutrino masses are given by Because S N is a SM singlet, it will not acquire its VEV from mixing with the SMH as in Eq. (2), but rather from the term µ S H † 0 H 0 S N , which gives, S N = µ S v 2 /M 2 S . Requiring electroweak values for the dimensionful parameter µ S = 100 GeV, and imposing constraints on singlet-Higgs mixing and active-sterile neutrino mixing leads to the allowed region in the (M S , M ν ) plane shown in Fig. 6. The gray regions are excluded due to mixings while the blue region predicts too heavy neutrinos 10 . The atmospheric mass splitting requires at least one neutrino to be heavier than 0.05 eV. A dashed black line corresponding to m ν = 0.01 eV is drawn to guide the reader. Notice that a sort of seesaw mechanism between the sterile neutrino and the singlet scalar is in place: the heavier is S, the smaller its VEV and thus the lighter is N (see arrows in Fig. 6). If instead we require the level-repulsion feedback No sharp prediction see Fig. 6 Type-II ∆ = µ v 2 Table III. Dormant Higgs bosons masses, for different neutrino mass generation mechanisms, assuming (1) the level-repulsion feedback on the Higgs is limited to (100 GeV) 2 for each neutrino, and (2) µ = 100 GeV for all mixings (see text for details).
Here, v = 175 GeV and g = 0.7. The requirement of a minimum mass for the heaviest neutrino typically translates into a low bound on the mass of some scalar. to the SMH mass to be small, no sharp prediction can be made about M S , as it can always be made arbitrarily large leading to pseudo-Dirac neutrinos.
In the case of type-II seesaw, in the absence of N R , the ∆ would get its VEV from the usual term µ ∆ H T 0 iσ 2 ∆ † H 0 . Neutrino masses would be given by m type-II ν = 0.42 eV g 0.7 µ ∆ 100 GeV such that electroweak values of µ ∆ would correspond to scalar triplets around the 10 8 GeV scale, as can be seen in Table III, "Case (2)" column. Restricting the level repulsion feedback on the Higgs potential leads to exactly the same mass constraints as in the Dirac neutrino case. Obtaining scalar triplet masses accessible at the LHC, M ∆ ∼ TeV would require µ ∆ ∼ 10 eV, as it happens in usual type-II seesaw scenarios. For the three mechanisms discussed, a summary of the masses of the scalar bound states responsible for the neutrino spectrum can be found in Table III (labeled as "Dirac", "Type-I" and "Type-II").

B. Lepton flavor violation
Any charged lepton flavor violating (CLFV) decay would immediately indicate the presence of new physics and as a consequence much effort has been put into detecting such effects. The result of this work places severe constraints on the branching ratios for the different CLFV decays of the charged leptons.
In the absence of mixing among dormant Higgses, the CLFV decays must be mediated by the Higgses associated to the lepton mass generation. As a result the CLFV transitions induced in our framework do not involve quarks. The relevant constraints are the radiative and the three-body decays are shown in Table IV.
The NP contribution to the decay of a charged lepton to another through the emission of a photon is governed, in the low-energy theory, by the dipole operator where i and j are flavor indices and C ij has inverse mass dimension. The decay width is given by In the Scalar Democracy, NP effects are induced by diagrams where a dormant Higgs propagates in the loop. The coefficients C ij depend only on the rotation matrices of the external fermions and the Higgs mass.
For concreteness, we will assume that the Dirac neutrino mass mechanism is at work (see Table III). The NP contribution to the coefficients of the dipole operators, from both neutral and charged Higgses, is given by to leading order in m 2 e /(M e ) 2 , where M e kl denotes the mass of the scalar that generates the (kl) entry of the corresponding leptonic mass matrix. Due to the smallness of neutrino masses, the neutrino Higgs contribution is completely irrelevant to the radiative decay. As an aside, we note that a quick estimate of the NP contribution to µ → µγ will show that even in the most optimistic scenario C 22 is several orders of magnitude too small to significantly influence the SM prediction for the muon g − 2.
The charged lepton decay to three lighter charged leptons is mediated at tree-level by neutral Higgses giving the effective operators of Eq. (40). The resulting decay width is Many of these decays only violate flavor by one unit, and so, in contrast to the NP contribution to mesonantimeson mixing, the operator does not vanish for all dormant Higgses in the limit where R e = 1 3 . Similarly to the quark sector, the only direct constraint on the lepton rotation matrices stem from the requirement V PMNS = L ν L † e . Therefore, we have the same freedom of choice in their parametrization. In contrast to the quark sector, the lepton observables do not significantly favor small rotation angles. For our benchmark points we have therefore chosen large rotation matrices to illustrate the case where NP will soon become detectable in LFV observables.

Benchmarks
We apply the PDG convention for the PMNS matrix [55] using the best-fit values of the three mixing angles and phase, δ CP [56]. The use of the tribimaximal (TBM) mixing matrix [57][58][59][60] is motivated from work on discrete flavor model and its structure is given by For each of the benchmark points we have determined the strongest individual limit on each e-type dormant Higgs by comparing the NP contribution to the decay bounds in Table IV. We have plotted the leading constraints on the individual Higgs masses resulting from our benchmark points in Fig. 7 together with the corresponding mass estimates resulting from Eq. (32). In the benchmarks almost all the Higgs mass constraints come from the experimental bounds on the muon decays-even for the dormant Higgses coupling primarily to the third generation leptons. This is due to the somewhat anarchic structure of the PMNS mixing matrix, meaning that all the Higgses mediate LFV muon decays. The PMNS mixing angle suppression is insufficient to compensate the weaker bounds on the LFV tau branching ratios. Despite the slightly stronger experimental bound on Br(µ → eγ) than Br(µ − → e − e + e − ), the loop suppression in the leading dormant Higgs contribution to the radiative decay is enough that almost all the H e Higgses get their strongest bound from the decay to three charged leptons. For at least some of the Higgses it is possible to engineer the rotation matrices in such a way as to avoid the experimental bound on Br(µ − → e − e + e − ).
Regardless of the anarchic flavor structure chosen for the benchmark points, Fig. 7 shows that there is only severe tension between the mass estimates of H e , H µ , H eµ and H µe and their corresponding constraints in Benchmark 2. We therefore expect that most of the parameter space (consistent with our assumptions) will pass the lepton-flavor constraints. We also anticipate that future CLFV experiments, with significantly improved sensitivities, will be important probes of this scenario.
The Higgses associated with neutrinos in the Dirac mass scenario, H ν , are only constrained by the radiative decays. In no case, regardless of the left-handed mixing, are the bounds on the H ν masses stronger than 5 TeV, much below what is expected from the fine-tuning argument. For this reason the corresponding constraints have not been displayed in Fig. 7.

VI. CONCLUSIONS
While our principle of Scalar Democracy lacks a detailed theoretical underpinning, it is nonetheless experimentally testable in the coming years. If some predictions are confirmed, such as the observation of H b , and possibly H τ , H ct and H c , with expected O(1) Yukawa couplings, it could reshape our view of the UV completion of the Standard Model.
We propose that a plethora of new scalar bosons exists in nature. We argue schematically that every fermion pair binds to form a boundstate scalar boson, due to a universal attractive interaction at a very high scale, Λ. Amongst many new states, including lepto-quarks, colored isodoublets and singlets, etc., this hypothesis implies the existence of a large number of sequentially more massive Higgs bosons.
We argue that the SM Higgs boson is the first of the sequence, and therefore must have a dynamical origin, essentially as a tt boundstate, but it is now part of the constellation of composite scalars that influence its mass. One immediate intrigue is that the top quark mass, predicted by the renormalization group fixed point [18,19], which is the prediction of top condensation models [6], now comes within a few percent of the observed value. This is expected to be improved by including the multi-Higgs boson masses and decouplings and is a first indicator that this could be a pathway to new UV physics.
The universal binding we invoke may be intimately associated with gravity, since strong scattering amplitudes may exist near the Planck-scale that are parameterized by d = 6 operators and scale as 1/M 2 Planck . Amongst these amplitudes we would expected a general four-fermion structure such as (g 2 /M 2 )(ψ i L ψ j R )(ψ i R ψ j L ), where i, j run over the conventional quark and lepton flavors and colors of the SM. In fact, since gravity does not distinguish between particle and anti-particle, we would expect similar operators with a Majorana-like structure. The operators must be gauge invariant, and with the usual SM fields this contains a subsector of a global SU(48) × U(1) chiral Lagrangian.
We therefore argue that there exists a sequence of Higgs doublets leading upwards to large mass scales. These are "dormant" isodoublets, each having only small "tadpole" VEVs that arise from mixing with the observed Higgs and that scale as 1/M 2 for heavy Higgs fields of mass M . Our main point in invoking gravity is to motivate that many Higgs doublets may exist in nature. They could also arise from strong gauge interactions. The masses of these scalars are lacking a theory-we view them at present as soft-symmetry-breaking, relevant operators that we insert by hand in accord with known phenomenolgy.
The virtue of the model is that the Yukawa interactions are "known," and all are of order unity. This is interesting as it transforms the problem of flavor dynamics in a fundamental way: we have only a single Yukawa coupling at Λ, and the derived couplings at low energies are subject to relatively small renormalization corrections. The leptons may see the most significant renormalization to a value of about 0.7 (smaller 'fine-structure" effects will be generated by the decoupling of heavy Higgs fields). The observed spectrum of quarks and leptons and their mixing angles can be codified in terms of the the mass-mixing problem in the extended Higgs sector.
Most interesting is the top-bottom system. The top mass and the SM Higgs calibrate the overall Yukawa coupling for all quarks, g = g top 1.
Here the first sequential dormant Higgs H b , couples essentially to bb, and has a positive mass, M 2 b . It has a mixing with the lightest Higgs H 0 causes it to acquire a small VEV that feeds the b-quark its mass. We can view the SM Higgs as having initially positive M 2 0 (or zero), and the effect of "level repulsion" by H b drives it to become tachyonic. To avoid significant fine tuning we would require µ 4 /M 2 b ∼ M 2 H , which leads to a soft "naturalness" upper limit on M H b (v weak /m b )(∼ 100 GeV) ∼ 3.5 TeV.
We discuss in some detail a particular search strategy to address the neutral component of H b = (h + b , h 0 b ), through the process pp → h 0 b → bb. This is available at the current LHC energy and luminosity, and possible upgrades. We find that the LHC already has the capability of excluding sequential H b with O(1) Yukawa coupling to bb in a mass range of order < ∼ 1 TeV (as in the model of [4]). Our conservative (rough) estimates indicate that the full range of H b masses may be accessible to the high luminosity and/or energy-doubled LHC, and should help justify the case for such future machines. Our estimates are preliminary and likely will bound above what can be done by improved detector based studies and more detailed deployment of cuts and search strategies. We strongly encourage our experimental colleagues to pursue this. We thus defer details of possible searches h + b → tb and H τ elsewhere.
The challenge to Scalar Democracy is to avoid unwanted d > 4 flavor transition operators. The rare processes we focus on are mainly neutral meson mixing, and we find that such processes are binding on the Higgs spectroscopy. However, avoiding these constraints in this model is possible since most of the new sequential Higgs bosons are very heavy. Conceivably, though we have not explored it, the extended Higgs sector may provide sources of observable rare processes in heavy meson decays. Discovery of the lighter doublets could thus provide impetus for partially unraveling the flavor problem.
We also discuss the leptons. Our model most naturally leads to a "neutrino Higgs seesaw" for neutrino Dirac masses in analogy to the quark-sector Higgs structure. However, we have additional fermionic boundstates that can, in principle, mix to develop tadpole VEVs and produce Majorana masses, realizing a type-I or type-II seesaw mechanism. Hence the double-β decay experiments will be important probes of the far UV Higgsology of this framework. Our model can also drive rare lepton number violation experiments, and we focus discussion on a subset of these, i → j + γ and → 3 . We certainly do not exclude other possible processes as probes of the leptonic Higgs system.
One view of the future evolution of fundamental physics with energy scale, argues that the couplings are small and asymptotically free, and that the theory "fades away" into a linear scale-invariance. We have arrived at a contrary point of view, that a rich spectrum of new scalar states lies immediately beyond current energy scales and is within reach of the LHC and its upgrade path. In part this is motivated by chirality: it is very unlikely that one can generate a small Yukawa coupling constant from zero (certainly not perturbatively). Hence the tiny g electron ∼ 10 −6 is most likely due to a power law suppression of a coupling that is of order unity, such as g electron ∼ g µ 2 e /M 2 e where g ∼ 1. The power-law suppression then demands new mass scales, such as µ and M , as realized in our model. We have argued that a natural way to achieve this "democratically" throughout the entire flavor system of the SM is with a grand enlargement of the scalar system. This flips the flavor problem: the lightest (heaviest) fermions are coupled to the heaviest (lightest) scalars. The resulting spectrum of fermions is due to the inverted spectrum of associated scalars. This underlies our "Scalar Democracy." It focuses urgent attention to the top-bottom (and perhaps the τ -ν τ subsystem) and is testable at the LHC and future upgrades. We urge our experimental colleagues and other theorists to take up the cause.