Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “astronomy data analysis”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Precise and accurate cosmology with CMB×LSS power spectra and bispectra

With the advent of a new generation of cosmological experiments that will provide high-precision measurements of the cosmic microwave background (CMB) and galaxies in the large-scale structure, it is pertinent to examine the potential of performing a joint analysis of multiple cosmological probes. In this paper, we study the cosmological information content contained in the one-loop power spectra and tree bispectra of galaxies cross-correlated with CMB lensing. We use the FFTLog method to compute angular correlations in spherical harmonic space, applicable for wide angles that can be accessed by forthcoming galaxy surveys. Here, we find that adding the bispectra and cross-correlations with CMB lensing offers a significant improvement in parameter constraints, including those on the total neutrino mass, M ν , and local non-Gaussianity amplitude, ƒ NL . In particular, our results suggest that the combination of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and CMB-S4 will be able to achieve σ(M ν )=42 meV from galaxy and CMB lensing correlations, and σ(M ν )=12 meV when further combined with the CMB temperature and polarization data, without any prior on the optical depth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Constraining galaxy–halo connection with high-order statistics

ABSTRACT We investigate using three-point statistics in constraining the galaxy–halo connection. We show that for some galaxy samples, the constraints on the halo occupation distribution parameters are dominated by the three-point function signal (over its two-point counterpart). We demonstrate this on mock catalogues corresponding to the Luminous red galaxies (LRGs), Emission-line galaxies (ELGs), and quasars (QSOs) targeted by the Dark Energy Spectroscopic Instrument (DESI) Survey. The projected three-point function for triangle sides less up to 20 h−1 Mpc measured from a cubic Gpc of data can constrain the characteristic minimum mass of the LRGs with a preci sion of 0.46 per cent. For comparison, similar constraints from the projected two-point function are 1.55 per cent. The improvements for the ELGs and QSOs targets are more modest. In the case of the QSOs, it is caused by the high shot-noise of the sample, and in the case of the ELGs, it is caused by the range of halo masses of the host haloes. The most time-consuming part of our pipeline is the measurement of the three-point functions. We adopt a tabulation method, proposed in earlier works for the two-point function, to significantly reduce the required compute time for the three-point analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The enigmatic Galactic Center excess: Spurious point sources and signal mismodeling

The Galactic Center GeV excess (GCE) has garnered great interest as a possible signal of either dark matter annihilation or some novel astrophysical phenomenon, such as a new population of gamma-ray emitting pulsars. In a companion paper, we showed that in a 10° radius region of interest (ROI) surrounding the Galactic Center, apparent evidence for GCE point sources (PSs) from non-Poissonian template fitting (NPTF) is actually an artifact of unmodeled north-south asymmetry of the GCE. In this work, we develop a simplified analytic description of how signal mismodeling can drive an apparent preference for a PS population, and demonstrate how the behavior pointed out in the companion paper also appears in simpler simulated datasets that contain no PS signals at all. We explore the generality of this behavior in the real gamma-ray data, and discuss the implications for past and future studies using NPTF techniques. While the drop in PS preference once north-south asymmetry is included is not ubiquitous in larger ROIs, we show that any overly-rigid signal model is expected to yield a spurious PS signal that can appear very convincing: as well as apparent significance comparable to what one would expect from a true PS population, the signal can exhibit stability against a range of variations in the analysis, and a source count function that is very consistent with previous apparent NPTF-based detections of a GCE PS population. This contrasts with previously studied forms of systematic mismodeling which are unlikely to mimic a PS population in the same way. Finally, in the light of this observation, and its explicit realization in the region where the GCE is brightest, we argue that a dominantly smooth origin for the GCE is not in tension with existing NPTF analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Misinterpreting modified gravity as dark energy: a quantitative study

Standard cosmological data analyses typically constrain simple phenomenological dark-energy parameters, for example the present-day value of the equation of state parameter, w 0 , and its variation with scale factor, w a . However, results from such an analysis cannot easily indicate the presence of modified gravity. Even if general relativity does not hold, experimental data could still be fit sufficiently well by a phenomenological w 0 w a CDM, unmodified-gravity model. Hence, it would be very useful to know if there are generic signatures of modified gravity in standard analyses. Here we present, for the first time to our knowledge, a quantitative mapping showing how modified gravity models look when (mis)interpreted within the standard unmodified-gravity analysis. Scanning through a broad space of modified-gravity (Horndeski) models, and assuming a near-future survey consisting of CMB, BAO, and SNIa observations, we report values of the best-fit set of cosmological parameters including (w 0 , w a ) that would be inferred if modified gravity were at work. We find that modified gravity models that can masquerade as standard gravity lead to very specific biases in standard-parameter spaces. Here, we also comment on implications for measurements of the amplitude of mass fluctuations described by the parameter S 8 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Pre-He White Dwarfs in Eclipsing Binaries. I. WASP 0131+28

We report the first BV light curves and high-resolution spectra of the post-mass transfer binary star WASP 0131+28 to study the absolute properties of extremely low-mass white dwarfs. From the observed spectra, the double-lined radial velocities were derived, and the effective temperature and rotational velocity of the brighter, more massive primary were found to be T{sub eff,1}=10,000±200 K and v{sub 1}sin i = 55 ± 10 km s{sup −1}, respectively. The combined analysis of the Transiting Exoplanet Survey Satellite (TESS) archive data and ours yielded the accurate fundamental parameters of the program target. The masses were derived to about 1.0% accuracy and the radii to 0.6%, or better. The secondary component's parameters of M{sub 2}=0.200±0.002 M{sub ⊙}, R{sub 2}=0.528±0.003 R{sub ⊙}, T{sub eff,2} = 11,186 ± 235 K, and L{sub 2}=3.9±0.3 L{sub ⊙} are in excellent agreement with the evolutionary sequence for a helium-core white dwarf of mass 0.203 M{sub ⊙}, and indicates that this star is halfway through the constant luminosity phase. The results presented in this article demonstrate that WASP 0131+28 is an EL CVn eclipsing binary in a thin disk, which is formed from the stable Roche-lobe overflow channel and composed of a main-sequence dwarf with a spectral type A0 and a pre-He white dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Atacama Cosmology Telescope: Summary of DR4 and DR5 Data Products and Data Access

Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg 2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Probing New physics with high-redshift quasars: axions and non-standard cosmology

The Hubble diagram of quasars, as candidates to “standardizable” candles, has been used to measure the expansion history of the Universe at late times, up to very high redshifts ( z ~ 7). It has been shown that this history, as inferred from the quasar dataset, deviates at ≳ 3 σ level from the concordance (ΛCDM) cosmology model preferred by the cosmic microwave background (CMB) and other datasets. In this article, we investigate whether new physics beyond ΛCDM (BΛCDM) or beyond the Standard Model (BSM) could make the quasar data consistent with the concordance model. We first show that an effective redshift-dependent relation between the quasar UV and X-ray luminosities, complementing previous phenomenological work in the literature, can potentially remedy the discrepancy. Such a redshift dependence can be realized in a BSM model with axion-photon conversion in the intergalactic medium (IGM), although the preferred parameter space is in tension with various other astrophysical constraints on axions, at a level depending on the specific assumptions made regarding the IGM magnetic field. We briefly discuss a variation of the axion model that could evade these astrophysical constraints. On the other hand, we show that models beyond ΛCDM such as one with a varying dark energy equation of state ( w CDM) or the phenomenological cosmographic model with a polynomial expansion of the luminosity distance, cannot alleviate the tension. The code for our analysis, based on emcee [1] and corner.py [2], is publicly available at github.com/ChenSun-Phys/high_z_candles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The VST ATLAS quasar survey I: Catalogue of photometrically selected quasar candidates

We present the VST ATLAS Quasar Survey, consisting of ~1229 000 quasar (QSO) candidates with 16 < g < 22.5 over ~4700 deg 2 . The catalogue is based on VST ATLAS+NEOWISE imaging surveys and aims to reach a QSO sky density of 130 deg –2 for $z$ < 2.2 and ~30 deg –2 for $z$ > 2.2. To guide our selection, we use X-ray/UV/optical/MIR data in the extended William Herschel Deep Field (WHDF) where we find a g < 22.5 broad-line QSO density of 269 ± 67 deg –2 , roughly consistent with the expected ~196 deg –2 . We find that ~25 per cent of our QSOs are morphologically classed as optically extended. Overall, we find that in these deep data, MIR, UV, and X-ray selections are ~70–90 per cent complete while X-ray suffers less contamination than MIR and UV. MIR is however more sensitive than X-ray or UV to $z$ > 2.2 QSOs at g < 22.5 and the S X (0.5 – 10 keV) > 1 X 10 -14 ergs cm –2 s –1 limit of eROSITA. We adjust the selection criteria from our previous 2QDES pilot survey and prioritize VST ATLAS candidates that show both UV and MIR excess, also selecting candidates initially classified as extended. We test our selections using data from DESI (which will be released in DR1) and 2dF to estimate the efficiency and completeness, and we use ANNz2 to determine photometric redshifts. Applying over the ~4700 deg 2 ATLAS area gives us ~ 917 000 $z$ < 2.2 QSO candidates of which 472 000 are likely to be $z$ < 2.2 QSOs, implying a sky density of ~100 deg –2 , which our WHDF analysis suggests will rise to at least 130 deg –2 when eROSITA X-ray candidates are included. At $z$ > 2.2, we find ~310() 000 candidates, of which 169 000 are likely to be QSOs for a sky density of ~36 deg –2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Statistics of Elsasser Increments in Solar Wind and Magnetohydrodynamic Turbulence

In this Letter we investigate the dependency with scale of the empirical probability distribution functions (PDF) of Elsasser increments using large sets of WIND data (collected between 1995 and 2017) near 1 au. The empirical PDF are compared to the ones obtained from high-resolution numerical simulations of steadily driven, homogeneous reduced MHD turbulence on a 2048 3 rectangular mesh. A large statistical sample of Alfvénic increments is obtained by using conditional analysis based on the solar wind average properties. The PDF tails obtained from observations and numerical simulations are found to have exponential behavior in the inertial range, with an exponential decrement that satisfies power laws of the form α l ∝ l –μ , where l is the scale size, with μ between 0.17 and 0.25 for observations and 0.43 for simulations. PDF tails were extrapolated assuming their exponential behavior extends to arbitrarily large increments in order to determine structure function scaling laws at very high orders. Our results point to potentially universal scaling laws governing the PDF of Elsasser increments and to an alternative approach to investigate high-order statistics in solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Theoretical Astrophysics at Fermilab

The Theoretical Astrophysics Group works on a broad range of topics ranging from string theory to data analysis in the Sloan Digital Sky Survey. The group is motivated by the belief that a deep understanding of fundamental physics is necessary to explain a wide variety of phenomena in the universe. During the three years 2001-2003 of our previous NASA grant, over 120 papers were written; ten of our postdocs went on to faculty positions; and we hosted or organized many workshops and conferences. Kolb and collaborators focused on the early universe, in particular and models and ramifications of the theory of inflation. They also studied models with extra dimensions, new types of dark matter, and the second order effects of super-horizon perturbations. S tebbins, Frieman, Hui, and Dodelson worked on phenomenological cosmology, extracting cosmological constraints from surveys such as the Sloan Digital Sky Survey. They also worked on theoretical topics such as weak lensing, reionization, and dark energy. This work has proved important to a number of experimental groups [including those at Fermilab] planning future observations. In general, the work of the Theoretical Astrophysics Group has served as a catalyst for experimental projects at Fennilab. An example of this is the Joint Dark Energy Mission. Fennilab is now a member of SNAP, and much of the work done here is by people formerly working on the accelerator. We have created an environment where many of these people made transition from physics to astronomy. We also worked on many other topics related to NASA s focus: cosmic rays, dark matter, the Sunyaev-Zel dovich effect, the galaxy distribution in the universe, and the Lyman alpha forest. The group organized and hosted a number of conferences and workshop over the years covered by the grant. Among them were:

Source record↗
DOE OSTI2023

Polarization and variability of compact sources measured in Planck time-ordered data

This paper introduces a new Planck Catalog of Polarized and Variable Compact Sources (PCCS-PV) comprising 153 sources, the majority of which are extragalactic. The data include both the total flux density and linear polarization measured by Planck with frequency coverage from 30 to 353 GHz, and temporal spacing ranging from days to years. We classify most sources as beamed, extragalactic radio sources; the catalog also includes several radio galaxies, Seyfert galaxies, and Galactic and Magellanic Cloud sources, including H II regions and planetary nebulae. An advanced extraction method applied directly to the multifrequency Planck time-ordered data, rather than the mission sky maps, was developed to allow an assessment of the variability of polarized sources. Our analysis of the time-ordered data from the Planck mission, tod2flux, allowed us to catalog the time-varying emission and polarization properties for these sources at the full range of polarized frequencies employed by Planck, 30-353 GHz. PCCS-PV provides the time-and frequency-dependent, polarized flux densities for all 153 sources. To illustrate some potential applications of the PCCS-PV, we conducted preliminary comparisons of our measurements of selected sources with published data from other astronomical instruments. In summary, we find general agreement between the Planck and the Institut de Radioastronomie Millimétrique (IRAM) polarization measurements as well as with the Metsähovi 37 GHz values at closely similar epochs. Here, these combined measurements also show the value of PCCS-PV results and the PCCS2 catalog for filling in missing spectral (or temporal) coverage and helping to define the spectral energy distributions of extragalactic sources. In turn, these results provide useful clues as to the physical properties of the sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DESI z ≳ 5 Quasar Survey. I. A First Sample of 400 New Quasars at z ~ 4.7–6.6

We report the first results of a high-redshift (z ≳ 5) quasar survey using the Dark Energy Spectroscopic Instrument (DESI). As a DESI secondary target program, this survey is designed to carry out a systematic search and investigation of quasars at 4.8 < z < 6.8. The target selection is based on the DESI Legacy Imaging Surveys (the Legacy Surveys) DR9 photometry, combined with the Pan-STARRS1 data and J-band photometry from public surveys. A first quasar sample has been constructed from the DESI Survey Validation 3 (SV3) and first-year observations until 2022 May. This sample includes more than 400 new quasars at redshift 4.7 ≤ z < 6.6, down to 21.5 magnitude (AB) in the z band, discovered from 35% of the entire target sample. Remarkably, there are 220 new quasars identified at z ≥ 5, more than one-third of existing quasars previously published at this redshift. The observations so far result in an average success rate of 23% at z > 4.7. The current spectral data set has already allowed analysis of interesting individual objects (e.g., quasars with damped Lyα absorbers and broad absorption line features), and statistical analysis will follow the survey's completion. A set of science projects will be carried out leveraging this program, including quasar luminosity function, quasar clustering, intergalactic medium, quasar spectral properties, intervening absorbers, and properties of early supermassive black holes. Additionally, a sample of 38 new quasars at z ~ 3.8–5.7 discovered from a pilot survey in the DESI SV1 is also published in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: cosmology from galaxy×CMB lensing and shear×CMB lensing

ABSTRACT We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence 〈δgκCMB〉, and source galaxy weak lensing shear with CMB lensing convergence 〈γκCMB〉. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of 〈δgκCMB〉 + 〈γκCMB〉 changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from 〈δgκCMB〉 and 〈γκCMB〉 are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in 〈δgκCMB〉 becomes subdominant to cosmic variance, preventing 〈δgκCMB〉 to further improve the constraints. This implies that to maximize the cosmological constraints from future 〈δgκCMB〉 + 〈γκCMB〉 analyses, we should focus more on the requirements on 〈γκCMB〉 instead of 〈δgκCMB〉. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Updated universal relations for tidal deformabilities of neutron stars from phenomenological equations of state

Equation of state (EOS) insensitive relations, so-called universal relations, between the neutron star (NS) compactness, its multipolar tidal deformability coefficients, and between the tidal parameters for binary systems are essential to break degeneracies in gravitational wave data analysis. Here, we validate and recalibrate these universal relations using a large set of almost 2 million phenomenological EOSs that are consistent with current observations. In doing so, we extend universal relations to a larger region of the EOS parameter space, most notably to softer EOSs and larger compactnesses. We show that waveform models that neglect higher-than-leading-order tidal deformations of the NSs accumulate as much as 3.5 radians of dephasing from 20Hz to merger. Further, we also perform a full Bayesian parameter estimation of the GW170817 data, and we compare the NS radius constraints produced using universal relations from the literature and the updated fits we propose here. We find that the new fits yield a NS radius that is smaller by about 500 meters. This difference is less than the statistical uncertainty on the radius at the signal-to-noise-ratio of GW170817, but it is significantly larger than the precision anticipated for next-generation detectors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Relation of Observable Stellar Parameters to Mass-loss Rate of AGB Stars in the LMC

Using the data set of Riebel et al. for 6889 pulsating AGB stars in the LMC, we have derived formulae for mass-loss rate as a function of luminosity and pulsation period, or luminosity and mass, in three ways, for each of five subsets of data: fundamental-mode oxygen-rich stars, first-overtone-mode oxygen-rich stars, fundamental-mode carbon stars, first-overtone-mode carbon stars, and extreme carbon stars. Using the distribution of the stars in period versus luminosity and mass versus luminosity, we are able to derive a power-law fit to the dependence of mass-loss rate on those quantities. This results in formulae that reproduce observed mass-loss rates and are in general agreement with the expectation from mass-loss models that the mass-loss rate is highly sensitive to luminosity, mass, and pulsation period. In the process of carrying out this analysis we have found radius–mass–luminosity relations and examined pulsation–mass–radius relations using published evolutionary and pulsation models. These allow us to derive mass and radius from the observed quantities luminosity and pulsation period. We also derived new mass-loss rate versus color relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Fermi-LAT Gamma-Ray Emission Discovered from the Composite Supernova Remnant B0453-685 in the Large Magellanic Cloud

Abstract We report the second extragalactic pulsar wind nebula (PWN) to be detected in the megaelectronvolt–gigaelectronvolt band by the Fermi-LAT, located within the Large Magellanic Cloud. The only other known PWN to emit in the Fermi band outside of the Milky Way is N157B, which lies to the west of the newly detected gamma-ray emission at an angular distance of 4°. Faint, pointlike gamma-ray emission is discovered at the location of the composite supernova remnant (SNR) B0453-685 with a ∼4 σ significance with energies ranging from 300 MeV–2 TeV. We present the Fermi-LAT data analysis of the new gamma-ray source, coupled with a detailed multiwavelength investigation to understand the nature of the observed emission. Combining the observed characteristics of the SNR and the physical implications from broadband modeling, we argue it is unlikely that the SNR is responsible for the gamma-ray emission. While the gamma-ray emission is too faint for a pulsation search, we try to distinguish between any pulsar and PWN component of SNR B0453-685 that could be responsible for the observed gamma-ray emission using semi-analytic models. We determine the most likely scenario is that the old PWN ( τ ∼ 14,000 yr) within B0453-685 has been impacted by the return of the SNR reverse shock with a possible substantial pulsar component below 5 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗