Models Of Interstellar Medium Research Articles

Neural network-based emulation of interstellar medium models

Context. The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too timeconsuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. Aims. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models for use in inference procedures. Methods. These emulators are defined with artificial neural networks (ANNs) with adapted architectures and are fitted using regression strategies instead of interpolation methods. The specificities inherent in ISM models need to be addressed to design and train adequate ANNs. Indeed, such models often predict numerous observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities and multistabilities. We propose applying five strategies to address these characteristics: (1) an outlier removal procedure; (2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; (3) a dimension reduction technique that enables us to adequately size the network architecture; (4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and (5) a dense architecture to ease the learning of simpler relations between line intensities and physical parameters. Results. We compare the proposed ANNs with four standard classes of interpolation methods, nearest-neighbor, linear, spline, and radial basis function (RBF), to emulate a representative ISM numerical model known as the Meudon PDR code. Combinations of the proposed strategies produce networks that outperform all interpolation methods in terms of accuracy by a factor of 2 in terms of the average error (reaching 4.5% on the Meudon PDR code) and a factor of 3 for the worst-case errors (33%). These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. Conclusions. This work will enable efficient inferences on wide-field multiline observations of the ISM.

Topological models to infer multiphase interstellar medium properties

Context.Spectroscopic observations of high-redshift galaxies slowly reveal the same complexity of the interstellar medium (ISM) as expected from resolved observations in nearby galaxies. While providing, in principle, a wealth of diagnostics concerning galaxy evolution, star formation, or the nature and influence of compact objects, such high-zspectra are often spatially and spectrally unresolved, and inferring reliable diagnostics represents a major obstacle. Bright, nearby, unresolved galaxies observed in the optical and infrared domains provide many constraints to design methods to infer ISM properties, but they have so far been limited to deterministic methods and/or with simple topological assumptions (e.g., single 1D model).Aims.It is urgent to build upon previous ISM multiphase and multicomponent methods by using a probabilistic approach that makes it possible to derive probability density functions for relevant parameters while also enabling a large number of free parameters with potential priors. The goal is to provide a flexible statistical framework that is agnostic to the model grid and that considers either a few discrete components defined by their parameter values and/or statistical distributions of parameters. In this paper, we present a first application with the objective to infer probability distributions of several physical parameters (e.g., the mass of H0, H2, escape fraction of ionizing photons, and metallicity) for the star-forming regions of the metal-poor dwarf galaxy I Zw 18 in order to confirm the low molecular gas content and high escape fraction of ionizing photons from Hiiregions.Methods.We present a Bayesian approach to model a suite of spectral lines using a sequential Monte Carlo method provided by the Python package PyMC which combines several concepts such as tempered likelihoods, importance sampling, and independent Metropolis-Hastings chains. The algorithm, provided by the associated code MULTIGRIS, accepts a few components which can be represented as sectors around one or several stellar clusters, or continuous (e.g., power-law, normal) distributions for any given parameter. We applied this approach to a grid of models calculated with the photoionization and photodissociation code Cloudy in order to produce topological models of I Zw 18.Results.The statistical framework we present makes it possible to consider a large number of spectroscopic tracers, with the extinction and systematic uncertainties as potential additional random variables. We applied this technique to the galaxy I Zw 18 in order to reproduce and go beyond previous topological models specifically tailored to this object. While our grid is designed for global properties of low-metallicity star-forming galaxies, we were able to calculate accurate values for the metallicity, number of ionizing photons, masses of ionized and neutral hydrogen, as well as the dust mass and the dust-to-gas mass ratio in I Zw 18. We find a relatively modest amount of H2(~105M⊙) which is predominantly CO-dark and traced by C+rather than C0. Nevertheless, more than 90% of the [Cii] emission is associated with the neutral atomic gas. Our models confirm the necessity to include an X-ray source with an inferred luminosity in good agreement with direct X-ray observations. Finally, we investigate the escape fraction of ionizing photons for different energy ranges. While the escape fraction for the main Hiiregion lies around 50–65%, we show that most of the soft X-ray photons are able to escape and may play a role in the ionization and heating of the circumgalactic or intergalactic medium.Conclusions.Multicomponent ISM models associate a complex enough distribution of matter and phases with a simple enough topological description to be constrained with probabilistic frameworks. Despite ignoring effects such as reflected light, the diffuse radiation field, or ionization by several non-cospatial sources, they remain well adapted to individual Hiiregions and to star-forming galaxies dominated by one or a few Hiiregions, and the improvement due to the combination of several components largely compensates for other secondary effects.

Simulations of black hole fueling in isolated and merging galaxies with an explicit, multiphase ISM

ABSTRACT We study gas inflows on to supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighbourhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions ≲3 × 103 M⊙. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM.

Thermal emission from bow shocks

The nearby, massive, runaway star ζ Ophiuchi has a large bow shock detected in optical and infrared light and, uniquely among runaway O stars, diffuse X-ray emission detected from the shocked stellar wind. Here we make the first detailed computational investigation of the bow shock of ζ Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare the detected X-ray emission with simulated emission maps. We reanalysed archival Chandra observations of the thermal diffuse X-ray emission from the shocked wind region of the bow shock, finding total unabsorbed X-ray flux in the 0.3–2keV band corresponding to a diffuse X-ray luminosity of LX = 2.33−1.54+1.12 × 1029 erg s−1, consistent with previous work. The diffuse X-ray emission arises from the region between the star and the bow shock. Three-dimensional magnetohydrodyanmic simulations were used to model the interaction of the star’s wind with a uniform interstellar medium (ISM) using a range of stellar and ISM parameters motivated by observational constraints. Synthetic infrared, Hα, soft X-ray, emission measure, and radio 6 GHz emission maps were generated from three simulations, for comparison with the relevant observations. Simulations where the space velocity of ζ Ophiuchi has a significant radial velocity produce infrared emission maps with the opening angle of the bow shock in better agreement with observations than for the case where motion is fully in the plane of the sky. All three simulations presented here have X-ray emission fainter than observed, in contrast to results for NGC 7635. The simulation with the highest pressure has the closest match to X-ray observations, with a flux level within a factor of two of the observational lower limit, and emission weighted temperature of log10(TA/K) = 6.4, although the morphology of the diffuse emission appears somewhat different. The observed X-ray emission is from a filled bubble that is brightest near the star, whereas simulations predict brightening towards the contact discontinuity as density increases. This first numerical study of the bow shock and wind bubble around ζ Ophiuchi uses a relatively simple model of a uniform ISM and ideal-magnetohydrodynamics, and can be used as a basis for comparing results from models incorporating more physical processes, or higher resolution simulations that may show more turbulent mixing.

A Mesofractal Model of Interstellar Cloudiness

The interstellar medium (ISM), serving as a background for the propagation of cosmic rays (CRs) and other information carriers, has a complex ragged structure. Being chaotically scattered over interstellar space, together with the magnetic field perturbations frozen in them, CRs are connected with each other by a network of magnetic field lines creating long-range correlations of a power-law type, similar to those observed in the spatial distribution of galaxies. These lines solving interstellar transfer problems require the choice of an ISM model, adequately and concisely representing their statistical properties. This article discusses one such model, the Uchaikin–Zolotarev model: a four-parameter approximation of the power spectrum spatial correlations, derived from the generalized Ornstein–Zernike equation. The numerical analysis confirmed that this approximation satisfactorily agrees with the numerical data obtained in the quasi-linear model of plasma turbulence.

Effects of Forcing on Shocks and Energy Dissipation in Interstellar and Intracluster Turbulences

Observations indicate that turbulence in the interstellar medium (ISM) is supersonic (M turb ≫ 1) and strongly magnetized (β ∼ 0.01–1), while in the intracluster medium (ICM) it is subsonic (M turb ≲ 1) and weakly magnetized (β ∼ 100). Here, M turb is the turbulent Mach number and β is the plasma beta. We study the properties of shocks induced in these disparate environments, including the distribution of the shock Mach number, M s , and the dissipation of the turbulent energy at shocks, through numerical simulations using a high-order, accurate code based on the weighted essentially nonoscillatory scheme. In particular, we investigate the effects of different modes of the forcing that drives turbulence: solenoidal, compressive, and a mixture of the two. In ISM turbulence, while the density distribution looks different with different forcings, the velocity power spectrum, P v , on small scales exhibits only weak dependence. Hence, the statistics of shocks depend weakly on forcing either. In the ISM models with M turb ≈ 10 and β ∼ 0.1, the fraction of the turbulent energy dissipated at shocks is estimated to be ∼15%, not sensitive to the forcing mode. In contrast, in ICM turbulence, P v as well as the density distribution show strong dependence on forcing. The frequency and average Mach number of shocks are greater for compressive forcing than for solenoidal forcing; so is the energy dissipation. The fraction of the ensuing shock dissipation is in the range of ∼10%–35% in the ICM models with M turb ≈ 0.5 and β ∼ 106. The rest of the turbulent energy should be dissipated through turbulent cascade.

Detecting FAST RADIO BURSTS in the MILKY WAY

ABSTRACT Fast radio bursts (FRBs) are highly energetic transient events with duration of order of microseconds to milliseconds and of unknown origin. They are known to lie at cosmological distances, through localization to host galaxies. Recently, an FRB-like event was seen from the Milky Way magnetar SGR 1935+2154 by the CHIME and Survey for Transient Astronomical Radio Emission 2 (STARE2) telescopes. This is the only magnetar that has produced FRB events in our Galaxy. Finding similar events in the Milky Way is of great interest to understanding FRB progenitors. Such events will be strongly affected by the turbulent interstellar medium (ISM) in the Milky Way, their intrinsic energy distribution, and their spatial locations within the plane of the Milky Way. We examine these effects using models for the distribution of electrons in the ISM to estimate the dispersion measure and pulse scattering of mock events, and a range of models for the spatial distribution and luminosity functions, including models motivated by the spatial distribution of the Milky Way’s magnetars. We evaluate the fraction of FRB events in the Milky Way that are detectable by STARE2 for a range of ISM models, spatial distributions, and burst luminosity functions. In all the models examined, only a fraction of burst events are detectable, mainly due to the scattering effects of the ISM. We find that Galactic Radio Explorer, a proposed all-sky experiment, could increase the detection rate of Milky Way FRB events by an order of magnitude, depending on assumptions made about the luminosity function and scale height of the FRBs.

LYRA I: Simulating the multi-phase ISM of a dwarf galaxy with variable energy supernovae from individual stars

Abstract We introduce the LYRA project, a new high resolution galaxy formation model built within the framework of the cosmological hydro-dynamical moving mesh code arepo. The model resolves the multi-phase interstellar medium down to 10 K. It forms individual stars sampled from the initial mass function (IMF), and tracks their lifetimes and death pathways individually. Single supernova (SN) blast waves with variable energy are followed within the hydrodynamic calculation to interact with the surrounding interstellar medium (ISM). In this paper, we present the methods and apply the model to a 1010 M⊙ isolated halo. We demonstrate that the majority of supernovae are Sedov-resolved at our fiducial gas mass resolution of 4 M⊙. We show that our SN feedback prescription self-consistently produces a hot phase within the ISM that drives significant outflows, reduces the gas density and suppresses star formation. Clustered SN play a major role in enhancing the effectiveness of feedback, because the majority of explosions occur in low density material. Accounting for variable SN energy allows the feedback to respond directly to stellar evolution. We show that the ISM is sensitive to the spatially distributed energy deposition. It strongly affects the outflow behaviour, reducing the mass loading by a factor of 2 − 3, thus allowing the galaxy to retain a higher fraction of mass and metals. LYRA makes it possible to use a comprehensive multi-physics ISM model directly in cosmological (zoom) simulations of dwarf and higher mass galaxies.

VALES

Context. Spatially resolved observations of the ionized and molecular gas are critical for understanding the physical processes that govern the interstellar medium (ISM) in galaxies. The observation of starburst systems is also important as they present extreme gas conditions that may help to test different ISM models. However, matched resolution imaging at ∼kpc scales for both ISM gas phases are usually scarce, and the ISM properties of starbursts still remain poorly understood. Aims. We aim to study the morpho-kinematic properties of the ionized and molecular gas in three dusty starburst galaxies at z = 0.12−0.17 to explore the relation between molecular ISM gas phase dynamics and the star-formation activity. Methods. We employ two-dimensional dynamical modelling to analyse Atacama Large Millimeter/submillimiter Array CO(1–0) and seeing-limited Spectrograph for INtegral Field Observations in the Near Infrared Paschen-α (Paα) observations, tracing the molecular and ionized gas morpho-kinematics at ∼kpc-scales. We use a dynamical mass model, which accounts for beam-smearing effects, to constrain the CO-to-H2 conversion factor and estimate the molecular gas mass content. Results. One starburst galaxy shows irregular morphology, which may indicate a major merger, while the other two systems show disc-like morpho-kinematics. The two disc-like starbursts show molecular gas velocity dispersion values comparable with those seen in local luminous and ultra luminous infrared galaxies but in an ISM with molecular gas fraction and surface density values in the range of the estimates reported for local star-forming galaxies. We find that these molecular gas velocity dispersion values can be explained by assuming vertical pressure equilibrium. We also find that the star-formation activity, traced by the Paα emission line, is well correlated with the molecular gas content, suggesting an enhanced star-formation efficiency and depletion times of the order of ∼0.1−1 Gyr. We find that the star-formation rate surface density (ΣSFR) correlates with the ISM pressure set by self-gravity (Pgrav) following a power law with an exponent close to 0.8. Conclusions. In dusty disc-like starburst galaxies, our data support the scenario in which the molecular gas velocity dispersion values are driven by the ISM pressure set by self-gravity and are responsible for maintaining the vertical pressure balance. The correlation between ΣSFR and Pgrav suggests that, in these dusty starbursts galaxies, the star-formation activity arises as a consequence of the ISM pressure balance.

Lyα Radiative Transfer: Monte Carlo Simulation of the Wouthuysen–Field Effect

Abstract A three-dimensional Monte Carlo Lyα radiative transfer (RT) code, named LaRT, is developed to study the Lyα RT and the Wouthuysen–Field (WF) effect. Using the code, we calculate the line profile of Lyα radiation within the multiphase interstellar medium (ISM), with a particular emphasis on gas at low densities. We show that the WF effect is in action: the central portion of the line profile tends to approach a small slice of the Planck function with a color temperature equal to the kinetic temperature of the gas, even in a system with an optical thickness as low as τ 0 ≈ 100–500. We also investigate the effects of the turbulent motion of the ISM on the emergent Lyα spectrum and color temperature. The turbulent motion broadens the emergent spectrum, as generally expected, but the color temperature is not affected by the turbulent motion in typical astrophysical environments. We utilize two multiphase ISM models, appropriate for the vicinity of the Sun, to calculate the 21 cm spin temperature of neutral hydrogen, including excitation via the Lyα resonant scattering. The first ISM model is a simple clumpy model, while the second is a self-consistent magnetohydrodynamics simulation model using the TIGRESS framework. Lyα photons originating from both H ii regions and the collisionally cooling gas are taken into account. We find that the Lyα radiation field is, in general, likely to be strong enough to bring the 21 cm spin temperature of the warm neutral medium close to the kinetic temperature. The escape fraction of Lyα in our ISM models is estimated to be ≈7%–20%.

Electron Density Structure of the Local Galactic Disk

Abstract Pulsar dispersion measures (DMs) have been used to model the electron density of the interstellar medium (ISM) in the Galactic disk as a plane-parallel medium, despite significant scatter in the DM-distance distribution and strong evidence for inhomogeneities in the ISM. We use a sample of pulsars with independent distance measurements to evaluate a model of the local ISM in the thick disk of the Galaxy that incorporates turbulent fluctuations, clumps, and voids in the electron density. The latter two components are required because about one-third of the lines of sight are discrepant from a strictly plane-parallel model. A likelihood analysis for smooth components of the model yields a scale height kpc and a mid-plane density n 0 = 0.015 ± 0.001 cm−3. The scatter in the DM-distance distribution is dominated by clumps and voids but receives significant contributions from a broad spectrum of density fluctuations, such as a Kolmogorov spectrum. The model is used to identify lines of sight with outlier values of DM. Three of these pulsars, J1614−2230, J1623−0908, and J1643−1224, lie behind known H ii regions, and the electron density model is combined with Hα intensity data to constrain the filling factors and other substructure properties of the H ii regions (Sh 2–7 and Sh 2–27). Several pulsars also exhibit enhanced DM fluctuations that are likely caused by an intersection of their lines of sight with the superbubble GSH 238+00+09.

ART2: a 3D parallel multiwavelength radiative transfer code for continuum and atomic and molecular lines

ABSTRACT ART2 is a 3D multiwavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The original ART2 has been extensively applied to hydrodynamics simulations to study panchromatic properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART2 can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [C ii], [N ii], [O iii], Ly α, and CO. These essential features, together with the multiphase ISM model and the adaptive grid, make ART2 a multipurpose code to study multiwavelength properties of a wide range of astrophysical systems from planetary discs to large-scale structures. To demonstrate the capability of the new ART2, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are vital for a broad array of studies. By enabling direct comparison between numerical simulations and multiband observations, ART2 provides a crucial theoretical framework for the understanding of existing and future surveys, and the synergy between multiband galaxy surveys and line intensity mappings. Therefore, ART2 is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

Clumpy galaxies in cosmological simulations: the effect of ISM model

ABSTRACT We study influence by models of interstellar medium (ISM) on properties of galaxies in cosmological simulations. We examine three models widely used in previous studies. The ISM models impose different equations of state on dense gas. Using zoom-in simulations, we demonstrate that switching the ISM models can control formation of giant clumps in massive discs at redshifts z ∼ 1–2, while their initial conditions and the other settings such as stellar feedback are unchanged. Thus, not only feedback but ISM models can also be responsible for clumpy morphologies of simulated galaxies. We find, however, that changing the ISM models hardly affects global properties of galaxies, such as the total stellar and gas masses, star formation rate, metallicity, and stellar angular momentum, irrespective of the significant difference of clumpiness; namely the ISM models only change clumpiness of discs. In addition, our approach provides a test to investigate impact by clump formation on the evolution of disc galaxies using the same initial conditions and feedback. We find that clump formation does not significantly alter the properties of galaxies and therefore could not be the causes of starburst or quenching.

No “Maunder Minimum” Candidates in M67: Mitigating Interstellar Contamination of Chromospheric Emission Lines

Abstract The solar analogs of M67 let us glimpse the probable behavior of the Sun on timescales surpassing the duration of human civilization. M67 can serve as a solar proxy because its stars share a similar age and composition with the Sun. Previous surveys of M67 observed that 15% of its Sun-like stars exhibited chromospheric activity levels below solar minimum, which suggest that these stars might be in activity-minimum states analogous to the Maunder Minimum. The activity diagnostic used, the HK index (relative intensities of the Ca ii H & K lines integrated over 1 Å bandpasses), was measured from low-resolution spectra (R ≈ 5000), as is traditional and suitable for nearby, bright stars. However, for stars beyond the Local Bubble, the interstellar medium (ISM) imprints absorption lines in spectra at Ca ii H & K, which negatively bias activity measurements when these lines fall within the HK index bandpass. I model the ISM clouds in the M67 foreground with high-resolution spectra of blue stragglers and solar analogs. I demonstrate that ISM absorption varies across the cluster and must be accounted for on a star-by-star basis. I then apply the ISM model to a solar spectrum and broaden it to the lower spectral resolution employed by prior surveys. Comparing HK indices measured from ISM-free and ISM-contaminated spectra, I find that all stars observed below solar minimum can be explained by this ISM bias. I conclude that there is no compelling evidence for Maunder Minimum candidates in M67 at this time.

Neutral ISM, Lyα, and Lyman-continuum in the Nearby Starburst Haro 11*

Abstract Star-forming galaxies are believed to be a major source of Lyman continuum (LyC) radiation responsible for reionizing the early universe. Direct observations of escaping ionizing radiation have however been sparse and with low escape fractions. In the local universe, only 10 emitters have been observed, with typical escape fractions of a few percent. The mechanisms regulating this escape need to be strongly evolving with redshift in order to account for the epoch of reionization. Gas content and star formation feedback are among the main suspects, known to both regulate neutral gas coverage and evolve with cosmic time. In this paper, we reanalyze Hubble Space Telescope (HST)-Cosmic Origins Spectrograph (COS) spectrocopy of the first detected local LyC leaker, Haro 11. We examine the connection between LyC leakage and Lyα line shape, and feedback-influenced neutral interstellar medium (ISM) properties like kinematics and gas distribution. We discuss the two extremes of an optically thin, density bounded ISM and a riddled, optically thick, ionization bounded ISM, and how Haro 11 fits into theoretical predictions. We find that the most likely ISM model is a clumpy neutral medium embedded in a highly ionized medium with a combined covering fraction of unity and a residual neutral gas column density in the ionized medium high enough to be optically thick to Lyα, but low enough to be at least partly transparent to LyC and undetected in Si ii. This suggests that star formation feedback and galaxy-scale interaction events play a major role in opening passageways for ionizing radiation through the neutral medium.

Thermal Pressure in the Cold Neutral Medium of Nearby Galaxies

Abstract Dynamic and thermal processes regulate the structure of the multiphase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ( ) of the diffuse, neutral gas. By combining [C ii] 158 μm, H i, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of ∼1 kiloparsec. We find a distribution of thermal pressures in the K cm−3 range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of P th/k ≈ 3600 K cm−3. These results are consistent with thermal pressure measurements using other observational methods. We find that increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral media could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about ∼30% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

PRESSURE EQUILIBRIUM BETWEEN THE LOCAL INTERSTELLAR CLOUDS AND THE LOCAL HOT BUBBLE

Three recent results related to the heliosphere and the local interstellar medium (ISM) have provided an improved insight into the distribution and conditions of material in the solar neighborhood. These are the measurement of the magnetic field outside of the heliosphere by Voyager 1, the improved mapping of the three-dimensional structure of neutral material surrounding the Local Cavity using extensive ISM absorption line and reddening data, and a sounding rocket flight which observed the heliospheric helium focusing cone in X-rays and provided a robust estimate of the contribution of solar wind charge exchange emission to the ROSAT All-Sky Survey 1/4?keV band data. Combining these disparate results, we show that the thermal pressure of the plasma in the Local Hot Bubble (LHB) is P/k = 10, 700?cm?3?K. If the LHB is relatively free of a global magnetic field, it can easily be in pressure (thermal plus magnetic field) equilibrium with the local interstellar clouds, eliminating a long-standing discrepancy in models of the local ISM.

Lyman-αemission properties of simulated galaxies: interstellar medium structure and inclination effects

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population.

The Interstellar Medium Surrounding the Sun

The Solar System is embedded in a flow of low-density, warm, and partially ionized interstellar material that has been sampled directly by in situ measurements of interstellar neutral gas and dust in the heliosphere. Absorption line data reveal that this interstellar gas is part of a larger cluster of local interstellar clouds, which is spatially and kinematically divided into additional small-scale structures indicating ongoing interactions. An origin for the clouds that is related to star formation in the Scorpius-Centaurus OB association is suggested by the dynamic characteristics of the flow. Variable depletions observed within the local interstellar medium (ISM) suggest an inhomogeneous Galactic environment, with shocks that destroy grains in some regions. Although photoionization models of the circumheliospheric ISM do an excellent job of reproducing the observed properties of the surrounding ISM, the unknown characteristics of the very low-density hot plasma filling the Local Bubble introduces uncertainty about the source of ionization and nature of cloud boundaries. Recent observations of small cold clouds provide new insight into the processes affecting the local ISM. A fuller understanding of the local ISM can provide insights into the past and future Galactic environment of the Sun, and deeper knowledge of the astrospheres of nearby stars.

ROBO: a model and a code for studying the interstellar medium

We present ROBO, a model and its companion code for the study of the interstellar medium (ISM). The aim is to provide an accurate description of the physical evolution of the ISM and to set the ground for an ancillary tool to be inserted in NBody-Tree-SPH (NB-TSPH) simulations of large scale structures in cosmological context or of the formation and evolution of individual galaxies. The ISM model consists of gas and dust. The gas chemical composition is regulated by a network of reactions that includes a large number of species (hydrogen and deuterium based molecules, helium, and metals). New reaction rates for the charge transfer in $\mathrm H^+$ and $\mathrm H_2$ collisions are presented. The dust contains the standard mixture of carbonaceous grains (graphite grains and PAHs) and silicates of which the model follows the formation and destruction by several processes. The model takes into account an accurate treatment of the cooling process, based on several physical mechanisms, and cooling functions recently reported in the literature. The model is applied to a wide range of the input parameters and the results for important quantities describing the physical state of the gas and dust are presented. The results are organized in a database suited to the artificial neural networks (ANNs). Once trained, the ANNs yield the same results obtained by ROBO, with great accuracy. We plan to develop ANNs suitably tailored for applications to NB-TSPH simulations of cosmological structures and/or galaxies.