Search

The intracluster medium (ICM) is the low-density diffuse magnetized plasma in galaxy clusters, which reaches virial temperatures of up to 10^8 K. Under these conditions, the plasma is weakly collisional and therefore has an anisotropic pressure tensor with respect to the local direction of the magnetic field. This triggers very fast, Larmor-scale, pressure-anisotropy-driven kinetic instabilities that alter magnetic field amplification. We study magnetic field amplification through a turbulent small-scale dynamo, including the effects of the kinetic instabilities, during the evolution of a typical massive galaxy cluster. A specific aim of this work is to establish a redshift limit from which a dynamo has to start to amplify the magnetic field up to equipartition with the turbulent velocity field at redshift z=0. We implemented 1D radial profiles for various plasma quantities for merger trees generated with the Modified GALFORM algorithm. We assume that turbulence is driven by successive mergers of dark matter halos and construct effective models for the Reynolds number Re_eff dependence on the magnetic field in three different magnetization regimes, including the effects of kinetic instabilities. The magnetic field growth rate is calculated for the different Re_eff models. The model results in a higher magnetic field growth rate at higher redshift. For all scenarios considered, to reach equipartition at z=0, the amplification of the magnetic field has to start at redshift z_start=1.5 and above. The time to reach equipartition can be significantly shorter, in cases with systematically smaller turbulent forcing scales, and for the highest Re_eff models. Merger trees are useful tools for studying the evolution of magnetic fields in weakly collisional plasmas, and could also be used to constrain the different stages of the dynamo that potentially could be observed by future radio telescopes., 15 pages, 9 figures

The stochastic signal detected by pulsar timing arrays (PTAs) has raised great interest in understanding its physical origin. Assuming the signal is a cosmological gravitational-wave background produced by overly large primordial density perturbations, we investigate the sound speed resonance effect with an oscillatory behavior using the combined PTA data from NANOGrav 15-yr data set, PPTA DR3, and EPTA DR2. We find that the stochastic signal can be explained by the induced gravitational waves sourced by the sound speed resonance mechanism, with the oscillation frequency $f_* \in [1.52, 4.67] \times 10^{-7}$Hz and the start time of oscillation $|\tau_0| \in [2.17, 95.7] \times 10^7$s., Comment: 16 pages, 2 figures, 1 table

We show that, whenever the perturbations of some field are excited during inflation by a physical process on sub-horizon scales, they unavoidably generate, even through gravitational interactions alone, a significant resonant IR cascade of power down to scales that are of the order of the horizon at that time (we denote these scales as near IR). We provide general analytic one-loop results for the enhancement of the IR power of the curvature perturbation generated by this effect, highlighting the role played by the resonance. We then study a number of examples in which the excited state is: (i) an isocurvature field, (ii) the curvature perturbation itself, (iii) a mixture of curvature and isocurvature fluctuations driven to an excited state by their coupled dynamics. In the cases shown, the cascade significantly modifies the near IR part of the power spectrum of the curvature perturbation with respect to the linear theory, indicating that this effect can impact the phenomenology associated with a variety of mechanisms considered in the literature, notably concerning primordial black holes and gravitational waves., 37 pages, 5 figures

This thesis treats the topic of CMB Spectral Distortions (SDs), which represent any deviation from a pure black body shape of the CMB energy spectrum. As such, they can be used to probe the inflationary, expansion and thermal evolution of the universe both within $\Lambda$CDM and beyond it. The currently missing observation of this rich probe of the universe makes of it an ideal target for future observational campaigns. In fact, while the $\Lambda$CDM signal guarantees a discovery, the sensitivity to a wide variety of new physics opens the door to an enormous uncharted territory. In light of these considerations, the thesis opens by reviewing the topic of CMB SDs in a pedagogical and illustrative fashion, aimed at waking the interest of the broader community. This introductory premise sets the stage for the first main contribution of the thesis to the field of SDs: their implementation in the Boltzmann solver CLASS and the parameter inference code MontePython. The CLASS+MontePython pipeline is publicly available, fast, it includes all sources of SDs within $\Lambda$CDM and many others beyond that, and allows to consistently account for any observational setup. By means of these numerical tools, the second main contribution of the thesis consists in showcasing the versatility and competitiveness of SDs for several cosmological models as well as for a number of different mission designs. Among others, the results cover features in the primordial power spectrum, primordial gravitational waves, non-standard dark matter properties, primordial black holes, primordial magnetic fields and Hubble tension. Finally, the manuscript is disseminated with (20) follow-up ideas that naturally extend the work carried out so far, highlighting how rich of unexplored possibilities the field of CMB SDs still is. The hope is that these suggestions will become a propeller for further interesting developments., Comment: PhD thesis. Pedagogical review of theory, experimental status and numerical tools (CLASS+MontePython) with broad overview of applications. Includes 20 original follow-up ideas

Invertible disformal transformations serve as a useful tool to explore ghost-free scalar-tensor theories. In this paper, we construct a generalization of invertible disformal transformations that involves arbitrary higher-order covariant derivatives of the scalar field. As a result, we obtain a more general class of ghost-free scalar-tensor theories than ever. Notably, our generalization is such that matter fields can be consistently coupled to these theories without introducing an unwanted extra degree of freedom in the unitary gauge., 9 pages

With the emersion of precise cosmology and the emergence of cosmic tensions, we are faced with the question of whether the simple model of cold dark matter needs to be extended and whether doing so can alleviate the tensions and improve our understanding of the properties of dark matter. In this study, we investigate one of the generalized models of dark matter so that the behavior of this dark matter changes according to the scale of $k$. In large scales (small $k$'s), the dark matter is cold, while it becomes warm for small scales (large $k$'s). This behavior is modeled phenomenologically for two different scenarios. We show that the $S_8$ tension can be alleviated, but the $H_0$ tension becomes milder while not too much., 8 pages, 6 figures. comments are welcome

Recently it has been shown that the cosmological dynamics of covariant $f(Q)$ gravity depend on different affine connections. In this paper, two specific $f(Q)$ models are investigated with SNe+CC+BAO+QSO observational data, and the spatial curvature of the universe is studied in covariant $f(Q)$ gravity. It is found that the parameter $\mathcal{X}$ characterizing affine connections significantly affects the behavior of the effective equation of state $w_Q$ and may drive it across the phantom divide line. A flat universe could be favored in covariant $f(Q)$ gravity rather than the $\Lambda$CDM model. These results imply some inertial effects of the universe change the cosmic dynamics and renew our cognition of the geometric structure of the universe., Comment: 14 pages, 6 figures, 3 tables

We study a new source of stochastic gravitational wave background (SGWB) from the final collapse of a string-wall network. In the context of $N_{\rm DW}=1$ axionic string-wall network, the final collapse of walls bounded by strings can release gravitational waves (GWs). This source is typically considered negligible due to its subdominance compared to GW emissions throughout the long-term evolution in the scaling regime. However, in some cases, a network can be driven outside of horizon by inflation and later re-enter horizon. Then, the network's final collapse after re-entering horzion becomes the dominant GW source and therefore cannot be neglected. Our caculation of the corresponding GW spectrum suggests it could potentially explain the nano-Hertz SGWB signal possibly detected by various Pulsar Timing Array experiments. In addition, with different parameter choices, the resultant GWs could be probed by various GW interferometry experiments., 5 pages, 3 figures

In addition to the Standard Model, the introduction of a singlet complex scalar field that acquires vacuum expectation value may give rise to a cosmologically stable pseudo-Nambu-Goldstone boson (pNGB), a suitable dark matter (DM) candidate. This work extends this scenario by including a second cosmological stable particle: a fermion singlet. The pNGB and the new fermion can be regarded as DM candidates simultaneously, both interacting with the Standard Model through Higgs portals via two non-degenerate Higgs bosons. We explore the thermal freeze-out of this scenario, with particular emphasis on the increasing yield of the pNGB before it completely decouples (recently called \textit{Bouncing DM}). We test the model under collider, relic abundance, and direct detection, and we explore the consequences of the yield bouncing on indirect detection observables today.

Dark Matter (DM) particles with sufficiently large cross sections may scatter as they travel through Earth's bulk. The corresponding changes in the DM flux give rise to a characteristic daily modulation signal in detectors sensitive to DM-electron interactions. Here, we report results obtained from the first underground operation of the DAMIC-M prototype detector searching for such a signal from DM with MeV-scale mass. A model-independent analysis finds no modulation in the rate of 1$e^-$ events with periods in the range 1-48 h. We then use these data to place exclusion limits on DM in the mass range [0.53, 2.7] MeV/c$^2$ interacting with electrons via a dark photon mediator. Taking advantage of the time-dependent signal we improve by $\sim$2 orders of magnitude on our previous limit obtained from the total rate of 1$e^-$ events, using the same data set. This daily modulation search represents the current strongest limit on DM-electron scattering via ultralight mediators for DM masses around 1 MeV/c$^2$.

We present the possibility that the seesaw mechanism and nonthermal leptogenesis can be probed via primordial non-Gaussianities in the context of a majoron curvaton model. Originating as a massless Nambu-Goldstone boson from the spontaneous breaking of the global baryon ($B$) minus lepton ($L$) number symmetry at a scale $v_{B-L}$, majoron becomes massive when it couples to a new confining sector through anomaly. Acting as a curvaton, majoron produces the observed red-tilted curvature power spectrum without relying on any inflaton contribution, and its decay in the post-inflationary era gives rise to a nonthermal population of right-handed neutrinos that participate in leptogenesis. A distinctive feature of the mechanism is the generation of observable non-Gaussianity, nontrivially linked to the scale of seesaw and leptogenesis. Identifying the viable parameter space, we show that the non-Gaussianity parameter $f_{\rm NL} \gtrsim \mathcal{O} (0.1)$ is produced for high-scale seesaw ($v_{B-L}$ at $\mathcal{O}(10^{14-17})$ GeV) and leptogenesis ($M_1 \gtrsim \mathcal{O}(10^6)$ GeV) where the latter represents the lightest right-handed neutrino mass. While the current bounds on local non-Gaussianity excludes some part of parameter space, the rest can be fully probed by future experiments like CMB-S4, LSST, and 21 cm tomography., 19 pages + references, 7 figures

We present a comprehensive study of Parker-type bounds on magnetic monopoles with arbitrary magnetic charge, including minicharged monopoles and magnetic black holes. We derive the bounds based on the survival of galactic magnetic fields, seed magnetic fields, as well as primordial magnetic fields. We find that monopoles with different magnetic charges are best constrained by different astrophysical systems: while monopoles with a Dirac charge are tightly constrained by seed galactic magnetic fields, minicharged monopoles are strongly constrained by primordial magnetic fields, and magnetic black holes by the density of dark matter. We also assess the viability of the various types of monopoles as dark matter, by studying whether they can cluster with galaxies hosting magnetic fields., 33 pages, 6 figures

We obtain a new sample of 1192 Type I quasars with the UV-optical, radio and X-ray wavebands coverage by combining \citet{Huang2022} and other matching data of SDSS-DR16 with FIRST, XMM-Newton, and Chandra Source Catalog, and a sample of 407 flat-spectrum radio-loud quasars (FSRLQs) of blazars from the Roma-BZCAT, which can be used to investigate their multi-band luminosity correlations and measure the luminosity distances of these Type I radio-loud quasars (RLQs) samples. We check the correlation between X-ray, UV-optical, and radio luminosity for various groupings of radio-quiet quasars (RQQs) and RLQs by parameterizing X-ray luminosity as a sole function of UV-optical or radio luminosity and as a joint function of UV-optical radio luminosity, which also can be employed to determine these cosmological distances. By Bayesian information criterion (BIC), the data suggest that the X-ray luminosity of RQQs is indirectly correlative with radio luminosity because of the connection between UV-optical and radio luminosity. But for RLQs, the X-Ray luminosity is directly related to radio luminosity, and the correlations between X-ray, optical/UV, and radio luminosity increase with the ratio of monochromatic luminosities logR. Meanwhile, we compare the results from RLQs with different UV-optical power law index ${\Gamma _{UV}}$, the goodness of fit for RLQs with ${\Gamma _{UV}}\le 1.6$ seems to be better. Finally, we apply a combination of Type I RLQs and SN Ia Pantheon to verify the nature of dark energy concerning whether or not its density deviates from the constant, and give the statistical results., Comment: arXiv admin note: text overlap with arXiv:2207.08390

Dark matter induced event rate in an Earth-based detector is predicted to show an annual modulation as a result of the Earth's orbital motion around the Sun. We searched for this modulation signature using the ionization signal of the DarkSide-50 liquid argon time projection chamber. No significant signature compatible with dark matter is observed in the electron recoil equivalent energy range above $40~{\rm eV_{ee}}$, the lowest threshold ever achieved in such a search., 8 pages, 4 figures

Majoron-like bosons would emerge from a supernova (SN) core by neutrino coalescence of the form $\nu\nu\to\phi$ and $\bar\nu\bar\nu\to\phi$ with 100 MeV-range energies. Subsequent decays to (anti)neutrinos of all flavors provide a flux component with energies much larger than the usual flux from the "neutrino sphere." The absence of 100 MeV-range events in the Kamiokande-II and Irvine-Michigan-Brookhaven signal of SN 1987A implies that less than 1% of the total energy was thus emitted and provides the strongest constraint on the Majoron-neutrino coupling of $g\lesssim 10^{-9}\,{\rm MeV}/m_\phi$ for $100~{\rm eV}\lesssim m_\phi\lesssim 100~{\rm MeV}$. It is straightforward to extend our new argument to other hypothetical feebly interacting particles., Comment: 5+8 pages, 2+6 figures. Matches published version

In this paper, we propose a novel mechanism in $\mathbb{T}^2$-inflation to enhance the power spectrum large enough to seed primordial black holes (PBHs) formation. To accomplish this, we consider the coupling function between the inflaton field and $\mathbb{T}^2= T_{\mu \nu}T^{\mu \nu}$ term. PBHs formed within this scenario can contribute partially or entirely to dark matter (DM) abundance. Furthermore, the amplification in the scalar power spectrum will concurrently produce significant scalar-induced gravitational waves (SIGWs) as a second-order effect. In addition, the energy spectrum associated with SIGWs can be compatible with the recent NANOGrav 15-year stochastic gravitational wave detection and fall into the sensitivity range of other forthcoming GW observatories., Comment: 8 pages, 7 figures, and one table

We provide a pedagogical approach to gravitational waves in cosmology with focus on gravitational wave signals related to primordial black holes. These lectures notes contain more details than one is able to present in the two two-hour lectures they are meant to and, as such, they should be thought as a complementary material. The main aim of these lectures is that, by the end, one obtains a certain degree of intuition on gravitational waves in cosmology and understands the basic features of scalar induced gravitational waves. We also highlight must-check properties of induced gravitational waves as well as current issues regarding secondary gravitational waves in cosmology. Throughout the lecture we provide exercises, supplementary information and activities with public codes to be ready to derive your own results., Lectures notes prepared for the ICCUB School 2023 on Primordial Black holes in the University of Barcelona. Comments and corrections are welcome

The Square Kilometre Array (SKA) will have the sensitivity to take the 3D light-cones of the 21 cm signal from the epoch of reionization. This signal, however, is highly non-Gaussian and can not be fully interpreted by the traditional statistic using power spectrum. In this work, we introduce the 3D ScatterNet that combines the normalizing flows with solid harmonic wavelet scattering transform, a 3D CNN featurizer with inductive bias, to perform implicit likelihood inference (ILI) from 21 cm light-cones. We show that 3D ScatterNet outperforms the ILI with a fine-tuned 3D CNN in the literature. It also reaches better performance than ILI with the power spectrum for varied light-cone effects and varied signal contaminations., 9 pages, 4 figures, 2 tables. Accepted to ICML 2023 Machine Learning for Astrophysics workshop. Comments and suggestions are welcome

The inner slope gammaDM of the dark matter (DM) density profile of cosmological halos carries information about the properties of DM and/or baryonic processes affecting the halo gravitational potential. Cold DM cosmological simulations predict steep inner slopes, gammaDM>~1. We test this prediction on the MACS J1206.2-0847 cluster at redshift z=0.44, whose DM density profile was claimed to be cored at the center. We determine the cluster DM density profile from 2 kpc from the cluster center to the virial radius (~2 Mpc), using the velocity distribution of ~500 cluster galaxies and the velocity dispersion profile of the Brightest Cluster Galaxy (BCG), obtained from VIMOS@VLT and MUSE@VLT data. We solve the Jeans equation of dynamical equilibrium using an upgraded version of the MAMPOSSt method. The total mass profile is modeled as a sum of a generalized-NFW profile that describes the DM component, allowing for a free inner slope of the density profile, a Jaffe profile that describes the BCG stellar mass component, and a non-parametric baryonic profile that describes the sum of the remaining galaxy stellar mass and of the hot intra-cluster gas mass. Our total mass profile is in remarkable agreement with independent determinations based on X-ray observations and strong lensing. We find gammaDM=0.7 -0.1 +0.2 (68% confidence levels), consistent with predictions from recent LambdaCDM cosmological numerical simulations., Submitted to ApJ on June, 1st 2023. 14 pages, 9 figures

Generative adversarial networks (GANs) are frequently utilized in astronomy to construct an emulator of numerical simulations. Nevertheless, training GANs can prove to be a precarious task, as they are prone to instability and often lead to mode collapse problems. Conversely, the diffusion model also has the ability to generate high-quality data without adversarial training. It has shown superiority over GANs with regard to several natural image datasets. In this study, we undertake a quantitative comparison between the denoising diffusion probabilistic model (DDPM) and StyleGAN2 (one of the most robust types of GANs) via a set of robust summary statistics from scattering transform. In particular, we utilize both models to generate the images of 21 cm brightness temperature mapping, as a case study, conditionally based on astrophysical parameters that govern the process of cosmic reionization. Using our new Fr\'echet Scattering Distance (FSD) as the evaluation metric to quantitatively compare the sample distribution between generative models and simulations, we demonstrate that DDPM outperforms StyleGAN2 on varied sizes of training sets. Through Fisher forecasts, we demonstrate that on our datasets, StyleGAN2 exhibits mode collapses in varied ways, while DDPM yields a more robust generation. We also explore the role of classifier-free guidance in DDPM and show the preference for a non-zero guidance scale only when the training data is limited. Our findings indicate that the diffusion model presents a promising alternative to GANs in the generation of accurate images. These images can subsequently provide reliable parameter constraints, particularly in the realm of astrophysics., Comment: 13 pages, 10 figures, 1 table. Submitted to MNRAS. Comments welcome

The Laser Interferometer Space Antenna (LISA) will detect gravitational waves (GWs) emitted by massive black hole binaries (MBHBs) in the low-frequency ($\sim$mHz) band. Low-mass lenses, such as dark-matter (DM) subhalos, have sizes comparable to the wavelength of these GWs. Encounters with these lenses produce wave-optics (WO) effects that alter waveform phase and amplitude. Thus, a single event with observable WO effects can be used to probe the lens properties. In this paper, we first compute the probability of observing WO effects in a model-agnostic way. We perform parameter estimation over approximately 1000 MBHBs with total mass, mass ratio, and redshift spanning the ranges relevant to LISA. We then calculate lensing rates using three semi-analytical models of MBHB populations. In both cases, we use a waveform model that includes merger, ringdown, and higher-order modes. We use two lens population models: the theory-based Press-Schechter halo mass function and an observation-based model derived from Sloan Digital Sky Survey, called the measured velocity function. We find that the probability of detecting WO effects can be as large as $\sim 3\%$, $\sim1.5\%$, and $\sim 1 \%$ at $1\sigma$, $3\sigma$, and $5\sigma$ confidence levels, respectively. The most optimistic MBHB population model yields $\sim 8$, $\sim 4$, and $\sim 3$ events at the same confidence levels, while the rates drop to $\sim 0.01$ in the more pessimistic scenarios. The most likely lens masses probed by LISA are in the range $(10^3, 10^8)\, M_{\odot}$, and the most probable redshifts are in the range $(0.3, 1.7)$. Therefore, LISA observations of WO effects can probe DM subhalos, complementing strong lensing and other observations., Comment: 24 pages, 15 figures. Comments are welcome

The density of relic neutrinos is expected to be enhanced due to clustering in our local neighbourhood at Earth. We introduce a novel analytical technique to calculate the neutrino overdensity, based on kinetic field theory. Kinetic field theory is a particle-based theory for cosmic structure formation and in this work we apply it for the first time to massive neutrinos. The gravitational interaction is expanded in a perturbation series and we take into account the first-order contribution to the local density of relic neutrinos. For neutrino masses that are consistent with cosmological neutrino mass bounds, our results are in excellent agreement with state-of-the-art calculations., Comment: 7 pages, 1 figure

Powerful radio sources associated with super-massive black holes are among the most luminous objects in the Universe, and are frequently recognized both as cosmological probes and active constituents in the evolution of galaxies. We present alignments between radio jets and cold molecular gas in the environment of distant radio galaxies, and show that the brightness of the radio synchrotron source can be enhanced by its interplay with the molecular gas. Our work is based on CO J>1 observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of three radio galaxies with redshifts in the range 1.4 < z < 2.1, namely MRC 0114-211 (z = 1.41), MRC 0156-252 (z = 2.02), and MRC 2048-272 (z = 2.05). These ALMA observations support previous work that found molecular gas out to 50 kpc in the circumgalactic environment, based on a CO(1-0) survey performed with the Australia Telescope Compact Array (ATCA). The CO emission is found along the radio axes but beyond the main radio lobes. When compared to a large sample of high-z radio galaxies from the literature, we find that the presence of this cold molecular medium correlates with an increased flux-density ratio of the main vs. counter lobe. This suggest that the radio lobe brightens when encountering cold molecular gas in the environment. While part of the molecular gas is likely related to the interstellar medium (ISM) from either the host or a companion galaxy, a significant fraction of the molecular gas in these systems shows very low excitation, with r$_{2-1/1-0}$ and r$_{3-2/1-0}$ values $\lesssim$0.2. This could be part of the circumgalactic medium (CGM)., Accepted for publication in ApJ (19 pages, 6 figures)

This paper presents an octree construction method, called Cornerstone, that facilitates global domain decomposition and interactions between particles in mesh-free numerical simulations. Our method is based on algorithms developed for 3D computer graphics, which we extend to distributed high performance computing (HPC) systems. Cornerstone yields global and locally essential octrees and is able to operate on all levels of tree hierarchies in parallel. The resulting octrees are suitable for supporting the computation of various kinds of short and long range interactions in N-body methods, such as Barnes-Hut and the Fast Multipole Method (FMM). While we provide a CPU implementation, Cornerstone may run entirely on GPUs. This results in significantly faster tree construction compared to execution on CPUs and serves as a powerful building block for the design of simulation codes that move beyond an offloading approach, where only numerically intensive tasks are dispatched to GPUs. With data residing exclusively in GPU memory, Cornerstone eliminates data movements between CPUs and GPUs. As an example, we employ Cornerstone to generate locally essential octrees for a Barnes-Hut treecode running on almost the full LUMI-G system with up to 8 trillion particles.

Within the extension of the $$\varLambda $$ Λ CDM model, allowing for the presence of neutrinos or warm dark matter, we develop the analytical cosmological perturbation theory. It covers all spatial scales where the weak gravitational field regime represents a valid approximation. Discrete particles – the sources of the inhomogeneous gravitational field – may be relativistic. Similarly to the previously investigated case of nonrelativistic matter, the Yukawa interaction range is naturally incorporated into the first-order scalar metric corrections.

We present a cross-correlation analysis between $1'$ resolution total intensity and polarization observations from the Atacama Cosmology Telescope (ACT) at 150 and 220 GHz and 15$''$ mid-infrared photometry from the Wide-field Infrared Survey Explorer (WISE) over 107 12.5$^\circ\times$12.5$^\circ$ patches of sky. We detect a spatially isotropic signal in the WISE$\times$ACT $TT$ cross power spectrum at 30$σ$ significance that we interpret as the correlation between the cosmic infrared background at ACT frequencies and polycyclic aromatic hydrocarbon (PAH) emission from galaxies in WISE, i.e., the cosmic PAH background. Within the Milky Way, the Galactic dust $TT$ spectra are generally well-described by power laws in $\ell$ over the range 10$^3 < \ell < $10$^4$, but there is evidence both for variability in the power law index and for non-power law behavior in some regions. We measure a positive correlation between WISE total intensity and ACT $E$-mode polarization at 1000$ < \ell \lesssim $6000 at $>$3$σ$ in each of 35 distinct $\sim$100 deg$^2$ regions of the sky, suggesting alignment between Galactic density structures and the local magnetic field persists to sub-parsec physical scales in these regions. The distribution of $TE$ amplitudes in this $\ell$ range across all 107 regions is biased to positive values, while there is no evidence for such a bias in the $TB$ spectra. This work constitutes the highest-$\ell$ measurements of the Galactic dust $TE$ spectrum to date and indicates that cross-correlation with high-resolution mid-infrared measurements of dust emission is a promising tool for constraining the spatial statistics of dust emission at millimeter wavelengths., 20 pages, 14 figures, submitted to ApJ

We study the reheating and leptogenesis in the case of a vector inflaton. We concentrate on particle production during the phase of oscillating background, especially gravitational production induced by the presence of non-minimal coupling imposed by an isotropic and homogeneous Universe. Including processes involving the exchange of graviton, we then extend our study to decay into fermions via direct or anomalous couplings. The necessity of non-minimal gravitational coupling and the gauge nature of couplings to fermions implies a much richer phenomenology than for a scalar inflaton., 24 pages, 4 figures, 13 tikzpicture Feynman diagrams, 2 tables

The astrophysical Stochastic Gravitational Wave Background (SGWB) originates from the mergers of compact binary objects in the observable Universe that are otherwise undetected as individual events, along with other sources such as supernovae, magnetars etc. The individual GW signal is time-varying over a time scale that depends on the chirp mass of the coalescing binaries. Another timescale that plays a role is the timescale at which the sources repeat, which depends on the merger rate. The combined effect of these two leads to a breakdown of the time-translation symmetry of the observed SGWB and leads to a correlation between different frequency modes in the signal covariance matrix of the SGWB. Using an ensemble of SGWB sky map due to binary black hole (BBH) coalescence, calculated using simulations of different astrophysical and primordial black hole mass distribution and merger rates, we show how the structure of the signal covariance matrix varies. This structure in the signal covariance matrix brings additional information about the sources on top of the power spectrum of the SGWB. We show that there is up to a factor of 2.5 improvement in the Figure of Merit (FoM) by using this additional information in comparison to only power spectrum estimation for the LIGO-Virgo-KAGRA (LVK) network of detectors with the design sensitivity noise with two years of observation time. The inclusion of the off-diagonal correlation in the covariance of the SGWB in the data analysis pipelines will be beneficial in the quest for the SGWB signal in LVK frequency bands as well as in lower frequencies (nano-Hz, mili-Hz, deci-Hz) and in getting an insight into the origin of the SGWB signal., 14 pages, 13 figures. Prepare for submission in MNRAS

The accurate modelling of the Point Spread Function (PSF) is of paramount importance in astronomical observations, as it allows for the correction of distortions and blurring caused by the telescope and atmosphere. PSF modelling is crucial for accurately measuring celestial objects' properties. The last decades brought us a steady increase in the power and complexity of astronomical telescopes and instruments. Upcoming galaxy surveys like Euclid and LSST will observe an unprecedented amount and quality of data. Modelling the PSF for these new facilities and surveys requires novel modelling techniques that can cope with the ever-tightening error requirements. The purpose of this review is three-fold. First, we introduce the optical background required for a more physically-motivated PSF modelling and propose an observational model that can be reused for future developments. Second, we provide an overview of the different physical contributors of the PSF, including the optic- and detector-level contributors and the atmosphere. We expect that the overview will help better understand the modelled effects. Third, we discuss the different methods for PSF modelling from the parametric and non-parametric families for ground- and space-based telescopes, with their advantages and limitations. Validation methods for PSF models are then addressed, with several metrics related to weak lensing studies discussed in detail. Finally, we explore current challenges and future directions in PSF modelling for astronomical telescopes., Comment: 63 pages, 14 figures. Submitted

With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries., White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html)

Time delays from strong gravitational lensing provide a one-step absolute distance measurement. Thus, they measure $H_0$ independently of all other probes. We first review the foundations and history of time-delay cosmography. Then, we illustrate the current state of the art by means of two recent case studies that have been real breakthroughs: i) the quadruply imaged quasar lensed by a galaxy-scale deflector RXJ1131$-$1231, for which spatially resolved stellar kinematics is available; ii) the multiply imaged supernova "Refsdal", the first with measured time delays, lensed by cluster MACS1149.5$+$2223. We conclude by discussing the exciting future prospects of time-delay cosmography in the coming decade., Invited chapter for the edited book "Hubble Constant Tension" (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024)

Self-interacting dark matter (SIDM) has been proposed as an alternative to the standard collisionless cold dark matter to explain the diversity of galactic rotation curves and core-cusp problems seen at small scales. Here, we estimate the constraints on SIDM for a sample of 11 relaxed galaxy groups with X-ray observations from Chandra and XMM-Newton. We fit the dark matter density distribution to the Einasto profile and use the estimated Einasto $\alpha$ parameter to constrain the SIDM cross-section, based on the empirical relation between the two, which was obtained in Eckert et al (2022). We obtain a non-zero central estimate for the cross-section per unit mass ($\sigma/m$) for seven groups, with the most precise estimate obtained for NGC 5044, given by $\sigma/m=0.165 \pm 0.025~\rm{cm^2/g}$, for dark matter velocity dispersion of about 300 km/sec. For the remaining four groups, we obtain 95% c.l. upper limits on $\sigma/m < 0.16-6.61~\rm{cm^2/g}$ with dark matter velocity dispersions between 200-500 km/sec, with the most stringent limit for our sample obtained for the group MKW 4, given by $\sigma/m< 0.16~\rm{cm^2/g}$ for dark matter velocity dispersion of about 350 km/sec., Comment: 14 pages, 22 figures. Accepted in Phys. Dark Universe

In this paper, we investigate an extension of standard $\Lambda$CDM model by allowing: a temporal evolution in the equation of state (EoS) of DM via Chevallier-Polarski-Linder parametrization, and the constant non-null sound speed. We also consider the properties of neutrinos, such as the effective neutrino mass and the effective number of neutrino species as free parameters. We derive the constraints on this scenario by using the data from the Planck-2018 cosmic microwave background (CMB), baryonic acoustic oscillations (BAO), the local value of the Hubble constant from the Hubble Space Telescope (HST), and some large scale structure (LSS) information from the abundance of galaxy clusters. We find constraints on the EoS and sound speed of DM very close to the null value in all cases thus no significant evidence can be stated beyond the standard CDM paradigm, and we conclude that the present observational data favor DM as a pressureless fluid. In all cases, we find the tighter upper bounds on the sum of neutrino masses. The most tight upper bound is $\sum m_\nu < 0.17$ eV at 95\% CL, imposed with the addition of HST prior in the analysis. We also observe the effects of neutrino properties on the extended DM parameters and also on other parameters. We find significantly lower mean values of $\sigma_8$ in all cases which are in agreement with the LSS measurements. Thus, the well-known $\sigma_8$ tension is reconciled in the considered model., Comment: 10 pages, 4 Figure panels. arXiv admin note: text overlap with arXiv:1901.07549

The Bimetric theory of gravity is an extension of General Relativity that describes a massive spin-$2$ particle in addition to the standard massless graviton. The theory is based on two dynamical metric tensors with their interactions constrained by requiring the absence of the so-called Boulware-Deser ghost. It has been realized that the quantum interactions of matter fields with gravity are bound to generate modifications to the standard Einstein-Hilbert action such as quadratic curvature terms. Such a quadratic Ricci scalar term is present in the so-called Starobinsky model which has been proven to be rather robust in its inflationary predictions. In the present article we study a generalization of the Starobinsky model within the bimetric theory and find that its inflationary behaviour stays intact while keeping all consistency requirements of the bimetric framework. The interpretation of the massive spin-2 particle as dark matter remains a viable scenario, as in standard Bigravity., 7 pages, 3 figures

Filaments are elongated structures that connect groups and clusters of galaxies and are visually the striking feature in cosmological maps. In the literature, typically filaments are defined only using galaxies, assuming that these are good tracers of the dark matter distribution, despite the fact that galaxies are a biased indicator. Here we apply the topological filament extractor DisPerSE to the predictions of the semi-analytic code GAEA to investigate the correspondence between the properties of $z=0$ filaments extracted using the distribution of dark matter and the distribution of model galaxies evolving within the same large-scale structure. We focus on filaments around massive clusters with a mass comparable to Virgo and Coma, with the intent of investigating the influence of massive systems and their feeding filamentary structure on the physical properties of galaxies. We apply different methods to compare the properties of filaments based on the different tracers and study how the sample selection impacts the extraction. Overall, filaments extracted using different tracers agree, although they never coincide totally. We also find that the number of filaments ending up in the massive clusters identified using galaxies distribution is typically underestimated with respect to the corresponding dark matter filament extraction., resubmitted to MNRAS after minor comments

Recently the international pulsar timing array collaboration has announced the first strong evidence for an isotropic gravitational wave background (GWB). We propose that rapid small oscillations (wiggles) in the Hubble parameter would trigger a resonance with the propagating gravitational waves, leaving novel signature in the GWB spectrum in the form of sharp resonance peaks. The proposed signal can appear at all frequency ranges and is common to continuous spectrum GWBs with arbitrary origin. Due to its resonant nature, the signal strength differs by a perturbation order depending on whether the GWB is primordial or not, which makes it a smoking gun for the primordial origin of the observed GWB. We show that a large part of the parameter space of such signal can be constrained by near future PTA observations, while fitting the signal template to the current NanoGrav 15yr data already hints an interesting feature near 15 nHz., 10 pages plut appendix, 4 figures

We explore the complementarity of direct detection (DD) and spallation source (SS) experiments for the study of sterile neutrino physics. We focus on the sterile baryonic neutrino model: an extension of the Standard Model that introduces a massive sterile neutrino with couplings to the quark sector via a new gauge boson. In this scenario, the inelastic scattering of an active neutrino with the target material in both DD and SS experiments gives rise to a characteristic nuclear recoil energy spectrum that can allow for the reconstruction of the neutrino mass in the event of a positive detection. We first derive new bounds on this model based on the data from the COHERENT collaboration on CsI and LAr targets, which we find do not yet probe new areas of the parameter space. We then assess how well future SS experiments will be able to measure the sterile neutrino mass and mixings, showing that masses in the range 15-50 MeV can be reconstructed. We show that there is a degeneracy in the measurement of the sterile neutrino mixing that substantially affects the reconstruction of parameters for masses of the order of 40 MeV. Thanks to their lower energy threshold and sensitivity to the solar tau neutrino flux, DD experiments allow us to partially lift the degeneracy in the sterile neutrino mixings and considerably improve its mass reconstruction down to 9 MeV. Our results demonstrate the excellent complementarity between DD and SS experiments in measuring the sterile neutrino mass and highlight the power of DD experiments in searching for new physics in the neutrino sector., 22 pages, 10 figures

The combined study of the flaring of Active Galactic Nuclei (AGN) at radio wavelengths and pc-scale jet kinematics with Very Long Baseline Interferometry (VLBI) has led to the view that i) the observed flares are associated with ejections of synchrotron blobs from the core, and ii) most of the flaring would follow a one-to-one correlation with the component ejection. Recent results have provided mounting evidence that the quasi-regular component injections into the relativistic jet may not be the only cause of the flux variability. We propose that AGN flux variability and jet morphology changes can both be of deterministic nature, i.e. having a geometric/kinetic origin linked to the time-variable Doppler beaming of the jet emission as its direction changes due to precession (and nutation). The physics of the underlying jet leads to shocks, instabilities, or to ejections of plasmoids. The appearance (morphology, flux, etc.) of the jet can, however, be strongly affected and modulated by precession. We demonstrate this modulating power of precession for OJ 287. For the first time, we show that the spectral state of the Spectral Energy Distribution (SED) can be directly related to the jet's precession phase. We model the SED evolution and reproduce the precession parameters. Further, we apply our precession model to eleven prominent AGN. We show that for OJ 287 precession seems to dominate the long-term variability ($\gtrsim 1\,{\rm yr}$) of the AGN flux, SED spectral state, and jet morphology, while stochastic processes affect the variability on short timescales ($\lesssim 0.2\,{\rm yr}$)., 48 pages, 26 figures, 14 tables; published in the Astrophysical Journal

We present a numerical analysis of the cosmological evolution of scalar field dark matter (SFDM) in the Boltzmann code $\texttt{CLASS}$, based on a dynamical system analysis of previous works. We show a detailed study of the evolution of the different dynamical variables, and in particular of the energy density and its corresponding linear perturbations. The numerical results are in good agreement with those of the original SFDM equations of motion, and have better accuracy than other approaches. In addition, we calculate the temperature and matter power spectra and discuss the reliability of their numerical results. We also give simple examples in which we can put constraints on the field mass using recent likelihoods incorporated in the Monte Carlo Markov Chain sampler $\texttt{MontePython}$., 20 pages, 12 figures

First-order phase transitions, which take place when the symmetries are predominantly broken (and masses are then generated) through radiative corrections, produce observable gravitational waves and primordial black holes. We provide a model-independent approach that is valid for large-enough supercooling to quantitatively describe these phenomena in terms of few parameters, which are computable once the model is specified. The validity of a previously-proposed approach of this sort is extended here to a larger class of theories. Among other things, we identify regions of the parameter space that correspond to the background of gravitational waves recently detected by pulsar timing arrays (NANOGrav, CPTA, EPTA, PPTA) and others that are either excluded by the observing runs of LIGO and Virgo or within the reach of future gravitational wave detectors. Furthermore, we find regions of the parameter space where primordial black holes produced by large over-densities due to such phase transitions can account for dark matter. Finally, it is shown how this model-independent approach can be applied to specific cases, including a phenomenological completion of the Standard Model with right-handed neutrinos and gauged $B-L$ undergoing radiative symmetry breaking., 32 pages, 9 figures, for the dataset see https://zenodo.org/record/8128176

The proliferation of topological defects is a common out-of-equilibrium phenomenon when a system is driven into a phase of broken symmetry. The Kibble-Zurek mechanism (KZM) provides a theoretical framework for the critical dynamics and generation of topological defects in such scenarios. One of the early applications of KZM is the estimation of heavy magnetic monopoles left behind by the cosmological phase transitions in the early universe. The scarcity of such relic monopoles, which contradicts the prediction of KZM, is one of the main motivations for cosmological inflationary theories. On the other hand, magnetic monopoles as emergent quasi-particles have been observed in spin ices, a peculiar class of frustrated magnets that remain disordered at temperatures well below the energy scale of exchange interaction. Here we study the annihilation dynamics of magnetic monopoles when spin ice is cooled to zero temperature in a finite time. Through extensive Glauber dynamics simulations, we find that the density of residual monopole follows a power law dependence on the annealing rate. A kinetic reaction theory that precisely captures the annihilation process from Monte Carlo simulations is developed. We further show that the KZM can be generalized to describe the critical dynamics of spin ice, where the exponent of the power-law behavior is determined by the dynamic critical exponent $z$ and the cooling protocol., 13 pages, 7 figures

Spectroscopic analysis of Ly$α$ damping wings of bright sources at $z>6$ is a promising way to measure the reionization history of the universe. However, the theoretical interpretation of the damping wings is challenging due to the inhomogeneous nature of the reionization process and the proximity effect of bright sources. In this Letter, we analyze the damping wings arising from the neutral patches in the radiative transfer cosmological simulation suite Cosmic Reionization on Computers (CROC). We find that the damping wing profile remains a tight function of volume-weighted neutral fraction $\left< x_{\rm HI} \right>_{\rm v}$, especially when $\left< x_{\rm HI} \right>_{\rm v}>0.5$, despite the patchy nature of reionization and the proximity effect. This small scatter indicates that with a well-measured damping wing profile, we could constrain the volume-weighted neutral fraction as precise as $Δ\left< x_{\rm HI} \right>_{\rm v} \lesssim 0.1$ in the first half of reionization., 4 pages, 4 figures, submitted to MNRAS Letter

Many symmetry breaking patterns in grand unified theories (GUTs) give rise to cosmic strings that eventually decay when pairs of GUT monopoles spontaneously nucleate along the string cores. These strings are known as metastable cosmic strings and have intriguing implications for particle physics and cosmology. In this article, we discuss the current status of metastable cosmic strings, with a focus on possible GUT embeddings and connections to inflation, neutrinos, and gravitational waves (GWs). The GW signal emitted by a network of metastable cosmic strings in the early universe differs, in particular, from the signal emitted by topologically stable strings by a suppression at low frequencies. Therefore, if the underlying symmetry breaking scale is close to the GUT scale, the resulting GW spectrum can be accessible at current ground-based interferometers as well as at future space-based interferometers, such as LISA, and at the same time account for the signal in the most recent pulsar timing data sets. Metastable cosmic strings thus nourish the hope that future GW observations might shed light on fundamental physics close to the GUT scale., 28 pages, 3 figures

The ability to compress observational data and accurately estimate physical parameters relies heavily on informative summary statistics. In this paper, we introduce the use of mutual information (MI) as a means of evaluating the quality of summary statistics in inference tasks. MI can assess the sufficiency of summaries, and provide a quantitative basis for comparison. We propose to estimate MI using the Barber-Agakov lower bound and normalizing flow based variational distributions. To demonstrate the effectiveness of our method, we compare three different summary statistics (namely the power spectrum, bispectrum, and scattering transform) in the context of inferring reionization parameters from mock images of 21~cm observations with Square Kilometre Array. We find that this approach is able to correctly assess the informativeness of different summary statistics and allows us to select the optimal set of statistics for inference tasks., Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics, comments welcome

We study how well perturbative forward modeling can constrain cosmological parameters compared to conventional analyses. We exploit the fact that in perturbation theory the field-level posterior can be computed analytically in the limit of small noise. In the idealized case where the only relevant parameter for the nonlinear evolution is the nonlinear scale, we argue that information content in this posterior is the same as in the $n$-point correlation functions computed at the same perturbative order. In the real universe other parameters can be important, and there are possibly enhanced effects due to nonlinear interactions of long and short wavelength fluctuations that can either degrade the signal or increase covariance matrices. We identify several different parameters that control these enhancements and show that for some shapes of the linear power spectrum they can be large. This leads to degradation of constraints in the standard analyses, even though the effects are not dramatic for a $\Lambda$CDM-like cosmology. The aforementioned long-short couplings do not affect the field-level inference which remains optimal. Finally, we show how in these examples calculation of the perturbative posterior motivates new estimators that are easier to implement in practice than the full forward modelling but lead to nearly optimal constraints on cosmological parameters., Comment: 33 pages, 1 figure

After the successful full-shape analyses of BOSS data using the Effective Field Theory of Large-Scale Structure, we investigate what upcoming galaxy surveys might achieve. We introduce a ``perturbativity prior" that ensures that loop terms are as large as theoretically expected, which is effective in the case of a large number of EFT parameters. After validating our technique by comparison with already-performed analyses of BOSS data, we provide Fisher forecasts using the one-loop prediction for power spectrum and bispectrum for two benchmark surveys: DESI and MegaMapper. We find overall great improvements on the cosmological parameters. In particular, we find that MegaMapper (DESI) should obtain at least a 12$\sigma$ ($2\sigma$) evidence for non-vanishing neutrino masses, bound the curvature $\Omega_k$ to 0.0012 (0.012), and primordial inflationary non-Gaussianities as follows: $f_{\text{NL}}^{\text{loc.}}$ to $\pm 0.26$ (3.3), $f_{\text{NL}}^{\text{eq.}}$ to $\pm16$ (92), $f_{\text{NL}}^{\text{orth.}}$ to $\pm 4.2$ (27). Such measurements would provide much insight on the theory of Inflation. We investigate the limiting factor of shot noise and ignorance of the EFT parameters., Comment: 38+14 pages, 12 figures, 8 tables

Recent observations of weak gravitational lensing surveys indicate a smoother Universe compared to the predictions of the Cosmic Microwave Background (CMB). This is known as $\sigma_8$ tension or $S_8$ tension, where $\sigma_8$ represents the present root-mean-square matter fluctuation averaged over a sphere of radius $8 h^{-1} \mathrm{Mpc}$ and $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}$. In this letter, we investigate a kind of general Dirac-Born-Infeld (DBI) Lagrangian referred as surface-type DBI (s-DBI) model. We have found that, up to the linear order, the constraints on the s-DBI model with CMB from Planck2018 and low-redshift probes (WL and GC) yield $S_8= 0.7685_{-0.0066}^{+0.0077}$ and $S_8=0.766_{-0.0376}^{+0.0471}$, respectively, which are not only self-consistent but also consistent with the values derived from most low-redshift probes. Furthermore, we provide an outlook for searching the non-linear effects of this model, which could be helpful to resolve other issues by Cold Dark Matter on small scales., Comment: Submitted, comments welcome

With focus on the cosmological evolution of linear perturbations of matter and geometry, we calculate the equivalent expressions to that of the Newtonian and Synchronous gauges within the framework of Unimodular Gravity, being these two gauges commonly used and implemented in Boltzmann codes. An important aspect of our analysis is the inclusion of the energy-momentum current violation, as well as its perturbations. Moreover, for the first time we demonstrate that it is possible to fix both gauges consistently, although as it has been already noticed in previous literature, neither of them is recovered in the sense of the dynamics given in General Relativity for matter and metric fluctuations. Specifically, we show that since the unimodular constraint at the level of linear perturbations lead to only one degree of freedom of scalar modes of metric fluctuations, the dynamics in Unimodular Gravity forces to keep the anisotropic stress in the Newtonian gauge, whereas the cold dark matter comoving frame can not be set in the Synchronous gauge. The physical implications on the density contrast of cold dark matter is reviewed, and the Sachs-Wolfe effect is obtained and compared with previous results in the literature of cosmological perturbations in Unimodular Gravity., 31 pages, new references added

We perform a Bayesian analysis of NANOGrav 15yr and IPTA DR2 pulsar timing residuals and show that the recently detected stochastic gravitational-wave background (SGWB) is compatible with a SGWB produced by bubble dynamics during a cosmological first-order phase transition. The timing data suggests that the phase transition would occur around QCD confinement temperature and would have a slow rate of completion. This scenario can naturally lead to the abundant production of primordial black holes (PBHs) with solar masses. These PBHs can potentially be detected by current and advanced gravitational wave detectors LIGO-Virgo-Kagra, Einstein Telescope, Cosmic Explorer, by astrometry with GAIA and by 21-cm survey., 5 pages, 4 figures + appendices

We infer the connection between the stellar mass of galaxies from the Subaru Hyper Suprime-Cam (HSC) survey, and their dark matter halo masses and its evolution in two bins of redshifts between $[0.3, 0.8]$. We use the measurements of the weak lensing signal of galaxies using background sources from the Year 1 shape catalog from the HSC survey. We bin galaxies in stellar mass with varying thresholds ranging from $8.6 \leq \log [ M_*/(h^{-2} {M_\odot})] \leq 11.2$ and use stringent cuts in the selection of source galaxies to measure the weak lensing signal. We model these measurements of the weak lensing signal together with the abundance of galaxies in the halo occupation distribution framework. We obtain constraints on the halo occupation parameters of central galaxies $M_{\rm min}$ and $\sigma_{\log M}$, which correspond to the halo mass at which central galaxies for each threshold sample reach half occupancy, and its scatter, respectively, along with parameters that describe the occupation of the satellite galaxies. The measurements of abundance and weak lensing individually constrain different degeneracy directions in the $M_{\rm min}$ and $\sigma_{\log M}$ plane, thus breaking the degeneracy in these parameters. We demonstrate that the weak lensing measurements are best able to constrain the average central halo masses, $\langle M_{\rm cen} \rangle$. We compare our measurements to those obtained using the abundance and clustering of these galaxies as well as the subhalo abundance matching measurements and demonstrate qualitative agreement. We find that the galaxy-dark matter connection does not vary significantly between redshift bins we explore in this study. Uncertainties in the photometric redshift of the lens galaxies imply that more efforts are required to understand the true underlying stellar mass-halo mass relation of galaxies and its evolution over cosmic epoch.