Researcher profile

Shmuel Bialy

Shmuel Bialy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 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.