Dark Matter Annihilation Research Articles

Complex scalar dark matter in a new gauged U(1) symmetry with kinetic and direct mixings

We propose a scalar dark matter model featuring a hidden gauge symmetry, denoted as U(1)X, with two complex scalars, Φ and S. In this framework, Φ spontaneously breaks the U(1)X gauge symmetry while S serves as a viable dark matter candidate. Particularly, the kinetic and direct mixings between the U(1)X and U(1)Y gauge groups provide a portal between dark matter and the Standard Model particles. These mixings offer a plausible explanation for the W boson mass anomaly observed by the Collider Detector at Fermilab Collaboration. We study the comprehensive phenomenological constraints of this model from colliders and dark matter detection experiments, including Z′ searches at the LHC, the 125 GeV Higgs boson measurements, the relic density of dark matter, and the indirect detection of dark matter annihilation. By randomly scanning the parameter space, we find that the regions where mZ′≳4750 GeV and mZ′≲4750 GeV for gx close to 1 remain viable and can be tested by future experiments. Published by the American Physical Society 2024

ToMCCA: a Toy Monte Carlo Coalescence Afterburner

Antinuclei in our Galaxy may stem either from annihilation or decay of dark matter, or from collisions of cosmic rays with the interstellar medium, which constitute the background of indirect dark matter searches. Understanding the formation mechanism of (anti)nuclei is crucial for setting limits on their production in space. Coalescence models, which describe the formation of light nuclei from final-state interaction of nucleons, have been widely employed in high-energy collisions. In this work, we introduce ToMCCA (Toy Monte Carlo Coalescence Afterburner), which allows for detailed studies of the nuclear formation processes without the overload of general-purpose event generators. ToMCCA contains parameterizations of the multiplicity dependence of the transverse momentum distributions of protons and of the baryon-emitting source size, extracted from ALICE measurements in pp collisions at s=5-13\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s} = 5 - 13$$\\end{document} TeV, as well as of the event multiplicity distributions, taken from the EPOS event generator. ToMCCA provides predictions of the deuteron transverse momentum distributions, with agreement of ∼5%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 5\\%$$\\end{document} with the experimental data. The results of ToMCCA show that the coalescence mechanism in pp collisions depends only on the event multiplicity, not on the collision system or its energy. This allows the model to be utilized for predictions at lower center-of-mass collision energies, which are the most relevant for the production of antinuclei from processes related to dark matter. This model can also be extended to heavier nuclei as long as the target nucleus wavefunction and its Wigner function are known.

Impact of inhomogeneous diffusion on secondary cosmic ray and antiproton local spectra

Recent γ-ray and neutrino observations seem to favor the consideration of non-uniform diffusion of cosmic rays (CRs) throughout the Galaxy.In this study, we investigate the consequences of spatially-dependent inhomogeneous propagation of CRs on the fluxes of secondary CRs and antiprotons detected at Earth. A comparison is made among different scenarios in search of potential features that may guide us toward favoring one over another in the near future. We also examine both the influence of inhomogeneous propagation in the production of secondary CRs from interactions with the gas, and the effects of this scenario on the local fluxes of antiprotons and light antinuclei produced as final products of dark matter annihilation. Our results indicate that the consideration of an inhomogeneous diffusion model could improve the compatibility of the predicted local antiproton flux with that of B, Be and Li, assuming only secondary origin of these particles. In addition, our model predicts a slightly harder local antiproton spectrum, making it more compatible with the high energy measurements of AMS-02.Finally, no significant changes are expected in the predicted local flux of antiprotons and antinuclei produced from dark matter among the different considered propagation scenarios.

Sommerfeld enhancement at NLO and the dark matter unitarity bound

We reexamine the consequences of perturbative unitarity on dark matter freeze-out when both Sommerfeld enhancement and bound state formation affect dark matter annihilations. At leading order the annihilation cross section is infrared dominated and the connection between the unitarity bound and the upper bound on the dark matter mass depends only on how the different partial waves are populated. We compute how this picture is modified at next-to-leading order with the goal of assigning a reliable theory uncertainty to the freeze-out predictions. We explicitly compute next-to-leading order corrections in a simple model with Abelian gauge interactions and provide an estimate of the theoretical uncertainty for the thermal masses of heavy electroweak n-plets. Along the way, we clarify the regularization and matching procedure necessary to deal with singular potentials in quantum mechanics with a calculable, relativistic UV completion. Published by the American Physical Society 2024

Surviving the waves: evidence for a dark matter cusp in the tidally disrupting Small Magellanic Cloud

ABSTRACT We use spectroscopic data for ${\sim }6000$ red giant branch stars in the Small Magellanic Cloud (SMC), together with proper motion data from Gaia Early Data Release 3, to build a mass model of the SMC. We test our Jeans mass modelling method (binulator + gravsphere) on mock data for an SMC-like dwarf undergoing severe tidal disruption, showing that we are able to successfully remove tidally unbound interlopers, recovering the dark matter density and stellar velocity anisotropy profiles within our 95 per cent confidence intervals. We then apply our method to real SMC data, finding that the stars of the cleaned sample are isotropic at all radii (at 95 per cent confidence) and that the inner dark matter density profile is dense, $\rho _{\rm DM}(150\ {\rm pc}) = 1.58_{-0.58}^{+0.80}\times 10^8 \ {\rm M}_{\odot }\, \rm kpc^{-3}$, consistent with a $\Lambda$ cold dark matter cusp. Our model gives a new estimate of the SMC’s total mass within 3 kpc $(M_{\rm tot} \le 3\ {\rm kpc})$ of $2.29\pm 0.46 \times 10^9 \ {\rm M}_{\odot }$. We also derive an astrophysical ‘J-factor’ of $18.99\pm 0.16$ GeV$^2$ cm$^{-5}$ and a ‘D-factor’ of $18.73\pm 0.04$ GeV$^2$ cm$^{-5}$, making the SMC a promising target for dark matter annihilation and decay searches. Finally, we combine our findings with literature measurements to test models in which dark matter is ‘heated up’ by baryonic effects. We find good qualitative agreement with the Di Cintio et al. model but we deviate from the Lazar et al. model at high $M_*/M_{200} > 10^{-2}$. We provide a new, analytical, density profile that reproduces dark matter heating behaviour over the range $10^{-4} < M_*/M_{200} < 10^{-1}$.

Bounds on ALP-mediated dark matter models from celestial objects

We have studied the signals from axion-like particles (ALPs) as dark matter mediators from celestial objects such as neutron stars, brown dwarfs or white dwarfs. We consider the accumulation of dark matter inside the celestial objects using the multiscatter capturing process. The production of ALP from the dark matter annihilation can escape the celestial object and decay into gamma-rays and neutrinos before reaching the Earth. We investigate our model using gamma-ray observations from Fermi and H.E.S.S. and neutrino observations from IceCube and ANTARES. The effective Lagrangian approach allows us to place constraints on the ALP-photon and ALP-fermion couplings. In the gamma-ray channel, our results are able to rule out the existence of ALP with mass up to ~ O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(10) GeV. On the other hand, the neutrino observations can be used to probe a higher mass range with ALP mass up to ~ O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(100) GeV.

High-energy neutrinos from millicharged dark matter annihilation in the Sun

Millicharged dark matter particles can be efficiently captured by the Sun, where they annihilate into tau leptons, leading to the production of high-energy neutrinos. In contrast to the Earth, the high temperature of the Sun suppresses the fraction of millicharged particles that are bound to nuclei, allowing for potentially high annihilation rates. We recast existing constraints from the IceCube Neutrino Observatory and use this information to place new limits on the fraction of the dark matter that is millicharged. This analysis excludes previously unexplored parameter space for masses of mχ∼(5–100) GeV, charges of qχ∼10−3–10−2, and fractional abundances as small as fDM∼10−5. Published by the American Physical Society 2024

Searching accretion-enhanced dark matter annihilation signals in the Galactic Centre

This study reanalyzes the detection prospects of dark matter (DM) annihilation signals in the Galactic Center, focusing on velocity-dependent dynamics within a spike density near the supermassive black hole (Sgr A⋆). We investigate three annihilation processes — p-wave, resonance, and forbidden annihilation — under semi-relativistic velocities, leveraging gamma-ray data from Fermi and DAMPE telescopes. Our analysis integrates a fermionic DM model with an electroweak axion-like particle (ALP) portal, exploring annihilation into two or four photons. Employing a comprehensive six-dimensional integration, we precisely calculate DM-induced gamma-ray fluxes near Sgr A⋆, incorporating velocity and positional dependencies in the annihilation cross-section and photon yield spectra. Our findings highlight scenarios of resonance and forbidden annihilation, where the larger ALP-DM-DM coupling constant Caχχ can affect spike density, potentially yielding detectable gamma-ray line spectra within Fermi and DAMPE energy resolution. We set upper limits for Caχχ across these scenarios, offering insights into the detectability and spectral characteristics of DM annihilation signals from the Galactic Center.

Revisiting the decoupling limit of the Georgi-Machacek model with a scalar singlet

We study the connection between collider and dark matter phenomenology in the singlet extension of the Georgi-Machacek model. In this framework, the singlet scalar serves as a suitable thermal dark matter (DM) candidate. Our focus lies on the region vχ< 1 GeV, where vχ is the common vacuum expectation value of the neutral components of the scalar triplets of the model. Setting bounds on the model parameters from theoretical, electroweak precision and LHC experimental constraints, we find that the BSM Higgs sector is highly constrained. Allowed values for the masses of the custodial fiveplets, triplets and singlet are restricted to the range 140 GeV <MH50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {M}_{H_5^0} $$\\end{document}< 350 GeV, 150 GeV <MH30\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {M}_{H_3^0} $$\\end{document}< 270 GeV and 145 GeV < MH< 300 GeV. The extended scalar sector provides new channels for DM annihilation into BSM scalars that allow to satisfy the observed relic density constraint while being consistent with direct DM detection limits. The allowed region of the parameter space of the model can be explored in the upcoming DM detection experiments, both direct and indirect. In particular, the possible high values of BR(H50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {H}_5^0 $$\\end{document} → γγ) can lead to an indirect DM signal within the reach of CTA. The same feature also provides the possibility of exploring the model at the High-Luminosity run of the LHC. In a simple cut-based analysis, we find that a signal of about 4σ significance can be achieved in final states with at least two photons for one of our benchmark points.

Search for dark matter annihilation to gamma-rays from SPT-SZ selected galaxy clusters

We search for dark matter annihilation from galaxy clusters in the energy range from 1–300 GeV using nearly 16 years of Fermi-LAT data. For this purpose, we use 350 galaxy clusters selected fromthe 2500 deg2 SPT-SZ survey. We model the dark matter distribution using the NFW profile for the main halo along with the Einasto profile for the substructure. The largest signal is seen for the cluster SPT-CL J2021-5257 with a significance of around 3σ. The best-fit dark matter mass and annihilation cross-section for this cluster are equal to (60.0 ± 11.8) GeV and ⟨σv⟩ = (6.0 ± 0.6) × 10-25 cm3 s-1 for the b̅ b annihilation channel. However, this central estimate is in conflict with the limits on annihilation cross-section from dwarf spheroidal galaxies, and hence cannot be attributed to dark matter annihilation. Three other clusters show significance between 2-2.5σ, whereas all the remaining clusters show null results. The most stringent 95% c.l. upper limit for the WIMP annihilation cross-section among all the clusters is from SPT-CL J0455-4159, viz. ⟨σv⟩ = 6.44 × 10-26 cm3 s-1 for m χ = 10 GeV and b b̅ annihilation channel.

Cosmic-ray propagation models elucidate the prospects for antinuclei detection

Tentative observations of cosmic-ray antihelium by the AMS-02 collaboration have re-energized the quest to use antinuclei to search for physics beyond the standard model. However, our transition to a data-driven era requires more accurate models of the expected astrophysical antinuclei fluxes. We use a state-of-the-art cosmic-ray propagation model, fit to high-precision antiproton and cosmic-ray nuclei (B, Be, Li) data, to constrain the antinuclei flux from both astrophysical and dark matter annihilation models. We show that astrophysical sources are capable of producing 𝒪(1) antideuteron events and 𝒪(0.1) antihelium-3 events over 15 years of AMS-02 observations. Standard dark matter models could potentially produce higher levels of these antinuclei, but showing a different energy-dependence. Given the uncertainties in these models, dark matter annihilation is still the most promising candidate to explain preliminary AMS-02 results. Meanwhile, any robust detection of antihelium-4 events would require more novel dark matter model building or a new astrophysical production mechanism.

Significant impact of Galactic dark matter particles on annihilation signals from Sagittarius analogues

We examine the gamma-ray signal from dark matter (DM) annihilation from analogues of the Sagittarius (Sgr) dwarf spheroidal galaxy in the Auriga cosmological simulations. For velocity-dependent annihilation cross sections, we compute emissions from simulated Sgr subhalos and from the Milky Way (MW) foreground. In addition to the annihilation signals from DM particles bound to Sgr, we consider for the first time the annihilation of DM particles bound to the MW that overlap spatially with Sgr. For p-wave models this contribution can enhance the signal by over an order of magnitude, while for d-wave models the enhancement can be over three orders of magnitude. For Sommerfeld and s-wave models, the corresponding emission does not significantly change. For the Sommerfeld model, the Sgr source can be visible above the MW foreground emission, while for s, p and d-wave models, the signal towards Sgr is most likely dominated by foreground MW emission. We interpret our results within the context of the observed gamma-ray emission from Sgr. We find that, given the background emission estimated from this region, the templates from simulations likely have spatial morphology that is too extended to explain the point-like emission that is observed.

Prospects for γ-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Galaxy clusters are expected to be both dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at γ-ray energies and are predicted to be sources of large-scale γ-ray emission due to hadronic interactions in the intracluster medium (ICM).In this paper, we estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse γ-ray emission from the Perseus galaxy cluster.We first perform a detailed spatial and spectral modelling of the expected signal for both the DM and the CRp components. For each case, we compute the expected CTA sensitivity accounting for the CTA instrument response functions. The CTA observing strategy of the Perseus cluster is also discussed.In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio X 500 within the characteristic radius R 500 down to about X 500 < 3 × 10-3, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index αCRp = 2.3. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure αCRp down to about ΔαCRp ≃ 0.1 and the CRp spatial distribution with 10% precision, respectively. Regarding DM, CTA should improve the current ground-based γ-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to ∼ 5, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with τ χ > 1027 s for DM masses above 1 TeV.These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario.

Singlet–doublet Dirac Fermion dark matter from Peccei–Quinn symmetry

Weakly interacting massive particles (WIMPs) and axions are arguably the most compelling dark matter (DM) candidates in the literature. Here, we consider a model where the Peccei–Quinn (PQ) symmetry solves the strong CP problem, generates radiatively Dirac neutrino masses and gives origin to a multicomponent dark sector. Specifically, scotogenic Dirac neutrino masses arise at one-loop level. The lightest fermionic mediator acts as the second DM candidate due to a residual [Formula: see text] symmetry resulting from the PQ symmetry breaking. The WIMP DM component resembles the well-known singlet–doublet fermion DM. While the lower WIMP dark mass region is usually excluded, our model reopens that portion of the parameter space (for DM masses below [Formula: see text][Formula: see text]GeV). Therefore, we perform a phenomenological analysis that addresses the constraints from direct searches of DM, neutrino oscillation data and charged lepton flavor violating processes. The model can be tested in future facilities where neutrino telescopes search for DM annihilation into Standard Model particles.

Dark branes for dark matter

We propose a setup for the origin of dark matter based on spacetime with a warped extra dimension and three branes: the Planck brane, the TeV brane, at a (few) TeV scale ρT, and a dark brane, at a (sub-)GeV scale ρ1≲100 GeV≪ρT. The Standard Model (SM) is localized in the TeV brane, thus solving the Higgs hierarchy problem, while the dark matter χ, a Dirac fermion with mass mχ&lt;ρ1, is localized in the dark brane. The radion, with mass mr&lt;mχ, interacts strongly [∼mχ/ρ1∼O(1)] with dark matter and very weakly (∼mfρ1/ρT2≪1) with the Standard Model matter f. The generic conflict between the bounds on its detection signatures and its proper relic abundance is avoided as dark matter annihilation is p-wave suppressed. The former is determined by its very weak interactions with the SM and the latter by its much stronger annihilation into radions. Therefore, there is a vast range in the dark matter’s parameter space where the correct relic abundance is achieved consistently with the existing bounds. Moreover, for the dark brane with ρ1≲3 GeV, a confinement/deconfinement first order phase transition, where the radion condensates, produces a stochastic gravitational wave background at the nanohertz frequencies, which can be identified with the signal detected by the Pulsar Timing Array (PTA) experiments. In the PTA window, for 0.15≲mχ≲2 GeV the relic abundance is reproduced and all constraints are satisfied. Published by the American Physical Society 2024

Gamma-ray signal from ZN≥3 dark matter-companion models

In Ref. [10], we proposed to replace the final dark matter (DM) particle in the semi-annihilation mode DM+DM→antiDM+Higgs boson with its ZN≥3 companion, thus reducing DM number density without DM-nucleon scattering. In this work, we study the indirect detection signals from DM annihilation, the Higgs boson pair with one of them from the companion decay being on- or off- shell, depending on the DM-companion mass splitting. We generate the photon spectrum by using PYTHIA8 and study the properties of the spectrum, to find that the hard part of the spectrum in our model is mainly shaped by the direct Higgs boson and thus does not differ much from that of the conventional semi-annihilation mode. Using the Fermi-LAT data of dwarf spheroidal galaxies, we derive the current limit of the DM annihilation cross section for DM+DM→companion⁎+Higgs boson, and for the relatively light DM, it reaches the typical thermal cross section. However, for the TeV scale DM, we have to rely on the Cherenkov Telescope Array, which is able to rule out the whole parameter space except for the coannihilation region.

Investigating higgsino dark matter in the semi-constrained NMSSM* *Supported by the National Natural Science Foundation of China (12275066, 11605123)

In this study, we explored the characteristics of higgsino-dominated dark matter (DM) within the semi-constrained Next-to-Minimal Supersymmetric Standard Model (scNMSSM), covering a mass range from hundreds of GeV to several TeV. We carefully analyzed the parameter space under existing theoretical and experimental constraints to confirm the viability of higgsino-dominated lightest supersymmetric particles (LSPs) with masses between 100 GeV and 4 TeV. Our study examined various DM annihilation mechanisms, emphasizing the significant role of coannihilation with the next-to-lightest supersymmetric particle (NLSP), which includes other higgsino-dominated particles such as and . We categorize the annihilation processes into three main classes: coannihilation, Higgs funnel annihilation, and coannihilation. Each class combines interactions with . Our results indicate that achieving the correct relic density in heavier higgsino LSPs requires a combination of coannihilation and Higgs funnel mechanisms. We also assessed the potential of future experiments, such as XENONnT, LUX-ZEPLIN (LZ), PandaX-xT, and the Cherenkov Telescope Array (CTA), to probe these DM scenarios through direct and indirect detections. In particular, future spin-independent DM detections may cover all samples with the correct DM relic density for GeV. Furthermore, future colliders such as the International Linear Collider (ILC) and Compact Linear Collider (CLIC) are expected to exceed the detection capabilities of current hadron colliders, especially for higher mass NLSPs. Notably, CLIC, which will operate at 3000 GeV, is anticipated to enable thorough investigation of all samples with insufficient DM relic density for GeV.

NLO thermal corrections to dark matter annihilation cross sections: A novel approach

Dark matter relic density has been increasingly accurately measured by successive generations of experiments. The Boltzmann equation determines the yields using the dark matter annihilation cross section as one of the inputs; the accurate computation of the latter including thermal contributions thus assumes importance. We report here the next-to-leading-order (NLO) thermal corrections to the cross sections for (Majorana) dark matter annihilation of standard model fermions, χχ→ff¯, via charged scalars. We use a novel approach, utilizing the technique of Grammer and Yennie, extended to thermal field theories, where the cancellation of soft infrared divergences occurs naturally. We present the NLO thermal cross sections in full detail both for the relativistic case and in the nonrelativistic limit. Our independent calculation verifies earlier results where the leading contribution at order O(T2) was shown to be proportional to the square of the fermion mass in the nonrelativistic limit, just as at leading order. We find that the O(T4) contributions have the same dependence on the fermion mass as well. Published by the American Physical Society 2024

Dynamical origin of neutrino masses and dark matter from a new confining sector

A dynamical mechanism, based on a confining non-Abelian dark symmetry, which generates Majorana masses for hyperchargeless fermions, is proposed. We apply it to the inverse seesaw scenario, which allows us to generate light neutrino masses from the interplay of TeV-scale pseudo-Dirac mass terms and a small explicit breaking of lepton number. A single generation of vectorlike dark quarks, transforming under a SU(3)D gauge symmetry, is coupled to a real singlet scalar, which serves as a portal between the dark quark condensate and three generations of heavy sterile neutrinos. Such a dark sector and the Standard Model (SM) are kept in thermal equilibrium with each other via sizable Yukawa couplings to the heavy neutrinos. In this framework, the lightest dark baryon, which has spin 3/2 and is stabilized at the renormalizable level by an accidental dark baryon number symmetry, can account for the observed relic density via thermal freeze-out from annihilations into the lightest dark mesons. These mesons, in turn, decay to heavy neutrinos, which produce SM final states upon decay. This model may be probed by next generation neutrino telescopes via neutrino lines produced from dark matter annihilations. Published by the American Physical Society 2024

Constraints on Ultraheavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations.

In this Letter we try to search for signals generated by ultraheavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible γ rays by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of the LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter that have low fluxes of astrophysical γ-ray background while having large amount of dark matter. By analyzing more than 700days of observational data at LHAASO, no significant dark matter signal from 1TeV to 1EeV is detected. Accordingly we derive the most stringent constraints on the ultraheavy dark matter annihilation cross section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived.