Showing 1–49 of 49 results for author: Bisbas, T

Turbulent Diffuse Molecular Media with Non-ideal Magnetohydrodynamics and Consistent Thermochemistry: Numerical Simulations and Dynamic Characteristics

Authors: Nannan Yue, Lile Wang, Thomas Bisbas, Donghui Quan, Di Li

Abstract: Turbulent diffuse molecular clouds can exhibit complicated morphologies caused by the interactions among radiation, chemistry, fluids, and fields. We performed full 3D simulations for turbulent diffuse molecular interstellar media, featuring time-dependent non-equilibrium thermochemistry co-evolved with magnetohydrodynamics (MHD). Simulation results exhibit the relative abundances of key chemical… ▽ More Turbulent diffuse molecular clouds can exhibit complicated morphologies caused by the interactions among radiation, chemistry, fluids, and fields. We performed full 3D simulations for turbulent diffuse molecular interstellar media, featuring time-dependent non-equilibrium thermochemistry co-evolved with magnetohydrodynamics (MHD). Simulation results exhibit the relative abundances of key chemical species (e.g., C, CO, OH) vary by more than one order of magnitude for the "premature" epoch of chemical evolution ($t\lesssim 2\times 10^5~{\rm yr}$). Various simulations are also conducted to study the impacts of physical parameters. Non-ideal MHD effects are essential in shaping the behavior of gases, and strong magnetic fields ($\sim 10~μ{\rm G}$) tend to inhibit vigorous compressions and thus reduce the fraction of warm gases ($T\gtrsim 10^2~{\rm K}$). Thermodynamical and chemical conditions of the gas are sensitive to modulation by dynamic conditions, especially the energy injection by turbulence. Chemical features, including ionization (cosmic ray and diffuse interstellar radiation), would not directly affect the turbulence power spectra. Nonetheless, their effects are prominent in the distribution profiles of temperatures and gas densities. Comprehensive observations are necessary and useful to eliminate the degeneracies of physical parameters and constrain the properties of diffuse molecular clouds with confidence. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 21 pages, 13 figures; re-submitted the revised manuscript to the Astrophysical Journal

Molecular gas excitation in the circumgalactic medium of MACS1931-26

Authors: L. Ghodsi, J. Zhou, P. Andreani, C. De Breuck, A. W. S. Man, Y. Miyamoto, T. G. Bisbas, A. Lundgren, Z. -Y. Zhang

Abstract: The evolution of galaxies is largely affected by exchanging material with their close environment, the circumgalactic medium (CGM). In this work, we investigate the CGM and the interstellar medium (ISM) of the bright central galaxy (BCG) of the galaxy cluster, MACS1931-26 at z~0.35. We detected [CI](2-1), CO(1-0), and CO(7-6) emission lines with the APEX 12-m and NRO 45-m telescopes. We complement… ▽ More The evolution of galaxies is largely affected by exchanging material with their close environment, the circumgalactic medium (CGM). In this work, we investigate the CGM and the interstellar medium (ISM) of the bright central galaxy (BCG) of the galaxy cluster, MACS1931-26 at z~0.35. We detected [CI](2-1), CO(1-0), and CO(7-6) emission lines with the APEX 12-m and NRO 45-m telescopes. We complemented these single-dish observations with CO(1-0), CO(3-2), and CO(4-3) ALMA interferometric data and inferred the cold molecular hydrogen physical properties. Using a modified large velocity gradient (LVG) model, we modelled the CO and CI emission of the CGM and BCG to extract the gas thermodynamical properties, including the kinetic temperature, the density, and the virialisation factor. Our study shows that the gas in the BCG is highly excited, comparable to the gas in local ultra luminous infrared galaxies (ULIRGs), while the CGM is likely less excited, colder, less dense, and less bound compared to the ISM of the BCG. The molecular hydrogen mass of the whole system derived using [CI](2-1) is larger than the mass derived from CO(1-0) in literature, showing that part of the gas in this system is CO-poor. Additional spatially resolved CI observations in both transitions, CO(1-0) and [CI](2-1), and the completion of the CO SLED with higher CO transitions are crucial to trace the different phases of the gas in such systems and constrain their properties. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 13 pages, 5 figures, accepted for publication by A&A

A case study of gas impacted by black-hole jets with the JWST: outflows, bow shocks, and high excitation of the gas in the galaxy IC5063

Authors: K. M. Dasyra, G. F. Paraschos, F. Combes, P. Patapis, G. Helou, M. Papachristou, J. A. Fernandez-Ontiveros, T. G. Bisbas, L. Spinoglio, L. Armus, M. Malkan

Abstract: We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and se… ▽ More We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and sensitivity in the infrared, now reveal that there are more than ten discrete regions with outflows, nearly doubling the number of such known regions. Outflows exist near the radio lobes, at the nucleus, in a biconical structure perpendicular to the jet, and in a bubble moving against the disk. In some of them, velocities above escape velocity are observed. Stratification is also observed, with higher ionization or excitation gas attaining higher velocities. More outflows and bow shocks, found further away from the nucleus than the radio lobes, in regions without significant radio emission, reveal the existence of past or weak radio jets that interacted with the interstellar medium. The coincidence of the bow shocks with the optical extended emission line region (EELR) suggests that the jets also contributed to the gas ionization. Maps of the H2 gas excitation temperature, T_ex, indicate that the molecular gas is most excited in regions with radio emission. There, T_ex is more than 100K higher than in the EELR interior. We argue that a combination of jet-related shocks and cosmic rays is likely responsible for this excess molecular gas excitation. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: Submitted to ApJ

The CO-dark molecular gas in the cold HI arc

Authors: Gan Luo, Di Li, Zhi-yu Zhang, Thomas G. Bisbas, Ningyu Tang, Lingrui Lin, Yichen Sun, Pei Zuo, Jing Zhou

Abstract: The CO-dark molecular gas (DMG), which refers to the molecular gas not traced by CO emission, is crucial for the evolution of the interstellar medium (ISM). While the gas properties of DMG have been widely explored in the Solar neighborhood, whether or not they are similar in the outer disk regions of the Milky Way is still not well understood. In this Letter, we confirm the existence of DMG towar… ▽ More The CO-dark molecular gas (DMG), which refers to the molecular gas not traced by CO emission, is crucial for the evolution of the interstellar medium (ISM). While the gas properties of DMG have been widely explored in the Solar neighborhood, whether or not they are similar in the outer disk regions of the Milky Way is still not well understood. In this Letter, we confirm the existence of DMG toward a cold HI arc structure at 13 kpc away from the Galactic center with both OH emission and HI narrow self-absorption (HINSA). This is the first detection of HINSA in the outer disk region, in which the HINSA fraction ($N_{\rm HINSA}$/$N_{\rm H_2}$ = 0.022$\pm$0.011) is an order of magnitude higher than the average value observed in nearby evolved dark clouds, but is consistent with that of the early evolutionary stage of dark clouds. The inferred H$_2$ column density from both extinction and OH emission ($N_{\rm H_2} \approx 10^{20}$ cm$^{-2}$) is an order of magnitude higher than previously estimated. Although the ISM environmental parameters are expected to be different between the outer Galactic disk regions and the Solar neighborhood, we find that the visual extinction ($A_{\rm V}$ = 0.19$\pm$0.03 mag), H$_2$-gas density ($n_{\rm H_2} = 91\pm46$ cm$^{-3}$), and molecular fraction (58\%$\pm$28\%) of the DMG are rather similar to those of nearby diffuse molecular clouds. The existence of DMG associated with the expanding HI supershell supports a scenario where the expansion of supershells may trigger the formation of molecular clouds within a crossing timescale of the shock wave ($\sim$10$^6$ yr). △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 9 pages, 5 figures, accepted by A&A Letter

Journal ref: A&A 685, L12 (2024)

$α$-enhanced Astrochemistry: the Carbon cycle in extreme galactic conditions

Authors: Thomas G. Bisbas, Zhi-Yu Zhang, Eda Gjergo, Ying-He Zhao, Gan Luo, Donghui Quan, Xue-Jian Jiang, Yichen Sun, Theodoros Topkaras, Di Li, Ziyi Guo

Abstract: Astrochemistry has been widely developed as a power tool to probe physical properties of the interstellar medium (ISM) in various conditions of the Milky Way (MW) Galaxy, and in near and distant galaxies. Most current studies conventionally apply linear scaling to all elemental abundances based on the gas-phase metallicity. However, these elements, including carbon and oxygen, are enriched differe… ▽ More Astrochemistry has been widely developed as a power tool to probe physical properties of the interstellar medium (ISM) in various conditions of the Milky Way (MW) Galaxy, and in near and distant galaxies. Most current studies conventionally apply linear scaling to all elemental abundances based on the gas-phase metallicity. However, these elements, including carbon and oxygen, are enriched differentially by stellar nucleosynthesis and the overall galactic chemical evolution, evident from $α$-enhancement in multiple galactic observations such as starbursts, high-redshift star-forming galaxies, and low-metallicity dwarfs. We perform astrochemical modeling to simulate the impact of an $α$-enhanced ISM gas cloud on the abundances of the three phases of carbon (C$^+$, C, CO) dubbed as `the carbon cycle'. The ISM environmental parameters considered include two cosmic-ray ionization rates ($ζ_{\rm CR}=10^{-17}$ and $10^{-15}\,{\rm s}^{-1}$), two isotropic FUV radiation field strengths ($χ/χ_0=1$ and $10^2$), and (sub-)linear dust-to-gas relations against metallicity, mimicking the ISM conditions of different galaxy types. In galaxies with [C/O] $<$ 0, CO, C and C$^+$ all decrease in both abundances and emission, though with differential biases. The low-$J$ CO emission is found to be the most stable tracer for the molecular gas, while C and C$^+$ trace H$_2$ gas only under limited conditions, in line with recent discoveries of [CI]-dark galaxies. We call for caution when using [CII]~$158μ$m and [CI](1-0) as alternative H$_2$-gas tracers for both diffuse and dense gas with non-zero [C/O] ratios. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: MNRAS accepted. 21 pages, 15 figures, comments welcome!

SUNRISE: The rich molecular inventory of high-redshift dusty galaxies revealed by broadband spectral line surveys

Authors: Chentao Yang, Alain Omont, Sergio Martín, Thomas G. Bisbas, Pierre Cox, Alexandre Beelen, Eduardo González-Alfonso, Raphaël Gavazzi, Susanne Aalto, Paola Andreani, Cecilia Ceccarelli, Yu Gao, Mark Gorski, Michel Guélin, Hai Fu, R. J. Ivison, Kirsten K. Knudsen, Matthew Lehnert, Hugo Messias, Sebastien Muller, Roberto Neri, Dominik Riechers, Paul van der Werf, Zhi-Yu Zhang

Abstract: Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar… ▽ More Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar at redshift $z=3.911$, and NCv1.143 (NC), a $z=3.565$ starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330-550 GHz. We report the detection of 38 and 25 emission lines in APM and NC, respectively. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the gas excitation analysis, we find that the physical properties and the chemical imprints of the ISM are different between them: the molecular gas is more excited in APM, exhibiting higher molecular-gas temperatures and densities compared to NC; the chemical abundances in APM are akin to the values of local AGN, showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or XDRs, while NC more closely resembles local starburst galaxies. The most significant differences are found in H2O, where the 448GHz H2O line is significantly brighter in APM, likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that at such high column densities, FUV radiation is less important in regulating the ISM, while CRs (X-rays/shocks) are the key players in shaping the abundance of the molecules and the initial conditions of star formation. Such deep spectral line surveys open a new window to study the physical and chemical properties of the ISM and the radiation field of galaxies in the early Universe. (abridged) △ Less

Submitted 22 October, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A. 39 pages (including 9 pages of Appendix). Updated after language editing

Journal ref: A&A 680, A95 (2023)

Abundance ratios of OH/CO and HCO+/CO as probes of the cosmic ray ionization rate in diffuse clouds

Authors: Gan Luo, Zhiyu Zhang, Thomas G. Bisbas, Di Li, Ping Zhou, Ningyu Tang, Junzhi Wang, Pei Zuo, Nannan Yue

Abstract: The cosmic-ray ionization rate (CRIR, $ζ_2$) is one of the key parameters controlling the formation and destruction of various molecules in molecular clouds. However, the current most commonly used CRIR tracers, such as H$_3^+$, OH$^+$, and H$_2$O$^+$, are hard to detect and require the presence of background massive stars for absorption measurements. In this work, we propose an alternative method… ▽ More The cosmic-ray ionization rate (CRIR, $ζ_2$) is one of the key parameters controlling the formation and destruction of various molecules in molecular clouds. However, the current most commonly used CRIR tracers, such as H$_3^+$, OH$^+$, and H$_2$O$^+$, are hard to detect and require the presence of background massive stars for absorption measurements. In this work, we propose an alternative method to infer the CRIR in diffuse clouds using the abundance ratios of OH/CO and HCO$^+$/CO. We have analyzed the response of chemical abundances of CO, OH, and HCO$^+$ on various environmental parameters of the interstellar medium in diffuse clouds and found that their abundances are proportional to $ζ_2$. Our analytic expressions give an excellent calculation of the abundance of OH for $ζ_2$ $\leq$10$^{-15}$ s$^{-1}$, which are potentially useful for modelling chemistry in hydrodynamical simulations. The abundances of OH and HCO$^+$ were found to monotonically decrease with increasing density, while the CO abundance shows the opposite trend. With high-sensitivity absorption transitions of both CO (1--0) and (2--1) lines from ALMA, we have derived the H$_2$ number densities ($n_{\rm H_2}$) toward 4 line-of-sights (LOSs); assuming a kinetic temperature of $T_{\rm k}=50\,{\rm K}$, we find a range of (0.14$\pm$0.03--1.2$\pm$0.1)$\times$10$^2$ cm$^{-3}$}. By comparing the observed and modelled HCO$^+$/CO ratios, we find that $ζ_2$ in our diffuse gas sample is in the { range of $1.0_{-1.0}^{+14.8}$ $\times$10$^{-16}- 2.5_{-2.4}^{+1.4}$ $\times$10$^{-15}$ s$^{-1}$. This is $\sim$2 times higher than the average value measured at higher extinction, supporting an attenuation of CRs as suggested by theoretical models. △ Less

Submitted 6 April, 2023; v1 submitted 24 February, 2023; originally announced February 2023.

Comments: 22 pages, 9 figures, accepted by ApJ

Cold gas disks in main-sequence galaxies at cosmic noon: Low turbulence, flat rotation curves, and disk-halo degeneracy

Authors: Federico Lelli, Zhi-Yu Zhang, Thomas G. Bisbas, Lingrui Lin, Padelis Papadopoulos, James M. Schombert, Enrico Di Teodoro, Antonino Marasco, Stacy S. McGaugh

Abstract: We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation cur… ▽ More We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation curves are flat out to $\sim$8 kpc contrary to previous results pointing to outer declines in the rotation speed $V_{\rm rot}$; (2) The intrinsic velocity dispersions are low ($σ_{\rm CO}\lesssim15$ km/s for CO and $σ_{\rm Hα}\lesssim37$ km/s for H$α$) and imply $V_{\rm rot}/σ_{\rm CO}\gtrsim17-22$ yielding no significant pressure support; (3) Mass models using HST images display a severe disk-halo degeneracy: models with inner baryon dominance and models with "cuspy" dark matter halos can fit the rotation curves equally well due to the uncertainties on stellar and gas masses; (4) Milgromian dynamics (MOND) can successfully fit the rotation curves with the same acceleration scale $a_0$ measured at $z\simeq0$. The question of the amount and distribution of dark matter in high-$z$ galaxies remains unsettled due to the limited spatial extent of the available kinematic data; we discuss the suitability of various emission lines to trace extended rotation curves at high $z$. Nevertheless, the properties of these two high-$z$ galaxies (high $V_{\rm rot}/σ_{\rm V}$ ratios, inner rotation curve shapes, bulge-to-total mass ratios) are remarkably similar to those of massive spirals at $z\simeq0$, suggesting weak dynamical evolution over more than 10 Gyr of the Universe's lifetime. △ Less

Submitted 25 May, 2023; v1 submitted 31 January, 2023; originally announced February 2023.

Comments: 18 pages, 11 figures, 4 tables, 2 appendices. Minor revisions to match the published version in Astronomy and Astrophysics

Journal ref: A&A 672, A106 (2023)

Dependence of Chemical Abundance on the Cosmic Ray Ionization Rate in IC 348

Authors: Gan Luo, Zhi-Yu Zhang, Thomas G. Bisbas, Di Li, Ningyu Tang, Junzhi Wang, Ping Zhou, Pei Zuo, Nannan Yue, Jing Zhou, Lingrui Lin

Abstract: Ions (e.g., H$_3^+$, H$_2$O$^+$) have been used extensively to quantify the cosmic-ray ionization rate (CRIR) in diffuse sightlines. However, measurements of CRIR in low-to-intermediate density gas environments are rare, especially when background stars are absent. In this work, we combine molecular line observations of CO, OH, CH, and HCO$^+$ in the star-forming cloud IC~348, and chemical models… ▽ More Ions (e.g., H$_3^+$, H$_2$O$^+$) have been used extensively to quantify the cosmic-ray ionization rate (CRIR) in diffuse sightlines. However, measurements of CRIR in low-to-intermediate density gas environments are rare, especially when background stars are absent. In this work, we combine molecular line observations of CO, OH, CH, and HCO$^+$ in the star-forming cloud IC~348, and chemical models to constrain the value of CRIR and study the response of the chemical abundances distribution. The cloud boundary is found to have an $A_{\rm V}$ of approximately 4 mag. From the interior to the exterior of the cloud, the observed $^{13}$CO line intensities drop by an order of magnitude. The calculated average abundance of $^{12}$CO (assuming $^{12}$C/$^{13}$C = 65) is (1.2$\pm$0.9) $\times$10$^{-4}$, which increases by a factor of 6 from the interior to the outside regions. The average abundance of CH (3.3$\pm$0.7 $\times$ 10$^{-8}$) is in good agreement with previous findings in diffuse and translucent clouds ($A_{\rm V}$ $<$ 5 mag). However, we did not find a decline in CH abundance in regions of high extinction ($A_{\rm V}\simeq$8 mag) as previously reported in Taurus. By comparing the observed molecular abundances and chemical models, we find a decreasing trend of CRIR as $A_{\rm V}$ increases. The inferred CRIR of $ζ_{cr}$ = (4.7$\pm$1.5) $\times$ 10$^{-16}$ s$^{-1}$ at low $A_{\rm V}$ is consistent with H$^+_3$ measurements toward two nearby massive stars. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 21 pages, 11 figures. Submitted to ApJ

PDFchem: A new fast method to determine ISM properties and infer environmental parameters using probability distributions

Authors: Thomas G. Bisbas, Ewine F. van Dishoeck, Chia-Yu Hu, Andreas Schruba

Abstract: Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a… ▽ More Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a large parameter space. Building on our earlier approach, we present PDFchem, a novel algorithm that models the cold ISM at moderate and large scales using functions connecting the quantities of the local ($A_{\rm V,eff}$) and the observed ($A_{\rm V,obs}$) visual extinctions, and the local number density, $n_{\rm H}$, with probability density functions (PDF) of $A_{\rm V,obs}$ on cloud scales typically tens-to-hundreds of pc as an input. For any given $A_{\rm V,obs}$-PDF, sampled with thousands of clouds, the algorithm instantly computes the average abundances of the most important species (HI, H$_2$, CII, CI, CO, OH, OH$^+$, H$_2$O$^+$, CH, HCO$^+$) and performs radiative transfer calculations to estimate the average emission of the most commonly observed lines ([CII]~$158μ$m, both [CI] fine-structure lines and the first five rotational transitions of $^{12}$CO). We examine two $A_{\rm V,obs}$-PDFs corresponding to a non star-forming and a star-forming ISM region, under a variety of environmental parameters combinations. These cover FUV intensities in the range of $χ/χ_0=10^{-1}-10^3$, cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}-10^{-13}\,{\rm s}^{-1}$ and metallicities in the range of $Z=0.1-2\,{\rm Z}_{\odot}$. PDFchem is fast, easy to use, reproduces the PDR quantities of the time-consuming hydrodynamical models and can be used directly with observed data to understand the evolution of the cold ISM chemistry. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 33 pages, 27 figures. Accepted in MNRAS. The algorithm can be found in: https://github.com/tbisbas/PDFchem. Comments welcome!

Cosmic-ray-induced H$_2$ line emission: Astrochemical modeling and implications for JWST observations

Authors: Brandt A. L. Gaches, Shmuel Bialy, Thomas G. Bisbas, Marco Padovani, Daniel Seifried, Stefanie Walch

Abstract: Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the inte… ▽ More Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the interstellar medium (ISM). Aims: We aim to investigate the impact of certain assumptions regarding the H$_2$ chemical models and ISM density distributions on the emission of cosmic-ray induced H$_2$ emission lines. This is of particular importance for utilizing observations of these lines with the James Webb Space Telescope to constrain the cosmic-ray ionization rate. Methods: We compare the predicted emission from cosmic-ray induced, ro-vibrationally excited H$_2$ emission lines for different one- and three-dimensional models with varying assumptions on the gas chemistry and density distribution. Results: We find that the model predictions of the H$_2$ line intensities for the (1-0)S(0), (1-0)Q(2), (1-0)O(2) and (1-0)O(4) transitions at 2.22, 2.41, 2.63 and 3.00 $μ$m, respectively, are relatively independent of the astro-chemical model and the gas density distribution when compared against the H$_2$ column density, making them robust tracer of the cosmic-ray ionization rate. Conclusions: We recommend the use of ro-vibrational H$_2$ line emission in combination with estimation of the cloud's H$_2$ column density, to constrain the ionization rate and the spectrum of low energy cosmic-rays. △ Less

Submitted 6 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted by A&A

Journal ref: A&A 664, A150 (2022)

The origin of the [CII]-deficit in a simulated dwarf galaxies starburst

Authors: Thomas G. Bisbas, Stefanie Walch, Thorsten Naab, Natalia Lahén, Rodrigo Herrera-Camus, Ulrich P. Steinwandel, Constantina M. Fotopoulou, Chia-Yu Hu, Peter H. Johansson

Abstract: We present [CII] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star-formation rate by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding HII regions and supernova (SN) explosions. For galaxies with properties similar… ▽ More We present [CII] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star-formation rate by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding HII regions and supernova (SN) explosions. For galaxies with properties similar to the modelled ones, we find that the [CII] emission remains optically thin throughout the merging process. We identify the Warm Neutral Medium ($3<\log T_{\rm gas}<4$ with $χ_{\rm HI}>2χ_{\rm H2}$) to be the primary source of [CII] emission ($\sim58\%$ contribution), although at stages when the HII regions are young and dense (during star cluster formation or SNe in the form of ionized bubbles) they can contribute $\gtrsim50\%$ to the total [CII] emission. We find that the [CII]/FIR ratio decreases due to thermal saturation of the [CII] emission caused by strong FUV radiation fields emitted by the massive star clusters, leading to a [CII]-deficit medium. We investigate the [CII]-SFR relation and find an approximately linear correlation which agrees well with observations, particularly those from the Dwarf Galaxy Survey. Our simulation reproduces the observed trends of [CII]/FIR versus $Σ_{\rm SFR}$ and $Σ_{\rm FIR}$, and it agrees well with the Kennicutt relation of SFR-FIR luminosity. We propose that local peaks of [CII] in resolved observations may provide evidence for ongoing massive cluster formation. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: To appear in ApJ. 24 pages, 17 figures. Comments welcome!

Linking pressure gradients with the stability of molecular clouds in galactic outflows

Authors: K. M. Dasyra, G. F. Paraschos, T. Bisbas, F. Combes, J. A. Fernandez-Ontiveros

Abstract: The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impa… ▽ More The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impact molecular clouds. Using a radiative transfer code that self-consistently performs astrochemical and thermal balance calculations based on the available gas heating sources, we found that mechanical heating and cosmic ray (CR) heating are fully capable of individually reproducing the data. In our best-fit model, CRs provide about 1/3rd of the dense gas heating at the radio lobes, emphasizing the role of this often neglected mechanism in heating the gas and potentially generating outflows. The gas temperature and density indicate that the jet passage leads to an increase of about 1 order of magnitude in the internal pressure Pi of molecular clouds (with Pi/k from 8*10^5 up to 7*10^6 K cm^-3), irrespective of the excitation mechanism. From the fluxes of [S II] and [N II] lines in VLT MUSE data, the external pressure Pe of molecular clouds increases in several regions enough to exceed Pi. This result leads us to conclude that we are observing the expansion of an ionized overpressurized cocoon that compresses molecular clouds and that could lead to their collapse. Some jet-impacted clouds, nonetheless, near pathways that the jet cleared have increased Pi and decreased Pe. They are likely to undergo evaporation of their outer layers. Part of the evaporated layers could mass load the outflow thanks to ram pressure from co-spatial ionized gas flows. The observed pressure changes thus suggest that both star formation enhancement and suppression could simultaneously occur. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: Submitted to Nature Astronomy

The impact of cosmic-ray attenuation on the carbon cycle emission in molecular clouds

Authors: Brandt A. L. Gaches, Thomas G. Bisbas, Shmuel Bialy

Abstract: Observations of the emission of the carbon cycle species (C, C+ CO) are commonly used to diagnose gas properties in the interstellar medium but are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, $ζ$, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds.… ▽ More Observations of the emission of the carbon cycle species (C, C+ CO) are commonly used to diagnose gas properties in the interstellar medium but are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, $ζ$, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds. However, astrochemical models commonly assume a constant cosmic-ray ionization rate in the clouds. We investigate the effect of cosmic-ray attenuation on the emission of carbon cycle species from molecular clouds, in particular the [CII] 158 $μ$m, [CI] 609 $μ$m and CO (J = 1 - 0) 115.27 GHz lines. We use a post-processed chemical model of diffuse and dense simulated molecular clouds and quantify the variation in both column densities and velocity integrated line emission of the carbon cycle with different cosmic-ray ionization rate models. We find that the abundances and column densities of carbon cycle species is significantly impacted by the chosen cosmic-ray ionization rate model: no single constant ionization rate can reproduce the abundances modelled with an attenuated cosmic-ray model. Further, we show that constant ionization rate models fail to simultaneously reproduce the integrated emission of the lines we consider, and their deviations from a physically derived cosmic-ray attenuation model is too complex to be simply corrected. We demonstrate that the two clouds we model exhibit a similar average $A_{\rm V, eff}$ -- $n_{\rm H}$ relationship, resulting in an average relation between the cosmic-ray ionization rate and density $ζ(n_{\rm H})$. We conclude by providing a number of implementation recommendations for CRs in astrochemical models, but emphasize the necessity for column-dependent cosmic-ray ionization rate prescriptions. △ Less

Submitted 11 January, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted to A&A. Major results in Figures 4, 6, 10, 13

Journal ref: A&A 658, A151 (2022)

Photodissociation Region Diagnostics Across Galactic Environments

Authors: Thomas G. Bisbas, Jonathan C. Tan, Kei E. I. Tanaka

Abstract: We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}$-$10^{-14}\,{\rm s}^{-1}$, far-UV intensities in the range of $G_0=1$-$10^3$ and metallicities in the range of $Z=0.1$-… ▽ More We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}$-$10^{-14}\,{\rm s}^{-1}$, far-UV intensities in the range of $G_0=1$-$10^3$ and metallicities in the range of $Z=0.1$-$2\,{\rm Z}_{\odot}$. The simulations also probe a range of densities and levels of turbulence, including cases where the gas has undergone recent compression due to cloud-cloud collisions. We examine: i) the column densities of carbon species across the cycle of CII, CI and CO, along with OI, in relation to the HI-to-H$_2$ transition; ii) the velocity-integrated emission of [CII]~$158μ$m, [$^{13}$CII]~$158μ$m, [CI]~$609μ$m and $370μ$m, [OI]~$63μ$m and $146μ$m, and of the first ten $^{12}$CO rotational transitions; iii) the corresponding Spectral Line Energy Distributions; iv) the usage of [CII] and [OI]~$63μ$m to describe the dynamical state of the clouds; v) the behavior of the most commonly used ratios between transitions of CO and [CI]; and vi) the conversion factors for using CO and CI as H$_2$-gas tracers. We find that enhanced cosmic-ray energy densities enhance all aforementioned line intensities. At low metallicities, the emission of [CII] is well connected with the H$_2$ column, making it a promising new H$_2$ tracer in metal-poor environments. The conversion factors of $X_{\rm CO}$ and $X_{\rm CI}$ depend on metallicity and the cosmic-ray ionization rate, but not on FUV intensity. In the era of ALMA, SOFIA and the forthcoming CCAT-prime telescope, our results can be used to understand better the behaviour of systems in a wide range of galactic and extragalactic environments. △ Less

Submitted 12 December, 2020; originally announced December 2020.

Comments: 33 pages, 24 figures, 3 tables. Submitted to MNRAS. Comments welcome!

Star Cluster Formation in Orion A

Authors: Wanggi Lim, Fumitaka Nakamura, Benjamin Wu, Thomas G. Bisbas, Jonathan C. Tan, Edward Chambers, John Bally, Shuo Kong, Peregrine McGehee, Dariusz C. Lis, Volker Ossenkopf-Okada, Álvaro Sánchez-Monge

Abstract: We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $Σ$ scatter of older YSOs can be sig… ▽ More We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $Σ$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$μ$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A. △ Less

Submitted 25 June, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: 16 pages, 11 figures, 4 tables, published in PASJ

The Astrochemical Impact of Cosmic Rays in Protoclusters II: CI-to-H$_2$ and CO-to-H$_2$ Conversion Factors

Authors: Brandt A. L. Gaches, Stella S. R. Offner, Thomas G. Bisbas

Abstract: We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have… ▽ More We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have $X_{\rm CO} = (2.5 \pm 1)\times10^{20}$ cm$^{-2}$(K km s$^{-1}$)$^{-1}$, consistent with Milky Way observations. We find that changing the cosmic ray ionization rate from external sources from the canonical $ζ\approx 10^{-17}$ to $ζ\approx 10^{-16}$ s$^{-1}$, which better represents observations in diffuse gas, reduces $X_{\rm CO}$ by 0.2 dex for clusters with surface densities below 3 g cm$^{-2}$. We show that embedded sources regulate $X_{\rm CO}$ and decrease its variance across a wide range of surface densities and star formation efficiencies. Our models reproduce the trends of a decreased $X_{\rm CO}$ in extreme cosmic ray environments. $X_{\rm CI}$ has been proposed as an alternative to $X_{\rm CO}$ due to its brightness at high redshifts. The inclusion of internal cosmic ray sources leads to 1.2 dex dispersion in $X_{\rm CI}$ ranging from $2\times10^{20} < X_{\rm CI} < 4\times10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We show that $X_{\rm CI}$ is highly sensitive to the underlying cosmic ray model. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 9 pages, 3 figures. ApJ Accepted

SILCC-Zoom: H$_2$ and CO-dark gas in molecular clouds -- The impact of feedback and magnetic fields

Authors: D. Seifried, S. Haid, S. Walch, E. M. A. Borchert, T. G. Bisbas

Abstract: We analyse the CO-dark molecular gas content of simulated molecular clouds from the SILCC-Zoom project. The simulations reach a resolution of 0.1 pc and include H$_2$ and CO formation, radiative stellar feedback and magnetic fields. CO-dark gas is found in regions with local visual extinctions $A_\text{V, 3D} \sim$ 0.2 - 1.5, number densities of 10 - 10$^{3}$ cm$^{-3}$ and gas temperatures of few… ▽ More We analyse the CO-dark molecular gas content of simulated molecular clouds from the SILCC-Zoom project. The simulations reach a resolution of 0.1 pc and include H$_2$ and CO formation, radiative stellar feedback and magnetic fields. CO-dark gas is found in regions with local visual extinctions $A_\text{V, 3D} \sim$ 0.2 - 1.5, number densities of 10 - 10$^{3}$ cm$^{-3}$ and gas temperatures of few 10 K - 100 K. CO-bright gas is found at number densities above 300 cm$^{-3}$ and temperatures below 50 K. The CO-dark gas fractions range from 40% to 95% and scale inversely with the amount of well-shielded gas ($A_\text{V, 3D}$ $\gtrsim$ 1.5), which is smaller in magnetised molecular clouds. We show that the density, chemical abundances and $A_\text{V, 3D}$ along a given line-of-sight cannot be properly determined from projected quantities. As an example, pixels with a projected visual extinction of $A_\text{V, 2D} \simeq$ 2.5 - 5 can be both, CO-bright or CO-dark, which can be attributed to the presence or absence of strong density enhancements along the line-of-sight. By producing synthetic CO(1-0) emission maps of the simulations with RADMC-3D, we show that about 15 - 65\% of the H$_2$ is in regions with intensities below the detection limit. Our clouds have $X_\text{CO}$-factors around 1.5 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ with a spread of up to a factor $\sim$ 4, implying a similar uncertainty in the derived total H$_2$ masses and even worse for individual pixels. Based on our results, we suggest a new approach to determine the H$_2$ mass, which relies on the availability of CO(1-0) emission and $A_\text{V, 2D}$ maps. It reduces the uncertainty of the clouds' overall H$_2$ mass to a factor of $\lesssim$ 1.8 and for individual pixels, i.e. on sub-pc scales, to a factor of $\lesssim$ 3. △ Less

Submitted 17 December, 2019; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: 20 pages, 20 figures, accepted for publication in MNRAS

The Astrochemical Impact of Cosmic Rays in Protoclusters I: Molecular Cloud Chemistry

Authors: Brandt A. L. Gaches, Stella S. R. Offner, Thomas G. Bisbas

Abstract: We present astrochemical photo-dissociation region models in which cosmic ray attenuation has been fully coupled to the chemical evolution of the gas. We model the astrochemical impact of cosmic rays, including those accelerated by protostellar accretion shocks, on molecular clouds hosting protoclusters. Our models with embedded protostars reproduce observed ionization rates. We study the imprint… ▽ More We present astrochemical photo-dissociation region models in which cosmic ray attenuation has been fully coupled to the chemical evolution of the gas. We model the astrochemical impact of cosmic rays, including those accelerated by protostellar accretion shocks, on molecular clouds hosting protoclusters. Our models with embedded protostars reproduce observed ionization rates. We study the imprint of cosmic ray attenuation on ions for models with different surface cosmic ray spectra and different star formation efficiencies. We find that abundances, particularly ions, are sensitive to the treatment of cosmic rays. We show the column densities of ions are under predicted by the `classic' treatment of cosmic rays by an order of magnitude. We also test two common chemistry approximations used to infer ionization rates. We conclude that the approximation based on the H$_3^+$ abundance under predicts the ionization rate except in regions where the cosmic rays dominate the chemistry. Our models suggest the chemistry in dense gas will be significantly impacted by the increased ionization rates, leading to a reduction in molecules such as NH$_3$ and causing H$_2$-rich gas to become [C II] bright. △ Less

Submitted 9 May, 2019; v1 submitted 6 May, 2019; originally announced May 2019.

Comments: Accepted to ApJ. The code used is public at https://uclchem.github.io/3dpdr.html

Simulating the atomic and molecular content of molecular clouds using probability distributions of physical parameters

Authors: Thomas G. Bisbas, Andreas Schruba, Ewine F. van Dishoeck

Abstract: Modern observations of the interstellar medium (ISM) in galaxies detect a variety of atomic and molecular species. The goal is to connect these observations to the astrochemical properties of the ISM. 3D hydro-chemical simulations attempt this but due to extreme computational cost, they have to rely on simplified chemical networks and are bound to individual case studies. We present an alternative… ▽ More Modern observations of the interstellar medium (ISM) in galaxies detect a variety of atomic and molecular species. The goal is to connect these observations to the astrochemical properties of the ISM. 3D hydro-chemical simulations attempt this but due to extreme computational cost, they have to rely on simplified chemical networks and are bound to individual case studies. We present an alternative approach which models the ISM at larger scales by an ensemble of pre-calculated 1D thermo-chemical photodissociation region (PDR) calculations that determine the abundance and excitation of atomic and molecular species. We adopt lognormal distributions of column density (\avpdf s) for which each column density is linked to a volume density as derived by hydrodynamical simulations. We consider two lognormal {\avpdf}s: a diffuse, low density medium with average visual extinction of $\overline{{\rm A}_V}=0.75\,{\rm mag}$ and dispersion of $σ=0.5$ and a denser giant molecular cloud with $\overline{{\rm A}_V}=4\,{\rm mag}$ and $σ=0.8$. We treat the UV radiation field, cosmic-ray ionization rate and metallicity as free parameters. We find that the low density medium remains fully HI- and CII-dominated under all explored conditions. The denser cloud remains almost always molecular (i.e. H$_2$-dominated) while its carbon phase (CO, CI and CII) is sensitive to the above free parameters, implying that existing methods of tracing H$_2$-rich gas may require adjustments depending on environment. Our numerical framework can be used to estimate the PDR properties of large ISM regions and quantify trends with different environmental parameters as it is fast, covers wide parameter space, and is flexible for extensions. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: 16 pages, 11 figures, resubmitted to MNRAS, referee's comments have been accounted for

The Interstellar Medium Properties of Heavily Reddened Quasars & Companions at z~2.5 with ALMA & JVLA

Authors: Manda Banerji, Gareth Jones, Jeff Wagg, Chris Carilli, Thomas Bisbas, Paul Hewett

Abstract: We study the interstellar medium (ISM) properties of three heavily reddened quasars at $z\sim2.5$ as well as three millimetre-bright companion galaxies near these quasars. New JVLA and ALMA observations constrain the CO(1-0), CO(7-6) and [CI]$^3$P$_2-^3$P$_1$ line emission as well as the far infrared to radio continuum. The gas excitation and physical properties of the ISM are constrained by compa… ▽ More We study the interstellar medium (ISM) properties of three heavily reddened quasars at $z\sim2.5$ as well as three millimetre-bright companion galaxies near these quasars. New JVLA and ALMA observations constrain the CO(1-0), CO(7-6) and [CI]$^3$P$_2-^3$P$_1$ line emission as well as the far infrared to radio continuum. The gas excitation and physical properties of the ISM are constrained by comparing our observations to photo-dissociation region (PDR) models. The ISM in our high-redshift quasars is composed of very high-density, high-temperature gas which is already highly enriched in elements like carbon. One of our quasar hosts is shown to be a close-separation ($<$2 arcsec) major merger with different line emission properties in the millimeter-bright galaxy and quasar components. Low angular resolution observations of high-redshift quasars used to assess quasar excitation properties should therefore be interpreted with caution as they could potentially be averaging over multiple components with different ISM conditions. Our quasars and their companion galaxies show a range of CO excitation properties spanning the full extent from starburst-like to quasar-like spectral line energy distributions. We compare gas masses based on CO, CI and dust emission, and find that these can disagree when standard assumptions are made regarding the values of $α_{\rm{CO}}$, the gas-to-dust ratio and the atomic carbon abundances. We conclude that the ISM properties of our quasars and their companion galaxies are diverse and likely vary spatially across the full extent of these complex, merging systems. △ Less

Submitted 31 May, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: Accepted for publication in MNRAS (18 pages, 10 figures and 5 tables)

The classical D-type expansion of spherical H II regions

Authors: Robin J. R. Williams, Thomas G. Bisbas, Thomas J. Haworth, Jonathan Mackey

Abstract: Recent numerical and analytic work has highlighted some shortcomings in our understanding of the dynamics of H II region expansion, especially at late times, when the H II region approaches pressure equilibrium with the ambient medium. Here we reconsider the idealized case of a constant radiation source in a uniform and spherically symmetric ambient medium, with an isothermal equation of state. A… ▽ More Recent numerical and analytic work has highlighted some shortcomings in our understanding of the dynamics of H II region expansion, especially at late times, when the H II region approaches pressure equilibrium with the ambient medium. Here we reconsider the idealized case of a constant radiation source in a uniform and spherically symmetric ambient medium, with an isothermal equation of state. A thick-shell solution is developed which captures the stalling of the ionization front and the decay of the leading shock to a weak compression wave as it escapes to large radii. An acoustic approximation is introduced to capture the late-time damped oscillations of the H II region about the stagnation radius. Putting these together, a matched asymptotic equation is derived for the radius of the ionization front which accounts for both the inertia of the expanding shell and the finite temperature of the ambient medium. The solution to this equation is shown to agree very well with the numerical solution at all times, and is superior to all previously published solutions. The matched asymptotic solution can also accurately model the variation of H II region radius for a time-varying radiation source. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: 8 pages, 5 figures. Authors' original version, revised version with minor changes accepted for publication in MNRAS

New places and phases of CO-poor/CI-rich molecular gas in the Universe

Authors: Padelis P. Papadopoulos, Thomas G. Bisbas, Zhiyu Zhang

Abstract: In this work we extend the work on the recently discovered role of Cosmic Rays (CRs) in regulating the average CO/$\rm H_2$ abundance ratio in molecular clouds (and thus their CO line visibility) in starburst galaxies, and find that it can lead to a CO-poor/CI-rich $\rm H_2 $ gas phase even in environments with Galactic or in only modestly enhanced CR backgrounds expected in ordinary star-forming… ▽ More In this work we extend the work on the recently discovered role of Cosmic Rays (CRs) in regulating the average CO/$\rm H_2$ abundance ratio in molecular clouds (and thus their CO line visibility) in starburst galaxies, and find that it can lead to a CO-poor/CI-rich $\rm H_2 $ gas phase even in environments with Galactic or in only modestly enhanced CR backgrounds expected in ordinary star-forming galaxies. Furthermore, the same CR-driven astro-chemistry raises the possibility of a widespread phase transition of molecular gas towards a CO-poor/CI-rich phase in: a) molecular gas outflows found in star-forming galaxies, b) active galactic nuclei (AGNs), and c) near synchrotron-emitting radio jets and the radio-loud cores of powerful radio galaxies. For main sequence galaxies we find that CRs can render some of their molecular gas mass CO-invisible, compounding the effects of low metallicities. Imaging the two fine structure lines of atomic carbon with resolution high enough to search beyond the CI/CO-bright line regions associated with central starbursts can reveal such a CO-poor/CI-rich molecular gas phase, provided that relative brightness sensitivity levels of $T_b$(CI $1-0$)/$T_b$(CO $J=1-0$)$\sim $0.15 are reached. The capability to search for such gas in the Galaxy is now at hand with the new high-frequency survey telescope HEAT deployed in Antarctica and future ones to be deployed in Dome A. ALMA can search for such gas in star-forming spiral disks, galactic molecular gas outflows and the CR-intense galactic and circumgalactic gas-rich environments of radio-loud objects. △ Less

Submitted 26 April, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: 11 pages, 5 figures, MNRAS accepted

The Inception of Star Cluster Formation Revealed by [CII] Emission Around an Infrared Dark Cloud

Authors: Thomas G. Bisbas, Jonathan C. Tan, Timea Csengeri, Benjamin Wu, Wanggi Lim, Paola Caselli, Rolf Guesten, Oliver Ricken, Denise Riquelme

Abstract: We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC… ▽ More We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation. △ Less

Submitted 7 March, 2018; v1 submitted 1 March, 2018; originally announced March 2018.

Comments: 5 pages, 4 figures, submitted to MNRAS Letters

Synthetic observations of star formation and the interstellar medium

Authors: Thomas J. Haworth, Simon C. O. Glover, Christine M. Koepferl, Thomas G. Bisbas, James E. Dale

Abstract: Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium… ▽ More Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium. We also discuss some possible directions for future research using synthetic observations. △ Less

Submitted 9 November, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 71 Pages, 16 figures. Revised following reviewers comments

GMC Collisions as Triggers of Star Formation. V. Observational Signatures

Authors: Thomas G. Bisbas, Kei E. I. Tanaka, Jonathan C. Tan, Benjamin Wu, Fumitaka Nakamura

Abstract: We present calculations of molecular, atomic and ionic line emission from simulations of giant molecular cloud (GMC) collisions. We post-process snapshots of the magneto-hydrodynamical simulations presented in an earlier paper in this series by Wu et al. (2017) of colliding and non-colliding GMCs. Using photodissociation region (PDR) chemistry and radiative transfer we calculate the level populati… ▽ More We present calculations of molecular, atomic and ionic line emission from simulations of giant molecular cloud (GMC) collisions. We post-process snapshots of the magneto-hydrodynamical simulations presented in an earlier paper in this series by Wu et al. (2017) of colliding and non-colliding GMCs. Using photodissociation region (PDR) chemistry and radiative transfer we calculate the level populations and emission properties of $^{12}$CO $J=1-0$, [CI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ at $609\,μ$m, [CII] $158\,μ$m and [OI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ transition at $63\,μ$m. From integrated intensity emission maps and position-velocity diagrams, we find that fine-structure lines, particularly the [CII] $158\,μ$m, can be used as a diagnostic tracer for cloud-cloud collision activity. These results hold even in more evolved systems in which the collision signature in molecular lines has been diminished. △ Less

Submitted 17 October, 2017; v1 submitted 21 June, 2017; originally announced June 2017.

Comments: 10 pages, 7 figures, accepted for publication in ApJ, comments welcome

The Interstellar Medium and Star Formation of Galactic Disks. I. ISM and GMC properties with Diffuse FUV and Cosmic Ray Backgrounds

Authors: Qi Li, Jonathan C. Tan, Duncan Christie, Thomas G. Bisbas, Benjamin Wu

Abstract: We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gra… ▽ More We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gradient designed to mimic the Milky Way. The effects of cosmic rays, including radial gradients in their heating and ionization rates, are also explored. The final simulations in this series achieve $4\:$pc resolution across the $\sim20\:$kpc global disk diameter, with heating and cooling followed down to temperatures of $\sim10\:$K. The disks are evolved for $300\:$Myr, which is enough time for the ISM to achieve a quasi-statistical equilibrium. In particular, the mass fraction of molecular gas stabilizes by $\sim$200 Myr. Additional global ISM properties are analysed. Giant molecular clouds (GMCs) are also identified and the statistical properties of their populations examined. GMCs are tracked as the disks evolve. GMC collisions, which may be a means of triggering star cluster formation, are counted and the rates compared with analytic models. Relatively frequent GMC collision rates are seen in these simulations and their implications for understanding GMC properties, including the driving of internal turbulence, are discussed. △ Less

Submitted 13 November, 2017; v1 submitted 12 June, 2017; originally announced June 2017.

Comments: Accepted by PASJ (cloud-cloud collision special issue)

Cosmic-ray induced destruction of CO in star-forming galaxies

Authors: Thomas G. Bisbas, Ewine F. van Dishoeck, Padelis P. Papadopoulos, László Szücs, Shmuel Bialy, Zhi-Yu Zhang

Abstract: We explore the effects of the expected higher cosmic ray (CR) ionization rates $ζ_{\rm CR}$ on the abundances of carbon monoxide (CO), atomic carbon (C), and ionized carbon (C$^+$) in the H$_2$ clouds of star-forming galaxies. The study of Bisbas et al. (2015) is expanded by: a) using realistic inhomogeneous Giant Molecular Cloud (GMC) structures, b) a detailed chemical analysis behind the CR-indu… ▽ More We explore the effects of the expected higher cosmic ray (CR) ionization rates $ζ_{\rm CR}$ on the abundances of carbon monoxide (CO), atomic carbon (C), and ionized carbon (C$^+$) in the H$_2$ clouds of star-forming galaxies. The study of Bisbas et al. (2015) is expanded by: a) using realistic inhomogeneous Giant Molecular Cloud (GMC) structures, b) a detailed chemical analysis behind the CR-induced destruction of CO, and c) exploring the thermal state of CR-irradiated molecular gas. CRs permeating the interstellar medium with $ζ_{\rm CR}$$\gtrsim 10\times$(Galactic) are found to significantly reduce the [CO]/[H$_2$] abundance ratios throughout the mass of a GMC. CO rotational line imaging will then show much clumpier structures than the actual ones. For $ζ_{\rm CR}$$\gtrsim 100\times$(Galactic) this bias becomes severe, limiting the utility of CO lines for recovering structural and dynamical characteristics of H$_2$-rich galaxies throughout the Universe, including many of the so-called Main Sequence (MS) galaxies where the bulk of cosmic star formation occurs. Both C$^+$ and C abundances increase with rising $ζ_{\rm CR}$, with C remaining the most abundant of the two throughout H$_2$ clouds, when $ζ_{\rm CR}\sim (1-100)\times$(Galactic). C$^+$ starts to dominate for $ζ_{\rm CR}$$\gtrsim 10^3\times$(Galactic). The thermal state of the gas in the inner and denser regions of GMCs is invariant with $T_{\rm gas}\sim 10\,{\rm K}$ for $ζ_{\rm CR}\sim (1-10)\times$(Galactic). For $ζ_{\rm CR}$$\sim 10^3\times$(Galactic) this is no longer the case and $T_{\rm gas}\sim 30-50\,{\rm K}$ are reached. Finally we identify OH as the key species whose $T_{\rm gas}-$sensitive abundance could mitigate the destruction of CO at high temperatures. △ Less

Submitted 24 March, 2017; originally announced March 2017.

Comments: 17 pages, 12 figures, accepted by ApJ

ALMA observations of atomic carbon in z~4 dusty star-forming galaxies

Authors: M. S. Bothwell, J. E. Aguirre, M. Aravena, M. Bethermin, T. G. Bisbas, S. C. Chapman, C. De Breuck, A. H. Gonzalez, T. R. Greve, Y. Hezaveh, J. Ma, M. Malkan, D. P. Marrone, E. J. Murphy, J. S. Spilker, M. Strandet, J. D. Vieira, A. Weiss

Abstract: We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift… ▽ More We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift range $2 < z < 5$. Using the [CI] line as a tracer of the molecular ISM, we find a mean molecular gas mass for SPT-DSFGs of $6.6 \times 10^{10}$ M$_{\odot}$. This is in tension with gas masses derived via low-$J$ $^{12}$CO and dust masses; bringing the estimates into accordance requires either (a) an elevated CO-to-H$_2$ conversion factor for our sample of $α_{\rm CO} \sim 2.5$ and a gas-to-dust ratio $\sim200$, or (b) an high carbon abundance $X_{\rm CI} \sim 7\times10^{-5}$. Using observations of a range of additional atomic and molecular lines (including [CI], [CII], and multiple transitions of CO), we use a modern Photodissociation Region code (3D-PDR) to assess the physical conditions (including the density, UV radiation field strength, and gas temperature) within the ISM of the DSFGs in our sample. We find that the ISM within our DSFGs is characterised by dense gas permeated by strong UV fields. We note that previous efforts to characterise PDR regions in DSFGs may have significantly underestimated the density of the ISM. Combined, our analysis suggests that the ISM of extreme dusty starbursts at high redshift consists of dense, carbon-rich gas not directly comparable to the ISM of starbursts in the local Universe. △ Less

Submitted 13 December, 2016; originally announced December 2016.

Comments: 21 pages, 12 figures. Accepted for publication in MNRAS

Photochemical-dynamical models of externally FUV irradiated protoplanetary discs

Authors: Thomas J. Haworth, Douglas Boubert, Stefano Facchini, Thomas G. Bisbas, Cathie J. Clarke

Abstract: There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simul… ▽ More There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their "critical radius". This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes. △ Less

Submitted 7 September, 2016; originally announced September 2016.

Comments: 14 Pages. Accepted for publication in MNRAS

Radiative transfer meets Bayesian statistics: where does a galaxy's [CII] emission come from?

Authors: Gioacchino Accurso, Amélie Saintonge, Thomas. G. Bisbas, Serena Viti

Abstract: The [CII] 158$μ$m emission line can arise in all phases of the ISM, therefore being able to disentangle the different contributions is an important yet unresolved problem when undertaking galaxy-wide, integrated [CII] observations. We present a new multi-phase 3D radiative transfer interface that couples Starburst99, a stellar spectrophotometric code, with the photoionisation and astrochemistry co… ▽ More The [CII] 158$μ$m emission line can arise in all phases of the ISM, therefore being able to disentangle the different contributions is an important yet unresolved problem when undertaking galaxy-wide, integrated [CII] observations. We present a new multi-phase 3D radiative transfer interface that couples Starburst99, a stellar spectrophotometric code, with the photoionisation and astrochemistry codes Mocassin and 3D-PDR. We model entire star forming regions, including the ionised, atomic and molecular phases of the ISM, and apply a Bayesian inference methodology to parametrise how the fraction of the [CII] emission originating from molecular regions, $f_{[CII],mol}$, varies as a function of typical integrated properties of galaxies in the local Universe. The main parameters responsible for the variations of $f_{[CII],mol}$ are specific star formation rate (sSFR), gas phase metallicity, HII region electron number density ($n_e$), and dust mass fraction. For example, $f_{[CII],mol}$ can increase from 60% to 80% when either $n_e$ increases from 10$^{1.5}$ to 10$^{2.5}$cm$^{-3}$, or SSFR decreases from $10^{-9.6}$ to $10^{-10.6}$ yr$^{-1}$. Our model predicts for the Milky Way that $f_{[CII],mol}$$=75.8\pm5.9$%, in agreement with the measured value of 75%. When applying the new prescription to a complete sample of galaxies from the Herschel Reference Survey (HRS), we find that anywhere from 60 to 80% of the total integrated [CII] emission arises from molecular regions. △ Less

Submitted 5 October, 2016; v1 submitted 12 July, 2016; originally announced July 2016.

Comments: 17 pages including appendices, 14 figures, 3 tables, Accepted for publication in MNRAS

External photoevaporation of protoplanetary discs in sparse stellar groups: the impact of dust growth

Authors: Stefano Facchini, Cathie J. Clarke, Thomas G. Bisbas

Abstract: We estimate the mass loss rates of photoevaporative winds launched from the outer edge of protoplanetary discs impinged by an ambient radiation field. We focus on mild/moderate environments (the number of stars in the group/cluster is N ~ 50), and explore disc sizes ranging between 20 and 250 AU. We evaluate the steady-state structures of the photoevaporative winds by coupling temperature estimate… ▽ More We estimate the mass loss rates of photoevaporative winds launched from the outer edge of protoplanetary discs impinged by an ambient radiation field. We focus on mild/moderate environments (the number of stars in the group/cluster is N ~ 50), and explore disc sizes ranging between 20 and 250 AU. We evaluate the steady-state structures of the photoevaporative winds by coupling temperature estimates obtained with a PDR code with 1D radial hydrodynamical equations. We also consider the impact of dust dragging and grain growth on the final mass loss rates. We find that these winds are much more significant than have been appreciated hitherto when grain growth is included in the modelling: in particular, mass loss rates > 1e-8 M_sun/yr are predicted even for modest background field strengths ( ~ 30 G_0) in the case of discs that extend to R > 150 AU. Grain growth significantly affects the final mass loss rates by reducing the average cross section at FUV wavelengths, and thus allowing a much more vigorous flow. The radial profiles of observable quantities (in particular surface density, temperature and velocity patterns) indicate that these winds have characteristic features that are now potentially observable with ALMA. In particular, such discs should have extended gaseous emission that is dust depleted in the outer regions, characterised by a non-Keplerian rotation curve, and with a radially increasing temperature gradient. △ Less

Submitted 27 January, 2016; originally announced January 2016.

Comments: 18 pages, 13 figures and 1 table. Accepted for publication in MNRAS

TORUS-3DPDR: A self-consistent code treating three-dimensional photoionization and photodissociation regions

Authors: T. G. Bisbas, T. J. Haworth, M. J. Barlow, S. Viti, T. J. Harries, T. Bell, J. A. Yates

Abstract: The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium eith… ▽ More The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations. △ Less

Submitted 16 September, 2015; originally announced September 2015.

Comments: 17 pages, 10 Figures, accepted for publication in MNRAS. Comments welcome

On the relative importance of different microphysics on the D-type expansion of galactic HII regions

Authors: Thomas J. Haworth, Tim J. Harries, David M. Acreman, Thomas G. Bisbas

Abstract: Radiation hydrodynamics (RHD) simulations are used to study many astrophysical phenomena, however they require the use of simplified radiation transport and thermal prescriptions to reduce computational cost. In this paper we present a systematic study of the importance of microphysical processes in RHD simulations using the example of D-type HII region expansion. We compare the simplest hydrogen-… ▽ More Radiation hydrodynamics (RHD) simulations are used to study many astrophysical phenomena, however they require the use of simplified radiation transport and thermal prescriptions to reduce computational cost. In this paper we present a systematic study of the importance of microphysical processes in RHD simulations using the example of D-type HII region expansion. We compare the simplest hydrogen-only models with those that include: ionisation of H, He, C, N, O, S and Ne, different gas metallicity, non-LTE metal line blanketed stellar spectral models of varying metallicity, radiation pressure, dust and treatment of photodissociation regions. Each of these processes are explicitly treated using modern numerical methods rather than parameterisation. In line with expectations, changes due to microphysics in either the effective number of ionising photons or the thermal structure of the gas lead to differences in D-type expansion. In general we find that more realistic calculations lead to the onset of D-type expansion at smaller radii and a slower subsequent expansion. Simulations of star forming regions using simplified microphysics are therefore likely overestimating the strength of radiative feedback. We find that both variations in gas metallicity and the inclusion of dust can affect the ionisation front evolution at the 10-20 per cent level over 500kyr, which could substantially modify the results of simplified 3D models including feedback. Stellar metallicity, radiation pressure and the inclusion of photodissociation regions are all less significant effects at the 1 per cent level or less, rendering them of minor importance in the modelling the dynamical evolution of HII regions. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 14 pages, 9 figures, accepted for publication in the MNRAS main journal

StarBench: The D-type expansion of an HII region

Authors: T. G. Bisbas, T. J. Haworth, R. J. R. Williams, J. Mackey, P. Tremblin, A. C. Raga, S. J. Arthur, C. Baczynski, J. E. Dale, T. Frostholm, S. Geen, T. Haugboelle, D. Hubber, I. T. Iliev, R. Kuiper, J. Rosdahl, D. Sullivan, S. Walch, R. Wuensch

Abstract: StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase o… ▽ More StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: 20 pages, 12 Figures, 4 Tables. Accepted for publication in MNRAS. Comments are welcome. Participation in future StarBench tests is also welcome

Effective destruction of CO by cosmic rays: implications for tracing H$_2$ gas in the Universe

Authors: Thomas G. Bisbas, Padelis P. Papadopoulos, Serena Viti

Abstract: We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems thro… ▽ More We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests. △ Less

Submitted 19 February, 2015; v1 submitted 14 February, 2015; originally announced February 2015.

Comments: Accepted for publication in ApJ, 29 pages

Astrochemical Correlations in Molecular Clouds

Authors: Brandt A. L. Gaches, Stella S. R. Offner, Erik W. Rosolowsky, Thomas G. Bisbas

Abstract: We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spati… ▽ More We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: Accepted by ApJ. 19 pages, 13 figures

A Photodissociation Region study of NGC 4038

Authors: T. G. Bisbas, T. A. Bell, S. Viti, M. J. Barlow, J. Yates, M. Vasta

Abstract: We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than thre… ▽ More We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\leχ\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609μm and [OI] 146μm fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609μm/[OI] 146μm is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value. △ Less

Submitted 6 June, 2014; originally announced June 2014.

Comments: 12 pages, 6 figures, 3 tables, to appear in MNRAS

An Alternative Accurate Tracer of Molecular Clouds: The "$X_{\rm CI}$-Factor"

Authors: Stella S. R. Offner, Thomas G. Bisbas, Tom A. Bell, Serena Viti

Abstract: We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $μ$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compar… ▽ More We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $μ$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compare the emission of CI and CO for model clouds irradiated by 1 and 10 times the average background and demonstrate that CI is a comparable or superior tracer of the molecular gas distribution for column densities up to $6 \times 10^{23}$ cm$^{-2}$. Our results hold for both reduced and full chemical networks. For our fiducial Galactic cloud we derive an average $X_{\rm CO}$ of $3.0\times 10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$s and $X_{\rm CI}$ of $1.1\times 10^{21}$ cm$^{-2}$K$^{-1}$km$^{-1}$s. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: 5 pages, 4 figures, 1 table, accepted to MNRAS Letters

Clumps and triggered star formation in ionised molecular clouds

Authors: S. Walch, A. P. Whitworth, T. G. Bisbas, R. Wunsch, D. A. Hubber

Abstract: Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is cha… ▽ More Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is characterised by its fractal dimension, D (with 2.0 <= D <= 2.8), and the variance of its (log-normal) density distribution, sigma_0^2 (with 0.36 <= sigma_0^2 <= 1.42). Our study focuses on the morphology of the swept-up cold gas and the distribution and statistics of the resulting star formation. If the fractal dimension is low, the border of the HII region is dominated by extended shell-like structures, and these break up into a small number of massive high-density clumps which then spawn star clusters; star formation occurs relatively quickly, and delivers somewhat higher stellar masses. Conversely, if the fractal dimension is high, the border of the HII region is dominated by a large number of pillars and cometary globules, which contain compact dense clumps and tend to spawn single stars or individual multiple systems; star formation occurs later, the stellar masses are somewhat lower, and the stars are more widely distributed. △ Less

Submitted 18 June, 2013; originally announced June 2013.

Comments: 12 pages, 11 figures, accepted for publication by MNRAS

Modeling the Atomic-to-Molecular Transition and Chemical Distributions of Turbulent Star-Forming Clouds

Authors: Stella S. R. Offner, Thomas G. Bisbas, Serena Viti, Thomas A. Bell

Abstract: We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions a… ▽ More We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches. △ Less

Submitted 25 April, 2013; originally announced April 2013.

Comments: 18 pages, 18 Figures, accepted to ApJ

3D-PDR: A new three-dimensional astrochemistry code for treating Photodissociation Regions

Authors: T. G. Bisbas, T. A. Bell, S. Viti, J. Yates, M. J. Barlow

Abstract: Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a… ▽ More Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPix-based ray-tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, ii) a uniform-density spherical cloud interacting with a two-component external radiation field, and iii) a cometary globule interacting with a plane-parallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment. △ Less

Submitted 4 December, 2012; v1 submitted 5 September, 2012; originally announced September 2012.

Comments: 19 pages, 11 figures, accepted for publication in MNRAS

A nebula in your computer: simulating the physics and chemistry of an HII region

Authors: Thomas G. Bisbas

Abstract: In this contribution we discuss about numerical modeling of nebulae. In particular we emphasize on the dynamical evolution of an HII region and the chemical structure of a Photodissociation region. We do this using the Smoothed Particles Hydrodynamics code SEREN and the recently developed astrochemistry code 3D-PDR, respectively. We show an example application by simulating a cometary globule usin… ▽ More In this contribution we discuss about numerical modeling of nebulae. In particular we emphasize on the dynamical evolution of an HII region and the chemical structure of a Photodissociation region. We do this using the Smoothed Particles Hydrodynamics code SEREN and the recently developed astrochemistry code 3D-PDR, respectively. We show an example application by simulating a cometary globule using these two codes. △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: 4 pages, 1 figure; to appear in The Labyrinth of Star Formation, (eds.) D. Stamatellos, S. Goodwin, and D. Ward-Thompson, Springer, in press

Dispersal of molecular clouds by ionising radiation

Authors: S. Walch, A. P. Whitworth, T. Bisbas, R. Wunsch, D. Hubber

Abstract: The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-no… ▽ More The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-normal density PDF. (i) As regards star formation, in the short term ionising feedback is positive, in the sense that star formation occurs much more quickly in gas that is compressed by the high pressure of the ionised gas. However, in the long term ionising feedback is negative, in the sense that most of the cloud is dispersed with an outflow rate of up to ~0.01M_sun/yr, on a timescale comparable with the sound-crossing time for the ionised gas (~1-2Myr), and triggered star formation is therefore limited to a few percent of the cloud's mass. (ii) As regards the morphology of the ionisation fronts (IFs) bounding the HII region and the systematics of outflowing gas, we distinguish two regimes. For low D<=2.2, the initial cloud is dominated by large-scale structures, so the neutral gas tends to be swept up into a few extended coherent shells, and the ionised gas blows out through a few large holes between these shells; we term these HII regions "shell-dominated". Conversely, for high D>=2.6, the initial cloud is dominated by small-scale structures, and these are quickly overrun by the advancing IF, thereby producing neutral pillars whilst the ionised gas blows out through a large number of small holes between the pillars; we term these HII regions "pillar-dominated". (iii) As regards the injection of bulk kinetic energy, by ~1Myr, the expansion of the HII region has delivered a rms velocity of ~6km/s; this represents less than 0.1% of the total energy radiated by the O7 star. △ Less

Submitted 27 June, 2012; originally announced June 2012.

Comments: 13 pages, 8 figures, 2 tables; submitted to MNRAS

Comparing simulations of ionisation triggered star formation and observations in RCW 120

Authors: S. Walch, A. Whitworth, T. Bisbas, D. A. Hubber, R. Wuensch

Abstract: Massive clumps within the swept-up shells of bubbles, like that surrounding the galactic HII region RCW 120, have been interpreted in terms of the Collect and Collapse (C&C) mechanism for triggered star formation. The cold, dusty clumps surrounding RCW 120 are arranged in an almost spherical shell and harbour many young stellar objects. By performing high-resolution, three-dimensional SPH simulati… ▽ More Massive clumps within the swept-up shells of bubbles, like that surrounding the galactic HII region RCW 120, have been interpreted in terms of the Collect and Collapse (C&C) mechanism for triggered star formation. The cold, dusty clumps surrounding RCW 120 are arranged in an almost spherical shell and harbour many young stellar objects. By performing high-resolution, three-dimensional SPH simulations of HII regions expanding into fractal molecular clouds, we investigate whether the formation of massive clumps in dense, swept-up shells necessarily requires the C&C mechanism. In a second step, we use RADMC-3D to compute the synthetic dust continuum emission from our simulations, in order to compare them with observations of RCW 120 made with APEX-LABOCA at 870 micron. We show that a distribution of clumps similar to the one seen in RCW 120 can readily be explained by a non-uniform initial molecular cloud structure. Hence, a shell-like configuration of massive clumps does not imply that the C&C mechanism is at work. Rather, we find a hybrid form of triggering, which combines elements of C&C and Radiatively Driven Implosion (RDI). In addition, we investigate the reliability of deriving clump masses from their 870 micron emission. We find that for clumps with more than 100 M_sun the observational estimates are accurate to within a factor of two and that, even at these long wavelengths, it is important to account for the radiative heating from triggered, embedded protostars. △ Less

Submitted 21 May, 2015; v1 submitted 15 September, 2011; originally announced September 2011.

Comments: 10 pages, 7 figures, submitted to MNRAS

Radiation driven implosion and triggered star formation

Authors: Thomas G. Bisbas, Richard Wunsch, Anthony P. Whitworth, David A. Hubber, Stefanie Walch

Abstract: We present simulations of initially stable isothermal clouds exposed to ionising radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN and an HEALPix-based photo-ionisation algorithm to simulate the propagation of the ionising radiation and the resulting dynamical ev… ▽ More We present simulations of initially stable isothermal clouds exposed to ionising radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN and an HEALPix-based photo-ionisation algorithm to simulate the propagation of the ionising radiation and the resulting dynamical evolution of the cloud. We find that the incident ionising flux, $Φ_{_{\rm LyC}}$, is the critical parameter determining the cloud evolution. At moderate fluxes, a large fraction of the cloud mass is converted into stars. As the flux is increased, the fraction of the cloud mass that is converted into stars and the mean masses of the individual stars both decrease. Very high fluxes simply disperse the cloud. Newly-formed stars tend to be concentrated along the central axis of the cloud (i.e. the axis pointing in the direction of the incident flux). For given cloud parameters, the time, $t_{_\star}$, at which star formation starts is proportional to $Φ_{_{\rm LyC}}^{-1/3}$. The pattern of star formation found in the simulations is similar to that observed in bright-rimmed clouds. △ Less

Submitted 18 May, 2011; originally announced May 2011.

Comments: 36 pages, 13 figures, accepted in ApJ

Radiation Driven Implosion and Triggered Star Formation

Authors: T. G. Bisbas, A. P. Whitworth, R. Wünsch, D. A. Hubber, S. Walch

Abstract: We present simulations of initially stable isothermal clouds exposed to ionizing radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas et al. (2009). We find that the incident io… ▽ More We present simulations of initially stable isothermal clouds exposed to ionizing radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas et al. (2009). We find that the incident ionizing flux is the critical parameter determining the evolution: high fluxes simply disperse the cloud, whereas low fluxes trigger star formation. We find a clear connection between the intensity of the incident flux and the parameters of star formation. △ Less

Submitted 16 July, 2010; originally announced July 2010.

Comments: 4 pages, 2 figures, conference proceedings, IAU Symposium 270 (eds. Alves, Elmegreen, Girart, Trimble)

Smoothed Particle Hydrodynamics simulations of expanding HII regions. I. Numerical methods and tests

Authors: T. G. Bisbas, R. Wunsch, A. P. Whitworth, D. A. Hubber

Abstract: We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of H… ▽ More We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of HEALPix to make the algorithm adaptive, so that fine angular resolution is invoked only where it is needed, and the computational cost is kept low. We present simulations of (i) the spherically symmetric expansion of an HII region inside a uniform-density, non--self-gravitating cloud; (ii) the spherically symmetric expansion of an HII region inside a uniform-density, self-gravitating cloud; (iii) the expansion of an off-centre HII region inside a uniform-density, non--self-gravitating cloud, resulting in rocket acceleration and dispersal of the cloud; and (iv) radiatively driven compression and ablation of a core overrun by an HII region. The new algorithm provides the means to explore and evaluate the role of ionizing radiation in regulating the efficiency and statistics of star formation. △ Less

Submitted 16 February, 2009; v1 submitted 13 February, 2009; originally announced February 2009.

Comments: 12 pages, 16 figures, simulation movies available at http://galaxy.ig.cas.cz/~richard/HIIregion/

Radiative transfer and the energy equation in SPH simulations of star formation

Authors: D. Stamatellos, A. P. Whitworth, T. Bisbas, S. Goodwin

Abstract: We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particle's heating and cooling. The method captures -- at minimal computational… ▽ More We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particle's heating and cooling. The method captures -- at minimal computational cost -- the effects of (i) the rotational and vibrational degrees of freedom of H2, H2 dissociation, H0 ionisation, (ii) opacity changes due to ice mantle melting, sublimation of dust, molecular lines, H-, bound-free and free-free processes and electron scattering; (iv) external irradiation; and (v) thermal inertia. The new algorithm reproduces the results of previous authors and/or known analytic solutions. The computational cost is comparable to a standard SPH simulation with a simple barotropic equation of state. The method is easy to implement, can be applied to both particle- and grid-based codes, and handles optical depths 0<tau<10^{11}. △ Less

Submitted 2 August, 2007; v1 submitted 1 May, 2007; originally announced May 2007.

Comments: Submitted to A&A, recommended for publication