2025-04-03 | | Total: 25
Halo bias links the statistical properties of the spatial distribution of dark matter halos to those of the underlying dark matter field, providing insights into clustering properties in both general relativity (GR) and modified-gravity scenarios such as f(R) models. While the primary halo mass-dependent bias has been studied in detailed, the secondary bias, which accounts for the additional dependencies on other internal halo properties, can offer a sensitive probe for testing gravity beyond the ΛCDM model. To quantify any potential deviations between ΛCDM and f(R) gravity models in halo clustering, at both the primary and secondary level, as well as in the distributions of halo properties in the cosmic web. Using N-body simulations of f(R) gravity models, we assess the scaling relations and the primary and secondary bias signals of halo populations on the basis of a halo-by-halo estimator of large-scale effective bias. Our analysis is performed using halo number density as the independent variable. The relative difference in the effective bias between the f(R) models and ΛCDM is sensitive, albeit slightly, to the power index of modified gravity. The largest deviations from GR are measured for low-mass halos, where the average bias at fixed number density decreases by up to 5\% for fixed scaling indices. We also show that the scaling relations for some environmental properties, including neighbour statistics, Mach number and local overdensity, exhibit small but non-negligible deviations (~3-5\%) from GR for a wide range of number densities. Our results also suggests that the properties of halos in sheets and voids show the largest departures from GR (> 10\% in some cases). In terms of secondary bias, we do not find any statistically significant deviations with respect to ΛCDM for any of the properties explored in this work.
Wave or fuzzy dark matter produced with high momenta behaves in many ways like hot particle dark matter while also possessing seemingly different phenomenology due to wave interference. We develop wave perturbation theory to show that white noise density fluctuations generated by the interference of high-momenta waves are gravitationally unstable in the usual way during matter domination above the free streaming scale and stabilize below the free streaming scale, much like the analogous effects for massive neutrinos in hot dark matter. We verify and illustrate these effects in the density power spectra of Newtonian Schrödinger-Poisson simulations. In the cosmological context, this would cause a gradual suppression of the initial white noise isocurvature perturbations below the free streaming scale at matter radiation equality, unlike cold dark matter isocurvature fluctuations, and virial stability of dark matter halos.
The distance duality relation (DDR) relates two independent ways of measuring cosmological distances, namely the angular diameter distance and the luminosity distance. These can be measured with baryon acoustic oscillations (BAO) and Type Ia supernovae (SNe Ia), respectively. Here, we use recent DESI DR1, Pantheon+, SH0ES and DES-SN5YR data to test this fundamental relation. We employ a parametrised approach and also use model-independent Generic Algorithms (GA), which are a machine learning method where functions evolve loosely based on biological evolution. When we use DESI and Pantheon+ data without Cepheid calibration or big bang nucleosynthesis (BBN), there is a 2σ violation of the DDR in the parametrised approach. Then, we add high-redshift BBN data and the low-redshift SH0ES Cepheid calibration. This reflects the Hubble tension since both data sets are in tension in the standard cosmological model ΛCDM. In this case, we find a significant violation of the DDR in the parametrised case at 6σ. Replacing the Pantheon+ SNe Ia data by DES-SN5YR, we find similar results. For the model-independent approach, we find no deviation in the uncalibrated case and a small deviation with BBN and Cepheids which remains at 1σ. This shows the importance of considering model-independent approaches for the DDR.
The standard cosmological model, ΛCDM, assumes isotropy on large cosmic scales. However, recent studies using galaxy cluster scaling relations reported an apparent H0 anisotropy at 5.4σ that could be attributed to large bulk flows extending beyond 500Mpc, in disagreement with ΛCDM. To quantify the statistical tension of the observational galaxy cluster data used in past studies with ΛCDM, we utilise the isotropic (2.8Gpc)3 run of the FLAMINGO (ΛCDM) simulations, the largest hydrodynamical cosmological simulation available to date. We create 1728 simulated lightcones and study the apparent level of anisotropy traced by X-ray and thermal Sunyaev-Zeldovich scaling relations in the same cluster sample selection and methodology as in Migkas et al. (2021, arXiv:2103.13904). We find the probability of such apparent anisotropies randomly emerging in cluster scaling relations within a ΛCDM universe to be 0.12%(3.2σ). The discrepancy goes up to ∼3.6σ when modelled as a bulk flow at z<0.1. We find that statistical noise accounts for over 80% of the anisotropy amplitude in each lightcone, with large peculiar velocities contributing less than 20%. We also show that anisotropy amplitudes are highly sensitive to the intrinsic scatter in the scaling relations, with tighter relations providing stronger constraints. Nevertheless, the tension between Migkas et al. (2021, arXiv:2103.13904) and ΛCDM persists, however, at a lower significance than previously reported.
We study the dynamics of Friedmann-Lemaître-Robertson-Walker models where a dark energy component with a quadratic equation of state (EoS) nonlinearly interacts with cold dark matter. Thus, two energy scales naturally come into play: ρ∗ is the scale at which the nonlinearity of the EoS becomes relevant; ρi is the energy scale around which the interaction starts to play an important role in the dynamics. Our focus is to understand whether there are parameter ranges for this system that can produce non-singular bouncing and emergent cosmologies for any initial condition. We complete a dynamical systems analysis, and find the parameter range such that trajectories always expand from a high energy non-singular de Sitter state. For flat and negative curvature models this de Sitter state is represented by a fixed point, the asymptotic past from which the universe emerges from. We find a subset of positive curvature models that during contraction get arbitrarily close to the de Sitter state, thus having a quasi-de Sitter bounce, then emerge from the bounce and expand, evolving toward spatial flatness. We find that the dimensionless parameter q≡ρ∗/ρi, which measures the relative strength of the nonlinear terms in the system, plays a crucial role in the topology of the phase space. When q<3, some trajectories expand toward a singularity, while others evolve toward a low energy cosmological constant at late-times, with a subset going through a decelerated matter dominate era before the final acceleration. When q>3, all trajectories are non-singular, and evolve toward a late-time cosmological constant. We find a subclass in this case in which all trajectories have at least one decelerated matter dominated phase, and accelerate at late-times.
The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-time inflationary physics, and underpinned by gravitation through general relativity. There have always been open questions about the soundness of the foundations of the standard model. However, recent years have shown that there may also be questions from the observational sector with the emergence of differences between certain cosmological probes. In this White Paper, we identify the key objectives that need to be addressed over the coming decade together with the core science projects that aim to meet these challenges. These discordances primarily rest on the divergence in the measurement of core cosmological parameters with varying levels of statistical confidence. These possible statistical tensions may be partially accounted for by systematics in various measurements or cosmological probes but there is also a growing indication of potential new physics beyond the standard model. After reviewing the principal probes used in the measurement of cosmological parameters, as well as potential systematics, we discuss the most promising array of potential new physics that may be observable in upcoming surveys. We also discuss the growing set of novel data analysis approaches that go beyond traditional methods to test physical models. [Abridged]
Time-delay distance measurements from strongly lensed quasars provide a robust, independent method for determining the Hubble constant (H0). This approach cross-checks H0 estimates from the distance ladder in the late universe and the cosmic microwave background in the early universe. However, the mass-sheet degeneracy in lensing models introduces systematic uncertainty, limiting precision. Dynamical modeling complements strong lensing by constraining the mass distribution with independent observational data. We develop a methodology and software framework for joint modeling of stellar kinematics and lensing data. Using simulated data for the lensed quasar RXJ1131−1131, we demonstrate that high-quality kinematic data can achieve ∼4% precision on H0. Through extensive modeling, we examine the impact of the presence of a supermassive black hole in the lens galaxy and potential systematic biases in kinematic data on H0 measurements. Our results show that imposing priors on black hole mass and orbital anisotropy, or excluding central kinematic bins, mitigates biases in H0 estimates. By testing on mock kinematic data with systematic biases, we highlight the need for sub-percent control of kinematic systematics, which is achievable with current technology. Additionally, we leverage GPU parallelization to accelerate Bayesian inference, reducing a previously month-long process by an order of magnitude. This pipeline offers significant potential for advancing cosmological and galaxy evolution studies with large datasets.
Quasars are effective tracers of the large-scale distribution of galaxies at high redshift thanks to their high luminosity and dedicated surveys. Previous studies have shown that quasars exhibit a bias similar to that of rich groups, indicating that quasar pairs may be linked to higher density environments, serving as protocluster proxies. In this work, we aim to characterize close quasar pairs within the same halo by identifying them in redshift space. We analyze pair-quasar cross correlations and CMB-derived lensing convergence profiles centered in these systems. We have identified 2777 pairs of quasars in the redshift range from 1.2 to 2.8. Quasar pairs exhibit a distribution of relative luminosities that differs from that of random pairs with the same redshift distribution, indicating that quasars in these systems are distinct from isolated ones. The cross-correlation between pairs and quasars shows a larger amplitude than the auto-correlation function of quasars, suggesting these systems are more strongly biased toward the large-scale mass distribution and reside in more massive halos. This is further supported by higher convergence CMB lensing profiles of pairs compared to isolated quasars with a similar redshift distribution. Our results indicate that quasar pairs are suitable precursors to present-day clusters of galaxies, unlike isolated quasars, which are associated with less dense environments.
Galaxy groups are essential for studying the distribution of matter on a large scale in redshift surveys and for deciphering the link between galaxy traits and their associated halos. In this work, we propose a widely applicable method for identifying groups through machine learning techniques in real space taking into account the impact of redshift distortion. Our methodology involves two neural networks: one is a classification model for identifying central galaxy groups, and the other is a regression model for predicting the mass of these groups. Both models input observable galaxy traits, allowing future applicability to real survey data. Testing on simulated datasets indicates our method accurately identifies over 92% of groups with Mvir≥1011h−1M⊙, with 80% achieving a membership completeness of at least 80%. The predicted group masses vary by less than 0.3 dex across different mass scales, even in the absence of a priori data. Our network adapts seamlessly to expand to sparse samples with a flux limit of mr<14, to high redshift samples at z=1.08, and to galaxy samples from the TNG300 hydrodynamical simulation without further training. Furthermore, the framework can easily adjust to real surveys by training on redshift distorted samples without needing parameter changes. Careful consideration of different observational effects in redshift space makes it promising that this method will be applicable to real galaxy surveys.
We present a hybrid method for reconstructing the primordial density from late-time halos and galaxies. Our approach involves two steps: (1) apply standard Baryon Acoustic Oscillation (BAO) reconstruction to recover the large-scale features in the primordial density field and (2) train a deep learning model to learn small-scale corrections on partitioned subgrids of the full volume. At inference, this correction is then convolved across the full survey volume, enabling scaling to large survey volumes. We train our method on both mock halo catalogs and mock galaxy catalogs in both configuration and redshift space from the Quijote 1(h−1Gpc)3 simulation suite. When evaluated on held-out simulations, our combined approach significantly improves the reconstruction cross-correlation coefficient with the true initial density field and remains robust to moderate model misspecification. Additionally, we show that models trained on 1(h−1Gpc)3 can be applied to larger boxes--e.g., (3h−1Gpc)3--without retraining. Finally, we perform a Fisher analysis on our method's recovery of the BAO peak, and find that it significantly improves the error on the acoustic scale relative to standard BAO reconstruction. Ultimately, this method robustly captures nonlinearities and bias without sacrificing large-scale accuracy, and its flexibility to handle arbitrarily large volumes without escalating computational requirements makes it especially promising for large-volume surveys like DESI.
Bootstrap techniques relying on the constraints imposed by Extended Galilean Invariance (EGI), have proved to be very useful in the context of perturbation theory of the Large Scale Structure (LSS). It has been formulated in both the Eulerian as well as Lagrangian space. While the Eulerian bootstrap formalism has been successfully applied to both tracer and matter kernels, the application of bootstrap methods in Lagrangian space has so far been restricted to matter. Up to third order, it has been shown that implementing EGI constraints in Eulerian space fully reproduces the bias expansion for tracers. Previous studies have demonstrated that time non-locality affects the bias expansion in a non-trivial way starting from fifth order. Motivated by this fact, we extend the bootstrap approach upto fifth order in both Eulerian and Lagrangian space and demonstrate that it fully captures the time non-local effects. For this, we generalize the Lagrangian bootstrap for tracers, and found that it agrees with the corresponding results obtained in Eulerian space. One of the major challenges in implementing EGI constraints in Eulerian space, is to systematically find out all the "spurious poles" and make them vanish. We have proposed a method that bypasses this difficulty making the procedure tractable at higher orders. From Lagrangian perspective, we have identified coefficients in the tracer kernel whose ratios are independent of tracer properties and may serve as direct probes of the underlying cosmology.
We present the weak-lensing mass calibration of 124 galaxy clusters and groups at redshift 0.1<z<0.8 in the first Data Release of the eROSITA All-Sky Survey (eRASS1) and constrain their BCG stellar-mass-to-halo-mass-and-redshift relation, using the data sets from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The cluster survey is conducted by the eROSITA X-ray telescope aboard the Spectrum-Roentgen-Gamma space observatory. The cluster sample is X-ray-selected and optically confirmed with a negligibly low contamination rate. On a basis of individual clusters, the shear profiles of 96 clusters are derived using the HSC Three-Year weak-lensing data, while the BCG stellar masses M⋆,BCG of 101 clusters are estimated using the SED template fitting to the HSC five-band photometry. The count rate (CR−M−z) and BCG stellar mass (M⋆,BCG−M−z) relations are simultaneously constrained in a population modelling. The CR−M−z relation reveals a mass trend of ∝M1.50+0.20−0.30, which is steeper than the self-similar prediction, and no deviation from the self-similar redshift scaling. We obtain the power-law indices of the mass (BBCG=0.25+0.11−0.15) and redshift (γBCG=0.87+0.92−0.78) scaling for the M⋆,BCG−M−z relation, and find a strong degeneracy between them. By adopting an informative prior on γBCG to break the degeneracy, we obtain a statistically consistent M⋆,BCG−M−z relation. Our results suggest that the BCG stellar mass at a fixed halo mass has remained stable with a moderate increase since redshift z≈0.8. This finding supports the picture of the ``rapid-then-slow'' BCG formation, where the majority of the stellar mass must have been assembled at much earlier cosmic time.
Galaxy and halo scaling relations, connecting a broad range of parameters, are well established from observations. The origin of many of these relations and their scatter is still a matter of debate. It remains a sizable challenge for models to simultaneously and self-consistently reproduce as many scaling relations as possible. We introduce the Magneticum Pathfinder hydrodynamical cosmological simulation suite, to date the suite that self-consistently covers the largest range in box volumes and resolutions. It is the only cosmological simulation suite that is tuned on the hot gas content of galaxy clusters instead of the stellar mass function. By assessing the successes and shortcomings of tuning to the hot gas component of galaxy clusters, we aim to further our understanding of the physical processes shaping the Universe. We analyze the importance of the hot and cold gas components for galaxy and structure evolution. We analyze 28 scaling relations, covering large-scale global parameters as well as internal properties for halos ranging from massive galaxy clusters down to galaxies, and show their predicted evolution from z=4 to z=0 in comparison with observations. These include the halo-to-stellar-mass and Kennicutt--Schmidt relations, the cosmic star formation rate density as well as the Fundamental Plane. Magneticum Pathfinder matches a remarkable number of the observed scaling relations from z=4 to z=0, including challenging relations like the number density of quiescent galaxies at cosmic dawn, the mass--size evolution, the mass--metallicity relation, the Magorrian relation, and the temperature--mass relation. We compile our data to allow for straightforward future comparisons. Galaxy properties and scaling relations arise naturally and the large scatter in observables at high redshift is crucial to distinguish the various galaxy formation models reproducing the z=0 relations.
We discuss a minimal renormalizable Pati-Salam theory based on the SU(4)C×SU(2)L×SU(2)R gauge group, with unification scale Higgs multiplets taken as SU(2)L and SU(2)R doublets, which lead to neutrino Dirac picture. Although a number of scalar particles could be light, even lying at the LHC energies, the unification scale is hopelessly out of reach in any foreseeable future. Moreover, phase transition in the early Universe leads to the production of magnetic monopoles and domain walls, both incompatible with the standard cosmological model. A small explicit breaking of the discrete left-right symmetry allows the domain walls to decay, and in the process possibly sweep away the monopoles, analogously to the previously discussed case of SU(5) grand unified theory. This leaves an important imprint of gravitational waves, within the reach of next generation searches, correlated with monopole detection and new light particles at collider energies. The theory has a dark matter candidate in the form of an inert scalar doublet, with a mass below TeV, which can further trigger electroweak baryogenesis.
We present a comprehensive analysis of the MIRI Extremely Red Object Virgil, a Lyman-α emitter at zspec=6.6379±0.0035 with the photometric properties of a Little Red Dot. Leveraging new JWST/MIRI imaging from the MIDIS and PAHSPECS programs, we confirm Virgil's extraordinary nature among galaxies in JADES/GOODS-South, exhibiting a strikingly red NIRCam-to-MIRI color (F444W − F1500W = 2.84±0.04~mag). Deep NIRSpec/PRISM spectroscopy from the OASIS program offers key insights into the host galaxy, revealing properties of an average star-forming galaxy during Cosmic Reionization, such as a subsolar metallicity, low-to-moderate dust content, and a relatively high ionization parameter and electron temperature. By estimating the star formation rate of Virgil from UV and Hα, we find evidence that the galaxy is either entering or fading out of a bursty episode. Although line-ratio diagnostics employed at high-z would classify Virgil as an Active Galactic Nucleus (AGN), this classification becomes ambiguous once redshift evolution is considered. Nonetheless, Virgil occupies the same parameter space as recently confirmed AGNs at similar redshifts. The new deep MIRI data at 15~μm reinforce the AGN nature of Virgil, as inferred from multiple spectral energy distribution (SED) fitting codes. Virgil's rising infrared SED and UV excess resemble those of Dust-Obscured Galaxies (DOGs) studied with Spitzer at Cosmic Noon, particularly blue-excess HotDOGs. Our results highlight the need for a multi-wavelength approach incorporating MIRI to uncover such extreme sources at z≳ and to shed light on the interplay between galaxy evolution and early black hole growth during Cosmic Reionization.
The Brans-Dicke scalar-tensor cosmological models are studied in both Einstein and Jordan frames, using hydrodynamical and self-interacting scalar field representations of the energy-momentum tensor, leading to the same background solutions. The main features of the corresponding cosmological scenarios are determined. In many cases, we identify a transition from decelerated to accelerated regimes. The properties of the self-interacting scalar field, including the corresponding potential, are determined. Finally, some possible realistic configurations are analyzed.
Within the framework of general U(1) scenario, we demonstrate that the ultra high energy neutrinos recently detected by KM3NeT could originate from a decaying right handed neutrino dark matter (DM), with a mass of 440 PeV. Considering DM production via freeze-in, we delineate the parameter space that satisfies the observed relic abundance and also lies within the reach of multiple gravitational wave detectors. Our study provides a testable new physics scenario, enabled by multi-messenger astronomy.
If a coupling between the inflaton and the Gauss-Bonnet term is introduced, many models of inflation that were ruled out by the most recent Planck data can be made viable again. The predictions for the scalar spectral index and tensor-to-scalar ratio are typically computed using the slow-roll approximation. In this paper we instead study the full equations of motion and determine the necessary initial conditions for reasonable inflation epoch. We derive the conditions under which the Friedmann equation admits positive solutions for the Hubble parameter. Then we study the possibility of the inflaton becoming trapped in a local potential minimum induced by the Gauss-Bonnet term. Finally we demonstrate the results on monomial potential models with a quadratic and a quartic potential and show that the slow-roll approximation becomes imprecise in the quartic case.
We investigate the influence of the cosmic web on galaxy properties in the IllustrisTNG simulations. To disentangle the effects of galaxy groups and cosmic filaments, we classify the cosmic web environment into four categories: group, group-dominated, filament-dominated, and field. By controlling for stellar mass, we reveal evident differences in specific star formation rates (sSFR), quenched fraction, gas fractions, local density, and stellar ages among central galaxies in different cosmic web environments, particularly for lower-mass galaxies. However, these differences largely diminish when the effect of local overdensity is further accounted for, indicating its dominant role. Additionally, we observe distinct differences in these properties among satellite galaxies across environments, mainly driven by stellar mass, halo mass, and overdensity. Notably, residual differences between satellites in field and filament-dominated region persist even after controlling for these factors, suggesting a stronger susceptibility of satellite galaxies to filaments compared to centrals. Our findings highlight the importance of differentiating between central and satellite to accurately assess the environmental effects of the cosmic web. Our analysis suggests that the relationship between galaxy properties and their distance from filaments arises from a combination of factors, including stellar and halo mass, groups, overdensity, and the intrinsic influence of the cosmic web. Additionally, we find that the effect of the cosmic web on galaxy properties is reduced at z=0.5, compared to z=0. Furthermore, central galaxies near thick filaments tend to exhibit slightly to moderately lower sSFR and cold gas fractions compared to those near thin filaments.
A perturbative strategy for inflation described by two-inflaton fields is developed using a mathematical analogy with the renormalization-group. Two small quantities, \alpha and \lambda, corresponding to standard slow-roll parameters are defined and systematic expansions of all inflationary quantities in terms of powers of \alpha and \lambda are found. No other slow-roll parameters are needed. To illustrate this perturbative method in the multi-field context, we adopt a simple two-inflaton model with quadratic potentials in the parameter range where both fields contribute similarly to the dynamics of inflation. The model, even though it is not a viable alternative for phenomenological description of the inflationary period, nicely illustrates subtleties of the perturbative approach. In particular, this method allows us to derive two independent gauge-invariant scalar perturbations that are conserved in the superhorizon limit, overcoming typical problems that emerge in multi-field inflation. Furthermore, it is possible to perform nonperturbative resummations that allow to study the model in a true multi-field regime. We derive tensor and scalar power spectra to the next-to-next-to-leading and next-to-leading orders, respectively, as well as their spectral indices. The hierarchy between the two scalar perturbations allows us to single out the dominant entry of the scalar power-spectrum matrix. Modifications due to the second inflaton occur already at the leading order. Finally, we explain why the quadratic two-inflaton model is not compatible with the present experimental constraints even though non-trivial corrections to scalar perturbations do emerge.
Steven Weinberg productive scientific life teaches us many things, one of the most important of which is the power of his example. This essay contains personal reminiscences and a speculation about gravitational wave propagation based on one of his very last papers -- as seems fitting, given his penchant for putting interesting physics into the essays he wrote for other luminaries over the years. (See arXiv:2502.10979 [physics.hist-ph] for a less personal summary of his main scientific accomplishments.)
Giant bulgeless disk galaxies, theoretically expected to be rare in the early Universe, have been confirmed by the James Webb Space Telescope (JWST) to exist as early as 2 billion years after the Big Bang. These morphologically extreme systems offer valuable insights into the physics of disk formation and the interplay between galaxies and their dark-matter halos. Using cosmological simulations, we identify analogs of such galaxies with stellar masses around 10^{11} M_\odot and half-light radii up to 6 kpc at z \sim 3 and characterize the factors that contribute to their formation. These galaxies form in young cosmic knots, populating host halos of high spin, low concentration, and spherical shapes. They feature dynamically coherent circum-galactic medium, as well as gas-rich, coherent mergers, which preserve their disk morphology and drive their large sizes. Interestingly, all the simulated giant disks harbor a compact, aligned inner disk, marginally resolvable in JWST images with a Sérsic index near unity. These findings highlight the environmental and structural conditions necessary for forming and sustaining giant bulgeless disks and provide a theoretical framework for interpreting JWST observations of extreme disk morphologies in the early Universe.
Through gravitational lensing, galaxy clusters can magnify supernovae (SNe) and create multiple images of the same SN. This enables measurements of cosmological parameters, which will be increasingly important in light of upcoming telescopic surveys. We study the prospects of detecting strongly lensed SNe in cluster fields with the Nancy Grace Roman Space Telescope (Roman)'s High Latitude Time Domain Survey (HLTDS) and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). We employed two approaches: one focusing on known multiply imaged galaxies behind clusters, along with the SN rates specific to those galaxies, and another based on the expected number of lensed SNe exploding in a given volume behind a galaxy cluster. We collected all the clusters in the literature that feature a well-constrained lens model and multiply imaged galaxies behind clusters with high-quality data for the lensed galaxies. This allowed us to determine the SN rate for each galaxy. We provide predictions for 46 clusters visible to the Vera C. Rubin Observatory, as well as for 9 observable by Roman's HLTDS, depending on whether the clusters fall within the survey's observing field. We predict that the number of multiply imaged SNe discovered by LSST in its first three years is 3.95 \pm 0.89 from the first approach or 4.94 \pm 1.02 from the second. For the HLTDS, the expected number of multiply imaged SNe ranges from 0.38 \pm 0.15 to 5.2 \pm 2.2, depending on the specific cluster observed, however, the fields to be targeted remain a matter of discussion. We conclude that LSST offers great prospects for detecting multiply imaged SNe. Our predictions are effectively lower limits, as we only considered the most massive and well-studied clusters. We provide a recommendation for HLTDS observing field selection, namely: either MACS J0553.4-3342 or Abell 1758a should be observed by the survey.
We investigate general relativistic effects on the photon spectrum emitted from decaying (or annihilating) particle dark matter in the halo surrounding a primordial black hole. The spectrum undergoes significant modification due to gravitational redshifts, which induces broadening as a result of the intense gravitational field near the black hole. This characteristic alteration in the photon spectrum presents a unique observational signature. Future observations of such spectral features may provide critical evidence for a mixed dark matter scenario, involving both primordial black holes and particle dark matter.
James Webb Space Telescope (JWST) has opened up a new chapter in infrared astronomy. Besides the discovery and a deeper understanding of various astrophysical sources, JWST can also uncover the non-gravitational nature of dark matter (DM). If DM is QCD axion or an eV-scale Axion-like particle (ALP), it can decay into two photons in the infrared band. This will produce a distinct line signature in the spectroscopic observations made by JWST. Using the latest NIRSpec IFU spectroscopic observations from JWST, we put the strongest bound on the photon coupling for QCD axion/ ALP DM in the mass range between 0.47 and 2.55 eV. In particular, we are able to probe a new mass range for ALP DM between \sim 0.47 eV to 0.78 eV beyond what can be probed by globular cluster observations. We constrain well-motivated and UV complete models of QCD axion and ALP DM, including predictions from some models derived from string theory and/ or various Grand Unification scenarios. Future JWST observations of DM-rich systems with a better understanding of the astrophysical and instrumental backgrounds can thus enable us to potentially discover QCD axion and ALP DM. The datasets used in this work are available at: https://dx.doi.org/10.17909/3e5f-nv69