Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feature selection”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Minimal cyclic behavior in sheared amorphous solids

Although jammed packings of soft spheres exist in potential energy landscapes with a vast number of minima, when subjected to cyclic shear they may revisit the same configurations repeatedly. Simple hysteretic spin models, in which particle rearrangements are represented by interacting spin flips called hysterons, capture many features of this periodic behavior. Yet it has been unclear to what extent individual rearrangements can be described by such binary objects and how such objects interact with one another. Using a particularly sensitive algorithm, we identify rearrangements in simulated jammed packings and select pairs of rearrangements that undo one another to create periodic cyclic behavior. We find that the rearrangement pairs surprisingly persist down to the smallest increments in strain, even in the smallest systems we can study. We explore the statistics of these rearrangement pairs and find that there is a relation between the amount of hysteresis and the energy drop and mean-square displacement of the particles; these results are inconsistent with the scaling found in models that treat rearrangements as localized buckling events. Finally, our analysis shows that there is no clean distinction between the particle motions that represent the identity of a single, individual rearrangement and the particle motions that lead to interactions between separated rearrangements or hysterons. These results offer insight into how complex systems such as amorphous solids can reach a limit cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Applicability of Unit Systems to High-Performance Computing Applications

Dimensional analysis is a key technique used to verify the soundness of scientific models. Most experts agree that engineering and scientific software would be made more reliable by integrating dimensional analysis in their type system. We explored how High Performance Computing (HPC) applications could integrate compile-time dimensional analysis. We started by investigating various implementation of unit systems for C++. Eventually, selecting the latest (and most advanced) one to apply to our test codes. We worked with code of increasing complexity, from a projectile trajectory calculation to the proxy-application Lulesh. This included our code, Springs-3D, which focuses on demonstrating language features while performing simple physic computations. Finally, our main contribution is a source-code analysis which extracts constraints on the dimension of all variables, functions, and constants in an application. This resulting system of equations is solved using the dimensions of a few of these objects. This analysis has the potential to greatly reduce the time spent performing dimensional analysis when refactoring application to use a representation of units.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗
DOE OSTI2022

Unitary Selective Coupled-Cluster Method

Simulating molecules using the Variational Quantum Eigensolver method is one of the promising applications for NISQ-era quantum computers. Designing an efficient ansatz to represent the electronic wave function is crucial in such simulations. Standard unitary coupled-cluster with singles and doubles (UCCSD) ansatz tends to have a large number of insignificant terms that do not lower the energy of the system. In this work, we present a unitary selective coupled-cluster method, a way to construct a unitary coupled-cluster ansatz iteratively using a selection procedure with excitations up to fourth order. This approach uses the electronic Hamiltonian matrix elements and the amplitudes for excitations already present in the ansatz to find the important excitations of higher order and to add them to the ansatz. The important feature of the method is that it systematically reduces the energy error with increasing ansatz size for a set of test molecules. {The main advantage of the proposed method is that the effort to increase the ansatz does not require any additional measurements on a quantum computer.}

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Selective Methane Oxidation by Heterogenized Iridium Catalysts

Oxidative methane (CH 4 ) carbonylation promises a direct route to the synthesis of value-added oxygenates such as acetic acid (CH 3 COOH). Here, in this paper, we report a strategy to realize oxidative CH 4 carbonylation through immobilized Ir complexes on an oxide support. Our immobilization approach not only enables direct CH 4 activation but also allows for easy separation and reutilization of the catalyst. Furthermore, we show that a key step, methyl migration, that forms a C–C bond, is sensitive to the electrophilicity of carbonyl, which can be tuned by a gentle reduction to the Ir centers. While the as-prepared catalyst that mainly featured Ir(IV) preferred CH 3 COOH production, a reduced catalyst featuring predominantly Ir(III) led to a significant increase of CH 3 OH production at the expense of the reduced yield of CH 3 COOH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Exploitation of Pore Structure for Increased CO 2 Selectivity in Type 3 Porous Liquids

CO 2 capture requires materials with high adsorption selectivity and an industrial ease of implementation. To address these needs, a new class of porous materials was recently developed that combines the fluidity of solvents with the porosity of solids. Type 3 porous liquids (PLs) composed of solvents and metal–organic frameworks (MOFs) offer a promising alternative to current liquid carbon capture methods due to the inherent tunability of the nanoporous MOFs. However, the effects of MOF structural features and solvent properties on CO 2 –MOF interactions within PLs are not well understood. Herein experimental and computational data of CO 2 gas adsorption isotherms were used to elucidate both solvent and pore structure influences on ZIF-based PLs. The roles of the pore structure including solvent size exclusion, structural environment, and MOF porosity on PL CO 2 uptake were examined. A comparison of the pore structure and pore aperture was performed using ZIF-8, ZIF-L, and amorphous-ZIF-8. Adsorption experiments here have verified our previously proposed solvent size design principle for ZIF-based PLs (1.8× ZIF pore aperture). Furthermore, the CO 2 adsorption isotherms of the ZIF-based PLs indicated that judicious selection of the pore environment allows for an increase in CO 2 selectivity greater than expected from the individual PL components or their combination. This nonlinear increase in the CO 2 selectivity is an emergent behavior resulting from the complex mixture of components specific to the ZIF-L + 2'-hydroxyacetophenone-based PL.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Consequences of Intrapore Liquids on Reactivity, Selectivity, and Stability for Aldol Condensation Reactions on Anatase TiO 2 Catalysts

Abstract This study provides evidence and mechanistic interpretations for the significant consequences of intrapore non‐polar liquids on acetone aldol condensation turnover rates, selectivity to primary dimer products, and catalyst stability for reactions at Lewis acid‐base site pairs on TiO 2 surfaces. These non‐polar liquids confer such benefits through the preferential stabilization of transition states (TS) for adsorption (“entry”) and desorption (“exit”) steps, which place their respective reactants or products within a solvating outer sphere environment. The extent to which non‐polar fluids ( n ‐heptane) form an intrapore liquid phase within TiO 2 voids was obtained from N 2 uptakes using established formalisms that consider the different molal volume, surface tension, and volatility between N 2 and n ‐heptane. Acetone condensation rates are limited by C−H activation, an “entry” step that forms bound prop‐1‐en‐2‐olates via a TS stabilized by intrapore liquids, leading to higher aldol condensation turnover rates as n ‐heptane pressure increases and active TiO 2 surfaces become increasingly immersed within a non‐polar liquid phase. These liquids solvate the late TS structures that mediate the desorption of primary C 6 condensation products even more effectively than those involved in prop‐1‐en‐2‐olate formation or in nucleophilic attack events that later form C−C bonds. Such preferential solvation favors desorption over C−C coupling events, thus disfavoring the formation of larger oligomers that become stranded at active sites, thus leading to much slower deactivation. Moreover, solvation by non‐polar liquids also leads to C 6 alkanones as the sole products formed in a single surface sojourn. These effects of a non‐polar dense phase circumvent the inherent stability, reactivity, and selectivity hurdles that have precluded practical aldol condensation catalysis on Lewis acid‐base pairs at oxide surfaces; these consequences are demonstrated here for TiO 2 catalysts, acetone aldol condensation reactions, and n ‐heptane as the non‐polar liquid but through strategies, concepts, and mechanistic features that extend to other systems. More generally, these observations and their mechanistic origins demonstrate how a contacting liquid preferentially solvates TS structures for elementary steps that involve either reactants arriving from or products entering into an outer sphere environment that contains a dense non‐polar phase.

Kadam, Shashikant A.↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The role of chalcogen vacancies for atomic defect emission in MoS 2

For two-dimensional (2D) layered semiconductors, control over atomic defects and understanding of their electronic and optical functionality represent major challenges towards developing a mature semiconductor technology using such materials. Here, we correlate generation, optical spectroscopy, atomic resolution imaging, and ab initio theory of chalcogen vacancies in monolayer MoS 2 . Chalcogen vacancies are selectively generated by in-vacuo annealing, but also focused ion beam exposure. The defect generation rate, atomic imaging and the optical signatures support this claim. We discriminate the narrow linewidth photoluminescence signatures of vacancies, resulting predominantly from localized defect orbitals, from broad luminescence features in the same spectral range, resulting from adsorbates. Vacancies can be patterned with a precision below 10 nm by ion beams, show single photon emission, and open the possibility for advanced defect engineering of 2D semiconductors at the ultimate scale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Robust unfolding of MeV x-ray spectra from filter stack spectrometer data

Here, we present an inversion method capable of robustly unfolding MeV x-ray spectra from filter stack spectrometer (FSS) data without requiring an a priori specification of a spectral shape or arbitrary termination of the algorithm. Our inversion method is based upon the perturbative minimization (PM) algorithm, which has previously been shown to be capable of unfolding x-ray transmission data, albeit for a limited regime in which the x-ray mass attenuation coefficient of the filter material increases monotonically with x-ray energy. Our inversion method improves upon the PM algorithm through regular smoothing of the candidate spectrum and by adding stochasticity to the search. With these additions, the inversion method does not require a physics model for an initial guess, fitting, or user-selected termination of the search. Instead, the only assumption made by the inversion method is that the x-ray spectrum should be near a smooth curve. Testing with synthetic data shows that the inversion method can successfully recover the primary large-scale features of MeV x-ray spectra, including the number of x-rays in energy bins of several-MeV widths to within 10%. Fine-scale features, however, are more difficult to recover accurately. Examples of unfolding experimental FSS data obtained at the Texas Petawatt Laser Facility and the OMEGA EP laser facility are also presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Evolution of nontrivial Fermi surface features in the band structures of the homologous members Pb 5 Bi 6 Se 14 and Pb 5 Bi 12 Se 23

High-quality single crystals of (PbSe) 5 (Bi 2 Se 3 ) 3m were grown and analyzed by nuclear magnetic resonance (NMR) spectroscopy. We report on 77 Se and 207 Pb NMR shifts and nuclear spin-lattice relaxation measurements in the naturally formed heterostructure homology (PbSe) 5 (Bi 2 Se 3 ) 3m with m = 1 (Pb 5 Bi 6 Se 14 ) and m = 2 (Pb 5 Bi 12 Se 23 ). A distinct site-specific contribution has been detected for both nuclei as a function of temperature, which reveals an electronic changeover from a semiconducting (Pb 5 Bi 6 Se 14 ) to a semimetalliclike (Pb 5 Bi 12 Se 23 ) system with nontrivial band structure features near the Fermi level. The temperature dependences of the relaxation rates are dominated by significant changes in the topology of energy dispersions accompanied with band edges and crossings in the region of the Fermi surface. These results, which interrogate nuclear spin interactions from selected atomic sites, clearly expose the effects of the added Bi 2 Se 3 layer on the crystal and electronic structure of (Pb 5 Bi 12 Se 23 ). These findings provide direct microscopic insight into the unconventional and dual nature of the electronic structure of these homologous thermoelectric and topologically nontrivial compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Muon Neutrino Reconstruction at ICARUS with Machine Learning

The ICARUS T600 LArTPC detector successfully ran for three years at the underground LNGS laboratories, providing a first sensitive search for LSND-like anomalous electron neutrino appearance in the CNGS beam. After a significant overhauling at CERN, the T600 detector has been placed in its experimental hall at Fermilab, fully commissioned, and the first events observed with full detector readout. Regular data-taking began in May 2021 with neutrinos from the Booster Neutrino Beam (BNB) and neutrinos six degrees off-axis from the Neutrinos at the Main Injector (NuMI). Modern developments in machine learning have allowed for the development of an end-to-end machine learning-based event reconstruction for ICARUS data. This reconstruction folds in 3D voxel-level feature extraction using sparse convolutional neural networks and particle clustering using graph neural networks to produce outputs suitable for physics analyses. This poster will summarize the performance of a high-purity and high-efficiency end-to-end machine learning-based selection of muon neutrinos from the BNB and highlight studies of electromagnetic shower reconstruction from a neutral pion selection.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Feature extraction for subtle anomaly detection using semi-supervised learning

The demand for automated and effective monitoring techniques has soared with the increased digitization of industrial monitoring systems. State-of-the-art machine learning methods are effectively detecting abrupt changes in system states. However, these methods lack comparable maturity in detecting subtle changes that may be signs of incipient faults. This manuscript argues that the current anomaly detection methods can be enhanced by exploring weak patterns to enable subtle variation detection. Specifically, the concept of semi-supervised learning is employed, with labels representing knowledge about some anomalous conditions of a system. The basic idea is to extract a candidate set of weak patterns discarded by state-of-the-art baselining algorithms. With few labeled anomalous data, the algorithm selects the weak patterns and allows for their possible fusion using the highest sensitivity to the labeled anomalies. Here, the method’s applicability is demonstrated using a representative pressurized water reactor (PWR) model simulated by Dymola.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Solvent selection criteria for temperature-resilient lithium–sulfur batteries

All-climate temperature operation capability and increased energy density have been recognized as two crucial targets, but they are rarely achieved together in rechargeable lithium (Li) batteries. Herein, we demonstrate an electrolyte system by using monodentate dibutyl ether with both low melting and high boiling points as the sole solvent. Its weak solvation endows an aggregate solvation structure and low solubility toward polysulfide species in a relatively low electrolyte concentration (2 mol L −1 ). These features were found to be vital in avoiding dendrite growth and enabling Li metal Coulombic efficiencies of 99.0%, 98.2%, and 98.7% at 23 °C, −40 °C, and 50 °C, respectively. Pouch cells employing thin Li metal (50 μm) and high-loading sulfurized polyacrylonitrile (3.3 mAh cm −2 ) cathodes (negative-to-positive capacity ratio = 2) output 87.5% and 115.9% of their room temperature capacity at −40 °C and 50 °C, respectively. This work provides solvent-based design criteria for a wide temperature range Li-sulfur pouch cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Review of the Effects of Polymer Binder Properties on Microstructure and Irreversible Volume Growth of Plastic Bonded Explosives Formulations

The rational design of effective polymeric binders for the formulation of plastic-bonded explosives (PBX) is challenging due to their inherent compositional complexity. The composites comprise irregularly shaped energetic material (EM) powders coated with low weight fractions of polymer via non-equilibrium processes such as slurry coating. Defects can deleteriously affect PBX stability and performance: nano- to micrometer-scale voids can act as loci for hot spots, lowering deflagration and detonation temperatures in unpredictable ways. Furthermore, some nominally desirable polymer properties are at odds with each other: e. g. good flow characteristics are desirable for coating and adhesion, but mechanical stiffness is needed to prevent deformation and cracking of PBX under mechanical stress. Good binder adhesion is critical, but the best means to predict and measure adhesion in PBX is not obvious. Experimental methods of determining binder adhesion on model surfaces may not capture polymer structural configurations relevant to deposition during coating. Molecular dynamics-based computational models have predicted key observables in PBX formulation, suggesting that they may be powerful tools for binder selection. In this review, primarily recent (~2006 and later) literature on polymeric binders for insensitive HE (IHE) is surveyed. We focus on how binder properties influence observable PBX properties as resistance to irreversible volume growth and void formation in PBX formulations mainly (but not exclusively) featuring 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the EM (e. g. PBX 9502, LX-17 and others). Conclusions from these studies yield useful guidelines for choosing HE binder candidates, as well as for the general design of highly-filled polymer composite materials. Finally, studies describing challenges in PBX formulation with the newer and more energetically dense high explosive, LLM-105, will be discussed.

adhesion↗
DOE OSTI2020

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗
DOE OSTI2023

Correlations between flow and transverse momentum in Xe + Xe and Pb + Pb collisions at the LHC with the ATLAS detector: A probe of the heavy-ion initial state and nuclear deformation

The correlations between flow harmonics v n for n = 2, 3, and 4 and mean transverse momentum [p T ] in 129 Xe + 129 Xe and 208 Pb + 208 Pb collisions at $\sqrt{s}$ = 5.44 and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are potentially sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from nonflow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and an event-activity selection based on particle production at very forward rapidity. The v n -[p T ] correlations show strong dependencies on centrality, harmonic number n, p T , and pseudorapidity range. Current models qualitatively describe the overall centralityand system-dependent trends but fail to quantitatively reproduce all features of the data. In central collisions, where models generally show good agreement, the v 2 -[p T ] correlations are sensitive to the triaxiality of the quadruple deformation. Comparison of the model with the Pb + Pb and Xe + Xe data confirms that the 129 Xe nucleus is a highly deformed triaxial ellipsoid that has neither a prolate nor oblate shape. This provides strong evidence for a triaxial deformation of the 129 Xe nucleus from high-energy heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 in the presence of a generalized incident local radiation field computed by Cloudy. 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-16Ry) 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 ≳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↗