Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 73 records · Page 4

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral solar irradiance on inclined surfaces: A fast Monte Carlo approach

Estimating spectral plane-of-array (POA) solar irradiance on inclined surfaces is an important step in the design and performance evaluation of both photovoltaic and concentrated solar plants. This work introduces a fast, line-by-line spectral, Monte Carlo (MC) radiative transfer model approach to simulate anisotropic distributions of shortwave radiation through the atmosphere as photon bundles impinge on inclined surfaces. This fast Monte Carlo approach reproduces the angular distribution of solar irradiance without the undesirable effects of spatial discretization and thus computes detailed POA irradiance values on surfaces at any orientation and also when surfaces are subjected to the anisotropic ground and atmospheric scattering. Polarization effects are also easily incorporated into this approach that can be considered as direct numerical simulation of the physics involved. Here, we compare our Monte Carlo radiative transfer model with the most widely used empirical transposition model, Perez4, under various conditions. The results show that the Perez4 model reproduces the more detailed Monte Carlo simulations with less than 10% deviation under clear skies for all relevant surface tilt and azimuth angles. When optically thin clouds are present, observed deviations are larger, especially when the receiving surface is strongly tilted. Deviations are also observed for large azimuth angle differences between the receiving surface and the solar position. When optically thick clouds are present, the two models agree within 15% deviation for nearly all surface orientation and tilt angles. The overall deviations are smaller when compared with cases for optically thin clouds. The Perez4 model performs very well (6.0% deviation) in comparison with the detailed MC simulations for all cases, thus validating its widespread use for practical solar applications. When detailed atmospheric profiles and cloud optical properties are available, the proposed fast Monte Carlo radiative model reproduces accurate spectral and angular POA irradiance levels for various atmospheric and cloud cover conditions, surface orientations, and different surface and ground properties.

42 ENGINEERING↗

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↗

A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model

Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model (Final Report)

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗

A detailed experimental and modeling comparison of molecular radiative heat loss in a spark-ignition engine

Radiative heat transfer has been chiefly considered negligible in internal combustion engines, except for Diesel engines where soot radiation was recognized as a significant radiative transfer source. Only more recently, detailed simulations have shown that molecular radiation can be substantial as well. In extension to this, molecular radiative heat transfer can reach detectable levels of about 5–10% of the total heat transfer in spark-ignited engines. The broadband radiative nature of the significant emitting molecules, carbon dioxide and water, makes it necessary to address whether radiative trapping plays a substantial role in either total heat loss or energy redistribution within the cylinder. An experimental setup that allows measurements of the infrared emissions at a 2-crank-angle-degree resolution was developed to determine the importance of radiative trapping by carbon dioxide and water. Measurements were conducted in the well-characterized and documented TCC-III engine at the University of Michigan. The engine operated on a stoichiometric propane/air mixture at a speed of 1,300 RPM. Large Eddy Simulations with added line-by-line photon Monte-Carlo simulations of the molecular radiation were an integral part of this study to plan and devise the measurements. Post-processing of the simulation data included an accurate representation of the experimental volume from where infrared signals are collected. Additionally, the experimental data were used for validation of the photon Monte-Carlo simulations. Joint analysis of experimental and simulated spectra allowed quantifying the significance of radiative trapping. Results show the significant role of radiative trapping when predicting radiative heat transfer in the TCC-III engine.

42 ENGINEERING↗

Photon (Non)Conservation in the Reduced Speed of Light Approximation and How to (Almost) Fix It

The "Reduced Speed of Light" (RSL) approximation is commonly used to speed up radiative transfer calculations in cosmological simulations. However, it has been shown previously that the RSL approximation leads to photon non-conservation when the radiation field is rapidly evolving in time. I show that these missing photons can be counted exactly for some numerical schemes. Adding them back into a simulation, however, is a much harder task. I show one example of such a scheme, which achieves sub-percent accuracy on simple tests. Unfortunately, the scheme performs much worse on semi-realistic simulations of cosmic reionization, leading to a faster overlap and significant errors in the point-wise comparison of the RSL radiation field with the reference simulation that maintains the full speed of light for the radiative transfer.

Gnedin, Nickolay Y. [Fermilab; Chicago U., KICP; C↗

Probing circular polarization and magnetic field structure in active galactic nuclei

The composition and magnetic field morphology of relativistic jets can be studied using circular polarization (CP) measurement. Recent three-dimensional relativistic magnetohydrodynamic (3D RMHD) simulations coupled with radiative transfer (RT) calculations make strong predictions about the level (and morphology) of the jet’s CP emission. These simulations show that the sign of CP and the electric vector position angle (EVPA) are both sensitive to the jet’s magnetic field morphology within the radio core. We probe this theory by exploring whether the jet’s radio core EVPA orientation is consistent with the observed sign of the core CP in deep full-track polarimetric observations. Based on a selection of sources from earlier MOJAVE observations, we aim to probe the nature of linear polarization (LP) and CP in the innermost regions of jets from a small sample of nine blazars. This sample includes sources that have exhibited: (i) positive CP; (ii) negative CP; or (iii) positive and negative CP simultaneously in the radio core region. By coupling deep polarimetric observations of a carefully selected sample of blazars with state-of-the-art RMHD and RT calculations, we hope to gain a deeper understanding of the physics of blazar jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗

Using soil library hyperspectral reflectance and machine learning to predict soil organic carbon: Assessing potential of airborne and spaceborne optical soil sensing

Soil organic carbon (SOC) is a key variable to determine soil functioning, ecosystem services, and global carbon cycles. Spectroscopy, particularly optical hyperspectral reflectance coupled with machine learning, can provide rapid, efficient, and cost-effective quantification of SOC. However, how to exploit soil hyperspectral reflectance to predict SOC concentration, and the potential performance of airborne and satellite data for predicting surface SOC at large scales remain relatively underknown. Here, this study utilized a continental-scale soil laboratory spectral library (37,540 full-pedon 350–2500 nm reflectance spectra with SOC concentration of 0–780 g·kg –1 across the US) to thoroughly evaluate seven machine learning algorithms including Partial-Least Squares Regression (PLSR), Random Forest (RF), K-Nearest Neighbors (KNN), Ridge, Artificial Neural Networks (ANN), Convolutional Neural Networks (CNN), and Long Short-Term Memory (LSTM) along with four preprocessed spectra, i.e. original, vector normalization, continuum removal, and first-order derivative, to quantify SOC concentration. Furthermore, by using the coupled soil-vegetation-atmosphere radiative transfer model, we simulated twelve airborne and spaceborne hyper/multi-spectral remote sensing data from surface bare soil laboratory spectra to evaluate their potential for estimating SOC concentration of surface bare soils. Results show that LSTM achieved best predictive performance of quantifying SOC concentration for the whole data sets (R 2 = 0.96, RMSE = 30.81 g·kg –1 ), mineral soils (SOC ≤ 120 g·kg –1 , R 2 = 0.71, RMSE = 10.60 g·kg –1 ), and organic soils (SOC > 120 g·kg –1 , R 2 = 0.78, RMSE = 62.31 g·kg –1 ). Spectral data preprocessing, particularly the first-order derivative, improved the performance of PLSR, RF, Ridge, KNN, and ANN, but not LSTM or CNN. We found that the SOC models of mineral and organic soils should be distinguished given their distinct spectral signatures. Finally, we identified that the shortwave infrared is vital for airborne and spaceborne hyperspectral sensors to monitor surface SOC. This study highlights the high accuracy of LSTM with hyperspectral/multispectral data to mitigate a certain level of noise (soil moisture <0.4 m 3 ·m –3 , green leaf area < 0.3 m 2 ·m –2 , plant residue <0.4 m 2 ·m –2 ) for quantifying surface SOC concentration. Forthcoming satellite hyperspectral missions like Surface Biology and Geology (SBG) have a high potential for future global soil carbon monitoring, while high-resolution satellite multispectral fusion data can be an alternative.

54 ENVIRONMENTAL SCIENCES↗

An asymmetric electron-scattering photosphere around optical tidal disruption events

We report a star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events have now been identified in optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three tidal disruption events. The continuum polarization appears independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code, we simulate the polarization degree as a function of disk mass and viewing angle and we show that the predicted levels are compatible with the observations for extended reprocessing envelopes of ~1,000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for tidal disruption event models, and opens an important new line of exploration in the study of tidal disruption events.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of Aerosol and Cloud on the Solar Energy Potential over the Central Gangetic Himalayan Region

We examine the impact of atmospheric aerosols and clouds on the surface solar radiation and solar energy at Nainital, a high-altitude remote location in the central Gangetic Himalayan region (CGHR). For this purpose, we exploited the synergy of remote-sensed data in terms of ground-based AERONET Sun Photometer and satellite observations from the MODerate Resolution Imaging Spectroradiometer (MODIS) and the Meteosat Second Generation (MSG), with radiative transfer model (RTM) simulations and 1 day forecasts from the Copernicus Atmosphere Monitoring Service (CAMS). Clouds and aerosols are one of the most common sources of solar irradiance attenuation and hence causing performance issues in the photovoltaic (PV) and concentrated solar power (CSP) plant installations. The outputs of RTM results presented with high accuracy under clear, cloudy sky and dust conditions for global horizontal (GHI) and beam horizontal irradiance (BHI). On an annual basis the total aerosol attenuation was found to be up to 105 kWh m –2 for the GHI and 266 kWh m –2 for BHI, respectively, while the cloud effect is much stronger with an attenuation of 245 and 271 kWh m –2 on GHI and BHI. The results of this study will support the Indian solar energy producers and electricity handling entities in order to quantify the energy and financial losses due to cloud and aerosol presence.

54 ENVIRONMENTAL SCIENCES↗

Axisymmetric Radiative Transfer Models of Kilonovae

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]

Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling CN Zeeman effect observations of the envelopes of a low-mass protostellar disc and a massive protostar

We use the POLARIS radiative transfer code to produce simulated circular polarization Zeeman emission maps of the cyanide (CN) J = 1–0 molecular line transition for two types of protostellar envelope magnetohydrodynamic simulations. Our first model is a low-mass disc envelope system (box length L = 200 au), and our second model is the envelope of a massive protostar (L = 104 au) with a protostellar wind and a CN-enhanced outflow shell. We compute the velocity-integrated Stokes I and V, as well as the implied V/I polarization percentage, for each detector pixel location in our simulated emission maps. Our results show that both types of protostellar environments are in principle accessible with current circular polarization instruments, with each containing swaths of envelope area that yield percentage polarizations that exceed the 1.8 percent nominal sensitivity limit for circular polarization experiments with the Atacama Large Millimeter/submillimeter Array. In both systems, high-polarization (≳1.8 percent) pixels tend to lie at an intermediate distance away from the central star and where the line-centre opacity of the CN emission is moderately optically thin (τLC ~ 0.1–1). Furthermore, our computed V/I values scale roughly with the density-weighted mean line-of-sight magnetic field strength, indicating that Zeeman observations can effectively diagnose the strength of envelope-scale magnetic fields. We also find that pixels with large V/I are preferentially co-located where the absolute value of the velocity-integrated V is also large, suggesting that locations with favourable percentage polarization are also favourable in terms of raw signal.

79 ASTRONOMY AND ASTROPHYSICS↗

Artificial Intelligence Assisted Inversion (AIAI): Quantifying the Spectral Features of 56 Ni of Type Ia Supernovae

Abstract Following our previous study of Artificial Intelligence Assisted Inversion (AIAI) of supernova analyses, we train a set of deep neural networks based on the 1D radiative transfer code TARDIS to simulate the optical spectra of Type Ia supernovae (SNe Ia) between 10 and 40 days after the explosion. The neural networks are applied to derive the mass of 56 Ni in velocity ranges above the photosphere for a sample of 124 well-observed SNe Ia in the TARDIS model context. A subset of the SNe have multi-epoch observations for which the decay of the radioactive 56 Ni can be used to test the AIAI quantitatively. The 56 Ni mass derived from AIAI using the observed spectra as inputs for this subset agrees with the radioactive decay rate of 56 Ni. AIAI reveals that a spectral signature near 3890 Å is related to the Niii4067Å line, and the 56 Ni mass deduced from AIAI is found to be correlated with the light-curve shapes of SNe Ia, with SNe Ia with broader light curves showing larger 56 Ni mass in the envelope above the photosphere. AIAI enables spectral data of SNe to be quantitatively analyzed under theoretical frameworks based on well-defined physical assumptions.

Astronomy & Astrophysics↗

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗

Spatial Decorrelation of Young Stars and Dense Gas as a Probe of the Star Formation–Feedback Cycle in Galaxies

The spatial decorrelation of dense molecular gas and young stars observed on ≲1 kpc scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. In this work, we explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC 300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC 300, and we show that this agreement is due to different aspects of feedback, including H 2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affects the correlation on ≲100 pc scales, while supernovae play a significant role on ≳ 100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per free fall time, ϵ ff , which provides a strong observational constraint on ϵ ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗