Next Article in Journal
The Formation of Fullerenes in Planetary Nebulae
Previous Article in Journal
AGBs, Post-AGBs and the Shaping of Planetary Nebulae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Conference Report

Challenges and Techniques for Simulating Line Emission

by
Karen P. Olsen
1,*,
Andrea Pallottini
2,3,
Aida Wofford
4,
Marios Chatzikos
5,
Mitchell Revalski
6,
Francisco Guzmán
5,
Gergö Popping
7,
Enrique Vázquez-Semadeni
8,
Georgios E. Magdis
9,
Mark L. A. Richardson
10,
Michaela Hirschmann
11 and
William J. Gray
12
1
School of Earth and Space Exploration, Arizona State University, 781 South Terrace Road, Tempe, AZ 85287, USA
2
Centro Fermi, Museo Storico della Fisica e Centro Studi e Ricerche “Enrico Fermi”, Piazza del Viminale 1, 00184 Roma, Italy
3
Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
4
Instituto de Astronomía, Universidad Nacional Autónoma de México, Unidad Académica en Ensenada, Km 103 Carr. Tijuana-Ensenada, Ensenada 22860, Mexico
5
Department of Physics & Astronomy, University of Kentucky, Lexington, KY 40506, USA
6
Department of Physics and Astronomy, Georgia State University, 25 Park Place, Suite 605, Atlanta, GA 30303, USA
7
Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117 Heidelberg, Germany
8
Instituto de Radioastronomía y Astrofísica, UNAM, Morelia, Michoacán 58089, Mexico
9
Cosmic DAWN Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen, Denmark
10
Department of Astrophysics, American Museum of Natural History, 79th Street at Central Park West, New York, NY 10024, USA
11
Sorbonne Universités, UPMC-CNRS, UMR7095, Institut d’Astrophysique de Paris, F-75014 Paris, France
12
CLASP, College of Engineering, University of Michigan, Ann Arbor, MI 48109, USA
*
Author to whom correspondence should be addressed.
Galaxies 2018, 6(4), 100; https://doi.org/10.3390/galaxies6040100
Submission received: 3 July 2018 / Revised: 24 August 2018 / Accepted: 1 September 2018 / Published: 20 September 2018

Abstract

:
Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of “Walking the Line”, a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.

1. Introduction

Line emission from the Interstellar Medium (ISM) of galaxies carries information that is crucial in understanding galaxy evolution. Observing line emission across the electromagnetic spectrum allows us to characterize the mass, composition, and chemical state of the ISM, as well as to trace galaxy properties such as star formation rate (SFR), metallicity and dynamics. For example, the emission from major cooling lines, such as H α or [C ii], is sensitive to the physical conditions (densities, radiation field) and dynamics of the ISM. In addition, emission lines work on all physical scales, from galaxy dynamics and inflows to turbulent and collapse motions in star-forming clouds and cores. By systematically comparing spectral-line signatures of different physical models, one can correctly identify the physical processes occurring in these regions. Furthermore, the emission from ionized interstellar gas contains particularly valuable information about the nature of the ionizing radiation sources in a galaxy. In fact, prominent optical emission lines are routinely used to estimate whether ionization is dominated by young massive stars (tracing SFR), an AGN or evolved, post-asymptotic giant branch (post-AGB) stars. Three of the most widely used line-ratio diagnostic “BPT” diagrams1, relate the [OIII]/H β ratio to the [NII]/H α , [SII]/H α and [OI]/H α ratios. These diagrams have proven useful in identifying the nature of the ionizing radiation in large samples of galaxies in the local Universe [2,3]. Complementary to line emission are the observations of absorption lines of the circumgalactic medium (CGM), which can give key information on the history of the feedback, in terms of chemical, ionization, and thermodynamical state of the outflowing/inflowing gas, that regulates the star formation process. Gas kinematics, from both emission and absorption, give information about large scale gas flows. Thus galactic outflows, from active galactic nuclei (AGN) and starbursts, can be combined with CGM absorption line observations, to study the star formation history, AGN activity history, and feedback processes that regulate both the evolution of the galaxy and its environment.
Looking back on the past three decades, the approach to creating synthetic observations of line emission has gone from simplified analytical modeling to complex simulations, increasing the reliability of the results. Since many important cooling lines emerge from the photodissociation regions (PDRs) of the ISM, modeling line emission from the PDRs has been an active field of research since the basic 1D PDR models came into place in the 80’s (e.g., [4,5,6,7,8]). This early modeling work included the reprocessing of starlight in the UV to infrared continuum by dust grains and polycyclic aromatic hydrocarbons (PAHs), and the stratified layers of increasingly photo-dissociated species as one moves through the neutral gas of a cloud towards the HII region. This picture is illustrated in the top panel of Figure 1. In the late 90s, the models presented in the 80’s [4,9] were used as basis for a modeling effort to create line ratio diagnostic plots of the [O i], [C ii], C i and CO FIR emission lines [10]. Later the modeling was improved with the online PDR Toolbox2 as a result [11,12].
Still, one of the main problems when comparing these models to actual observations is that the models are based on patches of gas with constant density, metallicity, and radiation field strength whereas any given galaxy (or region in a galaxy) will be a superposition of different gas states and radiation conditions as illustrated in Figure 2 (Ideally, each cloud should also have a radial gradient, but this feature has been omitted from the Figure for the sake of clarity). With improved simulations of galaxy (and cloud) formation, it became possible in the 21st century to use galaxy/cloud simulations as direct input to photoionization codes, that also calculate the line emission (e.g., [15,16,17,18,19,20,21,22,23]). This new approach gave way for a more realistic picture in which each region of the ISM is approximated by a combination of different gas conditions. Some of these simulations include radiative transfer and the non-equilibrium impact on the gas temperature of the local radiation field generated by nearby star formation (e.g., [24]), allowing for more accurate estimation of the nebular line emission. This method is being extended to simulations of the CGM, with increased resolution to capture the impact of outflows on accreting structure (see Section 4.4).
Improved numerical techniques such as those mentioned above have also meant a slight division in the field of line modeling, between groups would who specialize in creating the photoionization codes and groups that apply these codes to hydrodynamical simulations. The photoiononization codes have evolved to also calculate line emission and are kept updated as experimental values for collision strengths, chemical rate coefficients and photoelectric heating rates are being revised and databases of atomic and molecular line data increase in size (e.g., [25,26,27,28,29,30,31,32]).
Similar to the case of PDRs, significant progress has been achieved over the last few years in modelling nebular emission from ionized regions around young star clusters, AGN and post-AGN stellar populations in a full cosmological context. Yet, fully self-consistent models of this kind are currently limited by the performance of cosmological radiation-hydrodynamic simulations and insufficient spatial resolution on the scales of individual ionized regions around stars and active nuclei. To circumvent these limitations, some pioneer studies proposed the post-processing of cosmological hydrodynamic simulations and semi-analytic models with photoionization models to compute the cosmic evolution of nebular emission [33,34,35]. Some further improvement has been achieved recently by [36], not only accounting for the integrated nebular emission from young stars (as in [34,35]), but also from AGN and post-AGB stars, based on the star formation, AGN luminosity, gas density, and chemical enrichment histories of the simulated galaxies.
Across different aspects of simulating line emission, key questions remain to be answered:
  • What are the best emission lines to trace various ISM properties and ionizing sources in galaxies?
  • How can we use emission lines to trace feedback and ISM evolution with redshift?
  • How should absorption features be correctly interpreted?
  • Where do we stand in deriving sub-grid physics and comparing codes?
  • How do we coordinate our efforts?
To address such questions, the conference “Walking the Line” was held in Phoenix, Arizona on 14–16 March 2018. The main topic was Simulating Line Emission from Galaxies3 and the program4 included 35 min presentations from three invited speakers and 15 min presentations from 24 of the remaining 27 participants. In total, nine participants were PhD students and 1/3 of all participants were women. Most talks are now publicly available with video recordings online5. In addition, the conference had daily discussion sessions stimulating in-depth exploration of new topics and sparking new connections. This way, the meeting was an opportunity to discuss the challenges that need to be solved in order to make more realistic simulations and a fairer comparison with observations.
To take a specific case study where line emission has come to play an important role, one can consider the problem of measuring gas mass, which became a common thread for a large part of the workshop. Measuring the gas mass of galaxies, the fuel for star formation, as a function of cosmic time is a crucial component in understanding the formation and evolution of galaxies. Different groups have observed the gas content of main-sequence galaxies through cosmic time (the bulk of the galaxy population that is responsible for forming new stars at any epoch in cosmic history) either through their 12 CO (hereafter CO) emission or the dust continuum (e.g., [37,38,39,40,41,42,43,44,45,46,47]). Both approaches (CO and dust) rely on uncertain conversion factors between the observed luminosity and an estimated gas mass. ISM physical processes can further complicate the use of CO and dust continuum to reliably measure the gas mass of galaxies. For example, the contrast of the CO emission line and dust continuum against the cosmic microwave background becomes lower at higher redshifts (e.g., [23,48]) and under the influence of cosmic rays the CO molecule can be destroyed [49,50]. Reliable alternative measures of the gas mass in a well defined sample of galaxies are therefore essential to overcome the systematic uncertainties in the CO and dust conversion factors and to overcome the CMB contrast. Alternative options seen in the literature include for instance the emission from [CI] [51,52], [CII] [53], from H 2 O (e.g., [54]) and PAHs features (e.g., Cortzen et al. submitted 2018), and from optically thin isotopologues (e.g., [55]).The synergy between galaxy formation simulations, ISM chemistry simulations, and radiative transfer codes has the power to pave the way towards reliably measuring gas masses. This synergy allows for a controlled setting in which the conversion between sub-mm line/continuum strength and H 2 mass can be explored. This is important to (1) quantify under which physical conditions classical approaches to measure gas masses such as CO and dust continuum emission break down and (2) explore the robustness of other mid-IR and sub-mm emission lines such as PAHs, [CI], [CII] and H 2 O (although one could think of other examples) as a tracer of molecular hydrogen mass.
In this paper we summarize the topics covered in the conference, with particular attention to the break-out sessions of this workshop. We report our conclusions on the state of the art of modelling emission lines, by going to progressively larger physical scales (Section 2, Section 3 and Section 4), concluding by discussing the possible ways to move forward as a community (Section 5).
Note that the number of emission (and absorption) lines observed from the ISM of galaxies is ever increasing, and we did not attempt to cover them all with this workshop. However, to give the reader a quick overview of the lines that were discussed at the workshop and are commonly used to diagnose the ISM and/or CGM, we present in Table 1 a list of lines that will be mentioned in this paper.

2. On the Micro-Physical Level

2.1. Tools for Solving Level Populations and Line Excitation

Several software tools are available online for the calculation of line emission. While not covering all, the tools presented and discussed at this workshop are listed alphabetically in Table 2 for quick reference and comparison. Before comparing the tools, it is worth noting that they are built to achieve different levels of precision at the cost of computation time. As such, codes like LIME, MOLLIE, and RADMC-3D are very accurate in 3D, but slow compared to Cloudy and DESPOTIC which solve for chemistry and temperature in 1D and are intended for larger parameter studies that can be analyzed fast.
In terms of chemical networks, Cloudy is the most comprehensive tool with 625 species and a large range of permitted densities that allows for application to models on scales from clouds and HII regions to galaxies (e.g., [58,59,60]). DESPOTIC has been applied in post-processing to galaxy simulations mainly (e.g., [21,22,61]) but also to smaller regions such as the Galactic central molecular zone [62]. As mentioned above, Cloudy and DESPOTIC are restricted to problems in 1D whereas the remaining tools in Table 2 work in 3D. LIME, MOLLIE, and RADMC-3D do exact radiative transfer in 3D and are typically used on smaller scales where non-symmetric features such as filaments and turbulence become important (e.g., [63,64]). ART2 also does radiative transfer in 3D including the continuum from far-UV to radio wavelengths and incorporating the resonant line Ly α (e.g., [29,65,66]) (an updated version of the code including CO and some prominent far-infrared fine-structure lines such as CII, OI, OIII and NII will be out later this year). Finally, MAIHEM is unique in taking into account the non-equilibrium effects of turbulence on line emission (c.f. Section 3.3).

2.2. Correcting Heavy-Element Energy Levels

Predictions on line emission modeling depends on the reliability of fundamental atomic data. Software tools such as those discussed in the previous section use theoretical and experimental atomic and molecular rates and energies to calculate emissivities and line wavelengths respectively, mostly from external databases as CHIANTI [69,70] or NIST. Nowadays, theoretical energy levels of heavy many-electron ions disagree with experimental values, leading to incorrectly predicted wavelengths. This happens in all the energy and wavelength ranges. For instance, there are differences between the energies of Fe II reported by the NIST and CHIANTI databases, where several terms reported by NIST are missing in CHIANTI. Similarly, efforts to fit the spectra of X-ray satellite observations have accentuated the necessity of including transitions between levels with high main quantum number n in the theoretical calculations (e.g., [71,72,73,74]). A displacement of about 6 Å has been found in some X-ray lines due to theoretical and experimental differences for energy levels [74], and observation-theory intensity ratios as large as 10 have been reported for several lines in Fe XVI. It is believed that some of the differences could be explained by the blending of the lines with emissions from other ions [75]. Moreover, experimental measurements using ion traps (EBIT) have revealed significant differences with theoretical energies (see e.g., [74,76,77]). This situation has improved over the last decade. The high computational power reachable nowadays allows calculating energies and collision strengths for groups of ions with the same iso-sequence (e.g., [78,79,80,81]), or the same type of transitions for all ions [82]. As different codes calculate new sets of data as well as fit spectra (e.g., FAC [83], GRASP [84], AUTOSTRUCTURE [85]) and new data are stored in databases such as CHIANTI [69,70]) and HULLAC [73], there is better general agreement between theory and experiments (see for example [86]). However, there are still lines which remain unidentified, while others are poorly described because they are blended or have wrongly assigned transitions [87].
To illustrate this, two possible cases of disagreement between theory and experiment, both taken from CHIANTI (v7, [88]), are plotted in Figure 3. In the upper panel, the theoretical energies for Fe 16 + are, at most, roughly 1% off the experimental values. However, theory predicts many more levels that have not yet been observed. In the lower panel, the few levels predicted for Mg + are systematically lower than experimental values. Although the error for Mg + levels is around 3%, in the worst case, a line to the ground level can be shifted by as much as 25 Å.

2.3. Radiative Transfer of UV to Infrared

The reprocessing of UV light into infrared via dust grains is crucial for the energy balance in the ISM and thereby the line emission. In addition, it allows for a determination of the dust continuum emission which is as important an observable as line emission. However, radiative transport (RT) is a highly complex, non-local, and non-linear problem. It is considered one of the four grand challenges in computational astrophysics.6 Considerable progress has been achieved over the last 20 years in terms of carrying out exact 3D radiative transfer simulations of dust, as outlined in [92]. Yet, modeling dust radiative transfer in large hydrodynamical simulations is very computationally expensive, and it is usually ignored. A few pioneering studies that account “on the fly” for radiative transfer and the effects of dust are (e.g., [24,93,94,95,96,97,98]). RAMSES-RT [99,100] implements radiative transfer in a cosmological code, albeit without a full treatment for dust emission. Examples of other codes that transfer radiation on the fly include EMMA [101] in a cosmological setting, as well as ORION [102] and Athena [103] in an ISM setting. Generally speaking, the impact of dust on the energy balance of simulated galaxies is typically not accounted for when computed on the fly. Alternative approaches to modeling dust are therefore worth exploring.
Ideally, a full simulation of the hydrodynamical evolution of astrophysical clouds, including galaxies, needs to account for the propagation of radiation through the simulated volume, and for its interaction with matter. The latter requires detailed knowledge of the quantum structure of all ions, molecules, and dust present in the simulated gas, as well as keeping track of the chemical reactions that couple the various species. For a simulation with a few million cells or particles, the problem becomes extremely computationally expensive, and the non-local couplings induced by radiation render it practically intractable with the computational means available at present.
An alternate solution is to solve for the hydrodynamics and the microphysics separately, in post-processing. In this approach, a code such as those listed in Table 2 is required to simulate the microphysics. Cloudy [13] performs 1D spectral synthesis simulations of astrophysical clouds, accounting for all ions of elements up to Zn, a large number of molecular species, and all important microphysical processes. The propagation of radiation is done with the escape probability formalism, which is an expedient, and generally accurate approximation [104]. However, a few pathological cases exist for which escape probabilities are known to not give correct results [105]. Exact radiative transfer methods are required if a code is to be used in conjunction with hydrodynamical codes. The Cloudy team will seek to implement exact radiative transfer, eventually with full treatment of the interactions of dust and radiation. Even with these improvements, however, due to their 1D nature these simulations will not capture the time-varying three-dimensional radiative fields that exist in nature. An exact 3D treatment of radiative transfer will be needed to connect the microphysics and 3D hydrodynamical models.

3. Cloud-Scale Simulations

3.1. The Internal Density and Velocity Profile of Molecular Clouds

In calculating emission line strengths from molecular clouds, the adopted radial profiles of densities and radial velocity profile of the gas can drastically change the result. For example, at the scale of molecular clouds and their substructure, the so-called “infall profiles” are routinely used in identifying collapsing dense cores (see, e.g., [106]). These profiles consist of self-absorbed lines with a blue excess in moderately optically thick lines, in conjunction with single-peaked profiles of similar total width in optically thin lines. These profiles have been interpreted in terms of a dense molecular cloud core undergoing “inside-out collapse” [107], where the inner regions of the core are undergoing collapse, while the external envelope is in hydrostatic equilibrium, the two regions mediated an expanding rarefaction front. The assumed underlying radial velocity profile is essential for the generation of the lines. Recent studies (e.g., [108] Loughnane et al. 2018, submitted) have focused on the effect different assumed radial density and velocity profiles have on the emitted lines, and favored the collapse of marginally-unstable structures, characterized by “outside-in” velocity profiles rather than the canonical inside-out one. Thus, traditional line modeling has to assume a velocity profile, but there may be a degeneracy such that different combinations of density and velocity profiles may produce similar spectral lines. Synthetic observations of numerical simulations are thus essential in order to break the degeneracy, since they produce dynamically self-consistent density and velocity profiles, restricting the range of physically acceptable choices and thus facilitating the identification of the actual physical conditions that produce the signature spectral lines.

3.2. Simulating the Ionizing UV Field that Clouds Are Embedded in

The standard way to model the ionizing UV radiation, which ultimately determines the chemistry of the ISM, is to use stellar population synthesis (SPS) codes. In the SPS methodology, the integrated spectrum from populations of stars with given characteristics (IMF, star formation history, metallicity, etc.) is derived from stellar evolutionary tracks, which give the H-R diagram positions of stars of given masses with time, and from individual stellar spectra, which can be empirical or theoretical. If theoretical, they are computed with stellar atmosphere codes. Figure 4 shows examples of model predictions for populations of stars resulting from an instantaneous burst of star formation, which were computed with starburst997.
However, using SPS codes comes with a few cautions worth knowing. For example, SPS codes from independent research groups account for different stellar astrophysics (rotation, close-binarity, etc., see e.g., [109]) and predict different ionizing spectra (e.g., Figure 2 of [110]). In addition, the ionizing continuum from individual ionizing stars cannot be observed directly, because it is strongly absorbed in the ISM. This makes constraining model predictions of the ionizing spectrum of populations of stars difficult. The escape fractions of ionizing photons from a few distant galaxies have been estimated from direct detections of ionizing photons (e.g., [111]). However, these are insufficient to constrain the ionizing spectrum from massive star populations, because they are just escape fractions. Thus, we must rely on theory. However, there remain large uncertainties in massive star evolution (trajectories in H-R diagram and stellar atmospheres), particularly for post-MS phases and stars more massive than 20 M .

3.3. Direct Observations of Low-Metallicity Massive Stars

Massive ( M > 8 M ) extremely metal-poor stars ( Z / Z 0 . 1 ) are expected to be present in extremely metal-poor galaxies, i.e., galaxies with ionized-gas oxygen abundances of ≤solar/10. Since such galaxies are farther than ∼10 Mpc away, they are too distant to enable observations of their individual massive stars. This is a problem for understanding a variety of astrophysical objects which range from the first stars and galaxies which re-ionized the universe to the metal-poor stars that give rise to the formation of heavy binary black holes. Direct observations of individual, massive, extremely metal-poor stars will only be possible if telescopes such as LUVOIR8 and HabEx9, which are being proposed to NASA, ever see the light.

3.4. Implementing Turbulence and Shocks in Simulations of the ISM

In addition to the importance of modeling stars as described in the previous two subsections, another important source of heating and ionization is turbulence and shocks. The introduction of density gradients and additional sources of ionization through these processes can have a significant impact on the observed spectrum, and hence the emission line diagnostics. Examples of how different levels of turbulence affect the density structure of ISM gas is shown in Figure 5.
One solution to understanding the effects of turbulence is through the publicly available MAIHEM10 code. This FLASH (version 4.4) based code tracks ionizations, recombinations, and species by species radiative cooling for hydrogen, helium, carbon, nitrogen, oxygen, neon, sulfur, calcium, and iron [32,67,68]. By including a solenoidal stirring mechanism, turbulence is introduced into the gas that then propagates through all spatial scales.
At the meeting, there was an interest in applying MAIHEM to Baldwin-Phillips-Terlevich (BPT) diagrams [1,115] and other UV diagnostics, as shown in Figure 6, and other UV diagnostics. Specifically, turbulence in the gas might help explain the observed dispersion in some spatially-resolved BPT diagrams of nearby AGN [116], including a larger scatter seen in the neutral optical [O I] emission line as compared to [N II] and [S II].
In addition, the treatment of shocks is particularly complex as it requires understanding the physical conditions in the emission line gas before, during, and after the shock event. These events can cause collisional ionization and the resulting recombinations may release additional ionizing photons that interact with surrounding material, depending on the shock velocity. A new database of shock models calculated using the MAPPINGS V11 code [117] was presented at the workshop, and details can be found by consulting the documentation available on the Mexican Million Models database (3MdB) website12.

4. Galaxy-Scale Simulations

4.1. From Cloud to Galaxy Scale Simulations

Cloud scale simulations starting at masses of 10 4 10 6 M have shown us how gas structures emerge and evolve for different inflow rates from converging streams of diffuse gas, corresponding to large-scale velocity dispersions in galaxy simulations (e.g., [118,119,120,121,122,123,124]). However, these simulations are mostly made at CNM initial conditions, i.e., in a limited temperature range. Projects dealing with line emission on galaxy scales would benefit greatly from seeing the impact on small-scale simulations from different temperatures reflecting different radiation backgrounds and metallicities. The community could also benefit from more ISM simulation work on scales in between cloud and galaxy sizes (e.g., [125]) and also observations of individual molecular clouds in external galaxies (e.g., PHANGS survey13 and the LMT effort14).
When simulating line emission from objects extracted from analytical/cosmological simulations, it is often a challenge to identify the minimum amount of information you need in order to get a decent estimate of the emission. For example, some simulations calculate abundances of specific elements on the fly, whereas others only contain one number for the metallicity and an abundance table must be adopted. Including several elements separately instead of a fixed abundance table, has been shown to make a significant difference for modeling of the fine-structure [C ii] line [126].
Additionally, other possible uncertainties are present when considering galaxy simulations that include in their modeling both chemical networks (e.g., [127,128]) and radiative transfer (e.g., [129,130]); while such simulations can capture dynamical and thermodynamical effects due to photochemistry, they are not refined enough to calculate line emission, which is to be done in post-processing. However, in post-processing line calculations typically consider photoionization equilibrium and it is difficult to account for the past thermal and ionizational history of the simulated gas.
One of the more valuable conclusions from the discussions on galaxy-scales simulations, was the importance of simulating more than one emission line simultaneously. By simulating different lines, arising in different ISM phases, and comparing with observations, one ensures that the post-process recipes not only satisfy what is seen in one ISM phase, but is consistent across the entire galaxy.

4.2. Mapping Simulated Galaxy Samples to Observed Samples

An important step in simulating line emission from galaxies and from interstellar clouds is to compare with observed analogues of the model systems. Here, a number of issues continue to stand out, of which we discussed the following:
1
Observed SFRs come from tracers such as H α , UV, IR and radio sensitive to a time averaged SFR of 10 to a few hundred Myr. Model SFRs are often the instantaneous SFRs. The question here is whether models should use time averaged SFRs or produce continuum and line emission to measure SFRs like observers do? One of the better solutions is to do both and identify any possible biases that are introduced by using observational methods.
2
Model metallicities are often M metal / M gas | star , whereas observed metallicities are determined using a number of ionized emission lines and presented as e.g., 12 + log ( O / H ) . As with SFR, the best option at the moment might be to directly calculate 12 + log ( O / H ) if the oxygen abundance is tracked by the simulation used, although even this approach has problems as the observed 12 + log ( O / H ) will ultimately be luminosity-weighted which is hard to replicate in a simulation.
3
A realistic comparison between modeled and observed UV/optical emission lines (and continuum) requires the correct treatment of dust absorption within simulations. The amount of intervening dust in galaxy scale simulations between the photon source and the observer is ideally calculated self-consistently using dust chemical networks (e.g., [131,132,133,134]). However, often the amount of intervening dust is scaled as a function of the gas-phase metallicity by assuming a fixed grain size distribution. Finally, the geometry of the galaxy and the relative star-to-dust location is critical in correctly determining the observed properties (e.g., [135,136]), including sub-grid modeling of the dust properties of the birth-clouds of young stars.
4
Atomic hydrogen emission at 21cm is out of reach of current instrumentation beyond at z > 0 . 3 (SKA and SKA pathfinders such as APERTIF, MEERKAT, and ASKAP will be a remedy for this). This hampers the validation of the models that predict an evolution of atomic to molecular gas fraction with molecular gas dominating the gas budget within the optical scale of the galaxies at higher redshifts.
5
Analyzing a large sample of strongly star-forming SDSS galaxies with nebular He ii λ 4686 emission has shown that the luminosity of this line can only be reproduced with single bursts of star formation of 20 percent solar metallicity or higher, and ages of 4–5 Myr, when the extreme UV continuum is dominated by Wolf-Rayet stars [137]. He ii λ 1640 in the UV is 10 times stronger than the optical He ii λ 4686 line, and has been observed in broad (FWHM 1000 km s 1 ) and narrow emission for tens of dwarf star-forming galaxies also selected from SDSS. Reproducing the luminosity of the narrow nebular He ii λ 1640 emission from these galaxies has been challenging, even with state-of-the-art spectral synthesis models which combine the newest Charlot & Bruzual population synthesis models, which include very massive (300 M ) stars, with photoionization models, as described in [138]. This failure is reported for instance in [139]. In this workshop, Aida Wofford presented the case of one of the most metal poor nearby galaxies known, SBS 0335-052E, which has a metallicity of solar/20. None of the models which they tried were able to reproduce the luminosity of the He ii λ 1640 line. This is a problem because this line will be one of the only diagnostics of massive stars in future observations with large telescopes such as JWST (see a recent technical and scientific description in [140]), TMT [141], and e-ELT15. These telescopes will obtain rest-frame UV spectra for thousands of galaxies, at redshifts between 10 and 15, in the era of re-ionization when the first stars and galaxies formed.
6
How do AGN affect the comparison between observed and simulated galaxies? Comparing galaxies of similar SFR can become problematic, as the ionizing radiation and presence of a radio jet will enhance H α and radio emission that is typically used as a SFR diagnostic. In addition, the radiation will heat dust that is often used in mass estimations. See the following section for more on the effect of AGN and mass outflows. Although not directly discussed at this workshop, X-ray dominated regions (XDRs) are another important result of having an AGN present, and they must be modeled in order to reproduce certain emission lines. For example, this is extremely relevant for high-J CO lines, and it is in fact still not clear if high-J CO lines are influenced more by the presence of X-ray Photons or shocks [142]. Important theoretical work on modeling XDRs was published in 2005 [143] and Cloudy has also been used to model XDRs (e.g., [144]).

4.3. The Effects of AGN and Mass Outflows on Line Emission

The presence of an AGN can profoundly affect every aspect of observed emission line profiles including their velocities, widths, symmetry, and relative flux ratios. The AGN emits radiation anisotropically, and under the proper physical conditions can effectively couple with the ISM. This gives rise to outflows [145] of molecular and ionized gas that are likely radiatively driven. Understanding these outflows is critical as they may deliver feedback to the galaxy by clearing the inner regions of star forming gas, and assist in establishing the observed scaling relationships between galaxies and their supermassive black holes (SMBHs) [146,147,148,149].
At the workshop a subset of topics were covered, highlighting the physical mechanisms that couple the ISM to the AGN radiation field and affect the observed emission line properties. As an example, recent results from post-processing high-resolution, cosmological zoom-in simulations of massive galaxies with nebular emission models of galaxies and AGN were presented at the workshop [36,150]. One of the main results is that at least for massive galaxies, AGN-driven outflows seem to be a necessary component in order to reproduce the observed cosmic evolution of [OIII]/H β at fixed [NII]/H α via strongly regulating the star formation history, which controls nebular emission from young stars via the ionization parameter (see Figure 7).
There was also a particular focus on what can be gleaned from spatially resolved observations and models. In the local universe, spatially resolved BPT diagrams have proven their value in understanding how the physical conditions of the emitting gas change with distance from the SMBH, and can be used to constrain Cloudy photoionization models of ionized mass outflows [116]. A first theoretical modeling of such spatially resolved BPT diagrams, using cosmological zoom-in simulations of massive galaxies shows that central low ionization emission in massive galaxies is not necessarily caused by an AGN, but can be also due to the ionizing radiation coming from post-AGB stars (Hirschmann et al., in prep.). For some galaxies, direct signatures for AGN-driven outflows are also visible in nebular emission line maps. These studies have primarily focused on the optical nebular emission line gas; however, observations reveal that powerful outflows of molecular gas may equal or even dominate the feedback in some galaxies. A conundrum in this area has been understanding how molecular outflows remain stable over time, as the physical conditions in the gas suggest many molecules should be destroyed. New modeling with advanced chemistry networks has found that molecules can form in-situ within the outflow, providing a significant opacity source that helps to drive the molecular outflows to the observed velocities [151]. Additional work in this area is needed to fully understand how the presence of an AGN modifies the standard emission line diagnostics for determining properties such as dust temperature, metallicity, and the inferred star formation rates.
A special case combining the effect of AGN, going from the host galaxy to the surrounding CGM, is the high-redshift radio galaxy MRC 0943-242. S. Kolwa presented nuclear emission and CGM absorption observations of MRC 0943-242 (see Figure 8), whose emission is spatially integrated over a 0.6 arcsec (4.75 kpc) radius aperture centered at what is assumed to be an AGN. Hence, the lines reveal gas in the ISM in close proximity to the nucleus of the galaxy and ionized by the AGN radiation. The emitting gas is turbulent with a velocity width larger than 1000 km s 1 . Models and simulations are important in understanding the physical mechanisms behind the turbulence in the ionized gas caused by processes surrounding the AGN. Further, Figure 9 shows absorption from the gas of various column densities in the CGM. The kinematics of this gas will be studied soon with IFU observations. What impact the AGN has on heating or removing the galaxy gas may be uncovered by combining observations of the ISM emission excited by the AGN with emission from CGM absorbers possibly resulting from AGN outflows. Further, these observations may discriminate between different models of the origin and fate of CGM gas.

4.4. Simulating Line Absorption from the CGM

While the workshop focus was not specifically on absorption line profiles, or the CGM, several talks and break-out discussions emphasized the connection between modeling emission lines and absorption lines, and the value in collaborations between these two fields. In particular, the CGM can give further constraints on the nature of the galaxy that could be important for converting emission line observations to other characteristics of the galaxy, such as gas mass.
Simulations of the CGM and its corresponding absorption line profiles are progressing, even given the difficulty of resolving this medium. Observations of distant quasars with absorption from the CGM of foreground galaxies provide useful metrics to help improve simulations. The characteristics of these observations are described well in [152] and references therein. In particular, observations find an abundance of OVI in the CGM of galaxies of a range of masses, only well matched in simulations at high galaxy masses above 10 10 . 5 M [153]. Other simulations post-process their outputs with Trident [154] to produce mock absorption line profiles and compare them with quasar absorption line observations. The predicted CGM column densities and covering fractions of ions such as OVI depend on the models of star formation and stellar feedback (see Figure 10; Chuniaud et al. in prep). Star formation models influence the amount of metals produced in the simulations, while feedback models impact their ability to leave the galaxy and persist in the CGM.
Observations also find that these OVI lines typically have suprathermal line widths [159]. However, this could be due to blending of multiple components along the line of sight. Simulations could reveal the nature of these broad components, but this may depend on resolution. The dependence of the simulated CGM on the numerical resolution is a debated topic. It is unclear what are the scales of the cool gas in the CGM, which could range from tens of pc down to sub-pc scale (e.g., [160]). Increasing the resolution in the CGM leads to better resolving the central density peaks of outflowing and inflowing material in the CGM. The initially slightly higher density has faster cooling rates, which is amplified over time. Thus, simulations with better resolution in the CGM, such as in NEPHTHYS presented by M. L. A. Richardson, form more gas clumps, mostly in inflowing gas. Inflows are enriched from outflows, leading to enhanced cooling. However even the high resolution CGM simulations shown by Richardson et al. (nephthys) which have 100 physical pc resolution, may by numerically supported. Observed clumps in their simulations are roughly 10–20 cells across.
Many models of the CGM assume the ionizing radiation source to be the background UV reionization field, but the central galaxy can dominate the background over periods of time, including during AGN flaring or starburst episodes (e.g., [153,161]). Non-equilibrium treatment of the ionization state could make this effect be longer-lasting. Future work, including the production of look-up tables to convert simulation data to predicted ionization states, will need to consider this.
Finally, the CGM will be better constrained in the future with the ELT and other next-generation observatories that will be able to use Lyman-break galaxies, far more numerous than quasars, to have multiple line-of-sight absorption profiles per foreground galaxy. This will allow for spatially resolved study of the CGM, and give more instruction for combining different foreground galaxy column density to make general CGM column density profiles.

5. Discussion: How Can We as a Community Move Forward?

This workshop was the first of its kind (to our knowledge) to attempt to create a community of astronomers simulating line emission and working on the comparison with observations. In order to continue the sharing of knowledge within the community, we discussed ways in which to make the material of the workshop available to a larger group and encourage collaborations in the future. For this workshop, it was decided to use zenodo.org where all speakers had the option to upload their slides together with a recording of their talk (where available) to a community space16. This approach automatically provides a DOI for each upload for easy reference and we hope that future workshops of this kind will likewise attempt to make the presentations publicly available.
It should also be mentioned that the specific problems of simulating line emission from plasma are not only faced by the astronomical community. Within the field of plasma physics, collisional-radiative models are used to analyze plasma kinetics and spectra as well as to design plasma experiments in the laboratory. Yet, the two fields, i.e., that of astronomy and that of plasma physics, are currently far from fully benefiting from each other. We therefore wish to highlight corresponding established workshop series for non-LTE modeling, namely the “NLTE Code Comparison” workshops17 and the “Spectral Line Shapes in Plasmas” workshops18.

5.1. Getting Involved With the Community of Developers

As Table 2 shows, there are several software tools available for simulating line emission, each with their own group of developers behind. However, in order to achieve the best usage of these tools and adapting them to astrophysically relevant scenarios, it is crucial that users are able to connect with the developers to request support and updates, or even to adjust the underlying code themselves. Here we detail how this is done for two examples listed in Table 2.

5.1.1. The Cloudy Community

Small tutorials on using Cloudy have been uploaded to youtube by individuals not part of the Cloudy developer team and such video uploads are much appreciated by many and strongly encouraged by the Cloudy community. It was suggested to establish a Cloudy-branded youtube channel or youtube playlist for these videos. Also, many people are not aware that you can edit the databases used by Cloudy to for example include new species. The Cloudy quickstart guide has been noted to be an excellent tool for introducing new researchers and students to making predictions with the software. Maintaining it and keeping it up-to-date will be a priority. The Cloudy team is always looking for more participation on the Yahoo Groups forum. A goal is that questions and discussions will take place there among users, rather than solely between a user and cloudy developers. An interactive wiki has also been suggested as a place for users to post about specific cases or issues in their Cloudy research.

5.1.2. LIME

LIME resides on the code hosting platform GitHub (https://github.com/lime-rt/lime) from where it can be downloaded. A small team of active developers maintain and develop the package and minor releases appear regularly, approximately once per year. Users that experience problems when using LIME or have ideas to improve on the code, are welcome to create issues on GitHub where they will be answered by one or more developers. It is also possible to fork the code and develop it for specific needs. Documentation exists in the form of a written PDF manual as well as an issue tracker. Both are available on GitHub. A bi-annually summer school on Monte Carlo techniques in radiative transfer problems, aimed at PhD-students, is held at University of St Andrews School of Physics & Astronomy and this school offers both a lecture and tutorials on how to use LIME.

6. Conclusions

This paper reviews recent progress and identifies current problems in the field of simulating line emission from the ISM and CGM. Conclusions presented here are a result of discussion sessions at the workshop “Walking the Line” held in Phoenix Arizona on 14–16 March 2018. This international meeting was the first of its kind to gather astronomers specializing in the modeling of line emission and brought together 30 scientists from grad student to professor level. Most talks with follow-up questioning sessions have been made available online, but as a more digestible release of knowledge, we have condensed the outcome on the most relevant topics in this paper. Returning to the questions posed in the introduction, below are the main conclusions, focusing on the methods used to give the most reliable answers to each question.
What are the best emission lines to trace various ISM properties and ionized sources in galaxies? The best method to answer this question of course depends on the type and scale of the emitting source in the ISM considered. In order to provide the reader with an overview, though not complete, of tools used to calculate line emission either on-the-fly in a simulation or in post process, we listed the software tools that were discussed at the workshop. One outstanding issue that can affect these tools, mostly when considering line emission in the X-ray but also at longer wavelengths [162], is to correct theoretical heavy-element energy levels that do not agree with experimental values.
Key to deriving the line emission from any cloud is the UV field irradiating that cloud. However, our knowledge on the ionizing continuum of stars is severely hampered by the observational difficulties and thus currently restricted to theoretical predictions. Direct observations of individual, massive, low-metallicity stars will only be possible with telescopes such as LUVOIR and HabEx. Once the UV field is known, radiative transfer of UV into infrared light is a crucial part of any method used, especially when studying neutral or molecular regions. While still too computationally expensive to run alongside galaxy-scale simulations, work is ongoing to include it more regularly in post-process as part of any line emission simulation.
In order to compare simulated galaxies with observed ones, we need to put more thought into the choice of physical properties being compared. For example, a SFR that is modeled as instantaneous can hardly be compared to an observed SFR that is essentially derived as an average over 10 to a few-hundred Myr. Many other parameters such as the effects of an AGN or the inclination of a galaxy can distort an attempt to compare model galaxies with observed ones. Finally, it was concluded that the aim should always be to simulate more than one emission line simultaneously, in order to assure consistency across an the entire galaxy.
How can we use emission lines to trace feedback and ISM evolution with redshift? Feedback from AGN in the form of radiation and outflows has a profound effect on the nebular line emission. As an example, simulations of massive galaxies with and without AGN were presented in which the presence of an AGN could clearly be identified from the [OIII]/H β and [NII]/H α emission-line ratios. AGN-driven outflows thus strongly regulates the star formation history, at least for massive galaxies, which controls nebular emission from young stars via the ionization parameter. As an observational example, the high-redshift galaxy MRC 0943-242 is believed to contain an AGN, with line emission (and absorption) revealing ionized ISM in close proximity to the nucleus as well as outflows and high levels of turbulence.
In close relation to stellar and AGN feedback, turbulence and shocks are another important source of heating and ionization. Recently developed tools for the treatment of turbulence and shocks were presented to aid the line emission simulations in cases where necessitated.
On cloud-scale level, simulating the line emission from collapsing dense cores can reveal whether the molecular cloud is undergoing an “inside-out’’ or “outside-in’’ collapse. Work is undergoing to break this degeneracy by making synthetic observations of signature spectral lines for comparison with observations.
How should absorption features be correctly interpreted? Simulations of line absorption through the CGM are faced by its own set of problem, where for example the treatment of the ionizing field from the galaxy source can play a crucial role but has yet to be fully treated in a non-equilibrium manner. Another choice to be made when simulating the CGM is the level of spatial resolution needed, as better resolution leads generally to more gas clumps and a point of convergence must be found.
Where do we stand in deriving sub-grid physics and comparing codes? When studying ISM evolution with redshift, line emission must be derived from cosmological simulations which capture the evolution but comes at low spatial resolution compared to semi-analytical models. In addition, full chemistry and radiative transfer must typically be applied in post process. The uncertainties in the approach could be alleviated by more work on scales between galaxies and clouds, because simulations of the latter come at higher resolution and could hence be used as a benchmark on larger scales.
How do we coordinate our efforts? We hope to have set an example for more theoretical workshops of this kind to come, and look forward to seeing the results of the ongoing projects presented here. The developers of the software tools used to derive line emission, increasingly encourage their users to report problems or suggest improvements to their code via online platforms. With the list of caveats identified in this workshop, albeit long, we continue to believe that the future of simulating line emission is bright, especially when research groups in the field come together to find and solve common problems.

Author Contributions

The conference was led by K.O. who planned it with nine other members on the science organizing committee, including E.V.-S.; K.O. contributed to Section 1, Section 2.1, Section 5 and Section 6; A.P. contributed to the Section 1, Section 4.1 and Section 4.2; A.W. contributed to Section 3.2 and Section 3.3; M.C. contributed to Section 2.3 and Section 5; M.R. contributed to Section 3.4, Section 4.2 and Section 4.3; F.G. contributed to Section 2.2 and Section 5; G.P. contributed to the introduction and Section 4.2; E.V.-S. contributed to Section 3.1; G.E.M. contributed to the introduction; M.L.A.R. contributed to Section 4; M.H. contributed to Section 4.3, and W.G. contributed to Section 4.4.

Funding

The workshop of which the results are presented was in part sponsored by professor Rogier Windhorst and the School of Earth and Space Exploration at Arizona State University. M.R. gratefully acknowledges support from the National Science Foundation through the Graduate Research Fellowship Program (DGE-1550139). M.C. was supported by NASA through a grant (HST-AR-14556.001-A) from the Space Telescope Science Institute. F.G. acknowledges support by NSF (1412155). GEM acknowledges support from a research grant (13160) from Villum Fonden.

Acknowledgments

The authors thank the two anonymous referees for their thoughtful comments and suggestions. The authors thank all participants at the Walking the Line 2018 conference for presenting their work and stimulating discussions. The authors also thank Mark Krumholz, Livia Vallini, Christian Brinch, Robert Loughnane, Peter van Hoof, Moupiya Maji, William Gray, and Sthabile Kolwa for providing additional info and comments to this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baldwin, J.A.; Phillips, M.M.; Terlevich, R. Classification parameters for the emission-line spectra of extragalactic objects. Publ. ASP 1981, 93, 5–19. [Google Scholar] [CrossRef]
  2. Kewley, L.J.; Dopita, M.A.; Sutherland, R.S.; Heisler, C.A.; Trevena, J. Theoretical Modeling of Starburst Galaxies. Astrophys. J. 2001, 556, 121–140. [Google Scholar] [CrossRef] [Green Version]
  3. Kauffmann, G.; Heckman, T.M.; Tremonti, C.; Brinchmann, J.; Charlot, S.; White, S.D.M.; Ridgway, S.E.; Brinkmann, J.; Fukugita, M.; Hall, P.B.; et al. The host galaxies of active galactic nuclei. Mon. Not. Roy. Astron. Soc. 2003, 346, 1055–1077. [Google Scholar] [CrossRef] [Green Version]
  4. Tielens, A.G.G.M.; Hollenbach, D. Photodissociation regions. I–Basic model. II–A model for the Orion photodissociation region. Astrophys. J. 1985, 291, 722–754. [Google Scholar] [CrossRef]
  5. van Dishoeck, E.F.; Black, J.H. Comprehensive models of diffuse interstellar clouds—Physical conditions and molecular abundances. Astrophys. J. Suppl. 1986, 62, 109–145. [Google Scholar] [CrossRef]
  6. van Dishoeck, E.F.; Black, J.H. The Photodissociation of Interstellar Co/; Molecular Clouds, Milky-Way and External Galaxies; Lecture Notes in Physics; Dickman, R.L., Snell, R.L., Young, J.S., Eds.; Springer: Berlin, Germany, 1988; Volume 315, p. 168. [Google Scholar]
  7. van Dishoeck, E.F.; Black, J.H. The photodissociation and chemistry of interstellar CO. Astrophys. J. 1988, 334, 771–802. [Google Scholar] [CrossRef]
  8. Sternberg, A.; Dalgarno, A. The infrared response of molecular hydrogen gas to ultraviolet radiation–High-density regions. Astrophys. J. 1989, 338, 197–233. [Google Scholar] [CrossRef]
  9. Hollenbach, D.J.; Takahashi, T.; Tielens, A.G.G.M. Low-density photodissociation regions. Astrophys. J. 1991, 377, 192–209. [Google Scholar] [CrossRef]
  10. Kaufman, M.J.; Wolfire, M.G.; Hollenbach, D.J.; Luhman, M.L. Far-Infrared and Submillimeter Emission from Galactic and Extragalactic Photodissociation Regions. Astrophys. J. 1999, 527, 795–813. [Google Scholar] [CrossRef]
  11. Pound, M.W.; Wolfire, M.G. The Photo Dissociation Region Toolbox. In Astronomical Data Analysis Software and Systems XVII; Astronomical Society of the Pacific Conference Series; Argyle, R.W., Bunclark, P.S., Lewis, J.R., Eds.; Astronomical Society: San Francisco, CA, USA, 2008; Volume 394, p. 654. [Google Scholar]
  12. Kaufman, M.J.; Wolfire, M.G.; Hollenbach, D.J. [Si II], [Fe II], [C II], and H2 Emission from Massive Star-forming Regions. Astrophys. J. 2006, 644, 283–299. [Google Scholar] [CrossRef]
  13. Ferland, G.J.; Chatzikos, M.; Guzmán, F.; Lykins, M.L.; van Hoof, P.A.M.; Williams, R.J.R.; Abel, N.P.; Badnell, N.R.; Keenan, F.P.; Porter, R.L.; et al. The 2017 Release Cloudy. Rev. Mexicana Astron. Astrofis. 2017, 53, 385–438. [Google Scholar]
  14. Röllig, M.; Abel, N.P.; Bell, T.; Bensch, F.; Black, J.; Ferland, G.J.; Jonkheid, B.; Kamp, I.; Kaufman, M.J.; Le Bourlot, J.; et al. A photon dominated region code comparison study. Astron. Astrophys. 2007, 467, 187–206. [Google Scholar] [CrossRef] [Green Version]
  15. Bolatto, A.D.; Jackson, J.M.; Ingalls, J.G. A Semianalytical Model for the Observational Properties of the Dominant Carbon Species at Different Metallicities. Astrophys. J. 1999, 513, 275–286. [Google Scholar] [CrossRef] [Green Version]
  16. Röllig, M.; Ossenkopf, V.; Jeyakumar, S.; Stutzki, J.; Sternberg, A. [CII] 158 μm emission and metallicity in photon dominated regions. Astron. Astrophys. 2006, 451, 917–924. [Google Scholar] [CrossRef] [Green Version]
  17. Narayanan, D.; Kulesa, C.A.; Boss, A.; Walker, C.K. Molecular Line Emission from Gravitationally Unstable Protoplanetary Disks. Astrophys. J. 2006, 647, 1426–1436. [Google Scholar] [CrossRef] [Green Version]
  18. Popping, G.; Pérez-Beaupuits, J.P.; Spaans, M.; Trager, S.C.; Somerville, R.S. The nature of the ISM in galaxies during the star-formation activity peak of the Universe. Mon. Not. Roy. Astron. Soc. 2014, 444, 1301–1317. [Google Scholar] [CrossRef] [Green Version]
  19. Vallini, L.; Gallerani, S.; Ferrara, A.; Pallottini, A.; Yue, B. On the [CII]-SFR Relation in High Redshift Galaxies. Astrophys. J. 2015, 813, 36. [Google Scholar] [CrossRef]
  20. Olsen, K.P.; Greve, T.R.; Narayanan, D.; Thompson, R.; Toft, S.; Brinch, C. Simulator of Galaxy Millimeter/Submillimeter Emission (SÍGAME): The [C ii]-SFR Relationship of Massive z = 2 Main Sequence Galaxies. Astrophys. J. 2015, 814, 76. [Google Scholar] [CrossRef]
  21. Popping, G.; van Kampen, E.; Decarli, R.; Spaans, M.; Somerville, R.S.; Trager, S.C. Sub-mm emission line deep fields: CO and [C II] luminosity functions out to z = 6. Mon. Not. Roy. Astron. Soc. 2016, 461, 93–110. [Google Scholar] [CrossRef]
  22. Popping, G.; Narayanan, D.; Somerville, R.S.; Faisst, A.L.; Krumholz, M.R. The art of modeling CO, [CI], and [CII] in cosmological galaxy formation models. arXiv, 2018; arXiv:1805.11093. [Google Scholar]
  23. Vallini, L.; Pallottini, A.; Ferrara, A.; Gallerani, S.; Sobacchi, E.; Behrens, C. CO line emission from galaxies in the Epoch of Reionization. Mon. Not. Roy. Astron. Soc. 2018, 473, 271–285. [Google Scholar] [CrossRef]
  24. Rosdahl, J.; Katz, H.; Blaizot, J.; Kimm, T.; Michel-Dansac, L.; Garel, T.; Haehnelt, M.; Ocvirk, P.; Teyssier, R. The SPHINX Cosmological Simulations of the First Billion Years: The Impact of Binary Stars on Reionization. arXiv, 2018; arXiv:1801.07259. [Google Scholar] [CrossRef]
  25. Keto, E.; Rybicki, G.B.; Bergin, E.A.; Plume, R. Radiative Transfer and Starless Cores. Astrophys. J. 2004, 613, 355–373. [Google Scholar] [CrossRef] [Green Version]
  26. Brinch, C.; Hogerheijde, M.R. LIME—A flexible, non-LTE line excitation and radiation transfer method for millimeter and far-infrared wavelengths. Astron. Astrophys. 2010, 523, A25. [Google Scholar] [CrossRef]
  27. Keto, E.; Rybicki, G. Modeling Molecular Hyperfine Line Emission. Astrophys. J. 2010, 716, 1315–1322. [Google Scholar] [CrossRef]
  28. Dullemond, C.P.; Juhasz, A.; Pohl, A.; Sereshti, F.; Shetty, R.; Peters, T.; Commercon, B.; Flock, M. RADMC-3D: A Multi-Purpose Radiative Transfer Tool; Astrophysics Source Code Library: Houghton, MI, USA, 2012. [Google Scholar]
  29. Yajima, H.; Li, Y.; Zhu, Q.; Abel, T. ART2: coupling Lyα line and multi-wavelength continuum radiative transfer. Mon. Not. Roy. Astron. Soc. 2012, 424, 884–901. [Google Scholar] [CrossRef]
  30. Ferland, G.J.; Porter, R.L.; van Hoof, P.A.M.; Williams, R.J.R.; Abel, N.P.; Lykins, M.L.; Shaw, G.; Henney, W.J.; Stancil, P.C. The 2013 Release of Cloudy. arXiv, 2013; arXiv:1705.10877. [Google Scholar]
  31. Krumholz, M.R. DESPOTIC—A new software library to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds. Mon. Not. Roy. Astron. Soc. 2014, 437, 1662–1680. [Google Scholar] [CrossRef]
  32. Gray, W.J.; Scannapieco, E.; Kasen, D. Atomic Chemistry in Turbulent Astrophysical Media. I. Effect of Atomic Cooling. Astrophys. J. 2015, 801, 107. [Google Scholar] [CrossRef]
  33. Kewley, L.J.; Dopita, M.A.; Leitherer, C.; Davé, R.; Yuan, T.; Allen, M.; Groves, B.; Sutherland, R. Theoretical Evolution of Optical Strong Lines across Cosmic Time. Astrophys. J. 2013, 774, 100. [Google Scholar] [CrossRef]
  34. Orsi, Á.; Padilla, N.; Groves, B.; Cora, S.; Tecce, T.; Gargiulo, I.; Ruiz, A. The nebular emission of star-forming galaxies in a hierarchical universe. Mon. Not. Roy. Astron. Soc. 2014, 443, 799–814. [Google Scholar] [CrossRef] [Green Version]
  35. Shimizu, I.; Inoue, A.K.; Okamoto, T.; Yoshida, N. Nebular line emission from z = 7 galaxies in a cosmological simulation: Rest-frame UV to optical lines. Mon. Not. Roy. Astron. Soc. 2016, 461, 3563–3575. [Google Scholar] [CrossRef]
  36. Hirschmann, M.; Charlot, S.; Feltre, A.; Naab, T.; Choi, E.; Ostriker, J.P.; Somerville, R.S. Synthetic nebular emission from massive galaxies–I: Origin of the cosmic evolution of optical emission-line ratios. Mon. Not. Roy. Astron. Soc. 2017, 472, 2468–2495. [Google Scholar] [CrossRef]
  37. Aravena, M.; Carilli, C.; Daddi, E.; Wagg, J.; Walter, F.; Riechers, D.; Dannerbauer, H.; Morrison, G.E.; Stern, D.; Krips, M. Cold Molecular Gas in Massive, Star-forming Disk Galaxies at z = 1.5. Astrophys. J. 2010, 718, 177–183. [Google Scholar] [CrossRef]
  38. Daddi, E.; Bournaud, F.; Walter, F.; Dannerbauer, H.; Carilli, C.L.; Dickinson, M.; Elbaz, D.; Morrison, G.E.; Riechers, D.; Onodera, M.; et al. Very High Gas Fractions and Extended Gas Reservoirs in z = 1.5 Disk Galaxies. Astrophys. J. 2010, 713, 686–707. [Google Scholar] [CrossRef]
  39. Tacconi, L.J.; Genzel, R.; Neri, R.; Cox, P.; Cooper, M.C.; Shapiro, K.; Bolatto, A.; Bouché, N.; Bournaud, F.; Burkert, A.; et al. High molecular gas fractions in normal massive star-forming galaxies in the young Universe. Nature 2010, 463, 781–784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Tacconi, L.J.; Neri, R.; Genzel, R.; Combes, F.; Bolatto, A.; Cooper, M.C.; Wuyts, S.; Bournaud, F.; Burkert, A.; Comerford, J.; et al. Phibss: Molecular Gas Content and Scaling Relations in z ~ 1–3 Massive, Main-sequence Star-forming Galaxies. Astrophys. J. 2013, 768, 74. [Google Scholar] [CrossRef]
  41. Tacconi, L.J.; Genzel, R.; Saintonge, A.; Combes, F.; García-Burillo, S.; Neri, R.; Bolatto, A.; Contini, T.; Förster Schreiber, N.M.; Lilly, S.; et al. PHIBSS: Unified Scaling Relations of Gas Depletion Time and Molecular Gas Fractions. Astrophys. J. 2018, 853, 179. [Google Scholar] [CrossRef] [Green Version]
  42. Geach, J.E.; Smail, I.; Moran, S.M.; MacArthur, L.A.; Lagos, C.d.P.; Edge, A.C. On the Evolution of the Molecular Gas Fraction of Star-Forming Galaxies. Astrophys. J. Lett. 2011, 730, L19. [Google Scholar] [CrossRef]
  43. Magdis, G.E.; Daddi, E.; Sargent, M.; Elbaz, D.; Gobat, R.; Dannerbauer, H.; Feruglio, C.; Tan, Q.; Rigopoulou, D.; Charmandaris, V.; et al. The Molecular Gas Content of z = 3 Lyman Break Galaxies: Evidence of a Non-evolving Gas Fraction in Main-sequence Galaxies at z > 2. Astrophys. J. Lett. 2012, 758, L9. [Google Scholar] [CrossRef]
  44. Santini, P.; Maiolino, R.; Magnelli, B.; Lutz, D.; Lamastra, A.; Li Causi, G.; Eales, S.; Andreani, P.; Berta, S.; Buat, V.; et al. The evolution of the dust and gas content in galaxies. Astron. Astrophys. 2014, 562, A30. [Google Scholar] [CrossRef] [Green Version]
  45. Béthermin, M.; De Breuck, C.; Gullberg, B.; Aravena, M.; Bothwell, M.S.; Chapman, S.C.; Gonzalez, A.H.; Greve, T.R.; Litke, K.; Ma, J.; et al. An ALMA view of the interstellar medium of the z = 4.77 lensed starburst SPT-S J213242-5802.9. Astron. Astrophys. 2016, 586, L7. [Google Scholar] [CrossRef]
  46. Decarli, R.; Walter, F.; Aravena, M.; Carilli, C.; Bouwens, R.; da Cunha, E.; Daddi, E.; Elbaz, D.; Riechers, D.; Smail, I.; et al. ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies. arXiv, 2016; arXiv:1607.06771. [Google Scholar]
  47. Scoville, N.; Sheth, K.; Aussel, H.; Vanden Bout, P.; Capak, P.; Bongiorno, A.; Casey, C.M.; Murchikova, L.; Koda, J.; Álvarez-Márquez, J.; et al. ISM Masses and the Star formation Law at Z = 1 to 6: ALMA Observations of Dust Continuum in 145 Galaxies in the COSMOS Survey Field. Astrophys. J. 2016, 820, 83. [Google Scholar] [CrossRef]
  48. da Cunha, E.; Groves, B.; Walter, F.; Decarli, R.; Weiss, A.; Bertoldi, F.; Carilli, C.; Daddi, E.; Elbaz, D.; Ivison, R.; Maiolino, R.; Riechers, D.; Rix, H.W.; Sargent, M.; Smail, I. On the Effect of the Cosmic Microwave Background in High-redshift (Sub-)millimeter Observations. Astrophys. J. 2013, 766, 13. [Google Scholar] [CrossRef]
  49. Bisbas, T.G.; Papadopoulos, P.P.; Viti, S. Effective Destruction of CO by Cosmic Rays: Implications for Tracing H2 Gas in the Universe. Astrophys. J. 2015, 803, 37. [Google Scholar] [CrossRef]
  50. Glover, S.C.O.; Clark, P.C. Is atomic carbon a good tracer of molecular gas in metal-poor galaxies? Mon. Not. Roy. Astron. Soc. 2016, 456, 3596–3609. [Google Scholar] [CrossRef] [Green Version]
  51. Bothwell, M.S.; Aguirre, J.E.; Aravena, M.; Bethermin, M.; Bisbas, T.G.; Chapman, S.C.; De Breuck, C.; Gonzalez, A.H.; Greve, T.R.; Hezaveh, Y.; et al. ALMA observations of atomic carbon in z ~ 4 dusty star-forming galaxies. arXiv, 2016; arXiv:1612.04380. [Google Scholar]
  52. Popping, G.; Decarli, R.; Man, A.W.S.; Nelson, E.J.; Béthermin, M.; De Breuck, C.; Mainieri, V.; van Dokkum, P.G.; Gullberg, B.; van Kampen, E.; et al. ALMA reveals starburst-like interstellar medium conditions in a compact star-forming galaxy at z ~ 2 using [CI] and CO. Astron. Astrophys. 2017, 602, A11. [Google Scholar] [CrossRef] [Green Version]
  53. Zanella, A.; Daddi, E.; Magdis, G.; Diaz Santos, T.; Cormier, D.; Liu, D.; Cibinel, A.; Gobat, R.; Dickinson, M.; Sargent, M.; et al. The [C II] emission as a molecular 738 gas mass tracer in galaxies at low and high redshift. arXiv, 2018; arXiv:1808.10331. [Google Scholar]
  54. Liu, L.; Weiß, A.; Perez-Beaupuits, J.P.; Güsten, R.; Liu, D.; Gao, Y.; Menten, K.M.; van derWerf, P.; Israel, F.P.; Harris, A.; et al. HIFI Spectroscopy of H2O Submillimeter 741 Lines in Nuclei of Actively Star-forming Galaxies. Astrophys. J. 2017, 846, 5. [Google Scholar] [CrossRef]
  55. Cormier, D.; Bigiel, F.; Jiménez-Donaire, M.J.; Leroy, A.K.; Gallagher, M.; Usero, A.; Sandstrom, K.; Bolatto, A.; Hughes, A.; Kramer, C.; et al. Full-disc 13CO(1-0) mapping across nearby galaxies of the EMPIRE survey and the CO-to-H2 conversion factor. Mon. Not. Roy. Astron. Soc. 2018, 475, 3909–3933. [Google Scholar] [CrossRef]
  56. Sansonetti, C.J.; Kerber, F.; Reader, J.; Rosa, M.R. Characterization of the Far-ultraviolet Spectrum of Pt/Cr-Ne Hollow Cathode Lamps as Used on the Space Telescope Imaging Spectrograph on Board the Hubble Space Telescope. Astrophys. J. Suppl. 2004, 153, 555–579. [Google Scholar] [CrossRef]
  57. Leitherer, C.; Tremonti, C.A.; Heckman, T.M.; Calzetti, D. An Ultraviolet Spectroscopic Atlas of Local Starbursts and Star-forming Galaxies: The Legacy of FOS and GHRS. Astron. J. 2011, 141, 37. [Google Scholar] [CrossRef]
  58. Mittal, R.; O’Dea, C.P.; Ferland, G.; Oonk, J.B.R.; Edge, A.C.; Canning, R.E.A.; Russell, H.; Baum, S.A.; Böhringer, H.; Combes, F.; et al. Herschel observations of the Centaurus cluster–the dynamics of cold gas in a cool core. Mon. Not. Roy. Astron. Soc. 2011, 418, 2386–2402. [Google Scholar] [CrossRef] [Green Version]
  59. Pallottini, A.; Ferrara, A.; Gallerani, S.; Vallini, L.; Maiolino, R.; Salvadori, S. Zooming on the internal structure of z ≃ 6 galaxies. Mon. Not. Roy. Astron. Soc. 2017, 465, 2540–2558. [Google Scholar] [CrossRef]
  60. Xiao, L.; Stanway, E.R.; Eldridge, J.J. Emission-line diagnostics of nearby H II regions including interacting binary populations. Mon. Not. Roy. Astron. Soc. 2018, 477, 904–934. [Google Scholar] [CrossRef]
  61. Safranek-Shrader, C.; Krumholz, M.R.; Kim, C.G.; Ostriker, E.C.; Klein, R.I.; Li, S.; McKee, C.F.; Stone, J.M. Chemistry and radiative shielding in star-forming galactic discs. Mon. Not. Roy. Astron. Soc. 2017, 465, 885–905. [Google Scholar] [CrossRef]
  62. Ginsburg, A.; Henkel, C.; Ao, Y.; Riquelme, D.; Kauffmann, J.; Pillai, T.; Mills, E.A.C.; Requena-Torres, M.A.; Immer, K.; Testi, L.; et al. Dense gas in the Galactic central molecular zone is warm and heated by turbulence. Astron. Astrophys. 2016, 586, A50. [Google Scholar] [CrossRef]
  63. Peñaloza, C.H.; Clark, P.C.; Glover, S.C.O.; Shetty, R.; Klessen, R.S. Using CO line ratios to trace the physical properties of molecular clouds. Mon. Not. Roy. Astron. Soc. 2017, 465, 2277–2285. [Google Scholar] [CrossRef]
  64. Peñaloza, C.H.; Clark, P.C.; Glover, S.C.O.; Klessen, R.S. CO line ratios in molecular clouds: The impact of environment. Mon. Not. Roy. Astron. Soc. 2018, 475, 1508–1520. [Google Scholar] [CrossRef]
  65. Li, Y.; Hopkins, P.F.; Hernquist, L.; Finkbeiner, D.P.; Cox, T.J.; Springel, V.; Jiang, L.; Fan, X.; Yoshida, N. Modeling the Dust Properties of z ∼ 6 Quasars with ART2 - All-Wavelength Radiative Transfer with Adaptive Refinement Tree. Astrophys. J. 2008, 678, 41–63. [Google Scholar] [CrossRef]
  66. Yajima, H.; Li, Y.; Zhu, Q.; Abel, T.; Gronwall, C.; Ciardullo, R. Escape of Lyα and continuum photons from star-forming galaxies. Mon. Not. Roy. Astron. Soc. 2014, 440, 776–786. [Google Scholar] [CrossRef] [Green Version]
  67. Gray, W.J.; Scannapieco, E. Atomic Chemistry in Turbulent Astrophysical Media. II. Effect of the Redshift Zero Metagalactic Background. Astrophys. J. 2016, 818, 198. [Google Scholar] [CrossRef]
  68. Gray, W.J.; Scannapieco, E. The Effect of Turbulence on Nebular Emission Line Ratios. Astrophys. J. 2017, 849, 132. [Google Scholar] [CrossRef] [Green Version]
  69. Dere, K.P.; Landi, E.; Mason, H.E.; Monsignori Fossi, B.C.; Young, P.R. CHIANTI–an atomic database for emission lines. Astron. Astrophys. Suppl. 1997, 125, 149–173. [Google Scholar] [CrossRef]
  70. Del Zanna, G.; Dere, K.P.; Young, P.R.; Landi, E.; Mason, H.E. CHIANTI–An atomic database for emission lines. Version 8. Astron. Astrophys. 2015, 582, A56. [Google Scholar] [CrossRef]
  71. Brickhouse, N.S.; Dupree, A.K.; Edgar, R.J.; Liedahl, D.A.; Drake, S.A.; White, N.E.; Singh, K.P. Coronal Structure and Abundances of Capella from Simultaneous EUVE and ASCA Spectroscopy. Astrophys. J. 2000, 530, 387–402. [Google Scholar] [CrossRef]
  72. Kaastra, J.S.; Bleeker, J.A.M.; Mewe, R. X-ray diagnostics of supernova remnants. In UV and X-ray Spectroscopy of Astrophysical and Laboratory Plasmas; Yamashita, K., Watanabe, T., Eds.; Universal Academy Press: Tokyo, Japan, 1996; pp. 15–20. [Google Scholar]
  73. Bar-Shalom, A.; Oreg, J.; Klapisch, M. Recent developments in the SCROLL model. J. Quant. Spec. Radiat. Transf. 2000, 65, 43–55. [Google Scholar] [CrossRef]
  74. Del Zanna, G.; Berrington, K.A.; Mason, H.E. Benchmarking atomic data for astrophysics: Fe X. Astron. Astrophys. 2004, 422, 731–749. [Google Scholar] [CrossRef]
  75. Keenan, F.P.; Drake, J.J.; Aggarwal, K.M. An investigation of FeXVI emission lines in solar and stellar extreme-ultraviolet and soft X-ray spectra. Mon. Not. Roy. Astron. Soc. 2007, 381, 1727–1732. [Google Scholar] [CrossRef] [Green Version]
  76. Kato, D.; Sakaue, H.A.; Murakami, I.; Kato, T.; Nakamura, N.; Ohtani, S.; Yamamoto, N.; Watanabe, T. Electron-impact excitation of 2p53l → 2p53l’ line emission of Fe XVII. J. Phys. Conf. Ser. 2009, 163. [Google Scholar] [CrossRef]
  77. Beiersdorfer, P.; Träbert, E.; Lepson, J.K.; Brickhouse, N.S.; Golub, L. High-resolution Laboratory Measurements of Coronal Lines in the 198-218 Å Region. Astrophys. J. 2014, 788, 25. [Google Scholar]
  78. Liang, G.Y.; Whiteford, A.D.; Badnell, N.R. R-matrix inner-shell electron-impact excitation of the Na-like iso-electronic sequence. J. Phys. B At. Mol. Phys. 2009, 42, 225002. [Google Scholar] [CrossRef]
  79. Liang, G.Y.; Badnell, N.R. R-matrix electron-impact excitation data for the Ne-like iso-electronic sequence. Astron. Astrophys. 2010, 518, A64. [Google Scholar] [CrossRef]
  80. Wang, K.; Li, D.F.; Liu, H.T.; Han, X.Y.; Duan, B.; Li, C.Y.; Li, J.G.; Guo, X.L.; Chen, C.Y.; Yan, J. Systematic Calculations of Energy Levels and Transition Rates of C-like Ions with Z = 13–36. Astrophys. J. Suppl. 2014, 215, 26. [Google Scholar] [CrossRef]
  81. Del Zanna, G.; Badnell, N.R. Atomic data for astrophysics: Ni XII. Astron. Astrophys. 2016, 585, A118. [Google Scholar] [CrossRef]
  82. Liedahl, D.A.; Osterheld, A.L.; Goldstein, W.H. New calculations of Fe L-shell X-ray spectra in high-temperature plasmas. Astrophys. J. Lett. 1995, 438, L115–L118. [Google Scholar] [CrossRef]
  83. Gu, M.F. The flexible atomic code. Can. J. Phys. 2008, 86, 675–689. [Google Scholar] [CrossRef]
  84. Dyall, K.G.; Grant, I.P.; Johnson, C.T.; Parpia, F.A.; Plummer, E.P. GRASP: A general-purpose relativistic atomic structure program. Comput. Phys. Commun. 1989, 55, 425–456. [Google Scholar] [CrossRef]
  85. Badnell, N.R. Dielectric recombination of Fe(22+) and Fe(21+). J. Phys. B At. Mol. Phys. 1986, 19, 3827–3835. [Google Scholar] [CrossRef]
  86. Aggarwal, K.M.; Keenan, F.P. Energy levels, radiative rates and electron impact excitation rates for transitions in Fe XIV. Mon. Not. Roy. Astron. Soc. 2014, 445, 2015–2027. [Google Scholar] [CrossRef] [Green Version]
  87. Del Zanna, G.; Badnell, N.R. Atomic data and density diagnostics for S IV. Mon. Not. Roy. Astron. Soc. 2016, 456, 3720–3728. [Google Scholar] [CrossRef]
  88. Landi, E.; Del Zanna, G.; Young, P.R.; Dere, K.P.; Mason, H.E. CHIANTI–An Atomic Database for Emission Lines. XII. Version 7 of the Database. Astrophys. J. 2012, 744, 99. [Google Scholar] [CrossRef]
  89. Del Zanna, G.; Ishikawa, Y. Benchmarking atomic data for astrophysics: Fe XVII EUV lines. Astron. Astrophys. 2009, 508, 1517–1526. [Google Scholar] [CrossRef] [Green Version]
  90. Brown, G.V.; Beiersdorfer, P.; Liedahl, D.A.; Widmann, K.; Kahn, S.M. Laboratory Measurements and Modeling of the Fe XVII X-Ray Spectrum. Astrophys. J. 1998, 502, 1015–1026. [Google Scholar] [CrossRef]
  91. Kramida, A.; Ralchenko, Y.; Reader, J.; NIST ASD Team. NIST Atomic Spectra Database (ver. 5.2). National Institute of Standards and Technology: Gaithersburg, MD, USA, 2014. Available online: http://physics.nist.gov/asd (accessed on 28 November 2014).
  92. Steinacker, J.; Baes, M.; Gordon, K.D. Three-Dimensional Dust Radiative Transfer*. Annu. Rev. Astron. Astrophys. 2013, 51, 63–104. [Google Scholar] [CrossRef] [Green Version]
  93. Gnedin, N.Y.; Kaurov, A.A. Cosmic Reionization on Computers. II. Reionization History and Its Back-reaction on Early Galaxies. Astrophys. J. 2014, 793, 30. [Google Scholar] [CrossRef]
  94. So, G.C.; Norman, M.L.; Reynolds, D.R.; Wise, J.H. Fully Coupled Simulation of Cosmic Reionization. II. Recombinations, Clumping Factors, and the Photon Budget for Reionization. Astrophys. J. 2014, 789, 149. [Google Scholar] [CrossRef]
  95. Wise, J.H.; Demchenko, V.G.; Halicek, M.T.; Norman, M.L.; Turk, M.J.; Abel, T.; Smith, B.D. The birth of a galaxy–III. Propelling reionization with the faintest galaxies. Mon. Not. Roy. Astron. Soc. 2014, 442, 2560–2579. [Google Scholar] [CrossRef]
  96. Norman, M.L.; Reynolds, D.R.; So, G.C.; Harkness, R.P.; Wise, J.H. Fully Coupled Simulation of Cosmic Reionization. I. Numerical Methods and Tests. Astrophys. J. Suppl. 2015, 216, 16. [Google Scholar] [CrossRef]
  97. Pawlik, A.H.; Schaye, J.; Dalla Vecchia, C. Spatially adaptive radiation-hydrodynamical simulations of galaxy formation during cosmological reionization. Mon. Not. Roy. Astron. Soc. 2015, 451, 1586–1605. [Google Scholar] [CrossRef] [Green Version]
  98. Ocvirk, P.; Gillet, N.; Shapiro, P.R.; Aubert, D.; Iliev, I.T.; Teyssier, R.; Yepes, G.; Choi, J.H.; Sullivan, D.; Knebe, A.; et al. Cosmic Dawn (CoDa): The First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe. Mon. Not. Roy. Astron. Soc. 2016, 463, 1462–1485. [Google Scholar] [CrossRef]
  99. Rosdahl, J.; Blaizot, J.; Aubert, D.; Stranex, T.; Teyssier, R. RAMSES-RT: Radiation hydrodynamics in the cosmological context. Mon. Not. Roy. Astron. Soc. 2013, 436, 2188–2231. [Google Scholar] [CrossRef]
  100. Rosdahl, J.; Teyssier, R. A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT. Mon. Not. Roy. Astron. Soc. 2015, 449, 4380–4403. [Google Scholar] [CrossRef]
  101. Aubert, D.; Deparis, N.; Ocvirk, P. EMMA: An adaptive mesh refinement cosmological simulation code with radiative transfer. Mon. Not. Roy. Astron. Soc. 2015, 454, 1012–1037. [Google Scholar] [CrossRef]
  102. Krumholz, M.R.; Klein, R.I.; McKee, C.F.; Bolstad, J. Equations and Algorithms for Mixed-frame Flux-limited Diffusion Radiation Hydrodynamics. Astrophys. J. 2007, 667, 626–643. [Google Scholar] [CrossRef]
  103. Davis, S.W.; Stone, J.M.; Jiang, Y.F. A Radiation Transfer Solver for Athena Using Short Characteristics. Astrophys. J. Suppl. 2012, 199, 9. [Google Scholar] [CrossRef]
  104. Elitzur, M.; Asensio Ramos, A. A new exact method for line radiative transfer. Mon. Not. Roy. Astron. Soc. 2006, 365, 779–791. [Google Scholar] [CrossRef] [Green Version]
  105. Hubeny, I. From Escape Probabilities to Exact Radiative Transfer. In Spectroscopic Challenges of Photoionized Plasmas; Astronomical Society of the Pacific Conference Series; Ferland, G., Savin, D.W., Eds.; Astronomical Society: San Francisco, CA, USA, 2001; Volume 247, p. 197. [Google Scholar]
  106. Evans, N.J., II. Physical Conditions in Regions of Star Formation. Annu. Rev. Astron. Astrophys. 1999, 37, 311–362. [Google Scholar] [CrossRef]
  107. Shu, F.H. Self-similar collapse of isothermal spheres and star formation. Astrophys. J. 1977, 214, 488–497. [Google Scholar] [CrossRef]
  108. Keto, E.; Caselli, P.; Rawlings, J. The dynamics of collapsing cores and star formation. Mon. Not. Roy. Astron. Soc. 2015, 446, 3731–3740. [Google Scholar] [CrossRef]
  109. Eldridge, J.J.; Stanway, E.R.; Xiao, L.; McClelland, L.A.S.; Taylor, G.; Ng, M.; Greis, S.M.L.; Bray, J.C. Binary Population and Spectral Synthesis Version 2.1: Construction, Observational Verification, and New Results. Publ. Astron. Soc. Aust. 2017, 34, e058. [Google Scholar] [CrossRef]
  110. Wofford, A.; Charlot, S.; Bruzual, G.; Eldridge, J.J.; Calzetti, D.; Adamo, A.; Cignoni, M.; de Mink, S.E.; Gouliermis, D.A.; Grasha, K.; et al. A comprehensive comparative test of seven widely used spectral synthesis models against multi-band photometry of young massive-star clusters. Mon. Not. Roy. Astron. Soc. 2016, 457, 4296–4322. [Google Scholar] [CrossRef] [Green Version]
  111. Izotov, Y.I.; Schaerer, D.; Thuan, T.X.; Worseck, G.; Guseva, N.G.; Orlitová, I.; Verhamme, A. Detection of high Lyman continuum leakage from four low-redshift compact star-forming galaxies. Mon. Not. Roy. Astron. Soc. 2016, 461, 3683–3701. [Google Scholar] [CrossRef] [Green Version]
  112. Kroupa, P. On the variation of the initial mass function. Mon. Not. Roy. Astron. Soc. 2001, 322, 231–246. [Google Scholar] [CrossRef] [Green Version]
  113. Meynet, G.; Maeder, A.; Schaller, G.; Schaerer, D.; Charbonnel, C. Grids of massive stars with high mass loss rates. V. From 12 to 120 Msun_ at Z = 0.001, 0.004, 0.008, 0.020 and 0.040. Astron. Astrophys. Suppl. 1994, 103, 97–105. [Google Scholar]
  114. Habing, H.J. The interstellar radiation density between 912 A and 2400 A. Bull. Astron. Inst. Neth. 1968, 19, 421. [Google Scholar]
  115. Veilleux, S.; Osterbrock, D.E. Spectral classification of emission-line galaxies. Astrophys. J. Suppl. 1987, 63, 295–310. [Google Scholar] [CrossRef]
  116. Revalski, M.; Crenshaw, D.M.; Kraemer, S.B.; Fischer, T.C.; Schmitt, H.R.; Machuca, C. Quantifying Feedback from Narrow Line Region Outflows in Nearby Active Galaxies. I. Spatially Resolved Mass Outflow Rates for the Seyfert 2 Galaxy Markarian 573. Astrophys. J. 2018, 856, 46. [Google Scholar] [CrossRef] [Green Version]
  117. Sutherland, R.S.; Dopita, M.A. Effects of Preionization in Radiative Shocks. I. Self-consistent Models. Astrophys. J. Suppl. 2017, 229, 34. [Google Scholar] [CrossRef] [Green Version]
  118. Ballesteros-Paredes, J.; Vázquez-Semadeni, E.; Scalo, J. Clouds as Turbulent Density Fluctuations: Implications for Pressure Confinement and Spectral Line Data Interpretation. Astrophys. J. 1999, 515, 286–303. [Google Scholar] [CrossRef] [Green Version]
  119. Hennebelle, P.; Pérault, M. Dynamical condensation in a thermally bistable flow. Application to interstellar cirrus. Astron. Astrophys. 1999, 351, 309–322. [Google Scholar]
  120. Koyama, H.; Inutsuka, S.I. An Origin of Supersonic Motions in Interstellar Clouds. Astrophys. J. Lett. 2002, 564, L97–L100. [Google Scholar] [CrossRef]
  121. Audit, E.; Hennebelle, P. Thermal condensation in a turbulent atomic hydrogen flow. Astron. Astrophys. 2005, 433, 1–13. [Google Scholar] [CrossRef] [Green Version]
  122. Heitsch, F.; Burkert, A.; Hartmann, L.W.; Slyz, A.D.; Devriendt, J.E.G. Formation of Structure in Molecular Clouds: A Case Study. Astrophys. J. Lett. 2005, 633, L113–L116. [Google Scholar] [CrossRef]
  123. Vázquez-Semadeni, E.; Gómez, G.C.; Jappsen, A.K.; Ballesteros-Paredes, J.; González, R.F.; Klessen, R.S. Molecular Cloud Evolution. II. From Cloud Formation to the Early Stages of Star Formation in Decaying Conditions. Astrophys. J. 2007, 657, 870–883. [Google Scholar] [CrossRef] [Green Version]
  124. Banerjee, R.; Vázquez-Semadeni, E.; Hennebelle, P.; Klessen, R.S. Clump morphology and evolution in MHD simulations of molecular cloud formation. Mon. Not. Roy. Astron. Soc. 2009, 398, 1082–1092. [Google Scholar] [CrossRef] [Green Version]
  125. Dobbs, C.L.; Pringle, J.E.; Burkert, A. Giant molecular clouds: what are they made from, and how do they get there? Mon. Not. Roy. Astron. Soc. 2012, 425, 2157–2168. [Google Scholar] [CrossRef] [Green Version]
  126. Olsen, K.; Greve, T.R.; Narayanan, D.; Thompson, R.; Davé, R.; Niebla Rios, L.; Stawinski, S. SÍGAME Simulations of the [CII], [OI], and [OIII] Line Emission from Star-forming Galaxies at z ≃ 6. Astrophys. J. 2017, 846, 105. [Google Scholar] [CrossRef]
  127. Pallottini, A.; Ferrara, A.; Bovino, S.; Vallini, L.; Gallerani, S.; Maiolino, R.; Salvadori, S. The impact of chemistry on the structure of high-z galaxies. Mon. Not. Roy. Astron. Soc. 2017, 471, 4128–4143. [Google Scholar] [CrossRef] [Green Version]
  128. Capelo, P.R.; Bovino, S.; Lupi, A.; Schleicher, D.R.G.; Grassi, T. The effect of non-equilibrium metal cooling on the interstellar medium. Mon. Not. Roy. Astron. Soc. 2018, 475, 3283–3304. [Google Scholar] [CrossRef] [Green Version]
  129. Rosdahl, J.; Schaye, J.; Teyssier, R.; Agertz, O. Galaxies that shine: Radiation-hydrodynamical simulations of disc galaxies. Mon. Not. Roy. Astron. Soc. 2015, 451, 34–58. [Google Scholar] [CrossRef]
  130. Katz, H.; Kimm, T.; Sijacki, D.; Haehnelt, M.G. Interpreting ALMA observations of the ISM during the epoch of reionization. Mon. Not. Roy. Astron. Soc. 2017, 468, 4831–4861. [Google Scholar] [CrossRef]
  131. Grassi, T.; Bovino, S.; Haugbølle, T.; Schleicher, D.R.G. A detailed framework to incorporate dust in hydrodynamical simulations. Mon. Not. Roy. Astron. Soc. 2017, 466, 1259–1274. [Google Scholar] [CrossRef]
  132. Popping, G.; Somerville, R.S.; Galametz, M. The dust content of galaxies from z = 0 to z = 9. Mon. Not. Roy. Astron. Soc. 2017, 471, 3152–3185. [Google Scholar] [CrossRef] [Green Version]
  133. McKinnon, R.; Torrey, P.; Vogelsberger, M.; Hayward, C.C.; Marinacci, F. Simulating the dust content of galaxies: successes and failures. Mon. Not. Roy. Astron. Soc. 2017, 468, 1505–1521. [Google Scholar] [CrossRef]
  134. Aoyama, S.; Hou, K.C.; Hirashita, H.; Nagamine, K.; Shimizu, I. Cosmological simulation with dust formation and destruction. arXiv, 2018; arXiv:1802.04027. [Google Scholar] [CrossRef]
  135. Behrens, C.; Pallottini, A.; Ferrara, A.; Gallerani, S.; Vallini, L. Dusty galaxies in the Epoch of Reionization: simulations. Mon. Not. Roy. Astron. Soc. 2018, 477, 552–565. [Google Scholar] [CrossRef] [Green Version]
  136. Narayanan, D.; Conroy, C.; Dave, R.; Johnson, B.; Popping, G. A Theory for the Variation of Dust Attenuation Laws in Galaxies. arXiv, 2018; arXiv:1805.06905. [Google Scholar]
  137. Shirazi, M.; Brinchmann, J. Strongly star forming galaxies in the local Universe with nebular He IIλ4686 emission. Mon. Not. Roy. Astron. Soc. 2012, 421, 1043–1063. [Google Scholar] [CrossRef]
  138. Gutkin, J.; Charlot, S.; Bruzual, G. Modelling the nebular emission from primeval to present-day star-forming galaxies. Mon. Not. Roy. Astron. Soc. 2016, 462, 1757–1774. [Google Scholar] [CrossRef] [Green Version]
  139. Senchyna, P.; Stark, D.P.; Vidal-García, A.; Chevallard, J.; Charlot, S.; Mainali, R.; Jones, T.; Wofford, A.; Feltre, A.; Gutkin, J. Ultraviolet spectra of extreme nearby star-forming regions-approaching a local reference sample for JWST. Mon. Not. Roy. Astron. Soc. 2017, 472, 2608–2632. [Google Scholar] [CrossRef]
  140. Kalirai, J. Scientific discovery with the JamesWebb Space Telescope. Contemp. Phys. 2018, 59, 251–290. [Google Scholar] [CrossRef]
  141. Skidmore, W.; TMT International Science Development Teams; Science Advisory Committee. Thirty Meter Telescope Detailed Science Case: 2015. Res. Astron. Astrophys. 2015, 15, 1945. [Google Scholar] [CrossRef]
  142. Gallerani, S.; Ferrara, A.; Neri, R.; Maiolino, R. First CO(17-16) emission line detected in a z > 6 quasar. Mon. Not. Roy. Astron. Soc. 2014, 445, 2848–2853. [Google Scholar] [CrossRef]
  143. Meijerink, R.; Spaans, M. Diagnostics of irradiated gas in galaxy nuclei. I. A far-ultraviolet and X-ray dominated region code. Astron. Astrophys. 2005, 436, 397–409. [Google Scholar] [CrossRef]
  144. Mingozzi, M.; Vallini, L.; Pozzi, F.; Vignali, C.; Mignano, A.; Gruppioni, C.; Talia, M.; Cimatti, A.; Cresci, G.; Massardi, M. CO excitation in the Seyfert galaxy NGC 34: stars, shock or AGN driven? Mon. Not. Roy. Astron. Soc. 2018, 474, 3640–3648. [Google Scholar] [CrossRef]
  145. Crenshaw, D.M.; Kraemer, S.B.; George, I.M. Mass Loss from the Nuclei of Active Galaxies. Annu Rev. 2003, 41, 117–167. [Google Scholar] [CrossRef] [Green Version]
  146. Hopkins, P.F.; Hernquist, L.; Cox, T.J.; Di Matteo, T.; Martini, P.; Robertson, B.; Springel, V. Black Holes in Galaxy Mergers: Evolution of Quasars. Astrophys. J. 2005, 630, 705–715. [Google Scholar] [CrossRef] [Green Version]
  147. Kormendy, J.; Ho, L.C. Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies. Annu. Rev. Astron. Astrophys. 2013, 51, 511–653. [Google Scholar] [CrossRef] [Green Version]
  148. Batiste, M.; Bentz, M.C.; Raimundo, S.I.; Vestergaard, M.; Onken, C.A. Recalibration of the MBH-σ Relation for AGN. Astrophys. J. Lett. 2017, 838, L10. [Google Scholar] [CrossRef]
  149. Fiore, F.; Feruglio, C.; Shankar, F.; Bischetti, M.; Bongiorno, A.; Brusa, M.; Carniani, S.; Cicone, C.; Duras, F.; Lamastra, A.; et al. AGN wind scaling relations and the co-evolution of black holes and galaxies. Astron. Astrophys. 2017, 601, A143. [Google Scholar] [CrossRef] [Green Version]
  150. Choi, E.; Ostriker, J.P.; Naab, T.; Somerville, R.S.; Hirschmann, M.; Núñez, A.; Hu, C.Y.; Oser, L. Physics of Galactic Metals: Evolutionary Effects due to Production, Distribution, Feedback, and Interaction with Black Holes. arXiv, 2016; arXiv:1610.09389. [Google Scholar]
  151. Richings, A.J.; Faucher-Giguère, C.A. The origin of fast molecular outflows in quasars: Molecule formation in AGN-driven galactic winds. Mon. Not. Roy. Astron. Soc. 2018, 474, 3673–3699. [Google Scholar] [CrossRef]
  152. Tumlinson, J.; Peeples, M.S.; Werk, J.K. The Circumgalactic Medium. Annu. Rev. Astron. Astrophys. 2017, 55, 389–432. [Google Scholar] [CrossRef] [Green Version]
  153. Suresh, J.; Rubin, K.H.R.; Kannan, R.; Werk, J.K.; Hernquist, L.; Vogelsberger, M. On the OVI abundance in the circumgalactic medium of low-redshift galaxies. Mon. Not. Roy. Astron. Soc. 2017, 465, 2966–2982. [Google Scholar] [CrossRef]
  154. Hummels, C.B.; Smith, B.D.; Silvia, D.W. Trident: A Universal Tool for Generating Synthetic Absorption Spectra from Astrophysical Simulations. Astrophys. J. 2017, 847, 59. [Google Scholar] [CrossRef] [Green Version]
  155. Dubois, Y.; Teyssier, R. On the onset of galactic winds in quiescent star forming galaxies. Astron. Astrophys. 2008, 477, 79–94. [Google Scholar] [CrossRef]
  156. Padoan, P.; Nordlund, Å. The Star Formation Rate of Supersonic Magnetohydrodynamic Turbulence. Astrophys. J. 2011, 730, 40. [Google Scholar] [CrossRef]
  157. Federrath, C.; Klessen, R.S. The Star Formation Rate of Turbulent Magnetized Clouds: Comparing Theory, Simulations, and Observations. Astrophys. J. 2012, 761, 156. [Google Scholar] [CrossRef]
  158. Kimm, T.; Cen, R. Escape Fraction of Ionizing Photons during Reionization: Effects due to Supernova Feedback and Runaway OB Stars. Astrophys. J. 2014, 788, 121. [Google Scholar] [CrossRef]
  159. Werk, J.K.; Prochaska, J.X.; Cantalupo, S.; Fox, A.J.; Oppenheimer, B.; Tumlinson, J.; Tripp, T.M.; Lehner, N.; McQuinn, M. The COS-Halos Survey: Origins of the Highly Ionized Circumgalactic Medium of Star-Forming Galaxies. Astrophys. J. 2016, 833, 54. [Google Scholar] [CrossRef]
  160. McCourt, M.; Oh, S.P.; O’Leary, R.; Madigan, A.M. A characteristic scale for cold gas. Mon. Not. Roy. Astron. Soc. 2018, 473, 5407–5431. [Google Scholar] [CrossRef]
  161. Oppenheimer, B.D.; Segers, M.; Schaye, J.; Richings, A.J.; Crain, R.A. Flickering AGN can explain the strong circumgalactic O VI observed by COS-Halos. Mon. Not. Roy. Astron. Soc. 2018, 474, 4740–4755. [Google Scholar] [CrossRef]
  162. Del Zanna, G.; DeLuca, E.E. Solar Coronal Lines in the Visible and Infrared: A Rough Guide. Astrophys. J. 2018, 852, 52. [Google Scholar] [CrossRef] [Green Version]
1
“Baldwin, Phillips & Terlevich” (BPT) diagrams [1].
2
3
4
5
Slides from and video recordings of talks can be found at https://zenodo.org/communities/walk2018/.
6
E.g., the Grand Challenge Problems in Computational Astrophysics conference series at https://www.ipam.ucla.edu/.
7
8
The Large UV/Optical/IR Surveyor: https://asd.gsfc.nasa.gov/luvoir/.
9
Habitable Exoplanet Imaging Mission: https://www.jpl.nasa.gov/habex/.
10
11
12
13
14
15
16
17
18
Figure 1. (Top) Sketch of standard PDR structure, from innovative work in the 80’s [4] [©AAS. Reproduced with permission]. (Bottom) For comparison, an example of PDR structure computed this year with version C17 of Cloudy [13] for model V1 of [14], adapted to extend to A V = 100 . (Left panel) Densities of important atoms and molecules in the PDR. (Right panel) Local emissivities for some of the most important coolants.
Figure 1. (Top) Sketch of standard PDR structure, from innovative work in the 80’s [4] [©AAS. Reproduced with permission]. (Bottom) For comparison, an example of PDR structure computed this year with version C17 of Cloudy [13] for model V1 of [14], adapted to extend to A V = 100 . (Left panel) Densities of important atoms and molecules in the PDR. (Right panel) Local emissivities for some of the most important coolants.
Galaxies 06 00100 g001
Figure 2. Illustration of the difference between simulating line emission from a cell of constant density, temperature, metallicity etc., and using an ensemble of clouds superimposed to create a more realistic synthetic observation of a patch of the sky.
Figure 2. Illustration of the difference between simulating line emission from a cell of constant density, temperature, metallicity etc., and using an ensemble of clouds superimposed to create a more realistic synthetic observation of a patch of the sky.
Galaxies 06 00100 g002
Figure 3. Theoretical and experimental energies for levels obtained from the CHIANTI (v7, [88]) database. (Top) Fe 16 + : Observed energy levels (up to n = 5): [89]; observed energy levels (n = 6,7) from the wavelength measurements of [90]. (Bottom) Mg + : theoretical energies: [78]; observed energies from NIST: [91].
Figure 3. Theoretical and experimental energies for levels obtained from the CHIANTI (v7, [88]) database. (Top) Fe 16 + : Observed energy levels (up to n = 5): [89]; observed energy levels (n = 6,7) from the wavelength measurements of [90]. (Bottom) Mg + : theoretical energies: [78]; observed energies from NIST: [91].
Galaxies 06 00100 g003
Figure 4. Example of SEDs for stellar populations of increasing age and metallicity, calculated with the starburst99 tool (see the text). A Kroupa IMF [112] at masses 0 . 1 –100 M was used together with an instantaneous burst of star formation which happened at age = 0 Myr. The Geneva high-mass loss tracks [113] were used. The shaded area indicates the region in wavelength used to integrate the SED for obtaining the FUV flux in e.g., Habing units [114], setting the chemical state of the gas. We note that comparing these models of different metallicities at the same age can be misleading, because stars of different metallicities have different lifetimes, i.e., a star of 3 Myr with solar metallicity will not be in the same evolutionary stage as a star of 3 Myr with 0.1 times solar metallicity.
Figure 4. Example of SEDs for stellar populations of increasing age and metallicity, calculated with the starburst99 tool (see the text). A Kroupa IMF [112] at masses 0 . 1 –100 M was used together with an instantaneous burst of star formation which happened at age = 0 Myr. The Geneva high-mass loss tracks [113] were used. The shaded area indicates the region in wavelength used to integrate the SED for obtaining the FUV flux in e.g., Habing units [114], setting the chemical state of the gas. We note that comparing these models of different metallicities at the same age can be misleading, because stars of different metallicities have different lifetimes, i.e., a star of 3 Myr with solar metallicity will not be in the same evolutionary stage as a star of 3 Myr with 0.1 times solar metallicity.
Galaxies 06 00100 g004
Figure 5. An illustration of how the gas density structure changes with increasing levels of turbulence, σ , made with MAIHEM and presented in the talk by William Gray. At low σ the turbulence is subsonic and no shocks are formed creating a nearly uniform density distribution. At intermediate and high σ , shocks form and create regions of very high and low density.
Figure 5. An illustration of how the gas density structure changes with increasing levels of turbulence, σ , made with MAIHEM and presented in the talk by William Gray. At low σ the turbulence is subsonic and no shocks are formed creating a nearly uniform density distribution. At intermediate and high σ , shocks form and create regions of very high and low density.
Galaxies 06 00100 g005
Figure 6. The effects of turbulence on the classic [N II] and [S II] BPT diagrams. These figures are as presented at the workshop, see [68] for more details.
Figure 6. The effects of turbulence on the classic [N II] and [S II] BPT diagrams. These figures are as presented at the workshop, see [68] for more details.
Galaxies 06 00100 g006
Figure 7. Average [OIII]/H β emission-line ratio in bins of [NII]/H α for the star-forming subset (i.e., with log (BHAR/SFR) < 4 ) of the massive galaxies simulated with (right panel) and without AGN feedback (left panel) in different redshift ranges (connected stars of different colors). To guide the eye, SDSS data are shown in grey (adapted from [36]).
Figure 7. Average [OIII]/H β emission-line ratio in bins of [NII]/H α for the star-forming subset (i.e., with log (BHAR/SFR) < 4 ) of the massive galaxies simulated with (right panel) and without AGN feedback (left panel) in different redshift ranges (connected stars of different colors). To guide the eye, SDSS data are shown in grey (adapted from [36]).
Galaxies 06 00100 g007
Figure 8. From Kolwa et al. in prep. Line emission from gas in the central 0.6 arcsec (4.75 kpc) of MRC 0943-242. This gas is ionized by the local AGN radiation field. Line widths larger than 1000 km s 1 suggest a turbulent medium. For the top two figures the red profiles are best-fit single Gaussians. For the bottom two figures the doublets are fit with identical redshift. For the bottom two figures, the velocity scale of the bottom (top) horizontal axis is set to zero at the blue (orange) Gaussian centroid.
Figure 8. From Kolwa et al. in prep. Line emission from gas in the central 0.6 arcsec (4.75 kpc) of MRC 0943-242. This gas is ionized by the local AGN radiation field. Line widths larger than 1000 km s 1 suggest a turbulent medium. For the top two figures the red profiles are best-fit single Gaussians. For the bottom two figures the doublets are fit with identical redshift. For the bottom two figures, the velocity scale of the bottom (top) horizontal axis is set to zero at the blue (orange) Gaussian centroid.
Galaxies 06 00100 g008
Figure 9. From Kolwa et al. in prep. Line emission in the same spectra as Figure 8, however with clear components of absorption as well. The absorbing components’ centroids are marked with vertical lines. Absorbing components 1, 2a, and 2b are seen in both the Ly α (left) and CIV (right) emission. The best fit models (red line) result from the Gaussian emission with multiple Gaussian absorption profiles. For the right figure, the velocity scale of the bottom (top) horizontal axis is set to zero at the blue (orange) Gaussian centroid, and each centroid has its absorption component centroids marked with corresponding font color.
Figure 9. From Kolwa et al. in prep. Line emission in the same spectra as Figure 8, however with clear components of absorption as well. The absorbing components’ centroids are marked with vertical lines. Absorbing components 1, 2a, and 2b are seen in both the Ly α (left) and CIV (right) emission. The best fit models (red line) result from the Gaussian emission with multiple Gaussian absorption profiles. For the right figure, the velocity scale of the bottom (top) horizontal axis is set to zero at the blue (orange) Gaussian centroid, and each centroid has its absorption component centroids marked with corresponding font color.
Galaxies 06 00100 g009
Figure 10. CGM OVI covering fractions in a sub-L* simulated galaxy at z = 3 for different subgrid models. SF#FB# corresponds to whether the simulations uses star formation model 1 or 2 (SF1 and SF2), and stellar feedback model 1 or 2 (FB1 and FB2). SF1 uses a density threshold for triggering star formation. FB1 corresponds to feedback from a stellar population with a Salpeter IMF and the massive stars exploding with 10 51 erg per 10 M and the energy injected in the simulation matching the Sedov solution [155]. Both of these subgrid models depend on the resolution of the simulations (here 10 pc). SF2 uses a variable star formation efficiency for gas that is not supported against collapse by turbulence, following [156,157]. FB2 is consistent with a Chabrier IMF, and either the energy or momentum solution is injected depending on which is resolved by the simulations [158].
Figure 10. CGM OVI covering fractions in a sub-L* simulated galaxy at z = 3 for different subgrid models. SF#FB# corresponds to whether the simulations uses star formation model 1 or 2 (SF1 and SF2), and stellar feedback model 1 or 2 (FB1 and FB2). SF1 uses a density threshold for triggering star formation. FB1 corresponds to feedback from a stellar population with a Salpeter IMF and the massive stars exploding with 10 51 erg per 10 M and the energy injected in the simulation matching the Sedov solution [155]. Both of these subgrid models depend on the resolution of the simulations (here 10 pc). SF2 uses a variable star formation efficiency for gas that is not supported against collapse by turbulence, following [156,157]. FB2 is consistent with a Chabrier IMF, and either the energy or momentum solution is injected depending on which is resolved by the simulations [158].
Galaxies 06 00100 g010
Table 1. Overview of the lines discussed at the workshop and referred to in the remainder of this paper.
Table 1. Overview of the lines discussed at the workshop and referred to in the remainder of this paper.
NameTypeWavelength (s) (1)Tracer ofReference for Wavelength (s)
Ly α Recombination 1215 . 67 ÅIonized ISM[56]
C ivCE (2) 1548 . 19 , 1550 . 77 ÅStellar wind, ionized ISM[57]
O iii]CE (2) 1660 . 81 , 1666 . 15 ÅIonized ISM[57]
He iiRecombination 1640 . 42 ÅStellar wind, ionized ISM[57]
[C iii]CE (2) 1906 . 68 ÅISM[57]
C iii]CE (2) 1908 . 73 ÅIonized ISM[57]
H β Recombination 4861 . 36 ÅIonized ISMNIST (3)
[O iii]CE (2) 4958 . 91 , 5006 . 84 ÅIonized ISMNIST (3)
[O i]Recombination 6300 . 30 ÅIonized ISMNIST (3)
H α Recombination 6562 . 80 ÅIonized ISMNIST (3)
[N ii]CE (2) 6548 . 05 , 6583 . 45 ÅIonized ISMNIST (3)
[S ii]CE (2) 6716 . 44 , 6730 . 82 ÅIonized ISMNIST (3)
C iFine-structure 609 . 14 , 370 . 42 μ mAtomic and molecular gasLAMDA (4)
[C ii]Fine-structure 157 . 74 μ mAll ISMLAMDA (4)
[O i]Fine-structure 63 . 18 , 145 . 53 μ mAtomic and molecular gasLAMDA (4)
CORotational 2 . 6 , 1 . 3 , 0 . 87 mm …Molecular gasLAMDA (4)
( 1 ) We give air wavelengths above 2000 Å and below 20,000 Å, and vacuum wavelengths otherwise; ( 2 ) Collisionally excited line; ( 3 ) https://physics.nist.gov/PhysRefData/ASD/lines_form.html; ( 4 ) http://home.strw.leidenuniv.nl/~moldata/.
Table 2. Overview of the computational tools represented at the workshop.
Table 2. Overview of the computational tools represented at the workshop.
NameReferenceDensity Regime1D or 3DSpecies in Chemical Network for Calculating Ionization StatesExact Radiative Transfer Included?
ART2[29], Li 2018 (in prep) 10 2 10 7 cm−33DAtomic database of Cloudy, molecular database of LAMDAyes
Cloudy (1)[13] 10 6 10 15 cm−31D625 species (including atomic ions); CHIANTI, Stout, LAMDA databasesno (9)
DESPOTIC (2)[31]Cool atomic and molecular ISM1DC, O, H, and He, plus a super-species M that represents a composite of metallic elementsno (9)
LIME (3)[26] 10 2 10 12 cm−33DLAMDA databaseyes
MAIHEM (4)[32,67,68]Has been tested at 0 . 4 1200 cm−33D63 species (including atomic ions) (5)no
MOLLIE (6)[25,27] 10 2 10 12 cm−33D18 molecular species (7)yes
RADMC-3D (8)[28]No limits3DLAMDA database, but abundances can also be supplied by the useryes
( 1 ) https://www.nublado.org/; ( 2 ) https://bitbucket.org/krumholz/despotic/; ( 3 ) http://www.nbi.dk/~brinch/index.php?page=lime; ( 4 ) http://maihem.asu.edu/; ( 5 ) H, H + , He, He + , He 2 + , C-C 5 + , N-N 6 + , O-O 7 + , Ne-Ne 9 + , Na-Na 2 + , Mg-Mg 3 + , Si-Si 5 + , S-S 4 + , Ca-Ca 4 + , Fe-Fe 4 + , and electrons; ( 6 ) https://www.cfa.harvard.edu/~eketo/mollie; ( 7 ) CH 3 CN, NH 3 , N 2 H + , H 2 O (para + ortho), CO, CS, 13 CO, H 13 CO + , HCO + , C 17 O, C 18 O, HCN, H 13 CN, HC 15 N, N 2 D + , SiO, H 2 D + ; ( 8 ) http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/; ( 9 ) Escape probability formalism.

Share and Cite

MDPI and ACS Style

Olsen, K.P.; Pallottini, A.; Wofford, A.; Chatzikos, M.; Revalski, M.; Guzmán, F.; Popping, G.; Vázquez-Semadeni, E.; Magdis, G.E.; Richardson, M.L.A.; et al. Challenges and Techniques for Simulating Line Emission. Galaxies 2018, 6, 100. https://doi.org/10.3390/galaxies6040100

AMA Style

Olsen KP, Pallottini A, Wofford A, Chatzikos M, Revalski M, Guzmán F, Popping G, Vázquez-Semadeni E, Magdis GE, Richardson MLA, et al. Challenges and Techniques for Simulating Line Emission. Galaxies. 2018; 6(4):100. https://doi.org/10.3390/galaxies6040100

Chicago/Turabian Style

Olsen, Karen P., Andrea Pallottini, Aida Wofford, Marios Chatzikos, Mitchell Revalski, Francisco Guzmán, Gergö Popping, Enrique Vázquez-Semadeni, Georgios E. Magdis, Mark L. A. Richardson, and et al. 2018. "Challenges and Techniques for Simulating Line Emission" Galaxies 6, no. 4: 100. https://doi.org/10.3390/galaxies6040100

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop