Photoelectric Heating Research Articles

Influence of the nano-grain depletion in photon-dominated regions

Context. The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer irradiated part of PDRs. Aims. We highlight the influence of nano-grain depletion on the gas physics and chemistry in the Horsehead nebula, a prototypical PDR. Methods. We used a model for atomic and molecular gas in PDRs, the Meudon PDR code, using diffuse ISM-like dust and Horsehead-like dust to study the influence of nano-grain depletion on the gas physics and chemistry, focusing on the impact on photoelectric heating and H2 formation and, therefore, on the H2 gas lines. Results. We find that nano-grain depletion in the Horsehead strongly affects gas heating through the photoelectric effect and thus the gas temperature and the H2 formation, hence the H → H2 position. Consequently, the first four pure rotational lines of H2 (e.g. 0–0 S(0), S(1), S(2), and S(3)) vary by a factor of 2 to 14. The 0–0 S(3) line that is often underestimated in models is underestimated even more when taking nano-grain depletion into account due to the decrease in gas heating through the photoelectric effect. This strongly suggests that our understanding of the excitation of H2 and/or of heating processes in the Horsehead, and more generally in PDRs, is still incomplete. Conclusions. Nano-grain depletion in the outer part of the Horsehead has a strong influence on several gas tracers that will be prominent in JWST observations of irradiated clouds. We therefore need to take this depletion into account in order to improve our understanding of the Horsehead, and more generally PDRs, and to contribute to the optimal scientific return of the mission.

C ii] and CO Emission along the Bar and Counter-arms of NGC 7479*

Abstract We present new SOFIA [C ii] and ALMA CO J=1→0 observations of the nearby asymmetric barred spiral galaxy NGC 7479. The data, which cover the whole bar of the galaxy and the counter-arms visible in the radio continuum, are analyzed in conjunction with a wealth of existing visible, infrared, radio, and X-ray data. As in most normal galaxies, the [C ii] emission is generally consistent with emission from cooling gas excited by photoelectric heating in photodissociation regions. However, anomalously high [C ii]/CO ratios are seen at the two ends of the counter-arms. Both ends show shell-like structures, possibly bubbles, in Hα emission. In addition, the southern end has [C ii]-to-infrared emission ratios inconsistent with normal star formation. Because there is little H i emission at this location, the [C ii] emission probably originates in warm shocked molecular gas heated by the interaction of the radio jet forming the counter-arms with the interstellar medium in the galaxy. At two other locations, the high [C ii]/CO ratios provide evidence for the existence of patches of CO-dark molecular gas. The [C ii] and CO observations also reveal resolved velocity components along the bar. In particular, the CO emission can be separated into two components associated with gas along the leading edge of the bar and gas trailing the bar. The trailing gas component that amounts to approximately 40% of the gas around the bar region may be related to a minor merger.

The sensitivity of stellar feedback to IMF averaging versus IMF sampling in galaxy formation simulations

ABSTRACT Galaxy formation simulations frequently use initial mass function (IMF) averaged feedback prescriptions, where star particles are assumed to represent single stellar populations that fully sample the IMF. This approximation breaks down at high mass resolution, where stochastic variations in stellar populations become important. We discuss various schemes to populate star particles with stellar masses explicitly sampled from the IMF. We use Monte Carlo numerical experiments to examine the ability of the schemes to reproduce an input IMF in an unbiased manner while conserving mass. We present our preferred scheme which can easily be added to pre-existing star formation prescriptions. We then carry out a series of high-resolution isolated simulations of dwarf galaxies with supernovae (SNe), photoionization, and photoelectric heating to compare the differences between using IMF averaged feedback and explicitly sampling the IMF. We find that if SNe are the only form of feedback, triggering individual SNe from IMF averaged rates gives identical results to IMF sampling. However, we find that photoionization is more effective at regulating star formation when IMF averaged rates are used, creating more, smaller H ii regions than the rare, bright sources produced by IMF sampling. We note that the increased efficiency of the IMF averaged feedback versus IMF sampling is not necessarily a general trend and may be reversed depending on feedback channel, resolution and other details. However, IMF sampling is always the more physically motivated approach. We conservatively suggest that it should be used for star particles less massive than $\sim 500\, \mathrm{M_\odot }$.

The anatomy of a star-forming galaxy II: FUV heating via dust

ABSTRACT Far-ultraviolet (FUV) radiation greatly exceeds UV, supernovae (SNe), and winds in the energy budget of young star clusters but is poorly modelled in galaxy simulations. We present results of the first isolated galaxy disc simulations to include photoelectric heating of gas via dust grains from FUV radiation self-consistently, using a ray-tracing approach that calculates optical depths along the source–receiver sightline. This is the first science application of the TREVR radiative transfer algorithm. We find that FUV radiation alone cannot regulate star formation. However, FUV radiation produces warm neutral gas and is able to produce regulated galaxies with realistic scale heights. FUV is also a long-range feedback and is more important in the outer discs of galaxies. We also use the superbubble feedback model, which depends only on the SN energy per stellar mass, is more physically realistic than common, parameter-driven alternatives and thus better constrains SN feedback impacts. FUV and SNe together can regulate star formation without producing too much hot ionized medium and with less disruption to the interstellar medium compared to SNe alone.

Photoelectric heating effects on the evolution of luminous disc galaxies

ABSTRACT Photoelectric heating (PEH) influences the temperature and density of the interstellar medium (ISM), potentially also affecting star formation. PEH is expected to have a stronger effect on massive galaxies, as they host larger dust reservoirs compared to dwarf systems. Accordingly, in this paper, we study PEH effects in Milky Way-like galaxies using a smoothed particle hydrodynamics code, which self-consistently implements the evolution of the gas, dust, and interstellar radiation field. Dust evolution includes dust formation by stars, destruction by SNe, and growth in dense media. We find that PEH suppresses star formation due to the excess heating that reduces the ISM density. This suppression is seen across the entire range of gas fractions, star-formation recipes, dust models, and PEH efficiencies investigated by our code. The suppression ranges from negligible values to approximately a factor of five depending on the specific implementation. Galaxy models having higher gas fractions experience higher star-formation suppression. The adopted dust model also alters the extent of star-formation suppression. Moreover, when PEH is switched on, galaxy models show higher gas outflow rates and have higher loading factors, indicative of enhanced SNe feedback. In gas-rich models (i.e. a gas fraction of 0.5), we also find that PEH suppresses the formation of disc clumps via violent disc instabilities, and thus suppresses bulge formation via clump migration to the central regions.

Factories of CO-dark gas: molecular clouds with limited star formation efficiencies by far-ultraviolet feedback

ABSTRACT The star formation in molecular clouds is inefficient. The ionizing extreme-ultraviolet radiation (hν ≥ 13.6 eV) from young clusters has been considered as a primary feedback effect to limit the star formation efficiency (SFE). Here, we focus on the effects of stellar far-ultraviolet (FUV) radiation (6 eV ≤ hν ≤ 13.6 eV) during the cloud disruption stage. The FUV radiation may further reduce the SFE via photoelectric heating, and it also affects the chemical states of the gas that is not converted to stars (‘cloud remnants’) via photodissociation of molecules. We have developed a one-dimensional semi-analytical model that follows the evolution of both the thermal and chemical structure of a photodissociation region (PDR) during the dynamical expansion of an H ii region. We investigate how the FUV feedback limits the SFE, supposing that the star formation is quenched in the PDR where the temperature is above a threshold value (e.g. 100 K). Our model predicts that the FUV feedback contributes to reduce the SFEs for massive (Mcl ≳ 105 M⊙) clouds with low surface densities ($\Sigma _{\rm cl}\lesssim 100~{\rm M}_\odot \, {\rm pc}^{-2}$). Moreover, we show that a large part of the H2 molecular gas contained in the cloud remnants should be ‘CO-dark’ under the FUV feedback for a wide range of cloud properties. Therefore, the dispersed molecular clouds are potential factories of CO-dark gas, which returns into the cycle of the interstellar medium.

Softening of the Metagalactic Ionizing Background from Internal He ii Absorption in Quasars

Abstract Quasars and other active galactic nuclei (AGN) are significant contributors to the metagalactic ionizing background at redshifts z < 3. Recent Hubble Space Telescope (HST)/Cosmic Origins Spectrograph (COS) composite spectra of AGN find a harder flux distribution in the Lyman continuum, (α s = 1.41 ± 0.15), compared to previous studies. This index appears to be inconsistent with observed He ii/H i absorption ratios (η) in the Lyα forest. We explore the effects of internal AGN absorption in the He ii (4 ryd) continuum using an analytic source-function model of the ionizing background in which the emissivity (j ν ) arises from quasars, reprocessed by the opacity (κ ν ) of the intervening Lyα forest and distinct AGN escape fractions and at 1 ryd and 4 ryd, respectively. We also examine H i and He ii photoelectric heating from the reprocessed background, whose spectral index (α b > α s ) depends on α s , the H i column density slope β, and the ratio . We compare the model to Lyα absorption lines of He ii and H i toward the quasar HE 2347−4342. Internal AGN absorption with but would increase the index by Δα b ≈ 0.3–1.0, corresponding to η = 60–200 for β ≈ 1.5–1.6, in agreement with HST/COS observations at z ≈ 2.5–2.9. The observed range of ratios, η < 200, constrains α b < 3.4 and . Individual AGN with softer spectra, α s > 1.7, and more internal He ii absorption could produce a few absorbers with η > 300, in proximity to AGN transverse to the sight line.

Measuring the Heating and Cooling of the Interstellar Medium at High Redshift: PAH and [C ii] Observations of the Same Star-forming Galaxies at z ∼ 2

Abstract Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ∼ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C ii] 158 μm fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C ii] in six z ∼ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C ii] emission, and constrain their relationship as a function of total infrared luminosity (L IR) and IR surface density. [C ii] emission is detected in one galaxy at high signal-to-noise (34σ), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C ii] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C ii]/L IR and PAH/L IR observed in the literature. However, the ratio of [C ii] to PAH emission at z ∼ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O i] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.

Spitzer’s perspective of polycyclic aromatic hydrocarbons in galaxies

Polycyclic aromatic hydrocarbon (PAH) molecules are abundant and widespread throughout the Universe, as revealed by their distinctive set of emission bands at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 μm, which are characteristic of their vibrational modes. They are ubiquitously seen in a wide variety of astrophysical regions, ranging from planet-forming disks around young stars to the interstellar medium of the Milky Way and other galaxies out to high redshifts at z ≳ 4. PAHs profoundly influence the thermal budget and chemistry of the interstellar medium by dominating the photoelectric heating of the gas and controlling the ionization balance. Here I review the current state of knowledge of the astrophysics of PAHs, focusing on their observational characteristics obtained from the Spitzer Space Telescope and their diagnostic power for probing the local physical and chemical conditions and processes. Special attention is paid to the spectral properties of PAHs and their variations revealed by the Infrared Spectrograph onboard Spitzer across a much broader range of extragalactic environments (for example, distant galaxies, early-type galaxies, galactic halos, active galactic nuclei and low-metallicity galaxies) than was previously possible with the Infrared Space Observatory or any other telescope facilities. Also highlighted is the relation between the PAH abundance and the galaxy metallicity established for the first time by Spitzer. Spitzer revealed the power of astrophysical polycyclic aromatic hydrocarbon molecules to probe the local physical and chemical conditions and processes, for example, establishing the relation between their abundance and galaxy metallicity for the first time.

Global simulations of galactic discs: violent feedback from clustered supernovae during bursts of star formation

ABSTRACT A suite of idealized, global, gravitationally unstable, star-forming galactic disc simulations with 2 pc spatial resolution, performed with the adaptive mesh refinement code ramses, is used in this paper to predict the emergent effects of supernova feedback. The simulations include a simplified prescription for the formation of single stellar populations of mass $\sim 100 \, {\rm M}_{\odot }$, radiative cooling, photoelectric heating, an external gravitational potential for a dark matter halo and an old stellar disc, self-gravity, and a novel implementation of supernova feedback. The results of these simulations show that gravitationally unstable discs can generate violent supersonic winds with mass-loading factors η ≳ 10, followed by a galactic fountain phase. These violent winds are generated by highly clustered supernovae exploding in dense environments created by gravitational instability, and they are not produced in simulation without self-gravity. The violent winds significantly perturb the vertical structure of the disc, which is later re-established during the galactic fountain phase. Gas resettles into a quasi-steady, highly turbulent disc with volume-weighted velocity dispersion $\sigma \gt 50 \, {\rm km\, s}^{-1}$. The new configuration drives weaker galactic winds with a mass-loading factor η ≤ 0.1. The whole cycle takes place in ≤10 dynamical times. Such high time variability needs to be taken into account when interpreting observations of galactic winds from starburst and post-starburst galaxies.

Disentangling Planets from Photoelectric Instability in Gas-rich Optically Thin Dusty Disks

Abstract Structures in circumstellar disks such as gaps and rings are often attributed to planets. This connection has been difficult to show unequivocally, as other processes may also produce these features. In particular, a photoelectric instability (PEI) has been proposed, operating in gas-rich optically thin disks, that generates structures predicted by planet–disk interactions. We examine the question of how to disentangle the planetary effects on disk structure from the effects of the PEI. We use the Pencil Code to perform 2D global hydrodynamical models of the dynamics of gas and dust in a thin disk with and without planetary perturbers. Photoelectric heating is modeled with an equation of state where pressure is proportional to dust surface density. The drag force on grains and its backreaction on the gas are included. Analyzing the situation without PEI, we find that gas–dust interactions alter the shape of the planetary gap from the dust-free case when the local dust-to-gas ratio ε approaches unity. This result also applies to primordial disks, because dust drifting inward accumulates at the edge of the planetary gap, and any initial dust-to-gas ratio eventually achieves ε = 1 if the dust reservoir is sufficient. We find a result particular to high dust-to-gas ratio disks as well: as dust drifts inward, the dust front becomes a sharp transition, and the backreaction triggers the Rossby wave instability. When PEI is included, we find that it obscures structures induced by planets unless the planet’s mass is sufficiently large to carve a noticeable gap. Specifically, the instability generates arcs and rings of regular spacing: a planet is discernible when it carves a dust gap wider than the wavelength of the PEI.

Using [C ii] 158 μm Emission from Isolated ISM Phases as a Star Formation Rate Indicator

Abstract The brightest observed emission line in many star-forming galaxies is the 158 μm line, making it detectable up to z ∼ 7. In order to better understand and quantify the emission as a tracer of star formation, the theoretical ratio between the 205 μm emission and the 158 μm emission has been employed to empirically determine the fraction of emission that originates from the ionized and neutral phases of the interstellar medium (ISM). Sub-kiloparsec measurements of the 158 μm and 205 μm lines in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the emission line as a star formation rate (SFR) indicator and a better understanding of the deficit are pursued. emission is also compared to polycyclic aromatic hydrocarbon (PAH) emission in these regions to compare photoelectric heating from PAH molecules to cooling by in the neutral and ionized phases of the ISM. We find that the emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the emission originating in the ionized phase of the ISM as an SFR indicator for the normal star-forming galaxies included in this sample.

Radiative stellar feedback in galaxy formation: Methods and physics

ABSTRACT Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly understood. We explore this using high-resolution, multifrequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultrafaint dwarf through Milky Way mass scales, including H+He photoionization; photoelectric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically thin) and moments-based M1 (exact when optically thick). The most important RFB channels on galaxy scales are photoionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (∼1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photoionization heating from the UV background suppressing accretion. In MW-mass galaxies, metagalactic backgrounds have negligible effects; but local photoionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other ‘rapid’ FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyperdense, bound star clusters. This makes star formation more violent and ‘bursty’ when SNe explode in these hyperclustered objects: thus, including RFB ‘smoothes’ SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, $\sim 10{{\ \rm per\ cent}}$ of RP in massive galaxies): absorption occurs primarily in ‘normal’ GMCs with AV ∼ 1.

Herschel spectroscopy of massive young stellar objects in the Magellanic Clouds

ABSTRACT We present Herschel Space Observatory Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE FTS) spectroscopy of a sample of 20 massive Young Stellar Objects (YSOs) in the Large and Small Magellanic Clouds (LMC and SMC). We analyse the brightest far-infrared (far-IR) emission lines, that diagnose the conditions of the heated gas in the YSO envelope and pinpoint their physical origin. We compare the properties of massive Magellanic and Galactic YSOs. We find that [O i] and [C ii] emission, that originates from the photo-dissociation region associated with the YSOs, is enhanced with respect to the dust continuum in the Magellanic sample. Furthermore the photoelectric heating efficiency is systematically higher for Magellanic YSOs, consistent with reduced grain charge in low metallicity environments. The observed CO emission is likely due to multiple shock components. The gas temperatures, derived from the analysis of CO rotational diagrams, are similar to Galactic estimates. This suggests a common origin to the observed CO excitation, from low-luminosity to massive YSOs, both in the Galaxy and the Magellanic Clouds. Bright far-IR line emission provides a mechanism to cool the YSO environment. We find that, even though [O i], CO, and [C ii] are the main line coolants, there is an indication that CO becomes less important at low metallicity, especially for the SMC sources. This is consistent with a reduction in CO abundance in environments where the dust is warmer due to reduced ultraviolet-shielding. Weak H2O and OH emission is detected, consistent with a modest role in the energy balance of wider massive YSO environments.

Photoevaporation of Molecular Gas Clumps Illuminated by External Massive Stars: Clump Lifetimes and Metallicity Dependence

Abstract We perform a suite of 3D radiation hydrodynamics simulations of photoevaporation of molecular gas clumps illuminated by external massive stars. We study the fate of solar-mass clumps and derive their lifetimes by varying the gas metallicity over a range of . Our simulations incorporate radiation transfer of far- and extreme-ultraviolet photons and follow atomic/molecular line cooling and dust–gas collisional cooling. Nonequilibrium chemistry is coupled with the radiative transfer and hydrodynamics in a self-consistent manner. We show that radiation-driven shocks compress gas clumps to have a volume that is set by the pressure equilibrium with the hot ambient gas. Radiative cooling enables metal-rich clumps to condense and have small surface areas where photoevaporative flows are launched. For our fiducial setup with an O-type star at a distance of 0.1 pc, the resulting photoevaporation rate is as small as for metal-rich clumps, but it is larger for metal-poor clumps that have larger surface areas. The clumps are continuously accelerated away from the radiation source by the so-called rocket effect and can travel over ∼1 pc within the lifetime. We also study the photoevaporation of clumps in a photodissociation region. Photoelectric heating is inefficient for metal-poor clumps that contain a smaller amount of grains, and thus they survive for over 105 yr. We conclude that the gas metallicity strongly affects the clump lifetime and thus determines the strength of feedback from massive stars in star-forming regions.

Thermal Phases of the Neutral Atomic Interstellar Medium from Solar Metallicity to Primordial Gas

Abstract We study the thermal structure of the neutral atomic (H i) interstellar medium across a wide range of metallicities, from supersolar down to vanishing metallicity, and for varying UV intensities and cosmic-ray (CR) ionization rates. We calculate self-consistently the gas temperature and species abundances (with a special focus on the residual H2), assuming a thermal and chemical steady state. For solar metallicity, Z′ ≡ 1, we recover the known result that there exists a pressure range over which the gas is multiphased, with the warm (∼104 K, warm neutral medium (WNM)) and cold (∼100 K, cold neutral medium (CNM)) phases coexisting at the same pressure. At a metallicity Z′ ≈ 0.1, the CNM is colder (compared to Z′ = 1) due to the reduced efficiency of photoelectric heating. For Z′ ≲ 0.1, CR ionization becomes the dominant heating mechanism and the WNM-to-CNM transition shifts to ever-increasing pressure/density as the metallicity is reduced. For metallicities Z′ ≲ 0.01, H2 cooling becomes important, lowering the temperature of the WNM (down to ≈600 K), and smoothing out the multiphase phenomenon. At vanishing metallicities, H2 heating becomes effective and the multiphase phenomenon disappears entirely. We derive analytic expressions for the critical densities for the warm-to-cold phase transition in the different regimes, and the critical metallicities for H2 cooling and heating. We discuss potential implications on the star formation rates of galaxies and self-regulation theories.

AGN-driven quenching of satellite galaxies

ABSTRACT We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 haloes with masses ranging from 1012 to 1013.4$\rm {M}_{\odot}$ at z = 0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Strömgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z = 2. The mean gas-rich fraction of satellites at z = 0 decreases from 15 per cent in the noAGN simulation to 5 per cent in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z = 1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

Carbon radio recombination lines from gigahertz to megahertz frequencies towards Orion A

Context. The combined use of carbon radio recombination lines (CRRLs) and the 158μm-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158μm-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15–230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion A (M 42). We complement these observations with SOFIA observations of the 158μm-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C+/C/CO interface. By modeling the emission of the 158μm-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure >5 × 107 K cm−3 and a radiation field G0 ≈ 105 close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158μm-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C+/C/CO interface.

Simulating an isolated dwarf galaxy with multichannel feedback and chemical yields from individual stars

In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium.

Radiation Hydrodynamics Simulations of Photoevaporation of Protoplanetary Disks. II. Metallicity Dependence of UV and X-Ray Photoevaporation

Abstract We perform a suite of radiation hydrodynamics simulations of photoevaporating disks, varying the metallicity in a wide range of . We follow the disk evolution for over ∼5000 years by solving hydrodynamics, radiative transfer, and nonequilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultraviolet (UV) radiation only. At subsolar metallicities in the range of , the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. In contrast, the photoevaporation rate decreases at even lower metallicities of , because dust–gas collisional cooling remains efficient compared to far-UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at . Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far-UV radiation is strengthened by the X-rays and enhances the disk photoevaporation.