Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mock”

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 1,063 records · Page 59

Reconstruction of Galaxy Star Formation Histories through SED Fitting:The Dense Basis Approach

We introduce the dense basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M {sub *}, SFR, and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis star formation histories (SFHs) chosen by comparing the goodness-of-fit of mock galaxy SEDs to the goodness-of-reconstruction of their SFHs. We train and validate the method using a sample of realistic SFHs at z = 1 drawn from stochastic realizations, semi-analytic models, and a cosmological hydrodynamical galaxy formation simulation. The method is then applied to a sample of 1100 CANDELS GOODS-S galaxies at 1 10{sup 9.5} M{sub ⊙}. About 40% of the CANDELS galaxies have SFHs whose maximum occurs at or near the epoch of observation. The dense basis method is scalable and offers a general approach to a broad class of data-science problems.

79 ASTRONOMY AND ASTROPHYSICS↗

The Thermal Proximity Effect: A New Probe of the He ii Reionization History and Quasar Lifetime

Despite decades of effort, the timing and duration of He ii reionization and the properties of the quasars believed to drive it are still not well constrained. We present a new method to study both via the thermal proximity effect—the heating of the intergalactic medium (IGM) around quasars when their radiation doubly ionizes helium. We post-process hydrodynamical simulations with 1D radiative transfer and study how the thermal proximity effect depends on the He ii fraction, x{sub HeII,0}, which prevailed in the IGM before the quasar turned on, and the quasar lifetime t{sub Q}. We find that the amplitude of the temperature boost in the quasar environment depends on x{sub HeII,0}, with a characteristic value of ΔT≃10{sup 4} K for x{sub HeII,0}=1.0, whereas the size of the thermal proximity zone is sensitive to t{sub Q}, with typical sizes of ≃100 cMpc for t{sub Q}=10{sup 8} yr. This temperature boost increases the thermal broadening of H i absorption lines near the quasar. We introduce a new Bayesian statistical method based on measuring the Lyα forest power spectrum as a function of distance from the quasar, and demonstrate that the thermal proximity effect should be easily detectable. For a mock data set of 50 quasars at z≃4, we predict that one can measure x{sub HeII,0} to an (absolute) precision ≈0.04 and t{sub Q} to a precision of ≈0.1 dex. By applying our formalism to existing high-resolution Lyα forest spectra, one should be able to reconstruct the He ii reionization history, providing a global census of hard photons in the high-z universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting the Bulge–Halo Conspiracy. I. Dependence on Galaxy Properties and Halo Mass

We carry out a systematic investigation of the total mass density profile of massive (log M{sub star}/M{sub ⊙}≳11.3) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies with measured stellar masses, effective radii, and Sérsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical ΛCDM profile. Our main results are as follows. (1) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius γ′ is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (2) Models characterized by a Salpeter or variable initial mass function (IMF) and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier IMF and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (3) Currently available data on the mass density profiles of very massive galaxies (log M{sub star}/M{sub ⊙}≳12), with M{sub halo}≳3×10{sup 14} M{sub ⊙}, favor instead models with a stellar profile flatter than a Sérsic one in the very inner regions (r ≲ 3–5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ∼2 or ≲1.3, respectively, higher than those typically predicted by N-body numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Survival of Massive Star-forming Galaxies in Cluster Cores Drives Gas-phase Metallicity Gradients: The Effects of Ram Pressure Stripping

Recent observations of galaxies in a cluster at z = 0.35 show that their integrated gas-phase metallicities increase with decreasing cluster-centric distance. To test whether ram pressure stripping (RPS) is the underlying cause, we use a semianalytic model to quantify the “observational bias” that RPS introduces into the aperture-based metallicity measurements. We take integral field spectroscopy of local galaxies, remove gas from their outer galactic disks via RPS, and then conduct mock slit observations of cluster galaxies at z = 0.35. Our RPS model predicts a typical cluster-scale metallicity gradient of −0.03 dex/Mpc. By removing gas from the outer galactic disks, RPS introduces a mean metallicity enhancement of +0.02 dex at a fixed stellar mass. This gas removal and subsequent quenching of star formation preferentially removes low-mass cluster galaxies from the observed star-forming population. As only the more massive star-forming galaxies survive to reach the cluster core, RPS produces a cluster-scale stellar mass gradient of −0.05log(M{sub ∗}/M{sub ⊙})/Mpc. This mass segregation drives the predicted cluster-scale metallicity gradient of −0.03 dex/Mpc. However, the effects of RPS alone cannot explain the higher metallicities measured in cluster galaxies at z = 0.35. We hypothesize that additional mechanisms including steep internal metallicity gradients and self-enrichment due to gas strangulation are needed to reproduce our observations at z = 0.35.

79 ASTRONOMY AND ASTROPHYSICS↗

From Neutron Star Observables to the Equation of State. II. Bayesian Inference of Equation of State Pressures

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. Here, we present a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parameterized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, ρ{sub sat}, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ∼30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass–radius relation can lead to biases of nearly 1 km in radius, toward larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmological Constraints from the Redshift Dependence of the Volume Effect Using the Galaxy 2-point Correlation Function across the Line of Sight

We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line of sight, ξ(r{sub ⊥}), as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a redshift-dependent scaling in the galaxy distribution. This geometrical distortion can be observed as a redshift-dependent rescaling in the measured ξ(r{sub ⊥}). We test this methodology using a sample of 1.75 billion mock galaxies at redshifts 0, 0.5, 1, 1.5, and 2, drawn from the Horizon Run 4 N-body simulation. The shape of ξ(r{sub ⊥}) can exhibit a significant redshift evolution when the galaxy sample is analyzed under a cosmology differing from the true, simulated one. Other contributions, including the gravitational growth of structure, galaxy bias, and the redshift space distortions, do not produce large redshift evolution in the shape. We show that one can make use of this geometrical distortion to constrain the values of cosmological parameters governing the expansion history of the universe. This method could be applicable to future large-scale structure surveys, especially photometric surveys such as DES and LSST, to derive tight cosmological constraints. This work is a continuation of our previous works as a strategy to constrain cosmological parameters using redshift-invariant physical quantities.

79 ASTRONOMY AND ASTROPHYSICS↗

Robust Cross-correlation-based Measurement of Clump Sizes in Galaxies

Stars form in molecular complexes that are visible as giant clouds (∼10{sup 5−6} M{sub ⊙}) in nearby galaxies and as giant clumps (∼10{sup 8−9} M{sub ⊙}) in galaxies at redshifts z≈1–3. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g., Hα) can be recovered statistically using the two-point correlation function (2PCF) if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at ∼20% accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical measurements of Hα clump size agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

79 ASTRONOMY AND ASTROPHYSICS↗

X-Ray Eclipses of Active Galactic Nuclei

X-ray variation is a ubiquitous feature of active galactic nuclei (AGNs); however, its origin is not well understood. In this paper, we show that the X-ray flux variations in some AGNs, and correspondingly, the power spectral densities (PSDs) of the variations, may be interpreted as being caused by absorptions of eclipsing clouds or clumps in the broad-line region (BLR) and the dusty torus. By performing Monte Carlo simulations for a number of plausible cloud models, we systematically investigate the statistics of the X-ray variations resulting from the cloud eclipsing and the PSDs of the variations. For these models, we show that the number of eclipsing events can be significant and the absorption column densities due to those eclipsing clouds can be in the range from 10{sup 21} to 10{sup 24} cm{sup −2}, leading to significant X-ray variations. We find that the PSDs obtained from the mock observations for the X-ray flux and the absorption column density resulting from these models can be described by a broken double power law, similar to those directly measured from observations of some AGNs. The shape of the PSDs depend strongly on the kinematic structures and the intrinsic properties of the clouds in AGNs. We demonstrate that the X-ray eclipsing model can naturally lead to a strong correlation between the break frequencies (and correspondingly the break timescales) of the PSDs and the masses of the massive black holes (MBHs) in the model AGNs, which can be well consistent with the one obtained from observations. Future studies of the PSDs of the AGN X-ray (and possibly also the optical-UV) flux and column density variations may provide a powerful tool to constrain the structure of the BLR and the torus and to estimate the MBH masses in AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of Local Environment and Stellar Mass on Galaxy Quenching and Morphology at 0.5 < z < 2.0

We study galactic star formation activity as a function of environment and stellar mass over 0.5 < z < 2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, log(M/M{sub ⊙})>9 (9.5) at z = 1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision (σ{sub z}/(1+z)≲0.02) of ZFOURGE, accurately recovers galaxies in low- and high-density environments. We quantify the environmental quenching efficiency and show that at z>0.5, it depends on galaxy stellar mass, demonstrating that the effects of quenching related to (stellar) mass and environment are not separable. In high-density environments, the mass and environmental quenching efficiencies are comparable for massive galaxies (log(M/M{sub ⊙})≳10.5) at all redshifts. For lower-mass galaxies (log(M/M){sub ⊙})≲10), the environmental quenching efficiency is very low at z≳1.5, but increases rapidly with decreasing redshift. Environmental quenching can account for nearly all quiescent lower-mass galaxies (log(M/M{sub ⊙})∼9--10), which appear primarily at z≲1.0. The morphologies of lower-mass quiescent galaxies are inconsistent with those expected of recently quenched star-forming galaxies. Some environmental process must transform the morphologies on similar timescales as the environmental quenching itself. The evolution of the environmental quenching favors models that combine gas starvation (as galaxies become satellites) with gas exhaustion through star formation and outflows (“overconsumption”), and additional processes such as galaxy interactions, tidal stripping, and disk fading to account for the morphological differences between the quiescent and star-forming galaxy populations.

79 ASTRONOMY AND ASTROPHYSICS↗

Physical Origins of Gas Motions in Galaxy Cluster Cores: Interpreting Hitomi Observations of the Perseus Cluster

The Hitomi X-ray satellite has provided the first direct measurements of the plasma velocity dispersion in a galaxy cluster. It finds a relatively “quiescent” gas with a line-of-sight velocity dispersion σ{sub v,los}≃160 km s{sup −1}, at 30–60 kpc from the cluster center. This is surprising given the presence of jets and X-ray cavities that indicates on-going activity and feedback from the active galactic nucleus (AGN) at the cluster center. Using a set of mock Hitomi observations generated from a suite of state-of-the-art cosmological cluster simulations, and an isolated but higher resolution simulation of gas physics in the cluster core, including the effects of cooling and AGN feedback, we examine the likelihood of Hitomi detecting a cluster with the observed velocities. As long as the Perseus has not experienced a major merger in the last few gigayears, and AGN feedback is operating in a “‘gentle” mode, we reproduce the level of gas motions observed by Hitomi. The frequent mechanical AGN feedback generates net line-of-sight velocity dispersions ∼100--200 km s{sup −1}, bracketing the values measured in the Perseus core. The large-scale velocity shear observed across the core, on the other hand, is generated mainly by cosmic accretion such as mergers. We discuss the implications of these results for AGN feedback physics and cluster cosmology and progress that needs to be made in both simulations and observations, including a Hitomi re-flight and calorimeter-based instruments with higher spatial resolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-negative Matrix Factorization for Self-calibration of Photometric Redshift Scatter in Weak-lensing Surveys

Photo- z error is one of the major sources of systematics degrading the accuracy of weak-lensing cosmological inferences. Zhang et al. proposed a self-calibration method combining galaxy–galaxy correlations and galaxy–shear correlations between different photo- z bins. Fisher matrix analysis shows that it can determine the rate of photo- z outliers at a level of 0.01%–1% merely using photometric data and do not rely on any prior knowledge. In this paper, we develop a new algorithm to implement this method by solving a constrained nonlinear optimization problem arising in the self-calibration process. Based on the techniques of fixed-point iteration and non-negative matrix factorization, the proposed algorithm can efficiently and robustly reconstruct the scattering probabilities between the true- z and photo- z bins. The algorithm has been tested extensively by applying it to mock data from simulated stage IV weak-lensing projects. We find that the algorithm provides a successful recovery of the scatter rates at the level of 0.01%–1%, and the true mean redshifts of photo- z bins at the level of 0.001, which may satisfy the requirements in future lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

The Baryonic Collapse Efficiency of Galaxy Groups in the RESOLVE and ECO Surveys

We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at the group-integrated cold baryonic mass M{sub bary}{sup cold} ∼ 10{sup 11} M{sub ⊙}. The SAM, however, has significantly fewer groups at the transition mass ∼10{sup 11} M{sub ⊙} and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with a slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ∼2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M{sub halo}∼10{sup 11.4−12} M{sub ⊙}, which we label “nascent groups.” Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of Dark Matter Halo Mass from Dynamics of Satellite Galaxies

We propose that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using nine dark matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo, but has an intrinsic dependence on the halo formation history. We build this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of nine tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim 40 \% $. The accuracy can be improved to within ~25% if 30 tracers are used. Yet, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim 20 \% $ that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase-space information of other kinematic tracers is used.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Method to Measure the Post-reionization Ionizing Background from the Joint Distribution of Ly α and Ly β Forest Transmission

The amplitude of the ionizing background that pervades the intergalactic medium (IGM) at the end of the epoch of reionization provides a valuable constraint on the emissivity of the sources that reionized the universe. While measurements of the ionizing background at lower redshifts rely on a simulation-calibrated mapping between the photoionization rate and the mean transmission of the Lyα forest, at z ≳ 6 the IGM becomes increasingly opaque and transmission arises solely in narrow spikes separated by saturated Gunn–Peterson troughs. In this regime, the traditional approach of measuring the average transmission over large ~50 Mpc/h regions is less sensitive and suboptimal. In addition, the five times smaller oscillator strength of the Lyβ transition implies that the Lyβ forest is considerably more transparent at z ≳ 6, even in the presence of contamination by foreground z ~ 5 Lyα forest absorption. We present a novel statistical approach to analyze the joint distribution of transmission spikes in the cospatial z ~ 6 Lyα and Lyβ forests. Our method relies on approximate Bayesian computation (ABC), which circumvents the necessity of computing the intractable likelihood function describing the highly correlated Lyα and Lyβ transmission. We apply ABC to mock data generated from a large-volume hydrodynamical simulation combined with a state-of-the-art model of ionizing background fluctuations in the post-reionization IGM and show that it is sensitive to higher IGM neutral hydrogen fractions than previous techniques. As a proof of concept, we apply this methodology to a real spectrum of a z = 6.54 quasar and measure the ionizing background from 5.4 ≤ z ≤ 6.4 along this sightline with ~0.2 dex statistical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep Learning Identifies High-$z$ Galaxies in a Central Blue Nugget Phase in a Characteristic Mass Range

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a convolutional neural network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution— pre-BN, BN, and post-BN—and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of ~0.15 Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z$ = 1–3, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, 10 9.2–10.3 $M_⊙$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of z ~ 2.3 Cosmic Voids from 3D Ly α Forest Tomography in the COSMOS Field

We present the most distant detection of cosmic voids ( z ~ 2.3) and the first detection of three-dimensional voids in the Lyα forest. We used a 3D tomographic map of the absorption with an effective comoving spatial resolution of 2.5 h -1 Mpc and a volume of 3.15 &times; 10 5 h 3 Mpc 3 , which was reconstructed from moderate-resolution Keck I/LRIS spectra of 240 background Lyman-break galaxies and quasars in a 0.16 deg 2 footprint in the COSMOS field. Voids were detected using a spherical overdensity finder calibrated from hydrodynamical simulations of the intergalactic medium (IGM). This allows us to identify voids in the IGM corresponding to voids in the underlying matter density field, yielding a consistent volume fraction of voids in both data (19.5%) and simulations (18.2%). We fit excursion set models to the void radius function and compare the radially averaged stacked profiles of large voids ( r > 5 h -1 Mpc) to stacked voids in mock observations and the simulated density field. Comparing with 432 coeval galaxies with spectroscopic redshifts in the same volume as the tomographic map, we find that the tomography-identified voids are underdense in galaxies by 5.95σ compared to random cells.

79 ASTRONOMY AND ASTROPHYSICS↗

The Splashback Feature around DES Galaxy Clusters: Galaxy Density and Weak Lensing Profiles

Splashback refers to the process of matter that is accreting onto a dark matter halo reaching its first orbital apocenter and turning around in its orbit. The clustercentric radius at which this process occurs, r sp, defines a halo boundary that is connected to the dynamics of the cluster. A rapid decline in the halo profile is expected near r sp. We measure the galaxy number density and weak lensing mass profiles around redMaPPer galaxy clusters in the first-year Dark Energy Survey (DES) data. For a cluster sample with mean M 200m mass ≈2.5 × 10 14 M ⊙ , we find strong evidence of a splashback-like steepening of the galaxy density profile and measure r sp = 1.13 ± 0.07 h –1 Mpc, consistent with the earlier Sloan Digital Sky Survey measurements of More et al. and Baxter et al. Moreover, our weak lensing measurement demonstrates for the first time the existence of a splashback-like steepening of the matter profile of galaxy clusters. We measure r sp = 1.34 ± 0.21 h –1 Mpc from the weak lensing data, in good agreement with our galaxy density measurements. For different cluster and galaxy samples, we find that, consistent with ΛCDM simulations, r sp scales with R 200m and does not evolve with redshift over the redshift range of 0.3–0.6. We also find that potential systematic effects associated with the redMaPPer algorithm may impact the location of r sp . In conclusion, we discuss the progress needed to understand the systematic uncertainties and fully exploit forthcoming data from DES and future surveys, emphasizing the importance of more realistic mock catalogs and independent cluster samples.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Measurement of the Temperature–density Relation of the IGM from Voigt Profile Fitting

We decompose the Lyman-α (Lyα) forest of an extensive sample of 75 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-b cutoff in the b- distribution set primarily by the thermal state of the intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range 2.0 ≤ z ≤ 3.4, which allows us to measure the evolution of the IGM temperature-density (T = T 0 (ρ/ρ 0 ) γ-1 ) relation parameters T 0 and γ. We calibrate our measurements against mock Lyα forest data generated using 26 hydrodynamic simulations with different thermal histories from the THERMAL suite, also encompassing different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of T 0 (γ) shows a suggestive peak (dip) at z = 2.9 (z = 3). Our measured evolution of T 0 and γ is generally in good agreement with previous determinations in the literature. Both the peak in the evolution of T 0 at z = 2.8, as well as the high temperatures T 0 ≃ 15,000-20,000 K that we observe at 2.4 < z < 3.4, strongly suggest that a significant episode of heating occurred after the end of H I reionization, which was most likely the cosmic reionization of He II.

79 ASTRONOMY AND ASTROPHYSICS↗