Search

We report on the results of performing computational gravitational backreaction on cosmic string loops taken from a network simulation. The principal effect of backreaction is to smooth out small-scale structure on loops, which we demonstrate by various measures including the average loop power spectrum and the distribution of kink angles on the loops. Backreaction does lead to self-intersections in most cases, but these are typically small. An important effect discussed in prior work is the rounding off of kinks to form cusps, but we find that the cusps produced by that process are very weak and do not significantly contribute to the total gravitational-wave radiation of the loop. We comment briefly on extrapolating our results to loops as they would be found in nature., Comment: 26 pages, 11 figures. Power spectrum data available via Zenodo: https://doi.org/10.5281/zenodo.14037540

We calculate the two-loop corrections in curvature perturbation power spectrum in models of single field inflation incorporating an intermediate USR phase employed for PBHs formation. Among the total eleven one-particle irreducible Feynman diagrams, we calculate the corrections from the ``double scoop" two-loop diagram involving two vertices of quartic Hamiltonians. We demonstrate that the fractional two-loop correction in power spectrum scales like the square of the fractional one-loop correction. We confirm our previous findings that the loop corrections become arbitrarily large when the transition from the USR phase to the final attractor phase is very sharp. This suggests that in order for the analysis to be under perturbative control against the loop corrections, one requires a mild transition with a long enough relaxation period towards the final attractor phase., Comment: 21 pages, 2 figures. arXiv admin note: text overlap with arXiv:2403.03841

Fixed target experiments where beam electrons are focused upon a thin target have shown great potential for probing new physics, including the sub-GeV dark matter (DM) paradigm. However, a signal in future experiments such as the light dark matter experiment (LDMX) would require an independent validation to assert its DM origin. To this end, we propose to combine LDMX and next generation DM direct detection (DD) data in a four-step analysis strategy, which we here illustrate with Monte Carlo simulations. In the first step, the hypothetical LDMX signal (i.e. an excess in the final state electron energy and transverse momentum distributions) is $\textit{recorded}$. In the second step, a DM DD experiment operates with increasing exposure to test the DM origin of the LDMX signal. Here, LDMX and DD data are simulated. In the third step, a posterior probability density function (pdf) for the DM model parameters is extracted from the DD data, and used to $\textit{predict}$ the electron recoil energy and transverse momentum distributions at LDMX. In the last step, $\textit{predicted}$ and $\textit{recorded}$ electron recoil energy and transverse momentum distributions are compared in a chi-square test. We present the results of this comparison in terms of a threshold exposure that a DD experiment has to operate with to assert whether $\textit{predicted}$ and $\textit{recorded}$ distributions $\textit{can}$ be statistically dependent. We find that this threshold exposure grows with the DM particle mass, $m_\chi$. It varies from 0.012 kg-year for a DM mass of $m_\chi=4$ MeV to 1 kg-year for $m_\chi=25$ MeV, which is or will soon be within reach.

This study explores a two-component dark matter model in which one component, heavier dark matter, annihilates into a lighter dark matter. The lighter dark matter is expected to generate detectable signals in detectors due to its enhanced momentum, enabling direct detection even for MeV-scale dark matter. We investigate the effectiveness of directional direct detections, especially the nuclear emulsion detector NEWSdm, in verifying these boosted dark matter particles through nuclear recoil. In particular, we focus on light nuclei, such as protons and carbon, as suitable targets for this detection method due to their high sensitivity to MeV-scale dark matter. By modeling the interactions mediated by a dark photon in a hidden U(1)$_D$ gauge symmetry framework, we calculate the expected dark matter flux and scattering rates for various detector configurations. Our results show that nuclear emulsions have the potential to yield distinct, direction-sensitive dark matter signals from the Galactic center, providing a new way to probe low-mass dark matter parameter spaces that evade conventional detection methods., Comment: 17 pages, 10 figures, 3 tables

We investigate the impact of beam chromaticity, i.e., the frequency dependence of the beam window function, on cosmological and astrophysical parameter constraints from CMB power spectrum observations. We show that for future high-resolution CMB measurements it is necessary to include a color-corrected beam for each sky component with a distinct spectral energy distribution. We introduce a formalism able to easily implement the beam chromaticity in CMB power spectrum likelihood analyses and run a case study using a Simons Observatory (SO) Large Aperture Telescope-like experimental setup and within the public SO software stack. To quantify the impact, we assume that beam chromaticity is present in simulated spectra but omitted in the likelihood analysis. We find that, for passbands of fractional width $\Delta \nu/\nu \sim 0.2$, neglecting this effect leads to significant biases, with astrophysical foreground parameters shifting by more than $2\sigma$ and cosmological parameters by significant fractions of the error., Comment: 10 pages, 6 figures

A sharp step on a chaotic potential can enhance primordial curvature fluctuations on smaller scales to the $\mathcal{O}(10^{-2})$ to form primordial black holes (PBHs). The present study discusses an inflationary potential with a sharp step that results in the formation of PBHs in four distinct mass ranges. Also this inflationary model allows the separate consideration of observable parameters $n_s$ and $r$ on the CMB scale from the physics at small scales, where PBHs formation occur. In this work we computed the fractional abundance of PBHs ($f_{PBH}$) using the GLMS approximation of peak theory and also the Press-Schechter (PS) formalism. In the two typical mass windows, $10^{-13}M_\odot$ and $10^{-11}M_\odot$, $f_{PBH}$ calculated using the GLMS approximation is nearly equal to 1 and that calculated via PS is of $10^{-3}$. In the other two mass windows $1M_\odot$ and $6M_\odot$, $f_{PBH}$ obtained using GLMS approximation is 0.01 and 0.001 respectively, while $f_{PBH}$ calculated via PS formalism yields $10^{-5}$ and $10^{-6}$. The results obtained via GLMS approximation are found to be consistent with observational constraints. A comparative analysis of $f_{PBH}$ obtained using the GLMS perspective and the PS formalism is also included.

We revisit unflavoured leptogenesis in the seesaw model applying recent improvements in the computation of $CP$-conserving and $CP$-violating equilibration rates. These are relevant for the relativistic regime of the sterile Majorana fermions and the dynamics of the Standard Model particles acting as spectator processes. In order to probe the regime of large (${\cal O}(10^2)$) washout parameters, we add $\Delta L = 2$ washout processes, which we derive in the CTP-formalism. We then perform a parameter scan of the final baryon asymmetry and find a constraint $m_{\text{lightest}} \lesssim 0.15 \, \rm eV$ on the absolute neutrino mass scale, which is slightly less stringent than previously reported bounds obtained without the aforementioned improvements. The relaxation of the bounds is mainly due to partially equilibrated spectator fields, which protect part of the asymmetry from washout and lead to larger final asymmetries., Comment: 14 pages, 4 figures

We explore the potential of gravitational waves (GWs) to probe the pre-BBN era of the early universe, focusing on the effects of energy injection. Specifically, we examine a hidden sector alongside the Standard Model that undergoes a strong first-order phase transition (FOPT), producing a GW signal. Once the phase transition has completed, energy injection initiates reheating in the hidden sector, which positions the hidden sector field so that additional phase transitions can occur. This can result in a total of three distinct phase transitions with a unique three-peak GW spectrum. Among these transitions, the first and third are of the standard type, while the intermediate second transition is inverted, moving from a broken to an unbroken phase. Using polynomial potentials as a framework, we derive analytical relations among the phase transition parameters and the resulting GW spectrum. Our results indicate that the second and third transitions generate GWs with higher amplitudes than the first, with a peak frequency ratio differing by up to an order of magnitude. This three-peak GW spectrum is detectable by upcoming facilities such as LISA, BBO, and UDECIGO. Notably, the phenomenon is robust across various potentials and model parameters, suggesting that hidden sector GWs provide a powerful tool for exploring new physics scenarios in the pre-BBN era.

Early dark energy solutions to the Hubble tension introduce an additional scalar field which is frozen at early times but becomes dynamical around matter-radiation equality. In order to alleviate the tension, the scalar's share of the total energy density must rapidly shrink from $\sim 10\%$ at the onset of matter domination to $\ll 1\%$ by recombination. This typically requires a steep potential that is imposed $\textit{ad hoc}$ rather than emerging from a concrete particle physics model. Here, we point out an alternative possibility: a homogeneous scalar field coupled quadratically to a cosmological background of light thermal relics (such as the Standard Model neutrino) will acquire an effective potential which can reproduce the dynamics necessary to alleviate the tension. We identify the relevant parameter space for this "thermo-coupled" scenario and study its unique phenomenology at the background level, including the back-reaction on the neutrino mass. Follow-up numerical work is necessary to determine the constraints placed on the model by early-time measurements., Comment: 6 pages, 3 figures

The mass of galaxy clusters is a critical quantity for probing cluster cosmology and testing theories of gravity, but its measurement could be biased given assumptions are inevitable. In this paper, we employ and compare two mass proxies for galaxy clusters: thermodynamics of the intracluster medium and kinematics of member galaxies. We select 22 galaxy clusters from the cluster catalog in the first SRG/eROSITA All-Sky Survey (eRASS1) that have sufficient optical and near-infrared observations. We generate multi-band images in the energy range of (0.3, 7) keV for each cluster, and derive their temperature profiles, gas mass profiles and hydrostatic mass profiles using a parametric approach that does not assume dark matter halo models. With spectroscopically confirmed member galaxies collected from multiple surveys, we numerically solve the spherical Jeans equation for their dynamical mass profiles. Our results quantify the correlation between dynamical mass and line-of-sight velocity dispersion with an rms scatter of 0.14 dex. We find the two mass proxies lead to roughly the same total mass, with no observed systematic bias. As such, the $\sigma_8$ tension is not specific to hydrostatic mass or weak lensing shears, but also appears with galaxy kinematics. We also compare our hydrostatic masses with the latest weak lensing masses inferred with scaling relations. The comparison shows the weak lensing mass is significantly higher than our hydrostatic mass by $\sim$110%. This might explain the significantly larger value of $\sigma_8$ from the latest measurement using eRASS1 clusters than almost all previous estimates in the literature. Finally, we test the radial acceleration relation (RAR) established in disk galaxies. We confirm the missing baryon problem in the inner region of galaxy clusters using three independent mass proxies for the first time., Comment: 14 pages, 9 figures, 2 tables, with 12-page appendix for additional figures and tables. Accepted for publication in A&A

We present the first set of high-resolution, hydrodynamical cosmological simulations of galaxy formation in a Fuzzy Dark Matter (FDM) framework. These simulations were performed with a new version of the GASOLINE2 code, known as FUZZY-GASOLINE, which can simulate quantum FDM effects alongside a comprehensive baryonic model that includes metal cooling, star formation, supernova feedback, and black hole physics, previously used in the NIHAO simulation suite. Using thirty zoom-in simulations of galaxies with halo masses in the range $10^9 \lesssim M_{\text{halo}}/M_{\odot} \lesssim 10^{11}$, we explore how the interplay between FDM quantum potential and baryonic processes influences dark matter distributions and observable galaxy properties. Our findings indicate that both baryons and low-mass FDM contribute to core formation within dark matter profiles, though through distinct mechanisms: FDM-induced cores emerge in all haloes, particularly within low-mass systems at high redshift, while baryon-driven cores form within a specific mass range and at low redshift. Despite these significant differences in dark matter structure, key stellar observables such as star formation histories and velocity dispersion profiles remain remarkably similar to predictions from the Cold Dark Matter (CDM) model, making it challenging to distinguish between CDM and FDM solely through stellar observations., Comment: 13 pages, 8 figures, 6 tables. Paper submitted to MNRAS

In the year 1924, a paper by Carl Wirtz appeared in Astronomische Nachrichten, entitled 'De Sitter's cosmology and the radial motion of spiral galaxies'. This paper and its author remained largely unnoticed by the community, but it seems to be the first cosmological interpretation of the redshift of galaxies as a time dilation effect and the expansion of the Universe. Edwin Hubble knew Wirtz' publications quite well. The modern reader would find Wirtz' own understanding diffuse and contradictory in some aspects, but that reflected the early literature on nebulae, to which he himself made important contributions. The 100th anniversary provides a good opportunity to present an English transcription to the community, which can be found in the appendix. This anniversary also provokes to ask for the present status of cosmology which many authors see in a crisis. From an observational viewpoint it shall be illustrated that until today there is no consistent/convincing understanding of how the Universe evolved., Comment: 9 pages, no figures accepted by Astronomische Nachrichten

The large-scale structure of the Universe and its evolution over time contains an abundance of cosmological information. One way to unlock this is by measuring the density and momentum power spectrum from the positions and peculiar velocities of galaxies, and fitting the cosmological parameters from these power spectrum. In this paper, we will explore the cross power spectrum between the density and momentum fields of galaxies. We derive the estimator of the density-momentum cross power spectrum multipoles. The growth rate of the large-scale-structure, $f\sigma_8$ is measured from fitting the combined density monopole, momentum monopole and cross dipole power spectrum. The estimators and models of power spectrum as well as our fitting method have been tested using mock catalogues, and we find that they perform well in recovering the fiducial values of the cosmological parameters of the simulations, and we also find that the errors of the parameters can be largely reduced by including the cross-power spectrum in the fit. We measure the auto-density, auto-momentum and cross power spectrum using the Sloan Digital Sky Survey Data Release 14 peculiar velocity catalogue. The fit result of the growth rate $f\sigma_8$ is $f\sigma_8=0.413^{+0.050}_{-0.058}$ at effective redshift $z_{\mathrm{eff}}=0.073$, and our measurement is consistent with the prediction of the $\Lambda$ Cold Dark Matter cosmological model assuming General Relativity., Comment: 15 pages, 8 figures, 2 tables. Publish in ApJ

The 21 cm emission from neutral hydrogen surveys holds great potential as a valuable method for exploring the large-scale structure of the Universe. In this paper, we forecast for the cross-correlation between the Thermal Sunyaev-Zel'dovich (SZ) fluctuations as probed by the Planck satellite, and fluctuations in the HI brightness temperature as probed by the ground-based Five-hundred-meter Aperture Spherical Telescope (FAST), to trace the connection between galaxy clusters and the HI large-scale structure. Assuming that the measurement is limited by instrumental noise rather than by foreground, we estimate the potential detectability of the cross-correlation signal and their improvement in the measurement of the HI cosmic density, the hydrostatic mass bias parameter, and the universal pressure profile (UPP) parameters. We obtain a constraint on the cosmic neutral hydrogen density parameter significantly to $\sigma(\Omega_{\rm HI}) = 1.0 \times 10^{-6}$. We also find that the average halo masses contributing to the ${{\rm HI}-y}$ cross-power spectrum in the one-halo regime is $\sim 1.5\times 10^{14} M_{\odot}$. Our results also show that the HI-SZ cross-correlation has great potential to probe the distribution of neutral hydrogen (HI) within halos at low redshift., Comment: 18 pages, 4 figures, 3 tables, accepted in the Astrophysical Journal

Strongly lensed type Ia supernovae (SNe Ia) provide a unique cosmological probe to address the Hubble tension problem in cosmology. In addition to the sensitivity of the time delays to the value of the Hubble constant, the transient and standard candle nature of SNe Ia also enable valuable joint constraints on the model of the lens and the cosmological parameters. The upcoming Legacy Survey of Space and Time (LSST) with the Vera C. Rubin Observatory is expected to increase the number of observed SNe Ia by an order of magnitude in ten years of its lifetime. However, finding such systems in the LSST data is a challenge. In this work, we revisit the color-magnitude (CM) diagram used previously as a means to identify lensed SNe Ia and extend the work further as follows. We simulate LSST-like photometric data ($rizy$~bands) of lensed SNe Ia and analyze it in the CM parameter space. We find that a subset of lensed SNe Ia are redder compared to unlensed SNe Ia at a given magnitude, both in the rising and falling phases of their light curves and for SNe up to $z=3$. We propose a modified selection criterion based on these new results. We show that the contamination coming from the unlensed core-collapse (CC) SNe is negligible, whereas a small fraction of lensed CC SNe types Ib and Ic may get selected by this criterion as potential lensed SNe. Finally, we demonstrate that our criterion works well on a wide sample of observed unlensed SNe Ia, a handful of known multiply-imaged lensed SNe systems, and a representative sample of observed super-luminous supernovae.

It was hypothesized in the literature that some physical parameters may be time-evolving and the astrophysical data can serve as a probe. Recently, James Webb Space Telescope (JWST) have released its early observations. In this work, we select the JWST spectroscopic observations of the high redshift ($z>7.1$) galaxies with strong [OIII] ($\lambda=4959$ \AA \,and $5007$ \AA \,in the rest frame) emission lines to constraint the evolution of the fine structure constant ($\alpha$). With the spectra from two galaxies at redshifts of $7.19$ and $8.47$, the deviation of $\alpha$ to its fiducial value is found to be as small as $0.44^{+8.4+1.7}_{-8.3-1.7} \times 10^{-4}$ and $-10.0^{+18+1.5}_{-18-1.5} \times 10^{-4}$, respectively (the first error is statistical and the latter is systematic). The combination of our results with the previous data reveals that $\frac{1}{\alpha} \frac{d \alpha}{dt} = 0.30^{+4.5}_{-4.5} \times 10^{-17}~{\rm yr^{-1}}$. Clearly, there is no evidence for a cosmic evolution of $\alpha$. The prospect of further constraining the time evolution of $\alpha$ is also discussed. The scalar field of dark energy is hypothesized to drive the acceleration of the universe's expansion through an interaction with the electromagnetic field. By integrating the observational data of the fine-structure constant variation, $\frac{\Delta\alpha}{\alpha}(z)$, we have established a stringent upper limit on the coupling strength between dark energy and electromagnetism. Our analysis yields $\zeta \leq 3.92 \times 10^{-7}$ at the 95\% confidence level, representing the most stringent bound to date.

Pulsar timing array experiments have recently found evidence for a stochastic gravitational wave (GW) background, which induces correlations among pulsar timing residuals described by the Hellings and Downs (HD) curve. Standard calculations of the HD correlation and its variance assume an isotropic background. However, for a background of astrophysical origin, we expect a higher GW spectral density in directions with higher galaxy number densities. In a companion paper, we have developed a theoretical formalism to account for the anisotropies arising from large-scale galaxy clustering, leading to a new contribution to the variance of the HD correlation. In this subsequent work, we provide numerical results for this novel effect. We consider a GW background resulting from mergers of supermassive black hole binaries, and relate the merger number density to the overdensity of galaxies. We find that anisotropies due to large-scale galaxy clustering lead to a standard deviation of the HD correlation at most at percent level, remaining well below the standard contributions to the HD variance. Hence, this kind of anisotropies in the GW source distribution does not represent a substantial contamination to the correlations of timing residuals in present and future PTA surveys. Suitable statistical methods to extract the galaxy clustering signal from PTA data will be investigated in the future., Comment: 20 pages, 6 figures

We investigate the formation of bound states of non-relativistic dark matter particles subject to long-range interactions through radiative capture. The initial scattering and final bound states are described by Coulomb potentials with different strengths, as relevant for non-abelian gauge interactions or theories featuring charged scalars. For bound states with generic quantum numbers $n$ and $\ell$, we provide closed-form expressions for the bound-state formation (BSF) cross sections of monopole, dipole and quadrupole transitions, and of arbitrary multipole order when $\ell=n-1$. This allows us to investigate in detail a strong enhancement of BSF that occurs for initial states in a repulsive potential. For $\ell=n-1\gg 1$, we show that the BSF cross section for each single bound state violates the perturbative unitarity bound in the vicinity of a certain critical initial velocity, and provide an interpretation in terms of a smooth matching of classical trajectories. When summing the BSF cross section over all possible bound states in the final state, this leads to a unitarity violation below a certain velocity, but within the validity range of the weakly coupled non-relativistic description. We identify an effectively strong interaction as the origin of this unitarity violation, which is caused by an "anomalously" large overlap of scattering and bound-state wave functions in Coulomb potentials of different strength., Comment: 38p

Axions and axion-like particles can be probed through gravitational waves indirectly, often referred to as "audible axions". The usual concept of audible axion relies on the coupling between the axions and the gauge fields. Here we consider an axion-like mechanism with coupling to the Nieh-Yan term. This interaction leads to the direct and efficient production of gravitational waves during the radiation-dominated era, originating from the tachyonic instability of the gravitational perturbations with the Nieh-Yan term. We calculate the energy spectral density of the chiral gravitational wave background and the comoving energy density of axion-like fields. Based on the numerical results, we explore the parameter space of axion masses and decay constants for detectable gravitational wave signals, either in pulsar timing arrays or space-based gravitational wave detections., Comment: 11 pages, 9 figures, 1 table

New insights from the Dark Energy Spectroscopic Instrument (DESI) 2024 baryon acoustic oscillations (BAO) data, in conjunction with cosmic microwave background (CMB) and Type Ia supernova (SN) data, suggest that dark energy may not be a cosmological constant. In this work, we investigate the cosmological implications of holographic dark energy (HDE) and interacting holographic dark energy (IHDE) models, utilizing CMB, DESI BAO, and SN data. By considering the combined DESI BAO and SN data, we determine that in the IHDE model, the parameter $c > 1$ and the dark-energy equation of state $w$ does not cross $-1$ at the $1\sigma$ confidence level, whereas in the HDE model, it marginally falls below this threshold. Upon incorporating CMB data, we observe that in the HDE model, the parameter $c < 1$ and $w$ crosses $-1$ at a level beyond $10\sigma$. Conversely, for the IHDE model, the likelihood of $w$ crossing $-1$ is considerably diminished, implying that the introduction of interaction within the HDE model could potentially resolve or mitigate the cosmic big rip conundrum. Furthermore, our analysis reveals that the HDE and IHDE models are statistically as viable as the $\Lambda$CDM model when assessing Bayesian evidence with DESI BAO data combined with SN data. However, when CMB data are added, the HDE and IHDE models are significantly less favored compared to the $\Lambda$CDM model. Our findings advocate for further exploration of the HDE and IHDE models using forthcoming, more precise late-universe observations., Comment: 13 pages, 5 figures

We analyse the evolution of the reduced density matrix of inflationary perturbations, coupled to a heavy entropic field via the leading-order term within the Effective Field Theory of Inflation, for two nearly de Sitter backgrounds. We perform a full quantum treatment of the open system and derive a Fokker-Planck equation to describe decoherence and the entanglement structure of the adiabatic perturbations. We find that exotic phenomena, such as recoherence and transient negative growth of entanglement entropy, appearing for the attractor solution, are absent for the non-attractor background. We comment on the relationship of these to the non-Markovian nature of the system. Finally, we generalise to the case where a few e-folds of ultra-slow roll evolution are sandwiched between phases of slow-roll inflation to find its (memory) effects on the curvature perturbation., Comment: 42 pages, 14 figures

Jets launched by supermassive black holes transport relativistic leptons, magnetic fields, and atomic nuclei from the centres of galaxies to their outskirts and beyond. These outflows embody the most energetic pathway by which galaxies respond to their Cosmic Web environment. Studying black hole feedback is an astrophysical frontier, providing insights on star formation, galaxy cluster stability, and the origin of cosmic rays, magnetism, and heavy elements throughout the Universe. This feedback's cosmological importance is ultimately bounded by the reach of black hole jets, and could be sweeping if jets travel far at early epochs. Here we present the joint LOFAR-uGMRT-Keck discovery of a black hole jet pair extending over $7$ megaparsecs -- the largest galaxy-made structure ever found. The outflow, seen $7.5$ gigayears into the past, spans two-thirds of a typical cosmic void radius, thus penetrating voids at ${\sim}95\%$ probability. This system demonstrates that jets can avoid destruction by magnetohydrodynamical instabilities over cosmological distances, even at epochs when the Universe was 15 to 7 times denser than it is today. Whereas previous record-breaking outflows were powered by radiatively inefficient active galactic nuclei, this outflow is powered by a radiatively efficient active galactic nucleus, a type common at early epochs. If, as implied, a population of early void-penetrating outflows existed, then black hole jets could have overwritten the fields from primordial magnetogenesis. This outflow shows that energy transport from supermassive black holes operates on scales of the Cosmic Web and raises the possibility that cosmic rays and magnetism in the intergalactic medium have a non-local, cross-void origin., Comment: 30 pages, 10 figures, 1 table, published in Nature. This arXiv version pre-dates the peer review process

$F(R)$ models for dark energy generally exhibit a weak curvature singularity, which can be cured by adding an $R^2$ term. This correction allows for a unified description of primordial and late-time accelerated expansions. However, most existing models struggle to achieve this, as they become unstable over certain negative ranges of the Ricci scalar, where either the first or second derivative of $F(R)$ turns negative. These instabilities may disrupt the post-inflationary evolution when the Ricci scalar oscillates about the vacuum state after the $R^2$ inflation. In this work, we introduce a new model-building to guarantee global stability, i.e., the first and second derivatives are positive for all real Ricci scalars. By extending the idea from Appleby and Battye, we demonstrate that viable models can be constructed by imposing a positive, bounded first derivative of $F(R)$ with a sigmoid shape. As examples, we first reformulate and generalize the original Appleby-Battye model. Then, we propose a new dark energy model, which successfully explains the acceleration of cosmic expansion and passes local gravity tests., Comment: 31 pages, 3 figures

Axion-photon oscillation effect provides a possible explanation for the presence of very-high-energy (VHE) $\gamma$-ray signals from distant sources. In this work, we propose a model-dependent method to select possible sources that may give sufficient constraints on the axion parameters. We investigate such effect in the spectra of active galactic nuclei (AGN) B2 2234+28A and 3C 454.3 based on data obtained from Fermi Large Area Telescope (Fermi-LAT) and MAGIC U.L. We utilize the Markov Chain Monte Carlo method to fit the axion parameters, yielding a result of $g_{a\gamma}=3.05^{+0.51}_{-0.31} \times 10^{-11}$ GeV$^{-1}$ for the axion-photon coupling strength and $m_{a}=5.25^{+2.35}_{-2.65} \times 10^{-8} $ eV for the axion mass. We also perform 95\% confidence level (CL) constraints to set an upper limit for $g_{a\gamma}$., Comment: 11 pages, 10 figures, 1 table

The nature of dark matter is poorly constrained on subgalactic scales. Alternative models to cold dark matter, such as warm dark matter or self-interacting dark matter, could produce very different dark haloes on these scales. One of the few known dark haloes smaller than a galaxy was discovered in the triple source plane strong lens system J0946+1006. Previous studies have found that this structure is much more concentrated than expected in $\Lambda$CDM, but have assumed the dark halo is at the same redshift as the main deflector ($z_{\rm main}=0.222$). In this paper, we fit for the redshift of this dark halo. We reconstruct the first two sources in the system using a forward modelling approach, allowing for additional complexity from multipole perturbations. We find that the perturber redshift is $z_{\rm halo} = {0.207}^{+0.019}_{-0.019}$, and lower bounds on the evidence strongly prefer a subhalo over a line-of-sight structure. Whilst modelling both background sources does not improve constraints on the redshift of the subhalo, it breaks important degeneracies affecting the reconstruction of multipole perturbations. We find that the subhalo is a more than $5\sigma$ outlier from the $\Lambda$CDM $v_{\rm max}$-$r_{\rm max}$ relation and has a steep profile with an average slope of $\gamma_{\rm 2D} = {-1.81}^{+0.15}_{-0.11}$ for radii between $0.75-1.25$ kpc. This steep slope might indicate dark matter self-interactions causing the subhalo to undergo gravothermal collapse; such collapsed haloes are expected to have $\gamma_{\rm 2D} \approx -2$., Comment: 16 pages, 10 figures, 2 tables

The Cosmological Principle is a cornerstone of the standard model of cosmology and shapes how we view the Universe and our place within it. It is imperative, then, to devise multiple observational tests which can identify and quantify possible violations of this foundational principle. One possible method of probing large-scale anisotropies involves the use of weak gravitational lensing. We revisit this approach in order to analyse the imprint of late-time anisotropic expansion on cosmic shear. We show that the cross-correlation of shear $E$- and $B$-modes on large scales can be used to constrain the magnitude (and possibly direction) of anisotropic expansion. We estimate the signal to noise for multipoles $10\lesssim \ell\lesssim 100$ that is achievable by a Euclid-like survey. Our findings suggest that such a survey could detect the $E$-$B$ signal for reasonable values of the late-time anisotropy parameter.

Determining the spatial curvature $\Omega_K$ of the Universe has long been crucial in cosmology. In practice, this effort is often entangled with assumptions of dark energy. A combination of distance ($D_{\rm M}$, $D_{\rm L}$) and expansion rate ($H(z)$) measurements can break this degeneracy. However, fitting against discrete data points requires parameterizations of distance and expansion rate as functions of redshifts, which often induces cosmological model dependence. In this work, we propose a new dark energy model-independent parameterization of the cosmological comoving radial distance $\chi$. Fitting data combining distance ($D_{\rm M}$, $D_{\rm L}$) and Hubble parameter (or equivalently $D_H$) measurements, we are then able to obtain $\Omega_K$ in a dark energy model-independent manner. We test this parameterization and the associated fitting scheme with mock data generated with a wide range of fiducial dark energy equations of state ($-1.3

We demonstrate that a broad class of modular inflation models predicts the emergence of new physics within an energy range of approximately $10^{15}\, \mathrm{GeV}$ to $10^{17} \, \mathrm{GeV}$. This prediction arises by comparing the moduli-dependent species scale with observational constraints on inflation. Specifically, we illustrate this within the context of $SL(2, \mathbb{Z})$-modular inflation models by re-expressing inflationary observables in terms of the species scale. We further discuss the implications of this approach for generic Calabi-Yau threefolds, showing that this reformulation allows us to directly constrain the fundamental parameters related to the geometry of extra dimensions, specifically the second Chern numbers., Comment: 6 pages, 1 figure

We critically assess the impact of significant dipole and large-scale anisotropies on galaxy clustering signals, with a focus on radio continuum surveys. Our study reveals that these anisotropies -- resulting from intrinsic cosmological effects and/or observational systematics -- profoundly influence the two-point correlation function (2PCF) and angular power spectrum ($C_\ell$). Notably, large-scale anisotropies can obscure or simulate non-Gaussianity signals, complicating the extraction of precise cosmological information. The results emphasize that it is crucial to address systematics and rigorously mask the dipole and its surrounding multipoles to obtain accurate cosmological constraints. This approach is essential for extracting cosmological results from clustering signals, particularly for future surveys such as SKA, DESI, and LSST, to ensure the precision and reliability of cosmological analyses., Comment: 5 pages, 2 figures

The relationship between galaxies and their supermassive black holes (SMBHs) is an area of active research. One way to investigate this is to compare parsec-scale jets formed by SMBHs with the projected shape of their kiloparsec-scale host galaxies. We analyse Very Long Baseline Interferometry (VLBI) images of Active Galactic Nuclei (AGN) and optical images of their host galaxies. We compare the inner-jet position angle in VLBI-detected radio sources with the optical shapes of galaxies as measured by several large optical surveys. In total 6273 galaxy-AGN pairs were found. We carefully account for the systematics of the cross-matched sources and find that Dark Energy Spectroscopic Instrument Legacy Imaging Surveys data (DESI LS) is significantly less affected by them. Using DESI LS, with which 5853 galaxy-AGN pairs were cross-matched, we find a weak but significant alignment signal (with a p-value $\leqslant$ 0.01) between the parsec-scale AGN jet and the kpc projected minor axis of the optical host galaxy in sources with measured spectroscopic redshifts. Our results show that the observed source properties are connected over 3 orders of magnitude in scale. This points towards an intimate connection between the SMBH, their host galaxies and their subsequent evolution., Comment: 37 pages, 15 figures, published in Nature Astronomy: https://www.nature.com/articles/s41550-024-02407-4

The formation of supermassive black holes (SMBHs) in the early Universe is a subject of significant debate. In this study, we examine whether non-evaporating primordial black holes (PBHs) can offer a solution. We establish initial constraints on the range of PBH masses that correspond to Hawking radiation (HR) effective temperatures in the range needed to avoid the fragmentation of primordial gas into smaller, stellar-mass black holes. We also investigate the specific intensity of the HR from non-evaporating PBHs and compare it with the critical radiation needed for direct collapse black holes (DCBHs). We show that HR from non-evaporating PBHs cannot serve as the heating mechanism to facilitate the formation of the seeds for the SMBHs we observe in the high-redshift Universe unless, perhaps, the PBHs within the relevant mass range comprise a significant fraction of dark matter and are significantly clustered towards the center of the primordial halo., Comment: 16 pages, 6 figures. Project based on an undergradute thesis, to be submitted soon to the Open Journal of Astrophysics. Comments welcome

Two notable anomalies in radio observations -- the excess radiation in the Rayleigh-Jeans tail of the cosmic microwave background, revealed by ARCADE2, and the twice-deeper absorption trough of the global 21cm line, identified by EDGES -- remain unresolved. These phenomena may have a shared origin, as the enhancement of the 21cm absorption trough could arise from excess heating. We investigate this scenario through the framework of axion-like particles (ALPs), showing that the resonant conversion of ALPs into photons can produce a photon abundance sufficient to resolve both anomalies simultaneously. Our model naturally explains the observed radio excess between 0.4 and 10GHz while also enhances the 21cm absorption feature at 78MHz. Furthermore, it predicts a novel power-law scaling of the radio spectrum above 0.5GHz and an additional absorption trough below 30MHz, which could be verified through cross-detection in upcoming experiments., Comment: 11 pages, 4 figures

How much cosmological information can we reliably extract from existing and upcoming large-scale structure observations? Many summary statistics fall short in describing the non-Gaussian nature of the late-time Universe in comparison to existing and upcoming measurements. In this article we demonstrate that we can identify optimal summary statistics and that we can link them with existing summary statistics. Using simulation based inference (SBI) with automatic data-compression, we learn summary statistics for galaxy catalogs in the context of cosmological parameter estimation. By construction these summary statistics do not require the ability to write down an explicit likelihood. We demonstrate that they can be used for efficient parameter inference. These summary statistics offer a new avenue for analyzing different simulation models for baryonic physics with respect to their relevance for the resulting cosmological features. The learned summary statistics are low-dimensional, feature the underlying simulation parameters, and are similar across different network architectures. To link our models, we identify the relevant scales associated to our summary statistics (e.g. in the range of modes between $k= 5 - 30 h/\mathrm{Mpc}$) and we are able to match the summary statistics to underlying simulation parameters across various simulation models., Comment: 44 pages, 14 figures

The concentration-virial mass (c-M) relation is a fundamental scaling relation within the standard cold dark matter ($\Lambda$CDM) framework well established in numerical simulations. However, observational constraints of this relation are hampered by the difficulty of characterising the properties of dark matter haloes. Recent comparisons between simulations and observations have suggested a systematic difference of the c-M relation, with higher concentrations in the latter. In this work, we undertake detailed comparisons between simulated galaxies and observations of a sample of strong-lensing galaxies. We explore several factors of the comparison with strong gravitational lensing constraints, including the choice of the generic dark matter density profile, the effect of radial resolution, the reconstruction limits of observed versus simulated mass profiles, and the role of the initial mass function in the derivation of the dark matter parameters. Furthermore, we show the dependence of the c-M relation on reconstruction and model errors through a detailed comparison of real and simulated gravitational lensing systems. An effective reconciliation of simulated and observed c-M relations can be achieved if one considers less strict assumptions on the dark matter profile, for example, by changing the slope of a generic NFW profile or focusing on rather extreme combinations of stellar-to-dark matter distributions. A minor effect is inherent to the applied method: fits to the NFW profile on a less well-constrained inner mass profile yield slightly higher concentrations and lower virial masses., Comment: 9 pages, 9 figures, 2 tables. Accepted for publication in A&A

We perform a state-of-the-art study of the cosmological phase transitions of the real-scalar extended Standard Model. We carry out a broad scan of the parameter space of this model at next-to-next-to-leading order in powers of couplings. We use effective field theory to account for the necessary higher-order resummations, and to construct consistent real and gauge-invariant gravitational wave predictions. Our results provide a comprehensive account of the convergence of perturbative predictions for the gravitational wave signals in this model. For the majority of the parameter points in our study, we observe apparent convergence. While leading and next-to-leading order predictions of the gravitational wave amplitude typically suffer from relative errors between $10$ and $10^4$, at next-to-next-to-leading order the typical relative errors are reduced to between $0.5$ and $50$. Nevertheless, for those parameter points predicting the largest signals, potentially observable by future gravitational wave observatories, the validity of the perturbative expansion is in doubt., Comment: 25 pages, 3 figures

Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of model fitting. Here we follow an alternative approach based on the spectral variance in a model-independent way. We build upon the work presented in Sharbaf et al. (2023) that studied the covariance of high quality SDSS continuum-subtracted spectra in a relatively narrow range of velocity dispersion (100-150 km/s). Here the same analysis is applied to synthetic data from the EAGLE and Illustris TNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral variance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age-metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with AGN feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasise that this methodology only relies on the spectral variance to assess whether simulations provide a true representation of galaxy formation., Comment: Submitted to MNRAS. 27 pages including supplementary material. Comments welcome

The growth of ionized hydrogen bubbles in the intergalactic medium around early luminous objects is a fundamental process during the Epoch of Reionization (EoR). In this study, we analyze bubble sizes and their evolution using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation throughout a large (95.5 cMpc)^3 volume of the Universe. Analogous to the accretion and merger tree histories employed in galaxy formation simulations, we characterize the growth and merger rates of ionized bubbles by focusing on the spatially-resolved redshift of reionization. By tracing the chronological expansion of bubbles, we partition the simulation volume and construct a natural ionization history. We identify three distinct stages of ionized bubble growth: (1) initial slow expansion around the earliest ionizing sources seeding formation sites, (2) accelerated growth through percolation as bubbles begin to merge, and (3) rapid expansion dominated by the largest bubble. Notably, we find that the largest bubble emerges by z=9-10, well before the midpoint of reionization. This bubble becomes dominant during the second growth stage, and defines the third stage by rapidly expanding to eventually encompass the remainder of the simulation volume and becoming one of the few bubbles actively growing. Additionally, we observe a sharp decline in the number of bubbles with radii around ~10 cMpc compared to smaller sizes, indicating a characteristic scale in the final segmented bubble size distribution. Overall, these chronologically sequenced spatial reconstructions offer crucial insights into the physical mechanisms driving ionized bubble growth during the EoR and provide a framework for interpreting the structure and evolution of reionization itself., Comment: 17 pages, 23 figures. Comments welcome. Please visit www.thesan-project.com for more details

We show that the anisotropies in the spectrum of gravitational waves induced by scalar modes after the end of inflation in canonical, single-field models are completely determined by the tilt of the scalar and tensor power spectra. The latter contains information about anisotropies produced due to the propagation of the tensor modes in an inhomogeneous Universe, whereas the former represents the anisotropies generated at the time of production and arise only when non-Gaussian corrections to the angular power spectrum are considered. Our proof takes into account all scalar interactions in the cubic inflaton Lagrangian., Comment: 8 pages

Stochastic inflation, together with the $\Delta N$ formalism, provides a powerful tool for estimating the large-scale behaviour of primordial fluctuations. We construct a numerical code to capture the non-perturbative statistics of such fluctuations and test our code to obtain the exponential non-Gaussian tail of the curvature perturbations. We provide a numerical algorithm to compute the non-perturbative curvature power spectrum and apply it to slow-roll (SR) and ultra-slow-roll (USR) single-field models of inflation. For the USR case, we successfully reproduce the peak in the power spectrum, which agrees to a certain accuracy with the perturbative power spectrum. We highlight some important differences between non-perturbative and perturbative approaches that may suggest the inconsistency of the $\Delta N$ formalism at the transition stages between attractor and non-attractor regimes., Comment: 25 pages, 9 figures

We investigate the effects of ultralight axions (ULAs) on the differential brightness temperature fluctuations of neutral hydrogen at the post-reionization stage. Unique structure suppression features under the influence of ULAs are studied from angular correlations of the neutral hydrogen signals observed at different frequencies. Moreover, ULAs can behave like dark energy or dark matter at different mass regimes, implying that the fraction of the ULA energy density would be correlated with the sum of neutrino masses and the parameters of the dark energy equation of state for dark matter and dark energy like ULAs, respectively. We have explored the parameter space for the ULA physics and possible degeneracies among dark energy, neutrinos, and ULAs using the theoretical angular correlation functions expected from future intensity mapping surveys., Comment: 13 pages, 13 figures

Non-perturbative gravitational effects induce explicit global symmetry breaking terms within axion models. These exponentially suppressed terms in the potential give a mass contribution to the axion-like particles (ALPs). In this work we investigate this scenario with a scalar field charged under a global $U(1)$ symmetry and having a non-minimal coupling to gravity. Given the exponential dependence, the ALP can retain a mass spanning a wide range, which can act as a dark matter component. We specify pre-inflationary and post-inflationary production mechanisms of these ALPs, with the former from the misalignment mechanism and the latter from both the misalignment and cosmic-string decay. We identify the allowed parameter ranges that explain the dark matter abundance for both a general inflation case and a case where the radial mode scalar drives inflation, each in metric and Palatini formalisms. We show that the ALP can be the dominant component of the dark matter in a wide range of its mass, $m_{a} \in [10^{-21}~\mathrm{eV},\, \mathrm{TeV}]$, depending on the inflationary scenario and the $U(1)$ breaking scale. These results indicate that ALPs can be responsible for our dark matter abundance within a setup purely from non-perturbative gravitational effects., Comment: 25 pages, 5 figures, 1 table

The standard cosmological model is sufficiently well constrained that precise estimates can be provided for the redshift of various physically defined times in the chronology of the Universe. For example, it is well known that matter-radiation equality, recombination and reionisation happen at redshifts of around 3000, 1000 and 10, respectively, and these can be specified more precisely by fitting to data. What is less well known are the times in years (and their uncertainties) for these and other epochs in the history of the Universe. Here we provide precise time determinations for six epochs in cosmological history within the standard model, using data from the Planck satellite. Our main results are illustrated in a figure., Comment: 11 pages, 2 figures, 2 tables, data available at https://doi.org/10.5281/zenodo.14054334

We consider a scalar field theory with a Minkowski false vacuum and an unbounded (or very deep) true vacuum. We show compelling evidence that an AdS bubble of vanishing total energy, embedded in asymptotically flat spacetime, generically undergoes a spherical collapse which leads to a space-like curvature singularity after the formation of trapped surfaces and apparent horizons. The crunch singularity, which is hided behind an apparent horizon, occurs before the true vacuum is reached, and the existence of a lower bound of the scalar field potential is not a necessary condition for its formation., Comment: 7+6 pages, 7+7 figures

We investigate the formation of primordial black hole (PBH) based on numerical relativity simulations and peak theory as well as the corresponding scalar induced gravitational wave (SIGW) signals in the presence of \emph{logarithmic non-Gaussianities} which has recently been confirmed in a wide class of inflation models. Through numerical calculations, we find certain parameter spaces of the critical thresholds for the type A PBH formation and reveal a maximum critical threshold value. We also find that there is a region where no PBH is produced from type II fluctuations contrary to a previous study. We then confirm that SIGW signals originated from the logarithmic non-Gaussianity are detectable in the Laser Interferometer Space Antenna if PBHs account for whole dark matter. Finally, we discuss the SIGW interpretation of the nHz stochastic gravitational wave background reported by the recent pulsar timing array observations. We find that PBH overproduction is a serious problem for most of the parameter space, while this tension might still be alleviated in the non-perturbative regime., Comment: 20 pages, 11 figures

This study investigates the formation of primordial black holes (PBHs) resulting from the collapse of adiabatic fluctuations with large amplitudes and non-Gaussianity. Ref. \cite{Uehara:2024yyp} showed that fluctuations with large amplitudes lead to the formation of type B PBHs, characterized by the existence of the bifurcating trapping horizons, distinct from the more common type A PBHs without a bifurcating trapping horizon. We focus on the local type non-Gaussianity characterized by the curvature perturbation $\zeta$ given by a function of a Gaussian random variable $\zeta_{\rm G}$ as $\beta\zeta=-\ln(1-\beta \zeta_{\rm G})$ with a parameter $\beta$. Then we examine how the non-Gaussianity influences the dynamics and the type of PBH formed, particularly focusing on type II fluctuations, where the areal radius varies non-monotonically with the coordinate radius. Our findings indicate that, for $\beta>-2$, the threshold for distinguishing between type A and type B PBHs decreases with increasing $\beta$ similarly to the threshold for black hole formation. Additionally, for large positive values of $\beta$, the threshold for type B PBHs approaches that for type II fluctuations. We also find that, for a sufficiently large negative value of $\beta\lesssim-4.0$, the threshold value is in the type II region of $\mu$, i.e., there are fluctuations of type II that do not form black holes. Lastly, we calculate the PBH mass for several values of $\beta$. Then we observe that the final mass monotonically increases with the initial amplitude within the parameter region of type A PBHs, which differs from previous analytical expectations., Comment: 20pages, 14figures

We investigate the global structure of the recently discovered family of $SL(2,\mathbb{Z})$-invariant potentials describing inflationary $\alpha$-attractors. These potentials have an inflationary plateau consisting of the fundamental domain and its images fully covering the upper part of the Poincar\'e half-plane. Meanwhile, the lower part of the half-plane is covered by an infinitely large number of ridges, which, at first glance, are too sharp to support inflation. However, we show that this apparent sharpness is just an illusion created by hyperbolic geometry, and each of these ridges is physically equivalent to the inflationary plateau in the upper part of the Poincar\'e half-plane., Comment: 32 pages 17 Figures

The cosmic velocity field is an unbiased probe of the total matter distribution but is challenging to measure directly at intermediate and high redshifts. The large-scale velocity field imprints a signal in the cosmic microwave background (CMB) through the kinetic Sunyaev-Zeldovich (kSZ) effect. We perform the first 3d reconstruction of the large-scale velocity field from the kSZ effect by applying a quadratic estimator to CMB temperature maps and the 3d positions of galaxies. We do so by combining CMB data from the fifth data release of the Atacama Cosmology Telescope (in combination with Planck) and a spectroscopic galaxy sample from the Sloan Digital Sky Survey. We then measure the galaxy-velocity cross-power spectrum and detect the presence of the kSZ signal at a signal-to-noise ratio of 7.2$\sigma$. Using this galaxy-velocity cross-correlation alone, we constrain the amplitude of local primordial non-Gaussianity finding $f_{\rm NL}=-90^{+210}_{-350}$. This pathfinder measurement sets the stage for joint galaxy-CMB kSZ constraints to significantly enhance the $f_{\rm NL}$ information obtained from galaxy surveys through sample variance cancellation.

The standard $\Lambda$CDM cosmological model informed by cosmic microwave background (CMB) anisotropies makes a precise prediction for the growth of matter density fluctuations over cosmic time on linear scales. A variety of cosmological observables offer independent and complementary ways of testing this prediction, but results have been mixed, with many constraints on the amplitude of structure $S_8$ being 2-3$\sigma$ lower than the expectation from Planck primary CMB anisotropies. It is currently unclear whether these discrepancies are due to observational systematics, non-linearities and baryonic effects or new physics. We review how gravitational lensing of the CMB has and will continue to provide insights into this problem, including through tomographic cross-correlations with galaxy surveys over cosmic time., Comment: This is a review based on a talk I was invited to give at the Royal Society Meeting "Challenging the standard cosmological model" (London, April 2024). This matches the version accepted for publication in "Philosophical Transactions of the Royal Society A" (Accepted Oct 2024). 35 pages, 4 figures

In 2018 the EDGES experiment claimed the first detection of the global cosmic 21cm signal, which featured an absorption trough centered around $z \sim 17$ with a depth of approximately -500mK. This amplitude is deeper than the standard prediction (in which the radio background is determined by the cosmic microwave background) by a factor of two and potentially hints at the existence of a radio background excess. While this result was obtained by fitting the data with a phenomenological flattened-Gaussian shape for the cosmological signal, here we develop a physical model for the inhomogeneous radio background sourced by the first galaxies hosting population III stars. Star formation in these galaxies is quenched at lower redshifts due to various feedback mechanisms, so they serve as a natural candidate for the excess radio background hinted by EDGES, without violating present day measurements by ARCADE2. We forward-model the EDGES sky temperature data, jointly sampling our physical model for the cosmic signal, a foreground model, and residual calibration errors. We compare the Bayesian evidences obtained by varying the complexity and prior ranges for the systematics. We find that the data is best explained by a model with seven log-polynomial foreground terms, and that it requires calibration residuals. Interestingly, the presence of a cosmic 21cm signal with a non-standard depth is decisively disfavored. This is contrary to previous EDGES analysis in the context of extra radio background models, serving as a caution against using a ''pseudo-likelihood'' built on a model (flattened Gaussian) that is different from the one being used for inference. We make our simulation code and associated emulator publicly-available., Comment: 17 pages, 9 figures, 1 table, version submitted to A&A

In the present work, evaporation of a black hole immersed in a de Sitter environment is considered. Vaidya-de Sitter spacetime is used to model the process in a scenario of accelerated expansion of the Universe. The role of observers is highlighted in the development and Hayward thermodynamics for non stationary geometries is employed in the description of the compact objects. The results of the proposed dynamical model are compared with the usual description based on stationary geometries, focusing on primordial black holes (PBHs). It is found how the timescale of evaporation depends on the choice of a cosmological observer. It may differ substantially from the treatment based on stationary models for black holes. In particular, the standard assertion that there is a fixed initial mass just below $10^{15} \, \text{g} \sim 10^{-18} M_\odot$ for the PBHs which are ending their evaporation process today is imprecise, even when possible quantum corrections at the late stages are not considered. Deviations from this prediction appear when the evaporation is measured with respect to the cosmological time., Comment: 16 pages, 10 figures