Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERSTELLAR RADIATION”

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

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Analysis of Polarized Dust Emission Using Data from the First Flight of SPIDER

Using data from the first flight of Spider and from the Planck High Frequency Instrument, we probe the properties of polarized emission from interstellar dust in the Spider observing region. Component-separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses of diffuse Galactic dust emission spanning the full Spider region demonstrate (i) a spectral energy distribution that is broadly consistent with a modified-blackbody (MBB) model with a spectral index of β d = 1.45 ± 0.05 (1.47 ± 0.06) for E (B)-mode polarization, slightly lower than that reported by Planck for the full sky; (ii) an angular power spectrum broadly consistent with a power law; and (iii) no significant detection of line-of-sight polarization decorrelation. Tests of several modeling uncertainties find only a modest impact (~10% in σ r ) on Spider's sensitivity to the cosmological tensor-to-scalar ratio. The size of the Spider region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature T d = 19.6 K, an analysis of two independent subregions of that field results in inferred values of β d = 1.52 ± 0.06 and β d = 1.09 ± 0.09, which are inconsistent at the 3.9σ level. Furthermore, a joint analysis of Spider and Planck 217 and 353 GHz data within one subregion is inconsistent with a simple MBB at more than 3σ, assuming a common morphology of polarized dust emission over the full range of frequencies. This evidence of variation may inform the component-separation approaches of future cosmic microwave background polarization experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Anatomy of an Unusual Edge-on Protoplanetary Disk. II. Gas Temperature and a Warm Outer Region

We present high-resolution {sup 12}CO and {sup 13}CO 2–1 ALMA observations, as well as optical and near-infrared spectroscopy, of the highly inclined protoplanetary disk around SSTC2D J163131.2–242627. The spectral type we derive for the source is consistent with a 1.2 M {sub ⊙} star inferred from the ALMA observations. Despite its massive circumstellar disk, we find little to no evidence for ongoing accretion on the star. The CO maps reveal a disk that is unusually compact along the vertical direction, consistent with its appearance in scattered light images. The gas disk extends about twice as far away as both the submillimeter continuum and the optical scattered light. CO is detected from two surface layers separated by a midplane region in which CO emission is suppressed, as expected from freeze-out in the cold midplane. We apply a modified version of the tomographically reconstructed distribution method presented by Dutrey et al. to derive the temperature structure of the disk. We find a temperature in the CO-emitting layers and the midplane of ∼33 K and ∼20 K at R < 200 au, respectively. Outside of R > 200 au, the disk’s midplane temperature increases to ∼30 K, with a nearly vertically isothermal profile. The transition in CO temperature coincides with a dramatic reduction in the submicron and submillimeter emission from the disk. We interpret this as interstellar UV radiation providing an additional source of heating to the outer part of the disk.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

The Effect of Bars on the Ionized ISM: Optical Emission Lines from Milky Way Analogs

Gas interior to the bar of the Milky Way has recently been shown to be the closest example of a low-ionization (nuclear) emission region (LI(N)ER) in the universe. To better understand the nature of this gas, a sample of face-on galaxies with integral field spectroscopy is used to study the ionized gas conditions of 240 barred and 250 nonbarred galaxies, focusing on those that are most similar to the Milky Way. Strong optical line emission of [N ii] λ6584, Hα, [O iii] λ5007, and Hβ are used to diagnose the dominant ionization mechanisms of gas across galaxies and the Galaxy via Baldwin-Phillips-Terlevich diagrams. Barred galaxies show a strong suppression of star formation and an increase in composite and LI(N)ER-like spectra in their inner regions when compared with similar nonbarred counterparts. This effect is lessened in galaxies of very low (log 10 (M * /M ⊙ ) ≲ 10.4) or very high (log 10 (M * /M ⊙ ) ≳ 11.1) total stellar mass. Bar masks from Galaxy Zoo:3D show the bar's nonaxisymmetric effect on the ionized gas and help predict the face-on distribution of ionized gas conditions near the bar of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The star formation, dust, and abundance of galaxies with unWISE -CIB cross-correlations

Here, the cosmic infrared background (CIB) is the accumulated infrared (IR) radiation mainly from interstellar dust heated up by early stars. In this work, we measure the cross-correlation between galaxies from the unWISE catalog and the CIB maps from the Planck satellite to simultaneously constrain the cosmic star formation rate (SFR), dust spectral energy distribution (SED), and the halo occupation distribution (HOD). The unWISE galaxy catalog is divided into three tomographic bins centered at z ~ 0.6, 1.1, 1.5, and the CIB maps are at 353, 545, and 857 GHz. We measure the cross-correlations between these galaxy samples and CIB maps and get a 194σ signal within an angular scale 100 < ℓ < 2000, from which we constrain two CIB halo models from previous literature and one new model. The SFR, SED, and HOD model parameters are constrained consistently among the three models. Specifically, the dust temperature at z = 0 is constrained T 0 = $21.14^{+1.02}_{-1.34}$ K, which is slightly lower than T 0 = 24.4±1.9 K measured by the Planck collaboration. The halo mass that gives the most efficient star formation is around 10 11.79 $^{+0.73}_{-0.86}$ M ⊙ . From the model parameters, combined with the SFR density at z = 0 synthesized from multi-wavelength observations, we break the degeneracy between SED and SFR and recover the cosmic star formation history that is consistent with multi-wavelength surveys. We also constrain the graybody SED model in agreement with previous measurements from infrared flux stacking. From the HOD constraints, we derive an increasing trend of galaxy linear bias along redshifts that agrees with the results from cross- and auto-correlation with unWISE galaxies. This study indicates the power of using CIB-galaxy cross-correlation to study star formation, dust, and abundance of galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Diagnosis of 3D magnetic field and mode composition in MHD turbulence with Y -parameter

ABSTRACT Magnetic fields are crucial in numerous astrophysical processes within the interstellar medium (ISM). However, the detailed determination of magnetic field geometry is notoriously challenging. Based on the modern magnetohydrodynamic (MHD) turbulence theory, we introduce a novel statistical technique, the ‘Y-parameter’, to decipher the magnetic field inclination in the ISM and identify dominant turbulence modes. The Y-parameter, calculated as the ratio of anisotropies of different Stokes parameter combinations, displays contrasting trends with the mean-field inclination angle in Alfvénic and compressible turbulence modes. A Y-parameter value around 1.5 ± 0.5 provide a statistical boundary to determine the dominant MHD turbulence modes. We have discovered specific correlations between the Y-parameter value and the inclination angle that unveil the dominant turbulence mode. This methodology, when applied to future radio polarization surveys such as LOFAR and SKA, promises to significantly enhance our knowledge of 3D magnetic field in the ISM and improve our understanding of interstellar turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol , aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from Bmagpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Search for a diffuse flux of photons with energies above tens of PeV at the Pierre Auger Observatory

Diffuse photons of energy above 0.1 PeV, produced through the interactions between cosmic rays and either interstellar matter or background radiation fields, are powerful tracers of the distribution of cosmic rays in the Galaxy. Furthermore, the measurement of a diffuse photon flux would be an important probe to test models of super-heavy dark matter decaying into gamma-rays. In this work, we search for a diffuse photon flux in the energy range between 50 PeV and 200 PeV using data from the Pierre Auger Observatory. For the first time, we combine the air-shower measurements from a 2 km 2 surface array consisting of 19 water-Cherenkov surface detectors, spaced at 433 m, with the muon measurements from an array of buried scintillators placed in the same area. Using 15 months of data, collected while the array was still under construction, we derive upper limits to the integral photon flux ranging from 13.3 to 13.8 km -2 sr -1 yr -1 above tens of PeV. We extend the Pierre Auger Observatory photon search program towards lower energies, covering more than three decades of cosmic-ray energy. This work lays the foundation for future diffuse photon searches: with the data from the next 10 years of operation of the Observatory, this limit is expected to improve by a factor of ∼20.

cosmic rays detectors↗
DOE OSTI2024

Synchrotron-driven Instabilities in Relativistic Plasmas of Arbitrary Opacity

Recent work has shown that synchrotron emission from relativistic plasmas leads the electron distribution to form an anisotropic ring in momentum space, which can be unstable to both kinetic and hydrodynamic instabilities. Fundamental to these works was the assumption that the plasma was optically thin, allowing all emitted radiation to escape. Here, we examine the behavior of these instabilities as the plasma becomes more optically thick. To do this, we extend a recently developed Fokker–Planck operator for synchrotron emission and absorption in mildly relativistic plasmas to ultrarelativistic plasmas. For a given set of plasma parameters, photons emitted by higher-energy electrons tend to be higher frequency, and thus more easily escape the plasma. As a result, the ratio of the photon emission rate (radiative drag) to absorption rate (radiative diffusion) for a given electron is extremely energy dependent. Given this behavior, we determine the critical parameters that control the opacity, and show how the plasma gradually transitions to become more isotropic and stable at higher opacity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Interstellar Extinction, Polarization, and Grain Alignment in the Sh 2-185 (IC 59 and IC 63) Region

Optical and infrared continuum polarization from the interstellar medium is driven by radiative processes aligning the grains with the magnetic field. While a quantitative, predictive theory of radiative alignment torques (RATs) exists and has been extensively tested, several parameters of the theory remain to be fully constrained. In a recent paper, Medan & Andersson showed that the polarization efficiency (and therefore grain alignment efficiency) at different locations in the wall of the Local Bubble (LB) could be modeled as proportional to the integrated light intensity from the surrounding stars and OB associations. Here we probe that relationship at high radiation field intensities by studying the extinction and polarization in the two reflection nebulae IC 59 and IC 63 in the Sh 2-185 H ii region, illuminated by the B0 IV star γ Cassiopeia. We combine archival visual polarimetry with new seven-band photometry in the Vilnius system, to derive the polarization efficiency from the material. We find that the same linear relationship seen in the LB wall also applies to the Sh 2-185 region, strengthening the conclusion from the earlier study.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

H α emission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Constraint on Primordial B-modes from the First Flight of the Spider Balloon-borne Telescope

We present the first linear polarization measurements from the 2015 long-duration balloon flight of Spider, which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in the Spider bands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis of Spider and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and r parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, finding r < 0.11 and r < 0.19, respectively. Roughly half the uncertainty in r derives from noise associated with the template subtraction. New data at 280 GHz from Spider's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
DOE OSTI2020

From Interstellar Cloud to Star to Laboratory: Frontier HEDP Studies of Magnetized Colliding Plasma Flows with Strong Radiative Cooling (Final Report)

High-speed flows of gas are the norm in astrophysics where exploding stars, winds from the regions around black holes and streams of galactic material create “supersonic flows”. As the flows collide, strong shock waves are formed where gas is compressed and heated. The regions behind these shock waves are of great interest to scientists. New stars can form in these shock regions in the context of galaxies. In the context of planets like our Earth, the shocks waves come from the interaction of the solar wind with our planets’ magnetic field causing Aurora and placing astronauts and satellites in jeopardy. Thus, understanding supersonic flows of magnetized plasmas is essential for progress in many fields of astronomy and space physics (relevant to planets and the Earth). Remarkably, we can also create these kinds of flows in the laboratory using Z-pinches and laser driven experiments that include magnetic fields. The exploration of these forms of “High Energy Density Physics” (HEDP) flows is also critical in order to achieve Inertial Confinement Fusion, which is a step towards developing a sustainable source of Fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗