BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas G. Bisbas

Thomas G. Bisbas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.EP

Catalog footprint

What is connected

15works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photodissociation Region Diagnostics Across Galactic Environments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Star Cluster Formation in Orion A

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photochemical-dynamical models of externally FUV irradiated protoplanetary discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Astrochemical Correlations in Molecular Clouds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation driven implosion and triggered star formation

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.

0 citations0 reviewsNo code linkedOpen access