Strong evidence for an isotropic, Gaussian gravitational wave background (GWB) has been found by multiple pulsar timing arrays (PTAs). The GWB is expected to be sourced by a finite population of supermassive black hole binaries (SMBHBs) emitting in the PTA sensitivity band, and astrophysical inference of PTA data sets suggests a GWB signal that is at the higher end of GWB spectral amplitude estimates. However, current inference analyses make simplifying assumptions, such as modeling the GWB as Gaussian, assuming that all SMBHBs only emit at frequencies that are integer multiples of the total observing time, and ignoring the interference between the signals of different SMBHBs. In this paper, we build analytical and numerical models of an astrophysical GWB without the above approximations, and compare the statistical properties of its induced PTA signal to those of a signal produced by a Gaussian GWB. We show that interference effects introduce non-Gaussianities in the PTA signal, which are currently unmodeled in PTA analyses.

Ultralight dark matter may couple quadratically to Standard Model particles. Such quadratic interactions give rise to both coherent and stochastic signals in pulsar timing array (PTA) observations. In this work, we characterize these signals, including the effects of dark matter propagation in a finite-density medium, and assess the sensitivity of current and upcoming PTA observations to their detection. For coherent signals, we find that the sensitivity of current PTA observations competes with and sometimes exceeds that of other probes, such as equivalence principle tests and atomic clocks. For stochastic signals, we find that PTA sensitivities underperform equivalence principle constraints for both existing and upcoming PTA data sets.

With pulsar timing arrays (PTAs) having observed a gravitational wave background (GWB) at nanohertz frequencies, the focus of the field is shifting towards determining and characterizing its origin. While the primary candidate is a population of GW-emitting supermassive black hole binaries (SMBHBs), many other cosmological processes could produce a GWB with similar spectral properties as have been measured. One key argument to help differentiate an SMBHB GWB from a cosmologically sourced one is its level of anisotropy; a GWB sourced by a finite population will likely exhibit greater anisotropy than a cosmological GWB through finite source effects (“shot noise”) and potentially large-scale structure. Current PTA GWB anisotropy detection methods often use the frequentist PTA optimal statistic for its fast estimation of pulsar pair correlations and relatively low computational overhead compared to spatially-correlated Bayesian analyses. However, there are critical limitations with the status quo approach. In this paper, we improve this technique by incorporating three recent advancements: accounting for covariance between pulsar pairwise estimates of correlated GWB power; the per-frequency optimal statistic to dissect the GWB across the spectrum; and constructing null-hypothesis statistical distributions that include cosmic variance. By combining these methods, our new pipeline can localize GWB anisotropies to specific frequencies, through which anisotropy detection prospects — while impacted by cosmic variance — are shown to improve in our simulations from a pp-value of ∼0.2∼0.2 in a broadband search to ∼0.01∼0.01 in the per-frequency search. Our methods are already incorporated in community-available code and ready to deploy on forthcoming PTA datasets.

Recent findings from several Pulsar Timing Array (PTA) collaborations point to the existence of a Gravitational Wave Background (GWB) at nanohertz frequencies. A key next step towards characterizing this signal and identifying its origin is to map the sky distribution of its power. Several strategies have been proposed to reconstruct this distribution using PTA data. In this work, we compare these different strategies to determine which one is best suited to detect GWB anisotropies of different topologies. We find that, for both localized and large-scale anisotropies, reconstruction methods based on pixel and radiometer maps are the most promising. However, in both scenarios, even the optimistically large anisotropic signals discussed in this work remain challenging to detect with near-future PTA sensitivities. For example, we find that for a GWB hotspot contributing to $80\mathrm{\%}$of the GWB power in the second frequency bin, detection probabilities reach at most $\mathcal{O}(10\mathrm{\%})$ for a PTA with noise properties comparable with the ones of the upcoming IPTA third data release. Finally, we consider the fundamental limitations that cosmic variance poses to these kinds of searches by deriving the smallest deviations from isotropy that could be detected by an idealized PTA with no experimental or pulsar noise.

Several Pulsar Timing Array (PTA) collaborations have recently found evidence for a Gravitational Wave Background (GWB) by measuring the perturbations that this background induces in the time-of-arrivals of pulsar signals. These perturbations are expected to be correlated across different pulsars and, for isotropic GWBs, the expected values of these correlations (obtained by averaging over different GWB realizations) are a simple function of the pulsars’ angular separations, known as the Hellings-Downs (HD) correlation function. On the other hand, anisotropic GWBs would induce deviations from these HD correlations in a way that can be used to search for anisotropic distributions of the GWB power. However, even for isotropic GWBs, interference between GW sources radiating at overlapping frequencies induces deviations from the HD correlation pattern, an effect known in the literature as “cosmic variance”. In this work, we study the impact of cosmic variance on PTA anisotropy searches. We find that the fluctuations in cross-correlations related to cosmic variance can lead to the miss-classification of isotropic GWBs as anisotropic, leading to a false detection rate of ~50% for frequentist anisotropy searches. We also observe that cosmic variance complicates the reconstruction of the GWB sky map, making it more challenging to resolve bright GW hotspots, like the ones expected to be produced from a Supermassive Black Hole Binaries population. These findings highlight the need to refine anisotropy search techniques to improve our ability to reconstruct the GWB sky map and accurately assess the significance of any isotropy deviations we might find in it.

Several pulsar timing array (PTA) collaborations have recently found evidence for a gravitational wave background (GWB) permeating our Galaxy. The origin of this background is still unknown. Indeed, while the gravitational wave emission from inspiraling supermassive black hole binaries (SMBHB) is the primary candidate for its origin, several cosmological sources have also been proposed. One distinctive feature of SMBHB-generated backgrounds is the presence of GWB anisotropies stemming from the binaries distribution in the local Universe. However, none of the anisotropy searches performed to date reported a detection. In this work, we show that the lack of anisotropy detection is not currently in tension with a SMBHB origin of the background. We accomplish this by calculating the probability for present and future PTAs to observe deviations from an isotropic GWB. We find that a PTA with the noise characteristics of the NANOGrav 15-year data set had only a 2%–11% probability of detecting deviations from isotropy in an SMBHB-generated GWB, depending on the properties of the SMBHB population. However, we estimate that for the IPTA DR3 data set these probabilities will increase to 4%–28%, putting more pressure on the SMBHB interpretation in case of a null detection. We also identify SMBHB populations that are more likely to produce significant deviations from isotropy. This information could be used together with the spectral properties of the GWB to characterize the SMBHB population.

Pulsar timing array (PTA) searches for gravitational waves (GWs) aim to detect a characteristic correlation pattern in the timing residuals of galactic millisecond pulsars. This pattern is described by the PTA overlap reduction function (ORF) ${\mathrm{\Gamma }}_{ab}\left({\xi }_{ab}\right)$, which is known as the Hellings–Downs (HD) curve in general relativity (GR). In theories of modified gravity, the HD curve often receives corrections. Assuming, e.g. a subluminal GW phase velocity, one finds a drastically enhanced ORF in the limit of small angular separations between pulsar a and pulsar b in the sky, ${\xi }_{ab}\to 0$. In particular, working in harmonic space and performing an approximate resummation of all multipole contributions, the auto correlation coefficient Γ${}_{aa}$ seems to diverge. In this paper, we confirm that this divergence is unphysical and provide an exact and analytical expression for Γ${}_{aa}$ in dependence of the pulsar distance La and the GW phase velocity $v_{\mathrm{p}\mathrm{h}}$vph​. In the GR limit and assuming a large pulsar distance, our expression reduces to ${\mathrm{\Gamma }}_{aa}=1$. In the case of subluminal phase velocity, we show that the regularization of the naive divergent result is a finite-distance effect, meaning that Γ${}_{aa}$ scales linearly with fLa, where f is the GW frequency. For superluminal phase velocity (subluminal group velocity), which is relevant in the case of massive gravity, we correct an earlier analytical result for Γ${}_{ab}$b​. Our results pave the way for fitting modified-gravity theories with nonstandard phase velocity to PTA data, which requires a proper understanding of the auto correlation coefficient Γ${}_{aa}$.

Metric perturbations induced by ultralight dark matter (ULDM) fields have long been identified as a potential target for pulsar timing array (PTA) observations. Previous works have focused on the coherent oscillation of metric perturbations at the characteristic frequency set by the ULDM mass. In this work, we show that ULDM fields source low-frequency stochastic metric fluctuations and that these low-frequency fluctuations can produce distinctive detectable signals in PTA data. Using the NANOGrav 12.5-yr dataset and synthetic datasets mimicking present and future PTA capabilities, we show that the current and future PTA observations provide the strongest probe of ULDM density within the Solar System for masses in the range of ${10}^{–18}\mathrm{eV}-{10}^{–16}\mathrm{eV}$.

Pulsar timing arrays have found evidence for a low-frequency gravitational-wave background (GWB). Assuming that the GWB is produced by supermassive black hole binaries (SMBHBs), the next gravitational-wave (GW) signals astronomers anticipate are continuous waves (CWs) from single SMBHBs and their associated GWB anisotropy. The prospects for detecting CWs and anisotropy are highly dependent on the astrophysics of SMBHB populations. Thus, information from single sources can break degeneracies in astrophysical models and place much more stringent constraints than the GWB alone. We simulate and evolve SMBHB populations, model their GWs, and calculate their anisotropy and detectability. We investigate how varying components of our semianalytic model, including the galaxy stellar mass function, the SMBH–host galaxy relation (M${}_{BH}$​–M${}_{bulge}$​), and the binary evolution prescription, impact the expected detections. The CW occurrence rate is greatest for few total binaries, high SMBHB masses, large scatter in M${}_{BH}$​–M${}_{bulge}$​, and long hardening times. The occurrence rate depends most on the binary evolution parameters, implying that CWs offer a novel avenue to probe binary evolution. The most detectable CW sources are in the lowest frequency bin for a 16.03 yr PTA, have masses from ∼10${}^9$ to 10${}^{10}$M${}_{\odot }$​, and are ∼1 Gpc away. The level of anisotropy increases with frequency, with the angular power spectrum over multipole modes ℓ varying in low-frequency C${}_{\mathrm{\ell }>0}$/C${}_0$​ from ∼5 × 10${}^{\text{−}3}$ to ∼2 × 10${}^{\text{−}1}$, depending on the model; typical values are near current upper limits. Observing this anisotropy would support SMBHB models for the GWB over cosmological models, which tend to be isotropic.

Single phonon excitations, with energies in the 1–100 meV range, are a powerful probe of light dark matter (DM). Utilizing effective field theory, we derive a framework to compute DM absorption rates into single phonons starting from general DM-electron, proton, and neutron interactions. We apply the framework to a variety of DM models: Yukawa coupled scalars, axionlike particles with derivative interactions, and vector DM coupling via gauge interactions or Standard Model electric and magnetic dipole moments. We find that GaAs or ${\mathrm{Al}}_2{\mathrm{O}}_3$ targets can set powerful constraints on a $U(1{)}_{B–L}$ model, and targets with electronic spin ordering are similarly sensitive to DM coupling to the electron magnetic dipole moment. Lastly, we make the code, phonodark-abs (an extension of the existing phonodark code which computes general DM–single phonon scattering rates), publicly available.

This is a lightweight manual for PTArcade, a wrapper of ENTERPRISE and ceffyl that allows for easy implementation of new-physics searches in PTA data. In this manual, we describe how to get PTArcade installed (either on your local machine or an HPC cluster). We discuss how to define a stochastic or deterministic signal and how PTArcade implements these signals in PTA-analysis pipelines. Finally, we show how to handle and analyze the PTArcade output using a series of utility functions that come together with PTArcade.

The 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.

We present a recast in different benchmark models of the recent CMS search that uses the end cap muon detector system to identify displaced showers produced by decays of long-lived particles (LLPs). The exceptional shielding provided by the steel between the stations of the muon system drastically reduces the Standard Model background that limits other existing ATLAS and CMS searches. At the same time, by using the muon system as a sampling calorimeter, the search is sensitive to LLPs energies rather than masses. We show that, thanks to these characteristics, this new search approach is sensitive to LLPs masses even lighter than a GeV, and can be complementary to proposed and existing dedicated LLP experiments.

Pulsar timing arrays (PTAs) are exceptionally sensitive detectors in the frequency band $\mathrm{nHz}\lesssim f\lesssim \mu \mathrm{Hz}$. Ultralight dark matter (ULDM), with mass in the range ${10}^{–23}\mathrm{eV}\lesssim m_{\varphi }\lesssim {10}^{–20}\mathrm{eV}$, is one class of DM models known to generate signals in this frequency window. While purely gravitational signatures of ULDM have been studied previously, in this work we consider two signals in PTAs which arise in the presence of direct couplings between ULDM and ordinary matter. These couplings induce variations in fundamental constants, i.e., particle masses and couplings. These variations can alter the moment of inertia of pulsars, inducing pulsar spin fluctuations via conservation of angular momentum, or induce apparent timing residuals due to reference clock shifts. By using mock data mimicking current PTA datasets, we show that PTA experiments outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. In the case of coupling to quarks or photons, we find that PTAs and atomic clocks set similar constraints. Additionally, we discuss how future PTAs can further improve these constraints, and detail the unique properties of these signals relative to the previously studied effects of ULDM on PTAs.

Direct detection experiments for light (sub-GeV) dark matter are making enormous leaps in reaching previously unexplored theory space. The need for accurate characterizations of target responses has led to a growing interplay between particle and condensed matter physics. This white paper summarizes recent progress on direct detection calculations that utilize state-of-the-art numerical tools in condensed matter physics and effective field theory techniques. These new results provide the theoretical framework for interpreting ongoing and planned experiments using electronic and collective excitations, and for optimizing future searches.

Semiconductors with $\mathcal{O}(\mathrm{meV})$ band gaps have been shown to be promising targets to search for sub-MeV mass dark matter (DM). In this paper we focus on a class of materials where such narrow band gaps arise naturally as a consequence of spin-orbit coupling (SOC). Specifically, we are interested in computing DM-electron scattering and absorption rates in these materials using state-of-the-art density functional theory techniques. To do this, we extend the DM interaction rate calculation to include SOC effects which necessitates a generalization to spin-dependent wave functions. We apply our new formalism to calculate limits for several DM benchmark models using an example ${\text{ZrTe}}_5$ target and show that the inclusion of SOC can substantially alter projected constraints.

We revisit the calculation of bosonic dark matter absorption via electronic excitations. Working in an effective field theory framework and consistently taking into account in-medium effects, we clarify the relation between dark matter and photon absorption. As is well-known, for vector (dark photon) and pseudoscalar (axion-like particle) dark matter, the absorption rates can be simply related to the target material’s optical properties. However, this is not the case for scalar dark matter, where the dominant contribution comes from a different operator than the one contributing to photon absorption, which is formally next-to-leading-order and does not suffer from in-medium screening. It is therefore imperative to have reliable first-principles numerical calculations and/or semi-analytic modeling in order to predict the detection rate. We present updated sensitivity projections for semiconductor crystal and superconductor targets for ongoing and proposed direct detection experiments.

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5-year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black hole mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

Models of Dark Matter (DM) can leave unique imprints on the Universe’s small scale structure by boosting density perturbations on small scales. We study the capability of Pulsar Timing Arrays to search for, and constrain, subhalos from such models. The models of DM we consider are ordinary adiabatic perturbations in ΛCDM, QCD axion miniclusters, models with early matter domination, and vector DM produced during inflation. We show that ΛCDM, largely due to tidal stripping effects in the Milky Way, is out of reach for PTAs. Axion miniclusters may be within reach, although this depends crucially on whether the axion relic density is dominated by the misalignment or string contribution. Models where there is matter domination with a reheat temperature below 1 GeV may be observed with future PTAs. Lastly, vector DM produced during inflation can be detected if it is lighter than 10${}^{\text{−}16}$−16 GeV. We also make publicly available a Python Monte Carlo tool for generating the PTA time delay signal from any model of DM substructure.

Collective excitations in condensed matter systems, such as phonons and magnons, have recently been proposed as novel detection channels for light dark matter. We show that excitation of (i) optical phonon polaritons in polar materials in an O(1  T) magnetic field (via the axion-photon coupling), and (ii) gapped magnons in magnetically ordered materials (via the axion wind coupling to the electron spin), can cover the difficult-to-reach O(1–100)  meV mass window of QCD axion dark matter with less than a kilogram-year exposure. Finding materials with a large number of optical phonon or magnon modes that can couple to the axion field is crucial, suggesting a program to search for a range of materials with different resonant energies and excitation selection rules; we outline the rules and discuss a few candidate targets, leaving a more exhaustive search for future work. Ongoing development of single photon, phonon, and magnon detectors will provide the key for experimentally realizing the ideas presented here.

Dirac materials, because of their small O(meV) band gap, are a promising target for dark photon-mediated scattering and absorption of light dark matter. In this paper, we characterize the daily modulation rate of dark matter interacting with a Dirac material due to anisotropies in their crystal structure. We show that daily modulation is an O(1) fraction of the total rate for dark matter scattering in the Dirac material ZrTe5. When present, the modulation is dominated by the orientation of the material’s dielectric tensor with respect to the dark matter wind and is maximized when the crystal is oriented such that the dark matter wind is completely aligned with the largest and smallest components of the dielectric tensor at two different times of the day. Because of the large modulation, any putative dark matter scattering signal could be rapidly verified or ruled out by changing the orientation of the crystal with respect to the wind and observing how the daily modulation pattern changes.

Searching for new physics in large data sets needs a balance between two competing effects—signal identification vs background distortion. In this work, we perform a systematic study of both single variable and multivariate jet tagging methods that aim for this balance. The methods preserve the shape of the background distribution by either augmenting the training procedure or the data itself. Multiple quantitative metrics to compare the methods are considered, for tagging 2-, 3-, or 4-prong jets from the QCD background. This is the first study to show that the data augmentation techniques of Planing and PCA based scaling deliver similar performance as the augmented training techniques of Adversarial NN and uBoost, but are both easier to implement and computationally cheaper.

The relic cosmological abundance of stable or long-lived, charge neutral, colored particles gets reduced by up to 4 orders of magnitude by annihilations that occur after QCD confinement. We compute the abundance and the cosmological bounds on relic gluinos. The same postconfinement effect strongly enhances coannihilations with a lighter dark matter particle, provided that their mass difference is below a few giga-electron volts. Charged colored particles (such as stops) can instead form baryons, which can be (quasi)stable in some models.

We introduce the gluequark Dark Matter candidate, an accidentally stable bound state made of adjoint fermions and gluons from a new confining gauge force. Such scenario displays an unusual cosmological history where perturbative freeze-out is followed by a non-perturbative re-annihilation period with possible entropy injection. When the gluequark has electroweak quantum numbers, the critical density is obtained for masses as large as PeV. Independently of its mass, the size of the gluequark is determined by the confinement scale of the theory, leading at low energies to annihilation rates and elastic cross sections which are large for particle physics standards and potentially observable in indirect detection experiments.

We introduce the gluequark Dark Matter candidate, an accidentally stable bound state made of adjoint fermions and gluons from a new confining gauge force. Such scenario displays an unusual cosmological history where perturbative freeze-out is followed by a non-perturbative re-annihilation period with possible entropy injection. When the gluequark has electroweak quantum numbers, the critical density is obtained for masses as large as PeV. Independently of its mass, the size of the gluequark is determined by the confinement scale of the theory, leading at low energies to annihilation rates and elastic cross sections which are large for particle physics standards and potentially observable in indirect detection experiments.

We explore the possibility that dark matter (DM) is the lightest hadron made of two stable color octet Dirac fermions Q. The cosmological DM abundance is reproduced for MQ≈12.5  TeV, compatibly with direct searches (the Rayleigh cross section, suppressed by 1/MQ^6, is close to present bounds), indirect searches (enhanced by QQ+Q¯Q¯→QQ¯+QQ¯ recombination), and with collider searches (where Q manifests as tracks, pair produced via QCD). Hybrid hadrons, made of Q and of standard model quarks and gluons, have large QCD cross sections, and do not reach underground detectors. Their cosmological abundance is 105 times smaller than DM, such that their unusual signals seem compatible with bounds. Those in the Earth and stars sank to their centers; the Earth crust and meteorites later accumulate a secondary abundance, although their present abundance depends on nuclear and geological properties that we cannot compute from first principles.

Dark Matter might be an accidentally stable baryon of a new confining gauge interaction. We extend previous studies exploring the possibility that the DM is made of dark quarks heavier than the dark confinement scale. The resulting phenomenology contains new unusual elements: a two-stage DM cosmology (freeze-out followed by dark condensation), a large DM annihilation cross section through recombination of dark quarks (allowing to fit the positron excess). Light dark glue-balls are relatively long lived and give extra cosmological effects, DM itself can remain radioactive.

We present generic formulæ for computing how Sommerfeld corrections together with bound-state formation affects the thermal abundance of Dark Matter with non-abelian gauge interactions. We consider DM as a fermion 3plet (wino) or 5plet under SU(2)L. In the latter case bound states raise to 11.5 TeV the DM mass required to reproduce the cosmological DM abundance and give indirect detection signals such as (for this mass) a dominant γ-line around 70 GeV. Furthermore, we consider DM co-annihilating with a colored particle, such as a squark or a gluino, finding that bound state effects are especially relevant in the latter case.