Spectral Line Observations Research Articles

Magritte, a modern software library for spectral line radiative transfer

Spectral line observations are an indispensable tool to remotely probe the physical and chemical conditions throughout the universe. Modelling their behaviour is a computational challenge that requires dedicated software. In this paper, we present the first long-term stable release of Magritte, an open-source software library for line radiative transfer. First, we establish its necessity with two applications. Then, we introduce the overall design strategy and the application/programmer interface (API). Finally, we present three key improvements over previous versions: (1) an improved re-meshing algorithm to efficiently coarsen the spatial discretisation of a model; (2) a variation on Ng-acceleration, a popular acceleration-of-convergence method for non-LTE line transfer; and, (3) a semi-analytic approximation for line optical depths in the presence of large velocity gradients.

Long-term Integration Ability of the Submillimeter Wave Astronomy Satellite (SWAS) Spectral Line Receivers

The Submillimeter Wave Astronomy Satellite (SWAS) was the first space telescope capable of high spectral resolution observations of terahertz spectral lines. We have investigated the integration ability of its two receivers and spectrometer during five and a half years of on-orbit operation. The C i , O2, H2O, and 13CO spectra taken toward all observed Galactic sources were analyzed. The present results are based on spectra with a total integration time of up to 2.72 × 104 hr (≃108 s). The noise in the spectra is generally consistent with that expected from the radiometer equation, without any sign of approaching a noise floor. This noise performance reflects the extremely stable performance of the passively cooled front end as well as other relevant components in the SWAS instrument throughout its mission lifetime.

Observation of spectral lines in the exceptional GRB 221009A

Observation of spectral lines in the exceptional GRB 221009A

Mass Flows in Expanding Coronal Loops

An expansion of the cross-sectional area directly impacts the mass flow along a coronal loop and significantly alters the radiative and hydrodynamic evolution of that loop as a result. Previous studies have found that an area expansion from the chromosphere to the corona significantly lengthens the cooling time of the corona and appears to suppress draining from the corona. In this work, we examine the fluid dynamics to understand how the mass flow rate, the energy balance, and the cooling and draining timescales are affected by a nonuniform area. We find that in loops with moderate or large expansion (cross-sectional area expansion factors of 2, 3, 10, 30, 100 from the photosphere to the apex), impulsive heating, for either direct thermal heating or electron beam heating, induces a steady flow into the corona, so that the coronal density continues to rise during the cooling phase, whereas a uniform loop drains during the cooling phase. The induced upflow carries energy into the corona, balancing the losses from thermal conduction, and continues until thermal conduction weakens enough so that it can no longer support the radiative losses of the transition region. As a result, the plasma cools primarily radiatively until the onset of catastrophic collapse. The speed and duration of the induced upflow both increase in proportion to the rate of area expansion. We argue that observations of blueshifted spectral lines, therefore, could place a constraint on a loop’s area expansion.

The ALMA View of Positive Black Hole Feedback in the Dwarf Galaxy Henize 2–10

Henize 2–10 is a dwarf starburst galaxy hosting a ∼106 M ⊙ black hole (BH) that is driving an ionized outflow and triggering star formation within the central ∼100 pc of the galaxy. Here, we present Atacama Large Millimeter/submillimeter Array continuum observations from 99 to 340 GHz, as well as spectral line observations of the molecules CO (1–0, 3–2), HCN (1–0, 3–2), and HCO+ (1–0, 3–2), with a focus on the BH and its vicinity. Incorporating centimeter-wave radio measurements from the literature, we show that the spectral energy distribution of the BH is dominated by synchrotron emission from 1.4 to 340 GHz, with a spectral index of α ≈ − 0.5. We analyze the spectral line data and identify an elongated molecular gas structure around the BH with a velocity distinct from the surrounding regions. The physical extent of this molecular gas structure is ≈130 pc × 30 pc and the molecular gas mass is ∼106 M ⊙. Despite an abundance of molecular gas in this general region, the position of the BH is significantly offset from the peak intensity, which may explain why the BH is radiating at a very low Eddington ratio. Our analysis of the spatially resolved line ratio between CO J = 3–2 and J = 1–0 implies that the CO gas in the vicinity of the BH is highly excited, particularly at the interface between the BH outflow and the regions of triggered star formation. This suggests that the cold molecular gas is being shocked by the bipolar outflow from the BH, supporting the case for positive BH feedback.

A radical transition in the post-main-sequence system U Equulei

Context. U Equ is an unusual maser-hosting infrared source discovered in the 1990s. It was tentatively classified as a post-asymptotic giant branch (post-AGB) star with a unique optical spectrum displaying rare emission and absorption features from molecular gas at a temperature of about 500 K. In 2022, we serendipitously discovered that its optical spectrum had drastically changed since the last observations in the 1990s. Aims. We aim to characterize the drastic change in the spectrum and analyze the photometric behavior of the object since 1989. Methods. Optical high-resolution spectra of U Equ from the Southern African Large Telescope were supplemented by archival data and near-infrared photometry from the Nordic Optical Telescope. New spectral line observations with the Effelsberg 100 m radio telescope and Atacam Large Millimeter Array are presented. Radiative transfer modeling of multiple epoch spectral energy distributions was performed. Results. No circumstellar molecular features are present in the contemporary optical spectra of U Equ. Nonphotospheric absorption and emission from neutral and ionized species dominate the current spectrum. Some of the observed features indicate an outflow with a projected terminal velocity of 215 km s−1. Broad H&K lines of [Ca II] indicate a photosphere of spectral type F or similar. For the first time, we find SiO J = 1−0 υ=1 maser emission in U Equ. Our collected photometric measurements show that the source has monotonically increased its optical and near-infrared fluxes since about the beginning of this century and continues to do so. The current rise in the optical regime is about 1 mag. Spectral energy distributions at different epochs show dusty circumstellar material that is very likely arranged in a highly inclined disk. Adopting a distance of 4 kpc, informed by the Gaia parallax of U Equ, we find that the source luminosity is about 104 L⊙. This luminosity has likely increased by a factor of a few in the last decades, which is most probably related to the drastic change in the optical circumstellar spectrum of the object. Conclusions. The object has changed considerably in the past three decades, either due to geometrical reconfiguration of the circumstellar medium, evolutionary changes in the central star, or owing to an accretion event that started very recently in the system. Observationally, U Equ appears to resemble category 0 of disk-hosting post-AGB stars reported previously, especially the post-common-envelope binary HD 101584. It is uncertain whether the drastic spectral change and the associated optical and mid-infrared rise in brightness witnessed in U Equ are common in post-AGB stars, but this radical change may be related to the real-time onset of the evolution of the system into a planetary nebula. We find that the post-AGB star V576 Car has undergone a similar transformation as U Equ in the past few decades, which means that the phenomenon is not extremely rare.

Off-nuclear H2O maser and dense molecular gas in NGC 1068

ABSTRACT The results of high-resolution spectral-line observations of dense molecular gas are presented towards the nuclear region of the type 2 Seyfert galaxy NGC 1068. MERLIN observations of the 22 GHz H2O maser were made for imaging the known off-nuclear maser emission at radio jet component located about 0.3 arcsec north-east of the radio nucleus in the galaxy. High angular resolution ALMA observations have spatially resolved the molecular gas emissions of HCN and HCO+ in this region. The off-nuclear maser spots are found to nearly overlap with a ring-like molecular gas structure and are tracing an evolving shock-like structure, which appears to be energized by interaction between the radio jet and circumnuclear medium. The scenario of the dynamic jet–ISM interaction is further supported by a systematic shift of the centroid velocities of the off-nuclear maser features over a period of 35 yr. The enhanced integrated flux ratios of the HCN to HCO+ line emission features at component C suggest a kinetic temperature Tk ≳ 300 K and an H2 density of ≳ 106 cm−3, which are conditions where water masers may be formed. The diagnostics of the masering action in this jet–ISM interaction region is exemplary for galaxies hosting off-nuclear H2O maser emission.

A global view on star formation: The GLOSTAR Galactic plane survey

Context. Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. Aims. We aim to investigate the properties of molecular gas structures on different linear scales with the 4.8 GHz formaldehyde (H2CO) absorption line in Cygnus X. Methods. As part of the GLOSTAR Galactic plane survey, we performed large-scale (7º×3º) simultaneous H2CO (11,0–11,1) spectral line and radio continuum imaging observations toward Cygnus X at λ ~6 cm with the Karl G. Jansky Very Large Array and the Effelsberg 100 m radio telescope. We used auxiliary HI, 13CO (1–0), dust continuum, and dust polarization data for our analysis. Results. Our Effelsberg observations reveal widespread H2CO (11,0–11,1) absorption with a spatial extent of ≳50 pc in Cygnus X for the first time. On large scales of 4.4 pc, the relative orientation between the local velocity gradient and the magnetic field tends to be more parallel at H2 column densities of ≳1.8×1022 cm−2. On the smaller scale of 0.17 pc, our VLA+Effelsberg combined data reveal H2CO (11,0–11,1) absorption only towards three bright HII regions. Our observations demonstrate that H2CO (11,0–11,1) is optically thin in general. The kinematic analysis supports the assertion that molecular clouds generally exhibit supersonic motions on scales of 0.17−4.4 pc. We show a non-negligible contribution of the cosmic microwave background radiation to the extended absorption features in Cygnus X. Our observations suggest that H2CO (11,0–11,1) can trace molecular gas with H2 column densities of ≳5 × 1021 cm−2 (i.e., AV ≳ 5). The ortho-H2CO fractional abundance with respect to H2 has a mean value of 7.0 × 10−10. A comparison of the velocity dispersions on different linear scales suggests that the velocity dispersions of the dominant −3 km s−1 velocity component in the prominent DR21 region are nearly identical on scales of 0.17−4.4 pc, which deviates from the expected behavior of classic turbulence.

The rich molecular environment of the luminous blue variable star AFGL 2298

Context. Luminous blue variable (LBV) stars represent a short-lived stage in the late evolution of the most massive stars. Highly unstable, LBVs exhibit dense stellar winds and episodic eruptions that produce complex circumstellar nebulae, the study of which is crucial for properly constraining the impact of these sources at a Galactic scale from a structural, dynamical, and chemical perspective. Aims. We aim to investigate the molecular environment of AFGL 2298, an obscured Galactic LBV that hosts a highly structured circumstellar environment with hints of multiple mass-loss events in the last few 104 a. Methods. We present spectral line observations of AFGL 2298 at 1 and 3 mm performed with the IRAM 30 m radio telescope. Results. We report the detection of several carbon- and nitrogen-bearing species (CO, 13CO, C18O, C17O, HCO+, HCN, HNC, H13CO+, CN, N2H+, and C2H) in the surroundings of AFGL 2298. We identified three velocity components that clearly stand out from the Galactic background. The morphology, kinematics, masses, and isotopic ratios, together with a comparative study of the fractional abundances, lead us to suggest that two of these components (36 and 70 km s−1) have a stellar origin. The other component (46 km s−1) most likely traces swept-up interstellar material, and probably also harbours a photon-dominated region. Conclusions. We provide the first inventory of the circumstellar molecular gas around AFGL 2298. Our results are compatible with the hypothesis of former mass-loss events produced before the one that created the infrared nebula. The chemistry of this LBV suggests the presence of ejected stellar material, and also swept up gas. These findings will help us to better understand the mass-loss history of this class of evolved massive stars, which is important given that they heavily influence the overall chemical evolution of the Galaxy.

Wide-field CO isotopologue emission and the CO-to-H2 factor across the nearby spiral galaxy M101

Carbon monoxide (CO) emission constitutes the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H2 conversion factor, α12CO(1−0), links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of α12CO(1−0) across all environments within and across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using spectral line observations of CO and its isotopologues mapped across a nearby galaxy, we can constrain the molecular gas conditions and link them to a variation in α12CO(1−0). Here, we present new, wide-field (10 × 10 arcmin2) IRAM 30-m telescope 1 mm and 3 mm line observations of 12CO, 13CO, and C18O across the nearby, grand-design, spiral galaxy M101. From the CO isotopologue line ratio analysis alone, we find that selective nucleosynthesis and changes in the opacity are the main drivers of the variation in the line emission across the galaxy. In a further analysis step, we estimated α12CO(1−0) using different approaches, including (i) via the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) local thermal equilibrium (LTE) based measurements using the optically thin 13CO(1–0) intensity. We find an average value of ⟨α12CO(1 − 0)⟩ = 4.4 ± 0.9 M⊙ pc−2 (K km s−1)−1 across the disk of the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based α12CO(1−0) values are lower by a factor of 2–3 across the disk relative to the dust-based result. Accounting for α12CO(1−0) variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy’s center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if α12CO(1−0) variations are not accounted for.

A CO funnel in the Galactic centre: Molecular counterpart of the northern Galactic chimney

We report the discovery of a velocity coherent, funnel-shaped 13CO emission feature in the Galactic centre (GC) using data from the SEDIGISM survey. The molecular cloud appears as a low-velocity structure (VLSR = [ − 3.5, + 3.5] km s−1) with an angular extent of 0.95° ×1°, extending towards positive Galactic latitudes. The structure is offset from Sgr A* towards negative Galactic longitudes; it spatially and morphologically correlates well with the northern lobe of the 430 pc radio bubble, believed to be the radio counterpart of the multiwavelength GC chimney. Spectral line observations in the frequency range of 85–116 GHz have been carried out using the IRAM 30-m telescope towards 12 positions along the funnel-shaped emission. We examine the 12C/13C isotopic ratios using various molecules and their isotopologues. The mean 12C/13C isotope ratio (30.6 ± 2.9) is consistent with the structure located within inner 3 kpc of the Galaxy and possibly in the GC. The velocity of the molecular funnel is consistent with previous radio recombination line measurements of the northern lobe of radio bubble. Our multiwavelength analysis suggests that the funnel-shaped structure extending over 100 pc above the Galactic plane is the molecular counterpart of the northern GC chimney.

SKA Science Data Challenge 2: analysis and results

ABSTRACT The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarize the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterize 233 245 neutral hydrogen (H i) sources in a simulated data product representing a 2000 h SKA-Mid spectral line observation from redshifts 0.25–0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, ‘reproducibility awards’ were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy – which combined predictions from two independent machine learning techniques to yield a 20 per cent improvement in overall performance – underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical data sets.

ALMA Uncovers Highly Filamentary Structure toward the Sgr E Region

We report on the discovery of linear filaments observed in the CO(1-0) emission for a ∼2′ field of view toward the Sgr E star-forming region, centered at (l, b) = (358.°720, 0.°011). The Sgr E region is thought to be at the turbulent intersection of the “far dust lane” associated with the Galactic bar and the Central Molecular Zone (CMZ). This region is subject to strong accelerations, which are generally thought to inhibit star formation, yet Sgr E contains a large number of H ii regions. We present 12CO(1-0), 13CO(1-0), and C18O(1-0) spectral line observations from the Atacama Large Millimeter/submillimeter Array and provide measurements of the physical and kinematic properties for two of the brightest filaments. These filaments have widths (FWHMs) of ∼0.1 pc and are oriented nearly parallel to the Galactic plane, with angles from the Galactic plane of ∼2°. The filaments are elongated, with lower-limit aspect ratios of ∼5:1. For both filaments, we detect two distinct velocity components that are separated by about 15 km s−1. In the C18O spectral line data, with ∼0.09 pc spatial resolution, we find that these velocity components have relatively narrow (∼1–2 km s−1) FWHM line widths when compared to other sources toward the Galactic center. The properties of these filaments suggest that the gas in the Sgr E complex is being “stretched,” as it is rapidly accelerated by the gravitational field of the Galactic bar while falling toward the CMZ, a result that could provide insights into the extreme environment surrounding this region and the large-scale processes that fuel this environment.

CASA on the Fringe—Development of VLBI Processing Capabilities for CASA

New functionality to process Very Long Baseline Interferometry (VLBI) data has been implemented in the CASA package. This includes two new tasks to handle fringe fitting and VLBI-specific amplitude calibration steps. Existing tasks have been adjusted to handle VLBI visibility data and calibration meta-data properly. With these updates, it is now possible to process VLBI continuum and spectral line observations in CASA. This article describes the development and implementation, and presents an outline for the workflow when calibrating European VLBI Network or Very Long Baseline Array data in CASA. Though the CASA VLBI functionality has already been vetted extensively as part of the Event Horizon Telescope data processing, in this paper we compare results for the same data set processed in CASA and AIPS. We find identical results for the two packages and conclude that CASA in some cases performs better, though it cannot match AIPS for single-core processing time. The new functionality in CASA allows for easy development of pipelines or Jupyter notebooks, and thus contributes to raising VLBI data processing to present day standards for accessibility, reproducibility, and reusability.

Formation of H2CO and CH3OH in the Protoplanetary Disks

It is known that the complex organic molecules in planetary systems may be formed by the organic chemistry of gases and particles in protoplanetary disks. Among these molecules, H2CO and CH3OH are potential precursors of essential organic compounds like glycolaldehyde. We calculate the physical and chemical structures of the protoplanetary disks around LkCa 15 and MWC 480 stars, compare the obtained column densities with the observed values. Our results mostly show reasonable agreements with observations. We identify the formation of H2CO to be through the both gas-phase and dust surface reactions, while CH3OH is dominantly formed by surface chemistry in the original astrolabe and then released into the gas-phase by surface exothermic chemical reactions in the LkCa 15 disk. We studied the synthesis of complex organic molecules (NH2CHO and CH3OCHO) in the protoplanetary disk to predict the observable species and Molecular abundance. In the disks around LkCa 15 and MWC 480 stars, we suggest that the strong spectral line observation of NH2CHO belonging to band 4 may obtain good results with ALMA, the observation of CH3OCHO may be challenging.

Probing fragmentation with ALMA continuum and spectral line observations of the dense clumps in the ℓ = 224° region

ABSTRACT We report observations, performed with the Atacama Large Millimeter/submillimeter Array (ALMA), of 1 mm dust continuum emission and molecular line emission in 13CO(2–1) and C18O(2–1), towards a sample of starless and protostellar clumps selected from a region, towards the ℓ = 224° field, of the Herschel Infrared GALactic Plane Survey (Hi-GAL). Using the ALMA images and a source extraction algorithm we have analysed the small-scale (∼1000 AU) structure of the clumps and their population of cores (or fragments). We find in general multiple cores in each Hi-GAL clump, both in the continuum and spectral lines, but we do not find a dominant fragmentation mode and the morphologies are very different among the various sources. Our results suggest that during the transition phase from clump to core, those sources with a higher core formation efficiency are also associated with parent clumps that are more likely to convert a higher fraction of their initial mass into a single or a few cores. We were able to obtain a core mass function, or CoMF, covering masses in the range ∼2 × 10−3 to ∼1 M⊙ for the C18O cores, and ∼4 × 10−2 to ∼10 M⊙ for the continuum cores. We find that the CoMF in our sample is much shallower than the higher mass ($\gtrsim 1$ M⊙) IMF, thus indicating that while approaching the final phase of fragmentation the mass function does not resemble the IMF more closely.

First looks at solar active regions with ALMA

During the first few years of observing the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), the scientific community has acquired a number of observational datasets targeting various structures in active regions, including sunspot umbra and penumbra, active region pores, and plages. In this paper we review the results obtained from the extensive analysis of these interferometric millimeter data, together with the coordinated observations from IRIS, SDO, IBIS, and Hinode, that reveal information on the chromospheric thermal structure above active regions and properties of small-scale heating events near magnetic field concentrations. We discuss the properties of waves (especially the three-minute oscillations) in sunspots, plage, and network. We speculate how high-resolution millimeter data can supplement spectral line observations in the visible and UV and can improve chromospheric spectroscopic inversions. We identify challenges in the interpretation of the millimeter continuum emission due to the complex, non-local and time-dependent processes that determine the electron density through the chromosphere. Finally we overview the prospects for future active regions observations with ALMA during the ascending phase of the solar cycle.

Synthetic Next Generation Very Large Array line observations of a massive star-forming cloud

Context. Studies of the interstellar medium and the pre-stellar cloud evolution require spectral line observations that have a high sensitivity and high angular and velocity resolution. Regions of high-mass star formation are particularly challenging because of line-of-sight confusion, inhomogeneous physical conditions, and potentially very high optical depths. Aims. We wish to quantify to what accuracy the physical conditions within a massive star-forming cloud can be determined from observations. We are particularly interested in the possibilities offered by the Next Generation Very Large Array (ngVLA) interferometer. Methods. We used data from a magnetohydrodynamic simulation of star formation in a high-density environment. We concentrated on the study of a filamentary structure that has physical properties similar to a small infrared-dark cloud. We produced synthetic observations for spectral lines observable with the ngVLA and analysed these to measure column density, gas temperature, and kinematics. Results were compared to ideal line observations and the actual 3D model. Results. For a nominal cloud distance of 4kpc, ngVLA provides a resolution of ~0.01 pc even in its most compact configuration. For abundant molecules, such as HCO+, NH3, N2H+, and CO isotopomers, cloud kinematics and structure can be mapped down to subarcsecond scales in just a few hours. For NH3, a reliable column density map could be obtained for the entire 15″ × 40″ cloud, even without the help of additional single-dish data, and kinetic temperatures are recovered to a precision of ~1 K. At higher frequencies, the loss of large-scale emission becomes noticeable. The line observations are seen to accurately trace the cloud kinematics, except for the largest scales, where some artefacts appear due to the filtering of low spatial frequencies. The line-of-sight confusion complicates the interpretation of the kinematics, and the usefulness of collapse indicators based on the expected blue asymmetry of optically thick lines is limited. Conclusions. The ngVLA will be able to provide accurate data on the small-scale structure and the physical and chemical state of star-forming clouds, even in high-mass star-forming regions at kiloparsec distances. Complementary single-dish data are still essential for estimates of the total column density and the large-scale kinematics.

NaCl and KCl in Io’s Atmosphere

We present the first comprehensive study of NaCl and KCl gases in Io’s atmosphere in order to investigate their characteristics and infer properties of Io’s volcanoes and subsurface magma chambers. In this work, we compile all past spectral line observations of NaCl and KCl in Io's atmosphere from the Atacama Large Millimeter/submillimeter Array and use atmospheric models to constrain the physical properties of the gases on several dates between 2012 and 2018. NaCl and KCl appear to be largely spatially confined, and for observations with high spectral resolution, the temperatures are high (∼500–1000 K), implying a volcanic origin. The ratio of NaCl:KCl was found to be ∼5–6 in 2015 June and ∼3.5–10 in 2016 June, which is consistent with predictions based on observations of Io's extended atmosphere and less than half the Na/K ratio in chondrites. Assuming these gases are volcanic in origin, these ratios imply a magma temperature of ∼1300 K, such that the magma will preferentially outgas KCl over NaCl.

Low-energy cosmic rays: regulators of the dense interstellar medium

Low energy cosmic rays (up to the GeV energy domain) play a crucial role in the physics and chemistry of the densest phase of the interstellar medium. Unlike interstellar ionising radiation, they can penetrate large column densities of gas, and reach molecular cloud cores. By maintaining there a small but not negligible gas ionisation fraction, they dictate the coupling between the plasma and the magnetic field, which in turn affects the dynamical evolution of clouds and impacts on the process of star and planet formation. The cosmic-ray ionisation of molecular hydrogen in interstellar clouds also drives the rich interstellar chemistry revealed by observations of spectral lines in a broad region of the electromagnetic spectrum, spanning from the submillimetre to the visual band. Some recent developments in various branches of astrophysics provide us with an unprecedented view on low energy cosmic rays. Accurate measurements and constraints on the intensity of such particles are now available both for the very local interstellar medium and for distant interstellar clouds. The interpretation of these recent data is currently debated, and the emerging picture calls for a reassessment of the scenario invoked to describe the origin and/or the transport of low energy cosmic rays in the Galaxy.