Search NASA⌕ Search

SEARCH · Search NASA

Results for “Strong gravitational lensing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Projected Cosmological Constraints from Strongly Lensed Supernovae with the Roman Space Telescope

One of the primary mission objectives of the Roman Space Telescope is to investigate the nature of dark energy with a variety of methods. Observations of Type I supernovae (SNe Ia) will be one of the principal anchors of the Roman cosmology program through traditional luminosity distance measurements. This SNe Ia cosmology program can provide another valuable cosmological probe, without altering the strategy of the mission: time delay cosmography with gravitationally lensed supernova (SN). In this work, we forecast lensed SN cosmology constraints with the Roman Space Telescope, while providing useful tools for future work. Using the anticipated characteristics of the Roman SNe Ia survey, we have constructed mock catalogs of expected resolved lensing systems, as well as strongly lensed Type Ia and core-collapse (CC) SN light curves, including microlensing effects. We predict Roman will find ~11 lensed SNe Ia and ~20 CCSNe, depending on the survey strategy. Next, we estimate the time delay precision obtainable with Roman (Ia: ~2 days, CC: ~3 days), and use a Fisher matrix analysis to derive projected constraints on H 0 ,Ω m , and the dark energy equation of state, w, for each SNe Ia survey strategy. Here, a strategy optimized for the discovery of high-redshift SNe Ia is preferred when considering the constraints possible from both SNe Ia and lensed SN cosmology, also delivering ~1.5 times more lensed SNe than other proposed survey strategies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Assessing the sources of reionization: a spectroscopic case study of a 30× lensed galaxy at z ~ 5 with Ly α , C IV, Mg ιι, and [Ne ιιι]

We present a detailed spectroscopic analysis of a galaxy at z ≃ 4.88 that is, by chance, magnified ~30× by gravitational lensing. Only three sources at z ≳ 5 are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width ${{\rm C\, \small {IV}}}\ \lambda 1549\, \mathring{\rm A}$ emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionization. Using ultraviolet (UV) nebular line ratio diagnostics, Very Large Telescope (VLT)/X-shooter observations rule out strong active galactic nuclei (AGN) activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of ${[{\rm Ne\, \small {III}}]}\ \lambda 3870\, \mathring{\rm A}$ and use the [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] line ratio to constrain the ionization parameter and gas-phase metallicity. Closely related to the commonly used [${\rm O\, \small {III}}$]/[${\rm O\, \small {II}}$] ratio, our [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] measurement shows this source is similar to local ‘Green Pea’ galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the fundamental metallicity relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of ${{\rm Mg\, \small {II}}}\ \lambda 2796\, \mathring{\rm A}$, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting ${\rm Mg\, \small {II}}$ absorption. Strong ${\rm Mg\, \small {II}}$ emission has been observed in most z ~ 0 Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionizing radiation. In conclusion, this strongly lensed galaxy, observed just $300\, \mathrm{Myr}$ after reionization ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the James Webb Space Telescope (JWST)/Extremely Large Telescope (ELT) era.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

DESI Strong Lens Foundry. I. HST Observations and Modeling with GIGA-Lens

We present the Dark Energy Spectroscopic Instrument (DESI) Strong Lens Foundry. We discovered ∼3500 new strong gravitational lens candidates in the DESI Legacy Imaging Surveys using residual neural networks (ResNet). We observed a subset (51) of our candidates using the Hubble Space Telescope (HST). Except for one ambiguous case, we have confirmed 50 of the 51 candidates to be strong lenses. We also briefly describe spectroscopic follow-up observations by DESI and Keck NIRES programs. From this very rich data set, a number of studies will be carried out, including evaluating the quality of the ResNet search candidates and lens modeling. In this paper, we present our initial effort in these directions. In particular, as a demonstration, we present the lens model for DESI-165.4754−06.0423, with imaging data from HST, and lens and source redshifts from DESI and Keck NIRES, respectively. In this effort, we have applied a fully forward-modeling Bayesian approach (GIGA-Lens), using multiple GPUs, to a strong lens with HST data, and achieved statistical convergence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Geometry versus growth: Internal consistency of the flat ΛCDM model with KiDS-1000

Here, we carry out a multi-probe self-consistency test of the flat Lambda Cold Dark Matter (ΛCDM) model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of ΛCDM cosmological parameters. This extended model is constrained by a combination of weak gravitational lensing measurements from the Kilo-Degree Survey, galaxy clustering signatures extracted from Sloan Digital Sky Survey campaigns and the Six-Degree Field Galaxy Survey, and the angular baryon acoustic scale and the primordial scalar fluctuation power spectrum measured in Planck cosmic microwave background (CMB) data. For both the weak lensing data set individually and the combined probes, we find strong consistency between the geometry and growth parameters, as well as with the posterior of standard ΛCDM analysis. In the non-split analysis, for which one single set of parameters was used, tension in the amplitude of matter density fluctuations as measured by the parameter S 8 persists at around 3σ, with a 1.5% constraint of S 8 = 0.776$_{–0.008}^{+0.016}$ for the combined probes. We also observe a less significant preference (at least 2σ) for higher values of the Hubble constant, H 0 = 70.5$_{–1.5}^{+0.7}$ km s –1 Mpc –1 , as well as for lower values of the total matter density parameter Ω m = 0.289$_{–0.005}^{+0.007}$ compared to the full Planck analysis. Including the subset of the CMB information in the probe combination enhances these differences rather than alleviate them, which we link to the discrepancy between low and high multipoles in Planck data. Our geometry versus growth analysis does not yet yield clear signs regarding whether the origin of the discrepancies lies in ΛCDM structure growth or expansion history but holds promise as an insightful test for forthcoming, more powerful data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Identification of more than 40 gravitationally magnified stars in a galaxy at redshift 0.725

Strong gravitational magnification enables the detection of faint background sources and allows researchers to resolve their internal structures and even identify individual stars in distant galaxies. Highly magnified individual stars are useful in various applications, including studies of stellar populations in distant galaxies and constraining dark matter structures in the lensing plane. However, these applications have been hampered by the small number of individual stars observed, as typically one or a few stars are identified from each distant galaxy. Here, we report the discovery of more than 40 microlensed stars in a single galaxy behind Abell 370 at redshift of 0.725 (dubbed ‘the Dragon arc’) when the Universe was half of its current age, using James Webb Space Telescope observations with the time-domain technique. These events were found near the expected lensing critical curves, suggesting that these are magnified stars that appear as transients from intracluster stellar microlenses. Through multi-wavelength photometry, we constrained their stellar types and found that many of them are consistent with red giants or supergiants magnified by factors of hundreds. Furthermore, this finding reveals a high occurrence of microlensing events in the Dragon arc and demonstrates that time-domain observations by the James Webb Space Telescope could lead to the possibility of conducting statistical studies of high-redshift stars.

Early universe↗
DOE OSTI2024

Astronomaly at scale: searching for anomalies amongst 4 million galaxies

ABSTRACT Modern astronomical surveys are producing data sets of unprecedented size and richness, increasing the potential for high-impact scientific discovery. This possibility, coupled with the challenge of exploring a large number of sources, has led to the development of novel machine-learning-based anomaly detection approaches, such as astronomaly. For the first time, we test the scalability of astronomaly by applying it to almost 4 million images of galaxies from the Dark Energy Camera Legacy Survey. We use a trained deep learning algorithm to learn useful representations of the images and pass these to the anomaly detection algorithm isolation forest, coupled with astronomaly’s active learning method, to discover interesting sources. We find that data selection criteria have a significant impact on the trade-off between finding rare sources such as strong lenses and introducing artefacts into the data set. We demonstrate that active learning is required to identify the most interesting sources and reduce artefacts, while anomaly detection methods alone are insufficient. Using astronomaly, we find 1635 anomalies among the top 2000 sources in the data set after applying active learning, including eight strong gravitational lens candidates, 1609 galaxy merger candidates, and 18 previously unidentified sources exhibiting highly unusual morphology. Our results show that by leveraging the human–machine interface, astronomaly is able to rapidly identify sources of scientific interest even in large data sets.

Astronomy & Astrophysics↗
DOE OSTI2022

Cosmological analysis of three-dimensional BOSS galaxy clustering and Planck CMB lensing cross correlations via Lagrangian perturbation theory

In this work, we present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.78 and $σ_8$ = 0.743 ± 0.043. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.77 and $σ_8$ = 0.707 ± 0.035, in tension with CMB and cosmic shear constraints. The combination of Ω m and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($\textit{z}$ ≃ 0.38) and at large angular scales, which show a 20% deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Free-floating or Wide-orbit Planet in the Microlensing Event OGLE-2019-BLG-0551

High-cadence observations of the Galactic bulge by the microlensing surveys led to the discovery of a handful of extremely short-timescale microlensing events that can be attributed to free-floating or wide-orbit planets. Here, we report the discovery of another strong free-floating planet candidate, which was found from the analysis of the gravitational microlensing event OGLE-2019-BLG-0551. The light curve of the event is characterized by a very short duration (≲3 days) and a very small amplitude (≲0.1 mag). From modeling of the light curve, we find that the Einstein timescale, t{sub E}=0.381±0.017 day, is much shorter, and the angular Einstein radius, θ{sub E}=4.35±0.34 μas, is much smaller than those of typical lensing events produced by stellar-mass lenses (t{sub E}∼20 days, θ{sub E}∼0.3 mas), indicating that the lens is very likely to be a planetary-mass object. We conduct an extensive search for possible signatures of a companion star in the light curve of the event, finding no significant evidence for the putative host star. For the first time, we also demonstrate that the angular Einstein radius of the lens does not depend on blending in the low-magnification events with strong finite source effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Self-interacting dark matter and small-scale gravitational lenses in galaxy clusters

Recently, Meneghetti et al. reported an excess of small-scale gravitational lenses in galaxy clusters. We study its implications for self-interacting dark matter (SIDM) compared with standard cold dark matter (CDM). We design controlled N-body simulations that incorporate observational constraints. The presence of early-type galaxies in cluster substructures can deepen gravitational potential and reduce tidal mass loss. Both scenarios require a relatively high baryon concentration in the substructure to accommodate the lensing measurements, and their tangential caustics are similar. The SIDM substructure can experience gravothermal collapse and produce a steeper density profile than its CDM counterpart, leading to a larger radial galaxy-galaxy strong-lensing cross section, although this effect is hard to observe. Our results indicate that SIDM can provide a unified explanation to small-scale lenses in galaxy clusters and stellar motion in dwarf galaxies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on ultralight axion dark matter through galaxy cluster number counts

Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range 10 −32 eV to 10 −24 eV. We put upper bounds on the ultralight axion relic density Ω a in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around m a = 10 −27 eV with Ω a < 0.0035 and m a = 10 −26 eV with Ω a < 0.0079; both are at a 95% confidence level. When combined with cosmic microwave background probes, these bounds are tightened to Ω a < 0.0030 in the m a = 10 −27 eV mass bin and Ω a < 0.0058 in the m a = 10 −26 eV mass bin, with both at a 95% confidence level. This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.

cosmological parameters↗
DOE OSTI2021

Stacked CMB lensing and ISW signals around superstructures in the DESI Legacy Survey

The imprints of large-scale structures on the Cosmic Microwave Background (CMB) can be studied via the CMB lensing and Integrated Sachs–Wolfe (ISW) signals. In particular, the stacked ISW signal around supervoids has been claimed in several works to be anomalously high. In this study, we find cluster and void superstructures using four tomographic redshift bins with 0 < z < 0.8 from the DESI Legacy Survey and measure the stacked CMB lensing and ISW signals around them. To compare our measurements with ΛCDM model predictions, we construct a mock catalogue with matched galaxy number density and bias and apply the same photo-z uncertainty as the data. The consistency between the mock and the data is verified via the stacked galaxy density profiles around the superstructures and their quantity. The corresponding lensing convergence and ISW maps are then constructed and compared. The stacked lensing signal agrees with data well except at the highest redshift bin in density peaks, where the mock prediction is significantly higher, by approximately a factor of 1.3. The stacked ISW signal is generally consistent with the mock prediction. We do not obtain a significant signal from voids, A ISW = -0.10 ± 0.69, and the signal from clusters, A ISW = 1.52 ± 0.72, is at best weakly detected. However, these results are strongly inconsistent with previous claims of ISW signals at many times the level of the ΛCDM prediction. We discuss the comparison of our results with past work in this area and investigate possible explanations for this discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Velocity Dispersions of Massive Quiescent Galaxies from Weak Lensing and Spectroscopy

We use MMT spectroscopy and deep Subaru Hyper Suprime-Cam (HSC) imaging to compare the spectroscopic central stellar velocity dispersion of quiescent galaxies with the effective dispersion of the dark matter halo derived from the stacked lensing signal. The spectroscopic survey (the Smithsonian Hectospec Lensing Survey) provides a sample of 4585 quiescent galaxy lenses with measured line-of-sight central stellar velocity dispersion ($\sigma_{\rm SHELS}$) that is more than 85\% complete for $R < 20.6$, $D_{n}4000> 1.5$ and $M_{\star} > 10^{9.5}{\rm M}_{\odot}$.The median redshift of the sample of lenses is 0.32. We measure the stacked lensing signal from the HSC deep imaging. The central stellar velocity dispersion is directly proportional to the velocity dispersion derived from the lensing $\sigma_{\rm Lens}$, $\sigma_{\rm Lens} = (1.05\pm0.15)\sigma_{\rm SHELS}+(-21.17\pm35.19)$. The independent spectroscopic and weak lensing velocity dispersions probe different scales, $\sim3$kpc and $\gtrsim$ 100 kpc, respectively, and strongly indicate that the observable central stellar velocity dispersion for quiescent galaxies is a good proxy for the velocity dispersion of the dark matter halo. We thus demonstrate the power of combining high-quality imaging and spectroscopy to shed light on the connection between galaxies and their dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Addressing Experimental Challenges in Probing Dark Energy with Accuracy and Precision with the Rubin Observatory Legacy Survey of Space & Time (LSST)

With support from this award, the Stanford group under the leadership of PI Patricia Burchat focused on pixel-level algorithms that are key to the investigation of dark energy with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The work was carried out through strong engagement in the DOE-supported LSST Dark Energy Science Collaboration (DESC), for which the PI served as Deputy Spokesperson (07/2019 - 07/2021) and Technical Coordinator (07/2021 - 07/2023). Members of the Burchat group are Full DESC members and are among scientists who have been accepted to join the Rubin Observatory System Integration, Test, and Commissioning (SIT-Com) effort. The specific intellectual focus of the group was the optimization of gravitational lensing as a probe of the evolution of the distribution of dark matter in space and time to learn about the nature of dark energy. The Stanford group pursued an integrated research program to understand in detail, and then remediate or calibrate, systematic effects that could limit the sensitivity of the LSST to cosmic shear as a dark energy probe -- particularly those related to our understanding of the atmospheric point spread function (PSF), challenges associated with blended galaxy images, and impacts of residual camera artifacts. Key to this work was developing, enhancing, and validating the fidelity of simulations that play a key role in the development of our understanding of subtle systematic biases and then using these simulations to develop techniques to characterize and mitigate biases, and associated algorithms for implementation in the LSST or DESC science pipelines. The group’s work included contributions to the open-source GalSim simulation package, which is widely used in the community; simulations and studies of spatial correlations of atmospheric PSFs; studies of temporal correlations of the atmospheric PSF; impacts of chromatic PSFs on cosmic shear calibration; continued studies of anisotropic Gaussian process interpolation of PSF models across the focal plane; and measuring and mitigating impacts of residual camera artifacts on cosmic shear. A primary impact of this work is that it increases the scientific discovery potential of the DOE and NSF investments in the Rubin Observatory Legacy Survey of Space \& Time, by addressing challenges in interpreting astronomical images through more accurate descriptions of the blurring of images due to the atmosphere. Specifically, this work impacts our ability to understand the fundamental nature of dark energy.

Burchat, Patricia Rose↗
DOE OSTI2024

LoVoCCS. II. Weak Lensing Mass Distributions, Red-sequence Galaxy Distributions, and Their Alignment with the Brightest Cluster Galaxy in 58 Nearby X-Ray-luminous Galaxy Clusters

The Local Volume Complete Cluster Survey is an ongoing program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts 0.03 < z < 0.12 and X-ray luminosities in the 0.1–2.4 keV band L X500c > 10 44 erg s −1 ) with the Dark Energy Camera, capturing data in the u, g, r, i, z bands with a 5σ point source depth of approximately 25th–26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19° within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32° within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Bayesian approach to strong lens finding in the era of wide-area surveys

ABSTRACT The arrival of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), Euclid-Wide and Roman wide-area sensitive surveys will herald a new era in strong lens science in which the number of strong lenses known is expected to rise from $\mathcal {O}(10^3)$ to $\mathcal {O}(10^5)$. However, current lens-finding methods still require time-consuming follow-up visual inspection by strong lens experts to remove false positives which is only set to increase with these surveys. In this work, we demonstrate a range of methods to produce calibrated probabilities to help determine the veracity of any given lens candidate. To do this we use the classifications from citizen science and multiple neural networks for galaxies selected from the Hyper Suprime-Cam survey. Our methodology is not restricted to particular classifier types and could be applied to any strong lens classifier which produces quantitative scores. Using these calibrated probabilities, we generate an ensemble classifier, combining citizen science, and neural network lens finders. We find such an ensemble can provide improved classification over the individual classifiers. We find a false-positive rate of 10−3 can be achieved with a completeness of 46 per cent, compared to 34 per cent for the best individual classifier. Given the large number of galaxy–galaxy strong lenses anticipated in LSST, such improvement would still produce significant numbers of false positives, in which case using calibrated probabilities will be essential for population analysis of large populations of lenses and to help prioritize candidates for follow-up.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗
DOE OSTI2020

Signatures of self-interacting dark matter on cluster density profile and subhalo distributions

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Therefore, our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints (σ/m ≲ 2 cm 2 /g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

79 ASTRONOMY AND ASTROPHYSICS↗