A generalized Higgs potential with two degenerate minima for a dark QCD matter scenario

We consider the Higgs potential in generalizations of the Standard Model. The possibility of the potential to develop two almost degenerate minima is explored. This would imply that QCD matter at two distinct sets of quark masses is relevant for astrophysics and cosmology. If in the exotic minimum the QCD matter ground state is electromagnetically neutral, dark matter may consist of QCD matter and antimatter in bubbles of the Higgs field. We predict an abundance of γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} rays in the few MeV region as messengers of dark matter regions in space. In addition the ratio of dark matter to normal matter is expected to show a time dependence.

The purpose of our Letter is to discuss a possible alternative of such scenarios based on exotic QCD matter. In a recent work the authors presented a detailed study suggesting that QCD matter depends crucially on the Higgs field [17][18][19]. Within the SM the quark masses in QCD are proportional to the Higgs field. As a consequence, changing its ground state value does change the quark masses in QCD, however, in a manner that keeps all quark-mass ratios fixed. In [19] a possible first order transition along the Higgs field trajectory was discussed. It is compatible with current QCD lattice simulations of the baryon ground state masses, but should be scrutinized by further dedicated QCD lattice studies.
In Fig. 1 we show our prediction of the baryon masses along the Higgs trajectory [17][18][19]. The bands in the plot provide an estimate of uncertainties based on our Fit 1 and Fit 2 a e-mail: m.lutz@gsi.de (corresponding author) scenarios as discussed in [19]. At fixed ratio m s /m = 26 the masses are plotted as functions of the strange quark mass.
The key observation is that within a critical region of the Higgs field, baryonic matter and antimatter are composed from Λ andΛ particles rather than from nucleons and antinucleons. This follows from the relation M Λ < M N , which holds at a specific range of the strange quark mass. We point out that our dark-matter scenario does not rely necessarily on a first order transition. Since Λ particles are electromagnetically neutral such matter does not radiate and therefore appears dark. Since the Higgs sector of the SM drives a possible electroweak phase transition and underlies baryogenesis models (see e.g. [20][21][22][23]) it is important to explore exotic Higgs sector generalizations of the SM in more detail. Our dark-matter scenario should not be confused with the quarknugget scenario proposed by Witten a long time ago [24], in which dark matter would consist of deeply bound strange objects in the non-exotic Higgs phase.

The Higgs potential
We consider the Higgs sector of the SM [20,23,[25][26][27][28][29][30][31][32][33][34]. At tree-level the Higgs potential in the SM can be expressed in terms of two parameters only with the complex doublet Higgs field H , the Higgs mass parameter M h 125.2 GeV and the vacuum expectation value v 246.2 GeV of the Higgs field in its physical vacuum state [35]. The value of v plays a decisive role in the QCD part of the SM since all quark masses are proportional to v. In this work we are interested in the Higgs potential at H † H ≤ v 2 , where it is known that even loop corrections in the SM are sizeable (see e.g. [25,27,[29][30][31]36,37]). Since the Higgs potential will be affected in most extensions of the SM we follow here a phenomenological path where we explore the consequence of a fine-tuned potential with two degenerate minima. An effective field theory approach that implies two degenerate minima would require at least (H † H ) 3 and (H † H ) 4 operators. Consider the specific form with v a the vacuum expectation value of the Higgs field at the exotic minimum. By construction, the model potential (2) has two degenerate minima. At its physical one it recovers the empirical mass of the Higgs. The ratio r = v a /v determines the strange quark mass in the exotic minimum. The puzzle with (2) is that it may be unnatural in the size of its dimension-full operators. However, we may recast the problem by considering loop corrections (see e.g [25,27,30]). In the presence of multi-loop effects we may use the phenomenological ansatz where the particular form of the log term with the parameter γ is taken from [31]. There the value γ = 0.1 is used. We note that the Higgs sector is the least controlled part of the SM and therefore may be subject to significant model modifications. According to Fit 1 and Fit 2 we expect dark QCD matter in the range 0.39 < r < 0.57 and 0.39 < r < 0.54 respectively. The critical values are close to those as derived in [19] on the unphysical trajectory where m u + m d is kept constant. In Fig. 2 we plot the effective potentials of (1-3) as a function of   (2) and (3). It is compared with the tree-level potential of the SM (1) eter γ we can efficiently dial the magnitude of the Higgs potential close to the origin. In the vicinity of the two local minima we find a rather mild dependence on the form of our parametrization. The polynomial ansatz (2) or the log form (3) lead to almost indistinguishable results. We emphasize that both models are compatible with empirical constraints on the Higgs potential as discussed in [38]. For instance at r = 0.45 we extract from (2) and (3) the range 3.8 < κ λ < 6 for the three Higgs coupling constant κ λ . This is well compatible with the empirical 2-σ interval −5.0 < κ λ < 12.1 from ATLAS [38]. Our estimate excludes the SM value κ λ = 1.
One may object to such a fine-tuned Higgs potential. However, we wish to recall that there are ample cases in physics in which a system is driven by fine-tuned dynamical assumptions. In particular the SM itself has various fine-tuning issues already. At this stage of the development we would not worry too much. Rather, we discuss in some detail the consequences of a possible dark QCD matter scenario.
Dark matter is believed to account for approximately 85% of the matter in the universe (see e.g. [39]). How would we arrive at such a ratio in our scenario? We may assume here that the universe reheats above the electroweak scale after inflation. Alternative scenarios, where for instance the energy scale of inflation is around the electroweak scale (see e.g. [40]), are not considered here. Then the universe is baryon matter dominated already when the electrocweak symmetry is restored with H = (0, 0) (see e.g. [31,[41][42][43][44]). As the temperature lowers further our effective Higgs potential should take over. Here it is important to estimate the size of thermal effects. In the limiting case of sufficiently large T they are proportional to T 2 V (H ) (see [45]). Contrary to the SM, the finite temperature effect on the generalized Higgs potential is not monotonic in the Higgs field. From its functional form as modeled in Fig. 2 we expect the two local minima to merge to one global minimum. This is so since our two minima imply multiple sign changes in V (H ). Thus, at this stage of the cosmic evolution, matter sits in a 'conventional' Higgs field condensate. Only as we further lower the temperature, the additional exotic minimum will turn visible. We do not see any strong hint that in the electroweak era of the cosmic evolution any seed of dark matter is formed.
At the QCD scale T ∼ 1 GeV thermal effects from quarks and gluons will dominate the temperature effects in the Higgs potential. Here the competition of the two possible Higgs phases is more intricate. Despite our assumption on the almost degeneracy of the two Higgs phases at zero temperature, there will be a significant asymmetry from the dynamics of quarks and gluons. It is safe to assume a vanishing baryon chemical potential and consider the difference of the pressure densities in the two available local minima. The exotic Higgs phase wins here, since it is characterized by quark masses that are about twice as small as those in the normal phase. This is illustrated in Fig. 3, where the difference in the pressure densities of the two Higgs phases is shown as a function of the temperature T . At T = 1 GeV that difference of about Δp QC D 28 GeV/fm 3 is dominated by the presence of the heavy quark. If we ignored the charm quark the corresponding value would be Δp QC D 0.28 GeV/fm 3 . The computation follows from a non-interaction Fermi gas ansatz for the quark-gluon-plasma phase of QCD, which may serve as a rough estimate. Note that, the contributions from gluon degrees of freedom cancel out in the considered difference, at least approximatly [46,47]. Thus the results are determined by the quark masses in the two phases and the temperature only. In the limit of large temperatures we find, Δp Eventually the formation of such exotic bubbles will stop as the temperature is further reduced. We expect this to happen after QCD changed its degrees of freedom from quarks and gluons to hadrons at temperatures T < 150 MeV. Thus, altogether we suggest that dark baryonic matter is formed via a first-order phase transition taking place in the quark-gluonplasma era of the cosmic evolution. To reach the target of 15% ordinary to dark matter ratio appears realistic, however, requires detailed knowledge of QCD dynamics during that phase transition. A significant estimate of that ratio requires a quantitative study of the evolution process. That is beyond the scope of the current work.
If the exotic minimum in the Higgs potential is slightly metastable at zero temperature, we expect a scenario where the vacuum shows bubbles with dark QCD matter inside, but normal QCD matter outside. Inside the bubbles the matter or antimatter ground states consist of Λ orΛ particles, however with exotic properties as shown by Fig. 1.
Let us explore the stability of a possible Higgs bubble. Since the boundary of such a bubble stores a significant amount of energy, there is a tendency that such a bubble shrinks or even collapses. From Fig. 2 we estimate the energy density Higgs 1.1 × 10 9 GeV/fm 3 (4) from the Higgs potential taking in-between its two minima. We now assume a Higgs bubble with spherical geometry characterized by a radius R and a surface thickness d. That implies the total surface energy where the second term in (5) follows from the kinetic term of the Higgs field. From Fig. 2 we can read off the change of the Higgs field across the surface with (Δv) 2 v 2 /8. We estimate the bulk energy by a free-Fermi gas approximation with the gravitational constant G 6.709 × 10 −39 GeV −2 and N the total number of Λ's in the Higgs bubble. Their Fermi momentum is denoted by k F with ρ = k 3 F /(3 π 2 ), where ρ specifies the dark-matter density in the bubble. We parameterize a supposedly small difference in the vacuum energy densities at the two Higgs minima by Δ > 0, where we assume the dark-matter vacuum to be slightly disfavored.
A Higgs bubble can be stable, provided that it encloses a sufficient amount of dark matter. We can make this more quantitative by a minimization of its energy E = E bulk + Fig. 4 The radius and Fermi momentum as a function of N at fixed Δ . Note that N Sun M Sun /M Nucleon 1.2 × 10 57 E surface with respect to the surface thickness d and the radius R at a fixed value of the total number of Λ's in the bubble. From this we find the two relations, This implies that at given k F it follows that Δ must be smaller than a critical value, as to keep the dark-matter bubble stable. We checked that all second derivatives are positive so that with (7) we have at least a local minimum of the dark-matter system (5, 6). In Fig. 4 we show the radius, R, and Fermi momentum, k F , of the Higgs bubble as a function of N at various fixed values of Δ . Within the range 10 31 < N < 10 57 the value of k F < 250 MeV is small enough to justify our free-Fermi gas approximation. We expect our results to hold at the qualitative level.
It is left to check whether such a dark matter Higgs bubble is stable with respect to a decay into a more conventional object consisting out of normal baryonic matter. A useful quantity to consider is the energy per particle in the bubble, (E surface + E bulk )/N , with d and R as given in (7). In Fig. 5 we show such a dependence at various fixed values of Δ . The energy per particle is significantly smaller than the free nucleon mass and therefore, at least in the region 10 31 < N < 10 57 , where the effects from gravity are not dominating the system yet, there is no phase space available for such a decay. As an example consider N = M sun /M N 1.2 × 10 57 for which it follows R 16 km and k F 250 MeV at Δ = 0. The particle density with ρ 0.07 fm −3 would be sufficiently small as to justify the application of the free-Fermi gas approximation. As we further increase N > N Sun the dark matter will turn more dense and will be sensitive to the equation of state of the Λ particles in the Higgs bubble. It would be important to establish the latter from QCD and to derive mass limits for the dark matter bubbles. Here the Tolman-Oppenheimer-Volkoff equation has to be supplemented by a suitable boundary condition at the Higgs bubble surface. This may open the possibility for the existence of massive compact objects, with properties distinct to those of neutron stars and/or black holes of unconventionally small masses [48]. In a meson-exchange phenomenology, the repulsive omega-meson-exchange process is expected to dominate the short-range interaction of Λ particles in the Higgs bubble. Due to the approximate isospin conservation, pion-and rhomeson exchange processes are suppressed. It remains the etameson exchange and the two-pion exchange contributions, which may bring in some weak intermediate-range attractive forces [49][50][51].
While there appears to be a rather weak net attraction at the physical point [51][52][53][54] available studies suggest a sizeable quark-mass dependence thereof [49,[53][54][55]. We conclude that at the exotic Higgs minimum, that comes at much smaller up and down quark masses, there is little evidence to expect this weak attraction to survive. Whether and how massive dark-matter clusters form depends on the subtle balance of the gravitational force and the short-range strong interactions in the Higgs bubble.
We conclude that in any case the typical dark-matter density in a Higgs bubble should be significantly larger than the density of a cold interstellar medium, which is characterized by a baryon-number density smaller than about 10 6 /cm 3 . In this context we discuss the so-called Bullet Cluster [56,57]. While the radial velocity distributions of stars inside a galaxy or data on gravitational lensing effects (see e.g. [8,39,[58][59][60]) put constraints on the dark-matter distributions in and outside galaxies, more significant information on the possible nature of dark matter is set by the observation of collisions of galaxy clusters [56,57]. It is found that in such a collision there is no direct hint pointing at any sizeable interaction of dark matter with ordinary matter [57]. In this context we have to discuss how a Higgs bubble interacts with protons from the intergalactic hot gas. The relative velocity of the two colliding galaxies in [56] is of the order of 4500 km/s. An intergalactic gas of temperature T 6 keV implies a typical proton velocity of about 1300 km/s. Thus most of the protons from the gas do not have sufficient kinetic energy to invade the bubble. In turn there will be no strong interaction effects visible. Second we need to consider the case where Higgs bubbles from the two galaxies collide. The chance that this happens depends on the typical size of such bubbles, which is not well constrained at this stage. It depends on the details of the Higgs potential, in particular the size of Δ term, and a cosmological model. The smaller the typical size of the Higgs bubbles, the smaller the likelihood that such a process turns relevant in a galaxy merger event. Even if two bubbles start to overlap, we would expect that the two bubbles merge into a larger one, since this reduces the energy stored in their surface. The residual interaction of the Lambda particles with kinetic energies of at most a few MeV should be dominated by elastic processes. In turn we do not see any reason to expect a strong visible effect of the dark-matter component in such a galaxy collision event.
Last, we turn to a most interesting process where a sufficiently energetic cosmic proton tries to enter a dark-matter region in space with relative velocity, v p . Note that, depending on the energy such a proton may be even trapped inside the dark-matter bubble and therefore the ratio of dark matter to normal matter is expected to show a time dependence in our dark-matter scenario [61]. According to Fig. 1 the nucleon mass inside the bubble is only up 10 MeV = ΔM N larger than its mass outside the bubble. Thus on the way into the bubble the nucleon has to either transfer momentum to the Higgs bubble and/or radiate photons. Such a Bremsstrahlung spectrum should be limited to γ rays with energies less than that 10 MeV. Here a crucial parameter is the acceleration, a (c 2 /γ 2 p ) ΔM N /(d M N ) 8 × 10 23 (c/γ 2 p )/s, of the proton across the Higgs bubble surface, since its total radiation power is proportional to a 2 γ 4 p with γ p = (1 − v 2 p /c 2 ) −1/2 . To this extent our Higgs bubbles glim with a characteristic spectrum which depends on the details of the Higgs potential.
So far we made a rough estimate for the flux of sub MeV photons. To be explicit, we assumed a dark matter bubble of radius 16 km moving with 4500 km/s through an intergalactic gas of 6 keV temperature inside the bullet cluster. For this, the integrated X-ray flux arriving at a detector on the earth is 10 −48 W/cm 2 . It is instructive to confront this value with the integrated flux of 5.6 × 10 −19 W/cm 2 for X-rays with (0.1-2.4) keV from the Bullet Cluster [56,57]. Our corresponding estimate for integrated X-ray flux from dark-matter bubbles is less than 10 −36 W/cm 2 , i.e. down by 17 orders of magnitude. A similar value we predict from the photon emission of a corresponding single dark-matter bubble in the Milky Way, for which we estimate the rate to be smaller than 10 −38 W/cm 2 .
It is not very likely that such dark-matter photons can be detected with satellite-based detectors like e-Astrogram or AMEGO [62][63][64]. In particular, we note that so far there is gamma ray data in the (1-10) MeV region available with quite large uncertainties only.

Summary and conclusions
We constructed a phenomenological Higgs potential with two degenerate local minima. It was argued that such a generalization of the SM may lead to dark QCD matter that lives in bubbles of the Higgs field, with normal QCD matter outside and dark QCD matter inside. Within the bubbles we expect exotic Λ andΛ particles, that are formed by QCD at unconventionally small up, down and strange quark masses. We predict an abundance of γ rays in the few MeV region as messengers of dark-matter regions in space. In addition the ratio of dark matter to normal matter is expected to show a time dependence.
It would be interesting to further scrutinize the dark QCD matter scenario proposed here. With current QCD lattice techniques it is possible to substantiate or rule out such a scenario by further studies of the strange quark-mass dependence of the nucleon and Λ baryon masses. It would be important to establish a more fundamental framework in which such an exotic Higgs potential is implied.

Data Availability Statement
This manuscript has associated data in a data repository. [Authos' comment: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .