Search

The quantum dynamics of carbon chains through H2 and He collisions in the interstellar medium (ISM) is an important step toward accurate modeling of their abundance in non-local thermodynamic equilibrium conditions. The C5 ( Σ g + 1 ) molecule is the longest pure carbon chain detected in the ISM to date. While He collisions are computationally easy to perform, the collision with much more abundant H2 is both complicated and computationally demanding. Using templates for approximating p-H2 collisional rates, such as scaling He rates and using a reduced 4D → 2D potential energy surface (PES), has limited applicability. On the other hand, any such approximation does not exist for o-H2. Therefore, a full rotational dynamics of C5 with both p- and o-H2 is performed considering both molecules as rigid-rotors. The PES is calculated using CCSD(T)-F12a/AVTZ, and a neural network fitting model has been carefully chosen to strictly obey spectroscopic accuracy and augment the PES. The augmented PES is then expanded into radial terms using the bispherical harmonics function, and close coupling calculations have been done to get the cross sections and, subsequently, rate coefficients for various rotational transitions of C5. [ABSTRACT FROM AUTHOR]

The pure rotational and rovibrational spectra of the ν27 –NH2 torsion of cyclopropylamine (CPA) in the far-infrared region were measured with a high-resolution Fourier transform infrared coupled to a synchrotron. The complex spectra reflect the presence of both trans and gauche conformers. Analysis of the pure rotational spectra (34–64 cm−1) yielded accurate rotational and centrifugal distortion constants of the ground and first two torsional excited states of trans-CPA. The fundamental, hot bands and weak overtones were identified and assigned in the 200–550 cm−1 range. Global analysis of over 19 000 transitions provides accurate energy levels of the torsional polyads up to vT = 3. The torsional levels and their rotational constants were in agreement with the theoretical results from quasiadiabatic channel reaction path Hamiltonian (RPH) calculations, emphasizing the need for molecular-specific theoretical treatments for large amplitude motions. Tunneling components of the torsional fundamental of gauche-CPA were assigned based on the RPH results and symmetry considerations, differing from previous experimental and theoretical work. This comprehensive spectroscopic characterization of CPA is crucial for its potential detection in the interstellar medium as a precursor to complex prebiotic molecules, providing essential data for future astronomical searches and advancing our understanding of nitrogen-containing organic molecules in space. [ABSTRACT FROM AUTHOR]

The rotational dynamics of isocyanogen (CNCN) is studied for its collision with para (p-) and ortho (o-) hydrogen (H2) in the temperature range of 1–100 K. These temperatures correspond to the cold dense molecular clouds in the interstellar medium where molecular hydrogen is the primary collider. An ab initio 4D potential energy surface (PES) is constructed keeping the two molecules under rigid rotor approximation. The PES is generated using the CCSD(T)-F12b/AVTZ level of theory. The 4D PES is further fitted into a neural network (NN) model, which can augment the surface and account for missing data points within spectroscopic accuracy. This NN-fitted PES is then expanded over a bispherical harmonics function to get radial terms, which are expressed into analytic functions. Thereafter, the cross sections (σ) are computed for rotational transitions of CNCN (j → j′) using the close-coupling and centrifugal sudden methods for both p-H2 (jc = 0) and o-H2 (jc = 1) collision till 194 cm−1. In addition, p-H2 (jc = 0, 2) cross sections are also computed using the centrifugal sudden approximation method. The collisional rates are achieved by taking the Boltzmann distribution of σ over the translational energy of H2 till 100 K. Finally, the CNCN–H2 rates are compared to CNCN–He and NCCN–H2 collisional rates. Comparing even and odd transitions for the CNCN–H2 rates show a propensity toward higher rates for even transitions especially for o-H2 collisions considering low-order transitions. [ABSTRACT FROM AUTHOR]

Chloromethane, CH3Cl, is the first organohalogen molecule to be detected in the interstellar medium. Using the recently generated accurate three-dimensional potential energy surface (3D-PES) for the weakly bound CH335Cl–He complex, we deduced that of CH337Cl–He. Both 3D-PESs were incorporated later into dynamical computations code for nuclear motions treatments of each isotopologue colliding with He. We determined the state-to-state (de-)excitation cross sections for transitions among the 105 lowest rotational levels of both ortho-CH335Cl and ortho-CH337Cl (up to jkɛ = 206±). For total energies up to 100 cm−1, we used the quantum close-coupling methodology, and we adopted the coupled state approach for higher energies (up to 2500 cm−1). Rate coefficients for temperatures ranging from 5 to 150 K are then deduced after averaging these cross sections over a Maxwell–Boltzmann velocity distribution. Computations show that isotopic effects are weak, validating the use of collision rates of ortho-CH335Cl for ortho-CH337Cl and vice-versa. However, we expect that our results will be helpful to accurately determine the abundances of these isotopologues in the interstellar medium and, more generally, for understanding the complex and still not fully elucidated chemistry of chlorine containing molecular species there. [ABSTRACT FROM AUTHOR]

Structure determination is necessary to identify unknown organic molecules, such as those in natural products, forensic samples, the interstellar medium, and laboratory syntheses. Rotational spectroscopy enables structure determination by providing accurate 3D information about small organic molecules via their moments of inertia. Using these moments, Kraitchman analysis determines isotopic substitution coordinates, which are the unsigned |x|, |y|, |z| coordinates of all atoms with natural isotopic abundance, including carbon, nitrogen, and oxygen. While unsigned substitution coordinates can verify guesses of structures, the missing +/− signs make it challenging to determine the actual structure from the substitution coordinates alone. To tackle this inverse problem, we develop Kreed (Kraitchman REflection-Equivariant Diffusion), a generative diffusion model that infers a molecule's complete 3D structure from only its molecular formula, moments of inertia, and unsigned substitution coordinates of heavy atoms. Kreed's top-1 predictions identify the correct 3D structure with near-perfect accuracy on large simulated datasets when provided with substitution coordinates of all heavy atoms with natural isotopic abundance. Accuracy decreases as fewer substitution coordinates are provided, but is retained for smaller molecules. On a test set of experimentally measured substitution coordinates gathered from the literature, Kreed predicts the correct all-atom 3D structure in 25 of 33 cases, demonstrating experimental potential for de novo 3D structure determination with rotational spectroscopy. [ABSTRACT FROM AUTHOR]

There is an ongoing debate on the apparent depletion of sulfur in the interstellar medium (ISM) compared to its universal abundance; therefore, the investigation of sulfurous compounds at low temperatures is of utmost importance. This work aims to study thioacetamide, H3C–C(=S)–NH2, in low-temperature inert Ar and para-H2 matrices by IR spectroscopy. The samples have been exposed to various sources of irradiation, such as Lyman-α or laser UV photons as well as energetic electrons. Using different host materials enabled assessing the matrix's impact on precursor decomposition. The response of the molecule to different types of irradiation has also been evaluated. The existence of three main decomposition channels were deduced: formation of (i) CH3, CH4, and HNCS; (ii) H2S and H2C=C=NH; and (iii) NH3 and H2C=C=S. The H3C–CN and H3C–NC isomers of H2C=C=NH could also be identified. Secondary products such as HNC and HCN were also detected in the quantum solid para-H2 in contrast to the more rigid Ar matrix. The listed decomposition products have been observed in the ISM, with the exception of H2C=C=NH and H3C–NC. The results point to the potential sensitivity of the precursor molecule to energetic radiation in space environments. Finally, the findings of this work will serve as a foundation for future irradiation experiments using the astrochemically more relevant pure thioacetamide ice. [ABSTRACT FROM AUTHOR]

With the advent of JWST ice observations, dedicated studies on the formation reactions of detected molecules are becoming increasingly important. One of the most interesting molecules in interstellar ice is CO2. Despite its simplicity, the main formation reaction considered, CO + OH → CO2 + H through the energetic HOCO* intermediate on ice dust, is subject to uncertainty because it directly competes with the stabilization of HOCO as a final product, which is formed through energy dissipation of HOCO* to the water ice. When energy dissipation to the surface is effective during the reaction, HOCO can be a dominant product. In this study, we experimentally demonstrate that the major product of the reaction is indeed not CO2, but rather the highly reactive radical HOCO. The HOCO radical can later evolve into CO2 through H-abstraction reactions, but these reactions compete with additional reactions, leading to the formation of carboxylic acids (R-COOH). Our results highlight the importance of HOCO chemistry and encourage further exploration of the chemistry of this radical. [ABSTRACT FROM AUTHOR]

We assess the impact of quasi-stellar objects (QSOs) on the high-redshift (⁠|$z\gt 4$|⁠) intergalactic medium using Monte Carlo realizations of QSO populations and the He  |$\rm \scriptstyle III$| regions they generate, applied to the Sherwood–Relics simulations, allowing for uncertainties in the QSO luminosity function, its evolution, and QSO spectra and ages. While QSO luminosity functions based on optical–infrared selection are unable to reproduce the broadening H  |$\rm \scriptstyle I$| Ly  |$\alpha$| optical depth distributions at |$z\gt 5$|⁠ , much broader distributions are found for the higher numbers of QSOs based on X-ray selection, suggesting a large QSO contribution to the ultraviolet background at |$z\gt 5$| may offer an alternative to late reionization models to account for the broad H  |$\rm \scriptstyle I$| Ly  |$\alpha$| optical depth distributions. Realizations using QSOs based on the higher QSO counts also much better recover the measured pixel flux autocorrelation function at |$z\gt 5$|⁠. The He  |$\rm \scriptstyle III$| regions from QSO sources according to both types of luminosity function suppress the pixel flux power spectrum on small scales, |$k\gt 0.02\, \mathrm{s\, km^{-1}}$|⁠ , while enhancing it on larger, both by amounts of up to tens of per cent at |$z\gt 4$|⁠ , with the magnitude increasing with the intergalactic He  |$\rm \scriptstyle III$| filling factor and the boost in temperature within the He  |$\rm \scriptstyle III$| regions. [ABSTRACT FROM AUTHOR]

A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multipath propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present an improved method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, |$\alpha$|⁠ , as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report the measurements of the scatter-broadening time and |$\alpha$| from analysis of PSR J1643 |$-$| 1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and |$\alpha$| was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of this technique to PSR J1643 |$-$| 1224. [ABSTRACT FROM AUTHOR]

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be |$R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$|⁠. The mass-normalized rate is |$R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$| and the volumetric rate is |$R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$|⁠. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included. [ABSTRACT FROM AUTHOR]

Accurate astrometric and photometric measurements from Gaia have led to the construction of 3D dust extinction maps which can now be used for estimating the integrated extinctions of Galactic sources located within 5 kpc. These maps based on optical observations may not be reliable for use in the ultraviolet (UV) which is more sensitive to reddening. Past studies have focussed on studying UV extinction using main-sequence stars but lack comparison with 3D dust maps. White dwarfs with well-modelled hydrogen-dominated (DA) atmospheres provide an advantage over main-sequence stars affected by magnetic activity. In this work, we study the variation of UV extinction with 3D dust maps utilizing Hubble Space Telescope (HST) and Galaxy Evolution Explorer (GALEX) observations of DA white dwarfs located within 300 pc. We used HST COS spectroscopic data of 76 sightlines to calculate the optical extinction from Si  ii column densities and validate our results with the kinematic model predictions of the local interstellar medium. Also, we combined GALEX and Gaia photometric observations of 1158 DA white dwarfs to study UV reddening by comparing observed and modelled colour–colour relations. We calculated GALEX non-linearity corrections and derived reddening coefficients [ |$R(NUV-G)=6.52\pm 1.53$| and |$R(FUV-G)=6.04\pm 2.41$| ] considering their variations with optical extinction (⁠|$\rm{A_{V}}\lt 0.1$|  mag), and found them to be in good agreement with known extinction laws. HST analysis suggests a positive bias of 0.01–0.02 mag in the optical extinction from 3D maps depending on the Galactic latitude. These results independently confirm the validity of 3D dust maps to de-redden the optical and UV observations of white dwarfs. [ABSTRACT FROM AUTHOR]

Being the most prominent H  i line, Ly |$\alpha$| permeates the cosmic web in emission. Despite its potential as a cosmological probe, its detection on large scales remains elusive. We present a new methodology to perform Ly |$\alpha$| intensity mapping with broad-band optical images, by cross-correlating them with Ly |$\alpha$| forest data using a custom one-parameter estimator. We also develop an analytical large-scale Ly |$\alpha$| emission model with two parameters (average luminosity |$\langle L_{\rm Ly\alpha } \rangle$| and bias |$b_{\rm e}$|⁠) that respects observational constraints from quasar (QSO) luminosity functions. We compute a forecast for Dark Energy Camera Legacy Survey (DECaLS)/Beijing–Arizona Sky Survey (BASS) g -band images cross-correlated with Dark Energy Spectroscopic Instrument (DESI) Ly |$\alpha$| forest data, setting guidelines for reducing images into Ly |$\alpha$| intensity maps. Given the transversal scales of our cross-correlation (26.4 arcmin, |$\sim$| 33 cMpc h−1), our study effectively integrates Ly |$\alpha$| emission over all the cosmic volume inside the DESI footprint at |$2.2 \lt z \lt 3.4$| (the g -band Ly |$\alpha$| redshift range). Over the parameter space (⁠|$\langle L_{\rm Ly\alpha } \rangle$|⁠ , |$b_{\rm e}$|⁠) sampled by our forecast, we find a 3 |$\sigma$| of large-scale structure in Ly |$\alpha$| likely, with a probability of detection of 23.95 per cent for DESI–DECaLS/BASS, and 54.93 per cent for a hypothetical DESI phase-II with twice as much Ly |$\alpha$| QSOs. Without a detection, we derive upper bounds on |$\langle L_{\rm Ly\alpha } \rangle$| competitive with optimistic literature estimates [ |$2.3 \pm 1 \times 10^{\rm 41}$| erg s−1 (cMpc |$^3$|⁠)−1 for DESI, and |$\sim$| 35  per cent lower for its hypothetical phase-II]. Extrapolation to the DESI-Rubin overlap shows that a detection of large-scale structure with Ly |$\alpha$| intensity mapping using next-generation imaging surveys is certain. Such detection would allow constraining |$\langle L_{\rm Ly\alpha } \rangle$|⁠ , and explore the constraining power of Ly |$\alpha$| intensity mapping as a cosmological probe. [ABSTRACT FROM AUTHOR]

We study the properties of long gamma-ray burst (GRB) host galaxies using a statistical modelling framework derived to model damped Lyman- |$\alpha$| absorbers (DLAs) in quasar spectra at high redshift. The distribution of |$N_{\rm H\, {\small I}}$| for GRB-DLAs is |$\sim$| 10 times higher than what is found for quasar-DLAs at similar impact parameters. We interpret this as a temporal selection effect due to the short-lived GRB progenitor probing its host at the onset of a starburst where the interstellar medium may exhibit multiple overdense regions. Owing to the larger |$N_{\rm H\, {\small I}}$|⁠ , the dust extinction is larger with 29 per cent of GRB-DLAs exhibiting |$A(V)\gt 1$|  mag in agreement with the fraction of 'dark bursts'. Despite the differences in |$N_{\rm H\, {\small I}}$| distributions, we find that high-redshift |$2 \lt z \lt 3$| quasar- and GRB-DLAs trace the luminosity function of star-forming host galaxies in the same way. We propose that their differences may arise from the fact that the galaxies are sampled at different times in their star formation histories, and that the absorption sightlines probe the galaxy haloes differently. Quasar-DLAs sample the full H  i cross-section, whereas GRB-DLAs sample only regions hosting cold neutral medium. Previous studies have found that GRBs avoid high-metallicity galaxies (⁠|$\sim$| 0.5  |$Z_{\odot }$|⁠). Since at these redshifts galaxies on average have lower metallicities, our sample is only weakly sensitive to such a threshold. Lastly, we find that the modest detection rate of cold gas (H |$_2$| or C  i) in GRB spectra can be explained mainly by a low volume filling factor of cold gas clouds and to a lesser degree by destruction from the GRB explosion itself. [ABSTRACT FROM AUTHOR]

Understanding the formation of carbonyl sulfide (OCS) in interstellar ices is key to constrain the sulfur chemistry in the interstellar medium (ISM), since it is the only ice S-bearing molecule securely detected thus far. Two general pathways for OCS formation have been proposed: sulfurization of CO (CO + S) and oxidation of CS (CS + O), but their relative contribution in interstellar ices remains unconstrained. We have evaluated the contribution of both pathways to OCS formation upon energetic processing in isotopically labelled CO |$_2$|⁠ :CS |$_2$| and CO:CS |$_2$| ice samples at 7 |$-$| 50 K. Our results indicated that formation of OCS through the CS + O pathway was more favourable than through the CO + S pathway, as previously suggested by theoretical calculations. In addition, its relative contribution increased at higher temperatures. Therefore, this pathway could play a role in the ice formation of OCS, especially in warm regions where CO is expected to be preferentially in the gas phase. At the same time, we have explored the chemistry of CS |$_2$| -bearing, CO |$_2$| -, CO-, and also H |$_2$| O-rich ices, that could be relevant to the sulfur interstellar chemistry. We observed formation of a variety of S-bearing products in addition to OCS, including SO |$_2$|⁠ , C |$_3$| S |$_2$|⁠ , and S |$_2$|⁠. However, a significant fraction of sulfur was not detected at the end of the experiments, and could be locked in long, undetectable sulfur allotropes, one of the potential carriers of the missing sulfur in the dense ISM. [ABSTRACT FROM AUTHOR]

Strong low-ionization transitions like the Mg  ii |$\lambda 2796,2804$|  Å doublet are believed to produce a detectable 'metal-line forest', if metals pollute the neutral intergalatic medium (IGM). We measure the autocorrelation of the Mg  ii forest transmission using 10 ground-based |$z \ge 6.80$| quasar spectra probing the redshift range |$5.96 \lt z_{\rm Mg\, {\small II}} \lt 7.42$| (⁠|$z_{\rm Mg\, {\small II}, median} = 6.47$|⁠). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (⁠|$\Delta v = 768~{\rm km\, s^{-1}}$|⁠) arising from discrete absorbers in the circumgalactic medium of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a seminumerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95 per cent credibility upper limit of |$[{\rm Mg/H}] \lt -3.73$| at |$z_{\rm Mg\, {\small II},median} = 6.47$|⁠ , assuming uninformative priors on [Mg/H] and the IGM neutral fraction |$x_{\rm {H\, {\small I}}}$|⁠. Splitting the data into low- z (⁠|$5.96 \lt z_{\rm Mg\, {\small II}} \lt 6.47$|⁠ ; |$z_{\rm Mg\, {\small II},median} = 6.235$|⁠) and high- z (⁠|$6.47 \lt z_{\rm Mg\, {\small II}} \lt 7.42$|⁠ ; |$z_{\rm Mg\, {\small II},median} = 6.72$|⁠) subsamples again yields null detections and 95 per cent upper limits of |$[{\rm Mg/H}] \lt -3.75$| and |$\lt -3.45$|⁠ , respectively. These first measurements set the stage for making the Mg ii forest an emerging tool to precisely constrain the Universe reionization and enrichment history. [ABSTRACT FROM AUTHOR]

Using a combination of Hubble Space Telescope (HST), JWST, and Atacama Large Millimeter/submillimeter Array (ALMA) data, we perform spatially resolved spectral energy distributions (SED) fitting of fourteen 4 < z < 6 ultraviolet (UV)-selected main-sequence galaxies targeted by the ALMA Large Program [C ii ] Resolved ISM in Star-forming Galaxies. We consistently model the emission from stars and dust in ∼0.5–1 kpc spatial bins to obtain maps of their physical properties. We find no offsets between the stellar masses (M *) and star formation rates (SFRs) derived from their global emission and those from adding up the values in our spatial bins, suggesting there is no bias of outshining by young stars on the derived global properties. We show that ALMA observations are important to derive robust parameter maps because they reduce the uncertainties in L dust (hence, AV and SFR). Using these maps, we explore the resolved star-forming main sequence for z ∼ 5 galaxies, finding that this relation persists in typical star-forming galaxies in the early Universe. We find less obscured star formation where the M * (and SFR) surface densities are highest, typically in the central regions, contrary to the global relation between these parameters. We speculate this could be caused by feedback driving gas and dust out of these regions. However, more observations of IR luminosities with ALMA are needed to verify this. Finally, we test empirical SFR prescriptions based on the UV+IR and [C ii ] line luminosity, finding they work well at the scales probed (approximately kiloparsec). Our work demonstrates the usefulness of joint HST-, JWST-, and ALMA-resolved SED modeling analyses at high redshift. [ABSTRACT FROM AUTHOR]

Progenitor stars of long gamma-ray bursts (GRBs) could be surrounded by a significant and complex nebula structure lying at a parsec-scale distance. After the initial release of energy from the GRB jet, the jet will interact with this nebula environment. We show here that for a large, plausible parameter space region, the interaction between the jet blast wave and the wind termination (reverse) shock is expected to be weak, and may be associated with a precursor emission. As the jet blast wave encounters the contact discontinuity separating the shocked wind and the shocked interstellar medium, we find that a bright flash of synchrotron emission from the newly formed reverse shock is produced. This flash is expected to be observed at around ∼100 s after the initial explosion and precursor. Such a delayed emission thus constitutes a circumburst medium (CBM) phase in a GRB, having a physically distinct origin from the preceding prompt phase and the succeeding afterglow phase. The CBM phase emission may thus provide a natural explanation for bursts observed to have a precursor followed by an intense, synchrotron-dominated main episode that is found in a substantial minority, ∼10% of GRBs. A correct identification of the emission phase is thus required to infer the properties of the flow and of the immediate environment around GRB progenitors. [ABSTRACT FROM AUTHOR]

We investigate the long-term persistence of the time series of energetic neutral atom (ENA) fluxes recorded by the Interstellar Boundary Explorer (IBEX). ENAs provide global information about the outer heliosphere and its interactions with the very local interstellar medium. To avoid the IBEX Ribbon, here, we focus our analysis solely on the polar regions N60°–N90° and S60°–S90°. Each year, IBEX takes two measurements of every pixel in the sky. We make use of the whole set of 14 yr of IBEX data and adhere to the correct time order for the construction of the flux time series. We examine in detail both the trend and the fluctuations of these time series. Using modern methods of time series analysis and persistence characterization, we show that the time series (i) trend is influenced by the solar cycle, (ii) persistence can be established independently of the presence of this trend, and (iii) statistical properties of the fluctuations differs between north and south, pointing to the existence of anisotropy and thus a north–south asymmetry. [ABSTRACT FROM AUTHOR]

We present the first Chandra X-ray Observatory (CXO) catalog of "pulsar X-ray filaments," or "misaligned outflows." These are linear, synchrotron-radiating features powered by ultrarelativistic electrons and positrons that escape from bow shock pulsars. The filaments are misaligned with the (large) pulsar velocity, distinguishing them from the pulsar wind nebula (PWN) trail, which is also often visible in CXO ACIS images. Spectral fits and morphological properties are extracted for five secure filaments and three candidates using a uniform method. We present a search of archival CXO data for linear diffuse features; the known examples are recovered and a few additional weak candidates are identified. We also report on a snapshot CXO ACIS survey of pulsars with properties similar to the filament producers, finding no new filaments but some diffuse emission, including one PWN trail. Finally, we provide an updated model for the pulsar properties required to create filaments in light of these new observations. [ABSTRACT FROM AUTHOR]

The unique biosignature of life on Earth is the homochirality of organic compounds, such as amino acids, proteins, and sugars. High-energy spin-polarized (spinup or spindown) electrons (SPEs) from the β decay of radioactive nuclei have been proposed as a source of symmetry breaking, leading to homochirality; however, their exact role is much debated. Here, we propose magnetically aligned dust grains as a new source of SPEs, due to the photoemission of electrons having aligned spins by the Barnett effect. For the interstellar UV radiation field of strength G UV, we find that the SPE emission rate is Γ pe SPE ∼ 10 − 14 G UV electrons per second per H, the fraction of spin-polarized to total photoelectrons is ∼10%, and the SPE yield (photoelectron number per UV photon) can reach ∼1%, using the modern theory of grain alignment. SPEs emitted from aligned grains could play an important role in chiral-induced spin-selectivity-driven reduction chemistry in the icy grain mantles, producing an enantiomer excess of chiral molecules formed on the grain mantle. Finally, we suggest magnetically aligned grains could directly impact the enantio-selectivity through the chiral-induced spin-selective adsorption effect and exchange interaction. We estimate the disalignment of electron spins and depolarization by elastic scattering using the Mott theory and find that these effects are negligible for low-energy SPEs, so that the spins of SPEs remain well aligned during their journey through dust grains and the gas. Our proposed mechanism might explain the chiral asymmetry of prebiotic molecules detected in numerous comets, asteroids, and meteorites. [ABSTRACT FROM AUTHOR]

Polycyclic aromatic hydrocarbons (PAHs) are organic molecules containing adjacent aromatic rings. Infrared emission bands show that PAHs are abundant in space, but only a few specific PAHs have been detected in the interstellar medium. We detected 1-cyanopyrene, a cyano-substituted derivative of the related four-ring PAH pyrene, in radio observations of the dense cloud TMC-1, using the Green Bank Telescope. The measured column density of 1-cyanopyrene is ~1.52 × 1012 cm−2 , from which we estimate that pyrene contains up to 0.1% of the carbon in TMC-1. This abundance indicates that interstellar PAH chemistry favors the production of pyrene. We suggest that some of the carbon supplied to young planetary systems is carried by PAHs that originate in cold molecular clouds. [ABSTRACT FROM AUTHOR]

We present the largest Galactic neutral hydrogen H i absorption survey to date, utilizing the Australian SKA Pathfinder Telescope at an unprecedented spatial resolution of 30 arcsec. This survey, GASKAP-H i , unbiasedly targets 2714 continuum background sources over 250 square degrees in the direction of the Magellanic Clouds, a significant increase compared to a total of 373 sources observed by previous Galactic absorption surveys across the entire Milky Way. We aim to investigate the physical properties of cold (CNM) and warm (WNM) neutral atomic gas in the Milky Way foreground, characterized by two prominent filaments at high Galactic latitudes (between |$-45^{\circ }$| and |$-25^{\circ }$|⁠). We detected strong H i absorption along 462 lines of sight above the 3 |$\sigma$| threshold, achieving an absorption detection rate of 17 per cent. GASKAP-H i 's unprecedented angular resolution allows for simultaneous absorption and emission measurements to sample almost the same gas clouds along a line of sight. A joint Gaussian decomposition is then applied to absorption-emission spectra to provide direct estimates of H i optical depths, temperatures, and column densities for the CNM and WNM components. The thermal properties of CNM components are consistent with those previously observed along a wide range of Solar neighbourhood environments, indicating that cold H i properties are widely prevalent throughout the local interstellar medium. Across our region of interest, CNM accounts for |$\sim$| 30 per cent of the total H i gas, with the CNM fraction increasing with column density towards the two filaments. Our analysis reveals an anticorrelation between CNM temperature and its optical depth, which implies that CNM with lower optical depth leads to a higher temperature. [ABSTRACT FROM AUTHOR]

With numerical simulations that employ adaptive ray-tracing (ART) for radiative transfer at the same time as evolving gas magnetohydrodynamics, thermodynamics, and photochemistry, it is possible to obtain a high-resolution view of ultraviolet (UV) fields and their effects in realistic models of the multiphase interstellar medium. Here, we analyze results from TIGRESS-NCR simulations, which follow both far-UV (FUV) wavelengths, important for photoelectric heating and polycyclic aromatic hydrocarbon excitation, and the Lyman continuum (LyC), which photoionizes hydrogen. Considering two models, representing solar neighborhood and inner-galaxy conditions, we characterize the spatial distribution and time variation of UV radiation fields, and quantify their correlations with gas. We compare four approximate models for the FUV to simulated values to evaluate alternatives when full ART is infeasible. By convolving FUV radiation with density, we produce mock maps of dust emission. We introduce a method to calibrate mid-IR observations, for example from JWST, to obtain high-resolution gas surface density maps. We then consider the LyC radiation field, finding most of the gas exposed to this radiation to be in ionization–recombination equilibrium and to have a low neutral fraction. Additionally, we characterize the ionization parameter as a function of the environment. Using a simplified model of the LyC radiation field, we produce synthetic maps of emission measure (EM). We show that the simplified model can be used to extract an estimate of the neutral fraction of the photoionized gas and mean free path of ionizing radiation from observed EM maps in galaxies. [ABSTRACT FROM AUTHOR]

We report the first detection of the metal-bearing molecules sodium sulfide and magnesium sulfide and the tentative detection of calcium monoxide in the interstellar medium toward the Galactic center molecular cloud G+0.693-0.027. The derived column densities are (5.0 ± 1.1) × 1010 cm−2, (6.0 ± 0.6) × 1010 cm−2, and (2.0 ± 0.5) × 1010 cm−2, respectively. This translates into fractional abundances with respect to H2 of (3.7 ± 1.0) × 10−13, (4.4 ± 0.8) × 10−13, and (1.5 ± 0.4) × 10−13, respectively. We have also searched for other Na-, Mg-, and Ca-bearing species toward this source but none of them have been detected and thus we provide upper limits for their abundances. We discuss the possible chemical routes involved in the formation of these molecules containing metals under interstellar conditions. Finally, we compare the ratio between sulfur-bearing and oxygen-bearing molecules with and without metals, finding that metal-bearing sulfur molecules are much more abundant than metal-bearing oxygen ones, in contrast with the general trend found in the ratios between other nonmetal-oxygen- and sulfur-bearing molecules. This further strengthens the idea that sulfur may be a little depleted in G+0.693-0.027 as a result of the low-velocity shocks present in this source sputtering large amounts of material from dust grains. [ABSTRACT FROM AUTHOR]

The cold neutral medium (CNM) is where neutral atomic hydrogen (H i) is converted into molecular clouds, so the structure and kinematics of the CNM are key drivers of galaxy evolution. Here we provide new constraints on the vertical distribution of the CNM using the recently developed kinematic_scaleheight software package and a large catalog of sensitive H i absorption observations. We estimate the thickness of the CNM in the solar neighborhood to be σ z ∼ 50–90 pc, assuming a Gaussian vertical distribution. This is a factor of ∼2 smaller than typically assumed, indicating that the thickness of the CNM in the solar neighborhood is similar to that found in the inner Galaxy, consistent with recent simulation results. If we consider only structures with H i optical depths τ > 0.1 or column densities N (H i) > 1019.5 cm−2, which recent work suggests are thresholds for molecule formation, we find σ z ∼ 50 pc. Meanwhile, for structures with τ < 0.1 or column densities N (H i) < 1019.5 cm−2, we find σ z ∼ 120 pc. These thicknesses are similar to those derived for the thin- and thick-disk molecular cloud populations traced by CO emission, possibly suggesting that cold H i and CO are well mixed. Approximately 20% of CNM structures are identified as outliers, with kinematics that are not well explained by Galactic rotation. We show that some of these CNM structures—perhaps representing intermediate-velocity clouds—are associated with the Local Bubble wall. We compare our results to recent observations and simulations, and we discuss their implications for the multiphase structure of the Milky Way's interstellar medium. [ABSTRACT FROM AUTHOR]

Context. The scientific community employs complicated multiphysics simulations to understand the physics in solar, stellar, and interstellar media. These must be tested against known solutions to ensure their validity. Several well-known tests exist, such as the Sod shock tube test. However, a test for nonlinear diffusivity is missing. This problem is highly relevant in the solar atmosphere, where various events release energy that subsequently diffuses by Spitzer thermal conductivity. Aims. The aim is to derive an analytical solution for nonlinear diffusivity in 1D, 2D, and 3D, which allows for a nonzero background value. The solution is used to design a test for numerical solvers and study Spitzer conductivity in the solar atmosphere. Methods. There exists an ideal solution assuming zero background value. We performed an analytical first-order perturbation of this solution. The first-order solution was first tested against a dedicated nonlinear diffusion solver, whereupon it was used to benchmark the single- and multifluid radiative magnetohydrodynamics code Ebysus, used to study the Sun. The theory and numerical modeling were used to investigate the role of Spitzer conductivity in the transport of energy released in a nanoflare. Results. The derived analytical solution models nonlinear diffusivity accurately within its region of validity and approximately beyond. Various numerical schemes implemented in the Ebysus code is found to model Spitzer conductivity correctly. The energy from a representative nanoflare is found to diffuse 9 Mm within the first second of its lifetime due to Spitzer conductivity alone, strongly dependent on the electron density. Conclusions. The analytical first-order solution is a step forward in ensuring the physical validity of intricate simulations of the Sun. Additionally, since the derivation and argumentation are general, they can easily be followed to treat other nonlinear diffusion problems. [ABSTRACT FROM AUTHOR]

Context. Observations of ionic, atomic, or molecular lines are performed to improve our understanding of the interstellar medium (ISM). However, the potential of a line to constrain the physical conditions of the ISM is difficult to assess quantitatively, because of the complexity of the ISM physics. The situation is even more complex when trying to assess which combinations of lines are the most useful. Therefore, observation campaigns usually try to observe as many lines as possible for as much time as possible. Aims. We have searched for a quantitative statistical criterion to evaluate the full constraining power of a (combination of) tracer(s) with respect to physical conditions. Our goal with such a criterion is twofold. First, we want to improve our understanding of the statistical relationships between ISM tracers and physical conditions. Secondly, by exploiting this criterion, we aim to propose a method that helps observers to make their observation proposals; for example, by choosing to observe the lines with the highest constraining power given limited resources and time. Methods. We propose an approach based on information theory, in particular the concepts of conditional differential entropy and mutual information. The best (combination of) tracer(s) is obtained by comparing the mutual information between a physical parameter and different sets of lines. The presented analysis is independent of the choice of the estimation algorithm (e.g., neural network or χ2 minimization). We applied this method to simulations of radio molecular lines emitted by a photodissociation region similar to the Horsehead Nebula. In this simulated data, we considered the noise properties of a state-of-the-art single dish telescope such as the IRAM 30m telescope. We searched for the best lines to constrain the visual extinction, AVtot, or the ultraviolet illumination field, G0. We ran this search for different gas regimes, namely translucent gas, filamentary gas, and dense cores. Results. The most informative lines change with the physical regime (e.g., cloud extinction). However, the determination of the optimal (combination of) line(s) to constrain a physical parameter such as the visual extinction depends not only on the radiative transfer of the lines and chemistry of the associated species, but also on the achieved mean signal-to-noise ratio. The short integration time of the CO isotopologue J = 1 − 0 lines already yields much information on the total column density for a large range of (AVtot, G0) space. The best set of lines to constrain the visual extinction does not necessarily combine the most informative individual lines. Precise constraints on the radiation field are more difficult to achieve with molecular lines. They require spectral lines emitted at the cloud surface (e.g., [CII] and [CI] lines). Conclusions. This approach allows one to better explore the knowledge provided by ISM codes, and to guide future observation campaigns. [ABSTRACT FROM AUTHOR]

Context. The appearance of the jets and lobes of some radio galaxies makes it difficult to assign them to a known class of objects. This is often due to the activity of the central engine and/or interactions with the environment, as well as projection effects. Aims. We analyse the radio data for an apparently asymmetric radio source 4C 70.19, which is associated with the giant elliptical galaxy NGC 6048. The source shows distorted radio jets and lobes, one of which bends by 180°. The aim of our study is to explain the nature of the observed distortions. Methods. We used LOFAR, Effelsberg, and VLA radio data in a wide range of frequencies. At high frequencies, we also used radio polarimetry to study the properties of the magnetic fields. Additionally, we made use of optical, infra-red, and X-ray data. Results. Polarisation data suggest shearing of the magnetic fields at points where the jets bend. The low-frequency LOFAR map at 145 MHz, as well as the sensitive single-dish Effelsberg map at 8.35 GHz, reveal previously undetected diffuse emission around the source. The rotation measure (RM) derived from the polarimetric data allowed us to estimate the density of the medium surrounding the source, which agrees with typical densities of the intergalactic medium or the outer parts of insterstellar halos. Conclusions. We propose that the southern jet is bent in the same manner as the northern one, but that it is inclined to the sky plane. Both these bends are likely caused by the orbital motion within the galaxy group, as well as interactions with the intergalactic medium. Our analyses suggest that, despite its complex morphology, 4C 70.19 seems to be intrinsically symmetric with a physical extent of up to 600 kpc, and that the diffuse emission detected in our high-sensitivity maps is related to radio plumes that are expanding behind the source. [ABSTRACT FROM AUTHOR]

We present new λrest = 77 μm dust continuum observations from the Atacama Large Millimeter/submillimeter Array of HZ10 (CRISTAL-22). This dusty main sequence galaxy at z = 5.66 was observed as part of the [CII] Resolved Ism in STar-forming Alma Large program (CRISTAL). The high angular resolution of the ALMA Band 7 and new Band 9 data (∼0′′​​.4) reveals the complex structure of HZ10, which comprises two main components (HZ10-C and HZ10-W), along with a bridge-like dusty emission between them (i.e., "the bridge"). Using a modified blackbody function to model the dust spectral energy distribution (SED), we constrained the physical conditions of the interstellar medium (ISM) and its variations among the different components identified in HZ10. We find that HZ10-W (the more UV-obscured component) has an SED dust temperature of TSED ∼ 51.2 ± 13.1 K; this was found to be ∼5 K higher (which is statistically insignificant; i.e., less than 1σ) than that of the central component and previous global estimations for HZ10. Our new ALMA data allow us to reduce the uncertainties of global TSED measurements by a factor of ∼2.3, compared to previous studies. The HZ10 components have [CII]-to-far-infrared (FIR) luminosity ratios and FIR surface densities values that are consistent with local starburst galaxies. However, HZ10-W shows a lower [CII]/FIR ratio compared to the other two components (albeit still within the uncertainties), which may suggest a harder radiation field destroying polycyclic aromatic hydrocarbon associated with [CII] emission (e.g., active galactic nuclei or young stellar populations). While HZ10-C appears to follow the tight IRX-βUV relation seen in local UV-selected starburst galaxies and high-z star-forming galaxies, we find that both HZ10-W and the bridge depart from this relation and are well described by dust-screen models with holes in front of a hard UV radiation field. This suggests that the UV emission, which is likely coming from young stellar populations, is strongly attenuated in the "dustier" components of the HZ10 system. [ABSTRACT FROM AUTHOR]

The extraplanar diffuse ionised gas (eDIG) is a key component for understanding the feedback processes that connect galactic discs and their halos. In this paper, we present the second study of the Bidimensional Exploration of the warm-Temperature Ionised Gas (BETIS) project, the aim of which is to explore the possible ionisation mechanisms and characteristics of the eDIG. We use a sample of eight edge-on galaxies observed with the Multi-Unit Spectroscopic Explorer (MUSE) integral field spectrograph (IFS) and apply the methodology developed in the first paper of the BETIS project for obtaining binned emission line maps. We find that the vertical and radial profiles of the [N II]/Hα, [S II]/Hα, [O III]/Hβ, and [O I]/Hα ratios depict a complex ionisation structure within galactic halos – which is influenced by the spatial distribution of H II regions across the galactic plane as observed from our line of sight–, with Lyman continuum photon leakage from OB associations constituting the main ionisation source. Moreover, the electron temperature and S+/S ionisation ratio also exhibit a dependency on the distribution of H II regions within the galactic discs. Our analysis excludes low-mass, hot, and evolved stars (HOLMES) as viable candidates for secondary ionisation sources to elucidate the unusual behaviour of the line ratios at greater distances from the galactic midplane. In contrast, we ascertain that shocks induced in the interstellar medium by star formation(SF)-related feedback mechanisms represent a promising secondary ionisation source of the eDIG. We present a suite of models integrating ionisation mechanisms arising from fast shocks and photoionisation associated with star formation. When applied to the classical Baldwin–Phillips–Terlevich (BPT) diagrams, these models reveal that the ionisation budget of the eDIG ranges from 20% to 50% across our sample, with local variations of up to 20% within individual galaxy halos. This contribution correlates with the presence of filaments and other structural components observed within galaxy halos. The presence of shocks is additionally supported by the observation of a high density of high [O I]/Hα ratios, which is characteristic of shock-compressed ionised gas, and is likely induced by feedback from regions of intense SF within the galactic disc. These results demonstrate consistency across all galaxies analysed in this sample. [ABSTRACT FROM AUTHOR]

The Milky Way's Central Molecular Zone (CMZ) differs dramatically from our local solar neighbourhood, both in the extreme interstellar medium conditions it exhibits (e.g. high gas, stellar, and feedback density) and in the strong dynamics at play (e.g. due to shear and gas influx along the bar). Consequently, it is likely that there are large-scale physical structures within the CMZ that cannot form elsewhere in the Milky Way. In this paper, we present new results from the Atacama Large Millimeter/submillimeter Array (ALMA) large programme ACES (ALMA CMZ Exploration Survey) and conduct a multi-wavelength and kinematic analysis to determine the origin of the M0.8–0.2 ring, a molecular cloud with a distinct ring-like morphology. We estimate the projected inner and outer radii of the M0.8–0.2 ring to be 79″ and 154″, respectively (3.1 pc and 6.1 pc at an assumed Galactic Centre distance of 8.2 kpc) and calculate a mean gas density >104 cm−3, a mass of ~106 M⊙, and an expansion speed of ~20 km s−1, resulting in a high estimated kinetic energy (>1051 erg) and momentum (>107 M⊙ km s−1). We discuss several possible causes for the existence and expansion of the structure, including stellar feedback and large-scale dynamics. We propose that the most likely cause of the M0.8–0.2 ring is a single high-energy hypernova explosion. To viably explain the observed morphology and kinematics, such an explosion would need to have taken place inside a dense, very massive molecular cloud, the remnants of which we now see as the M0.8–0.2 ring. In this case, the structure provides an extreme example of how supernovae can affect molecular clouds. [ABSTRACT FROM AUTHOR]

Context. Accretion onto supermassive black holes at close to the Eddington rate is expected to drive powerful winds, which have the potential to majorly influence the properties of the host galaxy. Theoretical models of such winds can simultaneously explain observational correlations between supermassive black holes and their host galaxies, such as the M − σ relation, and the powerful multi-phase outflows that are observed in a number of active galaxies. Analytic models developed to understand these processes usually assume simple galaxy properties, namely spherical symmetry and a smooth gas distribution with an adiabatic equation of state. However, the interstellar medium in real galaxies is clumpy and cooling is important, complicating the analysis. Aims. We wish to determine how gas turbulence, uneven density distribution, and cooling influence the development of active galactic nucleus (AGN) wind-driven outflows and their global properties on kiloparsec scales. Methods. We calculated a suite of idealised hydrodynamical simulations of AGN outflows designed to isolate the effects of turbulence and cooling, both separately and in combination. All simulations initially consisted of a 1 kpc gas shell with an AGN in the centre. We measured the main outflow parameters – the velocity, the mass outflow rate (Ṁout), and the momentum (ṗoutc/LAGN) and energy (Ėout/LAGN) loading factors – as the system evolves over 1.2 Myr and estimated plausible observationally derived values. Results. We find that adiabatic simulations approximately reproduce the analytical estimates of outflow properties independently of the presence or absence of turbulence and clumpiness. Cooling, on the other hand, has a significant effect, reducing the outflow energy rate by one to two orders of magnitude in the smooth simulations and by up to one order of magnitude in the turbulent ones. The interplay between cooling and turbulence depends on AGN luminosity: in Eddington-limited AGN, turbulence enhances the coupling between the AGN wind and the gas, while in lower-luminosity simulations, the opposite is true. This mainly occurs because dense gas clumps are resilient to low-luminosity AGN feedback but get driven away by high-luminosity AGN feedback. The overall properties of multi-phase outflowing gas in our simulations qualitatively agree with observations of multi-phase outflows, although there are some quantitative differences. We also find that using 'observable' outflow properties leads to their parameters being underestimated by a factor of a few compared with real values. Conclusions. We conclude that the AGN wind-driven outflow model is capable of reproducing realistic outflow properties in close-to-realistic galaxy setups and that the M − σ relation can be established without efficient cooling of the shocked AGN wind. Furthermore, we suggest ways to improve large-scale numerical simulations by accounting for the effects of AGN wind. [ABSTRACT FROM AUTHOR]

Context. We investigate active galactic nuclei (AGN) feeding through the molecular gas (CO(2−1) emission) properties of the local Seyfert 1 galaxy NGC 4593, using Atacama Large Millimeter Array (ALMA) observations and other multi-wavelength data. Aims. Our study aims to understand the interplay between the AGN and the interstellar medium (ISM) in this galaxy, examining the role of the AGN in steering gas dynamics within its host galaxy, evaluating the energy injected into the ISM, and determining whether gas is inflowing or outflowing from the galaxy. Methods. After reducing the ALMA CO(2−1) images, we employed two models, 3D-BAROLO and DISCFIT, to construct a disc model and fit its emission to the ALMA data. Additionally, we used photometric data to build a spectral energy distribution (SED) and apply the CIGALE code to derive key physical properties of the AGN and its host. Results. Our analysis reveals a complex interplay within NGC 4593, including a clear rotational pattern, the influence of a non-axisymmetric bar potential, and a central molecular zone (CMZ)-like ring. We observe an outflow of CO(2−1) gas along the minor axis, at a distance of ∼220 pc from the nucleus. The total molecular gas mass is estimated to be 1 − 5 × 108 M⊙, with non-circular motions contributing 10%. Our SED analysis indicates an AGN fraction of 0.88 and a star formation rate (SFR) of 0.42 M⊙ yr−1. Conclusions. These findings highlight the complex dynamics in the centre of NGC 4593, which are significantly influenced by the presence of the AGN. The overall physical properties of this system suggest that the AGN has a substantial impact on the evolution of NGC 4593. [ABSTRACT FROM AUTHOR]

Resolved studies of the correlation between the radio and far-infrared (FIR) emission from galaxies at different frequencies can unveil the interplay between star formation and the relativistic interstellar medium (ISM). Thanks to the LOFAR LoTSS observations combined with VLA, Herschel, and WISE data, we study the role of cosmic rays and magnetic fields in the radio–FIR correlation on scales of ≳200 pc in the nearby galaxy IC 342. The thermal emission traced by the 22 μm emission, constitutes about 6%, 13%, and 30% of the observed radio emission at 0.14, 1.4, and 4.8 GHz, respectively, in star-forming regions and less in other parts. The nonthermal spectral index becomes flatter at frequencies lower than 1.4 GHz (αn = −0.51 ± 0.09, Sν ∝ ναn) than between 1.4 and 4.8 GHz (αn = −1.06 ± 0.19) on average, and this flattening occurs not only in star-forming regions but also in the diffuse ISM. The radio–FIR correlation holds at all radio frequencies; however, it is tighter at higher radio frequencies. A multi-scale analysis shows that this correlation cannot be maintained on small scales due to diffusion of cosmic ray electrons (CREs). The correlation breaks at a larger scale (≃320 pc) at 0.14 GHz than at 1.4 GHz (≃200 pc), indicating that the CREs traced at lower frequencies have diffused a longer path in the ISM. We find that the energy index of CREs becomes flatter in star-forming regions, in agreement with previous studies. Cooling of CREs due to the magnetic field is evident globally only after compensating for the effect of star formation activity that both accelerates CREs and amplifies magnetic fields. Compared with other nearby galaxies, we show that the smallest scale of the radio–FIR correlation is proportional to the propagation length of the CREs on which the ordered magnetic field has an important effect. [ABSTRACT FROM AUTHOR]

We present JWST Near-Infrared Spectrograph (NIRSpec) observations in integral field spectroscopic (IFS) mode of the galaxy GS5001 at redshift z = 3.47, the central member of a candidate protocluster in the GOODS-S field. The data cover a field of view (FoV) of 4″ × 4″ (∼30 × 30 kpc2) and were obtained as part of the Galaxy Assembly with NIRSpec IFS (GA-NIFS) GTO programme. The observations include both high (R ∼ 2700) and low (R ∼ 100) spectral resolution data, spanning the rest-frame wavelength ranges 3700–6780 Å and 1300–11850 Å, respectively. These observations enable the detection and mapping of the main optical emission lines from [O II]λλ3726, 29 to [S III]λ9531. We analysed the spatially resolved ionised gas kinematics and interstellar medium properties, including obscuration, gas metallicity, excitation, ionisation parameter, and electron density. In addition to the main galaxy (GS5001), the NIRSpec FoV covers a close companion in the south, with three sub-structures with velocities blueshifted by ∼ − 150 km s−1 with respect to GS5001, and another source in the north redshifted by ∼200 km s−1. Optical line ratio diagnostics indicate star formation ionisation and electron densities of ∼500 cm−3 across all sources in the FoV. The gas-phase metallicity in the main galaxy is 12+log(O/H) = 8.45 ± 0.04, and slightly lower in the companions (12+log(O/H)  =  8.34 − 8.42), consistent with the mass-metallicity relation at z ∼ 3. We find peculiar line ratios (high log[N II]/Hα = [−0.45, −0.3], low log[O III]/Hβ = [0.06, 0.10]) in the northern part of GS5001. These could be attributed to either higher metallicity, or to shocks resulting from the interaction of the main galaxy with the northern source. We identify a spatially resolved outflow in the main galaxy, traced by a broad symmetric component in Hα and [O III], with an extension of about 3 kpc. We find maximum outflow velocities of ∼400 km s−1, an outflow mass of (1.7 ± 0.4)×108 M⊙, a mass outflow rate of 23 ± 5 M⊙ yr−1, and a mass loading factor of 0.23. These properties are compatible with star formation being the driver of the outflow. Our analysis of these JWST NIRSpec IFS data therefore provides valuable, unprecedented insights into the interplay between star formation, galactic outflows, and interactions in the core of a z ∼ 3.5 candidate protocluster. [ABSTRACT FROM AUTHOR]

Extended Lyα emission is routinely found around single quasars across cosmic time. However, few studies have investigated how such emission changes in fields with physically associated quasar pairs, which should reside in dense environments and are predicted to be linked through intergalactic filaments. We present VLT/MUSE snapshot observations (45 minutes/source) to unveil extended Lyα emission on scales of the circumgalactic medium (CGM) around the largest sample of physically associated quasar pairs to date, encompassing eight pairs (14 observed quasars) at z ∼ 3 with an i-band magnitude between 18 and 22.75, corresponding to absolute magnitudes Mi(z = 2) between −29.6 and −24.9. The pairs are either at close (∼50–100 kpc, five pairs) or wide (∼450–500 kpc, three pairs) angular separation and have velocity differences of Δv ≤ 2000 km s−1. We detected extended emission around 12 of the 14 targeted quasars and investigated the luminosity, size, kinematics, and morphology of these Lyα nebulae. On average, they span about 90 kpc and are 2.8 × 1043 erg s−1 bright. Irrespective of the quasars' projected distance, the nebulae often (∼45%) extend toward the other quasar in the pair, leading to asymmetric emission whose flux-weighted centroid is at an offset position from any quasar location. We show that large nebulae are preferentially aligned with the large-scale structure, as traced by the direction between the two quasars, and conclude that the cool gas (104 K) in the CGM traces well the direction of cosmic web filaments. Additionally, the radial profile of the Lyα surface brightness around quasar pairs can be described by a power law with a shallower slope (∼−1.6) with respect to single quasars (∼−2), indicative of increased CGM densities out to large radii and/or an enhanced contribution from the intergalactic medium (IGM) due to the dense environments expected around quasar pairs. The sample presented in this study contains excellent targets for ultra-deep observations to directly study filamentary IGM structures in emission. This work demonstrates that a large snapshot survey of quasar pairs will pave the way to direct statistical study of the IGM. [ABSTRACT FROM AUTHOR]

Context. Feedback from stars in the form of radiation, stellar winds, and supernovae is crucial to regulating the star formation activity of galaxies. Dwarf galaxies are especially susceptible to these processes, making them an ideal test bed for studying the effects of stellar feedback in detail. Recent numerical models have aimed to resolve the interstellar medium (ISM) in dwarf galaxies with a very high resolution of several solar masses. However, when it comes to modeling the radiative feedback from stars, many models opt for simplified approaches instead of explicitly solving radiative transfer (RT) because of the computational complexity involved. Aims. We introduce the Realistic ISM modeling in Galaxy Evolution and Lifecycles (RIGEL) model, a novel framework to self-consistently model the effects of stellar feedback in the multiphase ISM of dwarf galaxies with explicit RT on a star-by-star basis. Methods. The RIGEL model integrates detailed implementations of feedback from individual massive stars into the state-of-the-art radiation-hydrodynamics code, AREPO-RT. It forms individual massive stars from the resolved multiphase ISM by sampling the initial mass function and tracks their evolution individually. The lifetimes, photon production rates, mass-loss rates, and wind velocities of these stars are determined by their initial masses and metallicities based on a library that incorporates a variety of stellar models. The RT equations are solved explicitly in seven spectral bins accounting for the infrared to He II ionizing bands, using a moment-base scheme with the M1 closure relation. The thermochemistry model tracks the nonequilibrium H, He chemistry as well as the equilibrium abundance of C I, C II, O I, O II, and CO in the irradiated ISM to capture the thermodynamics of all ISM phases, from cold molecular gas to hot ionized gas. Results. We evaluated the performance of the RIGEL model using 1 M⊙ resolution simulations of isolated dwarf galaxies. We found that the star formation rate (SFR) and interstellar radiation field (ISRF) show strong positive correlations with the metallicity of the galaxy. Photoionization and photoheating can reduce the SFR by an order of magnitude by removing the available cold, dense gas fuel for star formation. The presence of ISRF also significantly changes the thermal structure of the ISM. Radiative feedback occurs immediately after the birth of massive stars and rapidly disperses the molecular clouds within 1 Myr. As a consequence, radiative feedback reduces the age spread of star clusters to less than 2 Myr, prohibits the formation of massive star clusters, and shapes the cluster initial mass function to a steep power-law form with a slope of ∼ − 2. The mass-loading factor (measured at z = 1 kpc) of the fiducial galaxy has a median of ηM ∼ 50, while turning off radiative feedback reduces this factor by an order of magnitude. Conclusions. We demonstrate that RIGEL effectively captures the nonlinear coupling of early radiative feedback and supernova feedback in the multiphase ISM of dwarf galaxies. This novel framework enables the utilization of a comprehensive stellar feedback and ISM model in cosmological simulations of dwarf galaxies and various galactic environments spanning a wide dynamic range in both space and time. [ABSTRACT FROM AUTHOR]

Context. Near- to mid-infrared observations (from Spitzer and JWST) have revealed a hidden population of galaxies at redshift z = 3 − 6 called optically dark objects, which are believed to be massive and dusty star-formers. They contribute substantially to the cosmic star-formation rate (SFR) density at z ∼ 4 − 5 (up to 30 − 40%). Aims. While optically dark sources are widely recognized as a significant component of the stellar mass function, the history of their stellar mass assembly (and the evolution of their interstellar medium) remains unexplored. However, they are thought to be the progenitors of the more massive early-type galaxies found in present-day groups and clusters. It is thus important to examine the possible connection between dark sources and merging events in order to understand the environment in which they live. Methods. Here, we report our search for close companions in a sample of 19 optically-dark objects identified in the SMACS0723 JWST deep field. They were selected in the NIRCam F444W band and undetected below 2 μm. We restricted our analysis to the reddest (i.e., F277W–F444W > 1.3) and brightest (F444W < 26 mag) objects. Results. We identified KLAMA, an optically dark source showing a very close companion (angular distance < 0.5″). The spatially resolved SED fitting procedure indicates that all components lying within 1.5″ of who is it the dark source are indeed at z ∼ 5.7. Tidal features (leading to a whale-shaped morphology) corroborate the hypothesis that KLAMA is the most massive (log(M⋆/M⊙) > 10.3) and dusty (AV ∼ 3 at the core) system of an ongoing merger with a mass ratio of ∼10. Thus, around ten similar merger events would be required to double the stellar mass of KLAMA. Merging systems with properties similar to KLAMA are identified in the SERRA simulations, allowing us to reconstruct their stellar-mass assembly history and predict their molecular gas properties (in particular, the [CII] emission for the simulated system). Conclusions. The discovery of mergers within dark galaxies at the end of the Epoch of Reionization highlights the importance of conducting a statistical search for additional candidates in deep NIRCam fields. Such research will aid in our understanding of the significance of merging processes during the obscured phase of stellar-mass accumulation. [ABSTRACT FROM AUTHOR]

Context. Magnetohydrodynamic (MHD) turbulence plays a critical role in many key astrophysical processes, such as star formation, acceleration of cosmic rays, and heat conduction. However, its properties are still poorly understood. Aims. In this work, we explore how to extract the intermittency of compressible MHD turbulence from synthetic and real observations. Methods. We used three statistical methods, namely the probability distribution function, kurtosis, and scaling exponent of the multi-order structure function, to reveal the intermittency of MHD turbulence. Results. Our numerical results demonstrate that: (1) the synchrotron polarization intensity statistics can be used to probe the intermittency of magnetic turbulence, by which we can distinguish different turbulence regimes; (2) the intermittency of MHD turbulence is dominated by the slow mode in the sub-Alfvénic turbulence regime; and (3) the Galactic interstellar medium (ISM) in the low latitude region corresponds to the sub-Alfvénic and supersonic turbulence regime. Conclusions.We have successfully measured the intermittency of the Galactic ISM from synthetic and realistic observations. [ABSTRACT FROM AUTHOR]

Context. S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of Lj ∼ 5 × 1040 erg s−1. This implies that the accretion power to the black hole should be super-Eddington, of the order of Lacc ∼ Lj. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of LX ≈ 1037 erg s−1. Aims. We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. Methods. We first analyze multi-epoch X-ray observations of S26 obtained with XMM-Newton and model the super-Eddington disk and its wind. We then develop a jet model and study the particle acceleration and radiative processes that occur in shocks generated near the base of the jet and in its terminal region. Results. We find that the discrepancy between the jet and the apparent disk luminosities in S26 is caused by the complete absorption of the disk radiation by the wind ejected from the super-Eddington disk. The nonthermal X-rays are produced near the base of the jet, and the thermal X-rays are emitted in the terminal regions. The radio emission observed with the Australia Telescope Compact Array can be explained as synchrotron radiation produced at the reverse shock in the lobes. We also find that S26 can accelerate protons to PeV energies in both the inner jet and the lobes. The ultra-high energy protons accelerated in the lobes of S26 are injected into the interstellar medium with a total power of ∼1036 erg s−1. Conclusions. We conclude that S26 is a super-Eddington microquasar with a dense disk-driven wind that obscures the X-ray emission from the inner disk, and that the supercritical nature of the system allows the acceleration of cosmic rays to PeV energies. [ABSTRACT FROM AUTHOR]

Outflows from galaxies that are driven by active galactic nuclei and star formation activity spread magnetic fields into the intergalactic medium. The importance of this process can be assessed using cosmological magnetohydrodynamical numerical modeling of the baryonic feedback on the large-scale structure, such as that of IllustrisTNG simulations. We use the Faraday rotation measure data of the LOFAR Two-Metre Sky Survey (LoTSS) to test the IllustrisTNG baryonic feedback model. We show that the IllustrisTNG overpredicts the root mean square of the residual rotation measure in LoTSS, which suggests that pollution of the intergalactic medium by magnetized outflows from galaxies is less important than the estimate from IllustrisTNG. This fact supports the hypothesis that the volume-filling large-scale magnetic fields in voids of the large-scale structure are of cosmological origin. [ABSTRACT FROM AUTHOR]

The molecular-to-atomic gas ratio is crucial to our understanding of the evolution of the interstellar medium (ISM) in galaxies. We investigated the balance between the atomic (ΣHI) and molecular gas (ΣH2) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for H I and CO, respectively). We defined the molecular gas ratio as Rmol = ΣH2/ΣHI and measured how Rmol depends on local conditions in the galaxy disks using multiwavelength observations. We find that, depending on the galaxy, H I is detected at > 3σ out to 20 − 120 kpc in galactocentric radius (rgal). The typical radius at which ΣHI reaches 1 M⊙ pc−2 is rH I ≈ 22 kpc, which corresponds to 1 − 3 times the optical radius (r25). We note that, Rmol correlates best with the dynamical equilibrium pressure, PDE, among potential drivers studied, with a median correlation coefficient of ⟨ρ⟩ = 0.89. Correlations between Rmol and the star formation rate surface density, total gas surface density, stellar surface density, metallicity, and ΣSFR/PDE (a proxy for the combined effect of the UV radiation field and number density) are present but somewhat weaker. Our results also show a direct correlation between PDE and ΣSFR, supporting self-regulation models. Quantitatively, we measured similar scalings as previous works, and attribute the modest differences that we do find to the effect of varying resolution and sensitivity. At rgal ≳ 0.4r25, atomic gas dominates over molecular gas among our studied galaxies, and at the balance of these two gas phases (Rmol = 1), we find that the baryon mass is dominated by stars, with Σ* > 5 Σgas. Our study constitutes an important step in the statistical investigation of how local galaxy properties (stellar mass, star formation rate, or morphology) impact the conversion from atomic to molecular gas in nearby galaxies. [ABSTRACT FROM AUTHOR]

Radiation from stars and active galactic nuclei (AGN) plays an important role in galaxy formation and evolution, and profoundly transforms the intergalactic, circumgalactic, and interstellar medium (IGM, CGM, and ISM). On-the-fly radiative transfer (RT) has started being incorporated in cosmological simulations, but the complex evolving radiation spectra are often crudely approximated with a small number of broad bands with piece-wise constant intensity and a fixed photo-ionisation cross-section. Such a treatment is unable to capture the changes to the spectrum as light is absorbed while it propagates through a medium with non-zero opacity. This can lead to large errors in photo-ionisation and heating rates. In this work we present a novel approach of discretising the radiation field at discrete photon energies, at the edges of the typically used photo-ionising bands, in order to capture the power-law slope of the radiation field. In combination with power-law approximations for the photo-ionisation cross-sections, this model allows us to self-consistently combine radiation from sources with different spectra and accurately follow the ionisation states of primordial and metal species through time. The method is implemented in GASOLINE2 in connection with TREVR2, a fast reverse ray tracing algorithm with 풪(Nactive log2 N) scaling. We compare our new piece-wise power-law reconstruction to the piece-wise constant method in calculating the primordial chemistry photo-ionisation and heating rates under an evolving UV background (UVB) and stellar spectrum, and find that our method reduces errors significantly, by up to two orders of magnitude in the case of HeII ionisation. We apply our new spectral reconstruction method in RT post-processing of a cosmological zoom-in simulation, MUGS2 g1536, including radiation from stars and a live UVB, and find a significant increase in total neutral hydrogen (HI) mass in the ISM and the CGM due to shielding of the UVB and a low escape fraction of the stellar radiation. This demonstrates the importance of RT and an accurate spectral approximation in simulating the CGM-galaxy ecosystem. [ABSTRACT FROM AUTHOR]

The extended Lyα haloes (LAHs) have been found to be prevalent around high-redshift star-forming galaxies. However, the origin of the LAHs is still a subject of debate. The spatially resolved analysis of Lyα profiles provides an important diagnostic. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Lyα emitters (LAEs) at redshifts of 3 < z < 4 in the MUSE Extremely Deep Field. Our analysis reveals that, with respect to the Lyα line of the target LAE, the peak of the Lyα line at large distances becomes increasingly more blueshifted up to a projected distance of 60 kpc (≈3× virial radius), with a velocity offset of ≈250 km/s. This trend is evident in both the mean and median stacks, suggesting that it is a general property of our LAE sample, which typically has a Lyα luminosity of ≈1041.1 erg s−1. However, due to the absence of systemic redshift data, it remains unclear whether the Lyα line peak at large projected distances is less redshifted compared to the inner regions or truly blueshifted with respect to the systemic velocity. We explore various scenarios to explain the large-scale kinematics of the Lyα line. [ABSTRACT FROM AUTHOR]

Context. Growing evidence suggests that protoplanetary discs may be influenced by late stage infall from the interstellar medium (ISM). It remains unclear the degree to which infall shapes disc populations at ages ≳1 Myr. Aims. We explored possible spatial correlations between stellar accretion rates in the Lupus star-forming region, which would support the hypothesis that infall can regulate stellar accretion. Methods. We considered both the 'clustered' stars towards the centre of Lupus 3, and the 'distributed' stars that are more sparsely distributed across the Lupus complex. We took the observed accretion rates in the literature and explore spatial correlations. In particular, we tested whether the clustered stars exhibit a radial gradient in normalised accretion rates, and whether the distributed stars have spatially correlated accretion rates. Results. We found statistically significant correlations for both the clustered and distributed samples. The clustered sample exhibits higher accretion rates in the central region, consistent with the expected Bondi-Hoyle-Lyttleton accretion rate. Stars that are spatially closer among the distributed population also exhibit more similar accretion rates. These results cannot be explained by the stellar mass distribution for either sample. Age gradients are disfavoured, though not discounted, because normalised disc dust masses are not spatially correlated across the region. Conclusions. Spatially correlated stellar accretion rates within the Lupus star-forming region argue in favour of an environmental influence on stellar accretion, possibly combined with internal processes in the inner disc. Refined age measurements and searches for evidence of infalling material are potential ways to further test this finding. [ABSTRACT FROM AUTHOR]

The Square Kilometre Array mid-frequency array will enable high-redshift detections of neutral hydrogen (H i) emission in galaxies, providing important constraints on the evolution of cold gas in galaxies over cosmic time. Strong gravitational lensing will push back the H i emission frontier towards cosmic noon (⁠|$z\sim 2$|⁠), as has been done for all prominent spectral lines in the interstellar medium of galaxies. Chakraborty & Roy report a |$z=1.3$| H i emission detection towards the well-modelled, galaxy-scale gravitational lens, SDSS J0826+5630. We carry out H i source modelling of the system and find that their claimed H i magnification, |$\mu _{\rm H\, {\small I}} = 29 \pm 6$|⁠ , requires an H i disc radius of |$\lesssim\!\! 1.5$| kpc, which implies an implausible mean H i surface mass density in excess of |$\Sigma _{\rm H\, {\small I}} > 2000 \ \rm{ M}_\odot \, \rm{pc}^{-2}$|⁠. This is several orders of magnitude above the highest measured peak values (⁠|$\Sigma _{\rm H\, {\small I}} \sim 10~{\rm M}_\odot \, {\rm pc}^{-2}$|⁠), above which H i is converted into molecular hydrogen. Our re-analysis requires this to be the highest H i mass galaxy known (⁠|$M_{\rm H\, {\small I}} \sim 10^{11} \ \mathrm{ M}_\odot$|⁠), as well as strongly lensed, the latter having a typical probability of the order of 1 in 103–104. We conclude that the claimed detection is spurious. [ABSTRACT FROM AUTHOR]

We explore the relations between the gas-phase metallicity radial profiles (few hundred inner parsec) and multiple galaxy properties for 15 Seyfert galaxies from the active galactic nuclei integral field spectroscopy sample using optical integral field unit observations from Gemini Multi-Object Spectrograph (GMOS) and Multi Unit Spectroscopic Explorer (MUSE) processed archival data. The data were selected at |$z \lesssim 0.013$| within black hole mass range |$\left[6\lt \log \left(M_{\rm BH}/{\rm M_\odot } \right)\lt 9\right]$| with moderate 14–150-keV X-ray luminosities |$\left[42\, \lesssim \, \log L_X (\rm erg\, s^{-1})\, \lesssim \, 44\right]$|⁠. We estimated the gas-phase metallicity using the strong-line methods and found mean values for the oxygen-dependent (⁠|$Z \sim 0.75\, \mathrm{Z}_\odot$|⁠) and nitrogen-dependent (⁠|$Z \sim 1.14\, \mathrm{Z}_\odot$|⁠) calibrations. These estimates show excellent agreement with |$\Delta Z \approx 0.19$| and 0.18 dex between the mean values from the two strong-line calibrations for GMOS and MUSE, respectively, consistent with the order of metallicity uncertainty via the strong-line methods. We contend that our findings align with a scenario wherein local Seyferts have undergone seamless gas accretion histories, resulting in positive metallicity profiles over an extended period of time, thereby providing insights into galaxy evolution and the chemical enrichment or depletion of the universe. Additionally, we argue that metal-poor gas inflow from the local interstellar medium and accreted through the circumgalactic medium on to the galaxy systems regulates the star formation processes by diluting their central metallicity and inverting their gas-phase metallicity radial gradients, producing a more prominent anticorrelation between gas-phase metallicity and Eddington ratio. [ABSTRACT FROM AUTHOR]

Cosmological information is usually extracted from the Lyman- |$\alpha$| (Ly  |$\alpha$|⁠) forest correlations using only either large-scale information interpreted through linear theory or using small-scale information interpreted by means of expensive hydrodynamical simulations. A complete cosmological interpretation of the three-dimensional (3D) correlations at all measurable scales is challenged by the need of more realistic models including the complex growth of non-linear small scales that can only be studied within large hydrodynamical simulations. Past works were often limited by the trade-off between the simulated cosmological volume and the resolution of the low-density intergalactic medium from which the Ly  |$\alpha$| signal originates. We conduct a suite of hydrodynamical simulations of the intergalactic medium, including one of the largest Ly  |$\alpha$| simulations ever performed in terms of volume (640 |$\, h^{-1}\, \mathrm{Mpc}$|⁠), alongside simulations in smaller volumes with resolutions up to 25 |$\, h^{-1}\, \mathrm{kpc}$|⁠ , which will be further improved to show resolution convergence in future studies. We compare the 3D Ly  |$\alpha$| power spectra (⁠|$P_{3\mathrm{D},\alpha }$|⁠) predicted by those simulations to different non-linear models. The inferred Ly  |$\alpha$| bias and redshift space distortion parameters, |$b_\alpha$| and |$\beta _\alpha$| are in remarkable agreement with those measured in SDSS (Sloan Digital Sky Survey) and DESI (Dark Energy Spectroscopic Instrument) data. We find that, contrary to intuition, the convergence of large-scale modes of the |$P_{3\mathrm{D},\alpha }$|⁠ , which determines |$\beta _\alpha$|⁠ , is primarily influenced by the resolution of the simulation box through mode coupling, rather than the box size itself. Finally, we study the Baryon Acoustic Oscillation (BAO) signal encoded in |$P_{3\mathrm{D},\alpha }$|⁠. For the first time with a hydrodynamical simulation, we clearly detect the BAO signal; however, we only marginally detect its damping, associated with the non-linear growth of the structures. [ABSTRACT FROM AUTHOR]