We present a critical assessment of the SN1987A supernova cooling bound on axions and other light particles. Core-collapse simulations used in the literature to substantiate the bound omitted from the calculation the envelope exterior to the proto-neutron star (PNS). As a result, the only source of neutrinos in these simulations was, by construction, a cooling PNS. We show that if the canonical delayed neutrino mechanism failed to explode SN1987A, and if the pre-collapse star was rotating, then an accretion disk would form that could explain the late-time ($t\gtrsim5$ sec) neutrino events. Such accretion disk would be a natural feature if SN1987A was a collapse-induced thermonuclear explosion. Axions do not cool the disk and do not affect its neutrino output, provided the disk is optically-thin to neutrinos, as it naturally is. These considerations cast doubt on the supernova cooling bound.

The time delay measured between the images of gravitationally lensed quasars probes a combination of the angular diameter distance to the source-lens system and the mass density profile of the lens. Observational campaigns to measure such systems have reported a determination of the Hubble parameter H0 that shows significant tension with independent determination based on the cosmic microwave background (CMB) and large scale structure (LSS). We show that lens mass models that exhibit a cored component, coexisting with a cusp, probe a degenerate direction in the lens model parameter space, being an approximate mass sheet transformation. This family of lens models has not been considered by the cosmographic analyses. Once added to the model, the cosmographic error budget should become dependent on stellar kinematics uncertainties. We propose that a core component coexisting with a cusp could bring the lensing measurements of H0 to accord with the CMB/LSS value.

Measurements of the dynamical environment of supermassive black holes (SMBHs) are becoming abundant and precise. We use such measurements to look for ultralight dark matter (ULDM), which is predicted to form dense cores ("solitons") in the centre of galactic halos. We search for the gravitational imprint of an ULDM soliton on stellar orbits near Sgr A* and by combining stellar velocity measurements with Event Horizon Telescope imaging of M87*. Finding no positive evidence, we set limits on the soliton mass for different values of the ULDM particle mass $m$. The constraints we derive exclude the solitons predicted by a naive extrapolation of the soliton-halo relation, found in DM-only numerical simulations, for $2\times10^{-20}~{\rm eV}\lesssim m\lesssim8\times10^{-19}~{\rm eV}$ (from Sgr A*) and $m\lesssim4\times10^{-22}~{\rm eV}$ (from M87*). However, we present theoretical arguments suggesting that an extrapolation of the soliton-halo relation may not be adequate: in some regions of the parameter space, the dynamical effect of the SMBH could cause this extrapolation to over-predict the soliton mass by orders of magnitude.

Analytic arguments and numerical simulations show that bosonic ultra-light dark matter (ULDM) would form cored density distributions (`solitons') at the center of galaxies. ULDM solitons offer a promising way to exclude or detect ULDM by looking for a distinctive feature in the central region of galactic rotation curves. Baryonic contributions to the gravitational potential pose an obstacle to such analyses, being (i) dynamically important in the inner galaxy and (ii) highly non-spherical in rotation-supported galaxies, resulting in non-spherical solitons. We present an algorithm for finding the ground state soliton solution in the presence of stationary non-spherical background baryonic mass distribution. We quantify the impact of baryons on the predicted ULDM soliton in the Milky Way and in low surface-brightness galaxies from the SPARC database.

We study the ratio of neutrino-proton elastic scattering to inverse beta decay event counts, measurable in a scintillation detector like JUNO, as a key observable for identifying the explosion mechanism of a galactic core-collapse supernova. If the supernova is not powered by the core but rather, e.g., by collapse-induced thermonuclear explosion, then a prolonged period of accretion-dominated neutrino luminosity is predicted. Using 1D numerical simulations, we show that the distinct resulting flavour composition of the neutrino burst can be tested in JUNO with high significance, overcoming theoretical uncertainties in the progenitor star profile and equation of state.

Bosonic ultra-light dark matter (ULDM) would form cored density distributions at the center of galaxies. These cores, seen in numerical simulations, admit analytic description as the lowest energy bound state solution ("soliton") of the Schroedinger-Poisson equations. Numerical simulations of ULDM galactic halos found empirical scaling relations between the mass of the large-scale host halo and the mass of the central soliton. We discuss how the simulation results of different groups can be understood in terms of the basic properties of the soliton. Importantly, simulations imply that the energy per unit mass in the soliton and in the virialised host halo should be approximately equal. This relation lends itself to observational tests, because it predicts that the peak circular velocity, measured for the host halo in the outskirts of the galaxy, should approximately repeat itself in the central region. Contrasting this prediction to the measured rotation curves of well-resolved near-by galaxies, we show that ULDM in the mass range $m\sim (10^{-22}\div 10^{-21})$ eV, which has been invoked as a possible solution to the small-scale puzzles of $\Lambda$CDM, is in tension with the data. We suggest that a dedicated analysis of the Milky Way inner gravitational potential could probe ULDM up to $m\lesssim 10^{-19}$ eV.

Cosmic-ray anti-deuterium and anti-helium have long been suggested as probes of dark matter, as their secondary astrophysical production was thought extremely scarce. But how does one actually predict the secondary flux? Anti-nuclei are dominantly produced in pp collisions, where laboratory cross section data is lacking. We make a new attempt at tackling this problem by appealing to a scaling law of nuclear coalescence with the physical volume of the hadronic emission region. The same volume is probed by Hanbury Brown-Twiss (HBT) two-particle correlations. We demonstrate the consistency of the scaling law with systems ranging from central and off-axis AA collisions to pA collisions, spanning 3 orders of magnitude in coalescence yield. Extending the volume scaling to the pp system, HBT data allows us to make a new estimate of coalescence, that we test against preliminary ALICE pp data. For anti-helium the resulting cross section is 1-2 orders of magnitude higher than earlier estimates. The astrophysical secondary flux of anti-helium could be within reach of a five-year exposure of AMS02.

In recent years, space-born experiments have delivered new measurements of high energy cosmic-ray (CR) $\bar p$ and $e^+$. In addition, unprecedented sensitivity to CR composite anti-nuclei anti-d and anti-He is expected to be achieved in the near future. We report on the theoretical interpretation of these measurements. While CR antimatter is a promising discovery tool for new physics or exotic astrophysical phenomena, an irreducible background arises from secondary production by primary CR collisions with interstellar matter. Understanding this irreducible background or constraining it from first principles is an interesting challenge. We review the attempt to obtain such understanding and apply it to CR $\bar p,\, e^+,$ anti-d and anti-He. Based on state of the art Galactic cosmic ray measurements, dominated currently by the AMS-02 experiment, we show that: (i) CR $\bar p$ most likely come from CR-gas collisions; (ii) $e^+$ data is consistent with, and suggestive of the same secondary astrophysical production mechanism responsible for $\bar p$ and dominated by proton-proton collisions. In addition, based on recent accelerator analyses we show that the flux of secondary high energy anti-He may be observable with a few years exposure of AMS-02. We highlight key open questions, as well as the role played by recent and upcoming space and accelerator data in clarifying the origins of CR antimatter.

Secondary astrophysical production of $e^+$ and $\bar{p}$ cosmic rays is considered. Inclusive $\pi$, $K$, and $\bar{p}$ production cross sections in $pp$ collisions at large $\sqrt{s}$ are parametrised using recent experimental data at LHC energies. The astrophysical production rate ratio $Q_{e^+}/Q_{\bar{p}}$ is calculated for an input cosmic ray proton flux consistent with local measurements. At $10<E<100$$\sim$GeV the cosmic ray flux ratio $J_{e^+}/J_{\bar{p}}$ measured by AMS02 falls below the production rate ratio by about 50\%, while at high energy $E>100$$\sim$GeV the measured flux ratio coincides with the production rate ratio of the secondary source.

We revisit cosmic microwave background (CMB) constraints on primordial black hole dark matter. Spectral distortion limits from COBE/FIRAS do not impose a relevant constraint. Planck CMB anisotropy power spectra imply that primordial black holes with $m_{BH}\gtrsim 5~M_{\odot}$ are disfavored. However, this is susceptible to sizeable uncertainties due to the treatment of the black hole accretion process. These constraints are weaker than those quoted in earlier literature for the same observables.

Collapse-induced thermonuclear explosion (CITE) may explain core-collapse supernovae (CCSNe). We present a preliminary analysis of the neutrino signal predicted by CITE and compare it to the neutrino burst of SN1987A. For strong CCSNe, as SN1987A, CITE predicts a proto-neutron star (PNS) accretion phase, accompanied by the corresponding neutrino luminosity, that can last a few seconds and that is cut-off abruptly by black hole (BH) formation. The neutrino luminosity can later be revived by accretion disc emission after a dead time of few to a few ten seconds. In contrast, the neutrino mechanism for CCSNe predicts a shorter PNS accretion phase, followed by a slowly declining PNS cooling luminosity. We repeat statistical analyses used in the literature to interpret the neutrino mechanism, and apply them to CITE. The first 1-2 sec of the neutrino burst are equally compatible with CITE and with the neutrino mechanism. However, the data hints to a luminosity drop at t=2-3 sec, in some tension with the neutrino mechanism while being naturally attributed to BH formation in CITE. The occurrence of neutrino events at 5 sec in SN1987A suggests that the accretion disc formed by that time. We perform 2D numerical simulations, showing that CITE may be able to accommodate this disc formation time while reproducing the ejected $^{56}$Ni mass and ejecta kinetic energy within factors 2-3 of observations. We estimate the disc neutrino luminosity and show that it can roughly match the data. This suggests that direct BH formation is compatible with the neutrino burst of SN1987A. With current neutrino detectors, the neutrino burst of the next Galactic CCSN may give us front-row seats to the formation of an event horizon in real time. Access to phenomena near the event horizon motivates the construction of a few Megaton neutrino detector that should observe extragalactic CCSNe on a yearly basis.

The combination of S-matrix unitarity and the dynamics of thermal freeze-out for massive relic particles (denoted here simply by WIMPs) implies a lower limit on the density of such particles, that provide a (potentially sub-dominant) contribution to dark matter. This then translates to lower limits to the signal rates for a variety of techniques for direct and indirect detection of dark matter. For illustration, we focus on models where annihilation is s-wave dominated. We derive lower limits to the flux of gamma-rays from WIMP annihilation at the Galactic center; direct detection of WIMPs; energetic neutrinos from WIMP annihilation in the Sun; and the effects of WIMPs on the angular power spectrum and frequency spectrum of the cosmic microwave background radiation. The results suggest that a variety of dark-matter-search techniques may provide interesting avenues to seek new physics, even if WIMPs do not constitute all the dark matter. While the limits are quantitatively some distance from the reach of current measurements, they may be interesting for long-range planning exercises.

We study a model in which dark matter couples to the Standard Model through a dilaton of a sector with spontaneously broken approximate scale invariance. Scale invariance fixes the dilaton couplings to the Standard Model and dark matter fields, leaving three main free parameters: the symmetry breaking scale $f$, the dilaton mass $m_{\sigma}$, and the dark matter mass $m_{\chi}$. We analyze the experimental constraints on the parameter space from collider, direct and indirect detection experiments including the effect of Sommerfeld enhancement, and show that dilaton exchange provides a consistent, calculable framework for cold dark matter with $f,\,m_\sigma,\,m_\chi$ of roughly similar magnitude and in the range $\sim 1-10$ TeV. Direct and indirect detection experiments, notably future ground-based gamma ray and space-based cosmic ray measurements, can probe the model all the way to dark matter mass in the multi-TeV regime.

We show that Big Bang Nucleosynthesis (BBN) significantly constrains axion-like dark matter. The axion acts like an oscillating QCD $\theta$ angle that redshifts in the early universe, increasing the neutron-proton mass difference at neutron freeze-out. An axion-like particle that couples too strongly to QCD results in the underproduction of 4He during BBN and is thus excluded. The BBN bound overlaps with much of the parameter space that would be covered by proposed searches for time-varying neutron EDMs. The QCD axion does not couple strongly enough to affect BBN.

We derive constraints on elastic scattering between baryons and dark matter using the cosmic microwave background (CMB) data from the Planck satellite and the Lyman-alpha forest data from the Sloan Digital Sky Survey. Elastic scattering allows baryons and dark matter to exchange momentum, affecting the dynamics of linear density perturbations in the early Universe. We derive constraints to scattering cross sections of the form sigma \propto v^n, allowing for a wide range of velocity dependencies with n between -4 and 2. We improve and correct previous estimates where they exist, including velocity-independent cross section as well as dark matter millicharge and electromagnetic dipole moments. Lyman-alpha forest data dominates the constraints for n>-3, where the improvement over CMB data alone can be several orders of magnitude. Dark matter-baryon scattering cannot affect the halo mass function on mass scales M>10^{12} M_{solar}. Our results imply, model-independently, that a baryon in the halo of a galaxy like our own Milky Way, does not scatter from dark matter particles during the age of the galaxy.

We show that the recent AMS02 positron fraction measurement is consistent with a secondary origin for positrons, and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production (Katz et al 2009), obtained by neglecting radiative losses. This coincidence, which will be further tested by upcoming AMS02 data at higher energy, is a compelling indication for a secondary source. Within the secondary model the AMS02 data imply a cosmic ray propagation time in the Galaxy of < Myr and an average traversed interstellar matter density of order 1/cc, comparable to the density of the Milky Way gaseous disk, at a rigidity of 300 GV.

We show that dark matter annihilation around the time of recombination can lead to growing ionization fraction perturbations, that track the linear collapse of matter over-densities. This amplifies small scale cosmological perturbations to the free electron density by a significant amount compared to the usual acoustic oscillations. Electron density perturbations distort the CMB, inducing secondary non-gaussianity. We calculate the CMB bispectrum from recombination, that is marginally observable by Planck. Even though electron perturbations can be markedly boosted compared with the Standard Model prediction, the dark matter effect in the CMB bispectrum turns out to be small and will be difficult to disentangle from the standard model in the foreseeable future.

We compare the DAMA signal to the muon flux underground, which is annually modulated due to temperature variations in the stratosphere. We show that the muon flux at LNGS and the DAMA signal are tightly correlated. Different mechanisms were proposed in the literature by which muon-induced events may dominate the signal region in DAMA. We discuss simple statistical constraints on such mechanisms and show that the DAMA collaboration can falsify the muon hypothesis, if it is wrong, by reporting their annual baseline count rates.

We take a fresh look at high energy radioactive nuclei data reported in the 90's and at the positron data recently reported by PAMELA. Our aim is to study the model independent implications of these data for the propagation time scales of cosmic rays in the Galaxy. Considering radioactive nuclei, using decaying charge to decayed charge ratios -- the only directly relevant data available at relativistic energies -- we show that a rigidity independent residence time is consistent with observations. The data for all nuclei can be described by f_{s,i}=(t_i/100 Myr)^{0.7}, where f_{s,i} is the suppression of the flux due to decay and t_i is the observer frame lifetime for nucleus specie i. Considering positron measurements, we argue that the positron flux is consistent with a secondary origin. Comparing the positron data with radioactive nuclei at the same energy range, we derive an upper bound on the mean electromagnetic energy density traversed by the positrons, \bar U_T<1.25 eV/cm^3 at a rigidity of R=40 GV. Charge ratio measurements within easy reach of the AMS-02 experiment, most notably a determination of the Cl/Ar ratio extending up to R\sim100 GV, will constrain the energy dependence of the positron cooling time. Such constraints can be used to distinguish between different propagation scenarios, as well as to test the secondary origin hypothesis for the positrons in detail.

We present a critical analysis of the observational constraints on, and of the theoretical modeling of, aspects of cosmic ray (CR) generation and propagation in the Galaxy, which are relevant for the interpretation of recent positron and anti-proton measurements. We give simple, analytic, model independent expressions for the secondary pbar flux, and an upper limit for the secondary e+ flux, obtained by neglecting e+ radiative losses, e+/(e+ + e-)<0.2\pm0.1 up to ~300 GeV. These expressions are completely determined by the rigidity dependent grammage, which is measured from stable CR secondaries up to ~150 GeV/nuc, and by nuclear cross sections measured in the laboratory. pbar and e+ measurements, available up to ~100 GeV, are consistent with these estimates, implying that there is no need for new, non-secondary, pbar or e+ sources. The radiative loss suppression factor f_{s,e+} of the e+ flux depends on the e+ propagation in the Galaxy, which is not understood theoretically. A rough, model independent estimate of f_{s,e+} 1/3 can be obtained at a single energy, E\sim20 GeV, from unstable secondary decay and is found to be consistent with e+ measurements, including the positron fraction measured by PAMELA. We show that specific detailed models, that agree with compositional CR data, agree with our simple expressions for the e+ and pbar flux, and that the claims that the positron fraction measured by PAMELA requires new primary e+ sources are based on assumptions, that are not supported by observations. If PAMELA results are correct, they suggest that f_{s,e+} is slightly increasing with energy, which provides an interesting constraint on CR propagation models. We argue that measurements of the e+ to pbar ratio are more useful for challenging secondary production models than the positron fraction.

We explore, within the warped extra dimensional framework, the possibility of finding anti-matter signals in cosmic rays (CRs) from dark matter (DM) annihilation. Exchange of order 100 GeV radion, an integral part of our setup, generically results in Sommerfeld enhancement of the annihilation rate for TeV DM mass. No dark sector is required to obtain boosted annihilation cross sections. A mild hierarchy between the radion and DM masses can be natural due to the pseudo-Goldstone boson nature of the radion. Implications of Sommerfeld enhancement in warped grand unified theory (GUT) models, where proton stability implies a DM candidate, are studied. We show, via partially unified Pati-Salam group, how to incorporate a custodial symmetry for Z->b\bar b into the GUT framework such that a few TeV Kaluza-Klein (KK) mass scale is allowed by precision tests. The model with smallest fully unified SO(10) representation allows us to decouple the DM from the electroweak sector. Thus, a correct DM relic density is obtained and direct detection bounds are satisfied. Looking at robust CR observables, a possible future signal in the \bar p / p flux ratio is found. We show how to embed a similar custodial symmetry for the right handed tau, allowing it to be strongly coupled to KK particles. Such a scenario might lead to observed signal in CR positrons; however, the DM candidate in this case can not constitute all of the DM in the universe. Independently of the above, the strong coupling between KK particles and tau's can lead to striking LHC signals.

Measurements of flavor ratios of astrophysical neutrino fluxes are sensitive to the two yet unknown mixing parameters $\theta_{13}$ and $\delta$ through the combination $\sin\theta_{13}\cos\delta$. We extend previous studies by considering the possibility that neutrino fluxes from more than a single type of sources will be measured. We point out that, if reactor experiments establish a lower bound on $\theta_{13}$, then neutrino telescopes might establish an upper bound on $|\cos\delta|$ that is smaller than one, and by that prove that CP is violated in neutrino oscillations. Such a measurement requires several favorable ingredients to occur: (i) $\theta_{13}$ is not far below the present upper bound; (ii) The uncertainties in $\theta_{12}$ and $\theta_{23}$ are reduced by a factor of about two; (iii) Neutrino fluxes from muon-damped sources are identified, and their flavor ratios measured with accuracy of order 10% or better. For the last condition to be achieved with the planned km^3 detectors, the neutrino flux should be close to the Waxman-Bahcall bound. It motivates neutrino telescopes that are effectively about 10 times larger than IceCube for energies of O(100 TeV), even at the expense of a higher energy threshold.