We show that strong interactions violating the conservation of lepton numbers in the neutrino sector could significantly alter the standard low entropy picture for the pre-supernova collapsing core of a massive star. A rapid neutrino-antineutrino equilibration leads to entropy generation and enhanced electron capture and, hence, a lower electron fraction than in the standard model. This would affect the downstream core evolution, the prospects for a supernova explosion, and the emergent neutrino signal. If realized by lepton-number-violating neutrino self-interactions (LNV $\nu$SI), the relevant mediator mass and coupling ranges can be probed by future accelerator-based experiments.
Xianjin Shen, Hong-Li Liu, Zhiyuan Ren, Anandmayee Tej, Di Li, Hauyu Baobab Liu, Gary A. Fuller, Jinjin Xie, Sihan Jiao, Aiyuan Yang, Patrick M. Koch, Fengwei Xu, Patricio Sanhueza, Pham N. Diep, Nicolas Peretto, Ram K. Yadav, Busaba H. Kramer, Koichiro Sugiyama, Mark Rawlings, Chang Won Lee, et al (15) Filaments are believed to play a key role in high-mass star formation. We present a systematic study of the filaments and their hosting clumps in the G35 molecular complex using JCMT SCUBA-2 850 $\micron$ continuum data. We identified five clouds in the complex and 91 filaments within them, some of which form 10 hub-filament systems (HFSs), each with at least 3 hub-composing filaments. We also compiled a catalogue of 350 dense clumps, 183 of which are associated with the filaments. We investigated the physical properties of the filaments and clumps, such as mass, density, and size, and their relation to star formation. We find that the global mass-length trend of the filaments is consistent with a turbulent origin, while the hub-composing filaments of high line masses ($m_{\rm l}\,>$\u2009230\,$\mathrm{M_{\odot}~pc^{-1}}$) in HFSs deviate from this relation, possibly due to feedback from massive star formation. We also find that the most massive and densest clumps (R\,$>$\u20090.2\u2009pc, M\,$>35\,\mathrm{M_{\odot}}$, $\mathrm{\Sigma}>\,0.05\,\mathrm{g~cm^{-2}}$) are located in the filaments and in the hubs of HFS with the latter bearing a higher probability of occurrence of high-mass star-forming signatures, highlighting the preferential sites of HFSs for high-mass star formation. We do not find significant variation in the clump mass surface density across different evolutionary environments of the clouds, which may reflect the balance between mass accretion and stellar feedback.
Understanding how stars form, evolve and impact molecular clouds is key to understanding why star formation is such an inefficient process globally. In this paper, we use the infrared bright fraction, $f_\text{IRB}$ (the fraction of a given molecular cloud that appears bright against the 8 $\mu$m Milky Way background) as a proxy for time evolution to test how cloud properties change as star formation evolves. We apply this metric to 12,000 high-mass star-forming molecular clouds we identify using the Herschel-Hi-GAL survey between $|l|<70^\circ$ on the Milky Way plane. We find clouds are not static while forming stars. Instead, molecular clouds continuously gain mass while star formation progresses. By performing principal component analysis on the cloud properties, we find that they evolve down two paths distinguished by their mass gain. Most clouds (80%) gain four times more mass as a function of $f_\text{IRB}$. The remaining 20% experience an extreme period of growth, growing in mass by a factor of 150 on average and during this period, they initially gain mass fast enough to outpace their star formation. For all clouds, it is only after half their area becomes star forming that mass loss occurs. We expect stellar feedback and potentially galactic shear is responsible. By analysing cloud positions, we suggest that the rate of mass growth may be linked to the larger galactic environment. Altogether, these results have strong implications on how we assess star forming ability on cloud scales when assuming molecular cloud masses are fixed in time.
K. Pattle, P. S. Barry, A. W. Blain, M. Booth, R. A. Booth, D. L. Clements, M. J. Currie, S. Doyle, D. Eden, G. A. Fuller, M. Griffin, P. G. Huggard, J. D. Ilee, J. Karoly, Z. A. Khan, N. Klimovich, E. Kontar, P. Klaassen, A. J. Rigby, P. Scicluna, et al (14) In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre and millimetre community to determine their key priorities for both the near-term and long-term future of the field. We further performed detailed reviews of UK leadership in submillimetre/millimetre science and instrumentation. Our key strategic priorities are as follows: 1. The UK must be a key partner in the forthcoming AtLAST telescope, for which it is essential that the UK remains a key partner in the JCMT in the intermediate term. 2. The UK must maintain, and if possible enhance, access to ALMA and aim to lead parts of instrument development for ALMA2040. Our strategic priorities complement one another: AtLAST (a 50m single-dish telescope) and an upgraded ALMA (a large configurable interferometric array) would be in synergy, not competition, with one another. Both have identified and are working towards the same overarching science goals, and both are required in order to fully address these goals.
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, et al (30) 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.
V. S. Veena, W.-J. Kim, Alvaro Sanchez-Monge, P. Schilke, K. M. Menten, G. A. Fuller, M. C. Sormani, F. Wyrowski, W. E. Banda-Barragan, D. Riquelme, P. Tarrio, P. de Vicente The expanding molecular ring (EMR) manifests itself as a parallelogram in the position-velocity diagram of spectral line emission from the Central Molecular Zone (CMZ) surrounding the Galacic centre (GC). Using multiwavelength data, we investigate the gas kinematics, star formation activity, and the presence of shocked gas in a 200 pc long high velocity gas stream (V~ +150 km/s) with a double helix morphology named the helix stream, that is located 15-55 pc above the CMZ and is kinematically associated with the EMR/parallelogram. We carried out molecular line observations using the IRAM 30m, Yebes 40m, and APEX 12m telescopes. The detection of four rotational transitions of the SiO molecule indicate the presence of shocks. We derived the SiO column densities and abundances in different regions of the helix stream. The presence of protostellar clumps and a candidate HII region signify the ongoing star formation activity within the helix stream. The cloud is massive (2.5x10^6 M_sun) and highly turbulent. We find evidence of cloud-cloud collisions towards the eastern edge (l~1.3\deg), suggesting a dynamic interaction with the CMZ. An expanding shell is detected within the cloud with radius of 6.7 pc and an expansion velocity of 35 km/s. The shell might be powered by several supernovae or a single hypernova. The SiO abundance within the helix stream implies extensive shock processes occurring on large scales. The helical or cork-screw velocity structure of the helix stream indicates twisting and turning motions within the cloud. We propose that the helix stream is the continuation of the near side bar lane, that is overshooting after brushing the CMZ. Our findings carry profound implications for understanding star formation in extreme conditions and elucidate the intricate properties of gas and dust associated with nuclear inflows in barred spiral galaxies.
We calculate the extent to which collisionless dark matter impacts the stability of supermassive stars $(M\gtrsim10^4\,M_\odot)$. We find that, depending on the star's mass, a dark matter content in excess of ${\sim}1\%$ by mass throughout the entire star can raise the critical central density for the onset general relativistic instability, in some cases by orders of magnitude. We consider implications of this effect for the onset of nuclear burning and significant neutrino energy losses.
To study the impact of the initial effects of metallicity (i.e., the abundance of elements heavier than helium) on star formation and the formation of different molecular species, we searched for hot molecular cores in the sub-solar metallicity environment of the Large Magellanic Cloud (LMC). We conducted an ALMA Band 6 observations of 20 fields centered on young stellar objects (YSOs) distributed over the LMC in order to search for hot molecular cores in this galaxy. We detected a total of 65 compact 1.2 mm continuum cores in the 20 ALMA fields and analyzed their spectra with XCLASS software. The main temperature tracers are CH3OH and SO2, with more than two transitions detected in the observed frequency ranges. Other molecular lines with high detection rates in our sample are CS , SO, H13CO+, H13CN, HC15 N, and SiO. More complex molecules, such as HNCO, HDCO, HC3N, CH3CN, and NH2CHO, and multiple transitions of SO and SO2 isotopologues showed tentative or definite detection toward a small subset of the cores. According to the chemical richness of the cores and high temperatures from the XCLASS fitting, we report the detection of four hot cores and one hot core candidate. With one new hot core detection in this study, the number of detected hot cores in the LMC increases to seven. Six out of seven hot cores detected in the LMC to date are located in the stellar bar region of this galaxy. These six hot cores show emission from complex organic molecules (COMs), such as CH3OH, CH3CN, CH3OCHO, and CH3OCH3. The only known hot core in the LMC with no detection of COMs is located outside the bar region. The metallicity in the LMC presents a shallow gradient increasing from outer regions toward the bar. We suggest that the formation of hot molecular cores containing COMs ensues from the new generation of stars forming in the more metal-rich environment of the LMC bar.
D. J. Eden, Tie Liu, T.J.T. Moore, J. Di Francesco, G. Fuller, Kee-Tae Kim, Di Li, S.-Y. Liu, R. Plume, Ken'ichi Tatematsu, M.A. Thompson, Y. Wu, L. Bronfman, H.M. Butner, M.J. Currie, G. Garay, P.F. Goldsmith, N. Hirano, D. Johnstone, M. Juvela, et al (15) We have investigated the physical properties of Planck Galactic Cold Clumps (PGCCs) located in the Galactic Plane, using the JCMT Plane Survey (JPS) and the SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. By utilising a suite of molecular-line surveys, velocities and distances were assigned to the compact sources within the PGCCs, placing them in a Galactic context. The properties of these compact sources show no large-scale variations with Galactic environment. Investigating the star-forming content of the sample, we find that the luminosity-to-mass ratio (L/M) is an order of magnitude lower than in other Galactic studies, indicating that these objects are hosting lower levels of star formation. Finally, by comparing ATLASGAL sources that are associated or are not associated with PGCCs, we find that those associated with PGCCs are typically colder, denser, and have a lower L/M ratio, hinting that PGCCs are a distinct population of Galactic Plane sources.
This study presents a detailed analysis of the GAL045.804-0.356 massive star-forming clump. A high-angular resolution and sensitivity observations were conducted using MeerKAT at 1.28 GHz and ALMA interferometer at 1.3 mm. Two distinct centimetre radio continuum emissions (source A and source B) were identified within the clump. A comprehensive investigation was carried out on source A, the G45.804-0.355 star-forming region (SFR) due to its association with Extended Green Object (EGO), 6.7 GHz methanol maser and the spatial coincidence with the peak of the dust continuum emission at 870 $\mu$m. The ALMA observations revealed seven dense dust condensations (MM1 to MM7) in source A. The brightest ($S_{\rm \nu} \sim$ 87 mJy) and massive main dense core, MM1, was co-located with the 6.7 GHz methanol maser. Explorations into the kinematics revealed gas motions characterised by a velocity gradient across the MM1 core. Furthermore, molecular line emission showed the presence of an extended arm-like structure, with a physical size of 0.25 pc $\times$ 0.18 pc ($\sim$ 50000 au $\times$ 30000 au) at a distance of 7.3 kpc. Amongst these arms, two arms were prominently identified in both the dust continuum and some of the molecular lines. A blue-shifted absorption P-Cygni profile was seen in the H$_2$CO line spectrum. The findings of this study are both intriguing and new, utilising data from MeerKAT and ALMA to investigate the characteristics of the AGAL45 clump. The evidence of spiral arms, the compact nature of the EGO and $<$ 2 km s$^{-1}$ velocity gradient are all indicative of G45.804-0.355 being oriented face-on.
Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (C$\nu$B) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The C$\nu$B emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature range $T\sim 10\,{\rm MeV}$ to $T\sim10\,{\rm keV}$. This is the epoch where beyond Standard Model (BSM) physics may be influential in setting the light element abundances and the necessarily distorted fossil neutrino energy spectra. Here we use techniques honed in extensive CMB studies to analyze the C$\nu$B as calculated in detailed neutrino energy transport and nuclear reaction simulations. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and primordial abundance measurements to provide new insights into the C$\nu$B and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the C$\nu$B throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. The success of our compact description within the Standard Model motivates its use also in BSM scenarios. We demonstrate how observations of primordial light element abundances can be used to place constraints on the C$\nu$B energy spectrum, deriving response functions that can be applied for general C$\nu$B spectral distortions. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the C$\nu$B.
We analyse the radio-to-submillimetre spectral energy distribution (SED) for the central pseudobulge of NGC~1365 using archival data from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Very Large Array (VLA). This analysis shows that free-free emission dominates the continuum emission at 50--120~GHz and produces about 75 per cent of the 103~GHz continuum emission. However, the fraction of 103~GHz continuum emission originating from free-free emission varies significantly among different subregions in the pseudobulge, particularly for an outflow from the AGN on the eastern pseudobulge where the synchrotron emission produces half of the 103~GHz continuum emission. Free-free emission also dominates at 103~GHz within the central 400 pc diameter region, but this emission is associated with the AGN rather than star formation. The star formation rate (SFR) within the pseudobulge derived from the ALMA free-free emission is $8.9 \pm 1.1$~M$_\odot$~yr$^{-1}$. This is comparable to the SFR from the mid-infrared emission but higher than the SFR from the extinction-corrected H$\alpha$ line emission, mainly because the pseudobulge is heavily dust obscured. The 1.5 GHz emission yields a comparable SFR for the pseudobulge but may have lower SFRs within subregions of the pseudobulge because of the diffusion outside of these regions of the electrons producing the synchrotron radiation. We propose that applying a correction factor of 75 per cent to the 80--110~GHz continuum emission could provide valuable estimates of the free-free emission without performing any SED decomposition, which could derive extinction-free SFRs within 20 per cent accuracy.
We present a mechanism for producing a cosmologically-significant relic density of one or more sterile neutrinos. This scheme invokes two steps: First, a population of "heavy" sterile neutrinos is created by scattering-induced decoherence of active neutrinos; Second, this population is transferred, via sterile neutrino self-interaction-mediated scatterings and decays, to one or more lighter mass ($\sim 10\,{\rm keV}$ to $\sim 1\,{\rm GeV}$) sterile neutrinos that are far more weakly (or not at all) mixed with active species and could constitute dark matter. Dark matter produced this way can evade current electromagnetic and structure-based bounds, but may nevertheless be probed by future observations.
Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, et al (137) We report 850 $\mu$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity.
High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1 pc-scale observations from NOEMA and the IRAM 30-m telescope to trace dense cores with 2.8 mm continuum, and gas kinematics in C$^{18}$O, HCO$^+$, HNC, and N$_2$H$^+$ ($J$=1$-$0). We supplement our continuum sample with six IRDCs observed at 2.9 mm with ALMA, and examine the relationships between core- and clump-scale properties. We have developed a fully-automated multiple-velocity component hyperfine line-fitting code called mwydyn which we employ to trace the dense gas kinematics in N$_2$H$^+$ (1$-$0), revealing highly complex and dynamic clump interiors. We find that parsec-scale clump mass is the most important factor driving the evolution; more massive clumps are able to concentrate more mass into their most massive cores - with a log-normally distributed efficiency of around 9% - in addition to containing the most dynamic gas. Distributions of linewidths within the most massive cores are similar to the ambient gas, suggesting that they are not dynamically decoupled, but are similarly chaotic. A number of studies have previously suggested that clumps are globally collapsing; in such a scenario, the observed kinematics of clump centres would be the direct result of gravity-driven mass inflows that become ever more complex as the clumps evolve, which in turn leads to the chaotic mass growth of their core populations.
The neutrinos in the diffuse supernova neutrino background (DSNB) travel over cosmological distances and this provides them with an excellent opportunity to interact with dark relics. We show that a cosmologically-significant relic population of keV-mass sterile neutrinos with strong self-interactions could imprint their presence in the DSNB. The signatures of the self-interactions would be ``dips" in the otherwise smooth DSNB spectrum. Upcoming large-scale neutrino detectors, for example Hyper-Kamiokande, have a good chance of detecting the DSNB and these dips. If no dips are detected, this method serves as an independent constraint on the sterile neutrino self-interaction strength and mixing with active neutrinos. We show that relic sterile neutrino parameters that evade X-ray and structure bounds may nevertheless be testable by future detectors like TRISTAN, but may also produce dips in the DSNB which could be detectable. Such a detection would suggest the existence of a cosmologically-significant, strongly self-interacting sterile neutrino background, likely embedded in a richer dark sector.
The application of silicon monoxide (SiO) as a shock tracer arises from its propensity to occur in the gas phase as a result of shock-induced phenomena, including outflow activity and interactions between molecular clouds and expanding HII regions or supernova remnants. We searched for indications of shocks toward 366 massive star-forming regions by observing the ground rotational transition of SiO ($v=0$, $J=1-0$) at 43 GHz with the Korean VLBI Network (KVN) 21 m telescopes to extend our understanding on the origins of SiO in star-forming regions. We detected SiO emission toward 104 regions that consist of 57 IRDCs, 21 HMPOs, and 26 UCHIIs. The determined median SiO column density, $N$(SiO), and abundance, $X$(SiO), relative to $N$(H$_2$) are $8.12\times10^{12}$ cm$^{-2}$ and $1.28\times10^{-10}$, respectively. These values are similar to those obtained toward other star-forming regions and also consistent with predicted values from shock models with low-velocity shocks ($\lesssim$10 - 15 km s$^{-1}$). While the $X$(SiO) does not exhibit any strong correlation with the evolutionary stages of their host clumps, $L_{\rm SiO}$ is highly correlated with dust clump mass, and $L_{\rm SiO}/L_{\rm bol}$ also has a strong negative correlation with $T_{\rm dust}$. This shows that colder and younger clumps have high $L_{\rm SiO}/L_{\rm bol}$ suggestive of an evolutionary trend. This trend is not due to excess emission at higher velocities, such as SiO wing features, as the colder sources with high $L_{\rm SiO}/L_{\rm bol}$ ratios lack wing features. Comparing SiO emission with H$_2$O and Class I CH$_3$OH masers, we find a significant correlation between $L_{\rm SiO}$/$L_{\rm bol}$ and $L_{\rm CH_3OH}/L_{\rm bol}$ ratios, whereas no similar correlation is seen for the H$_2$O maser emission. This suggests a similar origin for the SiO and Class I CH$_3$OH emission in these sources.
A. Avison, G. A. Fuller, N. Asabre Frimpong, S. Etoka, M. Hoare, B.M. Jones, N. Peretto, A. Traficante, F. van der Tak, J.E. Pineda, M. Beltrán, F. Wyrowski, M. Thompson, S. Lumsden, Z. Nagy, T. Hill, S. Viti, F. Fontani, P. Schilke The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regions observed with ALMA at 1.3mm and explore the fragmentation, clustering and flux density properties of the sample. The TEMPO sample fields are each found to contain multiple fragments (between 2-15 per field). The flux density budget is split evenly (53%-47%) between fields where emission is dominated by a single high flux density fragment and those in which the combined flux density of fainter objects dominates. The fragmentation scales observed in most fields are not comparable with the thermal Jeans length, $\lambda_J$, being larger in the majority of cases, suggestive of some non-thermal mechanism. A tentative evolutionary trend is seen between luminosity of the clump and the `spectral line richness' of the TEMPO fields; with 6.7GHz maser associated fields found to be lower luminosity and more line rich. This work also describes a method of line-free continuum channel selection within ALMA data and a generalised approach used to distinguishing sources which are potentially star-forming from those which are not, utilising interferometric visibility properties.
We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 85.69 and 99.02 GHz continuum emission and H42$\alpha$ and H40$\alpha$ lines emission from the central 1~kpc of NGC 1808. These forms of emission are tracers of photoionizing stars but unaffected by dust obscuration that we use to test the applicability of other commonly star formation metrics. An analysis of the spectral energy distributions shows that free-free emission contributes about 60 to 90 per cent of the continuum emission in the 85-100 GHz frequency range, dependent on the region. The star formation rate (SFR) derived from the ALMA free-free emission is $3.1\pm0.3$~M$_\odot$~yr$^{-1}$. This is comparable to the SFRs measured from the infrared emission, mainly because most of the bolometric energy from the heavily obscured region is emitted as infrared emission. The radio 1.5~GHz emission yields a SFR 25 per cent lower than the ALMA value, probably because of the diffusion of the electrons producing the synchrotron emission beyond the star-forming regions. The SFRs measured from the extinction-corrected H$\alpha$ line emission are about 40 to 65 per cent of the SFR derived from the ALMA data, likely because this metric was not calibrated for high extinction regions. Some SFRs based on extinction-corrected ultraviolet emission are similar to those from ALMA and infrared data, but given that the ultraviolet terms in the extinction correction equations are very small, these metrics seem inappropriate to apply to this dusty starburst.
The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material toward the filaments within Barnard 5 (B5) and the infall from the envelope to the protostellar disk of the embedded protostar B5-IRS1. Our goal is to follow the flow of material from the larger, dense core scale, to the protostellar disk scale. We present new HC$_3$N line data from the NOEMA and 30m telescopes covering the coherence zone of B5, together with ALMA H$_2$CO and C$^{18}$O maps toward the protostellar envelope. We fit multiple Gaussian components to the lines so as to decompose their individual physical components. We investigate the HC$_3$N velocity gradients to determine the direction of chemically-fresh gas flow. At envelope scales, we use a clustering algorithm to disentangle the different kinematic components within H$_2$CO emission. At dense core scales, HC$_3$N traces the infall from the B5 region toward the filaments. HC$_3$N velocity gradients are consistent with accretion toward the filament spines plus flow along them. We found a $\sim2800$ au streamer in H$_2$CO emission which is blueshifted with respect to the protostar and deposits gas at outer disk scales. The strongest velocity gradients at large scales curve toward the position of the streamer at small scales, suggesting a connection between both flows. Our analysis suggests that the gas can flow from the dense core to the protostar. This implies that the mass available for a protostar is not limited to its envelope, and can receiving chemically-unprocessed gas after the main accretion phase.
We report the discovery of a velocity coherent, funnel shaped ^13CO emission feature in the Galactic centre (GC) using data from the SEDIGISM survey. The molecular cloud appears as a low velocity structure (V_LSR=[-3.5, +3.5] km/s) with an angular extent of 0.95\deg x 1\deg, extending toward positive Galactic latitudes. The structure is offset from Sgr A* toward negative Galactic longitudes and spatially and morphologically correlates well with the northern lobe of the 430 pc radio bubble, believed to be the radio counterpart of the multiwavelength GC chimney. Spectral line observations in the frequency range of 85-116 GHz have been carried out using the IRAM 30 metre telescope toward 12 positions along the funnel-shaped emission. We examine the ^12C/^13C isotopic ratios using various molecules and their isotopologues. The mean ^12C/^13C isotope ratio (30.6+-2.9) is consistent with the structure located within inner 3 kpc of the Galaxy and possibly in the GC. The velocity of the molecular funnel is consistent with previous radio recombination line measurements of the northern lobe of radio bubble. Our multiwavelength analysis suggests that the funnel shaped structure extending over 100 pc above the Galactic plane is the molecular counterpart of the northern GC chimney.
Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, et al (138) We present observations of polarized dust emission at 850 $\mu$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $\mu$G in the main starless core and up to $\sim$90$\pm$40 $\mu$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores.
The formation of stellar clusters dictates the pace at which galaxies evolve, and solving the question of their formation will undoubtedly lead to a better understanding of the Universe as a whole. While it is well known that star clusters form within parsec-scale over-densities of interstellar molecular gas called clumps, it is, however, unclear whether these clumps represent the high-density tip of a continuous gaseous flow that gradually leads towards the formation of stars, or a transition within the gas physical properties. Here, we present a unique analysis of a sample of 27 infrared dark clouds embedded within 24 individual molecular clouds that combine a large set of observations, allowing us to compute the mass and velocity dispersion profiles of each, from the scale of tens of parsecs down to the scale of tenths of a parsec. These profiles reveal that the vast majority of the clouds, if not all, are consistent with being self-gravitating on all scales, and that the clumps, on parsec-scale, are often dynamically decoupled from their surrounding molecular clouds, exhibiting steeper density profiles ($\rho\propto r^{-2}$) and flat velocity dispersion profiles ($\sigma\propto r^0$), clearly departing from Larson's relations. These findings suggest that the formation of star clusters correspond to a transition regime within the properties of the self-gravitating molecular gas. We propose that this transition regime is one that corresponds to the gravitational collapse of parsec-scale clumps within otherwise stable molecular clouds.
Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, et al (138) We present BISTRO Survey 850 \mum dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965).
As the early universe expands and cools the rates of the weak interactions that keep neutrinos in thermal equilibrium with the matter and the related rates of the reactions that inter-convert neutrons and protons decrease. Eventually, these rates fall below the expansion rate -- they freeze out. Likewise, the rates of the strong and electromagnetic nuclear reactions that build up and tear down nuclei, though fast enough to maintain equilibrium early on, slow down and ultimately lead to freeze out. Together these freeze out processes comprise the epoch of Big Bang Nucleosynthesis (BBN). The relics emerging from this early time include the light element abundances, for example of helium and deuterium, and a background of decoupled neutrinos, a "C$\nu$B" , roughly analogous to the Cosmic Microwave Background, the CMB. These fossil relics encode the history of the physics operating in the early universe. Consequently, BBN has emerged as a key tool for constraining new, beyond-standard-model (BSM) physics. BBN may become an even finer probe of BSM physics, given the anticipated higher precision in measurements of the primordial abundances of deuterium and helium afforded by the advent of large optical telescopes and Stage-4 CMB experiments. The latter experiments will also provide higher precision determinations of $N_{\rm eff}$, a measure of the relativistic energy density at the photon decoupling epoch and, hence, an important probe of the C$\nu$B.
A. Traficante, B. M. Jones, A. Avison, G. A. Fuller, M. Benedettini, D. Elia, S. Molinari, N. Peretto, S. Pezzuto, T. Pillai, K. L. J. Rygl, E. Schisano, R. J. Smith The formation mechanism of the most massive stars is far from completely understood. It is still unclear if the formation is core-fed or clump-fed, i.e. if the process is an extension of what happens in low-mass stars, or if the process is more dynamical such as a continuous, multi-scale accretion from the gas at parsec (or even larger) scales. In this context we introduce the SQUALO project, an ALMA 1.3 mm and 3 mm survey designed to investigate the properties of 13 massive clumps selected at various evolutionary stages, with the common feature that they all show evidence for accretion at the clump scale. In this work we present the results obtained from the 1.3 mm continuum data. Our observations identify 55 objects with masses in the range 0.4 <~ M <~ 309 M_sun, with evidence that the youngest clumps already present some degree of fragmentation. The data show that physical properties such as mass and surface density of the fragments and their parent clumps are tightly correlated. The minimum distance between fragments decreases with evolution, suggesting a dynamical scenario in which massive clumps first fragment under the influence of non-thermal motions driven by the competition between turbulence and gravity. With time gravitational collapse takes over and the fragments organize themselves into more thermally supported objects while continuing to accrete from their parent clump. Finally, one source does not fragment, suggesting that the support of other mechanisms (such as magnetic fields) is crucial only in specific star-forming regions.
Feng-Wei Xu, Ke Wang, Tie Liu, Paul F. Goldsmith, Qizhou Zhang, Mika Juvela, Hong-Li Liu, Sheng-Li Qin, Guang-Xing Li, Anandmayee Tej, Guido Garay, Leonardo Bronfman, Shanghuo Li, Yue-Fang Wu, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Ken'ichi Tatematsu, Zhiyuan Ren, Yong Zhang, L. Viktor Toth, et al (20) We present ALMA Band-3/7 observations towards "the Heart" of a massive hub-filament system (HFS) SDC335, to investigate its fragmentation and accretion. At a resolution of $\sim0.03$ pc, 3 mm continuum emission resolves two massive dense cores MM1 and MM2, with $383(^{+234}_{-120})$ $M_\odot$ (10-24% mass of "the Heart") and $74(^{+47}_{-24})$ $M_\odot$, respectively. With a resolution down to 0.01 pc, 0.87 mm continuum emission shows MM1 further fragments into six condensations and multi-transition lines of H$_2$CS provide temperature estimation. The relation between separation and mass of condensations at a scale of 0.01 pc favors turbulent Jeans fragmentation where the turbulence seems to be scale-free rather than scale-dependent. We use the H$^{13}$CO$^+$ (1-0) emission line to resolve the complex gas motion inside "the Heart" in position-position-velocity space. We identify four major gas streams connected to large-scale filaments, inheriting the anti-clockwise spiral pattern. Along these streams, gas feeds the central massive core MM1. Assuming an inclination angle of $45(\pm15)^{\circ}$ and a H$^{13}$CO$^+$ abundance of $5(\pm3)\times10^{-11}$, the total mass infall rate is estimated to be $2.40(\pm0.78)\times10^{-3}$ $M_\odot$ yr$^{-1}$, numerically consistent with the accretion rates derived from the clump-scale spherical infall model and the core-scale outflows. The consistency suggests a continuous, near steady-state, and efficient accretion from global collapse, therefore ensuring core feeding. Our comprehensive study of SDC335 showcases the detailed gas kinematics in a prototypical massive infalling clump and calls for further systematic and statistical analyses in a large sample.
T. Irabor, M.G. Hoare, M. Burton, W.D. Cotton, P. Diamond, S. Dougherty, S.P. Ellingsen, R. Fender, G.A. Fuller, S. Garrington, P.F. Goldsmith, J. Green, A.G. Gunn, J. Jackson, S. Kurtz, S.L. Lumsden J. Marti, I. McDonald, S. Molinari, T.J. Moore, M. Mutale, et al (14) We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$\sigma$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$\mu m$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane.
Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, et al (125) We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines.
Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, et al (137) We present 850 $\mu$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments.
Beyond Standard Model extensions of QCD could result in quark and gluon confinement occurring well above a temperature of $\sim$GeV. These models can also alter the order of the QCD phase transition. The enhanced production of primordial black holes (PBHs) that can accompany the change in relativistic degrees of freedom at the QCD transition therefore could favor the production of PBHs with mass scales smaller than the Standard Model QCD horizon scale. Consequently, and unlike PBHs associated with a standard GeV-scale QCD transition, such PBHs can account for all the dark matter abundance in the unconstrained asteroid-mass window. This links beyond Standard Model modifications of QCD physics over a broad range of unexplored temperature regimes ($\sim 10-10^3$ TeV) with microlensing surveys searching for PBHs. Additionally, we discuss implications of these models for gravitational wave experiments. We show that a first order QCD phase transition at $\sim7$ TeV is consistent with the Subaru Hyper-Suprime Cam candidate event, while a $\sim 70$ GeV transition is consistent with OGLE candidate events, and also could account for the claimed NANOGrav gravitational wave signal.
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 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^{-\Gamma}$, $\Gamma \simeq 0$. Abridged
Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, et al (140) We present and analyze observations of polarized dust emission at 850 $\mu$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical.
B. Westbrook, P. A. R. Ade, M. Aguilar, Y. Akiba, K. Arnold, C. Baccigalupi, D. Barron, D. Beck, S. Beckman, A. N. Bender, F. Bianchini, D. Boettger, J. Borrill, S. Chapman, Y. Chinone, G. Coppi, K. Crowley, A. Cukierman, T. de, R. Dünner, et al (73) We present on the status of POLARBEAR-2 A (PB2-A) focal plane fabrication. The PB2-A is the first of three telescopes in the Simon Array (SA), which is an array of three cosmic microwave background (CMB) polarization sensitive telescopes located at the POLARBEAR (PB) site in Northern Chile. As the successor to the PB experiment, each telescope and receiver combination is named as PB2-A, PB2-B, and PB2-C. PB2-A and -B will have nearly identical receivers operating at 90 and 150 GHz while PB2-C will house a receiver operating at 220 and 270 GHz. Each receiver contains a focal plane consisting of seven close-hex packed lenslet coupled sinuous antenna transition edge sensor bolometer arrays. Each array contains 271 di-chroic optical pixels each of which have four TES bolometers for a total of 7588 detectors per receiver. We have produced a set of two types of candidate arrays for PB2-A. The first we call Version 11 (V11) and uses a silicon oxide (SiOx) for the transmission lines and cross-over process for orthogonal polarizations. The second we call Version 13 (V13) and uses silicon nitride (SiNx) for the transmission lines and cross-under process for orthogonal polarizations. We have produced enough of each type of array to fully populate the focal plane of the PB2-A receiver. The average wirebond yield for V11 and V13 arrays is 93.2% and 95.6% respectively. The V11 arrays had a superconducting transition temperature (Tc) of 452 +/- 15 mK, a normal resistance (Rn) of 1.25 +/- 0.20 Ohms, and saturations powers of 5.2 +/- 1.0 pW and 13 +/- 1.2 pW for the 90 and 150 GHz bands respectively. The V13 arrays had a superconducting transition temperature (Tc) of 456 +/-6 mK, a normal resistance (Rn) of 1.1 +/- 0.2 Ohms, and saturations powers of 10.8 +/- 1.8 pW and 22.9 +/- 2.6 pW for the 90 and 150 GHz bands respectively.
In the multi-scale view of the star formation process the material flows from large molecular clouds down to clumps and cores. In this paradigm it is still unclear if it is gravity or turbulence that drives the observed supersonic non-thermal motions during the collapse, in particular in high-mass regions, and at which scales gravity becomes eventually dominant over the turbulence of the interstellar medium. To investigate this problem we have combined the dynamics of a sample of 70 micron-quiet clumps, selected to cover a wide range of masses and surface densities, with the dynamics of the parent filaments in which they are embedded. We observe a continuous interplay between turbulence and gravity, where the former creates structures at all scales and the latter takes the lead when a critical value of the surface density is reached, Sigma_th = 0.1 g cm^-2. In the densest filaments this transition can occur at the parsec, or even larger scales, leading to a global collapse of the whole region and most likely to the formation of the massive objects.
P. Nazari, J. D. Meijerhof, M. L. van Gelder, A. Ahmadi, E. F. van Dishoeck, B. Tabone, D. Langeroodi, N. F. W. Ligterink, J. Jaspers, M. T. Beltrán, G. A. Fuller, Á. Sánchez-Monge, P. Schilke No statistical study of COMs toward a large sample of high-mass protostars with ALMA has been carried out so far. We aim to study six N-bearing species: CH$_3$CN, HNCO, NH$_2$CHO, C$_2$H$_5$CN, C$_2$H$_3$CN and CH$_3$NH$_2$ in a large sample of high-mass protostars. From the ALMAGAL survey, 37 of the most line-rich hot molecular cores are selected. Next, we fit their spectra and find column densities and excitation temperatures of the above N-bearing species, in addition to CH$_3$OH. We (tentatively) detect CH$_3$NH$_2$ in $\sim32%$ of the sources. We find three groups of species when comparing their excitation temperatures: hot (NH$_2$CHO; Tex > 250 K), warm (C$_2$H$_3$CN, HN$^{13}$CO and CH$_{3}^{13}$CN; 100 K < Tex < 250 K) and cold species (CH$_3$OH and CH$_3$NH$_2$; Tex < 100 K). This temperature segregation reflects the trend seen in their sublimation temperatures and validates the idea of onion-like structure of COMs around protostars. Moreover, the molecules studied here show constant column density ratios across low- and high-mass protostars with scatter less than a factor $\sim3$ around the mean. The constant column density ratios point to a common formation environment of COMs or their precursors, most likely in the pre-stellar ices. The scatter around the mean of the ratios, although small, varies depending on the species considered. This spread can either have a physical origin (source structure, line or dust optical depth) or a chemical one. Formamide is most prone to the physical effects as it is tracing the closest regions to the protostars, whereas such effects are small for other species. Assuming that all molecules form in the pre-stellar ices, the scatter variations could be explained by differences in lifetimes or physical conditions of the pre-stellar clouds. If the pre-stellar lifetimes are the main factor, they should be similar for low- and high-mass protostars.
The physics of the mysterious and stealthy neutrino is at the heart of many phenomena in the cosmos. These particles interact with matter and with each other through the aptly named weak interaction. At typical astrophysical energies the weak interaction is some twenty orders of magnitude weaker than the electromagnetic interaction. However, in the early universe and in collapsing stars neutrinos can more than make up for their feeble interaction strength with huge numbers. Neutrinos can dominate the dynamics in these sites and set the conditions that govern the synthesis of the elements. Here we journey through the history of the discovery of these particles and describe their role in stellar evolution and collapse, the big bang, and multi-messenger astrophysics. Neutrino physics is at the frontier of elementary particle physics, nuclear physics, astrophysics and cosmology. All of these fields overlap in the neutrino story.
The deuteration of molecules forming in the ices such as methanol (CH$_3$OH) is sensitive to the physical conditions during their formation in dense cold clouds and can be probed through observations of deuterated methanol in hot cores. Observations with ALMA containing transitions of CH$_3$OH, CH$_2$DOH, CHD$_2$OH, $^{13}$CH$_3$OH, and CH$_3^{18}$OH are investigated. The column densities of CH$_2$DOH, CHD$_2$OH, and CH$_3$OH are determined for all sources, where the column density of CH$_3$OH is derived from optically thin $^{13}$C and $^{18}$O isotopologues. Consequently, the D/H ratio of methanol is derived taking statistical effects into account. Singly deuterated methanol (CH$_2$DOH) is detected toward 25 of the 99 sources in our sample of the high-mass protostars. Including upper limits, the $\rm (D/H)_{CH_3OH}$ ratio inferred from $N_\mathrm{CH_2DOH}/N_\mathrm{CH_3OH}$ was derived for 38 of the 99 sources and varies between $\sim10^{-3}-10^{-2}$. Including other high-mass hot cores from the literature, the mean methanol D/H ratio is $1.1\pm0.7\times10^{-3}$. This is more than one order of magnitude lower than what is seen for low-mass protostellar systems ($2.2\pm1.2\times10^{-2}$). Doubly deuterated methanol (CHD$_2$OH) is detected toward 11 of the 99 sources. Including upper limits for 15 sources, the $\rm (D/H)_{CH_2DOH}$ ratios derived from $N_\mathrm{CHD_2OH}/N_\mathrm{CH_2DOH}$ are more than two orders of magnitude higher than $\rm (D/H)_{CH_3OH}$ with an average of $2.0\pm0.8\times10^{-1}$ which is similar to what is found for low-mass sources. Comparison with literature GRAINOBLE models suggests that the high-mass prestellar phases are either warm ($>20$ K) or live shorter than the free-fall timescale. In contrast, for low-mass protostars, both a low temperature of $<15$ K and a prestellar phase timescale longer than the free-fall timescale are necessary.
The distribution of ultraviolet (UV) radiation field provides critical constraints on the physical environments of molecular clouds. Within 1 kpc of our solar system and fostering protostars of different masses, the giant molecular clouds in the Gould Belt present an excellent opportunity to resolve the UV field structure in star forming regions. We performed spectral energy distribution (SED) fitting of the archival data from the Herschel Gould Belt Survey (HGBS). Dust radiative transfer analysis with the DUSTY code were applied to 23 regions in 14 molecular complexes of the Gould Belt, resulting in the spatial distribution of radiation field in these regions. For 10 of 15 regions with independent measurements of star formation rate, their star formation rate and UV radiation intensity largely conform to a linear correlation found in previous studies.
H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, et al (145) Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
Hub-filament systems (HFSs) are potential sites of protocluster and massive star formation, and play a key role in mass accumulation. We report JCMT POL-2 850 $\mu$m polarization observations toward the massive HFS SDC13. We detect an organized magnetic field near the hub center with a cloud-scale "U-shape" morphology following the western edge of the hub. Together with larger-scale APEX 13CO and PLANCK polarization data, we find that SDC13 is located at the convergent point of three giant molecular clouds (GMCs) along a large-scale, partially spiral-like magnetic field. The smaller "U-shape" magnetic field is perpendicular to the large-scale magnetic field and the converging GMCs. We explain this as the result of a cloud-cloud collision. Within SDC13, we find that local gravity and velocity gradients point toward filament ridges and hub center. This suggests that gas can locally be pulled onto filaments and overall converges to the hub center. A virial analysis of the central hub shows that gravity dominates magnetic and kinematic energy. Combining large- and small-scale analyses, we propose that SDC13 is initially formed from a collision of clouds moving along the large-scale magnetic field. In the post-shock regions, after the initial turbulent energy has dissipated, gravity takes over and starts to drive the gas accretion along the filaments toward the hub center.
Ken'ichi Tatematsu, You-Ting Yeh, Naomi Hirano, Sheng-Yuan Liu, Tie Liu, Somnath Dutta, Dipen Sahu, Neal J. Evans II, Mika Juvela, Hee-Weon Yi, Jeong-Eun Lee, Patricio Sanhueza, Shanghuo Li, David Eden, Gwanjeong Kim, Chin-Fei Lee, Yuefang Wu, Kee-Tae Kim, L. Viktor T'oth, Minho Choi, et al (11) In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $\delta V$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $\delta V$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0).
CMB-HD Collaboration, Simone Aiola, Yashar Akrami, Kaustuv Basu, Michael Boylan-Kolchin, Thejs Brinckmann, Sean Bryan, Caitlin M. Casey, Jens Chluba, Sebastien Clesse, Francis-Yan Cyr-Racine, Luca Di Mascolo, Simon Dicker, Thomas Essinger-Hileman, Gerrit S. Farren, Michael A. Fedderke, Simone Ferraro, George M. Fuller, Nicholas Galitzki, Vera Gluscevic, et al (44) CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.
M. L. van Gelder, P. Nazari, B. Tabone, A. Ahmadi, E. F. van Dishoeck, M. T. Beltrán, G. A. Fuller, N. Sakai, Á. Sánchez-Monge, P. Schilke, Y.-L. Yang, Y. Zhang Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.70\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs.
Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, et al (132) We present 850 $\mu$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~\mu$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.
Ana K. Diaz-Rodriguez, Guillem Anglada, Guillermo Blázquez-Calero, Mayra Osorio, José F. Gómez, Gary A. Fuller, Robert Estalella, José M. Torrelles, Sylvie Cabrit, Luis F. Rodríguez, Charlène Lefèvre, Enrique Macías, Carlos Carrasco-González, Luis A. Zapata, Itziar de Gregorio-Monsalvo, Paul T. P. Ho We present VLA and ALMA observations of the close (0.3" = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii $\sim$12 and $\sim$9 au in dust, and $\sim$30 au in gas) with masses of $\sim$0.004-0.009 $M_{sun}$ for VLA 4A (the western component) and $\sim$0.009-0.030 $M_{sun}$ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms extending $\sim$500 au and a mass of $\sim$0.052 $M_{sun}$ appears to be in the earliest stages of formation. The dust emission is more compact and with a very high optical depth toward VLA 4B, while toward VLA 4A the dust column density is lower, allowing the detection of stronger molecular transitions. We infer rotational temperatures of $\sim$140 K, on scales of $\sim$30 au, across the whole source, and a rich chemistry. Molecular transitions typical of hot corinos are detected toward both protostars, being stronger toward VLA 4A, with several ethylene glycol transitions detected only toward this source. There are clear velocity gradients, that we interpret in terms of infall plus rotation of the circumbinary disk, and purely rotation of the circumstellar disk of VLA 4A. We measured orbital proper motions and determined a total stellar mass of 1 $M_{sun}$. From the molecular kinematics we infer the geometry and orientation of the system, and stellar masses of $\sim$0.26 $M_{sun}$ for VLA 4A and $\sim$0.60 $M_{sun}$ for VLA 4B.
Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified `local' Brink-Axel hypothesis for Gamow-Teller transitions for $pf$-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semi-converged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing to previous calculations of astrophysical rates.
A-Ran Lyo, Jongsoo Kim, Sarah Sadavoy, Doug Johnstone, David Berry, Kate Pattle, Woojin Kwon, Pierre Bastien, Takashi Onaka, James Di Francesco, Ji-Hyun Kang, Ray Furuya, Charles L. H. Hull, Motohide Tamura, Patrick M. Koch, Derek Ward-Thompson, Tetsuo Hasegawa, Thiem Hoang, Doris Arzoumanian, Chang Won Lee, et al (126) We present the results of simultaneous 450 $\mu$m and 850 $\mu$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $\mu$G in the central $\sim$0.12 pc region from 850 $\mu$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $\mu$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $\mu$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $\mu$m and $0.36 \pm 0.04$ at 850 $\mu$m, respectively. We think that the shallow slope at 850 $\mu$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects.
The hot and dense early Universe combined with the promise of high-precision cosmological observations provide an intriguing laboratory for Beyond Standard Model (BSM) physics. We simulate the early Universe to examine the effects of the decay of thermally populated sterile neutrino states into Standard Model products around the time of weak decoupling. These decays deposit a significant amount of entropy into the plasma as well as produce a population of high-energy out-of-equilibrium active neutrinos. As a result, we can constrain these models by their inferred value of $N_{\rm eff}$, the effective number of relativistic degrees of freedom. In this work, we explore a variety of models with $N_{\rm eff}{}$ values consistent with CMB observations, but with vastly different active neutrino spectra which will challenge the standard cosmological model, affect lepton capture rates on free nucleons, and may significantly affect Big Bang Nucleosynthesis (BBN).
S. Martín, J. G. Mangum, N. Harada, F. Costagliola, K. Sakamoto, S. Muller, R. Aladro, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, S. Aalto, E. Behrens, L. Colzi, K. L. Emig, G. A. Fuller, S. García-Burillo, T. R. Greve, C. Henkel, et al (15) We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We focus on the analysis of the spectra extracted from the $15''$ ($\sim255$ pc) resolution ALMA Compact Array data. We model the molecular emission assuming local thermodynamic equilibrium with 78 species detected. Additionally, multiple hydrogen and helium recombination lines are identified. Spectral lines contribute 5 to 36\% of the total emission in frequency bins of 50 GHz. We report the first extragalactic detections of C$_2$H$_5$OH, HOCN, HC$_3$HO, and several rare isotopologues. Isotopic ratios of carbon, oxygen, sulfur, nitrogen and silicon were measure with multiple species. Infrared pumped vibrationaly excited HCN, HNC, and HC$_3$N emission, originating in massive star formation locations, is clearly detected at low resolution, while we do not detect it for HCO$^+$. We suggest high temperature conditions in these regions driving a seemingly "carbon-rich" chemistry which may also explain the observed high abundance of organic species close to those in Galactic hot cores. The $L_{vib}/L_{IR}$ ratio is used as a proxy to estimate a $3\%$ contribution from proto super star cluster to the global infrared emission. Measured isotopic ratios with high dipole moment species agree with those within the central kiloparsec of the Galaxy, while those derived from $\rm^{13}C^{18}O$ are a factor of 5 larger, confirming the existence of multiple ISM components within NGC 253 with different degrees of nucleosynthesis enrichment. ALCHEMI provides a template for early Universe galaxies.
The physical processes behind the transfer of mass from parsec-scale clumps to massive-star-forming cores remain elusive. We investigate the relation between the clump morphology and the mass fraction that ends up in its most massive core (MMC) as a function of infrared brightness, i.e. a clump evolutionary tracer. Using ALMA 12 m and ACA we surveyed 6 infrared-dark hubs in 2.9mm continuum at $\sim$3" resolution. To put our sample into context, we also re-analysed published ALMA data from a sample of 29 high mass-surface density ATLASGAL sources. We characterise the size, mass, morphology, and infrared brightness of the clumps using Herschel and Spitzer data. Within the 6 newly observed hubs, we identify 67 cores, and find that the MMCs have masses between 15-911 $\mathrm{M}_{\odot}$ within a radius of 0.018-0.156 pc. The MMC of each hub contains 3-24% of the clump mass ($f_\mathrm{MMC}$), becoming 5-36% once core masses are normalised to the median core radius. Across the 35 clumps, we find no significant difference in the median $f_\mathrm{MMC}$ values of hub and non-hub systems, likely the consequence of a sample bias. However, we find that $f_\mathrm{MMC}$ is $\sim$7.9 times larger for infrared-dark clumps compared to infrared-bright ones. This factor increases up to $\sim$14.5 when comparing our sample of 6 infrared-dark hubs to infrared-bright clumps. We speculate that hub-filament systems efficiently concentrate mass within their MMC early on during its evolution. As clumps evolve, they grow in mass, but such growth does not lead to the formation of more massive MMCs.
Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings is a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the post-shock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. We also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star (PNS) convection, asymmetric neutrino emission, neutrino absorptions and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.
M. Massardi, F. Stoehr, G. J. Bendo, M. Bonato, J. Brand, V. Galluzzi, F. Guglielmetti, E. Liuzzo, N. Marchili, A. M. S. Richards, K. L. J. Rygl, F. Bedosti, A. Giannetti, M. Stagni, C. Knapic, M. Sponza, G. A. Fuller, T. W. B. Muxlow The Additional Representative Images for Legacy (ARI-L) project is a European Development project for ALMA Upgrade approved by the Joint ALMA Observatory (JAO) and the European Southern Observatory (ESO), started in June 2019. It aims to increase the legacy value of the ALMA Science Archive (ASA) by bringing the reduction level of ALMA data from Cycles 2-4 close to that of data from more recent Cycles processed for imaging with the ALMA Pipeline. As of mid-2021 more than 150000 images have been returned to the ASA for public use. At its completion in 2022, the project will have provided enhanced products for at least 70% of the observational data from Cycles 2-4 processable with the ALMA Pipeline. In this paper we present the project rationale, its implementation, and the new opportunities offered to ASA users by the ARI-L products. The ARI-L cubes and images complement the much limited number of archival image products generated during the data quality assurance stages (QA2), which cover only a small fraction of the available data for those Cycles. ARI-L imaging products are highly relevant for many science cases and significantly enhance the possibilities for exploiting archival data. Indeed, ARI-L products facilitate archive access and data usage for science purposes even for non-expert data miners, provide a homogeneous view of all data for better dataset comparisons and download selections, make the archive more accessible to visualization and analysis tools, and enable the generation of preview images and plots similar to those possible for subsequent Cycles.
Ken'ichi Tatematsu, Gwanjeong Kim, Tie Liu, Neal J. Evans II, Hee-Weon Yi, Jeong-Eun Lee, Yuefang Wu, Naomi Hirano, Sheng-Yuan Liu, Somnath Dutta, Dipen Sahu, Patricio Sanhueza, Kee-Tae Kim, Mika Juvela, L. Viktor T'oth, Orsolya Feh'er, Jinhua He, J. X. Ge, Siyi Feng, Minho Choi, et al (18) We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$\mu$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$\mu$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.
We investigate the kinematics of high mass protostellar objects within the high mass star forming region IRAS 19410+2336. We performed high angular resolution observations of 6.7-GHz methanol and 22 GHz water masers using the MERLIN (Multi-Element Radio Linked Interferometer Network) and e-MERLIN interferometers. The 6.7-GHz methanol maser emission line was detected within the $\sim$ 16--27 km s$^{-1}$ velocity range with a peak flux density $\sim$50 Jy. The maser spots are spread over $\sim$1.3 arcsec on the sky, corresponding to $\sim$2800 au at a distance of 2.16 kpc. These are the first astrometric measurements at 6.7 GHz in IRAS 19410+2336. The 22-GHz water maser line was imaged in 2005 and 2019 (the latter with good astrometry). Its velocities range from 13 to $\sim$29 km s$^{-1}$. The peak flux density was found to be 18.7 Jy and 13.487 Jy in 2005, and 2019, respectively. The distribution of the water maser components is up to 165 mas, $\sim$350 au at 2.16 kpc. We find that the Eastern methanol masers most probably trace outflows from the region of millimetre source mm1. The water masers to the West lie in a disc (flared or interacting with outflow/infall) around another more evolved millimetre source (13-s). The maser distribution suggests that the disc lies at an angle of 60$^{\circ}$ or more to the plane of the sky and the observed line of sight velocities then suggest an enclosed mass between 44 M$_{\odot}$ and as little as 11 M$_{\odot}$ if the disc is edge-on. The Western methanol masers may be infalling.
Jinjin Xie, Jingwen Wu, Gary A. Fuller, Nicolas Peretto, Zhiyuan Ren, Longfei Chen, Yaoting Yan, Guodong Li, Yan Duan, Jifeng Xia, Yongxiong Wang, Di Li We investigate the infall properties in a sample of 11 infrared dark clouds (IRDCs) showing blue-asymmetry signatures in HCO$^{+}$ J=1--0 line profiles. We used JCMT to conduct mapping observations in HCO$^{+}$ J=4--3 as well as single-pointing observations in HCO$^{+}$ J =3--2, towards 23 clumps in these IRDCs. We applied the HILL model to fit these observations and derived infall velocities in the range of 0.5-2.7 km s$^{-1}$, with a median value of 1.0 km s$^{-1}$, and obtained mass accretion rates of 0.5-14$\times$10$^{-3}$ Msun yr$^{-1}$. These values are comparable to those found in massive star forming clumps in later evolutionary stages. These IRDC clumps are more likely to form star clusters. HCO$^{+}$ J =3--2 and HCO$^{+}$ J =1--0 were shown to trace infall signatures well in these IRDCs with comparable inferred properties. HCO$^{+}$ J=4--3, on the other hand, exhibits infall signatures only in a few very massive clumps, due to smaller opacties. No obvious correlation for these clumps was found between infall velocity and the NH3/CCS ratio.
Deuterium represents the only bound isotope in the universe with atomic mass number $A=2$. Motivated by the possibility of other universes, where the strong force could be stronger, this paper considers the effects of bound diprotons and dineutrons on stars. We find that the existence of additional stable nuclei with $A=2$ has relatively modest effects on the universe. Previous work indicates that Big Bang Nucleosynthesis (BBN) produces more deuterium, but does not lead to catastrophic heavy element production. This paper revisits BBN considerations and confirms that the universe is left with an ample supply of hydrogen and other light nuclei for typical cosmological parameters. Using the $MESA$ numerical package, we carry out stellar evolution calculations for universes with stable diprotons, with nuclear cross sections enhanced by large factors $X$. This work focuses on $X=10^{15}-10^{18}$, but explores the wider range $X$ = $10^{-3}-10^{18}$. For a given stellar mass, the presence of stable diprotons leads to somewhat brighter stars, with the radii and photospheric temperatures roughly comparable to thoese of red giants. The central temperature decreases from the characteristic value of $T_c\approx1.5\times10^7$ K for hydrogen burning down to the value of $T_c\approx10^6$ K characteristic of deuterium burning. The stellar lifetimes are smaller for a given mass, but with the extended possible mass range, the smallest stars live for trillions of years, far longer than the current cosmic age. Finally, the enhanced cross sections allow for small, partially degenerate objects with mass $M_\ast=1-10M_J$ to produce significant steady-state luminosity and thereby function as stars.
Jinjin Xie, Gary A. Fuller, Di Li, Longfei Chen, Zhiyuan Ren, Jingwen Wu, Yan Duan, Junzhi Wang, Juan Li, Nicolas Peretto, Tie Liu, Zhiqiang Shen We present one of the first Shanghai Tian Ma Radio Telescope (TMRT) K Band observations towards a sample of 26 infrared dark clouds (IRDCs). We observed the (1,1), (2,2), (3,3), and (4,4) transitions of NH$_{3}$ together with CCS (2$_{1}$-1$_{0}$) and HC$_{3}$N $J\,$=2-1, simultaneously. The survey dramatically increases the existing CCS-detected IRDC sample from 8 to 23, enabling a better statistical study of the ratios of carbon-chain molecules (CCM) to N-bearing molecules in IRDCs. With the newly developed hyperfine group ratio (HFGR) method of fitting NH$_{3}$ inversion lines, we found the gas temperature to be between 10 and 18 K. The column density ratios of CCS to NH$_{3}$ for most of the IRDCs are less than 10$^{-2}$, distinguishing IRDCs from low-mass star-forming regions. We carried out chemical evolution simulations based on a three-phase chemical model NAUTILUS. Our measurements of the column density ratios between CCM and NH$_{3}$ are consistent with chemical evolutionary ages of $\lesssim$10$^{5}$ yr in the models. Comparisons of the data and chemical models suggest that CCS, HC$_{3}$N, and NH$_{3}$ are sensitive to the chemical evolutionary stages of the sources.
Chakali Eswaraiah, Di Li, Ray S. Furuya, Tetsuo Hasegawa, Derek Ward-Thompson, Keping Qiu, Nagayoshi Ohashi, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, David Berry, Yasuo Doi, Tao-Chung Ching, Shih-Ping Lai, Jia-Wei Wang, Patrick M. Koch, Jungmi Kwon, Woojin Kwon, Pierre Bastien, Doris Arzoumanian, et al (125) We have obtained sensitive dust continuum polarization observations at 850 $\mu$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $\mu$G, 44$\pm$16 $\mu$G, and 12$\pm$5 $\mu$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux.
N. Falstad, S. Aalto, S. König, K. Onishi, S. Muller, M. Gorski, M. Sato, F. Stanley, F. Combes, E. González-Alfonso, J. G. Mangum, A. S. Evans, L. Barcos-Muñoz, G. C. Privon, S. T. Linden, T. Díaz-Santos, S. Martín, K. Sakamoto, N. Harada, G. A. Fuller, et al (11) Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact and dusty nuclei. The intense infrared radiation arising from warm dust in these sources is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib), with the brightest emission found in compact obscured nuclei (CONs). We aim to establish how common CONs are in the local Universe, and whether their prevalence depends on the luminosity or other properties of the host galaxy. We have conducted an Atacama Large Millimeter/submillimeter Array (ALMA) survey of the rotational J=3-2 transition of HCN-vib in a sample of 46 far-infrared luminous galaxies. Compact obscured nuclei are identified in 38 percent of ULIRGs, 21 percent of LIRGs, and 0 percent of lower luminosity galaxies. We find no dependence on the inclination of the host galaxy, but strong evidence of lower IRAS 25 to 60 \mum flux density ratios (f25/f60) in CONs compared to the rest of the sample. Furthermore, we find that CONs have stronger silicate features (s9.7\mum) but similar PAH equivalent widths (EQW6.2\mum) compared to other galaxies. In the local Universe, CONs are primarily found in (U)LIRGs. High resolution continuum observations of the individual nuclei are required to determine if the CON phenomenon is related to the inclinations of the nuclear disks. The lower f25/f60 ratios in CONs as well as the results for the mid-infrared diagnostics investigated are consistent with large dust columns shifting the nuclear radiation to longer wavelengths, making the mid- and far-infrared "photospheres" significantly cooler than the interior regions. To assess the importance of CONs in the context of galaxy evolution, it is necessary to extend this study to higher redshifts where (U)LIRGs are more common.
C. Gieser, H. Beuther, D. Semenov, A. Ahmadi, S. Suri, T. Möller, M.T. Beltran, P. Klaassen, Q. Zhang, J.S. Urquhart, Th. Henning, S. Feng, R. Galván-Madrid, V. de Souza Magalhães, L. Moscadelli, S. Longmore, S. Leurini, R. Kuiper, T. Peters, K.M. Menten, et al (13) We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again.
Collapsing supermassive stars ($M \gtrsim 3 \times 10^4 M_{\odot}$) at high redshifts can naturally provide seeds and explain the origin of the supermassive black holes observed in the centers of nearly all galaxies. During the collapse of supermassive stars, a burst of non-thermal neutrinos is generated with a luminosity that could greatly exceed that of a conventional core collapse supernova explosion. In this work, we investigate the extent to which the neutrinos produced in these explosions can be observed via coherent elastic neutrino-nucleus scattering (CE$\nu$NS). Large scale direct dark matter detection experiments provide particularly favorable targets. We find that upcoming $\mathcal{O}(100)$ tonne-scale experiments will be sensitive to the collapse of individual supermassive stars at distances as large as $\mathcal{O}(10)$ Mpc.
This study investigate the effectiveness of using Deep Learning (DL) for the classification of planetary nebulae (PNe). It focusses on distinguishing PNe from other types of objects, as well as their morphological classification. We adopted the deep transfer learning approach using three ImageNet pre-trained algorithms. This study was conducted using images from the Hong Kong/Australian Astronomical Observatory/Strasbourg Observatory H-alpha Planetary Nebula research platform database (HASH DB) and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). We found that the algorithm has high success in distinguishing True PNe from other types of objects even without any parameter tuning. The Matthews correlation coefficient is 0.9. Our analysis shows that DenseNet201 is the most effective DL algorithm. For the morphological classification, we found for three classes, Bipolar, Elliptical and Round, half of objects are correctly classified. Further improvement may require more data and/or training. We discuss the trade-offs and potential avenues for future work and conclude that deep transfer learning can be utilized to classify wide-field astronomical images.
We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $\mu$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical with a mean mass per unit length, $M/L$, of ~80 M$_\odot$ pc$^{-1}$, after background subtraction, with local increases reaching values of ~150 M$_\odot$ pc$^{-1}$. This would require a magnetic field strength of ~500 $\mu$G to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500-7000 au (~0.03 pc) along the filaments. This equals ~2.0 times the radius of the flat inner region. We find an anti-correlation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints of the distribution of the dense gas within super-critical filaments showing pre- and protostellar activity.
D. Arzoumanian, R. Furuya, T. Hasegawa, M. Tahani, S. Sadavoy, C. L. H. Hull, D. Johnstone, P. M. Koch, S.-i. Inutsuka, Y. Doi, T. Hoang, T. Onaka, K. Iwasaki, Y. Shimajiri, T. Inoue, N. Peretto, P. André, P. Bastien, D. Berry, H.-R. V. Chen, et al (128) [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$\mu$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($\theta_{B}$). We derive the power spectrum of the intensity and $\theta_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $\theta_{B}$ varies coherently along the filaments. The observed power spectrum of $\theta_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $\theta_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $\theta_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars.
The $pp$-chain of nuclear reactions is the primary route for energy production in the Sun. The first step in that reaction sequence converts two protons to a deuterium nucleus with the emission of a positron and electron neutrino. This reaction is extremely slow because it is a weak interaction, and significantly, it involves quantum tunneling through the Coulomb barrier. Though the reaction rate can be calculated with high confidence in the Standard Model, it has not been measured at solar energies. If there exist interactions that are engendered by non-standard mediators then the rate of this reaction in the Sun could be altered. We probe such non-standard interactions by comparing calculations of solar evolution to the current solar system age in the presence and absence of the non-standard mediators. These reveal ranges of non-standard mediator mass and couplings that are inconsistent with measured properties of the Sun, including solar neutrino results. Our constraints on these non-standard parameters, in many cases overlapping those derived via other considerations, could be extended further with better confidence in the value of the metalicity of the Sun and the solar neutrino CNO flux. Intriguingly, our work reveals a degeneracy between the solar metalicity and the presence of the invoked non-standard mediators.
The IRDC SDC335.579-0.292 (SDC335) is a massive star-forming cloud found to be globally collapsing towards one of the most massive star forming cores in the Galaxy. SDC335 hosts three high-mass protostellar objects at early stages of their evolution and archival ALMA Cycle 0 data indicate the presence of at least one molecular outflow in the region. Observations of molecular outflows from massive protostellar objects allow us to estimate the accretion rates of the protostars as well as to assess the disruptive impact that stars have on their natal clouds. The aim of this work is to identify and analyse the properties of the protostellar-driven molecular outflows within SDC335 and use these outflows to help refine the properties of the protostars. We imaged the molecular outflows in SDC335 using new data from the ATCA of SiO and Class I CH$_3$OH maser emission (~3 arcsec) alongside observations of four CO transitions made with APEX and archival ALMA CO, $^{13}$CO (~1 arcsec), and HNC data. We introduced a generalised argument to constrain outflow inclination angles based on observed outflow properties. We used the properties of each outflow to infer the accretion rates on the protostellar sources driving them and to deduce the evolutionary characteristics of the sources. We identify three molecular outflows in SDC335, one associated with each of the known compact HII regions. The outflow properties show that the SDC335 protostars are in the early stages (Class 0) of their evolution, with the potential to form stars in excess of 50 M$_{\odot}$. The measured total accretion rate onto the protostars is $1.4(\pm 0.1) \times 10^{-3}$M$_{\odot}$ yr$^{-1}$, comparable to the total mass infall rate toward the cloud centre on parsec scales of 2.5$(\pm 1.0) \times 10^{-3}$M$_{\odot}$ yr$^{-1}$, suggesting a near-continuous flow of material from cloud to core scales. [abridged].
Nguyen Bich Ngoc, Pham Ngoc Diep, Harriet Parsons, Kate Pattle, Thiem Hoang, Derek Ward-Thompson, Le Ngoc Tram, Charles L. H. Hull, Mehrnoosh Tahani, Ray Furuya, Pierre Bastien, Keping Qiu, Tetsuo Hasegawa, Woojin Kwon, Yasuo Doi, Shih-Ping Lai, Simon Coude, David Berry, Tao-Chung Ching, Jihye Hwang, et al (129) We report the first high spatial resolution measurement of magnetic fields surrounding LkH$\alpha$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $\mu$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrm{\mu}$G. The mass-to-magnetic-flux ratio in critical value units, $\lambda\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$\alpha$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $\delta B/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $\alpha=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$\alpha$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend.
Hsi-Wei Yen, Patrick M. Koch, Charles L. H. Hull, Derek Ward-Thompson, Pierre Bastien, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Tao-Chung Ching, Eun Jung Chung, Simon Coude, James Di Francesco, Pham Ngoc Diep, Yasuo Doi, Chakali Eswaraiah, Sam Falle, Gary Fuller, Ray S. Furuya, Ilseung Han, et al (26) We compare the directions of molecular outflows of 62 low-mass Class 0 and I protostars in nearby (<450 pc) star-forming regions with the mean orientations of the magnetic fields on 0.05-0.5 pc scales in the dense cores/clumps where they are embedded. The magnetic field orientations were measured using the JCMT POL-2 data taken by the BISTRO-1 survey and from the archive. The outflow directions were observed with interferometers in the literature. The observed distribution of the angles between the outflows and the magnetic fields peaks between 15 and 35 degrees. After considering projection effects, our results could suggest that the outflows tend to be misaligned with the magnetic fields by 50+/-15 degrees in three-dimensional space and are less likely (but not ruled out) randomly oriented with respect to the magnetic fields. There is no correlation between the misalignment and the bolometric temperatures in our sample. In several sources, the small-scale (1000-3000 au) magnetic fields is more misaligned with the outflows than their large-scale magnetic fields, suggesting that the small-scale magnetic field has been twisted by the dynamics. In comparison with turbulent MHD simulations of core formation, our observational results are more consistent with models in which the energy densities in the magnetic field and the turbulence of the gas are comparable. Our results also suggest that the misalignment alone cannot sufficiently reduce the efficiency of magnetic braking to enable formation of the observed number of large Keplerian disks with sizes larger than 30-50 au.
The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for estimation of model parameters and neutrino flavor evolution histories. We combine data assimilation, ordinary differential equation solvers, and neural networks to craft an inference approach tailored for non-linear dynamical systems. Using this architecture, and a simple two-neutrino, two-flavor model, we test various optimization algorithms with the help of four experimental setups. We find that employing this new architecture, together with evolutionary optimization algorithms, accurately captures flavor histories in the four experiments. This work provides more options for extending inference techniques to large numbers of neutrinos.