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.

Low energy cosmic-rays (CRs) are responsible for gas heating and ionization of interstellar clouds, which in turn introduces coupling to Galactic magnetic fields. So far the CR ionization rate (CRIR) has been estimated using indirect methods, such as its effect on the abundances of various rare molecular species. Here we show that the CRIR may be constrained from line emission of H$_2$ rovibrational transitions, excited by CRs. We derive the required conditions for CRs to dominate line excitation, and show that CR-excited lines may be detected with the Very Large Telescope (VLT) over 8 hours integration. Our method, if successfully applied to a variety of clouds at different Galactic locations will provide improved constraints on the spectrum of low energy CRs and their origins.

Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters.

We study the thermal structure of the neutral atomic (H {\small I}) interstellar medium across a wide range of metallicities, from supersolar down to vanishing metallicity, and for varying UV intensities and cosmic-ray ionization rates. We calculate self-consistently the gas temperature and species abundances (with a special focus on the residual H$_2$), assuming thermal and chemical steady-state. For solar metallicity, $Z' \equiv 1$, we recover the known result that there exists a pressure range over which the gas is multiphased, with the warm ($\sim 10^4$ K, WNM) and cold ($\sim 100$ K, CNM) phases coexisting at the same pressure. At a metallicity $Z' \approx 0.1$, the CNM is colder (compared to $Z'=1$) due to the reduced efficiency of photoelectric heating. For $Z' \lesssim 0.1$, cosmic-ray ionization becomes the dominant heating mechanism and the WNM-to-CNM transition shifts to ever increasing pressure/density as the metallicity is reduced. For metallicities $Z' \lesssim 0.01$, H$_2$ cooling becomes important, lowering the temperature of the WNM (down to $\approx 600$ K), and smoothing out the multiphase phenomenon. At vanishing metallicities, H$_2$ heating becomes effective and the multiphase phenomenon disappears entirely. We derive analytic expressions for the critical densities for the warm-to-cold phase transition in the different regimes, and the critical metallicities for H$_2$ cooling and heating. We discuss potential implications on the star-formation rates of galaxies and self-regulation theories.

Active Galactic Nuclei (AGNs) emit substantial fluxes of high-energy electromagnetic radiation, and have therefore attracted some recent attention for their negative impact on galactic habitability. In this paper, we propose that AGNs may also engender the following beneficial effects: (i) prebiotic synthesis of biomolecular building blocks mediated by ultraviolet (UV) radiation, and (ii) powering photosynthesis on certain free-floating planets and moons. We also reassess the harmful biological impact of UV radiation originating from AGNs, and find that their significance could have been overestimated. Our calculations suggest that neither the positive nor negative ramifications stemming from a hypothetical AGN in the Milky Way are likely to affect putative biospheres in most of our Galaxy. On the other hand, we find that a sizable fraction of all planetary systems in galaxies with either disproportionately massive black holes ($\sim 10^{9-10}\,M_\odot$) or high stellar densities (e.g., compact dwarf galaxies) might be susceptible to both the beneficial and detrimental consequences of AGNs, with the former potentially encompassing a greater spatial extent than the latter.

`Oumuamua (1I/2017 U1) is the first object of interstellar origin observed in the Solar System. Recently, \citet{Micheli2018} reported that `Oumuamua showed deviations from a Keplerian orbit at a high statistical significance. The observed trajectory is best explained by an excess radial acceleration $\Delta a \propto r^{-2}$, where $r$ is the distance of `Oumuamua from the Sun. Such an acceleration is naturally expected for comets, driven by the evaporating material. However, recent observational and theoretical studies imply that `Oumuamua is not an active comet. We explore the possibility that the excess acceleration results from Solar radiation pressure. The required mass-to-area ratio is $(m/A)\approx 0.1$ g cm$^{-2}$. For a thin sheet this requires a thickness of $\approx 0.3-0.9$ mm. We find that although extremely thin, such an object would survive an interstellar travel over Galactic distances of $\sim 5$ kpc, withstanding collisions with gas and dust-grains as well as stresses from rotation and tidal forces. We discuss the possible origins of such an object. Our general results apply to any light probes designed for interstellar travel.

We present the detection analysis of a diffuse molecular cloud at z$_{abs}$=2.4636 towards the quasar SDSS J1513+0352(z$_{em}\,\simeq$ 2.68) observed with the X-shooter spectrograph(VLT). We measure very high column densities of atomic and molecular hydrogen, with log N(HI,H$_2$)$\simeq$21.8,21.3. This is the highest H$_2$ column density ever measured in an intervening damped Lyman-alpha system but we do not detect CO, implying log N(CO)/N(H$_2$) < -7.8, which could be due to a low metallicity of the cloud. From the metal absorption lines, we derive the metallicity to be Z $\simeq$ 0.15 Z$_{\odot}$ and determine the amount of dust by measuring the induced extinction of the background quasar light, A$_V$ $\simeq$ 0.4. We also detect Ly-$\alpha$ emission at the same redshift, with a centroid located at a most probable impact parameter of only $\rho\,\simeq$ 1.4 kpc. We argue that the line of sight is therefore likely passing through the ISM of a galaxy as opposed to the CGM. The relation between the surface density of gas and that of star formation seems to follow the global empirical relation derived in the nearby Universe although our constraints on the star formation rate and on the galaxy extent remain too loose to be conclusive. We study the transition from atomic to molecular hydrogen using a theoretical description based on the microphysics of molecular hydrogen. We use the derived chemical properties of the cloud and physical conditions (T$_k\,\simeq$90 K and n$\simeq$250 cm$^{-3}$ derived through the excitation of H$_2$ rotational levels and neutral carbon fine structure transitions to constrain the fundamental parameters that govern this transition. By comparing the theoretical and observed HI column densities, we are able to bring an independent constraint on the incident UV flux, which we find to be in agreement with that estimated from the observed star formation rate.

We investigate the metallicity dependence of HI surface densities in star-forming regions along many lines of sight within 70 nearby galaxies, probing kpc to 50 pc scales. We employ HI, SFR, stellar mass, and metallicity (gradient) measurements from the literature, spanning a wide range (5 dex) in stellar and gas mass and (1.6 dex) in metallicity. We consider metallicities as observed, or rescaled to match the mass-metallicity relation determined for SDSS galaxies. At intermediate to high metallicities (0.3-2 times solar), we find that the HI surface densities saturate at sufficiently large total gas surface density. The maximal HI columns vary approximately inversely with metallicity, and show little variation with spatial resolution, galactocentric radius, or among galaxies. In the central parts of massive spiral galaxies the HI gas is depressed by factors of 2. The observed behavior is naturally reproduced by metallicity dependent shielding theories for the HI-to-H2 transitions in star-forming galaxies. We show that the inverse scaling of the maximal HI columns with metallicity suggests that the area filling fraction of atomic-molecular complexes in galaxies is of order unity, and weakly dependent on metallicity.

We study the effect of density fluctuations induced by turbulence on the HI/H$_2$ structure in photodissociation regions (PDRs) both analytically and numerically. We perform magnetohydrodynamic numerical simulations for both subsonic and supersonic turbulent gas, and chemical HI/H$_2$ balance calculations. We derive atomic-to-molecular density profiles and the HI column density probability density function (PDF) assuming chemical equilibrium. We find that while the HI/H$_2$ density profiles are strongly perturbed in turbulent gas, the mean HI column density is well approximated by the uniform-density analytic formula of Sternberg et al. (2014). The PDF width depends on (a) the radiation intensity to mean density ratio, (b) the sonic Mach number and (c) the turbulence decorrelation scale, or driving scale. We derive an analytic model for the HI PDF and demonstrate how our model, combined with 21 cm observations, can be used to constrain the Mach number and driving scale of turbulent gas. As an example, we apply our model to observations of HI in the Perseus molecular cloud. We show that a narrow observed HI PDF may imply small scale decorrelation, pointing to the potential importance of subcloud-scale turbulence driving.

We explore the effects of the expected higher cosmic ray (CR) ionization rates $\zeta_{\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 $\zeta_{\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 $\zeta_{\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 $\zeta_{\rm CR}$, with C remaining the most abundant of the two throughout H$_2$ clouds, when $\zeta_{\rm CR}\sim (1-100)\times$(Galactic). C$^+$ starts to dominate for $\zeta_{\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 $\zeta_{\rm CR}\sim (1-10)\times$(Galactic). For $\zeta_{\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.

The process of atomic-to-molecular (HI-to-H$_2$) gas conversion is fundamental for molecular-cloud formation and star formation. 21 cm observations of the star-forming region W43 revealed extremely high HI column densities, of 120-180 M$_{\odot}$ pc$^{-2}$, a factor of 10-20 larger than predicted by HI-to-H$_2$ transition theories. We analyze the observed HI with an HI-to-H$_2$ transition theoretical model, and show that the theory-observation discrepancy cannot be explained by the intense radiation in W43, nor by variations of the assumed volume density or H$_2$ formation-rate coefficient. We show that the large observed HI columns are naturally explained by several ($9-22$) HI-to-H$_2$ transition layers, superimposed along the sightlines of W43. We discuss other possible interpretations such as a non-steady-state scenario, and inefficient dust absorption. The case of W43 suggests that HI thresholds reported in extra-galactic observations are probably not associated with a single HI-to-H$_2$ transition, but are rather a result of several transition layers (clouds) along the sightlines, beam-diluted with diffuse inter-cloud gas.

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 "$\alpha 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 $\alpha 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.

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.

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 $\alpha G$, and the ratio $I_{\rm UV}/n$, of the FUV intensity to hydrogen gas density. We find $\alpha 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.

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, $\zeta/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 $\zeta/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.

We demonstrate that high abundances of water vapor could have existed in extremely low metallicity ($10^{-3}$ solar) partially shielded gas, during the epoch of first metal enrichment of the interstellar medium of galaxies at high redshifts.