Search

We present a reconstruction of the line-of-sight motions of the local interstellar medium (ISM) based on the combination of a state-of-the-art model of the three-dimensional dust density distribution within 1.25 kpc from the Sun and the HI and CO line emission within Galactic latitudes $|b| < 5^{\circ}$. We use the histogram of oriented gradient (HOG) method to match the plane-of-the-sky 3D dust distribution across distances with the atomic and molecular line emission. We identify a significant correlation between the 3D dust model and the line emission. We employ this correlation to assign line-of-sight velocities to the dust across density channels and produce a face-on map of the local ISM radial motions with respect to the local standard of rest (LSR). We find that most of the material in the 3D dust model follows the large-scale pattern of Galactic rotation; however, we also report local departures from the rotation pattern with standard deviations of 12.1 and 6.1 km/s for the HI and CO line emission, respectively. The mean kinetic energy densities of these streaming motions are around 0.68 and 0.18 eV/cm$^{3}$ from either gas tracer. Assuming homogeneity and isotropy in the velocity field, these values are within a factor of a few of the total kinetic energy density. These kinetic energy values are roughly comparable to other energy densities, thus confirming the near-equipartition introduced by the feedback loops connecting the physical processes in the local ISM. Yet, we find energy and momentum overdensities of around a factor of ten concentrated in the Radcliffe Wave, the Split, and other local density structures. Although we do not find evidence of the Local Spiral Arm impact in these overdensities, their distribution suggests the influence of other large-scale effects that, in addition to supernova feedback, shape the energy distribution in the Solar neighborhood., Comment: 32 pages, 43 figures, submitted to Astronomy & Astrophysics (15NOV2024)

SPIDER is a balloon-borne instrument designed to map the cosmic microwave background at degree-angular scales in the presence of Galactic foregrounds. SPIDER has mapped a large sky area in the Southern Hemisphere using more than 2000 transition-edge sensors (TESs) during two NASA Long Duration Balloon flights above the Antarctic continent. During its first flight in January 2015, SPIDER observed in the 95 GHz and 150 GHz frequency bands, setting constraints on the B-mode signature of primordial gravitational waves. Its second flight in the 2022-23 season added new receivers at 280 GHz, each using an array of TESs coupled to the sky through feedhorns formed from stacks of silicon wafers. These receivers are optimized to produce deep maps of polarized Galactic dust emission over a large sky area, providing a unique data set with lasting value to the field. In this work, we describe the instrument's performance during SPIDER's second flight., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, JATIS

We use data from the ALMA Evolutionary Study of High Mass Protocluster Formation in the Galaxy (ALMAGAL) survey to study 100 ALMAGAL regions at $\sim$ 1 arsecond resolution located between $\sim$ 2 and 6 kpc distance. Using ALMAGAL $\sim$ 1.3mm line and continuum data we estimate flow rates onto individual cores. We focus specifically on flow rates along filamentary structures associated with these cores. Our primary analysis is centered around position velocity cuts in H$_2$CO (3$_{0,3}$ - 2$_{0,2}$) which allow us to measure the velocity fields, surrounding these cores. Combining this work with column density estimates we derive the flow rates along the extended filamentary structures associated with cores in these regions. We select a sample of 100 ALMAGAL regions covering four evolutionary stages from quiescent to protostellar, Young Stellar Objects (YSOs), and HII regions (25 each). Using dendrogram and line analysis, we identify a final sample of 182 cores in 87 regions. In this paper, we present 728 flow rates for our sample (4 per core), analysed in the context of evolutionary stage, distance from the core, and core mass. On average, for the whole sample, we derive flow rates on the order of $\sim$10$^{-4}$ M$_{sun}$yr$^{-1}$ with estimated uncertainties of $\pm$50%. We see increasing differences in the values among evolutionary stages, most notably between the less evolved (quiescent/protostellar) and more evolved (YSO/HII region) sources. We also see an increasing trend as we move further away from the centre of these cores. We also find a clear relationship between the flow rates and core masses $\sim$M$^{2/3}$ which is in line with the result expected from the tidal-lobe accretion mechanism. Overall, we see increasing trends in the relationships between the flow rate and the three investigated parameters; evolutionary stage, distance from the core, and core mass., Comment: 11 pages, 11 figures, accepted for publication in A&A

Using data from the first flight of SPIDER and from Planck HFI, we probe the properties of polarized emission from interstellar dust in the SPIDER observing region. Component separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses spanning the full SPIDER region demonstrate that i) the spectral energy distribution of diffuse Galactic dust emission is broadly consistent with a modified-blackbody (MBB) model with a spectral index of $\beta_\mathrm{d}=1.45\pm0.05$ $(1.47\pm0.06)$ for $E$ ($B$)-mode polarization, slightly lower than that reported by Planck for the full sky; ii) its angular power spectrum is broadly consistent with a power law; and iii) there is no significant detection of line-of-sight decorrelation of the astrophysical polarization. The size of the SPIDER region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature $T_\mathrm{d}=19.6$ K, an analysis of two independent sub-regions of that field results in inferred values of $\beta_\mathrm{d}=1.52\pm0.06$ and $\beta_\mathrm{d}=1.09\pm0.09$, which are inconsistent at the $3.9\,\sigma$ level. Furthermore, a joint analysis of SPIDER and Planck 217 and 353 GHz data within a subset of the SPIDER region is inconsistent with a simple MBB at more than $3\,\sigma$, assuming a common morphology of polarized dust emission over the full range of frequencies. These modeling uncertainties have a small--but non-negligible--impact on limits on the cosmological tensor-to-scalar ratio derived from the \spider dataset. The fidelity of the component separation approaches of future CMB polarization experiments may thus have a significant impact on their constraining power., Comment: 21 pages, 15 figures

We investigate the geometry of the magnetic field towards the Radcliffe Wave, a coherent 3-kpc-long part of the nearby Local Arm recently discovered via three-dimensional dust mapping. We use archival stellar polarization in the optical and new measurements in the near-infrared to trace the magnetic field as projected on the plane of the sky. Our new observations cover the portion of the structure that is closest to the Sun, between Galactic longitudes of 122$^\circ$ and 188$^\circ$. The polarization angles of stars immediately background to the Radcliffe Wave appear to be aligned with the structure as projected on the plane of the sky. The observed magnetic field configuration is inclined with respect to the Galactic disk at an angle of 18$^\circ$. This departure from a geometry parallel to the plane of the Galaxy is contrary to previous constraints from more distant stars and polarized dust emission. We confirm that the polarization angle of stars at larger distances shows a mean orientation parallel to the Galactic disk. We discuss implications of the observed morphology of the magnetic field for models of the large-scale Galactic magnetic field, as well as formation scenarios for the Radcliffe Wave itself., Comment: Submitted to A&A

The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $\sigma_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $\sigma_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $\mu$G. Finally, we discuss the implications of our results within the observational context in the CMZ.

Targeted therapeutic interventions utilizing low-inten\-sity ultrasound (LIUS) exhibit substantial potential for hindering the proliferation of cancer stem cells. This investigation introduces a multiscale model and computational framework to comprehensively explore the therapeutic LIUS on poroelastic tumor dynamics, thereby unraveling the intricacies of mechanotransduction mechanisms at play. Our model includes both macroscopic timescales encompassing days and rapid timescales spanning from microseconds to seconds, facilitating an in-depth comprehension of tumor behavior. We unveil the discerning suppression or reorientation of cancer cell proliferation and migration, enhancing a notable redistribution of cellular phases and stresses within the tumor microenvironment. Our findings defy existing paradigms by elucidating the impact of LIUS on cancer stem cell behavior. This endeavor advances our fundamental understanding of mechanotransduction phenomena in the context of LIUS therapy, thus underscoring its promising as a targeted therapeutic modality for cancer treatment. Furthermore, our results make a substantial contribution to the broader scientific community by shedding light on the intricate interplay between mechanical forces, cellular responses, and the spatiotemporal evolution of tumors. These insights hold the promising to promote a new perspective for the future development of pioneering and highly efficacious therapeutic strategies for combating cancer in a personalized manner.

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution data from The HI/OH/Recombination line survey of the Milky Way (THOR), we aim to constrain the kinematic and physical properties of the cold atomic hydrogen gas phase toward the inner Galactic plane. HI self-absorption (HISA) has proven to be a viable method to detect cold atomic hydrogen clouds in the Galactic plane. With the help of a newly developed self-absorption extraction routine (astroSABER), we build upon previous case studies to identify HI self-absorption toward a sample of Giant Molecular Filaments (GMFs). We find the cold atomic gas to be spatially correlated with the molecular gas on a global scale. The column densities of the cold atomic gas traced by HISA are usually of the order of $10^{20}\rm\,cm^{-2}$ whereas those of molecular hydrogen traced by $\rm^{13}CO$ are at least an order of magnitude higher. The HISA column densities are attributed to a cold gas component that accounts for a fraction of $\sim$5% of the total atomic gas budget within the clouds. The HISA column density distributions show pronounced log-normal shapes that are broader than those traced by HI emission. The cold atomic gas is found to be moderately supersonic with Mach numbers of a $\sim$few. In contrast, highly supersonic dynamics drive the molecular gas within most filaments., Comment: 41 pages, 28 figures, accepted for publication in A&A

Many experimental quantities show a power-law distribution $p(x)\propto x^{-\alpha}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $\alpha$ and by the extremes $x_\text{min}$ $x_\text{max}$ where the distribution extends. There are no mathematical tools that derive the three unknowns at the same time. In general, one estimates a set of $\alpha$ corresponding to different guesses of $x_\text{min}$ $x_\text{max}$. Then, the best set of values describing the observed data is selected a posteriori. In this paper, we present a tool that finds contextually the three parameters based on simple assumptions on how the observed values $x_i$ populate the unknown range between $x_\text{min}$ and $x_\text{max}$ for a given $\alpha$. Our tool, freely downloadable, finds the best values through a non-linear least-squares fit. We compare our technique with the maximum likelihood estimators for power-law distributions, both truncated and not. Through simulated data, we show for each method the reliability of the computed parameters as a function of the number $N$ of data in the sample. We then apply our method to observed data to derive: i) the slope of the core mass function in the Perseus star-forming region, finding two power-law distributions: $\alpha=2.576$ between $1.06\,M_{\sun}$ and $3.35\,M_{\sun}$, $\alpha=3.39$ between $3.48\,M_{\sun}$ and $33.4\,M_{\sun}$; ii) the slope of the $\gamma$-ray spectrum of the blazar J0011.4+0057, extracted from the Fermi-LAT archive. For the latter case, we derive $\alpha=2.89$ between 1,484~MeV and 28.7~GeV; then we derive the time-resolved slopes using subsets 200 photons each., Comment: Corrected few typos found at proof-reading stage. The only important modification is in Table 4 where "x_M not constrained" is now "x_M=40"

We present a comparison of the Milky Way's star formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $\Sigma_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method determines $\Sigma_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun. We find an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods for computing the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10\,Myr., Comment: Accepted for publication at Astronomy & Astrophysics (11SEP2023). 8 pages, 7 figures

We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) method to match the morphology in a 3D reconstruction of the dust density (3D dust) and the distribution of the gas tracers' emission. The result of the HOG analysis is a map of the relationship between the distances and radial velocities. The HOG comparison between the 3D dust and the HI emission indicates a morphological match at the distance of Taurus but an anti-correlation between the dust density and the HI emission, which uncovers a significant amount of cold HI within the Taurus MC. The HOG between the 3D dust and $^{12}$CO reveals a pattern in radial velocities and distances that is consistent with converging motions of the gas in the Taurus MC, with the near side of the cloud moving at higher velocities and the far side moving at lower velocities. This convergence of flows is likely triggered by the large-scale gas compression caused by the interaction of the Local Bubble and the Per-Tau shell, with Taurus lying at the intersection of the two bubble surfaces., Comment: 18 pages, 21 figures. Accepted for publication in Astronomy & Astrophysics (22MAY2023)

High-spatial-resolution HI observations have led to the realisation that the nearby (within few hundreds of parsecs) Galactic atomic filamentary structures are aligned with the ambient magnetic field. Enabled by the high quality data from the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope for the Galactic ASKAP HI (GASKAP-HI) survey, we investigate the potential magnetic alignment of the $\gtrsim 10\,{\rm pc}$-scale HI filaments in the Small Magellanic Cloud (SMC). Using the Rolling Hough Transform (RHT) technique that automatically identifies filamentary structures, combined with our newly devised ray-tracing algorithm that compares the HI and starlight polarisation data, we find that the HI filaments in the northeastern end of the SMC main body ("Bar" region) and the transition area between the main body and the tidal feature ("Wing" region) appear preferentially aligned with the magnetic field traced by starlight polarisation. Meanwhile, the remaining SMC volume lacks starlight polarisation data of sufficient quality to draw any conclusions. This suggests for the first time that filamentary HI structures can be magnetically aligned across a large spatial volume ($\gtrsim\,{\rm kpc}$) outside of the Milky Way. In addition, we generate maps of the preferred orientation of HI filaments throughout the entire SMC, revealing the highly complex gaseous structures of the galaxy likely shaped by a combination of the intrinsic internal gas dynamics, tidal interactions, and star formation feedback processes. These maps can further be compared with future measurements of the magnetic structures in other regions of the SMC., Comment: 24 pages, MNRAS accepted

We present a study of synthetic observations of polarized dust emission at 353 GHz as seen by an observer within a cavity in the interstellar medium (ISM). The cavity is selected from a magnetohydrodynamic simulation of the local ISM with time-dependent chemistry, star formation, and stellar feedback in form of supernova explosions with physical properties comparable to the Local Bubble ones. We find that the local density enhancement together with the coherent magnetic field in the cavity walls makes the selected candidate a translucent polarization filter to the emission coming from beyond its domains. This underlines the importance of studying the Local Bubble in further detail. The magnetic field lines inferred from synthetic dust polarization data are qualitatively in agreement with the all-sky maps of polarized emission at 353 GHz from the Planck satellite in the latitudes interval 15deg <= |b| <= 65deg. As our numerical simulation allows us to track the Galactic midplane only out to distances of 250 pc, we exclude the region |b|<15deg from our analysis. At large Galactic latitudes, our model exhibits a high degree of small-scale structures. On the contrary, the observed polarization pattern around the Galactic poles is relatively coherent and regular, and we argue that the global toroidal magnetic field of the Milky Way is important for explaining the data at |b| > 65deg. We show that from our synthetic polarization maps, it is difficult to distinguish between an open and a closed Galactic cap using the inferred magnetic field morphology alone.

We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr$^{-1}$, of which $1.7 \pm 0.6$~M$_{\odot}$~yr$^{-1}$ coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius $R_\mathrm{gal} \sim 5$~kpc, followed by an exponential decrease as $\log(\Sigma_\mathrm{SFR}/[\mathrm{M}_\odot~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}])=a\,R_\mathrm{gal}/[\mathrm{kpc}]+b $, with $a=-0.28 \pm 0.01$. In this regard, the fraction of SFR produced within and outside the Solar circle is 84\% and 16\%, respectively; the fraction corresponding to the far outer Galaxy ($R_\mathrm{gal} > 13.5$~kpc) is only 1\%. We also find that, for $R_\mathrm{gal}>3$~kpc, our data follow a power law as a function of density, similarly to the Kennicutt-Schmidt relation. Finally, we compare the distribution of the SFR density across the face-on Galactic plane and those of median parameters, such as temperature, luminosity/mass ratio and bolometric temperature, describing the evolutionary stage of Hi-GAL clumps. We found no clear correlation between the SFR and the clump evolutionary stage., Comment: To be published in ApJ

We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher column densities. For diffuse, quiescent, and low-mass SF clouds, we propose that the two log-normals arise from the atomic and molecular phase, respectively. For massive clouds, we suggest that the first log-normal is built up by turbulently mixed H2 and the second one by compressed (via stellar feedback) molecular gas. Nearly all clouds have two PLTs with slopes consistent with self-gravity, where the second one can be flatter or steeper than the first one. A flatter PLT could be caused by stellar feedback or other physical processes that slow down collapse and reduce the flow of mass toward higher densities. The steeper slope could arise if the magnetic field is oriented perpendicular to the LOS column density distribution. The first deviation point (DP), where the N-PDF turns from log-normal into a PLT, shows a clustering around values of a visual extinction of AV (DP1) around 2-5. The second DP, which defines the break between the two PLTs, varies strongly. Using the Delta-variance, we observe that the AV value, where the slope changes between the first and second PLT, increases with the characteristic size scale in the variance spectrum. We conclude that at low column densities, atomic and molecular gas is turbulently mixed, while at high column densities, the gas is fully molecular and dominated by self-gravity. The best fitting model N-PDFs of molecular clouds is thus one with log-normal low column density distributions, followed by one or two PLTs., Comment: Astronomy and Astrophysics in presss

Experiments deploying large arrays of transition-edge sensors (TESs) often require a robust method to monitor gain variations with minimal loss of observing time. We propose a sensitive and non-intrusive method for monitoring variations in TES responsivity using small square waves applied to the TES bias. We construct an estimator for a TES's small-signal power response from its electrical response that is exact in the limit of strong electrothermal feedback. We discuss the application and validation of this method using flight data from SPIDER, a balloon-borne telescope that observes the polarization of the cosmic microwave background with more than 2000 TESs. This method may prove useful for future balloon- and space-based instruments, where observing time and ground control bandwidth are limited., Comment: 7 pages, 3 figures; Proceedings of the 19th International Workshop on Low Temperature Detectors (LTD19); Minor updates to match published version

Modern CMB analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-$C_\ell$ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, $F_{\ell}$, that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, $J_{\ell\ell'}$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the SPIDER balloon-borne CMB experiment., Comment: 14 pages, 7 figures; updated to match published version

The atomic phase of the interstellar medium plays a key role in the formation process of molecular clouds. Due to the line-of-sight confusion in the Galactic plane that is associated with its ubiquity, atomic hydrogen emission has been challenging to study. Employing the high-angular resolution data from the THOR survey, we identify one of the largest, coherent, mostly atomic HI filaments in the Milky Way at the line-of-sight velocities around -54 km/s. The giant atomic filament "Maggie", with a total length of 1.2 kpc, is not detected in most other tracers, and does not show signs of active star formation. At a kinematic distance of 17 kpc, Maggie is situated below (by 500 pc) but parallel to the Galactic HI disk and is trailing the predicted location of the Outer Arm by 5-10 km/s in longitude-velocity space. The centroid velocity exhibits a smooth gradient of less than $\pm$3 km/s /10 pc and a coherent structure to within $\pm$6 km/s. The line widths of 10 km/s along the spine of the filament are dominated by non-thermal effects. After correcting for optical depth effects, the mass of Maggie's dense spine is estimated to be $7.2\times10^5\,M_{\odot}$. The mean number density of the filament is 4$\rm\,cm^{-3}$, which is best explained by the filament being a mix of cold and warm neutral gas. In contrast to molecular filaments, the turbulent Mach number and velocity structure function suggest that Maggie is driven by transonic to moderately supersonic velocities that are likely associated with the Galactic potential rather than being subject to the effects of self-gravity or stellar feedback. The column density PDF displays a log-normal shape around a mean of $N_{\rm HI} = 4.8\times 10^{20}\rm\,cm^{-2}$, thus reflecting the absence of dominating effects of gravitational contraction., Comment: 19 pages, 17 figures, accepted for publication in A&A

G0.253+0.016, commonly referred to as "the Brick" and located within the Central Molecular Zone, is one of the densest ($\approx10^{3-4}$ cm$^{-3}$) molecular clouds in the Galaxy to lack signatures of widespread star formation. We set out to constrain the origins of an arc-shaped molecular line emission feature located within the cloud. We determine that the arc, centred on $\{l_{0},b_{0}\}=\{0.248^{\circ}, 0.18^{\circ}\}$, has a radius of $1.3$ pc and kinematics indicative of the presence of a shell expanding at $5.2^{+2.7}_{-1.9}$ km s$^{-1}$. Extended radio continuum emission fills the arc cavity and recombination line emission peaks at a similar velocity to the arc, implying that the molecular and ionised gas are physically related. The inferred Lyman continuum photon rate is $N_{\rm LyC}=10^{46.0}-10^{47.9}$ photons s$^{-1}$, consistent with a star of spectral type B1-O8.5, corresponding to a mass of $\approx12-20$ M$_{\odot}$. We explore two scenarios for the origin of the arc: i) a partial shell swept up by the wind of an interloper high-mass star; ii) a partial shell swept up by stellar feedback resulting from in-situ star formation. We favour the latter scenario, finding reasonable (factor of a few) agreement between its morphology, dynamics, and energetics and those predicted for an expanding bubble driven by the wind from a high-mass star. The immediate implication is that G0.253+0.016 may not be as quiescent as is commonly accepted. We speculate that the cloud may have produced a $\lesssim10^{3}$ M$_{\odot}$ star cluster $\gtrsim0.4$ Myr ago, and demonstrate that the high-extinction and stellar crowding observed towards G0.253+0.016 may help to obscure such a star cluster from detection., Comment: 17 pages, 7 figures. Accepted for publication in MNRAS (October 15, 2021)

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

We present a study of the cold atomic hydrogen (HI) content of molecular clouds simulated within the SILCC-Zoom project. We produce synthetic observations of HI at 21 cm including HI self-absorption (HISA) and observational effects. We find that HI column densities, $N_\textrm{HI}$, of $\gtrsim$10$^{22}$ cm$^{-2}$ are frequently reached in molecular clouds with HI temperatures as low as $\sim$10 K. Hence, HISA observations assuming a fixed HI temperature tend to underestimate the amount of cold HI in molecular clouds by a factor of 3 - 10 and produce an artificial upper limit of $N_\textrm{HI}$ around 10$^{21}$ cm$^{-2}$. We thus argue that the cold HI mass in clouds could be a factor of a few higher than previously estimated. Also $N_\textrm{HI}$-PDFs obtained from HISA observations might be subject to observational biases and should be considered with caution. The underestimation of cold HI in HISA observations is due to both the large HI temperature variations and the effect of noise in regions of high optical depth. We find optical depths of cold HI around 1 - 10 making optical depth corrections essential. We show that the high HI column densities ($\gtrsim$10$^{22}$ cm$^{-2}$) can in parts be attributed to the occurrence of up to 10 individual HI-H$_2$ transitions along the line of sight. This is also reflected in the spectra, necessitating Gaussian decomposition algorithms for their analysis. However, also for a single HI-H$_2$ transition, $N_\textrm{HI}$ frequently exceeds 10$^{21}$ cm$^{-2}$, challenging 1D, semi-analytical models. This is due to non-equilibrium chemistry effects and the fact that HI-H$_2$ transition regions usually do not possess a 1-dimensional geometry. Finally, we show that the HI gas is moderately supersonic with Mach numbers of a few. The corresponding non-thermal velocity dispersion can be determined via HISA observations within a factor of $\sim$2., Comment: 20 pages including appendix, 24 figures, accepted for publication in MNRAS

We present a statistical study of the filamentary structure orientation in the CO emission observations obtained in the Milky Way Imaging Scroll Painting (MWISP) survey in the range $25.8\deg < l < 49.7\deg$, $|b| \leq 1.25\deg$, and $-100 < v_{\rm LSR} < 135$ km/s. We found that most of the filamentary structures in the $^{12}$CO and $^{13}$CO emission do not show a global preferential orientation either parallel or perpendicular to the Galactic plane. However, we found ranges in Galactic longitude and radial velocity where the $^{12}$CO and $^{13}$CO filamentary structures are parallel to the Galactic plane. These preferential orientations are different from those found for the HI emission. We consider this an indication that the molecular structures do not simply inherit these properties from parental atomic clouds. Instead, they are shaped by local physical conditions, such as stellar feedback, magnetic fields, and Galactic spiral shocks., Comment: Accepted for publication in Astronomy & Astrophysics (7 pages, 7 figures). arXiv admin note: text overlap with arXiv:2007.07285

Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation., Comment: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set., Comment: 18 pages, 11 figures

We present the first linear polarization measurements from the 2015 long-duration balloon flight of SPIDER, an experiment designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. Results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance; Galactic synchrotron emission is found to be negligible in the SPIDER bands. We employ two independent foreground-removal techniques in order to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, employing a joint analysis of SPIDER and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and $r$ parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman-Cousins and Bayesian constructions, finding $r<0.11$ and $r<0.19$, respectively. Roughly half the uncertainty in $r$ derives from noise associated with the template subtraction. New data at 280 GHz from SPIDER's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission., Comment: 29 pages, 13 figures

In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for the second Antarctic flight will incorporate three new 280 GHz receivers alongside three refurbished 95- and 150 GHz receivers from Spider's first flight. In this work we discuss the design and characterization of these new receivers, which employ over 1500 feedhorn-coupled transition-edge sensors. We describe pre-flight laboratory measurements of detector properties, and the optical performance of completed receivers. These receivers will map a wide area of the sky at 280 GHz, providing new information on polarized Galactic dust emission that will help to separate it from the cosmological signal., Comment: 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how much magnetic fields shape the star formation process or whether gravity is largely dominating. Here, we show that for the high-mass star-forming region G327.3 the magnetic field morphology appears to be dominantly shaped by the gravitational contraction of the central massive gas core where the star formation proceeds. We find that in the outer parts of the region, the magnetic field is directed toward the gravitational center of the region. Filamentary structures feeding the central core exhibit U-shaped magnetic field morphologies directed toward the gravitational center as well, again showing the gravitational drag toward the center. The inner part then shows rotational signatures, potentially associated with an embedded disk, and there the magnetic field morphology appears to be rotationally dominated. Hence, our results demonstrate that for this region gravity and rotation are dominating the dynamics and shaping the magnetic field morphology., Comment: 10 pages, 4 figures, accepted for the Astrophysical Journal, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

We report the first map of large-scale (10 pc in length) emission of millimeter-wavelength hydrogen recombination lines (mm-RRLs) toward the giant H II region around the W43-Main young massive star cluster (YMC). Our mm-RRL data come from the IRAM 30 m telescope and are analyzed together with radio continuum and cm-RRL data from the Karl G. Jansky Very Large Array and HCO$^{+}$ 1-0 line emission data from the IRAM 30 m. The mm-RRLs reveal an expanding wind-blown ionized gas shell with an electron density ~70-1500 cm$^{-3}$ driven by the WR/OB cluster, which produces a total Ly$\alpha$ photon flux of 1.5 x 10$^{50}$ s$^{-1}$. This shell is interacting with the dense neutral molecular gas in the W43-Main dense cloud. Combining the high spectral and angular resolution mm-RRL and cm-RRL cubes, we derive the two-dimensional relative distributions of dynamical and pressure broadening of the ionized gas emission and find that the RRL line shapes are dominated by pressure broadening (4-55 km s$^{-1}$) near the YMC and by dynamical broadening (8-36 km s$^{-1}$) near the shell's edge. Ionized gas clumps hosting ultra-compact H II regions found at the edge of the shell suggest that large-scale ionized gas motion triggers the formation of new star generation near the periphery of the shell., Comment: Published at https://iopscience.iop.org/article/10.3847/2041-8213/aa7d48/pdf

Molecular clouds, which harbor the birthplaces of stars, form out of the atomic phase of the interstellar medium (ISM). We aim to characterize the atomic and molecular phases of the ISM and set their physical properties into the context of cloud formation processes. We studied the cold neutral medium (CNM) by means of $\rm HI$ self-absorption (HISA) toward the giant molecular filament GMF20.0-17.9 and compared our results with molecular gas traced by $^{13}\rm CO$ emission. We fitted baselines of HISA features to $\rm HI$ emission spectra using first and second order polynomial functions. The CNM identified by this method spatially correlates with the morphology of the molecular gas toward the western region. However, no spatial correlation between HISA and $^{13}\rm CO$ is evident toward the eastern part of the filament. The distribution of HISA peak velocities and line widths agrees well with $^{13}\rm CO$ within the whole filament. The column density probability density functions (N-PDFs) of HISA (CNM) and $\rm HI$ emission (tracing both the CNM and the warm neutral medium, WNM) have a log-normal shape for all parts of the filament, indicative of turbulent motions as the main driver for these structures. The $\rm H_2$ N-PDFs show a broad log-normal distribution with a power-law tail suggesting the onset of gravitational contraction. The saturation of $\rm HI$ column density is observed at $\sim$25$\rm\,M_{\odot}\,pc^{-2}$. We conjecture that different evolutionary stages are evident within the filament. In the eastern region, we witness the onset of molecular cloud formation out of the atomic gas reservoir while the western part is more evolved, as it reveals pronounced $\rm H_2$ column density peaks and signs of active star formation., Comment: 21 pages, 19 PDF figures

We present a study of the filamentary structure in the emission from the neutral atomic hydrogen (HI) at 21 cm across velocity channels in the 40'' and 1.5-km/s resolution position-position-velocity cube resulting from the combination of the single-dish and interferometric observations in The HI/OH/Recombination (THOR) line survey. Using the Hessian matrix method in combination with tools from circular statistics, we find that the majority of the filamentary structures in the HI emission are aligned with the Galactic plane. Part of this trend can be assigned to long filamentary structures that are coherent across several velocity channels. However, we also find ranges of Galactic longitude and radial velocity where the HI filamentary structures are preferentially oriented perpendicular to the Galactic plane. These are located (i) around the tangent point of the Scutum spiral arm, $l \approx 28^{\circ}$ and $v_{\rm LSR}\approx 100$ km/s, (ii) toward $l \approx 45^{\circ}$ and $v_{\rm LSR}\approx 50$ km/s, (iii) around the Riegel-Crutcher cloud, and (iv) toward the terminal velocities. Comparison with numerical simulations indicates that the prevalence of horizontal filamentary structures is most likely the result of the large-scale dynamics and that vertical structures identified in (i) and (ii) may arise from the combined effect of supernova (SN) feedback and strong magnetic fields. The vertical filamentary structures in (iv) can be related to the presence of clouds from extra-planar HI gas falling back into the Galactic plane after being expelled by SNe. Our results indicate that a systematic characterization of the emission morphology toward the Galactic plane provides an unexplored link between the observations and the dynamical behaviour of the interstellar medium, from the effect of large-scale Galactic dynamics to the Galactic fountains driven by SNe., Comment: 28 pages. 37 figures. Accepted for publication in Astronomy & Astrophysics (09SEP2020)

We deduce a model for glioma invasion making use of DTI data and accounting for the dynamics of brain tissue being actively degraded by tumor cells via excessive acidity production, but also according to the local orientation of tissue fibers. Our approach has a multiscale character: we start with a microscopic description of single cell dynamics including biochemical and/or biophysical effects of the tumor microenvironment, translated on the one hand into cell stress and corresponding forces and on the other hand into receptor binding dynamics; these lead on the mesoscopic level to kinetic equations involving transport terms w.r.t. all kinetic variables and eventually, by appropriate upscaling, to a macroscopic reaction-diffusion equation for glioma density with multiple taxis, coupled to (integro-)differential equations characterizing the evolution of acidity and macro- and mesoscopic tissue. Our approach also allows for a switch between fast and slower moving regimes, %diffusion- and drift-dominated regimes, according to the local tissue anisotropy. We perform numerical simulations to investigate the behavior of solutions w.r.t. various scenarios of tissue dynamics and the dominance of each of the tactic terms, also suggesting how the model can be used to perform a numerical necrosis-based tumor grading or support radiotherapy planning by dose painting. We also provide a discussion about alternative ways of including cell level environmental influences in such multiscale modeling approach, ultimately leading in the macroscopic limit to (multiple) taxis.

We examine the role of the interstellar magnetic field to modulate the orientation of turbulent flows within the Taurus molecular cloud using spatial gradients of thin velocity slices of 12CO and 13CO antenna temperatures. Our analysis accounts for the random errors of the gradients that arise from the thermal noise of the spectra. The orientations of the vectors normal to the antenna temperature gradient vectors are compared to the magnetic field orientations that are calculated from Planck 353~GHz polarization data. These relative orientations are parameterized with the projected Rayleigh statistic and mean resultant vector. For 12CO, 28 percent and 39 percent of the cloud area exhibit strongly parallel or strongly perpendicular relative orientations respectively. For the lower opacity 13CO emission, strongly parallel and strongly perpendicular orientations are found in 7 percent and 43 percent of the cloud area respectively. For both isotopologues, strongly parallel or perpendicular alignments are restricted to localized regions with low levels of turbulence. If the relative orientations serve as an observational proxy to the Alfvenic Mach number then our results imply local variations of the Alfvenic Mach number throughout the cloud., Comment: 20 pages, 15 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows., Comment: 13 pages, 15 figures, Astronomy & Astrophysics in press, a higher resolution version can be found at http://www.mpia.de/homes/beuther/papers.html

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs., Comment: 17 pages, 13 figures, accepted by A&A

We present synthetic dust polarization maps of simulated molecular clouds (MCs) with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\bf{B}$) with respect to the MC sub-structures identified in density ($n$; 3D) and column density ($N$; 2D). The polarization maps are generated with the radiative transfer code POLARIS, including self-consistently calculated efficiencies for radiative torque alignment. The MCs are formed in two sets of 3D MHD simulations: in (i) colliding flows (CF), and (ii) the SILCC-Zoom simulations. In 3D, for the CF simulations with an initial field strength below $\sim$5 $\mu$G, $\bf{B}$ is oriented parallel or randomly with respect to the $n$-structures. For CF runs with stronger initial fields and all SILCC-Zoom simulations, which have an initial field strength of 3 $\mu$G, a flip from parallel to perpendicular orientation occurs at high densities of $n_\text{trans}$ $\simeq$ 10$^2$ - 10$^3$ cm$^{-3}$. We suggest that this flip happens if the MC's mass-to-flux ratio, $\mu$, is close to or below the critical value of 1. This corresponds to a field strength around 3 - 5 $\mu$G. In 2D, we use the Projected Rayleigh Statistics (PRS) to study the orientation of $\bf{B}$. If present, the flip in orientation occurs at $N_\text{trans}$ $\simeq$ 10$^{21 - 21.5}$ cm$^{-2}$, similar to the observed transition value from sub- to supercritical magnetic fields in the ISM. However, projection effects can reduce the power of the PRS method: Depending on the MC or LOS, the projected maps of the SILCC-Zoom simulations do not always show the flip, although expected from the 3D morphology. Such projection effects can explain the variety of recently observed field configurations, in particular within a single MC. Finally, we do not find a correlation between the observed orientation of $\bf{B}$ and the $N$-PDF., Comment: 20 pages, 12 figures, accepted for publication in MNRAS

Future mm-wave and sub-mm space missions will employ large arrays of multiplexed Transition Edge Sensor (TES) bolometers. Such instruments must contend with the high flux of cosmic rays beyond our atmosphere that induce "glitches" in bolometer data, which posed a challenge to data analysis from the Planck bolometers. Future instruments will face the additional challenges of shared substrate wafers and multiplexed readout wiring. In this work we explore the susceptibility of modern TES arrays to the cosmic ray environment of space using two data sets: the 2015 long-duration balloon flight of the SPIDER cosmic microwave background polarimeter, and a laboratory exposure of SPIDER flight hardware to radioactive sources. We find manageable glitch rates and short glitch durations, leading to minimal effect on SPIDER analysis. We constrain energy propagation within the substrate through a study of multi-detector coincidences, and give a preliminary look at pulse shapes in laboratory data., Comment: 9 pages, 6 figures, Proceedings of the 18th International Workshop on Low Temperature Detectors

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes., Comment: 22 pages, 25 figures, accepted for publication by A&A

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$n\alpha$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $\tau$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions., Comment: 22 pages, 13 figures, accepted to be published in Astronomy and Astrophysics

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys., Comment: Accepted for publication in Astrophysical Journal Letters

Using a method based on Faraday rotation measurements, Tahani et al. 2018 find the line-of-sight component of magnetic fields in Orion-A and show that their direction changes from the eastern side of this filamentary structure to its western side. Three possible magnetic field morphologies that can explain this reversal across the Orion-A region are toroidal, helical, and bow-shaped morphologies. In this paper we construct simple models to represent these three morphologies and compare them with the available observational data to find the most probable morphology(ies). To compare the observations with the models, we use probability values and a Monte-Carlo analysis to determine the most likely magnetic field morphology among these three morphologies. We find that the bow morphology has the highest probability values and that our Monte-Carlo analysis suggests that the bow morphology is more likely. We suggest that the bow morphology is the most likely and the most natural of the three morphologies that could explain a magnetic field reversal across the Orion-A filamentary structure (i.e., bow, helical and toroidal morphologies)., Comment: 11 pages, 15 Figures

The Orion-Eridanus superbubble has been blown by supernovae and supersonic winds of the massive stars in the Orion OB associations. The formation history and current structure of the superbubble are still poorly understood. It possibly consists of a combination of nested shells along the line of sight. We have investigated the composite structure of the Eridanus side of the superbubble in the light of a new decomposition of the atomic and molecular gas. We used HI and CO emission lines to separate coherent gas shells in space and velocity, and we studied their relation to the warm ionised gas probed in H$\alpha$ emission, to the hot plasma emitting X-rays, and to the magnetic fields traced by dust polarised emission. We also constrained the relative distances to the clouds using dust reddening maps and X-ray absorption. We used the dust polarisation data to estimate the plane-of-sky components of the magnetic field in several clouds and along the outer rim of the superbubble. Our gas decomposition has revealed several shells inside the superbubble that span distances from about 150 pc to 250 pc. One of these shells forms a nearly complete ring filled with hot plasma. Other shells likely correspond to the layers of swept-up gas that is compressed behind the expanding outer shock wave. We used the gas and magnetic-field data downstream of the shock to derive a shock expansion velocity close to 20 km/s. Taking the X-ray absorption by the gas into account, we find that the pressure of the hot plasma inside the superbubble exceeds that in the Local Bubble. It comprises a mix of hotter (3-9 MK) and cooler (0.3-1.2 MK) gas along the lines of sight. The magnetic field along the western and southern rims and in the approaching wall of the superbubble appears to be shaped and compressed by the ongoing expansion. We find plane-of-sky magnetic strengths ranging from 3 to 15 $\mu$G along the rim., Comment: accepted for publication in A&A