Search NASA⌕ Search

SEARCH · Search NASA

Results for “quasars: general”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

How does the radio enhancement of broad absorption line quasars relate to colour and accretion rate?

ABSTRACT The origin of radio emission in different populations of radio-quiet quasars is relatively unknown, but recent work has uncovered various drivers of increased radio-detection fraction. In this work, we pull together three known factors: optical colour (g − i), C iv distance (a proxy for L/LEdd), and whether or not the quasar contains broad absorption lines (BALQSOs) which signify an outflow. We use SDSS (Sloan Digital Sky Survey) DR14 spectra along with the LOFAR Two Metre Sky Survey Data Release 2 and find that each of these properties have an independent effect. BALQSOs are marginally more likely to be radio-detected than non-BALQSOs at similar colours and L/LEdd, moderate reddening significantly increases the radio-detection fraction and the radio detection increases with L/LEdd above a threshold for all populations. We test a widely used simple model for radio wind shock emission and calculate energetic efficiencies that would be required to reproduce the observed radio properties. We discuss interpretations of these results concerning radio-quiet quasars more generally. We suggest that radio emission in BALQSOs is connected to a different physical origin than the general quasar population since they show different radio properties independent of colour and C iv distance.

Petley, James W. (ORCID:0000000244960754)↗
DOE OSTI2025

A multifrequency study of sub-parsec jets with the Event Horizon Telescope

Context. The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN.Aims. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion.Methods. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132–1654, J0006–0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2–86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole.Results. Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factor δ ∝ r≤0.5. This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radius B ∝ r−3, hinting at jet acceleration or efficient magnetic energy dissipation.Key words: techniques: interferometric / galaxies: active / galaxies: jets / galaxies: nuclei / quasars: general / quasars: supermassive black holes⋆⋆ Deceased.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1987

Is optical Fe II emission related to the soft X-ray properties of quasars?

Radio-quiet quasars generally show broad, blended multiplets of Fe II emission in their optical and UV spectra. Radio-loud quasars also show UV Fe II emission, but their optical Fe II emission is generally weaker. No satisfactory theory connecting the generation of Fe II and radio emission has been found to explain this effect. A second, well-established distinction between the two clases of quasar is in their X-ray properties: radio-loud quasars are more X-ray luminous, and recent results have shown that they also have systematically flatter soft X-ray slopes. Here it is proposed that the second effect causes the first; i.e., that the primary factor controlling the optical Fe II emission is the soft X-ray spectrum. This proposition is supported by X-ray and optical data for nine quasars, which shows a correlation between the soft X-ray slope and the strength of the optical Fe II emission. One of these quasars (1803+676) is radio-quiet, and yet its optical spectrum shows no evidence for Fe II emission. This quasar is also unusual in that it has a flat X-ray spectrum. This further supports the proposal that the X-ray spectrum is important in determining the relative strengths of UV and optical Fe II emission.

Wilkes, Belinda J.↗
NASA NTRS1986

IRAS observations of active galactic nuclei

IRAS data has been analyzed for 35 radio galaxies and 42 radio-loud quasars, using the coadded survey database. The detection rate is 50 percent with more than 40 percent detected in more than one band. Quasars generally have smooth energy spectra over four or five decades in frequency, consistent with incoherent synchrotron emission from 1 cm to 1 micron. Radio galaxies by contrast have discontinuous spectra, with a large range in the spectral index at IRAS wavelengths. The energy budget shows that the majority of the power per unit bandwidth emerges in the infrared (1-100 microns). The fraction of bolometric luminosity emitted in the infrared is 20-30 percent for the quasars and 40-60 percent for the radio galaxies.

Impey, Chris↗
DOE OSTI2025

GA-NIFS and EIGER: A merging quasar host at z = 7 with an overmassive black hole

The James Webb Space Telescope is revolutionising our ability to understand the host galaxies and local environments of high-z quasars. Here we obtain a comprehensive understanding of the host galaxy of the z = 7.08 quasar J1120+0641 by combining NIRSpec integral field spectroscopy with NIRCam photometry of the host continuum emission. Our emission-line maps reveal that this quasar host is undergoing a merger with a bright companion galaxy. The quasar host and the companion have similar dynamical masses of ∼10 10 M ⊙ , suggesting that this is a major galaxy interaction. Through detailed quasar subtraction and SED fitting using the NIRCam data, we obtained an estimate of the host stellar mass of M ∗ = (3.0 −1.4 +2.5 ) × 10 9 M ⊙ , with M ∗ = (2.7 −0.5 +0.5 ) × 10 9 M ⊙ for the companion galaxy. Using the Hβ Balmer line, we estimated a virial black hole mass of M BH = (1.9 −1.1 +2.9 ) × 109 M ⊙ . Thus, J1120+0641 has an extreme black hole–stellar mass ratio of M BH /M ∗ = 0.63 −0.31 +0.54 , which is ∼3 dex larger than expected by the local scaling relations between black hole and stellar mass. J1120+0641 is powered by an overmassive black hole with the highest reported black hole–stellar mass ratio in a quasar host that is currently undergoing a major merger. These new insights highlight the power of JWST for measuring and understanding these extreme first quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Uncertainty quantification of the virial black hole mass with conformal prediction

Precise measurements of the black hole mass are essential to gain insight on the black hole and host galaxy co-evolution. A direct measure of the black hole mass is often restricted to nearest galaxies and instead, an indirect method using the single-epoch virial black hole mass estimation is used for objects at high redshifts. However, this method is subjected to biases and uncertainties as it is reliant on the scaling relation from a small sample of local active galactic nuclei. In this study, we propose the application of conformalized quantile regression (CQR) to quantify the uncertainties of the black hole predictions in a machine learning setting. We compare CQR with various prediction interval techniques and demonstrated that CQR can provide a more useful prediction interval indicator. In contrast to baseline approaches for prediction interval estimation, we show that the CQR method provides prediction intervals that adjust to the black hole mass and its related properties. That is it yields a tighter constraint on the prediction interval (hence more certain) for a larger black hole mass, and accordingly, bright and broad spectral line width source. Using a combination of neural network model and CQR framework, the recovered virial black hole mass predictions and uncertainties are comparable to those measured from the Sloan Digital Sky Survey. The code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The mass distribution of quasars in optical time-domain surveys

The determination of supermassive black hole (SMBH) masses is the key to understanding the host galaxy build-up and the SMBH mass assembly histories. The SMBH masses of non-local quasars are frequently estimated via the single-epoch virial black-hole mass estimators, which may suffer from significant biases. Here we demonstrate a new approach to infer the mass distribution of SMBHs in quasars by modelling quasar ultraviolet /optical variability. Here, our inferred black hole masses are systematically smaller than the virial ones by 0.3–0.6 dex; the ∼0.3 dex offsets are roughly consistent with the expected biases of the virial black-hole mass estimators. In the upcoming time-domain astronomy era, our methodology can be used to constrain the cosmic evolution of quasar mass distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

OzDES Reverberation Mapping Program: Stacking analysis with Hβ, Mg ii , and C iv

ABSTRACT Reverberation mapping is the leading technique used to measure direct black hole masses outside of the local Universe. Additionally, reverberation measurements calibrate secondary mass-scaling relations used to estimate single-epoch virial black hole masses. The Australian Dark Energy Survey (OzDES) conducted one of the first multi-object reverberation mapping surveys, monitoring 735 AGN up to z ∼ 4, over 6 years. The limited temporal coverage of the OzDES data has hindered recovery of individual measurements for some classes of sources, particularly those with shorter reverberation lags or lags that fall within campaign season gaps. To alleviate this limitation, we perform a stacking analysis of the cross-correlation functions of sources with similar intrinsic properties to recover average composite reverberation lags. This analysis leads to the recovery of average lags in each redshift-luminosity bin across our sample. We present the average lags recovered for the Hβ, Mg ii, and C iv samples, as well as multiline measurements for redshift bins where two lines are accessible. The stacking analysis is consistent with the Radius–Luminosity relations for each line. Our results for the Hβ sample demonstrate that stacking has the potential to improve upon constraints on the R–L relation, which have been derived only from individual source measurements until now.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Stacked reverberation mapping of high-redshift quasars in DESI. I. Feasibility analysis

The broad-line region of quasars has long been probed by reverberation mapping techniques that measure time lags between continuum and broad emission-line variations. Stacked reverberation mapping has been proposed as a less observationally expensive alternative to traditional methods. This ensemble approach also reduces biases from small-number statistics. The Dark Energy Spectroscopic Instrument (DESI) is conducting the most extensive spectroscopic survey of quasars to date. We create mock light curves emulating expected DESI quasar observations at redshifts $1.48\lt z\lt 5.2$ and luminosities $44.68 \le \log \lambda L_{1350 \mathring{\rm A}{}} / \mathrm{erg\, s^{-1}} \le 45.99$ to test stacked reverberation mapping feasibility using sparse spectroscopic data paired with well-sampled photometric data. The pipeline, using the lag estimation code JAVELIN (Just Another Vehicle for Estimating Lags In Nuclei), successfully recovers the simulated C IV lags within 1σ of the true values using spectroscopic light curves composed of only a few spectral epochs (2–10) with irregular cadences. We investigate how observational factors, including C IV flux error magnitude, number of stacked quasars, and spectral epoch count, affect performance. This work motivates a pathway for future stacked reverberation mapping projects with large-scale spectroscopic surveys of quasars having $\ge 2$ spectroscopic observations. Our results suggest an economical alternative for constraining and extending the radius–luminosity relation to higher redshifts and luminosities. Subsequently, this relation can be employed more reliably in single-epoch black hole mass measurements and quasar cosmology in these distant regimes.

quasars: general, quasars: supermassive black hole↗
NASA NTRS2011

High Resolution Rapid Response Observations of Compact Radio Sources with the Ceduna Hobart Interferometer (CHI)

Context. Frequent, simultaneous observations across the electromagnetic spectrum are essential to the study of a range of astrophysical phenomena including Active Galactic Nuclei. A key tool of such studies is the ability to observe an object when it flares i.e. exhibits a rapid and significant increase in its flux density. Aims. We describe the specific observational procedures and the calibration techniques that have been developed and tested to create a single baseline radio interferometer. that can rapidly observe a flaring object. This is the only facility that is dedicated to rapid high resolution radio observations of an object south of -30 degrees declination. An immediate application is to provide rapid contemporaneous radio coverage of AGN flaring at y-ray frequencies detected by the Fermi Gamma-ray Space Telescope. Methods. A single baseline interferometer was formed with radio telescopes in Hobart, Tasmania and Ceduna, South Australia. A software correlator was set up at the University of Tasmania to correlate these data. Results. Measurements of the flux densities of flaring objects can be made using our observing strategy within half an hour of a triggering event. These observations can be calibrated with amplitude errors better than 20%. Lower limits to the brightness temperatures of the sources can also be calculated using CHI. Key words. instrumentation:interferometers - galaxies:active - galaxies:jets - galaxies:nuclei quasars:general gamma rays:galaxies- 1.

Blanchard, Jay M.↗
NASA NTRS2016

Extremes of the Jet-Accretion Power Relation of Blazars, as Explored by Nustar

Hard X-ray observations are crucial to study the non-thermal jet emission from high-redshift, powerful blazars. We observed two bright z > 2 flat-spectrum radio quasars (FSRQs) in hard X-rays to explore the details of their relativistic jets and their possible variability. S50014+81 (at z = 3.366) and B0222+185 (at z = 2.690) have been observed twice by the Nuclear Spectroscopic Telescope Array (NuSTAR) simultaneously with Swift X-ray Telescope, showing different variability behaviors. We found that NuSTAR is instrumental to explore the variability of powerful high-redshift blazars, even when no gamma-ray emission is detected. The two sources have proven to have respectively the most luminous accretion disc and the most powerful jet among known blazars. Thanks to these properties, they are located at the extreme end of the jet-accretion disc relation previously found for gamma-ray detected blazars, to which they are consistent.

galaxies: active – quasars: general – quasars:↗
DOE OSTI2023

The cosmic web of X-ray active galactic nuclei seen through the eROSITA Final Equatorial Depth Survey (eFEDS)

Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA’s contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited (FX > 6.5 × 10 -15 erg s -1 cm -2 ) AGN sample at low redshift, 0.05 < z < 0.55. Two summary statistics, the clustering using spectra from SDSS-V and galaxy-galaxy lensing with imaging from HSC, are measured and interpreted with halo occupation distribution and abundance matching models. Both models successfully account for the observations. We obtain an exceptionally complete view of the AGN halo occupation distribution. The population of AGNs is broadly distributed among halos with a mean mass of 3.9 -2.4 +2.0 × 10 12 M ⊙ . This corresponds to a large-scale halo bias of b(z = 0.34) = 0.99 -0.10 +0.08 . The central occupation has a large transition parameter, σ log 10 (M) = 1.28 ± 0.2. The satellite occupation distribution is characterized by a shallow slope, α sat = 0.73 ± 0.38. We find that AGNs in satellites are rare, with f sat < 20%. Most soft X-ray-selected AGNs are hosted by central galaxies in their dark matter halo. A weak correlation between soft X-ray luminosity and large-scale halo bias is confirmed (3.3σ). We discuss the implications of environmental-dependent AGN triggering. This study paves the way toward fully charting, in the coming decade, the coevolution of X-ray AGNs, their host galaxies, and dark matter halos by combining eROSITA with SDSS-V, 4MOST, DESI, LSST, and Euclid data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

VarIabiLity seLection of AstrophysIcal sources iN PTF (VILLAIN): I. Structure function fits to 71 million objects

Light-curve variability is well-suited to characterizing objects in surveys with high cadence and a long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences between quasars (QSOs), stars, and galaxies. We relate variability and colour information in preparation for future surveys. We fit joint likelihoods to structure functions (SFs) of 71 million PTF light curves with a Markov chain Monte Carlo method. For each object, we assume a power-law SF and extract two parameters: the amplitude on timescales of one year, A, and a power-law index, γ. With these parameters and colours in the optical (Pan-STARRS1) and mid-infrared (WISE), we identify regions of parameter space dominated by different types of spectroscopically confirmed objects from SDSS. Candidate QSOs, stars, and galaxies are selected to show their parameter distributions. QSOs show high-amplitude variations in the R band, and the highest γ values. Galaxies have a broader range of amplitudes and their variability shows relatively little dependency on timescale. With variability and colours, we achieve a photometric selection purity of 99.3% for QSOs. Even though hard cuts in monochromatic variability alone are not as effective as seven-band magnitude cuts, variability is useful in characterizing object subclasses. Through variability, we also find QSOs that were erroneously classified as stars in the SDSS. We discuss perspectives and computational solutions in view of the upcoming LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Discovering strongly lensed quasar candidates with catalogue-based methods from DESI Legacy Surveys

The Hubble tension, revealed by a ~5σ discrepancy between measurements of the Hubble-Lemaitre constant among observations of the early and local Universe, is one of the most significant problems in modern cosmology. In order to better understand the origin of this mismatch, independent techniques to measure H 0 , such as strong lensing time delays, are required. Notably, the sample size of such systems is key to minimising the statistical uncertainties and cosmic variance, which can be improved by exploring the datasets of large-scale sky surveys such as Dark Energy Spectroscopic Instrument (DESI). We identify possible strong lensing time-delay systems within DESI by selecting candidate multiply imaged lensed quasars from a catalogue of 24 440 816 candidate QSOs contained in the ninth data release of the DESI Legacy Imaging Surveys (DESI-LS). Using a friend-of-friends-like algorithm on spatial co-ordinates, our method generates an initial list of compact quasar groups. This list is subsequently filtered using a measure of the similarity of colours among a group’s members and the likelihood that they are quasars. A visual inspection finally selects candidate strong lensing systems based on the spatial configuration of the group members. We identified 620 new candidate multiply imaged lensed quasars (101 grade-A, 214 grade-B, 305 grade-C). This number excludes 53 known spectroscopically confirmed systems and existing candidate systems identified in other similar catalogues. When available, these new candidates will be further checked by combining the spectroscopic and photometric data from DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The miniJPAS survey quasar selection – II. Machine learning classification with photometric measurements and uncertainties

Astrophysical surveys rely heavily on the classification of sources as stars, galaxies, or quasars from multiband photometry. Surveys in narrow-band filters allow for greater discriminatory power, but the variety of different types and redshifts of the objects present a challenge to standard template-based methods. In this work, which is part of a larger effort that aims at building a catalogue of quasars from the miniJPAS survey, we present a machine learning-based method that employs convolutional neural networks (CNNs) to classify point-like sources including the information in the measurement errors. We validate our methods using data from the miniJPAS survey, a proof-of-concept project of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) collaboration covering ∼1 deg 2 of the northern sky using the 56 narrow-band filters of the J-PAS survey. Due to the scarcity of real data, we trained our algorithms using mocks that were purpose-built to reproduce the distributions of different types of objects that we expect to find in the miniJPAS survey, as well as the properties of the real observations in terms of signal and noise. We compare the performance of the CNNs with other well-established machine learning classification methods based on decision trees, finding that the CNNs improve the classification when the measurement errors are provided as inputs. The predicted distribution of objects in miniJPAS is consistent with the putative luminosity functions of stars, quasars, and unresolved galaxies. Our results are a proof of concept for the idea that the J-PAS survey will be able to detect unprecedented numbers of quasars with high confidence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Probing quasar viewing angle with the variability structure function

Given the anisotropic emission from quasar accretion discs, their viewing angle affects estimates of the quasar luminosity, black hole mass, and Eddington ratio. Discs appear overluminous when viewed pole-on and underluminous when viewed at high inclination. In radio-quiet quasars, the viewing angle is usually unknown, although spectroscopic indicators have been proposed. Here, we use a recently discovered universality in the variability structure function (SF) of quasar light curves (LCs), where all quasars show the same SF when clocks run in units of orbital time-scale. As an offset from the mean relation can be caused by incorrect orbital time-scales and thus incorrect luminosities, we correlate these offsets with suggested inclination indicators. We derive SFs from NASA/Asteroid Terrestrial-impact Last Alert System (ATLAS) LCs spanning ~6 yr of observation, using a sample of 183 luminous quasars with measured H β lines as well as 753 quasars with C $\tiny{IV}$ and Mg $\tiny{II}$ lines. Starting from the proposed orientation indicators, we expect quasars with narrower H β lines and with more blueshifted C $\tiny{IV}$ lines to be viewed more pole-on and thus appear overluminous. In contrast, our SF analysis finds that presumed pole-on discs appear underluminous, consistently for both line indicators. We discuss possible explanations for the behaviour of quasars with highly blueshifted C $\tiny{IV}$ lines irrespective of inclination angle, including dusty outflows that might render the accretion disc underluminous and flatter disc temperature profiles with longer orbital time-scales than in thin-disc models but reach no satisfying conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

No redshift evolution in the rest-frame ultraviolet emission line properties of quasars from z = 1.5 to z = 4.0

We analyse the rest-frame ultraviolet (UV) spectra of 2531 high-redshift (3.5 < z < 4.0) quasars from the Sloan Digital Sky Survey DR16Q catalogue. In combination with previous work, we study the redshift evolution of the rest-frame UV line properties across the entire redshift range, 1.5 < z < 4.0. We improve the systemic redshift estimates at z > 3.5 using a cross-correlation algorithm that employs high signal-to-noise template spectra spanning the full range in UV emission line properties. We then quantify the evolution of C iv and He ii emission line properties with redshift. The increase in C iv blueshifts with cosmological redshift can be fully explained by the higher luminosities of quasars observed at high redshifts. We recover broadly similar trends between the He ii equivalent width and C iv blueshift at both 1.5 < z < 2.65 and 3.5 < z < 4.0 suggesting that the blueshift depends systematically on the spectral energy density (SED) of the quasar and there is no evolution in the SED over the redshift range 1.5 < z < 4.0. C iv blueshifts are highest when L/L Edd ≥ 0.2 and M BH ≥ 10 9 M ⊙ for the entire 1.5 < z < 4.0 sample. We find that luminosity matching samples as a means to explore the evolution of their rest-frame UV emission line properties is only viable if the samples are also matched in the M BH –L/L Edd plane. Quasars at z ≥ 6 are on average less massive and have higher Eddington-scaled accretion rates than their luminosity-matched counterparts at 1.5 < z < 4.0, which could explain the observed evolution in their UV line properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

JWST lensed quasar dark matter survey – II. Strongest gravitational lensing limit on the dark matter free streaming length to date

This is the second in a series of papers in which we use JWST Mid Infrared Instrument multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of nine systems. The warm dust region is compact and sensitive to perturbations by populations of haloes down to masses $\sim 10^6$ M$_{\odot }$. Using these warm dust flux-ratio measurements in combination with five previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third- and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhaloes. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6}\, {\rm M}_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly $\alpha$ forest and Milky Way satellite galaxies.

dark matter↗