Very heavy dark Skyrmions

A dark sector with a solitonic component provides a means to circumvent the problem of generically low annihilation cross sections of very heavy dark matter particles. At the same time, enhanced annihilation cross sections are necessary for indirect detection of very heavy dark matter components beyond 100 TeV. Non-thermally produced dark matter in this mass range could therefore contribute to the cosmic \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}γ-ray and neutrino flux above 100 TeV, and massive Skyrmions provide an interesting framework for the discussion of these scenarios. Therefore a Higgs portal and a neutrino portal for very heavy Skyrmion dark matter are discussed. The Higgs portal model demonstrates a dark mediator bottleneck, where limitations on particle annihilation cross sections will prevent a signal from the potentially large soliton annihilation cross sections. This problem can be avoided in models where the dark mediator decays. This is illustrated by the neutrino portal for Skyrmion dark matter.


Introduction
The fact that direct search experiments so far could not confirm a dark matter signal in the theoretically well motivated WIMP mass range between about 10 GeV and a few TeV creates increasing pressure to look for light dark matter particles or for very heavy dark matter as alternative explanations of the dark matter puzzle. Light dark matter models can be motivated through axions, dilatons, or moduli fields in string theory, and they will be tested by upcoming experiments.
On the other hand, superheavy dark matter with masses above 1 EeV = 10 9 GeV had been discussed extensively as a consequence of initial lack of observation of a Greisen-Zatsepin-Kuzmin (GZK) cutoff [1,2] in the ultra-highenergy cosmic ray spectrum [3][4][5][6]. The early pioneering paper on ultra-high-energy cosmic rays from superheavy dark matter was Hill's paper on monopolonium decay [7]. a e-mail: rainer.dick@usask.ca However, the discovery of a GZK cutoff in the meantime [8][9][10][11], the successful matching of the spectrum above 3 EeV in terms of nuclear components [12], and the increasingly stringent limits on the fluxes of ultra-high-energy neutrinos [13] and photons [14,15], indicate that superheavy dark matter, if it exists, will not be detectable in direct or indirect search experiments as we know them. This is not entirely unexpected, since we knew from the start that the unitarity limits on the annihilation cross section of superheavy dark matter make it an ideal candidate for practically secluded dark matter in terms of particle physics experiments [6].
However, there is a mass range between the WIMP mass range and the superheavy mass range that warrants further exploration. We are particularly interested in the very heavy dark matter mass range between about 100 TeV and 10 PeV = 10 7 GeV because of the possibility that there might be a detectable flux of very high-energy γ -rays between 100 TeV and a few PeV, and because IceCube has seen neutrinos with PeV scale energies [16,17]. A detectable contribution from dark matter annihilation in this energy range calls for solitonic enhancement of annihilation cross sections, because the indirect search limits in the TeV mass range are already encroaching on the thermal production limit [18]. This makes indirect dark matter signals from dark particle annihilation for higher masses very unlikely, as these signals can be expected to drop with dark matter mass M roughly in proportion to M −3 , i.e. faster than the decrease of the cosmic γ -ray flux with energy. The indirect search limits in the TeV mass range therefore pose the question whether there is any hope for a potentially detectable indirect dark matter γ -ray signal in the energy range beyond 100 TeV, which could be motivated by theory. Solitonic states can avoid this negative verdict on indirect signals from very heavy dark matter, because their annihilation cross sections are size limited σ ∝ L 2 rather than mass limited σ ∝ M −2 . Non-topological dark solitons could arise as e.g. as dark Qballs [19][20][21][22][23][24][25][26][27]. However, in the present paper we focus on Skyrmions as an example of topological dark solitons.
Sommerfeld enhancement provides another way to achieve high cross section values and potentially observable indirect signals for heavy dark matter if the dark matter particles participate in interactions which are long range compared to their Compton wavelength M −1 [28]. This applies especially to heavy dark matter which is weakly charged [29][30][31][32][33][34][35][36][37], or to heavy Majorana dark matter which can exchange scalar bosons of lower mass. However, we will see that the size-induced solitonic enhancement factors can reach levels of order 10 7 for Skyrmion couplings of order g V 0.1, and therefore the size effect alone will be sufficient to generate observable signals from very heavy dark Skyrmions. Therefore we develop very heavy dark Skyrmion models where the Skyrmions are not charged under long range gauge interactions, although charged very heavy Skyrmions are an interesting topic for further investigations.
The question for indirect signals from very heavy dark matter is a timely question to address, since IceCube is already exploring neutrinos in the PeV energy range while the Cherenkov telescope array (CTA) will start to explore cosmic γ -rays in the energy range beyond 100 TeV. Searches for indirect dark matter signals will likely remain the primary, if not the only option to explore the mass range beyond 100 TeV in the foreseeable future [38]; see also [39] for an excellent discussion of motivation and mass reach of next generation colliders.
The lower end of the heavy mass range poses interesting theoretical challenges which force us to deviate from standard dark matter theory when we cross this boundary. Indeed, it is well known that standard thermal creation of dark matter particles would imply overclosure of the universe if the dark matter particles are heavier than about 100 TeV [40]. This problem could eventually be ameliorated by dilution of dark matter through late time entropy production [41,42]. However, there is also another, less well known problem appearing at about 10-50 TeV, depending on the dark matter model: Standard thermal dark matter creation also requires increasingly strong dark matter coupling g D−SM ≡ g of the dark matter to Standard Model particles to produce heavier particles. Once we reach a mass of a few ten TeV, the required coupling to standard model particles needs to be of nonperturbative strength g D−SM ∼ 4π for thermal dark matter creation. The reason for this is that the thermally averaged annihilation cross section at temperature T [43], varies asymptotically like M −2 for large dark matter mass, whereas the required cross section σ v r (T f ) for thermal dark matter creation at the freeze-out temperature T f only varies logarithmically with dark matter mass. Every perturbative dark matter model leads to a prediction for σ (s) = g 2 × σ (s)| g=1 , and thus to a prediction for σ v G (M, is determined by cosmology. This leads to the prediction of the masscoupling relations for dark matter models and the quadratic decrease in σ v G (M, T f )| g=1 with mass leads to a linear increase in the required dark matter coupling g(M) for thermal dark matter creation of very heavy dark matter. Due to the perturbative approximation σ (s) = g 2 × σ (s)| g=1 , this reasoning cannot be extended for calculating the required dark matter couplings in the non-perturbative high-mass regime, but it tells us for each model a mass bound where perturbative dark matter theory breaks down. Contrary to the breakdown in perturbativity of dark matter models for high dark matter mass, the overclosure problem does not involve any perturbative calculation and only assumes that particle creation and annihilation is described by a scattering matrix [40]. However, both the calculation of g(M) and the overclosure problem assume thermal creation of dark matter. Therefore both the overclosure problem and the need for non-perturbative coupling can be avoided through non-thermal creation of very heavy dark matter. Nonthermal dark matter creation can be achieved in several scenarios during or immediately after inflation [44,45]. Gravitational production due to the rapidly evolving scale factor during inflation or near its end is a promising possibility [46][47][48][49][50][51]. Other mechanisms for very heavy dark matter production include the preheating [52,53] and the reheating phases [46,54,55] at the end of inflation, or resonant production due to an evolving effective mass term from couplings to the inflaton through Yukawa type couplings [56,57] or kinetic couplings [58]. These are different proposals using different physical models for non-thermal dark matter creation in the early universe. However, there is one common denominator that is worth emphasizing: None of these proposals needs to describe the creation dynamics through a scattering matrix and corresponding reaction cross sections, and the standard thermal freeze-out estimate X ∝ σ v −1 [59] for the remnant dark matter abundance does not apply. Indeed, all the possible mechanisms for very heavy dark matter production are inherently semi-classical, either through directly integrating coupled systems of Lagrangian evolution equations or by evolving Bogolubov coefficients in a rapidly evolving classical spacetime.
This apparent semi-classical aspect of very heavy dark matter generation matches nicely with the observation that a detectable indirect signal in the cosmic γ -ray and neutrino flux beyond 100 TeV would favor a solitonic component in very heavy dark matter. To elucidate this point, we note that the spectral flux of cosmic rays from annihilation of heavy particle dark matter of density n X (r) and dark anti-matter of density 1 n X (r), can be estimated to scale with dark matter mass M roughly like M −3 , since the densities scale like M −1 and the velocity weighted annihilation cross section for very heavy dark matter (i.e. M much larger than the top quark mass) scales according to v X X σ X X ∝ M −2 , while more fragmentation products could be expected in proportion to M at lower energies. Equation (15) below provides an explicit example for the asymptotic behavior of the velocity weighted annihilation cross section. Comparing the expected drop in cosmic ray flux from annihilation of dark particles with the fact that the combined spectral flux of cosmic rays in all particles scales like E −2.7 for E 3 PeV, tells us that indirect signals from very heavy particle dark matter would be buried deeply in the cosmic ray flux from astrophysical accelerators. On the other hand, the annihilation signal from solitonic dark matter of size L S and mass M S would scale like L 2 S M −1 S , and could therefore contribute at a detectable level to the cosmic ray flux above 100 TeV. For Skyrmion dark matter the enhancement of the annihilation cross section is of order This yields an enhancement of order 10 6 for a weak scale Skyrmion coupling g V 0.3. This is relevant for indirect dark matter searches beyond 100 TeV. A Higgs portal Skyrmion model is therefore introduced in Sect. 2.
Skyrmions can arise as a consequence of a first order phase transition, e.g. due to chiral symmetry breaking, and Campbell et al. have shown that this can create the correct Skyrmion abundance for dark matter [60]. In these cases, the Skyrmions are indeed the dominant form of energy at least at early stages after the phase transition. However, we will see that a Higgs portal coupling implies that the w bosons, into which the Skyrmions annihilate, will generically also contribute a sizable particle dark matter component if the w bosons are stable. This reduces the observational significance of the Higgs portal model for Skyrmions. To avoid this problem, we also discuss a model where the bosons w decay into νν and νν pairs. We will address this model, which was analyzed for WIMP scale dark matter coupling to very heavy right-handed neutrinos by Dudas et al. [61], as the ν 2 -portal for dark mat- 1 The flux is usually written with the thermally averaged velocity weighted annihilation cross section σ v (T ), see also Eq. (1), but heavy dark matter at low redshift is practically at zero temperature. Furthermore, a generic expression dN (E, 2M)/dE has been chosen for the fragmentation of an annihilation event with energy √ s = 2M instead of the customary 2dN (E, M)/dE, which assumes fragmentation of two initial jets of energy M. The reason is that annihilation of bound states or solitons might not necessarily produce only two jets. ter, to avoid confusion with the neutrino(-Higgs) portals proposed in [62,63] (see also [64]). In the neutrino(-Higgs) portal models, unstable dark fermions χ i couple to left-handed fermion doublets i in the Standard Model through the same couplings as the right-handed neutrinos, λ i j i ·H · χ j +h.c., whereH = · H * = (H 0 * , −H + * ) is the Higgs doublet in the complex conjugate fundamental representation of the electroweak gauge group SU(2), mapped back into the fundamental representation. These models and their generalizations to higher-dimensional operators attracted a lot of interest in recent years due to their possible relevance for PeV scale neutrinos; see [65][66][67][68][69][70] and the references therein (see e.g. [71][72][73] for discussion of the possible astrophysical sources of PeV neutrinos). Here we wish instead to couple scalar bosons w i to the Standard Model through λ i jk ν i,R c ·ν j,R w k +h.c., utilizing the fact that the symmetries of the standard model are compatible with Majorana terms for the right-handed neutrinos. The w particles in this model are not the dark matter but are generated from dark Skyrmion annihilation, and they will decay very fast into ν R ν R pairs or ν R ν R pairs, with the right-handed neutrinos then decaying into left-handed neutrinos and Higgs particles. This model has the virtue that Skyrmion annihilation into the w bosons cannot build up a competing dark matter component.
The Higgs portal Skyrmion model is introduced in Sect. 2 and the ν 2 neutrino portal is introduced in Sect. 3. Section 4 summarizes our conclusions.

A Higgs portal model for heavy Skyrmion dark matter
Skyrmions have always been part of the toolsets of nuclear theory and mathematical physics, and they have also become very popular in condensed matter physics. Their emergence in black hole physics [74,75] and as dark matter candidates [76][77][78] is a relatively recent development that should not surprise us, given that we must naturally encounter low-energy effective field theories in every field of physics. Motivated by techniques of low-energy hadron physics and the observation that we need enhanced annihilation cross sections to see indirect signals of very massive dark matter, we model the dark matter components as Skyrmion excitations of a w-field where the σ i are the Pauli matrices of a dark SU(2) isospin symmetry, w i (x) is a triplet of states with the lowest mass in the very massive dark sector above 10 TeV, andŵ = w/|w| is the corresponding unit vector. We couple the resulting dark Skyrmion model to the standard model through a Higgs portal coupling. This will give rise to a two-stage annihilation process of very massive dark Skyrmions into Standard Model states. Skyrme had proposed [79,80], and Witten et al. [81,82] had demonstrated, in the large N limit of SU(N ) gauge theory, that baryonic states can be realized as topological excitations of mesonic states, and this reasoning would also apply to a dark gauge theory sector. We will show that this observation is particularly relevant for very heavy dark matter, since it provides a means to enhance the reaction cross sections for very heavy dark matter to observable levels. However, we remain agnostic with respect to the question whether heavy Skyrmion dark matter indeed arises as an effective description of very heavy bound states in a dark gauge theory sector, or as a genuine solitonic excitation of a heavy scalar field. The point is that either way, the resulting enhancement of annihilation cross sections makes this an interesting target for indirect dark matter searches with masses exceeding 100 TeV.
The proposed Skyrme model for heavy dark matter states has a Lagrangian with a Skyrme part where the notation g V for the Skyrme term coupling reminds us that this coupling may arise from vector dominance due to a hidden gauge invariance [83][84][85][86] (see also [87]). The term L D−SM describes the coupling to the standard model. A natural guess is to implement Higgs mediation into standard model states. Higgs exchange had been suggested on several occasions as a minimal coupling mechanism between dark matter and the Standard Model [88][89][90][91][92], and a Higgs portal coupling for a very heavy dark Skyrmion model would take the form Here H is the electroweak Higgs doublet with vacuum expectation value H + H = v 2 h /2. The coupling (6) would not contribute to any invisible Higgs decay width since all the masses of the dark sector states, including the mass m w of the w i (x) fields, are assumed to be much larger than the Higgs mass. It is also not constrained by direct search experiments which are not sensitive to the mass range above 10 TeV. The w-Higgs coupling in unitary gauge is This model differs from little Higgs Skyrmion dark matter models [76,77] by not assuming that the Higgs field itself is part of the fields w i , which parametrize a coset in the little Higgs models. The model is also different and much less constrained than the model in Ref. [78], where the Higgs field was coupled to the kinetic term Tr(∂ μ U · ∂ μ U + ) and the non-linear fields U (x) were related to the electroweak gauge bosons. These models therefore also did not include the Higgs portal coupling (7).
The scale f w is a mass scale, but contrary to the hadronic Skyrme models, it is not a decay constant. Recall that the π mesons of hadronic physics can only decay because their constituents couple to the lighter lepton sector through the electroweak gauge bosons. There is no corresponding low mass dark matter sector included in the w+Skyrmion picture of very heavy dark matter, and the w-particles in the model (4,5) with the coupling (7) are stable up to annihilation into Standard Model states through the Higgs portal. Therefore the dark sector in this type of dark Skyrmion model will generically consist both of w particles and of their Skyrmion excitations. Skyrmion annihilation and thermal creation can generate w-particles which will also contribute to the dark matter if they do not annihilate sufficiently fast into Standard Model states. Nevertheless, the Skyrmions can initially be produced as the dominant dark matter component if they arise as a consequence of symmetry breaking during a first order phase transition [60], and this is the assumption used here. However, it is also intriguing to ask what happens if the w particles are not just a low-energy effective description of a condensate in the low-energy effective theory of a broken symmetry, but are thermally produced in the early universe before Skyrmions are generated in a phase transition. For the reasons alluded to above, we can perform a perturbative analysis of this question only if the mass m w of the w particles is not too large, and we will return to this question once we have assembled the pertinent cross sections.
The Hamiltonian for the Lagrange density (4) is where the summation over lower pairs of 4-indices is Euclidean, e.g. ∂ μ U · ∂ μ U + ≡ ∂ 0 U · ∂ 0 U + + ∇U · ∇U + . Skyrme had demonstrated that the Lagrangian L S generates solitons which are stable due to the topologically conserved current (with the convention 0123 = −1) which he had proposed to use as a baryon current. Of course, in our model W μ appears as a Skyrmion current in the very heavy dark sector. The winding number W = d 3 xW 0 measures how often the mapping U (x) in Eq. (3) (with the compactifying boundary condition lim |x|→∞ U (x) = 1) wraps compactified R 3 around the SU (2) group manifold S 3 . The mass and the length scale of |W | = 1 Skyrmions are M S = 73 f w /g V and L S = 1.2/(g V f w ) [87,93]. The solitonic nature of Skyrmions implies that they are not described as particle states in the Fock space of the theory, and Skyrmion-Skyrmion interactions are not described by the usual scattering matrix formalism. Instead, their interactions have to be studied through mathematical analysis and numerical integration of the equations of motion, and study of the evolution of the topological density W 0 [94][95][96][97][98]. For the annihilation of a Skyrmion S and an anti-Skyrmions S, this leads in particular to the interesting result of sudden onset of annihilation through emission of a few w quanta once the distance between the Skyrmion and the anti-Skyrmion is down to the size of a Skyrmion [94,95]. It is a classical soliton-soliton interaction effect and determined by soliton size, whence the -wave unitarity limit σ ≤ 4π(2 + 1)/k 2 on reaction cross sections from scattering matrices [40] does not apply. On the other hand, the underlying w particles (into which the Skyrmions decay upon annihilation in our adoption of the Skyrmion picture for very heavy dark matter), are the basic quantum excitations of the Hamiltonian (8). Annihilation of very heavy dark matter therefore proceeds in two stages. Size-limited annihilation of the heavy Skyrmions into w particles proceeds into quantum mechanical ww annihilation into Standard Model matter through the Higgs portal coupling (6).
The w particles as the basic field quanta in Eq. (4) can annihilate into Higgs bosons, fermions, and gauge bosons through the coupling L with N c = 1 for leptons and N c = 3 for quarks, and Here δ z = 1 for annihilation into Z bosons and δ z = 0 for annihilation into W + W − . The velocity weighted cross sections are vσ = 2 1 − (4m 2 w /s)σ (s) and the thermally averaged annihilation cross section at temperature T can be calculated using Eq. (1).
The asymptotic cross sections for √ s much larger than the standard model masses are (with δ W = 1 for annihilation and and therefore the net leading order total cross section into standard model states is With Eq. (15) at hand, we can ask for which very high masses m m top and couplings g wh the w particles could account for the dark matter in the universe in a scenario of perturbative thermal creation. The requirement of equality of the thermally averaged cross section (1) with the necessary cross section σ v r (T f ) for thermal particle creation yields the relation between g wh and m w shown in Fig. 1. The w particles could still be thermally created for even higher masses below the unitarity constraint, m w 100 TeV, but this would require non-perturbative coupling and could not be analyzed as easily any more. Non-thermal creation could be possible for even higher mass. However, in all these cases, the w particles would be assumed to be the dominant dark matter component with the Skyrmion solutions to (4) as excited states, and the standard unitarity limits on annihilation cross sections of elementary particles would rule out indirect dark matter signals.
Therefore, we will rely in the following on the assumption of Skyrmion production due to a first order phase transition, whence the Skyrmions dominate the energy density at least early after the phase transition [60], and the w fields only parametrize the condensate giving rise to the Skyrmions. In this case the w particles arise from Skyrmion annihilation, and it is a more intriguing question whether the Higgs portal Skyrmion model could provide an indirect dark matter signal. In scenarios of primary creation of the Skyrmions, the effective balance equation for the total w density n w = 3 i=1 n w i is because the Skyrmion annihilation term will be much larger than the thermal w production term. We use the classical cross section σ SS = π L 2 S for the Skyrmions, noting that the solitons are classical objects of size L S and using the observation that Skyrmions annihilate when their separation compares to their size [94,95]. We also note that they typically annihilate into a small number 2 N w 7 of w particles [94,95], where it is assumed that the mass ratio M S /m w is larger than 2. This theoretical result is in agreement with observations of low-energy p p annihilations. These annihilations produce 5 ± 1 pions, with pion multiplicities 3 ≤ n π ≤ 7 in 99% of measurements [99]. Skyrmion annihilation will therefore inject relativistic w particles with energies M S ≥ E w ≥ m w into the ambient w-particle density, and w particles will annihilate through the Higgs portal with an average center of mass energy squared s , which will satisfy 4m 2 w ≤ s ≤ 4M 2 S . Detailed balance in the equilibrium of very heavy dark matter states then determines the w density, The expected final products in this scenario will be predominantly photons, light leptons and neutrinos from meson decays, because the primary products from heavy Skyrmion annihilation will be Higgs particles and massive electroweak gauge bosons. The Higgs particles at high energies also predominantly decay into Z Z and W + W − , and the massive gauge bosons predominantly decay into hadrons, with a dom-inant meson component which will decay into photons, light leptons and neutrinos. The results of Ref. [100], although derived for lower annihilation energy 10 TeV, are generic for shower composition with annihilation energies well above the top quark mass. The reason is that cross sections into standard model final states simply scale with the relativistic boost factor γ like σ fi ∼ 4/s = E −2 = γ −2 M −2 once E is well above the top quark mass m top . This is a consequence of vdσ fi = V dN fi /T for transitions from an initial 2-particle state |i into a final n f -particle state | f with densities of collision partners V −1 and reaction rate dN fi /T . Once E m top and no more additional Standard Model channels can open up, the Lorentz factors V → V /γ and T → γ T determine the scaling with energy, since v → 1 in a fixed target frame and v → 2 in the center of mass frame. The followup single particle decays scale like d fi = dN fi /T ∼ 1/γ , and therefore dominate the shower formation after the initial annihilation channels of the dark matter particles. Therefore the branching ratios of the initial annihilation events and the subsequent decays are practically fixed once E m top . The CTA is expected to reach a maximum sensitivity vσ 10 −26 cm 3 /s after 500 h for a dark matter mass at the few TeV scale annihilating to W + W − [38], and then the sensitivity limit scales with mass approximately like M 3 for very high dark matter masses. This yields a 500 h sensitivity limit of order vσ 10 −25 cm 3 s −1 × (m/100 TeV) 3 at the 100 TeV mass scale and beyond, whereas the Skyrmion value at that scale is The comparison of very massive Skyrmion dark matter annihilation cross sections with the anticipated CTA sensitivity is displayed in Fig. 2.
Solitonic enhancement factors like the Skyrmion factor L 2 S M 2 S 7.7 × 10 3 g −4 V are certainly needed for dark matter detection at those mass scales.
However, the Higgs portal to dark Skyrmions demonstrates a dark mediator bottleneck. In terms of mass ratios, the lower bound on n w in Eq. (17) becomes For weak scale couplings and not excessively small values of v SS and m w /M S , this implies that n 2 w should at least be comparable to n S n S , such that the Skyrmions would not be the only dark matter component. The Skyrmion annihilation cross sections in Fig. 2 would then have to be reduced by a factor 4n S n S /[(2n S + n w )(2n S + n w )] in comparison to the anticipated CTA sensitivity. This should likely be a generic feature for soliton enhanced dark matter models. If the carrier field w is an inherently stable quantum field, we should always anticipate that w particles form a comparable or even dominant dark matter component for generic mass ratios and coupling constants, even when the solitonic states are the primary states to be created during a phase transition.
One might hope to avoid this dark mediator bottleneck for indirect signals from dark solitons with stable carrier fields w i (x) through Sommerfeld enhancement of the annihilation of the carrier fields. The hope would be that this might prevent build-up of a competing particle dark matter component, and enhance the soliton annihilation signal through fast followup annihilation of the produced w particles. However, this will not remove the bottleneck because of the hard unitarity limit σ ≤ 4π(2 + 1)/k 2 on the -wave annihilation cross sections of particle dark matter of momentum k. This implies the lower limit on the density ratios, and the required small Skyrmion coupling for an observable indirect signal will always imply a noticeable competing particle dark matter component even for small mass ratio m w /M S . Suppression of the very high-energy signal from annihilation of very massive dark Skyrmions due to a competing particle dark matter component appears at least generic, if not even inevitable, in models with stable carriers w.

The neutrino portal to Skyrmion dark matter
For 10 4 N w v SS v ww m 2 w /(g 4 V M 2 S ) 1, the Skyrmions will be a subdominant heavy dark matter component of little observational interest, while at the same time the annihilation cross section for very heavy w particles will be too small to be observable.
A very heavy Skyrmion model with unstable carrier fields w i (x) avoids this dark mediator bottleneck. Unstable w(x) requires breaking of the vector SU (2) in Eq. (4), but introducing a standard model coupling term proportional to a · w(x) with fixed dark isospin vector a will not do the trick. That would destabilize only the component w a (x) =â · w(x), but leave the perpendicular w fields stable, such that again an undetectable w component could dominate the dark matter sector. Stated differently, the dark vector SU(2) has to be broken completely if we wish to avoid any possible remnant dark w component. A natural way to implement such a scenario of unstable carrier field w for the Skyrmions is to use the broken flavor symmetry in the lepton sector through a Majorana type coupling to the right-handed neutrinos, where is the charge conjugate of the right-handed neutrino, and we use C = iγ 2 for the charge conjugation matrix in the Dirac or Weyl representations of the γ matrices. In keeping with the assumption of a very heavy dark sector, we also add a Majorana mass matrix with large eigenvalues, in addition to the standard Yukawa coupling for mass generation for the isospin 1/2 components of the left-handed lepton doublets i = (ν L ,i , e − L ,i ), Here H = (H + , H 0 ) is the Higgs doublet withH = (H 0 * , −H + * ), and v h = 246 GeV is the Higgs vacuum expectation value. This yields according to the standard seesaw mechanism mass eigenstates which are mostly right-handed in the highmass sector and left-handed in the low mass sector, and accounts for the fact that the right-handed neutrinos do not contradict the Planck limit on low mass neutrinos [101]; see e.g. [61,67]. The heavy right-handed neutrinos will remain in thermal equilibrium for temperatures above a few percent of their masses [102,103], and could contribute to dark matter if they are stable, or if the eigenvalues of the mixing matrix m D ∼ √ m ν M R are extremely small to ensure a long lifetime of the heavy right-handed neutrinos, m D 10 −14 eV. However, for generic eigenvalues of m D the right-handed neu-trinos will decay fast into left-handed neutrinos and Higgs particles, and the left-handed neutrinos could then be seen by the neutrino detectors.
The coupling (21) respects the gauge symmetries and Lorentz invariance of the Standard Model, but we cannot rotate the local dark isospin vector with a global vector SU(2) rotation of the w i (x) fields to make only a particular component decay. The coupling does violate lepton number conservation, but it does not require the right-handed neutrinos to be Majorana neutrinos. It only uses the fact that the local symmetries of the standard model are compatible with a Majorana mass matrix for the righthanded neutrinos. In that sense the w i (x) could be thought of as dynamical Majorana masses for the right-handed neutrinos, which could generate such mass terms through a Higgs mechanism. Note also that we can choose the Yukawa couplings symmetric in the neutrino indices, λ i jk = λ jik , since the Dirac structure of the vertices is antisymmetric, The Yukawa couplings can therefore be diagonalized with respect to the neutrino flavors or mass eigenstates, if we wish to do so. It makes sense to refer to the coupling (21) as the ν 2 portal to dark matter, which complements the neutrino(-Higgs) portal introduced in [62,63]. The coupling (21) yields with m w m ν i averaged wdecay rates i.e.
or an averaged lifetime The decay rates and lifetime contain the boost factor γ = 2M S / N w m w because the w particles are created with energy 2M S / N w from Skyrmion annihilation. The lifetime is very small even for very small Yukawa couplings and large boost factors. There will not be a build-up of a dark w component from Skyrmion annihilation in the neutrino portal model, just as there will not be a build-up of heavy right-handed neutrinos for generic mixing matrix m D , and the dark Skyrmions will be the only effective dark matter component in this model. The approximately monochromatic high-energy neutrinos from the two-step process SS → N w w → N w νν(νν) will further decay into left-handed neutrinos and Higgs particles, with a coupling strength relative to the right-handed neutrino production of order ⎛ MeV for neutrino masses m ν < 1 eV [104] and M R 10 TeV. This then corresponds to a lifetime estimate τ R 5 × 10 −20 s for the heavy righthanded neutrinos. The decay of the right-handed neutrinos is slow compared to the decay of the w-particles, but both components are extremely short-lived compared to the stable Skyrmions, which can only decrease in abundance due to annihilation.
Since every high-energy right-handed neutrino decays into a Higgs particle and a high-energy left-handed neutrino, we can estimate the resulting integrated flux of very highenergy neutrinos from Skyrmion annihilation, Evaluation for a Navarro-Frenk-White (NFW) halo [105,106], n S (r) = n S (r) = 4n S (r s )r 3 s r (r + r s ) 2 (27) with r = 8 kpc, a scale radius r s = 10.7 kpc and a halo mass M h = 7.3×10 11 M within 385 kpc for the Milky Way [107] yields the fluxes reported in Fig. 3. The overall diffuse neutrino flux between 250 TeV and 2.5 PeV can be estimated from the IceCube spectral neutrino flux per steradian [16] as j ν 3.3 × 10 −13 cm −2 s −1 . The model therefore appears to be compatible with IceCube observations and the assumption that the very high-energy neutrino flux is dominated by astrophysical sources.

Conclusions
Observation of indirect signals from stable dark matter in the very heavy mass range above 100 TeV requires strongly enhanced annihilation cross sections. Soliton dark matter or bound states can provide a solution to this problem. However, if the field (here a dark isotriplet w) which carries the solitons is stable, small particle annihilation cross sections σ ww ∼ m −2 w will create a dark mediator bottleneck which can prevent an observable indirect dark matter signal. This problem can be avoided if the dark mediator w is unstable, e.g. through the ν 2 portal (21). The neutrino signal for Skyrmion dark matter through the ν 2 portal could contribute at a subdominant, but potentially noticeable level to the flux of very high-energy neutrinos.