Source author record

Shmuel Bialy

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

12 published item(s)

preprint2022arXiv

Constraining the cosmic-ray ionization rate and their spectrum with NIR spectroscopy of dense clouds -- A test-bed for JWST

Low-energy cosmic-rays (CRs) control the thermo-chemical state and the coupling between gas and magnetic fields in dense molecular clouds, the sites of star-formation. However, current estimates of the low-energy CR spectrum ($E \lesssim 1$ GeV) and the associated CR ionization rate are highly uncertain. We apply, for the first time, a new method for constraining the CR ionization rate and the CR spectral shape using H$_2$ rovibrational lines from cold molecular clouds. Using the MMIRS instrument on the MMT, we obtained deep near-infrared (NIR) spectra in six positions within four dense cores, G150, G157, G163, G198, with column densities $N_{\rm H_2} \approx 10^{22}$ cm$^{-2}$. We derive 3$σ$ upper limits on the H$_2$ $(1-0)$S(0) line (2.22 $μ$m) brightness in the range $I = 5.9 \times 10^{-8}$ to $1.2 \times 10^{-7}$ erg cm$^{-2}$ s$^{-1}$ sr$^{-1}$ for the different targets. Using both an analytic model and a numerical model of CR propagation, we convert these into upper limits on the CR ionization rate in the clouds' interior, $ζ= 1.5$ to $3.6 \times 10^{-16}$ s$^{-1}$, and lower limits on the low-energy spectral slope of interstellar CR protons, $α= -0.97$ to $-0.79$. We show that while MMT was unable to detect the H$_2$ lines due to high atmospheric noise, JWST/NIRSpec will be able to efficiently detect the CR-excited H$_2$ lines, making it the ideal method for constraining the otherwise elusive low-energy CRs, shedding light on the sources and propagation modes of CRs.

preprint2022arXiv

Cosmic rays in molecular clouds probed by H$_{2}$ rovibrational lines -- Perspectives for the James Webb Space Telescope

Cosmic rays (CRs) at sub-TeV energies play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H$_{2}$. Their characterisation is important both for the interpretation of observations and for the development of theoretical models. The methods used so far for estimating the CR ionisation rate ($ζ$) in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. We refine and extend the method proposed by Bialy (2020) to estimate $ζ$ by observing rovibrational transitions of H$_{2}$ at near-infrared wavelengths, which are mainly excited by secondary CR electrons. Combining models of interstellar CR propagation and attenuation with the calculation of the expected secondary electron spectrum and updated H$_{2}$ excitation cross sections by electron collisions, we derive the intensity of the four H$_{2}$ rovibrational transitions observable in dense, cold gas: (1-0)O(2), (1-0)Q(2), (1-0)S(0), and (1-0)O(4). The proposed method allows the estimation of $ζ$ for a given observed line intensity and H$_{2}$ column density. We are also able to deduce the shape of the low-energy CR proton spectrum impinging upon the molecular cloud. We present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar CR proton spectrum. We comment on the capability of the James Webb Space Telescope to detect these near-infrared H$_{2}$ lines, making it possible to derive for the first time spatial variation of $ζ$ in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will be possible to test competing models of CR propagation and attenuation in the interstellar medium, as well as compare CR spectra in different Galactic regions.

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.

preprint2022arXiv

Star formation near the Sun is driven by expansion of the Local Bubble

For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust. However, the precise shape and extent of this shell, the impetus and timescale for its formation, and its relationship to nearby star formation have remained uncertain, largely due to low-resolution models of the local interstellar medium. Leveraging new spatial and dynamical constraints from the Gaia space mission, here we report an analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a scenario where the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's center beginning 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

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.

preprint2020arXiv

Study Turbulence and Probe Magnetic Field Using Gradients Technique: Application to HI-to-H2 Transition Regions

The atomic-to-molecular (HI-to-H$_2$) transition in photodissociation regions (PDRs) has been investigated over the last several decades through analytic and numerical modeling. However, classical PDR models typically assume uniform density gas, ignoring the turbulent nature of the interstellar medium. Recently, Bialy et al. (2017b, 2019) have presented a theoretical framework for studying the HI-to-H$_2$ in a realistic turbulent medium with a non-homogeneous density structure. Here we extend these turbulent-chemical models to explore the possibility of tracing the magnetic field direction in turbulent PDRs using the Gradients Technique. We utilize both subsonic and supersonic magnetohydrodynamic numerical simulations for chemical HI/H$_2$ balance calculations. We confirm that the density fluctuations induced by turbulence can disperse the distribution of H$_2$ and HI fraction. We find that the energy spectrum of moment maps gets shallower when the sonic Mach number MS increases. We explore the ability in magnetic field tracing of gradients of higher-order velocity centroids and compare their performance with that of traditional velocity centroid gradients (VCGs) and with intensity gradients (IGs). We find that the velocity gradients of the second-order centroids (VC$_2$Gs) are more accurate than VCGs and IGs in probing the magnetic field orientation.

preprint2020arXiv

The driving scale -- density decorrelation scale relation in a turbulent medium

Density fluctuations produced by supersonic turbulence are of great importance to astrophysical chemical models. A property of these density fluctuations is that the two point correlation function decreases with increasing scale separation. The relation between the density decorrelation length scale ($L_{\rm dec}$) and the turbulence driving scale ($L_{\rm drive}$) determines how turbulence affects the density and chemical structures in the interstellar medium (ISM), and is a key component for using observations of atomic and molecular tracers to constrain turbulence properties. We run a set of numerical simulations of supersonic magnetohydrodynamic turbulence, with different sonic Mach numbers ($\mathcal{M}_s=4.5, 7$) , and driven on varying scales (1/2.5, 1/5, 1/7) the box length. We derive the $L_{\rm dec}-L_{\rm drive}$ relation as a function of Mach number, driving scale, and the orientation of the line-of-sight (LOS) in respect to the magnetic-field. We find that the mean ratio $L_{\rm dec}/L_{\rm drive} = 0.19 \pm 0.10$, when averaged over snapshots, Mach numbers, driving lengths, and the three LOSs. For LOS parallel to the magnetic field the density structures are statistically smaller and the $L_{\rm dec}-L_{\rm drive}$ relation is tighter, with $L_{\rm dec}/L_{\rm drive} = 0.112 \pm 0.024$. We discuss our results in the context of using observations of chemical tracers to constrain the dominant turbulence driving scale in the ISM.

preprint2016arXiv

Analytic HI-to-H2 Photodissociation Transition Profiles

We present a simple analytic procedure for generating atomic-to-molecular (HI-to-H$_2$) density profiles for optically thick clouds illuminated by far-ultraviolet radiation. Our procedure is based on the analytic theory for the structure of one-dimensional HI/H$_2$ photon-dominated regions, presented by Sternberg et al. (2014). Depth-dependent HI and H$_2$ density fractions may be computed for arbitrary gas density, far-ultraviolet field intensity, and the metallicity dependent H$_2$ formation rate coefficient, and dust absorption cross section. We use our procedure to generate a set of HI-to-H$_2$ transition profiles for a wide range of conditions, from the weak- to strong-field limits, and from super-solar down to low metallicities. We show that if presented as functions of dust optical depth the HI and H$_2$ density profiles depend primarily on the Sternberg "$αG$ parameter" (dimensionless) that determines the dust optical depth associated with the total photodissociated HI column. We derive a universal analytic formula for the HI-to-H$_2$ transition points as a function of just $αG$. Our formula will be useful for interpreting emission-line observations of HI/H$_2$ interfaces, for estimating star-formation thresholds, and for sub-grid components in hydrodynamics simulations.

preprint2015arXiv

CO/H2, C/CO, OH/CO, and OH/O2 in Dense Interstellar Gas: From High Ionization to Low Metallicity

We present numerical computations and analytic scaling relations for interstellar ion-molecule gas phase chemistry down to very low metallicities ($ 10^{-3} \times$ solar), and/or up to high driving ionization rates. Relevant environments include the cool interstellar medium (ISM) in low-metallicity dwarf galaxies, early enriched clouds at the reionization and Pop-II star formation era, and in dense cold gas exposed to intense X-ray or cosmic-ray sources. We focus on the behavior for H$_2$, CO, CH, OH, H$_2$O and O$_2$, at gas temperatures $\sim 100$ K, characteristic of a cooled ISM at low metallicities. We consider shielded or partially shielded one-zone gas parcels, and solve the gas phase chemical rate equations for the steady-state "metal-molecule" abundances for a wide range of ionization parameters, $ζ/n$, and metallicties, $Z'$. We find that the OH abundances are always maximal near the H-to-H$_2$ conversion points, and that large OH abundances persist at very low metallicities even when the hydrogen is predominantly atomic. We study the OH/O$_2$, C/CO and OH/CO abundance ratios, from large to small, as functions of $ζ/n$ and $Z'$. Much of the cold dense ISM for the Pop-II generation may have been OH-dominated and atomic rather than CO-dominated and molecular.

preprint2015arXiv

HI-to-H2 transitions in the Perseus Molecular Cloud

We apply the Sternberg et al. (2014) theoretical model to analyze HI and H2 observations in the Perseus molecular cloud. We constrain the physical properties of the HI shielding envelopes and the nature of the HI-to-H2 transitions. Our analysis (Bialy et al. 2015) implies that in addition to cold neutral gas (CNM), less dense thermally-unstable gas (UNM) significantly contributes to the shielding of the H2 cores in Perseus.

preprint2015arXiv

HI-to-H2 Transitions in the Perseus Molecular Cloud

We use the Sternberg et al. (2014) theory for interstellar atomic to molecular (HI-to-H$_2$) conversion to analyze HI-to-H$_2$ transitions in five (low-mass) star-forming and dark regions in the Perseus molecular cloud, B1, B1E, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M$_{\odot}$ pc$^{-2}$ are consistent with HI-to-H$_2$ transitions dominated by HI-dust shielding in predominantly atomic envelopes. For each source, we constrain the dimensionless parameter $αG$, and the ratio $I_{\rm UV}/n$, of the FUV intensity to hydrogen gas density. We find $αG$ values from 5.0 to 26.1, implying characteristic atomic hydrogen densities 11.8 to 1.8 cm$^{-3}$, for $I_{\rm UV} \approx 1$ appropriate for Perseus. Our analysis implies that the dusty HI shielding layers are probably multiphased, with thermally unstable UNM gas in addition to cold CNM within the 21 cm kinematic radius.