Search

We study the implications of relaxing the requirement for ultralight axions to account for all dark matter in the Universe by examining mixed dark matter (MDM) cosmologies with axion fractions $f \leq 0.3$ within the fuzzy dark matter (FDM) window $10^{-25}$ eV $\lesssim m \lesssim 10^{-23}$ eV. Our simulations, using a new MDM gravity solver implemented in AxiREPO, capture wave dynamics across various scales with high accuracy down to redshifts $z\approx 1$. We identify halos with Rockstar using the CDM component and find good agreement of inferred halo mass functions (HMFs) and concentration-mass relations with theoretical models across redshifts $z=1-10$. This justifies our halo finder approach a posteriori as well as the assumptions underlying the MDM halo model AxionHMcode. Using the inferred axion halo mass - cold halo mass relation $M_{\text{a}}(M_{\text{c}})$ and calibrating a generalised smoothing parameter $\alpha$ to our MDM simulations, we present a new version of AxionHMcode. The code exhibits excellent agreement with simulations on scales $k< 20 \ h$ cMpc$^{-1}$ at redshifts $z=1-3.5$ for $f\leq 0.1$ around the fiducial axion mass $m = 10^{-24.5}$ eV $ = 3.16\times 10^{-25}$ eV, with maximum deviations remaining below 10%. For axion fractions $f\leq 0.3$, the model maintains accuracy with deviations under 20% at redshifts $z\approx 1$ and scales $k< 10 \ h$ cMpc$^{-1}$, though deviations can reach up to 30% for higher redshifts when $f=0.3$. Reducing the run-time for a single evaluation of AxionHMcode to below $1$ minute, these results highlight the potential of AxionHMcode to provide a robust framework for parameter sampling across MDM cosmologies in Bayesian constraint and forecast analyses., Comment: 20 pages, 10 figures, 5 Tables, comments welcome

Recent studies have found evidence for parity violation in the BOSS spectroscopic galaxy survey, with statistical significance as high as $7\sigma$. These analyses assess the significance of the parity-odd four-point correlation function (4PCF) with a statistic called $\chi^2$. This statistic is biased if the parity-even eight-point correlation function (8PCF) of the data differs from the mock catalogs. We construct new statistics $\chi^2_\times$, $\chi^2_{\mathrm{null}}$ that separate the parity violation signal from the 8PCF bias term, allowing them to be jointly constrained. Applying these statistics to BOSS, we find that the parity violation signal ranges from $0$ to $2.5\sigma$ depending on analysis choices, whereas the 8PCF bias term is $\sim 6\sigma$. We conclude that there is no compelling evidence for parity violation in BOSS. Our new statistics can be used to search for parity violation in future surveys, such as DESI, without 8PCF biases., Comment: 46 pages, 11 figures; comments welcome!

Dark matter as scalar particles consisting of multiple species is well motivated in string theory where axion fields are ubiquitous. A two-field fuzzy dark matter (FDM) model features two species of ultralight axion particles with different masses, $m_1 \neq m_2$, which is extended from the standard one-field model with $m_a \sim 10^{-22}\,{\rm eV}$. Here we perform numerical simulations to explore the properties of two-field FDM haloes. We find that the central soliton has a nested structure when $m_2 \gg m_1$, which is distinguishable from the generic flat-core soliton in one-field haloes. However, the formation of this nested soliton is subject to many factors, including the density fraction and mass ratio of the two fields. Finally, we study non-linear structure formation in two-field cosmological simulations with self-consistent initial conditions and find that the small-scale structure in two-field cosmology is also distinct from the one-field model in terms of DM halo counts and soliton formation time., Comment: 11 pages, 5 figures. Published version

Using a single gravitational lens system observed at $\lesssim5$ milli-arcsecond resolution with very long baseline interferometry (VLBI), we place a lower bound on the mass of the fuzzy dark matter (FDM) particle, ruling out $m_\chi \leq 4.4\times10^{-21}~\mathrm{eV}$ with a 20:1 posterior odds ratio relative to a smooth lens model. We generalize our result to non-scalar and multiple-field models, such as vector FDM, with $m_{\chi,\mathrm{vec}} > 1.4 \times 10^{-21}~\mathrm{eV}$. Due to the extended source structure and high angular resolution of the observation, our analysis is directly sensitive to the presence of granule structures in the main dark matter halo of the lens, which is the most generic prediction of FDM theories. A model based on well-understood physics of ultra-light dark matter fields in a gravitational potential well makes our result robust to a wide range of assumed dark matter fractions and velocity dispersions in the lens galaxy. Our result is competitive with other lower bounds on $m_\chi$ from past analyses, which rely on intermediate modelling of structure formation and/or baryonic effects. Higher resolution observations taken at 10 to 100 GHz could improve our constraints by up to 2 orders of magnitude in the future., Comment: 5 pages, 2 figures. Accepted in MNRAS Letters

Fuzzy dark matter (FDM) is a dark matter candidate consisting of ultra-light scalar particles with masses around $10^{-22} \mathrm{eV}/c^2$, a regime where cold bosonic matter behaves as a collective wave rather than individual particles. It has increasingly attracted attention due to its rich phenomenology on astrophysical scales, with implications for the small-scale tensions present within the standard cosmological model, $\Lambda$CDM. Although constraints on FDM are accumulating in many different contexts, very few have been verified by self-consistent numerical simulations. We present new large numerical simulations of cosmic structure formation with FDM, solving the full Schr\"odinger-Poisson (SP) equations using the AxiREPO code, which implements a pseudo-spectral numerical method. Combined with our previous simulations, they allow us to draw a four-way comparison of matter clustering, contrasting results (such as power spectra) for each combination of initial conditions (FDM vs. CDM) and dynamics (SP vs. $N$-body). By disentangling the impact of initial conditions and non-linear dynamics, we can gauge the validity of approximate methods used in previous works, such as ordinary $N$-body simulations with an FDM initial power spectrum. Due to the comparatively large volume achieved in our FDM simulations, we are able to measure the FDM halo mass function from full wave simulations for the first time, and compare to previous results obtained using analytic or approximate approaches. We find that, due to the cut-off of small-scale power in the FDM power spectrum, haloes are linked via continuous, smooth, and dense filaments throughout the entire simulation volume (unlike for the standard $\Lambda$CDM power spectrum), posing significant challenges for reliably identifying haloes. We also investigate the density profiles of these filaments and compare to their CDM counterparts., Comment: 18 pages, 13 figures; submitted to MNRAS

Consider the symmetric group $S_n$ acting as a reflection group on the polynomial ring $k[x_1, \ldots, x_n]$, where $k$ is a field such that Char$(k)$ does not divide $n!$. We use Higher Specht polynomials to construct matrix factorizations of the discriminant of this group action: these matrix factorizations are indexed by partitions of $n$ and respect the decomposition of the coinvariant algebra into isotypical components. The maximal Cohen-Macaulay modules associated to these matrix factorizations give rise to a noncommutative resolution of the discriminant and they correspond to the nontrivial irreducible representations of $S_n$. All our constructions are implemented in Macaulay2 and we provide several examples. We also discuss extensions of these results to Young subgroups of $S_n$., Comment: 23 pages, comments welcome!

In the fuzzy dark matter (FDM) model, gravitationally collapsed objects always consist of a solitonic core located within a virialised halo. Although various numerical simulations have confirmed that the collapsed structure can be described by a cored NFW like density profile, there is still disagreement about the relation between the core mass and the halo mass. To fully understand this relation, we have assembled a large sample of cored haloes based on both idealised soliton mergers and cosmological simulations with various box sizes. We find that there exists a sizeable dispersion in the core-halo mass relation that increases with halo mass, indicating that the FDM model allows cores and haloes to coexist in diverse configurations. We provide a new empirical equation for a core halo mass relation with uncertainties that can encompass all previously found relations in the dispersion, and emphasise that any observational constraints on the particle mass using a tight one-to-one core-halo mass relation should suffer from an additional uncertainty on the order of 50 % for halo masses $ \ge 10^9 (8 \times 10^{-23} eV/ (mc^2))^{3/2} M_\odot$. We suggest that tidal stripping may be one of the effects contributing to the scatter in the relation.

We show that for the family of complex reflection groups $G=G(m,p,2)$ appearing in the Shephard--Todd classification, the endomorphism ring of the reduced hyperplane arrangement $A(G)$ is a non-commutative resolution for the coordinate ring of the discriminant $\Delta$ of $G$. This furthers the work of Buchweitz, Faber and Ingalls who showed that this result holds for any true reflection group. In particular, we construct a matrix factorization for $\Delta$ from $A(G)$ and decompose it using data from the irreducible representations of $G$. For $G(m,p,2)$ we give a full decomposition of this matrix factorization, including for each irreducible representation a corresponding a maximal Cohen--Macaulay module. The decomposition concludes that the endomorphism ring of the reduced hyperplane arrangement $A(G)$ will be a non-commutative resolution., Comment: 31 pages

An ultra-light bosonic particle of mass around $10^{-22}\,\mathrm{eV}/c^2$ is of special interest as a dark matter candidate, as it both has particle physics motivations, and may give rise to notable differences in the structures on highly non-linear scales due to the manifestation of quantum-physical wave effects on macroscopic scales, which could address a number of contentious small-scale tensions in the standard cosmological model, $\Lambda$CDM. Using a spectral technique, we here discuss simulations of such fuzzy dark matter (FDM), including the full non-linear wave dynamics, with a comparatively large dynamic range and for larger box sizes than considered previously. While the impact of suppressed small-scale power in the initial conditions associated with FDM has been studied before, the characteristic FDM dynamics are often neglected; in our simulations, we instead show the impact of the full non-linear dynamics on physical observables. We focus on the evolution of the matter power spectrum, give first results for the FDM halo mass function directly based on full FDM simulations, and discuss the computational challenges associated with the FDM equations. FDM shows a pronounced suppression of power on small scales relative to cold dark matter (CDM), which can be understood as a damping effect due to 'quantum pressure'. In certain regimes, however, the FDM power can exceed that of CDM, which may be interpreted as a reflection of order-unity density fluctuations occurring in FDM. In the halo mass function, FDM shows a significant abundance reduction below a characteristic mass scale only. This could in principle alleviate the need to invoke very strong feedback processes in small galaxies to reconcile $\Lambda$CDM with the observed galaxy luminosity function, but detailed studies that also include baryons will be needed to ultimately judge the viability of FDM., Comment: 16 pages, 14 figures; matches published version except minor corrections to numbers in fig. 1 and typos

minimal-lagrangians is a Python program which allows one to specify the field content of an extension of the Standard Model of particle physics and, using this information, to generate the most general renormalizable Lagrangian that describes such a model. As the program was originally created for the study of minimal dark matter models with radiative neutrino masses, it can handle additional scalar or Weyl fermion fields which are $\mathrm{SU}(3)_{\mathrm{C}}$ singlets, $\mathrm{SU}(2)_{\mathrm{L}}$ singlets, doublets or triplets, and can have arbitrary $\mathrm{U}(1)_{\mathrm{Y}}$ hypercharge. It is also possible to enforce an arbitrary number of global $\mathrm{U}(1)$ symmetries (with $\mathbb{Z}_2$ as a special case) so that the new fields can additionally carry such global charges. In addition to human-readable and $\mathrm{\LaTeX}$ output, the program can generate SARAH model files containing the computed Lagrangian, as well as information about the fields after electroweak symmetry breaking (EWSB), such as vacuum expectation values (VEVs) and mixing matrices. This capability allows further detailed investigation of the model in question, with minimal-lagrangians as the first component in a tool chain for rapid phenomenological studies of "minimal" dark matter models requiring little effort and no unnecessary input from the user., Comment: 40 pages, 1 figure; version accepted by and published in CPC; code is available at https://gitlab.com/Socob/minimal-lagrangians

In this master's thesis, the minimal dark matter models with radiative neutrino masses introduced in arXiv:1308.3655 are examined both in general and, using the model T1-3-B ($\alpha = 0$) as a representative, in more detail. In the process, it is both shown how such models are built and a numerical analysis tool chain is established, automating as many steps as possible from the very beginning (deriving the Lagrangian from a given field content) to the very end (calculating observables such as the dark matter relic density, the direct detection cross section and the neutrino masses and mixing matrices). As an important result in laying out the mathematical groundwork, a convention for representations of SU(2) is established which, in this form, is not found in previous works. Another major result is the development of the program minimal-lagrangians, whose capabilities are described in the thesis. It is of great use both for the construction of Lagrangians for new models as well their verification or quick surveys. minimal-lagrangians can be applied to any model within the class it is designed to handle, extending its utility far beyond a one-time result that is only obtained for a single model. Finally, the first steps in the analysis of the model T1-3-B ($\alpha = 0$) are taken. It is shown that the model can reproduce both the singlet-doublet fermion dark matter model and the inert triplet model in the appropriate limits. Moreover, the neutrino mass matrix is found and it is demonstrated that the model can accommodate current constraints on the neutrino masses. Lastly, the behavior of the dark matter relic density as it depends on the parameter responsible for the generation of neutrino masses is investigated. Further in-depth study of the model T1-3-B ($\alpha = 0$) was performed in a separate publication (arXiv:1812.11133)., Comment: 226 pages, 32 figures; master's thesis at Westf\"alische Wilhelms-Universit\"at M\"unster

In these proceedings, we present a study of a combined singlet--doublet fermion and triplet scalar model for dark matter (DM). Together, these models form a simple extension of the Standard Model (SM) that can account for DM and explain the existence of neutrino masses, which are generated radiatively. However, this also implies the existence of lepton flavour violating (LFV) processes. In addition, this particular model allows for gauge coupling unification. The new fields are odd under a new $\mathbb{Z}_2$ symmetry to stabilise the DM candidate. We analyse the DM, neutrino mass and LFV aspects, exploring the viable parameter space of the model. This is done using a numerical random scan imposing successively the neutrino mass and mixing, relic density, Higgs mass, direct detection, collider and LFV constraints. We find that DM in this model is fermionic for masses below about 1 TeV and scalar above. We observe a high degree of complementarity between direct detection and LFV experiments, which should soon allow to fully probe the fermionic DM sector and at least partially the scalar DM sector., Comment: 4 pages, 1 figure; contribution to the 2019 EW session of the 54th Rencontres de Moriond (summary of arXiv:1812.11133)

We present a detailed study of a combined singlet-doublet fermion and triplet scalar model for dark matter. These models have only been studied separately in the past. Together, they form a simple extension of the Standard Model that can account for dark matter and explain the existence of neutrino masses, which are generated radiatively. This holds even if singlet-doublet fermions and triplet scalars never contribute simultaneously to the dark matter abundance. However, this also implies the existence of lepton flavour violating processes. In addition, this particular model allows for gauge coupling unification. The new fields are odd under a new $\mathbb{Z}_2$ symmetry to stabilise the dark matter candidate. We analyse the dark matter, neutrino mass and lepton flavour violation aspects both separately and in conjunction, exploring the viable parameter space of the model. This is done using a numerical random scan imposing successively the neutrino mass and mixing, relic density, Higgs mass, direct detection, collider and lepton flavour violation constraints. We find that dark matter in this model is fermionic for masses below about 1 TeV and scalar above. The narrow mass regions found previously for the two separate models are enlarged by their coupling. While coannihilations of the weak isospin partners are sizeable, this is not the case for fermions and scalars despite their often similar masses due to the relatively small coupling of the two sectors, imposed by the small neutrino masses. We observe a high degree of complementarity between direct detection and lepton flavour violation experiments, which should soon allow to fully probe the fermionic dark matter sector and at least partially the scalar dark matter sector., Comment: 24 pages, 12 figures; version accepted by and published in JHEP