Thermal dark matter models generically include the prompt creation of highly-concentrated dark matter cusps in the early Universe. Recent studies find that these cusps can survive to the present day, as long as they do not fall into extremely dense regions of baryonic structure. In this work, we build models of dark matter annihilation within the prompt cusps that reside in galaxy clusters, showing that they dominate the total $\gamma$-ray annihilation signal. Using 15 years of Fermi-LAT data, we find no evidence for a $\gamma$-ray excess from these sources, and set strong constraints on annihilating dark matter. These constraints generically rule out the thermal annihilation cross-section to the $b\bar{b}$ channel for dark matter masses below $\sim$200~GeV.

Direct searches for sub-GeV dark matter are limited by the intrinsic quantum properties of the target material. In this proof-of-concept study, we argue that this problem is particularly well suited for machine learning. We demonstrate that a simple neural architecture consisting of a variational autoencoder and a multi-layer perceptron can efficiently generate unique molecules with desired properties. In specific, the energy threshold and signal (quantum) efficiency determine the minimum mass and cross section to which a detector can be sensitive. Organic molecules present a particularly interesting class of materials with intrinsically anisotropic electronic responses and $\mathcal{O}$(few) eV excitation energies. However, the space of possible organic compounds is intractably large, which makes traditional database screening challenging. We adopt excitation energies and proxy transition matrix elements as target properties learned by our network. Our model is able to generate molecules that are not in even the most expansive quantum chemistry databases and predict their relevant properties for high-throughput and efficient screening. Following a massive generation of novel molecules, we use clustering analysis to identify some of the most promising molecular structures that optimise the desired molecular properties for dark matter detection.

Observations have shown that spatially extended "TeV halos" are a common (and potentially generic) feature surrounding young and middle-aged pulsars. However, their morphology is not understood. They are larger than the "compact" region where the stellar remnant dominates the properties of the interstellar medium, but smaller than expected in models of cosmic-ray diffusion through the standard interstellar medium. Several explanations have been proposed, but all have shortcomings. Here, we revisit a class of models where the cosmic-ray gradient produced by the central source induces a streaming stability that "self-confines" the cosmic-ray population. We find that previous studies significantly underpredicted the degree of cosmic-ray confinement and show that corrected models can significantly inhibit cosmic-ray diffusion throughout the TeV halo, especially when similar contributions from the coincident supernova remnant are included.

We compute for the first time baryon number violation at zero temperature from Higgs bubble collisions and find that it can be of the same order as that from thermal sphalerons in the symmetric phase at electroweak temperatures. We study the dependence of the rate of Chern--Simons number transitions on the shape of the scalar potential and on the Lorentz factor of the bubble walls at collision via large-scale (3+1)D lattice simulations of the Higgs doublet and SU(2) gauge fields. We estimate the resulting baryon asymmetry assuming some CP-violating source activated by the Higgs-field variation during the phase transition.

We report an analysis of the interstellar $\gamma$-ray emission in the third Galactic quadrant measured by the {Fermi} Large Area Telescope. The window encompassing the Galactic plane from longitude $210\arcdeg$ to $250\arcdeg$ has kinematically well-defined segments of the Local and the Perseus arms, suitable to study the cosmic-ray densities across the outer Galaxy. We measure no large gradient with Galactocentric distance of the $\gamma$-ray emissivities per interstellar H atom over the regions sampled in this study. The gradient depends, however, on the optical depth correction applied to derive the \HI\ column densities. No significant variations are found in the interstellar spectra in the outer Galaxy, indicating similar shapes of the cosmic-ray spectrum up to the Perseus arm for particles with GeV to tens of GeV energies. The emissivity as a function of Galactocentric radius does not show a large enhancement in the spiral arms with respect to the interarm region. The measured emissivity gradient is flatter than expectations based on a cosmic-ray propagation model using the radial distribution of supernova remnants and uniform diffusion properties. In this context, observations require a larger halo size and/or a flatter CR source distribution than usually assumed. The molecular mass calibrating ratio, $X_{\rm CO} = N({\rm H_{2}})/W_{\rm CO}$, is found to be $(2.08 \pm 0.11) \times 10^{20} {\rm cm^{-2} (K km s^{-1})^{-1}}$ in the Local-arm clouds and is not significantly sensitive to the choice of \HI\ spin temperature. No significant variations are found for clouds in the interarm region.

We perform a frequentist analysis using the standard profile likelihood method for clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). While Bayesian inferences for Effective Field Theory models of galaxy clustering can be highly sensitive to the choice of priors for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare Bayesian and frequentist constraints for the parameter set $\{\sigma_8, H_0, \Omega_{\rm{m}}, w_0, w_a\}$ when fitting to the full-shape of the power spectrum multipoles, the post-reconstruction Baryon Acoustic Oscillation (BAO) measurements, as well as external datasets from the CMB and type Ia supernovae measurements. Bayesian prior effects are very significant for the $w_0w_a$CDM model; while the $1 \sigma$ frequentist confidence intervals encompass the maximum a posteriori (MAP), the Bayesian credible intervals almost always exclude the maximum likelihood estimate (MLE) and the MAP - indicating strong prior volume projection effects - unless supernovae data are included. We observe limited prior effects for the $\Lambda$CDM model, due to the reduced number of parameters. When DESI full-shape and BAO data are jointly fit, we obtain the following $1\sigma$ frequentist confidence intervals for $\Lambda$CDM ($w_0w_a$CDM): $\sigma_8 = 0.867^{+0.048}_{-0.041} , \ H_0 = 68.91^{+0.80}_{-0.79} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.3038\pm0.0110$ ($\sigma_8 = 0.793^{+0.069}_{-0.048} , \ H_0 = 64.9^{+4.8}_{-2.8} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.369^{+0.029}_{-0.059}$ , $w_0 = -0.24^{+0.17}_{-0.64}$ , $w_a = -2.5^{+1.9}_{}$), corresponding to 0.7$\sigma$, 0.3$\sigma$, 0.7$\sigma$ (1.9$\sigma$, 3.4$\sigma$, 5.6$\sigma$, 5.5$\sigma$, 5.6$\sigma$) shifts between the MLE relative to the Bayesian posterior mean for $\Lambda$CDM ($w_0w_a$CDM) respectively.

The Euclid satellite will provide data on the clustering of galaxies and on the distortion of their measured shapes, which can be used to constrain and test the cosmological model. However, the increase in precision places strong requirements on the accuracy of the theoretical modelling for the observables and of the full analysis pipeline. In this paper, we investigate the accuracy of the calculations performed by the Cosmology Likelihood for Observables in Euclid (CLOE), a software able to handle both the modelling of observables and their fit against observational data for both the photometric and spectroscopic surveys of Euclid, by comparing the output of CLOE with external codes used as benchmark. We perform such a comparison on the quantities entering the calculations of the observables, as well as on the final outputs of these calculations. Our results highlight the high accuracy of CLOE when comparing its calculation against external codes for Euclid observables on an extended range of operative cases. In particular, all the summary statistics of interest always differ less than $0.1\,\sigma$ from the chosen benchmark, and CLOE predictions are statistically compatible with simulated data obtained from benchmark codes. The same holds for the comparison of correlation function in configuration space for spectroscopic and photometric observables.

In a Galactic core-collapse supernova (SN), axionlike particles (ALPs) could be emitted via the Primakoff process and eventually convert into $\gamma$ rays in the magnetic field of the Milky Way. From a data-driven sensitivity estimate, we find that, for a SN exploding in our Galaxy, the Fermi Large Area Telescope (LAT) would be able to explore the photon-ALP coupling down to $g_{a\gamma} \simeq 2 \times 10^{-13}\,$GeV$^{-1}$ for an ALP mass $m_a \lesssim 10^{-9}\,$eV. These values are out of reach of next generation laboratory experiments. In this event, the Fermi LAT would probe large regions of the ALP parameter space invoked to explain the anomalous transparency of the Universe to $\gamma$ rays, stellar cooling anomalies, and cold dark matter. If no $\gamma$-ray emission were to be detected, Fermi-LAT observations would improve current bounds derived from SN1987A by more than one order of magnitude.

Dark matter in the Milky Way may annihilate directly into gamma rays, producing a monoenergetic spectral line. Therefore, detecting such a signature would be strong evidence for dark matter annihilation or decay. We search for spectral lines in the Fermi Large Area Telescope observations of the Milky Way halo in the energy range 200 MeV to 500 GeV using analysis methods from our most recent line searches. The main improvements relative to previous works are our use of 5.8 years of data reprocessed with the Pass 8 event-level analysis and the additional data resulting from the modified observing strategy designed to increase exposure of the Galactic center region. We searched in five sky regions selected to optimize sensitivity to different theoretically-motivated dark matter scenarios and find no significant detections. In addition to presenting the results from our search for lines, we also investigate the previously reported tentative detection of a line at 133 GeV using the new Pass 8 data.

The Sun can efficiently capture leptophilic dark matter that scatters with free electrons. If this dark matter subsequently annihilates into leptonic states, it can produce a detectable neutrino flux. Using 10 years of Super-Kamiokande observations, we set constraints on the dark-matter/electron scattering cross-section that exceed terrestrial direct detection searches by more than an order of magnitude for dark matter masses below 100 GeV, and reach cross-sections as low as $\sim$4$\times$10$^{-41}$cm$^{-2}$.

We investigate a nonrelativistic model of tetraquarks, which are assumed to be compact and to consist of diquark-antidiquark pairs. We fit, for the first time, basically all currently known values for the measured masses of 45 mesons, including both charmed and bottom mesons, to the model and predict masses of tetraquarks as well as diquarks. In particular, we find masses of four axial-vector diquarks, i.e., $qc$, $cc$, $qb$, and $bb$, where $q = u,d$, and 24 ground-state tetraquarks, including both heavy-light tetraquarks ($qc\overline{qc}$ and $qb\overline{qb}$) and heavy tetraquarks ($cc\overline{cc}$ and $bb\overline{bb}$). In general, our results for the masses of $qb\overline{qb}$, $cc\overline{cc}$, and $bb\overline{bb}$ are largely comparable with other reported results, whereas our results for the masses of $qc\overline{qc}$ are slightly larger than what has been found earlier. Finally, we identify some of the obtained predictions for masses of tetraquarks with masses of experimental tetraquark candidates, and especially, we find that $\psi(4660)$, $Z_b(10610)$, and $Z_b(10650)$ could be described by the model.

We provide a comprehensive analysis of the phenomenology of axion-like particles (ALPs) produced in core-collapse supernovae (ccSNe) through interactions with electrons and muons, both of which have a non-negligible abundance in the SN plasma. We identify and calculate six significant ALP-production channels, two of which are loop-level processes involving photons. We then examine several observational constraints on the ALP-electron and ALP-muon parameter spaces. Those include the bounds on anomalous cooling, energy deposition, decay into photons, diffuse gamma rays, and the 511 keV line. Our results provide updated and robust constraints on ALP couplings to electrons and muons from an improved treatment of production and absorption processes. Furthermore, we quantify the uncertainties of the results by using three state-of-the-art supernova models based on two independent simulation codes, finding that constraints vary by factors of O(2-10).

We present a cosmic shear consistency analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy). By adopting three tiers of consistency metrics, we compare cosmological constraints between subsets of the KiDS-Legacy dataset split by redshift, angular scale, galaxy colour and spatial region. We also review a range of two-point cosmic shear statistics. With the data passing all our consistency metric tests, we demonstrate that KiDS-Legacy is the most internally consistent KiDS catalogue to date. In a joint cosmological analysis of KiDS-Legacy and DES Y3 cosmic shear, combined with data from the Pantheon+ Type Ia supernovae compilation and baryon acoustic oscillations from DESI Y1, we find constraints consistent with Planck measurements of the cosmic microwave background with $S_8\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.814^{+0.011}_{-0.012}$ and $\sigma_8 = 0.802^{+0.022}_{-0.018}$.

We report on a search for monochromatic $\gamma$-ray features in the spectra of galaxy clusters observed by the \emph{Fermi} Large Area Telescope. Galaxy clusters are the largest structures in the Universe that are bound by dark matter (DM), making them an important testing ground for possible self-interactions or decays of the DM particles. Monochromatic $\gamma$-ray lines provide a unique signature due to the absence of astrophysical backgrounds and are as such considered a smoking-gun signature for new physics. An unbinned joint likelihood analysis of the sixteen most promising clusters using five years of data at energies between 10 and 400 GeV revealed no significant features. For the case of self-annihilation, we set upper limits on the monochromatic velocity-averaged interaction cross section. These limits are compatible with those obtained from observations of the Galactic Center, albeit weaker due to the larger distance to the studied clusters.

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log\, M_{\star}/M_{\odot}>10.5) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at z>5.5. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at z>7.5 agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at z>3.5 halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

We study the capture of WIMP dark matter by the Sun in the non-relativistic effective theory of dark matter self-interactions. The aim is to assess the impact of self-interactions on the expected neutrino flux from the annihilation of WIMPs trapped in the Sun in a model independent manner. We consider all non-relativistic Galilean invariant self-interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to 1 for WIMPs of spin equal to 0, 1/2 and 1. We show that for interaction operators depending at most linearly on the momentum transfer, the WIMP-induced neutrino flux can be enhanced by several orders of magnitude compared to the same flux in absence of self-interactions. This is true even for standard values of the thermally averaged annihilation cross-section. This conclusion impacts the analysis of present and future observations performed at neutrino telescopes.

An excess of $\gamma$ rays in the data measured by the Fermi Large Area Telescope in the direction of the Galactic center has been reported in several publications. This excess, labeled as the Galactic center excess (GCE), is detected analyzing the data with different interstellar emission models, point source catalogs and analysis techniques. The characteristics of the GCE, recently measured with unprecedented precision, are all compatible with dark matter particles (DM) annihilating in the main halo of our Galaxy, even if other interpretations are still not excluded. We investigate the DM candidates that fit the observed GCE spectrum and spatial morphology. We assume a simple scenario with DM annihilating into a single channel but we inspect also more complicated models with two and three channels. We perform a search for a $\gamma$-ray flux from a list of 48 Milky Way dwarf spheroidal galaxies (dSphs) using state-of-the-art estimation of the DM density in these objects. Since we do not find any significant signal from the dSphs, we put upper limits on the annihilation cross section that result to be compatible with the DM candidate that fits the GCE. However, we find that the GCE DM signal is excluded by the AMS-02 $\bar{p}$ flux data for all hadronic and semi-hadronic annihilation channels unless the vertical size of the diffusion halo is smaller than 2 kpc -- which is in tension with radioactive cosmic ray fluxes and radio data. Furthermore, AMS-02 $e^+$ data rule out pure or mixed channels with a component of $e^+ e^-$. The only DM candidate that fits the GCE spectrum and is compatible with constraints obtained with the combined dSphs analysis and the AMS-02 $\bar{p}$ and $e^+$ data annihilates purely into $\mu^+\mu^-$, has a mass of 60 GeV and roughly a thermal cross section.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.