Showing 1–50 of 128 results for author: Soler, D

Probing the Physics of Star-Formation (ProPStar) III. No evidence for dissipation of turbulence down to 20 mpc (4 000 au) scale

Authors: Jaime E. Pineda, Juan D. Soler, Stella Offner, Eric W. Koch, Dominique M. Segura-Cox, Roberto Neri, Michael Kuffmeier, Alexei V. Ivlev, Maria Teresa Valdivia-Mena, Olli Sipilä, Maria Jose Maureira, Paola Caselli, Nichol Cunningham, Anika Schmiedeke, Caroline Gieser, Michael Chen, Silvia Spezzano

Abstract: Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is $k^{-5/3}$, while for supersonic turbulence it is $k^{-2}$. Aims. We aim to determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected ma… ▽ More Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is $k^{-5/3}$, while for supersonic turbulence it is $k^{-2}$. Aims. We aim to determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyze observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from $\sim$10 pc down to 20 mpc. The largest scales are covered with the low density gas tracer $^{13}$CO (1-0) obtained with single dish, the intermediate scales are covered with single-dish observations of the C$^{18}$O (3-2) line, while the smallest scales are covered in H$^{13}$CO$^+$ (1-0) and HNC (1-0) with a combination of NOEMA interferometer and IRAM 30m single dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9$\pm$0.1 and no evidence of dissipation. Moreover, there is evidence for the power-spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite from theoretical expectations. Conclusions. We show new possibilities of studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations. △ Less

Submitted 31 August, 2024; originally announced September 2024.

Comments: 6 pages (+2 appendix), 5 figures. Accepted to A&A Letters. GitHub repo with code available at https://github.com/jpinedaf/NGC1333_NOEMA_turbulence

In-Flight Performance of Spider's 280 GHz Receivers

Authors: Elle C. Shaw, P. A. R. Ade, S. Akers, M. Amiri, J. Austermann, J. Beall, D. T. Becker, S. J. Benton, A. S. Bergman, J. J. Bock, J. R. Bond, S. A. Bryan, H. C. Chiang, C. R. Contaldi, R. S. Domagalski, O. Doré, S. M. Duff, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, L. M. Fissel, A. A. Fraisse, K. Freese, M. Galloway , et al. (62 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, JATIS

Dynamical Accretion Flows -- ALMAGAL: Flows along filamentary structures in high-mass star-forming clusters

Authors: M. R. A. Wells, H. Beuther, S. Molinari, P. Schilke, C. Battersby, P. Ho, Á. Sánchez-Monge, B. Jones, M. B. Scheuck, J. Syed, C. Gieser, R. Kuiper, D. Elia, A. Coletta, A. Traficante, J. Wallace, A. J. Rigby, R. S. Klessen, Q. Zhang, S. Walch, M. T. Beltrán, Y. Tang, G. A. Fuller, D. C. Lis, T. Möller , et al. (25 additional authors not shown)

Abstract: 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 thes… ▽ More 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. △ Less

Submitted 16 August, 2024; v1 submitted 15 August, 2024; originally announced August 2024.

Comments: 11 pages, 11 figures, accepted for publication in A&A

Analysis of Polarized Dust Emission from the First Flight of the SPIDER Balloon-Borne Telescope

Authors: SPIDER Collaboration, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Doré, A. J. Duivenvoorden, H. K. Eriksen, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, A. E. Gambrel, N. N. Gandilo, K. Ganga, S. Gourapura, R. Gualtieri, J. E. Gudmundsson , et al. (45 additional authors not shown)

Abstract: 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 demon… ▽ More 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 $β_\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 $β_\mathrm{d}=1.52\pm0.06$ and $β_\mathrm{d}=1.09\pm0.09$, which are inconsistent at the $3.9\,σ$ 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\,σ$, 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. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 21 pages, 15 figures

The magnetic field of the Radcliffe Wave: starlight polarization at nearest approach to the Sun

Authors: G. V. Panopoulou, C. Zucker, D. Clemens, V. Pelgrims, J. D. Soler, S. E. Clark, J. Alves, A. Goodman, J. Becker Tjus

Abstract: 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… ▽ More 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. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Submitted to A&A

Attribute-Based Authentication in Secure Group Messaging for Distributed Environments

Authors: David Soler, Carlos Dafonte, Manuel Fernández-Veiga, Ana Fernández Vilas, Francisco J. Nóvoa

Abstract: Messaging Layer security (MLS) and its underlying Continuous Group Key Agreement (CGKA) protocol allows a group of users to share a cryptographic secret in a dynamic manner, such that the secret is modified in member insertions and deletions. Although this flexibility makes MLS ideal for implementations in distributed environments, a number of issues need to be overcome. Particularly, the use of d… ▽ More Messaging Layer security (MLS) and its underlying Continuous Group Key Agreement (CGKA) protocol allows a group of users to share a cryptographic secret in a dynamic manner, such that the secret is modified in member insertions and deletions. Although this flexibility makes MLS ideal for implementations in distributed environments, a number of issues need to be overcome. Particularly, the use of digital certificates for authentication in a group goes against the group members' privacy. In this work we provide an alternative method of authentication in which the solicitors, instead of revealing their identity, only need to prove possession of certain attributes, dynamically defined by the group, to become a member. Instead of digital certificates, we employ Attribute-Based Credentials accompanied with Selective Disclosure in order to reveal the minimum required amount of information and to prevent attackers from linking the activity of a user through multiple groups. We formally define a CGKA variant named Attribute-Authenticated Continuous Group Key Agreement (AA-CGKA) and provide security proofs for its properties of Requirement Integrity, Unforgeability and Unlinkability. We also provide guidelines for an integration of our construction in MLS. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: Submitted to IEEE Transactions on Dependable and Secure Computing This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

Faraday tomography of LoTSS-DR2 data: II. Multi-tracer analysis in the high-latitude outer Galaxy

Authors: Ana Erceg, Vibor Jelić, Marijke Haverkorn, Andrea Bracco, Lana Ceraj, Luka Turić, Juan D. Soler

Abstract: We conducted a follow-up study on the analysis of the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) mosaic in the high-latitude outer Galaxy presented in the first paper of this series. Here, we focus on the search for alignment between the magnetic field traced by dust, HI filaments, starlight optical linear polarisation, and linear depolarised structures (depolarisation canals) observ… ▽ More We conducted a follow-up study on the analysis of the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) mosaic in the high-latitude outer Galaxy presented in the first paper of this series. Here, we focus on the search for alignment between the magnetic field traced by dust, HI filaments, starlight optical linear polarisation, and linear depolarised structures (depolarisation canals) observed in low-frequency synchrotron polarisation. This alignment was previously found in several smaller fields observed with LOFAR, offering valuable insights into the nature of the interstellar medium and the 3D spatial distribution of the diffuse ionised medium. We aim to determine whether the alignment of the interstellar medium (ISM) phases observed through multiple tracers is a common occurrence or an exception. Additionally, in areas where depolarisation canals align with the magnetic field, we use starlight polarisation to constrain the distance to the structures associated with the observed canals. We employed the Rolling Hough Transform (RHT) and projected Rayleigh statistics (PRS) to identify and quantify the alignment between the different tracers. On the scale of the whole mosaic, we did not find any evidence of a universal alignment among the three tracers. However, in one particular area, the western region (Dec between $29^\circ$ and $70^\circ$ and RA between $\mathrm{7^h44^m}$ and $\mathrm{9^h20^m}$), we do find a significant alignment between the magnetic field, depolarisation canals, and HI filaments. Based on this alignment, we used the starlight polarisation of stars with known parallax distances to estimate that the minimum distance to the structures observed by LOFAR in this region lies within the range of 200 to 240 pc. We associate these structures with the edge of the Local Bubble. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 11 pages, 10 figures; accepted for publication in A&A; data products are available at https://data.fulir.irb.hr/islandora/object/irb:367

Journal ref: A&A, 687, A23 (2024)

Testing kinematic distances under a realistic Galactic potential

Authors: Glen H. Hunter, Mattia C. Sormani, Jan P. Beckmann, Eugene Vasiliev, Simon C. O. Glover, Ralf S. Klessen, Juan D. Soler, Noé Brucy, Philipp Girichidis, Junia Göller, Loke Ohlin, Robin Tress, Sergio Molinari, Ortwin Gerhard, Milena Benedettini, Rowan Smith, Patrick Hennebelle, Leonardo Testi

Abstract: Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify t… ▽ More Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify the systematic errors that arise from the kinematic distance method in the presence of a Galactic potential that is non-axisymmetric. We investigate how these errors differ in certain regions of the Galaxy and how they relate to the underlying dynamics. We perform 2D isothermal hydrodynamical simulation of the gas disk with the moving-mesh code Arepo, adding the capability of using an external potential provided by the Agama library for galactic dynamics. We introduce a new analytic potential of the Milky Way, taking elements from existing models and adjusting parameters to match recent observational constraints. We find significant errors in the kinematic distance estimate for gas close to the Sun, along sight lines towards the Galactic centre and anti-centre, and significant deviations associated with the Galactic bar. Kinematic distance errors are low within the spiral arms as gas resides close to local potential minima and the resulting line-of-sight velocity is close to what is expected for an axisymmetric potential. Interarm regions exhibit large deviations at any given Galactic radius. This is caused by the gas being sped up or slowed down as it travels into or out of the spiral arm. We are able to define 'zones of avoidance' in the lv-diagram, where the kinematic distance method is particularly unreliable and should only be used with caution. We report a power law relation between the kinematic distance error and the deviation of the project line-of-sight velocity from circular motion. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 22 figures, 6 tables, submitted to A&A

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: 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 f… ▽ More 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 $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_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}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

Cloud properties across spatial scales in simulations of the interstellar medium

Authors: Tine Colman, Noé Brucy, Philipp Girichidis, Simon C. O Glover, Milena Benedettini, Juan D. Soler, Robin G. Tress, Alessio Traficante, Patrick Hennebelle, Ralf S. Klessen, Sergio Molinari, Marc-Antoine Miville-Deschênes

Abstract: Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud… ▽ More Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion and virial state. We apply the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass slope of $\mathrm{d} N/\mathrm{d}\ln M\propto-1$ in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $σ\propto R^{0.75}$ for large clouds ($R\gtrsim3\,\mathrm{pc}$). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $α_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 22 pages, 16 figures, proposed for acceptance in A&A

Journal ref: A&A 686, A155 (2024)

Reinforcement learning to maximise wind turbine energy generation

Authors: Daniel Soler, Oscar Mariño, David Huergo, Martín de Frutos, Esteban Ferrer

Abstract: We propose a reinforcement learning strategy to control wind turbine energy generation by actively changing the rotor speed, the rotor yaw angle and the blade pitch angle. A double deep Q-learning with a prioritized experience replay agent is coupled with a blade element momentum model and is trained to allow control for changing winds. The agent is trained to decide the best control (speed, yaw,… ▽ More We propose a reinforcement learning strategy to control wind turbine energy generation by actively changing the rotor speed, the rotor yaw angle and the blade pitch angle. A double deep Q-learning with a prioritized experience replay agent is coupled with a blade element momentum model and is trained to allow control for changing winds. The agent is trained to decide the best control (speed, yaw, pitch) for simple steady winds and is subsequently challenged with real dynamic turbulent winds, showing good performance. The double deep Q- learning is compared with a classic value iteration reinforcement learning control and both strategies outperform a classic PID control in all environments. Furthermore, the reinforcement learning approach is well suited to changing environments including turbulent/gusty winds, showing great adaptability. Finally, we compare all control strategies with real winds and compute the annual energy production. In this case, the double deep Q-learning algorithm also outperforms classic methodologies. △ Less

Submitted 17 February, 2024; originally announced February 2024.

QKDNetSim+: Improvement of the Quantum Network Simulator for NS-3

Authors: David Soler, Iván Cillero, Carlos Dafonte, Manuel Fernández-Veiga, Ana Fernández-Vilas, Francisco J. Nóvoa

Abstract: The first Quantum Key Distribution (QKD) networks are currently being deployed, but the implementation cost is still prohibitive for most researchers. As such, there is a need for realistic QKD network simulators. The \textit{QKDNetSim} module for the network simulator NS-3 focuses on the representation of packets and the management of key material in a QKD network at the application layer. Althou… ▽ More The first Quantum Key Distribution (QKD) networks are currently being deployed, but the implementation cost is still prohibitive for most researchers. As such, there is a need for realistic QKD network simulators. The \textit{QKDNetSim} module for the network simulator NS-3 focuses on the representation of packets and the management of key material in a QKD network at the application layer. Although QKDNetSim's representation of a QKD network is insightful, some its components lack the depth that would allow the simulator to faithfully represent the behaviour of a real quantum network. In this work, we analyse QKDNetSim's architecture to identify its limitations, and we present an enhanced version of QKDNetSim in which some of its components have been modified to provide a more realistic simulation environment. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 15 pages, 5 figures, preprint submitted to SoftwareX

A Privacy-preserving key transmission protocol to distribute QRNG keys using zk-SNARKs

Authors: David Soler, Carlos Dafonte, Manuel Fernández-Veiga, Ana Fernández Vilas, Francisco J. Nóvoa

Abstract: High-entropy random numbers are an essential part of cryptography, and Quantum Random Number Generators (QRNG) are an emergent technology that can provide high-quality keys for cryptographic algorithms but unfortunately are currently difficult to access. Existing Entropy-as-a-Service solutions require users to trust the central authority distributing the key material, which is not desirable in a h… ▽ More High-entropy random numbers are an essential part of cryptography, and Quantum Random Number Generators (QRNG) are an emergent technology that can provide high-quality keys for cryptographic algorithms but unfortunately are currently difficult to access. Existing Entropy-as-a-Service solutions require users to trust the central authority distributing the key material, which is not desirable in a high-privacy environment. In this paper, we present a novel key transmission protocol that allows users to obtain cryptographic material generated by a QRNG in such a way that the server is unable to identify which user is receiving each key. This is achieved with the inclusion of Zero Knowledge Succinct Non-interactive Arguments of Knowledge (zk-SNARK), a cryptographic primitive that allow users to prove knowledge of some value without needing to reveal it. The security analysis of the protocol proves that it satisfies the properties of Anonymity, Unforgeability and Confidentiality, as defined in this document. We also provide an implementation of the protocol demonstrating its functionality and performance, using NFC as the transmission channel for the QRNG key. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 36 pages, 6 figures. Submitted to Computer Networks

The BLAST Observatory: A Sensitivity Study for Far-IR Balloon-borne Polarimeters

Authors: The BLAST Observatory Collaboration, Gabriele Coppi, Simon Dicker, James E. Aguirre, Jason E. Austermann, James A. Beall, Susan E. Clark, Erin G. Cox, Mark J. Devlin, Laura M. Fissel, Nicholas Galitzki, Brandon S. Hensley, Johannes Hubmayr, Sergio Molinari, Federico Nati, Giles Novak, Eugenio Schisano, Juan D. Soler, Carole E. Tucker, Joel N. Ullom, Anna Vaskuri, Michael R. Vissers, Jordan D. Wheeler, Mario Zannoni

Abstract: Sensitive wide-field observations of polarized thermal emission from interstellar dust grains will allow astronomers to address key outstanding questions about the life cycle of matter and energy driving the formation of stars and the evolution of galaxies. Stratospheric balloon-borne telescopes can map this polarized emission at far-infrared wavelengths near the peak of the dust thermal spectrum… ▽ More Sensitive wide-field observations of polarized thermal emission from interstellar dust grains will allow astronomers to address key outstanding questions about the life cycle of matter and energy driving the formation of stars and the evolution of galaxies. Stratospheric balloon-borne telescopes can map this polarized emission at far-infrared wavelengths near the peak of the dust thermal spectrum - wavelengths that are inaccessible from the ground. In this paper we address the sensitivity achievable by a Super Pressure Balloon (SPB) polarimetry mission, using as an example the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) Observatory. By launching from Wanaka, New Zealand, BLAST Observatory can obtain a 30-day flight with excellent sky coverage - overcoming limitations of past experiments that suffered from short flight duration and/or launch sites with poor coverage of nearby star-forming regions. This proposed polarimetry mission will map large regions of the sky at sub-arcminute resolution, with simultaneous observations at 175, 250, and 350 $μm$, using a total of 8274 microwave kinetic inductance detectors. Here, we describe the scientific motivation for the BLAST Observatory, the proposed implementation, and the forecasting methods used to predict its sensitivity. We also compare our forecasted experiment sensitivity with other facilities. △ Less

Submitted 23 May, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: Published in PASP

Journal ref: 2024 PASP 136 035003

Density distributions, magnetic field structures and fragmentation in high-mass star formation

Authors: H. Beuther, C. Gieser, J. D. Soler, Q. Zhang, R. Rao, D. Semenov, Th. Henning, R. Pudritz, T. Peters, P. Klaassen, M. T. Beltran, A. Palau, T. Moeller, K. G. Johnston, H. Zinnecker, J. Urquhart, R. Kuiper, A. Ahmadi, A. Sanchez-Monge, S. Feng, S. Leurini, S. E. Ragan

Abstract: 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… ▽ More 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. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Evidence for Spatial Separation of Galactic Dust Populations

Authors: Corwin Shiu, Steven J. Benton, Jeffrey P. Filippini, Aurélien A. Fraisse, William C. Jones, Johanna M. Nagy, Ivan L. Padilla, Juan D. Soler

Abstract: We present an implementation of a Bayesian mixture model using Hamiltonian Monte Carlo (HMC) techniques to search for spatial separation of Galactic dust populations. Utilizing intensity measurements from Planck High Frequency Instrument (HFI), we apply this model to high-latitude Galactic dust emission. Our analysis reveals a strong preference for a spatially-varying two-population dust model ove… ▽ More We present an implementation of a Bayesian mixture model using Hamiltonian Monte Carlo (HMC) techniques to search for spatial separation of Galactic dust populations. Utilizing intensity measurements from Planck High Frequency Instrument (HFI), we apply this model to high-latitude Galactic dust emission. Our analysis reveals a strong preference for a spatially-varying two-population dust model over a one-population dust model, when the latter must capture the total variance in the sky. Each dust population is well characterized by a single-component spectral energy distribution (SED) and accommodates small variations. These populations could signify two distinct components, or may originate from a one-component model with different temperatures resulting in different SED scalings. While no spatial information is built into the likelihood, our investigation unveils large-scale spatially coherent structures with high significance, pointing to a physical origin for the observed spatial variation. These results are robust to our choice of likelihood and of input data. Furthermore, this spatially varying two-population model is the most favored from Bayesian evidence calculations. Incorporating IRAS 100 $μ$m to constrain the Wein-side of the blackbody function, we find the dust populations differ at the 2.5$σ$ level in the spectral index ($β_d$) vs. temperature ($T_d$) plane. The presence of multiple dust populations has implications for component separation techniques frequently employed in the recovery of the cosmic microwave background. △ Less

Submitted 6 May, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJ

Turbulent Structure In Supernova Remnants G46.8-0.3 And G39.2-0.3 From THOR Polarimetry

Authors: Russell Shanahan, Jeroen Stil, Loren Anderson, Henrik Beuther, Paul Goldsmith, Ralf S. Klessen, Michael Rugel, Juan D. Soler

Abstract: We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Farada… ▽ More We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Faraday dispersion (sigma_phi). From THOR, we find median sigma_phi of 15.9 +/- 3.2 rad m^-2 for G46.8-0.3 and 17.6 +/- 1.6 rad m^-2 for G39.2-0.3. When comparing to polarization at 6cm, we find evidence for sigma_phi >~ 30 rad m^-2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium (ISM) using the SF of extra-Galactic radio sources (EGRS) and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Has been accepted by ApJ for publication. Figures 3 and 4 did do not render well when using an internet browser to view them, but when the pdf file is downloaded these figures look as they should

Cold atomic gas identified by HI self-absorption. Cold atomic clouds toward giant molecular filaments

Authors: J. Syed, H. Beuther, P. F. Goldsmith, Th. Henning, M. Heyer, R. S. Klessen, J. M. Stil, J. D. Soler, L. D. Anderson, J. S. Urquhart, M. R. Rugel, K. G. Johnston, A. Brunthaler

Abstract: 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 dat… ▽ More 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. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 41 pages, 28 figures, accepted for publication in A&A

Journal ref: A&A 679, A130 (2023)

The Orion-Taurus ridge: a synchrotron radio loop at the edge of the Orion-Eridanus superbubble

Authors: Andrea Bracco, Marco Padovani, Juan D. Soler

Abstract: Large-scale synchrotron loops are recognized as the main source of diffuse radio-continuum emission in the Galaxy at intermediate and high Galactic latitudes. Their origin, however, remains rather unexplained. Using a combination of multi-frequency data in the radio band of total and polarized intensities, for the first time in this letter, we associate one arc -- hereafter, the Orion-Taurus ridge… ▽ More Large-scale synchrotron loops are recognized as the main source of diffuse radio-continuum emission in the Galaxy at intermediate and high Galactic latitudes. Their origin, however, remains rather unexplained. Using a combination of multi-frequency data in the radio band of total and polarized intensities, for the first time in this letter, we associate one arc -- hereafter, the Orion-Taurus ridge -- with the wall of the most prominent stellar-feedback blown shell in the Solar neighborhood, namely the Orion-Eridanus superbubble. We traced the Orion-Taurus ridge using 3D maps of interstellar dust extinction and column-density maps of molecular gas, $N_{\rm H_2}$. We found the Orion-Taurus ridge at a distance of 400\,pc, with a plane-of-the-sky extent of $180$\,pc. Its median $N_{\rm H_2}$ value is $(1.4^{+2.6}_{-0.6})\times 10^{21}$ cm$^{-2}$. Thanks to the broadband observations below 100 MHz of the Long Wavelength Array, we also computed the low-frequency spectral-index map of synchrotron emissivity, $β$, in the Orion-Taurus ridge. We found a flat distribution of $β$ with a median value of $-2.24^{+0.03}_{-0.02}$ that we interpreted in terms of depletion of low-energy ($<$ GeV) cosmic-ray electrons in recent supernova remnants ($10^5$ - $10^6$ yrs). Our results are consistent with plane-of-the-sky magnetic-field strengths in the Orion-Taurus ridge larger than a few tens of $μ$G ($> 30 - 40 \,μ$G). We report the first detection of diffuse synchrotron emission from cold-neutral, partly molecular, gas in the surroundings of the Orion-Eridanus superbubble. This observation opens a new perspective to study the multiphase and magnetized interstellar medium with the advent of future high-sensitivity radio facilities, such as the C-Band All-Sky Survey and the Square Kilometre Array. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, L11 (2023)

A new tool to derive simultaneously exponent and extremes of power-law distributions

Authors: S. Pezzuto, A. Coletta, R. S. Klessen, E. Schisano, M. Benedettini, D. Elia, S. Molinari, J. D. Soler, A. Traficante

Abstract: Many experimental quantities show a power-law distribution $p(x)\propto x^{-α}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $α$ 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.… ▽ More Many experimental quantities show a power-law distribution $p(x)\propto x^{-α}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $α$ 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 $α$ 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 $α$. 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: $α=2.576$ between $1.06\,M_{\sun}$ and $3.35\,M_{\sun}$, $α=3.39$ between $3.48\,M_{\sun}$ and $33.4\,M_{\sun}$; ii) the slope of the $γ$-ray spectrum of the blazar J0011.4+0057, extracted from the Fermi-LAT archive. For the latter case, we derive $α=2.89$ between 1,484~MeV and 28.7~GeV; then we derive the time-resolved slopes using subsets 200 photons each. △ Less

Submitted 19 September, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 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"

Journal ref: MNRAS, 525, 4744 (2023)

A deep-learning approach to the 3D reconstruction of dust density and temperature in star-forming regions

Authors: Victor F. Ksoll, Stefan Reissl, Ralf S. Klessen, Ian W. Stephens, Rowan J. Smith, Juan D. Soler, Alessio Traficante, Leonardo Testi, Patrick Hennebelle, Sergio Molinari

Abstract: Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic du… ▽ More Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 $μ$m. We simplify the task by reconstructing the cloud structure along individual lines of sight and train a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and is able to predict full posterior distributions for the target dust properties. We test different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluate the predictive performance of these models on synthetic test data. Results: We report an excellent reconstruction performance for the 23-wavelengths cINN model, achieving median absolute relative errors of about 1.8% in $\log(n/m^{-3})$ and 1% in $\log(T_{dust}/K)$, respectively. We identify trends towards overestimation at the low end of the density range and towards underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. Limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 3.3% and 2.5% in $\log(n/m^{-3})$ and $\log(T_{dust}/K)$. Conclusions: This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and even compatible with a more realistically constrained wavelength coverage. △ Less

Submitted 9 February, 2024; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 38 pages, 25 figures, accepted for publication in Astronomy & Astrophysics on 26.01.2024

A comparison of the Milky Way's recent star formation revealed by dust thermal emission and high-mass stars

Authors: J. D. Soler, E. Zari, D. Elia, S. Molinari, C. Mininni, E. Schisano, A. Traficante, R. S. Klessen, S. C. O. Glover, P. Hennebelle, T. Colman, N. Frankel, T. Wenger

Abstract: We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\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… ▽ More We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\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 $Σ_{\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. △ Less

Submitted 18 September, 2023; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication at Astronomy & Astrophysics (11SEP2023). 8 pages, 7 figures

Journal ref: A&A 678, A95 (2023)

A panoptic view of the Taurus molecular cloud I. The cloud dynamics revealed by gas emission and 3D dust

Authors: J. D. Soler, C. Zucker, J. E. G. Peek, M. Heyer, P. F. Goldsmith, S. C. O. Glover, S. Molinari, R. S. Klessen, P. Hennebelle, L. Testi, T. Colman, M. Benedettini, D. Elia, C. Mininni, S. Pezzuto, E. Schisano, A. Traficante

Abstract: 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) metho… ▽ More 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. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 21 figures. Accepted for publication in Astronomy & Astrophysics (22MAY2023)

Journal ref: A&A 675, A206 (2023)

Highly-Entangled Polyradical Nanographene with Coexisting Strong Correlation and Topological Frustration

Authors: Shaotang Song, Andrés Pinar Solé, Adam Matěj, Guangwu Li, Oleksandr Stetsovych, Diego Soler, Huimin Yang, Mykola Telychko, Jing Li, Manish Kumar, Jiri Brabec, Libor Veis, Jishan Wu, Pavel Jelinek, Jiong Lu

Abstract: Open-shell benzenoid polycyclic aromatic hydrocarbons, known as magnetic nanographenes, exhibit unconventional p-magnetism arising from topological frustration or strong electronic-electron (e-e) interaction. Imprinting multiple strongly entangled spins into polyradical nanographenes creates a major paradigm shift in realizing non-trivial collective quantum behaviors and exotic quantum phases in o… ▽ More Open-shell benzenoid polycyclic aromatic hydrocarbons, known as magnetic nanographenes, exhibit unconventional p-magnetism arising from topological frustration or strong electronic-electron (e-e) interaction. Imprinting multiple strongly entangled spins into polyradical nanographenes creates a major paradigm shift in realizing non-trivial collective quantum behaviors and exotic quantum phases in organic quantum materials. However, conventional design approaches are limited by a single magnetic origin, which can restrict the number of correlated spins or the type of magnetic ordering in open-shell nanographenes. Here, we present a novel design strategy combing topological frustration and e-e interactions to fabricate the largest fully-fused open-shell nanographene reported to date, a 'butterfly'-shaped tetraradical on Au(111). We employed bond-resolved scanning tunneling microscopy and spin excitation spectroscopy to unambiguously resolve the molecular backbone and reveal the strongly correlated open-shell character, respectively. This nanographene contains four unpaired electrons with both ferromagnetic and anti-ferromagnetic interactions, harboring a many-body singlet ground state and strong multi-spin entanglement, which can be well described by many-body calculations. Furthermore, we demonstrate that the nickelocene magnetic probe can sense highly-correlated spin states in nanographene. The ability to imprint and characterize many-body strongly correlated spins in polyradical nanographenes not only presents exciting opportunities for realizing non-trivial quantum magnetism and phases in organic materials but also paves the way toward high-density ultrafast spintronic devices and quantum information technologies. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 17 pages, 4 figures

HI filaments as potential compass needles? Comparing the magnetic field structure of the Small Magellanic Cloud to the orientation of GASKAP-HI filaments

Authors: Y. K. Ma, N. M. McClure-Griffiths, S. E. Clark, S. J. Gibson, J. Th. van Loon, J. D. Soler, M. E. Putman, J. M. Dickey, M. -Y. Lee, K. E. Jameson, L. Uscanga, J. Dempsey, H. Dénes, C. Lynn, N. M. Pingel

Abstract: 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 alig… ▽ More 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. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 24 pages, MNRAS accepted

On the distribution of the Cold Neutral Medium in galaxy discs

Authors: Rowan J. Smith, Robin Tress, Juan D. Soler, Ralf S. Klessen, Simon C. O. Glover, Patrick Hennebelle, Sergio Molinari, Mordecai-Mark Mac Low, David Whitworth

Abstract: The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and… ▽ More The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total HI due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii.However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM doesn't extend out radially as far as HI, and the vertical scale height is smaller in the outer galaxy compared to HI with no flaring. The CNM column density scales with total midplane pressure and disappears from the gas phase below values of PT/kB =1000 K/cm3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt-Schmidt law. In the outer galaxy we produce realistic vertical velocity dispersions in the HI purely from galactic dynamics but our models do not predict CNM at the extremely large radii observed in HI absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii. △ Less

Submitted 19 July, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 13 pages, 20 figures, accepted by MNRAS

Synthetic dust polarization emission maps at 353 GHz for an observer placed inside a Local Bubble-like cavity

Authors: E. Maconi, J. D. Soler, S. Reissl, P. Girichidis, R. S. Klessen, P. Hennebelle, S. Molinari, L. Testi, R. J. Smith, M. C. Sormani, J. W. Teh, A. Traficante, .

Abstract: 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 f… ▽ More 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. △ Less

Submitted 19 June, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

A 3D View of Orion: I. Barnard's Loop

Authors: Michael M. Foley, Alyssa Goodman, Catherine Zucker, John C. Forbes, Ralf Konietzka, Cameren Swiggum, João Alves, John Bally, Juan D. Soler, Josefa E. Großschedl, Shmuel Bialy, Michael Y. Grudić, Reimar Leike, Torsten Ensslin

Abstract: Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Br… ▽ More Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Briceño-1 (OBP-B1) cluster in the middle of a large dust cavity. The large-scale gas in the region also appears to be expanding from a central point, originally proposed to be Orion X. OBP-B1 appears to serve as another possible center, and we evaluate whether Orion X or OBP-B1 is more likely to be the cause of the expansion. We find that neither cluster served as the single expansion center, but rather a combination of feedback from both likely propelled the expansion. Recent 3D dust maps are used to characterize the 3D topology of the entire region, which shows Barnard's Loop's correspondence with a large dust cavity around the OPB-B1 cluster. The molecular clouds Orion A, Orion B, and Orion $λ$ reside on the shell of this cavity. Simple estimates of gravitational effects from both stars and gas indicate that the expansion of this asymmetric cavity likely induced anisotropy in the kinematics of OBP-B1. We conclude that feedback from OBP-B1 has affected the structure of the Orion A, Orion B, and Orion $λ$ molecular clouds and may have played a major role in the formation of Barnard's Loop. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 25 pages, 10 figures. Submitted to ApJ

The Star Formation Rate of the Milky Way as seen by Herschel

Authors: D. Elia, S. Molinari, E. Schisano, J. D. Soler, M. Merello, D. Russeil, M. Veneziani, A. Zavagno, A. Noriega-Crespo, L. Olmi, M. Benedettini, P. Hennebelle, R. S. Klessen, S. Leurini, R. Paladini, S. Pezzuto, A. Traficante, D. J. Eden, P. G. Martin, M. Sormani, A. Coletta, T. Colman, R. Plume, Y. Maruccia, C. Mininni , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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(Σ_\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. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: To be published in ApJ

How magnetic field and stellar radiative feedback influences the collapse and the stellar mass spectrum of a massive star forming clump

Authors: Patrick Hennebelle, Ugo Lebreuilly, Tine Colman, Davide Elia, Gary Fuller, Silvia Leurini, Thomas Nony, Eugenio Schisano, Juan D. Soler, Alessio Traficante, Ralf S. Klessen, Sergio Molinari, Leonardo Testi

Abstract: In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiati… ▽ More In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiative feedback and magnetic field with spatial resolution down to 1 AU. Both ideal and non-ideal MHD runs are performed and various radiative feedback efficiencies are considered. We also develop analytical models that we confront to the numerical results. The sum of the luminosities produced by the stars in the calculations is computed and it compares well with the bolometric luminosities reported in observations of massive star forming clumps. The temperatures, velocities and densities are also found to be in good agreement with recent observations. The stellar mass spectrum inferred for the simulations is, generally speaking, not strictly universal and in particular varies with magnetic intensity. It is also influenced by the choice of the radiative feedback efficiency. In all simulations, a sharp drop in the stellar distribution is found at about $M_{min} \simeq$ 0.1 M$_\odot$, which is likely a consequence of the adiabatic behaviour induced by dust opacities at high densities. As a consequence, when the combination of magnetic and thermal support is not too large, the mass distribution presents a peak located at 0.3-0.5 M$_\odot$. When magnetic and thermal support are large, the mass distribution is better described by a plateau, i.e. $d N / d \log M \propto M^{-Γ}$, $Γ\simeq 0$. Abridged △ Less

Submitted 22 October, 2022; originally announced October 2022.

Comments: accepted for publication in A&A

Journal ref: A&A 668, A147 (2022)

Properties of atomic hydrogen gas in the Galactic plane from THOR 21-cm absorption spectra: a comparison with the high latitude gas

Authors: Arghyadeep Basu, Nirupam Roy, Henrik Beuther, Jonas Syed, Jürgen Ott, Juan D. Soler, Jeroen Stil, Michael R. Rugel

Abstract: The neutral hydrogen 21 cm line is an excellent tracer of the atomic interstellar medium in the cold and the warm phases. Combined 21 cm emission and absorption observations are very useful to study the properties of the gas over a wide range of density and temperature. In this work, we have used 21 cm absorption spectra from recent interferometric surveys, along with the corresponding emission sp… ▽ More The neutral hydrogen 21 cm line is an excellent tracer of the atomic interstellar medium in the cold and the warm phases. Combined 21 cm emission and absorption observations are very useful to study the properties of the gas over a wide range of density and temperature. In this work, we have used 21 cm absorption spectra from recent interferometric surveys, along with the corresponding emission spectra from earlier single dish surveys to study the properties of the atomic gas in the Milky Way. In particular, we focus on a comparison of properties between lines of sight through the gas disk in the Galactic plane and high Galactic latitude lines of sight through more diffuse gas. As expected, the analysis shows a lower average temperature for the gas in the Galactic plane compared to that along the high latitude lines of sight. The gas in the plane also has a higher molecular fraction, showing a sharp transition and flattening in the dust - gas correlation. On the other hand, the observed correlation between 21 cm brightness temperature and optical depth indicates some intrinsic difference in spin temperature distribution and a fraction of gas in the Galactic plane having intermediate optical depth (for 0.02 < $τ$ < 0.2) but higher spin temperature, compared to that of the diffuse gas at high latitude with the same optical depth. This may be due to a small fraction of cold gas with slightly higher temperature and lower density present on the Galactic plane. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 6 pages, 5 figures, submitted and accepted for publication in MNRAS

Magnetic fields do not suppress global star formation in low metallicity dwarf galaxies

Authors: David J. Whitworth, Rowan J. Smith, Ralf S. Klessen, Mordecai-Mark Mac Low, Simon C. O. Glover, Robin Tress, Rudiger Pakmor, Juan D. Soler

Abstract: Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a wea… ▽ More Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution AREPO simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10$^{-6} μ$G and different sub-solar metallicities, and one starts with a saturated field of 10$^{-2} μ$G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1\% and 10\% of solar values. We find that the magnetic field has little impact on the global star formation rate, which is in tension with some previously published results. We further find that the initial field strength has little impact on the global star formation rate. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation. △ Less

Submitted 12 January, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 19 pages, 19 figures, accepted for publication in MNRAS

The signature of large scale turbulence driving on the structure of the interstellar medium

Authors: Tine Colman, Jean-François Robitaille, Patrick Hennebelle, Marc-Antoine Miville-Deschênes, Noé Brucy, Ralf S. Klessen, Simon C. O. Glover, Juan D. Soler, Davide Elia, Alessio Traficante, Sergio Molinari, Leonardo Testi

Abstract: The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stell… ▽ More The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stellar feedback. Large scale external turbulent driving of various intensities is added to mimic galactic driving mechanisms. We analyse the resulting column density maps with a technique called Multi-scale non-Gaussian segmentation that separates the coherent structures and the Gaussian background. This effectively discriminates between the various simulations and is a promising method to understand the ISM structure. In particular the power spectrum of the coherent structures flattens above 60 pc when turbulence is driven only by stellar feedback. When large-scale driving is applied, the turn-over shifts to larger scales. A systematic comparison with the Large Magellanic Cloud (LMC) is then performed. Only 1 out of 25 regions has a coherent power spectrum which is consistent with the feedback-only simulation. A detailed study of the turn-over scale leads us to conclude that regular stellar feedback is not enough to explain the observed ISM structure on scales larger than 60 pc. Extreme feedback in the form of supergiant shells likely plays an important role but cannot explain all the regions of the LMC. If we assume ISM structure is generated by turbulence, another large scale driving mechanism is needed to explain the entirety of the observations. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 15 pages, 11 figures, to be published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 3, August 2022

The Galactic dynamics revealed by the filamentary structure in atomic hydrogen emission

Authors: Juan D. Soler, Marc-Antoine Miville-Deschênes, Sergio Molinari, Ralf S. Klessen, Patrick Hennebelle, Leonardo Testi, Naomi M. McClure-Griffiths, Henrik Beuther, Davide Elia, Eugenio Schisano, Alessio Traficante, Philipp Girichidis, Simon C. O. Glover, Rowan J. Smith, Mattia Sormani, Robin Treß

Abstract: We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyo… ▽ More We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyond 10 kpc and up to roughly 18 kpc from the Galactic center display HI filamentary structures predominantly parallel to the Galactic plane. For regions at lower Galactocentric radii, we found that the HI filaments are mostly perpendicular or do not have a preferred orientation with respect to the Galactic plane. We interpret these results as the imprint of supernova feedback in the inner Galaxy and Galactic rotation in the outer Milky Way. We found that the HI filamentary structures follow the Galactic warp and that they highlight some of the variations interpreted as the effect of the gravitational interaction with satellite galaxies. In addition, the mean scale height of the filamentary structures is lower than that sampled by the bulk of the HI emission, thus indicating that the cold and warm atomic hydrogen phases have different scale heights in the outer galaxy. Finally, we found that the fraction of the column density in HI filaments is almost constant up to approximately 18 kpc from the Galactic center. This is possibly a result of the roughly constant ratio between the cold and warm atomic hydrogen phases inferred from the HI absorption studies. Our results indicate that the HI filamentary structures provide insight into the dynamical processes shaping the Galactic disk. Their orientations record how and where the stellar energy input, the Galactic fountain process, the cosmic ray diffusion, and the gas accretion have molded the diffuse interstellar medium in the Galactic plane. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 33 pages. 39 figures. Accepted for publication in Astronomy & Astrophysics (06MAY2022)

Report number: AA/2022/43334

Journal ref: A&A 662, A96 (2022)

From Bubbles and Filaments to Cores and Disks: Gas Gathering and Growth of Structure Leading to the Formation of Stellar Systems

Authors: Jaime E. Pineda, Doris Arzoumanian, Philippe André, Rachel K. Friesen, Annie Zavagno, Seamus D. Clarke, Tsuyoshi Inoue, Che-Yu Chen, Yueh-Ning Lee, Juan D. Soler, Michael Kuffmeier

Abstract: The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented… ▽ More The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50-100 pc scales through parsec-scale molecular filaments down to $<$0.1 pc dense cores and $<$1000 au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation. △ Less

Submitted 23 March, 2023; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 42 pages, 20 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Small changes to match proofs and the updated format

In-flight gain monitoring of SPIDER's transition-edge sensor arrays

Authors: J. P. Filippini, A. E. Gambrel, A. S. Rahlin, E. Y. Young, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Dore, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, A. A. Fraisse, K. Freese, M. Galloway, N. N. Gandilo, K. Ganga, R. Gualtieri , et al. (45 additional authors not shown)

Abstract: 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 respons… ▽ More 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. △ Less

Submitted 16 June, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 7 pages, 3 figures; Proceedings of the 19th International Workshop on Low Temperature Detectors (LTD19); Minor updates to match published version

Journal ref: Journal of Low Temperature Physics (2022)

A Simulation-Based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Authors: J. S. -Y. Leung, J. Hartley, J. M. Nagy, C. B. Netterfield, J. A. Shariff, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Doré, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, A. E. Gambrel , et al. (45 additional authors not shown)

Abstract: 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-… ▽ More 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. △ Less

Submitted 21 April, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 14 pages, 7 figures; updated to match published version

Journal ref: ApJ 928(2):109, 2022

The "Maggie" filament: Physical properties of a giant atomic cloud

Authors: J. Syed, J. D. Soler, H. Beuther, Y. Wang, S. Suri, J. D. Henshaw, M. Riener, S. Bialy, S. Rezaei Kh., J. M. Stil, P. F. Goldsmith, M. R. Rugel, S. C. O. Glover, R. S. Klessen, J. Kerp, J. S. Urquhart, J. Ott, N. Roy, N. Schneider, R. J. Smith, S. N. Longmore, H. Linz

Abstract: 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 M… ▽ More 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. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 19 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 657, A1 (2022)

A wind-blown bubble in the Central Molecular Zone cloud G0.253+0.016

Authors: J. D. Henshaw, M. R. Krumholz, N. O. Butterfield, J. Mackey, A. Ginsburg, T. J. Haworth, F. Nogueras-Lara, A. T. Barnes, S. N. Longmore, J. Bally, J. M. D. Kruijssen, E. A. C. Mills, H. Beuther, D. L. Walker, C. Battersby, A. Bulatek, T. Henning, J. Ott, J. D. Soler

Abstract: 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… ▽ More 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. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 17 pages, 7 figures. Accepted for publication in MNRAS (October 15, 2021)

Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953

Authors: C. Gieser, H. Beuther, D. Semenov, S. Suri, J. D. Soler, H. Linz, J. Syed, Th. Henning, S. Feng, T. Möller, A. Palau, J. M. Winters, M. T. Beltrán, R. Kuiper, L. Moscadelli, P. Klaassen, J. S. Urquhart, T. Peters, S. N. Longmore, Á. Sánchez-Monge, R. Galván-Madrid, R. E. Pudritz, K. G. Johnston

Abstract: 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 t… ▽ More 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. △ Less

Submitted 14 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 657, A3 (2022)

On the accuracy of HI observations in molecular clouds -- More cold HI than thought?

Authors: D. Seifried, H. Beuther, S. Walch, J. Syed, J. D. Soler, P. Girichidis, R. Wünsch

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: 20 pages including appendix, 24 figures, accepted for publication in MNRAS

The filamentary structures in the CO emission toward the Milky Way disk

Authors: J. D. Soler, H. Beuther, J. Syed, Y. Wang, Th. Henning, S. C. O. Glover, R. S. Klessen, M. C. Sormani, M. Heyer, R. J. Smith, J. S. Urquhart, J. Yang, Y. Su, X. Zhou

Abstract: 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°< l < 49.7°$, $|b| \leq 1.25°$, 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 p… ▽ More 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°< l < 49.7°$, $|b| \leq 1.25°$, 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. △ Less

Submitted 25 June, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: Accepted for publication in Astronomy & Astrophysics (7 pages, 7 figures). arXiv admin note: text overlap with arXiv:2007.07285

Report number: AA/2021/41327

Journal ref: A&A 651, L4 (2021)

Fragmentation and kinematics in high-mass star formation: CORE-extension targeting two very young high-mass star-forming regions

Authors: H. Beuther, C. Gieser, S. Suri, H. Linz, P. Klaassen, D. Semenov, J. M. Winters, Th. Henning, J. D. Soler, J. S. Urquhart, J. Syed, S . Feng, T. Moeller, M. T. Beltran, A. Sanchez-Monge, S. N. Longmore, T. Peters, J. Ballesteros-Paredes, P. Schilke, L. Moscadelli, A. Palau, R. Cesaroni, S. Lumsden, R. Pudritz, F. Wyrowski , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 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

Journal ref: A&A 649, A113 (2021)

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

Authors: A. E. Gambrel, A. S. Rahlin, X. Song, C. R. Contaldi, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, N. N. Gandilo, R. Gualtieri, J. E. Gudmundsson, M. Halpern , et al. (42 additional authors not shown)

Abstract: 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 method… ▽ More 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. △ Less

Submitted 24 May, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: 18 pages, 11 figures

A Constraint on Primordial $B$-Modes from the First Flight of the SPIDER Balloon-Borne Telescope

Authors: SPIDER Collaboration, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Doré, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, A. E. Gambrel, N. N. Gandilo, K. Ganga, R. Gualtieri, J. E. Gudmundsson , et al. (46 additional authors not shown)

Abstract: 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 test… ▽ More 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. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: 29 pages, 13 figures

Design and pre-flight performance of SPIDER 280 GHz receivers

Authors: E. C. Shaw, P. A. R. Ade, S. Akers, M. Amiri, J. Austermann, J. Beall, D. T. Becker, S. J. Benton, A. S. Bergman, J. J. Bock, J. R. Bond, S. A. Bryan, H. C. Chiang, C. R. Contaldi, R. S. Domagalski, O. Doré, S. M. Duff, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, L. M. Fissel, A. A. Fraisse, K. Freese, M. Galloway , et al. (57 additional authors not shown)

Abstract: 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 th… ▽ More 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. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

Journal ref: "Design and pre-flight performance of SPIDER 280 GHz receivers," Proc. SPIE 11453, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, 114532F (13 December 2020)

The Balloon-Borne Large Aperture Submillimeter Telescope Observatory

Authors: Ian Lowe, Gabriele Coppi, Peter A. R. Ade, Peter C. Ashton, Jason E. Austermann, James Beall, Susan Clark, Erin G. Cox, Mark J. Devlin, Simon Dicker, Bradley J. Dober, Valentina Fanfani, Laura M. Fissel, Nicholas Galitzki, Jiansong Gao, Brandon Hensley, Johannes Hubmayr, Steven Li, Zhi-Yun Li, Nathan P. Lourie, Peter G. Martin, Philip Mauskopf, Federico Nati, Giles Novak, Giampaolo Pisano , et al. (8 additional authors not shown)

Abstract: The BLAST Observatory is a proposed superpressure balloon-borne polarimeter designed for a future ultra-long duration balloon campaign from Wanaka, New Zealand. To maximize scientific output while staying within the stringent superpressure weight envelope, BLAST will feature new 1.8m off-axis optical system contained within a lightweight monocoque structure gondola. The payload will incorporate a… ▽ More The BLAST Observatory is a proposed superpressure balloon-borne polarimeter designed for a future ultra-long duration balloon campaign from Wanaka, New Zealand. To maximize scientific output while staying within the stringent superpressure weight envelope, BLAST will feature new 1.8m off-axis optical system contained within a lightweight monocoque structure gondola. The payload will incorporate a 300L $^4$He cryogenic receiver which will cool 8,274 microwave kinetic inductance detectors (MKIDs) to 100mK through the use of an adiabatic demagnetization refrigerator (ADR) in combination with a $^3$He sorption refrigerator all backed by a liquid helium pumped pot operating at 2K. The detector readout utilizes a new Xilinx RFSOC-based system which will run the next-generation of the BLAST-TNG KIDPy software. With this instrument we aim to answer outstanding questions about dust dynamics as well as provide community access to the polarized submillimeter sky made possible by high-altitude observing unrestricted by atmospheric transmission. The BLAST Observatory is designed for a minimum 31-day flight of which 70$\%$ will be dedicated to observations for BLAST scientific goals and the remaining 30$\%$ will be open to proposals from the wider astronomical community through a shared-risk proposals program. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: Presented at SPIE Ground-based and Airborne Telescopes VIII, December 13-18, 2020

Characterization, deployment, and in-flight performance of the BLAST-TNG cryogenic receiver

Authors: Ian Lowe, Peter A. R. Ade, Peter C. Ashton, Jason E. Austermann, Gabriele Coppi, Erin G. Cox, Mark J. Devlin, Bradley J. Dober, Valentina Fanfani, Laura M. Fissel, Nicholas Galitzki, Jiansong Gao, Samuel Gordon, Christopher E. Groppi, Gene C. Hilton, Johannes Hubmayr, Jeffrey Klein, Dale Li, Nathan P. Lourie, Hamdi Mani, Philip Mauskopf, Christopher McKenney, Federico Nati, Giles Novak, Giampaolo Pisano , et al. (8 additional authors not shown)

Abstract: The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter polarimeter designed to map interstellar dust and galactic foregrounds at 250, 350, and 500 microns during a 24-day Antarctic flight. The BLAST-TNG detector arrays are comprised of 918, 469, and 272 MKID pixels, respectively. The pixels are formed from two orthogonally oriented, crossed, linear-… ▽ More The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter polarimeter designed to map interstellar dust and galactic foregrounds at 250, 350, and 500 microns during a 24-day Antarctic flight. The BLAST-TNG detector arrays are comprised of 918, 469, and 272 MKID pixels, respectively. The pixels are formed from two orthogonally oriented, crossed, linear-polarization sensitive MKID antennae. The arrays are cooled to sub 300mK temperatures and stabilized via a closed cycle $^3$He sorption fridge in combination with a $^4$He vacuum pot. The detectors are read out through a combination of the second-generation Reconfigurable Open Architecture Computing Hardware (ROACH2) and custom RF electronics designed for BLAST-TNG. The firmware and software designed to readout and characterize these detectors was built from scratch by the BLAST team around these detectors, and has been adapted for use by other MKID instruments such as TolTEC and OLIMPO. We present an overview of these systems as well as in-depth methodology of the ground-based characterization and the measured in-flight performance. △ Less

Submitted 12 January, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: Presented at SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, December 13-18, 2020

Journal ref: Proc. SPIE 11453, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, 1145304 (13 December 2020)

In-flight performance of the BLAST-TNG telescope platform

Authors: Gabriele Coppi, Peter A. R. Ade, Peter C. Ashton, Jason E. Austermann, Erin G. Cox, Mark J. Devlin, Bradley J. Dober, Valentina Fanfani, Laura M. Fissel, Nicholas Galitzki, Jiansong Gao, Samuel Gordon, Christopher E. Groppi, Gene C. Hilton, Johannes Hubmayr, Jeffrey Klein, Dale Li, Nathan P. Lourie, Ian Lowe, Hamdi Mani, Philip Mauskopf, Christopher McKenney, Federico Nati, Giles Novak, Giampaolo Pisano , et al. (8 additional authors not shown)

Abstract: The Next Generation Balloon-Borne Large Aperture Submillimeter Telescope (BLAST-TNG) was a unique instrument for characterizing the polarized submillimeter sky at high-angular resolution. BLAST-TNG flew from the Long Duration Balloon Facility in Antarctica in January 2020. Despite the short flight duration, the instrument worked very well and is providing significant information about each subsyst… ▽ More The Next Generation Balloon-Borne Large Aperture Submillimeter Telescope (BLAST-TNG) was a unique instrument for characterizing the polarized submillimeter sky at high-angular resolution. BLAST-TNG flew from the Long Duration Balloon Facility in Antarctica in January 2020. Despite the short flight duration, the instrument worked very well and is providing significant information about each subsystem that will be invaluable for future balloon missions. In this contribution, we discuss the performance of telescope and gondola. △ Less

Submitted 2 January, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: Submitted to SPIE Astronomical Telescopes + Instrumentation, Ground-based and Airborne Telescopes VIII

Journal ref: Proc. SPIE 11445, Ground-based and Airborne Telescopes VIII, 2020

Gravity and rotation drag the magnetic field in high-mass star formation

Authors: H. Beuther, J. D. Soler, H. Linz, Th. Henning, C. Gieser, R. Kuiper, W. Vlemmings, P. Hennebelle, S. Feng, R. Smith, A. Ahmadi

Abstract: 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 muc… ▽ More 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. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 10 pages, 4 figures, accepted for the Astrophysical Journal, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html