We assess the tension between DAMA/LIBRA and the latest dark matter annual modulation results from the ANAIS-112 and COSINE-100 NaI experiments, under a range of hypotheses ranging from physical to general parameterisations. We find that, in the most physically-motivated cases, the tension between DAMA and these other NaI experiments exceeds 5$\sigma$. Lowering the tension to reasonable values requires significant tuning, such as overfitting with large numbers of free parameters, and opposite-sign modulation between recoil signals on sodium versus iodine.

Dark matter fermions interacting via attractive fifth forces mediated by a light mediator can form dark matter halos in the very early universe. We show that bound systems composed of these halos are capable of generating gravitational wave (GW) signals detectable today, even when the individual halos are very light. The Yukawa force dominates the dynamics of these halo binaries, rather than gravity. As a result, large GW signals can be produced at initially extremely high frequencies, which are then redshifted to frequency bands accessible to current or future GW observatories. In addition, the resulting GW signals carry distinctive features that enable future observations to distinguish them from conventional ones. Notably, even if only a tiny fraction of dark matter experiences strong fifth-force interactions, such effects provide a new avenue to discover self-interacting dark matter through GW observations.

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from $^{39}$Ar beta decays and is suppressed using pulseshape discrimination (PSD). We use two types of PSD algorithm: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the charge of a single photoelectron, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulseshape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected.

The knowledge of scintillation quenching of $\alpha$-particles plays a paramount role in understanding $\alpha$-induced backgrounds and improving the sensitivity of liquid argon-based direct detection of dark matter experiments. We performed a relative measurement of scintillation quenching in the MeV energy region using radioactive isotopes ($^{222}$Rn, $^{218}$Po and $^{214}$Po isotopes) present in trace amounts in the DEAP-3600 detector and quantified the uncertainty of extrapolating the quenching factor to the low-energy region.

We present a framework for relating gravitational wave (GW) sources to the astrophysical properties of spectroscopic galaxy samples. We show how this can enable using clustering measurements of gravitational wave (GW) sources to infer the relationship between the GW sources and the astrophysical properties of their host galaxies. We accomplish this by creating mock GW catalogs from the spectroscopic Sloan Digital Sky Survey (SDSS) DR7 galaxy survey. We populate the GWs using a joint host-galaxy probability function defined over stellar mass, star formation rate (SFR), and metallicity. This probability is modeled as the product of three broken power-law distributions, each with a turnover point motivated by astrophysical processes governing the relation between current-day galaxy properties and BBH mergers, such as galaxy quenching and BBH delay time. Our results show that GW bias is most sensitive to host-galaxy probability dependence on stellar mass, with increases of up to $\sim O (10)\%$ relative to galaxy bias as the stellar mass pivot scale rises. We also find a notable relationship between GW bias and SFR: when the host-galaxy probability favors low-SFR galaxies, the GW bias significantly increases. In contrast, we observe no strong correlation between GW bias and metallicity. These findings suggest that the spatial clustering of GW sources is primarily driven by the stellar mass and SFR of their host galaxies and shows how GW bias measurements can inform models of the host-galaxy probability function.

Sub-GeV dark matter (DM) particles produced via thermal freeze-out evade many of the strong constraints on heavier DM candidates but at the same time face a multitude of new constraints from laboratory experiments, astrophysical observations and cosmological data. In this work we combine all of these constraints in order to perform frequentist and Bayesian global analyses of fermionic and scalar sub-GeV DM coupled to a dark photon with kinetic mixing. For fermionic DM, we find viable parameter regions close to the dark photon resonance, which expand significantly when including a particle-antiparticle asymmetry. For scalar DM, the velocity-dependent annihilation cross section evades the strongest constraints even in the symmetric case. Using Bayesian model comparison, we show that both asymmetric fermionic DM and symmetric scalar DM are preferred over symmetric fermionic DM due to the reduced fine-tuning penalty. Finally, we explore the discovery prospects of near-future experiments both in the full parameter space and for specific benchmark points. We find that the most commonly used benchmark scenarios are already in tension with existing constraints and propose a new benchmark point that can be targeted with future searches.

Dark matter accumulates in the center of the Earth as the planet plows through the dark matter halo in the Milky Way. Possible annihilation of dark matter to Standard Model particles can be probed in indirect dark matter searches. Among possible messengers, neutrinos are uniquely ideal as they can escape dense regions. Therefore, neutrino telescopes, with their large volume and broad energy exposures, offer new opportunities to search for dark matter signals from the center of the Earth. However, such studies have been restricted to dark matter masses below $\sim$ PeV as the Earth becomes opaque to very-high-energy neutrinos. In this study, we demonstrate that neutrino telescopes operating at TeV-PeV energies can probe very heavy dark matter particles if they annihilate to tau neutrinos or tau leptons. Here, we report upper limits on the spin-independent dark matter-nucleon cross section for masses between $10^5$ GeV and $10^{10}$ GeV by using 7.5 years of IceCube high-energy starting event observations. Our results motivate detailed analyses in IceCube and other upcoming neutrino telescopes in the Northern Hemisphere.

Many of the tensions in cosmological models of the Universe lie in the low mass, low velocity regime. Probing this regime requires a statistically significant sample of galaxies with well measured kinematics and robustly measured uncertainties. WALLABY, as a wide area, untargetted HI survey is well positioned to construct this sample. As a first step towards this goal we develop a framework for testing kinematic modelling codes in the low resolution, low $S/N$, low rotation velocity regime. We find that the WALLABY Kinematic Analysis Proto-Pipeline (WKAPP) is remarkably successful at modelling these galaxies when compared to other algorithms, but, even in idealized tests, there are a significant fraction of false positives found below inclinations of $\approx 40^{\circ}$. We further examine the 11 detections with rotation velocities below $50~\kms$ in the WALLABY pilot data releases. We find that those galaxies with inclinations above $40^{\circ}$ lie within $1-2~\sigma$ of structural scaling relations that require reliable rotation velocity measurements, such as the baryonic Tully Fisher relation. Moreover, the subset that have consistent kinematic and photometric inclinations tend to lie nearer to the relations than those that have inconsistent inclination measures. This work both demonstrates the challenges faced in low-velocity kinematic modelling, and provides a framework for testing modelling codes as well as constructing a large sample of well measured low rotation models from untargetted surveys.

Dark matter particles with Planck-scale mass ($\simeq10^{19}\text{GeV}/c^2$) arise in well-motivated theories and could be produced by several cosmological mechanisms. Using a blind analysis of data collected over a 813 d live time with DEAP-3600, a 3.3 t single-phase liquid argon-based dark matter experiment at SNOLAB, a search for supermassive dark matter was performed, looking for multiple-scatter signals. No candidate signal events were observed, leading to the first direct detection constraints on Planck-scale mass dark matter. Leading limits constrain dark matter masses between $8.3\times10^{6}$ and $1.2\times10^{19} \text{GeV}/c^2$, and cross sections for scattering on $^{40}$Ar between $1.0\times10^{-23}$ and $2.4\times10^{-18} \text{cm}^2$. These are used to constrain two composite dark matter models.

We determine the upper limit on the mass of the lightest neutrino from the most robust recent cosmological and terrestrial data. Marginalizing over possible effective relativistic degrees of freedom at early times ($N_\mathrm{eff}$) and assuming normal mass ordering, the mass of the lightest neutrino is less than 0.037 eV at 95% confidence; with inverted ordering, the bound is 0.042 eV. These results improve upon the strength and robustness of other recent limits and constrain the mass of the lightest neutrino to be barely larger than the largest mass splitting. We show the impacts of realistic mass models, and different sources of $N_\mathrm{eff}$.

We study the electrodynamics of a kinetically mixed dark photon cloud that forms through superradiance around a spinning black hole, and design strategies to search for the resulting multimessenger signals. A dark photon superradiance cloud sources a rotating dark electromagnetic field which, through kinetic mixing, induces a rotating visible electromagnetic field. Standard model charged particles entering this field initiate a transient phase of particle production that populates a plasma inside the cloud and leads to a system which shares qualitative features with a pulsar magnetosphere. We study the electrodynamics of the dark photon cloud with resistive magnetohydrodynamics methods applicable to highly magnetized plasma, adapting techniques from simulations of pulsar magnetospheres. We identify turbulent magnetic field reconnection as the main source of dissipation and electromagnetic emission, and compute the peak luminosity from clouds around solar-mass black holes to be as large as $10^{43}$ erg/s for open dark photon parameter space. The emission is expected to have a significant X-ray component and is potentially periodic, with period set by the dark photon mass. The luminosity is comparable to the brightest X-ray sources in the Universe, allowing for searches at distances of up to hundreds of Mpc with existing telescopes. We discuss observational strategies, including targeted electromagnetic follow-ups of solar-mass black hole mergers and targeted continuous gravitational wave searches of anomalous pulsars.

A new approach is presented to compute entropy for massless scalar quantum fields. By perturbing a skewed correlation matrix composed of field operator correlation functions, the mutual information is obtained for disjoint spherical regions of size $r$ at separation $R$, including an expansion to all orders in $r/R$. This approach also permits a perturbative expansion for the thermal field entropy difference in the small temperature limit ($T \ll 1/r$).

A new comprehensive study on the Cs${_2}$ZrCl${_6}$ (CZC) crystal scintillating properties under different types of irradiation was performed over a wide temperature range from 5 to 300 K. The light yield (LY) at room temperature (RT), measured under irradiation by 662 keV $\gamma$ quanta of $^{137}$Cs, was evaluated to be 53,300$\pm$4,700 photons/MeV corresponding to approximately 71% of its estimated absolute value. The maximum light emission was observed in the temperature interval 135-165 K, where the LY reached 56,900 photons/MeV and 19,700 photons/MeV for $\gamma$ quanta and $\alpha$ particles, respectively. The quenching factor (QF) for $\alpha$ particles increases smoothly from QF = 0.30 at RT to QF = 0.36 at 135 K. The shape of scintillation pulses induced by $\alpha$ particles is characterized by three time-constants (0.3, 2.5 and 11.8 $\mu$s at RT), whereas the average pulse of $\gamma$ induced events is characterized by two time-constants (1.3 and 11.5 $\mu$s at RT). At the same time, scintillating properties and pulse-shape discrimination capability of the CZC exhibit an acute deterioration at temperatures below 135 K. The optimal operating conditions to maximize the scintillating performance of undoped CZC crystals are discussed.

$^{180m}$Ta is the longest-lived metastable state presently known. Its decay has not been observed yet. In this work, we report a new result on the decay of \mTa obtained with a $2015.12$-g tantalum sample measured for $527.7$ d with an ultra-low background HPGe detector in the STELLA laboratory of the Laboratori Nazionali del Gran Sasso (LNGS), in Italy. Before the measurement, the sample has been stored deep-underground for ten years, resulting in subdominant background contributions from cosmogenically activated $^{182}$Ta. We observe no signal in the regions of interest and set half-life limits on the process for the two channels EC and $\beta^-$: $T_{1/2,~\mathrm{EC}} > 1.6 \times 10^{18}$ yr and $T_{1/2,~\beta^-} > 1.1\times 10^{18}$ yr ($90$\% C.\,I.), respectively. We also set the limit on the $\gamma$ de-excitation / IC channel: $T_{1/2,~\mathrm{IC}} > 4.1 \times 10^{15}$ yr ($90$\% C.\,I.). These are, as of now, the most stringent bounds on the decay of $^{180m}$Ta worldwide. Finally, we test the hypothetical scenarios of de-excitation of $^{180m}$Ta by cosmological Dark Matter and constrain new parameter space for strongly-interacting dark-matter particle with mass up to $10^5$ GeV.

Scalar boson stars have attracted attention as simple models for exploring the nonlinear dynamics of a large class of ultra compact and black hole mimicking objects. Here, we study the impact of interactions in the scalar matter making up these stars. In particular, we show the pivotal role the scalar phase and vortex structure play during the late inspiral, merger, and post-merger oscillations of a binary boson star, as well as their impact on the properties of the merger remnant. To that end, we construct constraint satisfying binary boson star initial data and numerically evolve the nonlinear set of Einstein-Klein-Gordon equations. We demonstrate that the scalar interactions can significantly affect the inspiral gravitational wave amplitude and phase, and the length of a potential hypermassive phase shortly after merger. If a black hole is formed after merger, we find its spin angular momentum to be consistent with similar binary black hole and binary neutron star merger remnants. Furthermore, we formulate a mapping that approximately predicts the remnant properties of any given binary boson star merger. Guided by this mapping, we use numerical evolutions to explicitly demonstrate, for the first time, that rotating boson stars can form as remnants from the merger of two non-spinning boson stars. We characterize this new formation mechanism and discuss its robustness. Finally, we comment on the implications for rotating Proca stars.

Neutron stars provide a compelling testing ground for gravity, nuclear dynamics, and physics beyond the Standard Model, and so it will be useful to locate the neutron stars nearest to Earth. To that end, we revisit pulsar distance estimates extracted from the dispersion measure of pulsar radio waves scattering on electrons. In particular, we create a new electron density map for the local kiloparsec by fitting to parallax measurements of the nearest pulsars, which complements existing maps that are fit on the Galactic scale. This ``near-Earth'' electron density map implies that pulsars previously estimated to be 100-200 pc away may be as close as tens of parsecs away, which motivates a parallax-based measurement campaign to follow-up on these very-near candidate pulsars. Such nearby neutron stars would be valuable laboratories for testing fundamental physics phenomena, including several late-stage neutron star heating mechanisms, using current and forthcoming telescopes. We illustrate this by estimating the sensitivities of the upcoming Extremely Large Telescope and Thirty Meter Telescope to neutron stars heated by dark matter capture.

The half-life of $^{39}$Ar is measured using the DEAP-3600 detector located 2 km underground at SNOLAB. In 2016-2020, DEAP-3600 used a target mass of (3269 $\pm$ 24) kg of liquid argon distilled from the atmosphere in a direct-detection dark matter search. Such an argon mass also enables direct measurements of argon isotope properties. The decay of $^{39}$Ar in DEAP-3600 is the dominant source of triggers by two orders of magnitude, ensuring high statistics and making DEAP-3600 well-suited for measuring this isotope's half-life. Use of the pulse-shape discrimination technique in DEAP-3600 allows for powerful discrimination between nuclear recoils and electron recoils, resulting in the selection of a clean sample of $^{39}$Ar decays. Observing over a period of 3.4 years, the $^{39}$Ar half-life is measured to be $(302 \pm 8_{\rm stat} \pm 6_{\rm sys})$ years. This new direct measurement suggests that the half-life of $^{39}$Ar may be significantly longer than the accepted value, with potential implications for measurements using this isotope's half-life as input.

The flavor composition of TeV--PeV astrophysical neutrinos, i.e., the proportion of neutrinos of different flavors in their flux, is a versatile probe of high-energy astrophysics and fundamental physics. Because flavor identification is challenging and the number of detected high-energy astrophysical neutrinos is limited, so far measurements of the flavor composition have represented an average over the range of observed neutrino energies. Yet, this washes out the potential existence of changes in the flavor composition with energy and weakens our sensitivity to the many models that posit them. For the first time, we measure the energy dependence of the flavor composition, looking for a transition from low to high energies. Our present-day measurements, based on the 7.5-year public sample of IceCube High-Energy Starting Events (HESE), find no evidence of a flavor transition. The observation of HESE and through-going muons jointly by next-generation neutrino telescopes Baikal-GVD, IceCube-Gen2, KM3NeT, P-ONE, TAMBO, and TRIDENT may identify a flavor transition around 200TeV by 2030. By 2040, we could infer the flavor composition with which neutrinos are produced with enough precision to establish the transition from neutrino production via the full pion decay chain at low energies to muon-damped pion decay at high energies.

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.