We examine velocity-dependent dark matter annihilation in subhalos using a sample of six Milky Way-like galaxies from the Aurgia simulation suite. We quantify the enhancement in the annihilation rate in subhalos when including the contribution from particles in the smooth component of the halo that overlap with the subhalos. The enhancement in the annihilation rate scales with the smooth component of the host halo dark matter density, and is evident for subhalos over the resolvable mass range. Maximal enhancement factors are ti 48 for p-wave models, and ti 37, 000 for d-wave models. For p and d-wave annihilation models, ti 13 and ti 6 subhalos, respectively, across all six host halos have emission from dark matter annihilation in their direction that is above the foreground emission from the smooth dark matter component, and would therefore be resolvable as sources. Such subhalos with the most significant enhancement factors tend to be on the lower end of the mass range and located closer to the center of the host galaxy. We provide a prescription to calculate the enhancement for subhalos as a function of distance from the Galactic center, and use this to examine the impact on dark matter limits from a couple of example dwarf spheroidals. We show that, including the enhancement factors, limits from individual dwarf spheroidals are at a cross section scale that may approach those derived from the Galactic center.

We study the dark matter (DM) annihilation signals from the Large Magellanic Cloud (LMC) and the impact of the LMC on the DM annihilation signals from the Milky Way (MW) halo, using a MW-LMC analogue from the Auriga magneto-hydrodynamical simulations. We find that the gamma-ray signals from DM annihilation from the LMC rises above the MW foreground by a factor of greater than 100 for the s-wave velocity-independent annihilation model, as well as for the Sommerfeld, p-wave, and d-wave velocity-dependent models. We derive upper limits on the annihilation cross section of DM particles in the LMC using Fermi-LAT data for all velocity-dependent cross section models. Bounds for d-wave annihilation are more stringent by similar to 4-6 orders of magnitude relative to previous bounds from dwarf galaxies, and for p-wave emission our bounds are similar to 2-3 orders of magnitude more stringent. We also demonstrate that the impact of the LMC on the DM annihilation signals from the MW halo is greatest for the p-wave and d-wave models towards the outer MW halo, while the impact is minimal for Sommerfeld and s-wave models. The LMC boosts the DM density and velocity distribution in the outer MW halo, both by bringing in high-speed DM particles and by accelerating the DM particles of the MW, affecting the DM annihilation signals from the MW for the p-wave and d-wave models.

Orbits of Milky Way (MW) and M31 satellites in cosmological simulations have been shown to match analytic orbits computed in fine-tuned models of the extended potential of the MW/M31 halos. Using a sample of Local Group (LG) analogs from the Illustris TNG100 simulations, we investigate whether the same is true for LG satellites near the periphery of the LG turnaround radius (∼ 1 Mpc), but outside the virial radii of the MW and M31. From both dark matter-only (TNG100-Dark) and full hydrodynamical (TNG100-1) simulations, we track the orbits of dwarf galaxy analogs and compare them with orbits resulting from an analytic model of the LG potential, which is dominated by the dark matter halos of the MW and M31, along with, where appropriate, their most massive satellites. We find that the reconstructed orbits of the selected outer LG satellites agree with the simulated orbits to within ≲ 10% over the last 6 Gyr, with the scatter increasing with lookback time and reaching ∼ 40% by 6 Gyr. The fact that orbits of outer LG dwarfs track the analytical models suggests that any excess mass in the outer regions of our LG analogs or the tidal field surrounding the LG does not affect the orbits of outer LG satellites. Residual differences between the analytic and simulated orbits are likely due to a combination of mass accretion for the satellites and MW/M31, the non-spherical nature of the MW/M31 potential, and the prescription for dynamical friction in the analytical models.

We investigate the capabilities of upcoming kiloton-scale neutrino detectors, such as Hyper-Kamiokande, in determining the primary cosmic ray spectrum. These detectors provide full-sky coverage and long-term monitoring, unlike traditional satellite and balloon experiments that measure cosmic ray flux at specific altitudes and locations. By analyzing the atmospheric neutrino flux generated by cosmic ray interactions, we demonstrate that future detectors can differentiate between various cosmic ray models with high statistical significance, even when accounting for uncertainties in neutrino cross sections and hadronic interactions. We introduce a technique for reconstructing the primary cosmic ray spectrum using neutrino measurements, which reduces the flux uncertainty from approximately 20% to about 7%. We then show that Hyper-K has the potential to increase sensitivity to neutrino oscillation parameters, such as sin^2θ_23, by a factor of 2. Our results highlight the complementary role of neutrino detectors in cosmic ray physics and their critical importance for precision measurements in particle astrophysics.

The classic model of the Local Group (LG) is that of two dominant constituents, the Milky Way and M31, first separating and then detaching from the Hubble flow, leading to a nearly radial approaching orbit. This simple model has been confronted by new measurements of the 3D M31 kinematics, by cosmological simulations, and by theoretical understanding of the impact of massive substructures such as the Large Magellanic Cloud. This article explores the consequences of new observations and theory on the determination of the mass and dynamics of the LG. The M31 tangential velocity measurement and contribution from the cosmological constant both increase the implied timing mass of the LG to be ∼ 5 × 10^12 M_⊙. Timing mass estimates for the LG tend to be larger than the sum of the Milky Way and M31 halo masses, and larger than independent LG mass estimators. Precision future kinematics have the potential to explore the origin of this difference, shed light on dark matter in the LG, the origin of its angular momentum, and possibly even local values of cosmological parameters.

We present results over an 11-year Solar cycle of cosmic antiprotons based on 1.1×10^{6} events in the rigidity range from 1.00 to 41.9 GV. The p[over ¯] fluxes exhibit distinct properties. The magnitude of the p[over ¯] flux temporal variation is significantly smaller than those of p, e^{-}, and e^{+}. A hysteresis between the p[over ¯] fluxes and the p fluxes is observed, whereas the p[over ¯] and e^{-} fluxes show a linear correlation. With a model-independent analysis, we found a universal relation between the shape of the rigidity spectrum and the magnitude of flux temporal variation over an 11-year Solar cycle for both positively and negatively charged particles. The simultaneous results on p[over ¯] and p, e^{-}, and e^{+} provide unique information for understanding particle transport in the Solar System as a function of mass, charge, and spectral shape.

We report the properties of precision time structures of cosmic nuclei He, Li, Be, B, C, N, and O fluxes over an 11-year solar cycle from May 2011 to November 2022 in the rigidity range from 1.92 to 60.3 GV. The nuclei fluxes show similar but not identical time variations with amplitudes decreasing with increasing rigidity. In particular, below 3.64 GV the Li, Be, and B fluxes, and below 2.15 GV the C, N, and O fluxes, are significantly less affected by solar modulation than the He flux. We observe that these differences in solar modulation are linearly correlated with the differences in the spectral indices of the cosmic nuclei fluxes. This shows, in a model-independent way, that solar modulation of galactic cosmic nuclei depends on their spectral shape. In addition, solar modulation differences due to nuclei velocity dependence on the mass-to-charge ratio (A/Z) are not observed.

Using Gaia Data Release 3 data, we examine the kinematics of the central core of the Sagittarius (Sgr) dwarf spheroidal galaxy using data that include proper motions and line-of-sight velocities for member stars in addition to their projected positions. We extract a sample of bright stars that are high-probability members of Sgr. We model the distances to these stars, which is the only missing phase-space component measurement from our 5D sample, highlighting how their corresponding uncertainties propagate to affect the kinematics. Using line-of-sight velocity data only, which are not affected by the distance uncertainties, and assuming a Jeans-based equilibrium analysis, we obtain a velocity anisotropy of beta(a) = -2.24 +/- 1.99, which implies a system with tangentially biased orbits. With the full 5D data, we project that the data will significantly improve upon measurements of the log-slope of the dark matter density profile and the stellar velocity anisotropy. Tests with mock distance data show an improvement of anisotropy errors of approximately an order of magnitude, and log-slope at the half-light radius of approximately half an order of magnitude.

We analyze the prospects for studying scalar non-standard interactions (SNSI) using the neutrino burst from a Galactic supernova. SNSI modify the resonant flavor conversion and, correspondingly, the neutronization burst signal, and may be identifiable in future multi-tonne-scale experiments such as DUNE. We show that, in the presence of SNSI, neutrinos propagating out of the dense supernova environment acquire a density-squared-dependent contribution to their mass-squared differences, which in turn modifies the energy levels of the neutrino mass eigenstates. This phenomenon is not present in less dense environments like the Earth or the Sun. For a given mass ordering, supernova neutrinos can improve the sensitivity to SNSI parameters by up to four orders of magnitude compared to that achievable with solar or terrestrial neutrino sources.

The recent KM3NeT observation of an O(100 PeV) event KM3-230213A is puzzling because IceCube with much larger effective area times exposure has not found any such events. We propose a novel solution to this conundrum in terms of dark matter (DM) scattering in the Earth's crust. We show that intermediate dark-sector particles that decay into muons are copiously produced when high-energy (∼100 PeV) DM propagates through a sufficient amount of Earth overburden. The same interactions responsible for DM scattering in Earth also source the boosted DM flux from a high-luminosity blazar. We address the non-observation of similar events at IceCube via two examples of weakly coupled long-lived dark sector scenarios that satisfy all existing constraints. We calculate the corresponding dark sector cross sections, lifetimes and blazar luminosities required to yield one event at KM3NeT, and also predict the number of IceCube events for these parameters that can be tested very soon. Our proposed DM explanation of the event can also be distinguished from a neutrino-induced event in future high-energy neutrino flavor analyses, large-scale DM direct detection experiments, as well as at future colliders.

The orbital superposition method originally developed by Schwarzschild (1979) is used to study the dynamics of growth of a black hole and its host galaxy, and has uncovered new relationships between the galaxy's global characteristics. Scientists are specifically interested in finding optimal parameter choices for this model that best match physical measurements along with quantifying the uncertainty of such procedures. This renders a statistical calibration problem with multivariate outcomes. In this article, we develop a Bayesian method for calibration with multivariate outcomes using orthogonal bias functions thus ensuring parameter identifiability. Our approach is based on projecting the posterior to an appropriate space which allows the user to choose any nonparametric prior on the bias function(s) instead of having to model it (them) with Gaussian processes. We develop a functional projection approach using the theory of Hilbert spaces. A finite-dimensional analogue of the projection problem is also considered. We illustrate the proposed approach using a BART prior and apply it to calibrate the Schwarzschild model illustrating how a multivariate approach may resolve discrepancies resulting from a univariate calibration.

We consider the nuclear absorption of dark matter as an alternative to the typical indirect detection search channels of dark matter decay or annihilation. In this scenario, an atomic nucleus transitions to an excited state by absorbing a pseudoscalar dark matter particle and promptly emits a photon as it transitions back to its ground state. The nuclear excitation of carbon and oxygen in the Galactic Center would produce a discrete photon spectrum in the O (10) MeV range that could be detected by gamma-ray telescopes. Using the BIGSTICK large-scale shell-model code, we calculate the excitation energies of carbon and oxygen. We constrain the dark matter-nucleus coupling for current COMPTEL data, and provide projections for future experiments AMEGO-X, e-ASTROGAM, and GRAMS for dark matter masses from similar to 10 to 30 MeV. We find the excitation process to be very sensitive to the dark matter mass and find that the future experiments considered would improve constraints on the dark matter-nucleus coupling within an order of magnitude.

The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on Ar-40 and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called "brems flipping," as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage.

PandaX-4T and XENONnT have recently reported the first measurement of nuclear recoils induced by the 8B solar neutrino flux, through the coherent elastic neutrino-nucleus scattering (CEvNS) channel. As long anticipated, this is an important milestone for dark matter searches as well as for neutrino physics. This measurement means that these detectors have reached exposures such that searches for low mass, less than or similar to 10 GeV dark matter cannot be analyzed using the background-free paradigm going forward. It also opens a new era for these detectors to be used as neutrino observatories. In this paper we assess the sensitivity of these new measurements to new physics in the neutrino sector. We focus on neutrino nonstandard interactions (NSI) and show that-despite the still moderately low statistical significance of the signals- these data already provide valuable information. We find that limits on NSI from PandaX-4T and XENONnT measurements are comparable to those derived using combined COHERENT CsI and LAr data. Furthermore, they provide sensitivity to pure z flavor parameters that are not accessible using stoppedpion or reactor sources. With larger exposures and consequently improvements of statistical uncertainties, forthcoming data from these experiments will provide important, novel results for CEvNS-related physics.

This paper introduces a novel track-length extension fitting algorithm for measuring the kinetic energies of inelastically interacting particles in liquid argon time projection chambers (LArTPCs). The algorithm finds the most probable offset in track length for a track-like object by comparing the measured ionization density as a function of position with a theoretical prediction of the energy loss as a function of the energy, including models of electron recombination and detector response. The algorithm can be used to measure the energies of particles that interact before they stop, such as charged pions that are absorbed by argon nuclei. The algorithm's energy measurement resolutions and fractional biases are presented as functions of particle kinetic energy and number of track hits using samples of stopping secondary charged pions in data collected by the ProtoDUNE-SP detector, and also in a detailed simulation. Additional studies describe the impact of the dE/dx model on energy measurement performance. The method described in this paper to characterize the energy measurement performance can be repeated in any LArTPC experiment using stopping secondary charged pions.

We use new measurements of the M31 proper motion to examine the Milky Way (MW) - M31 orbit and angular momentum. For Local Group (LG) mass consistent with measured values, and assuming the system evolves in isolation, we show a wide range of orbits is possible. We compare to a sample of LG-like systems in the Illustris simulation and find that ∼ 13% of these pairs have undergone a pericentric passage. Using the simulated sample, we examine how accurately an isolated, two-body model describes the MW-M31 orbit, and show that ∼ 10% of the analogues in the simulation are well-modeled by such an orbit. Systems that evolve in isolation by this definition are found to have a lower rate of major mergers and, in particular, have no major mergers since z ≈ 0.3. For all systems, we find an increase in the orbital angular momentum, which is fairly independent of the merger rate and is possibly explained by the influence of tidal torques on the LG. Given the likely quiet recent major merger history of the MW, it is plausible that the isolated two-body model appropriately describes the orbit, though recent evidence for a major merger in M31 may complicate this interpretation.

The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmic ray events collected in the spring of 2021. We use this sample to demonstrate the imaging performance of the charge and light readout systems as well as the signal correlations between the two. We also report argon purity and detector uniformity measurements, and provide comparisons to detector simulations.

We consider the extraction of parameters of dark matter halos hosting ultra-faint dwarf galaxies, in the case where there are only O(10) identified member stars with measured line-of-sight velocities. This scenario is likely to be increasingly common, as upcoming newly discovered dwarf galaxies in the Milky Way, by e.g. the Rubin Observatory, will likely (at least initially) have only a few identified members. Assuming an NFW dark matter profile, equilibrium modeling likely can only robustly extract one halo parameter (ρ_s r_s), but the scale radius itself will typically be unconstrained. In these cases, the results obtainable from Jeans modeling can be well replicated by a simple scaling relation motivated by the half-light mass estimator. As a application, we examine the recently discovered stellar system Ursa Major III, which has been optimistically assessed to have the largest J-factor of any known object. We suggest that, because of the presence of outlier stars, the J-factor obtained from modeling of Ursa Major III is likely inflated, as it is inconsistent with the half-light mass estimator, while removal of the outliers will leave the J-factor unconstrained from below.

Micro-pattern gaseous detectors (MPGDs) are a class of technologies that enable the full three-dimensional spatial reconstruction of ionisation tracks from nuclear and electron recoils in gas. Anticipating near-future 30 m^3-scale time projection chambers with MPGD-based readout, we forecast the sensitivity of such directionally-sensitive low-energy recoil detectors to neutrino interactions beyond the Standard Model. We work in the framework of neutrino non-standard interactions (NSIs), and calculate the combined recoil energy-angle distributions of the electron recoil signal generated by solar neutrinos in atmospheric-pressure He:CF_4 gas. We estimate the expected exclusion limits that such an experiment could place on various NSI parameters, as well as the mass and coupling of a new light mediator that interacts with electrons and neutrinos. We find that with an achievable background reduction of around a factor of ten from current estimates for a 30 m^3 optical-readout detector using this gas mixture, sensitivity to NSI parameters would already approach Borexino's sensitivity. Directionality also allows for event-by-event neutrino energy reconstruction, which would provide a means to resolve some parameter degeneracies present in the modified neutrino cross section in this formalism. Our results strongly motivate the development of small-scale directionally-sensitive gas detectors for neutrino physics.