Electroweak-Skyrmion as topological dark matter

We show the existence of a nontrivial topological configuration of the Higgs field in the Standard Model with the Skyrme term. It is shown that the current upper bound of the mass of the topological object is about 34 TeV. We discuss the impact of the existence of the topological object on cosmology.


Introduction
The nature of the dark matter remains a mystery of the Universe today. Numbers of candidates for the dark matter have been proposed based on models beyond the Standard Model (SM). Here, we show the existence of a dark matter candidate in the SM with a minimal addition of a higher dimensional operator. The idea is analogous to that of the Skyrme model [1], in which nucleons are described by nontrivial toplogical configurations of the Nambu-Goldstone (NG)-boson fields (pions) associated with the chiral symmetry breaking triggered by the quantum chromodynamics (QCD). Since the structure of the global symmetry breaking in the electroweak symmetry breaking (EWSB) sector of the SM is exactly the same as that of the chiral symmetry breaking in QCD, the argument by Skyrme in the QCD chiral Lagrangian is directly applied to the electroweak (EW) chiral Lagrangian by a simple scale-up of the relevant energy scale. However, there is one important difference between the low-energy physics of QCD and the EWSB sector of the SM: that is the existence of a light scalar particle, the Higgs boson. Actually, from the generalized analysis of the Skyrme model, it is known that the existence of a scalar resonance would generate the operator which destabilize the Skyrmion configuration when it is integrated out from the theory [2,3]. However, that does not immediately mean the non-existence of topological objects in the SM. In order to clarify the (non-)existence of such an object, it is necessary to study the system with the scalar field included as a dynamical degree of freedom. In this paper, we derive the equations which should be satisfied by the topological configuration in a coupled system of the NG and the Higgs fields, and show that a non-trivial topological configuration actually exists in the system. The mass of the topological object, which we call the EW-Skyrmion, is related to the magnitude of the coefficient of the O(p 4 ) term in the Lagrangian. Therefore the experimental bound on that coefficient places the upper bound on the mass of the EW-Skyrmion. We show JHEP07(2016)037 that the current experimental constraints give an upper bound of about 34 TeV, which is expected to be significantly lowered near future. The impact of the existence of the EW-Skyrmion on cosmology will be discussed, and we show that constraint from the direct detection dark matter experiments give a lower bound of O(1) TeV.
The paper is organized as follows. In section 2, we give a brief review of the original Skyrme model. Then, in section 3, we introduce a light scalar (higgs) into the system, and derive the coupled equations which should be satisfied by the topological object. It will be shown that a non-trivial topological configuration actually exists in the system. In section 4, upper bound on the mass of the EW-Skyrmion will be derived, and in section 5, cosmological impact of the existence of the EW-Skyrmion will be discussed.

Skyrme model
In order to provide a theoretical foundation and to introduce various notations, before going to the system which includes a scalar particle, we review the scaled-up version of the Skyrme model in this section. The Lagrangian is obtained by simply replacing the pion decay constant f π in the original Skyrme model Lagrangian by the EW scale v EW : and [A, B] ≡ AB − BA. π(x) is the NG field which parametrizes the coset space of SU(2) L ⊗ SU(2) R /SU(2) V . The second term of the Lagrangian (2.1) is often called the Skyrme term, and e (which should not be confused with the electric charge) is the numerical parameter which defines the magnitude of the coefficient of the Skyrme term. Skyrmion is described as a static field configuration which has the following "hedgehog" shape: F (r) can be determined by requiring that the total energy of the system is minimized. For the discussion below, it is convenient to introduce a dimensionless coordinater which is defined as the radial coordinate r normalized by the reference scale R: Then the energy of the static configuration can be written as Note that the above equation does not explicitly depend on either e or v EW . ThereforeF (r) takes the same shape for any values of e and v EW , and the total energy is proportional to v EW /e as one can see from the expression in eq. (2.7). Field configurations in this system are characterized by the topological winding number which is identified as the baryon number in the original Skyrme model. When the system has B = n (n: integer),F (r) as a solution of eq. (2.9) takes the following boundary values: The latter is required so that the system has a finite energy. The numerical solution of F (r) for the case of n = 1 is shown as the dashed curve in figure 1. The total energy, or Skyrmion mass M , can be obtained by substituting this solution into eq. (2.7) as (2.12)

Skyrmion in the EW chiral Lagrangian with a light scalar
In this section, we extend the scaled-up version of the Skyrme model described in the previous section so that a light scalar degree of freedom is included in the system. 1 We introduce a scalar field h(x) by linearizing the non-linear field U (x) in the kinetic term of the Skyrme-model Laglangian (2.1), and add the potential of h(x) in the following way: where the potential of h(x) is defined as follows: In the limit of e → ∞, the Lagrangian defined above is nothing but the SM after replacing derivatives with covariant ones. The newly added Skyrme term describes anomalous quartic JHEP07(2016)037 interactions among EW gauge bosons. Phenomenologically, we know only the upper bound of the coefficient 1/16e 2 from the experiments of the EW gauge boson scattering. As a possible origin, it has been known that the exchange of the vector resonance effectively generates the Skyrme term with the appropriate signature to support the Skyrmions. In general, h(x) and U (x) do not have to form a single linear field, thus the factor in front of the first term of the Lagrangian (3.1) could be any function of h(x)/f with f being a scale which is not necessarily related to v EW . Also a similar factor can be multiplied to the Skyrme term as well. There could be other types of higher order terms in the Lagrangian, and the form of V (h(x)) does not have to be limited to the one shown in eq. (3.2) as well. The choice of the Lagrangian above is to make the study tractable and to make the difference from the SM Lagrangian as minimal as possible: only the difference between the Lagrangian discussed in this section and that of the SM is the existence of the Skyrme term.
We take the same ansatz for the form of the static configuration of U (x) as shown in eq. (2.4). As for h(x), we assume that the static solution, h 0 (x), is spherically symmetric: With these ansatz, the total energy of the system takes the following form: Here,φ(r) is defined by φ(r) = φ(rR) ≡φ(r), and we replaced the parameter λ by the mass of the scalar m h by using the relation λ = m 2 h /(2v 2 EW ). The energy is minimized whenF (r) andφ(r) satisfy the following coupled equations: EWr 2 2φ(r) + 3φ(r) 2 +φ(r) 4 = 0. (3.6) We numerically solved the above equations, and found that the system always has topologically non-trivial field configurations. We call this topological object as the EW-Skyrmion.
In figure 1 The blue curve asymptotically approaches to that value in the limit of large m h /(ev EW ). This is an expected behavior considering the fact that taking the value of m h /(ev EW ) infinity corresponds to taking the decoupling limit of the scalar particle with e and v EW being fixed. Meanwhile, when m h /(ev EW ) is taken to be smaller, the normalized JHEP07(2016)037 mass of the EW-Skyrmion is reduced, and it approaches to about half of the value of the decoupling limit. In the next section, we use this behavior to place an upper bound on the mass of the EW-Skyrmion from currently available experimental data.
In the quantum theory of Skyrmion, the angular momentum of the configurations is quantized and the topological object is chosen to be either spin integer or half-integer, i.e., bosons or fermions. The ansatz of the shape of the solution as well as the coefficient of the Wess-Zumino-Witten (WZW) term, if exists, give relations among the topological charge, the angular momentum, and the gauge quantum number such as the electric charge [5]. Hereafter, we simply choose the spin-0 and neutral solution as the lowest energy state that corresponds to the absence of the WZW term. (See ref. [6,7] for general discussion of the quantum numbers of the Skyrmion in the context of little Higgs models.)

Experimental constraint
and experimental constraints on the coefficients α 4,5 can be obtained from the scattering data of EW gauge bosons [9,10]. The right hand side of eq. (4.1) can be rewritten as (4.2) where {A, B} ≡ AB + BA. The first term is nothing but the Skyrme term, while the second term is another O(p 4 ) term, which we assumed to be absent in the current study. Therefore, to obtain the constraint on the coefficient of the Skyrme term from 95% CL allowed region on the α 4 -α 5 plane, one should look at the one-dimentional projection of the allowed region onto the α 4 + α 5 = 0. By reading the constraint on α 5 on the line of α 5 = −α 4 , both ref. [9] and ref. [10] give a similar lower bound, α 5 −0.4, which can be translated to the upper bound on the coefficient of the Skyrme term as follows:

Cosmology
In this section, we discuss observational constraints on the EW-Skyrmion as a dark matter. We first discuss the possibility of the thermal production of the EW-Skyrmion, in which case, the relic density Ω S is related to the pair annihilation cross section σ A as follows: Here, we assume that the geometric cross section (∼ πR 2 , where R = 1/(e v EW ) is a typical radius of the EW-Skyrmion) is a reasonable estimate of the annihilation cross section, and substitute σ A v rel = π/(e v EW ) 2 to obtain a rough estimate of the right hand side of eq. (5.1). By requiring that the relic density of the EW-Skyrmion does not exceed the current dark matter density, Ω DM h 2 ≈ 0.1, we obtain the following constraint on the parameter e: e 150. (5.2) When e takes the maximum value, e 150, in which case the relic density of the dark matter is explained solely by the EW-Skyrmion, the mass of the EW-Skyrmion is about M = 60 GeV. The EW-Skyrmion with such a small mass is excluded by the direct detection experiment of the dark matter, such as the LUX experiment [12], unless the coupling relevant to the direct detection is extremely small, which is unlikely in the case of the current scenario. If the mass of the EW-Skyrmion is larger, the constraint from the direct detection experiments become weaker, however, in that case, the amount of the contribution of the EW-Skyrmion to the current dark matter density becomes smaller. Also, when the mass of the EW-Skyrmion is too large, say M = 10 TeV, the freeze-out temperature of the EW Skyrmion becomes much higher than the EW scale. Since the EW-Skyrmion

JHEP07(2016)037
exists only in the broken phase of the EW symmetry, it is not self-consistent to discuss the EW-Skyrmion beeing in the thermal equilibrium at the temperature much higher than the critical temperature of the EW phase transition.
From the above considerations, it is unlikely that the thermal relic component of the EW-Skyrmion becomes the dominant source of the dark matter of the current Universe. Since the EW-Skyrmion has the conserved quantum number (topological winding number), the asymmetry of the quantum number in the full description of the theory may have been produced in the history of the Universe [13] and that can remain today as the energy density carried by the EW-Skyrmion. The mechanism how such a production occurs depends on the UV physics which realizes the low-energy effective picture of the EW-Skyrmion as well as the history of the inflationally Universe. Instead of trying to specify those, we discuss the consistency of the asymmetry scenario with existing data and future experimental prospects.
In order to discuss the constraint from the direct detection dark matter experiments, we need to know how the EW-Skyrmion interacts with SM particles, especially with the Higgs. Calculating couplings requires the precise definition of the UV physics and the proper treatment of the extended object. Here, instead of trying to be precise, we adopt a simple assumption that the EW-Skyrmion has a low-energy effective coupling to the Higgs doublet with the form L eff = −2κ|S| 2 |H 2 |, where H is the Higgs doublet and S is the EW-Skyrmion field. (See, ref. [14] for analyses of the equivalent effective model.) The purpose here is to give a rough idea of the cosmological impact of the EW-Skyrmion. With this assumption, the spin-independent elastic EW-Skyrmion-nucleon cross section as a function of the EW-Skyrmion mass M can be expressed as follows: 4) where m N is the nucleon mass and f is the form factor, which is taken to be 0.3 here. The upper bound on the WIMP-nucleon cross section given by the LUX experiment [12] can be applied to the above cross section, and we obtain M 1.5 TeV when we assume κ = 1.0. Taking larger (smaller) value of κ places stronger (weaker) constraint: when the coupling is taken to be κ = π (0.5), the bound becomes M 3.5 TeV (1 TeV). Combining this lower bound with the upper bound obtained in the previous section, we conclude that the EW-Skyrmion with mass between 1.5 TeV and 34 TeV, with certain amount of uncertainty, is consistent with current experimental bounds, while explaining the dark matter of the Universe. As was mentioned in the previous section, LHC RUN2 will give stronger upper bound, M 4 TeV, which will probe the most of the allowed region obtained here. The sensitivity of the direct detection dark matter experiments will be also improved, therefore, there is a possibility that we will find a hint of the dark matter in near future.

Summary
In this paper, we formulated coupled equations for a topologically stable object in the system of the electroweak chiral Lagrangian with a light scalar boson. We applied those JHEP07(2016)037 to the Higgs sector of the Standard Model with minimal addition of the O(p 4 ) (Skyrme) term, and showed that a non-trivial solution exists in the system. From the experimental constraint on the magnitude of the coefficient of the Skyrme term, upper bound on the mass of the topological object (EW-Skyrmion) was derived as M 34 TeV. We discussed observational constraint on the EW-Skyrmion as a dark matter, and showed that nonthermally produced EW-Skyrmion which is heavier than O(1) TeV is consistent with the direct detection dark matter experiments. Both upper and lower bounds derived in this paper will be improved by LHC RUN2 and future direct detection dark matter experiments, and currently allowed mass region of the EW-Skyrmion will be entirely probed near future.