Showing 1–50 of 132 results for author: Whitworth, A P

ALMA-IMF XIV: Free-Free Templates Derived from H$41α$ and Ionized Gas Content in Fifteen Massive Protoclusters

Authors: Roberto Galván-Madrid, Daniel J. Díaz-González, Frédérique Motte, Adam Ginsburg, Nichol Cunningham, Karl M. Menten, Mélanie Armante, Mélisse Bonfand, Jonathan Braine, Timea Csengeri, Pierre Dell'Ova, Fabien Louvet, Thomas Nony, Rudy Rivera-Soto, Patricio Sanhueza, Amelia M. Stutz, Friedrich Wyrowski, Rodrigo H. Álvarez-Gutiérrez, Tapas Baug, Sylvain Bontemps, Leonardo Bronfman, Manuel Fernández-López, Antoine Gusdorf, Atanu Koley, Hong-Li Liu , et al. (3 additional authors not shown)

Abstract: We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen i… ▽ More We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41α$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41α$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $ΔV_\mathrm{dyn}$ with $D$ and $Q_0$. $ΔV_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

ALMA-IMF IX: Catalog and Physical Properties of 315 SiO Outflow Candidates in 15 Massive Protoclusters

Authors: A. P. M. Towner, A. Ginsburg, P. Dell'Ova, A. Gusdorf, S. Bontemps, T. Csengeri, R. Galván-Madrid, F. K. Louvet, F. Motte, P. Sanhueza, A. M. Stutz, J. Bally, T. Baug, H. R. V. Chen, N. Cunningham, M. Fernández-López, H. -L. Liu, X. Lu, T. Nony, M. Valeille-Manet, B. Wu, R. H. Álvarez-Gutiérrez, M. Bonfand, J. Di Francesco, Q. Nguyen-Luong , et al. (2 additional authors not shown)

Abstract: We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,00… ▽ More We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,000 years, yielding median mass, momentum, and energy rates of $\dot{M}$ = 10$^{-4.4}$ M$_{\odot}$ yr$^{-1}$, $\dot{P}$ = 10$^{-3.2}$ M$_{\odot}$ km/s yr$^{-1}$, and $\dot{E}$ = 1 L$_{\odot}$. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (~3$-$5$σ$), and no correlations above 3$σ$ with clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass - which has been previously described in the literature - may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15-30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force-bolometric luminosity relationship in the literature, and our data extend this relationship up to L $\geq$ 10$^6$ L$_{\odot}$ and $\dot{P}$ $\geq$ 1 M$_{\odot}$ km/s yr$^{-1}$. Our lack of correlation with clump L/M is inconsistent with models of protocluster formation in which all protostars start forming at the same time. △ Less

Submitted 27 October, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 46 pages, 14 figures, 10 tables. This publication has an associated Zenodo entry, which can be found here: https://zenodo.org/records/8350595

Line emission from filaments in molecular clouds

Authors: F. D. Priestley, D. Arzoumanian, A. P. Whitworth

Abstract: Filamentary structures are often identified in column density maps of molecular clouds, and appear to be important for both low- and high-mass star formation. Theoretically, these structures are expected to form in regions where the supersonic cloud-scale turbulent velocity field converges. While this model of filament formation successfully reproduces several of their properties derived from colu… ▽ More Filamentary structures are often identified in column density maps of molecular clouds, and appear to be important for both low- and high-mass star formation. Theoretically, these structures are expected to form in regions where the supersonic cloud-scale turbulent velocity field converges. While this model of filament formation successfully reproduces several of their properties derived from column densities, it is unclear whether it can also reproduce their kinematic features. We use a combination of hydrodynamical, chemical and radiative transfer modelling to predict the emission properties of these dynamically-forming filaments in the $^{13}$CO, HCN and N$_2$H$^+$ $J=1-0$ rotational lines. The results are largely in agreement with observations; in particular, line widths are typically subsonic to transonic, even for filaments which have formed from highly supersonic inflows. If the observed filaments are formed dynamically, as our results suggest, no equilibrium analysis is possible, and simulations which presuppose the existence of a filament are likely to produce unrealistic results. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 9 pages, 9 figures. MNRAS accepted

ALMA-IMF. V. Prestellar and protostellar core populations in the W43 cloud complex

Authors: T. Nony, R. Galvan-Madrid, F. Motte, Y. Pouteau, N. Cunningham, F. Louvet, A. M. Stutz, B. Lefloch, S. Bontemps, N. Brouillet, A. Ginsburg, I. Joncour, F. Herpin, P. Sanhueza, T. Csengeri, A. P. M. Towner, M. Bonfand, M. Fernández-López, T. Baug, L. Bronfman, G. Busquet, J. Di Francesco, A. Gusdorf, X. Lu, F. Olguin , et al. (2 additional authors not shown)

Abstract: The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and… ▽ More The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $α=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $α=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores. △ Less

Submitted 15 March, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

Journal ref: A&A 674, A75 (2023)

Do simulated molecular clouds look like real ones?

Authors: F. D. Priestley, P. C. Clark, A. P Whitworth

Abstract: Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work… ▽ More Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work has discovered a tight, often close-to-linear relationship between the integrated intensity of molecular lines and the total column density of the cloud material. We combine magnetohydrodynamical simulations, time-dependent chemistry, and radiative transfer to produce synthetic molecular line observations of model clouds. We find similarly tight correlations between line intensity and column density to those observed, although the linear behaviour is only seen in isolated (as opposed to colliding) model clouds. This linear relationship is not due to optically thin emission; all lines investigated have high optical depths, and the increase in integrated intensity with column density is due to higher velocity dispersion along the line of sight. Overall, the idealised models commonly used in the literature appear to be reasonably accurate representations of real molecular clouds. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 10 pages, 10 figures. MNRAS accepted. Data publicly available at http://cloudzoo.astro.cf.ac.uk/Downloads/202211/

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: 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, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

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

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Differences in chemical evolution between isolated and embedded prestellar cores

Authors: F. D. Priestley, A. P. Whitworth, E. Fogerty

Abstract: Models of prestellar cores often assume that the cores are isolated from their environment - material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. We investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor-Ebert spheres,… ▽ More Models of prestellar cores often assume that the cores are isolated from their environment - material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. We investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor-Ebert spheres, and show that the density of the ambient medium has a large impact on the resulting chemical properties of the cores. Models embedded in high-density, low-temperature surroundings have greatly enhanced abundances of several molecules, such as CO and CS, compared to models with more diffuse surroundings, corresponding to relatively isolated cores. The predicted intensities and profile shapes of molecular lines are also affected. The density of the ambient medium has a stronger effect on the chemical evolution than whether the cores are initially in or out of equilibrium. This suggests that the impact of environment cannot be neglected when modelling chemistry in prestellar cores; the results of these models are highly sensitive to the assumptions made about the core surroundings. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 7 pages, 5 figures. MNRAS accepted

RJ-plots: An improved method to classify structures objectively

Authors: Seamus D. Clarke, Sarah E. Jaffa, Anthony P. Whitworth

Abstract: The interstellar medium is highly structured, presenting a range of morphologies across spatial scales. The large data sets resulting from observational surveys and state-of-the-art simulations studying these hierarchical structures means that identification and classification must be done in an automated fashion to be efficient. Here we present RJ-plots, an improved version of the automated morph… ▽ More The interstellar medium is highly structured, presenting a range of morphologies across spatial scales. The large data sets resulting from observational surveys and state-of-the-art simulations studying these hierarchical structures means that identification and classification must be done in an automated fashion to be efficient. Here we present RJ-plots, an improved version of the automated morphological classification technique J-plots developed by Jaffa et al. (2018). This method allows clear distinctions between quasi-circular/elongated structures and centrally over/under-dense structures. We use the recent morphological SEDIGISM catalogue of Neralwar et al. (2022) to show the improvement in classification resulting from RJ-plots, especially for ring-like and concentrated cloud types. We also find a strong correlation between the central concentration of a structure and its star formation efficiency and dense gas fraction, as well as a lack of correlation with elongation. Furthermore, we use the accreting filament simulations of Clarke et al. (2020) to highlight a multi-scale application of RJ-plots, finding that while spherical structures become more common at smaller scales they are never the dominant structure down to $r\sim0.03$ pc. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 11 pages, 9 figures, accepted by MNRAS

The widths of magnetised filaments in molecular clouds

Authors: F. D. Priestley, A. P. Whitworth

Abstract: Filaments are an ubiquitous feature of molecular clouds, and appear to play a critical role in assembling the material to form stars. The dominant filaments are observed to have a rather narrow range of widths around $\sim 0.1$ pc, and to be preferentially aligned perpendicularly to the direction of the local magnetic field. We have previously argued that the observed filament widths can be explai… ▽ More Filaments are an ubiquitous feature of molecular clouds, and appear to play a critical role in assembling the material to form stars. The dominant filaments are observed to have a rather narrow range of widths around $\sim 0.1$ pc, and to be preferentially aligned perpendicularly to the direction of the local magnetic field. We have previously argued that the observed filament widths can be explained if filaments are formed by converging, mildly supersonic flows, resulting from large-scale turbulent motions in the parent molecular cloud. Here we demonstrate that the introduction of a magnetic field perpendicular to the filament long axis does not greatly alter this conclusion, as long as the mass-to-flux ratio is supercritical. The distribution of widths for supercritical magnetised filaments formed via this mechanism is peaked at slightly higher values, and is slightly broader, than for non-magnetised filaments, but still reproduces the basic properties of the width distributions derived from far-infrared observations of molecular clouds. In contrast, subcritical filaments have width distributions with a fundamentally different shape, and typically have much larger widths than those observed. Both subcritical and supercritical filaments are consistent with the observed lack of correlation between filament widths and filament surface densities. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 9 pages, 8 figures. MNRAS accepted

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software p… ▽ More The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

Journal ref: A&A 664, A26 (2022)

Protostellar Outflows: a window to the past

Authors: P. F. Rohde, S. Walch, D. Seifried, A. P. Whitworth, S. D. Clarke

Abstract: During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink… ▽ More During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink particles, and the sink particles launch episodic outflows using a subgrid model. The Optics algorithm is used to identify individual episodic bullets within the outflows. The parameters of the overall outflow and the individual bullets are then used to estimate the age and energetics of the outflow, and the accretion events that triggered it; and to evaluate how reliable these estimates are, if observational uncertainties and selection effects (like inclination) are neglected. Of the commonly used methods for estimating outflow ages, it appears that those based on the length and speed of advance of the lobe are the most reliable in the early phases of evolution, and those based on the width of the outflow cavity and the speed of advance are most reliable during the later phases. We describe a new method that is almost as accurate as these methods, and reliable throughout the evolution. In addition we show how the accretion history of the protostar can be accurately reconstructed from the dynamics of the bullets if each lobe contains at least two bullets. The outflows entrain about ten times more mass than originally ejected by the protostar. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: 20 pages, 19 figures, published in MNRAS

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

Authors: F. Motte, S. Bontemps, T. Csengeri, Y. Pouteau, F. Louvet, A. M. Stutz, N. Cunningham, A. López-Sepulcre, N. Brouillet, R. Galván-Madrid, A. Ginsburg, L. Maud, A. Men'shchikov, F. Nakamura, T. Nony, P. Sanhueza, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, G. Busquet, E. Chapillon, D. Díaz-González, M. Fernández-López, A. E. Guzmán , et al. (39 additional authors not shown)

Abstract: The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2… ▽ More The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 27 pages, 9 figures, accepted for publication by A&A

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

The JCMT BISTRO Survey: Evidence for Pinched Magnetic Fields in Quiescent Filaments of NGC 1333

Authors: Yasuo Doi, Kohji Tomisaka, Tetsuo Hasegawa, Simon Coudé, Doris Arzoumanian, Pierre Bastien, Masafumi Matsumura, Mehrnoosh Tahani, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, James Di Francesco, Yoshito Shimajiri, Ray S. Furuya, Jungmi Kwon, Motohide Tamura, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Thiem Hoang, Florian Kirchschlager, Jihye Hwang, Chakali Eswaraiah, Patrick M. Koch, Anthony P. Whitworth, Kate Pattle , et al. (11 additional authors not shown)

Abstract: We investigate the internal 3D magnetic structure of dense interstellar filaments within NGC 1333 using polarization data at $850 μ\mathrm{m}$ from the $B$-fields In STar-forming Region Observations survey at the James Clerk Maxwell Telescope. Theoretical models predict that the magnetic field lines in a filament will tend to be dragged radially inward (i.e., pinched) toward the central axis due t… ▽ More We investigate the internal 3D magnetic structure of dense interstellar filaments within NGC 1333 using polarization data at $850 μ\mathrm{m}$ from the $B$-fields In STar-forming Region Observations survey at the James Clerk Maxwell Telescope. Theoretical models predict that the magnetic field lines in a filament will tend to be dragged radially inward (i.e., pinched) toward the central axis due to the filament's self-gravity. We study the cross-sectional profiles of the total intensity ($I$) and polarized intensity (PI) of dust emission in four segments of filaments unaffected by local star formation that are expected to retain a pristine magnetic field structure. We find that the filaments' FWHM in PI are not the same as those in $I$, with two segments being appreciably narrower in PI (FWHM ratio $\simeq 0.7-0.8$) and one segment being wider (FWHM ratio $\simeq 1.3$). The filament profiles of the polarization fraction ($P$) do not show a minimum at the spine of the filament, which is not in line with an anticorrelation between $P$ and $I$ normally seen in molecular clouds and protostellar cores. Dust grain alignment variation with density cannot reproduce the observed $P$ distribution. We demonstrate numerically that the $I$ and PI cross-sectional profiles of filaments in magnetohydrostatic equilibrium will have differing relative widths depending on the viewing angle. The observed variations of FWHM ratios in NGC 1333 are therefore consistent with models of pinched magnetic field structures inside filaments, and especially if they are magnetically near-critical or supercritical. △ Less

Submitted 8 December, 2021; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJL

Investigating variations in the dust emissivity index in the Andromeda galaxy

Authors: G. Athikkat-Eknath, S. A. Eales, M. W. L. Smith, A. Schruba, K. A. Marsh, A. P. Whitworth

Abstract: Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index ($β$). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use $^{12}$CO(J=1-0) observations taken by the Combined Array for Research in… ▽ More Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index ($β$). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use $^{12}$CO(J=1-0) observations taken by the Combined Array for Research in Millimeter Astronomy (CARMA) and dust maps derived from \textit{Herschel} images, both with an angular resolution of 8" and spatial resolution of 30 pc, to study variations in $β$ across an area of $\approx$ 18.6 kpc$^2$ in M31. We extract sources, which we identify as molecular clouds, by applying the astrodendro algorithm to the $^{12}$CO and dust maps, which as a byproduct allows us to compare continuum emission from dust and CO emission as alternative ways of finding molecular clouds. We then use these catalogues to investigate whether there is evidence that $β$ is different inside and outside molecular clouds. Our results confirm the radial variations of $β$ seen in previous studies. However, we find little difference between the average $β$ inside molecular clouds compared to outside molecular clouds, in disagreement with models which predict an increase of $β$ in dense environments. Finally, we find some clouds traced by dust with very little CO which may be either clouds dominated by atomic gas or clouds of molecular gas that contain little CO. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

The origin of a universal filament width in molecular clouds

Authors: F. D. Priestley, A. P. Whitworth

Abstract: Filamentary structures identified in far-infrared observations of molecular clouds are typically found to have full-widths at half-maximum $\sim\!0.1$ pc. However, the physical explanation for this phenomenon is currently uncertain. We use hydrodynamic simulations of cylindrically-symmetric converging flows to show that the full-width at half-maximum of the resulting filament's surface density pro… ▽ More Filamentary structures identified in far-infrared observations of molecular clouds are typically found to have full-widths at half-maximum $\sim\!0.1$ pc. However, the physical explanation for this phenomenon is currently uncertain. We use hydrodynamic simulations of cylindrically-symmetric converging flows to show that the full-width at half-maximum of the resulting filament's surface density profile, FWHM$_Σ$, is closely related to the location of the accretion shock, where the inflow meets the boundary of the filament. For inflow Mach Number, ${\cal M}$, between 1 and 5, filament FWHM$_Σ$s fall in the range $0.03$ pc $\lesssim$ FWHM$_Σ \lesssim 0.3$ pc, with higher ${\cal M}$ resulting in narrower filaments. A large sample of filaments, seen at different evolutionary stages and with different values of ${\cal M}$, naturally results in a peaked distribution of FWHM$_Σ$s similar in shape to that obtained from far-infrared observations of molecular clouds. However, unless the converging flows are limited to ${\cal M} \lesssim 3$, the peak of the distribution of FWHM$_Σ$s is below the observed $\sim 0.1$ pc. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 11 pages, 12 figures. MNRAS accepted

Molecular line signatures of cloud-cloud collisions

Authors: F. D. Priestley, A. P. Whitworth

Abstract: Collisions between interstellar gas clouds are potentially an important mechanism for triggering star formation. This is because they are able to rapidly generate large masses of dense gas. Observationally, cloud collisions are often identified in position-velocity (PV) space through bridging features between intensity peaks, usually of CO emission. Using a combination of hydrodynamical simulation… ▽ More Collisions between interstellar gas clouds are potentially an important mechanism for triggering star formation. This is because they are able to rapidly generate large masses of dense gas. Observationally, cloud collisions are often identified in position-velocity (PV) space through bridging features between intensity peaks, usually of CO emission. Using a combination of hydrodynamical simulations, time-dependent chemistry, and radiative transfer, we produce synthetic molecular line observations of overlapping clouds that are genuinely colliding, and overlapping clouds that are just chance superpositions. Molecules tracing denser material than CO, such as NH$_3$ and HCN, reach peak intensity ratios of $0.5$ and $0.2$ with respect to CO in the `bridging feature' region of PV space for genuinely colliding clouds. For overlapping clouds that are just chance superpositions, the peak NH$_3$ and HCN intensities are co-located with the CO intensity peaks. This represents a way of confirming cloud collisions observationally, and distinguishing them from chance alignments of unrelated material. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 7 pages, 4 figures, MNRAS accepted

Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay

Authors: Richard Wünsch, Stefanie Walch, František Dinnbier, Daniel Seifried, Sebastian Haid, Andre Klepitko, Anthony P. Whitworth, Jan Palouš

Abstract: The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium. In this paper we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code FLASH, as… ▽ More The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium. In this paper we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code FLASH, as a module of the tree solver developed by Wünsch et al. However, the method itself is independent of the host code and can be implemented in any grid based or particle based hydrodynamics code. A key advantage of TreeRay is that its computational cost is independent of the number of sources, making it suitable for simulations with many point sources (e.g. massive star clusters) as well as simulations where diffuse emission is important. A very efficient communication and tree-walk strategy enables TreeRay to achieve almost ideal parallel scalings. TreeRay can easily be extended with sub-modules to treat radiative transfer at different wavelengths and to implement related physical processes. Here, we focus on ionising (EUV) radiation and use the On-the-Spot approximation to test the method and its parameters. The ability to set the tree solver time step independently enables the speedy calculation of radiative transfer in a multi-phase interstellar medium, where the hydrodynamic time step is typically limited by the sound speed of the hot gas produced in stellar wind bubbles or supernova remnants. We show that complicated simulations of star clusters with feedback from multiple massive stars become feasible with TreeRay. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 26 pages, 21 figures, accepted by MNRAS

Feedback from an O star formed in a filament

Authors: Anthony P Whitworth, Felix D Priestley

Abstract: We explore a simple semi-analytic model for what happens when an O star (or cluster of O stars) forms in an isolated filamentary cloud. The model is characterised by three configuration parameters: the radius of the filament, R_FIL, the mean density of H_2 in the filament, n_FIL, and the rate at which the O star emits ionising photons, Ndot_LyC. We show that for a wide range of these configuration… ▽ More We explore a simple semi-analytic model for what happens when an O star (or cluster of O stars) forms in an isolated filamentary cloud. The model is characterised by three configuration parameters: the radius of the filament, R_FIL, the mean density of H_2 in the filament, n_FIL, and the rate at which the O star emits ionising photons, Ndot_LyC. We show that for a wide range of these configuration parameters, ionising radiation from the O star rapidly erodes the filament, and the ionised gas from the filament disperses into the surroundings. Under these circumstances the distance from the O star to the ionisation front (IF) is given approximately by L ~ 5.2 pc [R_FIL/0.2pc]^-1/6 [n_FIL/10^4cm^-3]^-1/3 [Ndot_LyC/10^49s^-1]^1/6 [t/Myr]^2/3, and we derive similar simple power-law expressions for other quantities, for example the rate at which ionised gas boils off the filament, and the mass of the shock-compressed layer (SCL) that is swept up behind the IF. We show that a very small fraction of the ionising radiation is expended locally, and a rather small amount of molecular gas is ionised and dispersed. We discuss some features of more realistic models, and the extent to which they might modify or invalidate the predictions of this idealised model. In particular we show that, for very large R_FIL and/or large n_FIL and/or low Ndot_LyC, continuing accretion onto the filament might trap the ionising radiation from the O star, slowing the erosion of the filament even further. △ Less

Submitted 30 April, 2021; originally announced April 2021.

Comments: 10 pages, 3 figures, accepted by MNRAS

A PPMAP analysis of the filamentary structures in Ophiuchus L1688 and L1689

Authors: A. D. P. Howard, A. P. Whitworth, M. J. Griffin, K. A. Marsh, M. W. L. Smith

Abstract: We use the PPMAP (Point Process MAPping) algorithm to re-analyse the \textit{Herschel} and SCUBA-2 observations of the L1688 and L1689 sub-regions of the Ophiuchus molecular cloud. PPMAP delivers maps with high resolution (here $14''$, corresponding to $\sim 0.01\,{\rm pc}$ at $\sim 140\,{\rm pc}$), by using the observations at their native resolutions. PPMAP also delivers more accurate dust optic… ▽ More We use the PPMAP (Point Process MAPping) algorithm to re-analyse the \textit{Herschel} and SCUBA-2 observations of the L1688 and L1689 sub-regions of the Ophiuchus molecular cloud. PPMAP delivers maps with high resolution (here $14''$, corresponding to $\sim 0.01\,{\rm pc}$ at $\sim 140\,{\rm pc}$), by using the observations at their native resolutions. PPMAP also delivers more accurate dust optical depths, by distinguishing dust of different types and at different temperatures. The filaments and prestellar cores almost all lie in regions with $N_{\rm H_2}\gtrsim 7\times 10^{21}\,{\rm cm}^{-2}$ (corresponding to $A_{_{\rm V}}\gtrsim 7$). The dust temperature, $T$, tends to be correlated with the dust opacity index, $β$, with low $T$ and low $β$ tend concentrated in the interiors of filaments. The one exception to this tendency is a section of filament in L1688 that falls -- in projection -- between the two B stars, S1 and HD147889; here $T$ and $β$ are relatively high, and there is compelling evidence that feedback from these two stars has heated and compressed the filament. Filament {\sc fwhm}s are typically in the range $0.10\,{\rm pc}$ to $0.15\,{\rm pc}$. Most filaments have line densities in the range $25\,{\rm M_{_\odot}\,pc^{-1}}$ to $65\,{\rm M_{_\odot}\,pc^{-1}}$. If their only support is thermal gas pressure, and the gas is at the canonical temperature of $10\,{\rm K}$, the filaments are highly supercritical. However, there is some evidence from ammonia observations that the gas is significantly warmer than this, and we cannot rule out the possibility of additional support from turbulence and/or magnetic fields. On the basis of their spatial distribution, we argue that most of the starless cores are likely to disperse (rather than evolving to become prestellar). △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 22 pages, 18 figures, Accepted for publication by MNRAS March 2021

Synthetic line and continuum observations of simulated turbulent clouds: the apparent widths of filaments

Authors: F. D. Priestley, A. P. Whitworth

Abstract: Filamentary structures are ubiquitous in observations of real molecular clouds, and also in simulations of turbulent, self-gravitating gas. However, making comparisons between observations and simulations is complicated by the difficulty of estimating volume-densities observationally. Here, we have post-processed hydrodynamical simulations of a turbulent isothermal molecular cloud, using a full ti… ▽ More Filamentary structures are ubiquitous in observations of real molecular clouds, and also in simulations of turbulent, self-gravitating gas. However, making comparisons between observations and simulations is complicated by the difficulty of estimating volume-densities observationally. Here, we have post-processed hydrodynamical simulations of a turbulent isothermal molecular cloud, using a full time-dependent chemical network. We have then run radiative transfer models to obtain synthetic line and continuum intensities that can be compared directly with those observed. We find that filaments have a characteristic width of $\,\sim\!0.1 \, {\rm pc}$, both on maps of their true surface density, and on maps of their $850 \, {\rm μm}$ dust-continuum emission, in agreement with previous work. On maps of line emission from CO isotopologues, the apparent widths of filaments are typically several times larger because the line intensities are poorly correlated with the surface density. On maps of line emission from dense-gas tracers such as N$_2$H$^+$ and HCN, the apparent widths of filaments are $\lesssim 0.1 \, {\rm pc}$. Thus, current observations of molecular-line emission are compatible with the universal $0.1 \, {\rm pc}$ filament width inferred from ${\it Herschel}$ observations, provided proper account is taken of abundance, optical-depth, and excitation considerations. We find evidence for $\sim 0.4 \, {\rm km \, s^{-1}}$ radial velocity differences across filaments. These radial velocity differences might be a useful indicator of the mechanism by which a filament has formed or is forming, for example the turbulent cloud scenario modelled here, as against other mechanisms such as cloud-cloud collisions. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: 11 pages, 11 figures. MNRAS accepted 06/10/20

The impact of episodic outflow feedback on stellar multiplicity and the star formation efficiency

Authors: P. F. Rohde, S. Walch, S. D. Clarke, D. Seifried, A. P. Whitworth, A. Klepitko

Abstract: The accretion of material onto young protostars is accompanied by the launching of outflows. Observations show that accretion, and therefore also outflows, are episodic. However, the effects of episodic outflow feedback on the core-scale are not well understood. We have performed 88 Smoothed Particle Hydrodynamic simulations of turbulent dense $1 \, \mathrm{M}_{\odot}$ cores, to study the influenc… ▽ More The accretion of material onto young protostars is accompanied by the launching of outflows. Observations show that accretion, and therefore also outflows, are episodic. However, the effects of episodic outflow feedback on the core-scale are not well understood. We have performed 88 Smoothed Particle Hydrodynamic simulations of turbulent dense $1 \, \mathrm{M}_{\odot}$ cores, to study the influence of episodic outflow feedback on the stellar multiplicity and the star formation efficiency (SFE). Protostars are represented by sink particles, which use a sub-grid model to capture stellar evolution, inner-disc evolution, episodic accretion and the launching of outflows. By comparing simulations with and without episodic outflow feedback, we show that simulations with outflow feedback reproduce the binary statistics of young stellar populations, including the relative proportions of singles, binaries, triples, etc. and the high incidence of twin binaries with $q\geq 0.95$; simulations without outflow feedback do not. Entrainment factors (the ratio between total outflowing mass and initially ejected mass) are typically $\sim 7\pm 2$, but can be much higher if the total mass of stars formed in a core is low and/or outflow episodes are infrequent. By decreasing both the mean mass of the stars formed and the number of stars formed, outflow feedback reduces the SFE by about a factor of 2 (as compared with simulations that do not include outflow feedback). △ Less

Submitted 24 September, 2020; originally announced September 2020.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. 19 pages, 12 figures

On the emergent System Mass Function: the contest between accretion and fragmentation

Authors: Paul C. Clark, Anthony P. Whitworth

Abstract: We propose a new model for the evolution of a star cluster's System Mass Function (SMF). The model involves both turbulent fragmentation and competitive accretion. Turbulent fragmentation creates low-mass seed proto-systems (i.e. single and multiple protostars). Some of these low-mass seed proto-systems then grow by competitive accretion to produce the high-mass power-law tail of the SMF. Turbulen… ▽ More We propose a new model for the evolution of a star cluster's System Mass Function (SMF). The model involves both turbulent fragmentation and competitive accretion. Turbulent fragmentation creates low-mass seed proto-systems (i.e. single and multiple protostars). Some of these low-mass seed proto-systems then grow by competitive accretion to produce the high-mass power-law tail of the SMF. Turbulent fragmentation is relatively inefficient, in the sense that the creation of low-mass seed proto-systems only consumes a fraction, $\sim 23\%$ (at most $\sim 50\%$), of the mass available for star formation. The remaining mass is consumed by competitive accretion. Provided the accretion rate onto a proto-system is approximately proportional to its mass ($dm/dt \propto m$), the SMF develops a power-law tail at high masses with the Salpeter slope ($\sim -2.3$). If the rate of supply of mass accelerates, the rate of proto-system formation also accelerates, as appears to be observed in many clusters. However, even if the rate of supply of mass decreases, or ceases and then resumes, the SMF evolves homologously, retaining the same overall shape, and the high-mass power-law tail simply extends to ever higher masses until the supply of gas runs out completely. The Chabrier SMF can be reproduced very accurately if the seed proto-systems have an approximately log-normal mass distribution with median mass $\sim 0.11 {\rm M}_{_\odot}$ and logarithmic standard deviation $σ_{\log_{10}(M/{\rm M}_\odot)}\sim 0.47$). △ Less

Submitted 22 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS. 13 pages, 4 figures, 1 appendix

The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333

Authors: Yasuo Doi, Tetsuo Hasegawa, Ray S. Furuya, Simon Coudé, Charles L. H. Hull, Doris Arzoumanian, Pierre Bastien, Michael Chun-Yuan Chen, James di Francesco, Rachel Friesen, Martin Houde, Shu-ichiro Inutsuka, Steve Mairs, Masafumi Matsumura, Takashi Onaka, Sarah Sadavoy, Yoshito Shimajiri, Mehrnoosh Tahani, Kohji Tomisaka, Chakali Eswaraiah, Patrick M. Koch, Kate Pattle, Chang Won Lee, Motohide Tamura, David Berry , et al. (113 additional authors not shown)

Abstract: We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary… ▽ More We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multiscale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~1 pc and remains continuous from the scales of filaments (~0.1 pc) to that of protostellar envelopes (~0.005 pc or ~1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data toward the filament network. △ Less

Submitted 21 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 33 pages, 21 figures, accepted for publication in ApJ

The apparent anti-correlation between the mass opacity of interstellar dust and the surface-density of interstellar gas

Authors: F. D. Priestley, A. P. Whitworth

Abstract: Recent analyses of ${\it Herschel}$ observations suggest that in nearby disc galaxies the dust mass opacity at $500 \, {\rm μm}$, $κ_{500}$, decreases with increasing gas surface density, $Σ_{\rm ISM}$ (Clark et al. 2019). This apparent anti-correlation between $κ_{500}$ and $Σ_{\rm ISM}$ is opposite to the behaviour expected from theoretical dust evolution models; in such models, dust in denser,… ▽ More Recent analyses of ${\it Herschel}$ observations suggest that in nearby disc galaxies the dust mass opacity at $500 \, {\rm μm}$, $κ_{500}$, decreases with increasing gas surface density, $Σ_{\rm ISM}$ (Clark et al. 2019). This apparent anti-correlation between $κ_{500}$ and $Σ_{\rm ISM}$ is opposite to the behaviour expected from theoretical dust evolution models; in such models, dust in denser, cooler regions (i.e. regions of increased $Σ_{\rm ISM}$) tends to grow and therefore to have increased $κ_{500}$. We show, using a toy model, that the presence of a range of dust temperatures along the line of sight can lead to spuriously low estimated values of $κ_{500}$. If in regions of higher $Σ_{\rm ISM}$ the range of dust temperatures extends to lower values (as seems likely), the magnitude of this effect may be sufficient to explain the apparent anti-correlation between $κ_{500}$ and $Σ_{\rm ISM}$. Therefore there may not be any need for spatial variation in the intrinsic dust properties that run counter to theoretical expectations. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 5 pages, 5 figures. MNRASL accepted 14/02/20

Characterising lognormal fractional-Brownian-motion density fields with a Convolutional Neural Network

Authors: M. L. Bates, A. P. Whitworth, O. D. Lomax

Abstract: In attempting to quantify statistically the density structure of the interstellar medium, astronomers have considered a variety of fractal models. Here we argue that, to properly characterise a fractal model, one needs to define precisely the algorithm used to generate the density field, and to specify -- at least -- three parameters: one parameter constrains the spatial structure of the field; on… ▽ More In attempting to quantify statistically the density structure of the interstellar medium, astronomers have considered a variety of fractal models. Here we argue that, to properly characterise a fractal model, one needs to define precisely the algorithm used to generate the density field, and to specify -- at least -- three parameters: one parameter constrains the spatial structure of the field; one parameter constrains the density contrast between structures on different scales; and one parameter constrains the dynamic range of spatial scales over which self-similarity is expected (either due to physical considerations, or due to the limitations of the observational or numerical technique generating the input data). A realistic fractal field must also be noisy and non-periodic. We illustrate this with the exponentiated fractional Brownian motion (xfBm) algorithm, which is popular because it delivers an approximately lognormal density field, and for which the three parameters are, respectively, the power spectrum exponent, $β$, the exponentiating factor, ${\cal S}$, and the dynamic range, ${\cal R}$. We then explore and compare two approaches that might be used to estimate these parameters: Machine Learning and the established $Δ$-Variance procedure. We show that for $2\leqβ\leq 4$ and $0\leq{\cal S}\leq 3$, a suitably trained Convolutional Neural Network is able to estimate objectively both $β$ (with root-mean-square error $ε_{_β}\sim 0.12$) and ${\cal S}$ (with $ε_{_{\cal S}}\sim 0.29$). $\;Δ$-variance is also able to estimate $β$, albeit with a somewhat larger error ($ε_{_β}\sim 0.17$) and with some human intervention, but is not able to estimate ${\cal S}$. △ Less

Submitted 13 February, 2020; originally announced February 2020.

The Hi-GAL catalogue of dusty filamentary structures in the Galactic Plane

Authors: E. Schisano, S. Molinari, D. Elia, M. Benedettini, L. Olmi, S. Pezzuto, A. Traficante, M. Brescia, S. Cavuoti, A. M. di Giorgio, S. J. Liu, T. J. T. Moore, A. Noriega-Crespo, G. Riccio, A. Baldeschi, U. Becciani, N. Peretto, M. Merello, F. Vitello, A. Zavagno, M. T. Beltrán, L. Cambrésy, D. J. Eden, G. Li Causi, M. Molinaro , et al. (5 additional authors not shown)

Abstract: The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structur… ▽ More The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 38 pages, 29 figures, 3 appendices

L1495 Revisited: A PPMAP View of a Star-Forming Filament

Authors: A. D. P. Howard, A. P. Whitworth, K. A. Marsh, S. D. Clarke, M. J. Griffin, M. W. L. Smith, O. D. Lomax

Abstract: We have analysed the Herschel and SCUBA-2 dust continuum observations of the main filament in the Taurus L1495 star forming region, using the Bayesian fitting procedure PPMAP. (i) If we construct an average profile along the whole length of the filament, it has fwhm $\simeq 0.087\pm 0.003\,{\rm pc};\;$, but the closeness to previous estimates is coincidental. (ii) If we analyse small local section… ▽ More We have analysed the Herschel and SCUBA-2 dust continuum observations of the main filament in the Taurus L1495 star forming region, using the Bayesian fitting procedure PPMAP. (i) If we construct an average profile along the whole length of the filament, it has fwhm $\simeq 0.087\pm 0.003\,{\rm pc};\;$, but the closeness to previous estimates is coincidental. (ii) If we analyse small local sections of the filament, the column-density profile approximates well to the form predicted for hydrostatic equilibrium of an isothermal cylinder. (iii) The ability of PPMAP to distinguish dust emitting at different temperatures, and thereby to discriminate between the warm outer layers of the filament and the cold inner layers near the spine, leads to a significant reduction in the surface-density, $\varSigma$, and hence in the line-density, $μ$. If we adopt the canonical value for the critical line-density at a gas-kinetic temperature of $10\,{\rm K}$, $μ_{_{\rm CRIT}}\simeq 16\,{\rm M_{_\odot}\,pc^{-1}}$, the filament is on average trans-critical, with ${\barμ}\sim μ_{_{\rm CRIT}};\;$ local sections where $μ>μ_{_{\rm CRIT}}$ tend to lie close to pre-stellar cores. (iv) The ability of PPMAP to distinguish different types of dust, i.e. dust characterised by different values of the emissivity index, $β$, reveals that the dust in the filament has a lower emissivity index, $β\leq1.5$, than the dust outside the filament, $β\geq 1.7$, implying that the physical conditions in the filament have effected a change in the properties of the dust. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 16 pages, 9 figures, Accepted for publication by MNRAS August 2019

The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

Authors: Simon Coudé, Pierre Bastien, Martin Houde, Sarah Sadavoy, Rachel Friesen, James Di Francesco, Doug Johnstone, Steve Mairs, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Derek Ward-Thompson, David Berry, Michael Chun-Yuan Chen, Jason Fiege, Erica Franzmann, Jennifer Hatchell, Kevin Lacaille, Brenda C. Matthews, Gerald H. Moriarty-Schieven, Andy Pon, Philippe André, Doris Arzoumanian, Yusuke Aso , et al. (96 additional authors not shown)

Abstract: We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We th… ▽ More We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $μ$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $μ$G and a criticality criterion $λ_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component. △ Less

Submitted 23 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal

The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $ρ$ Ophiuchus C

Authors: Junhao Liu, Keping Qiu, David Berry, James Di Francesco, Pierre Bastien, Patrick M. Koch, Ray S. Furuya, Kee-Tae Kim, Simon Coudé, Chang Won Lee, Archana Soam, Chakali Eswaraiah, Di Li, Jihye Hwang, A-Ran Lyo, Kate Pattle, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Derek Ward-Thompson, Tao-Chung Ching, Zhiwei Chen, Qilao Gu, Dalei Li, Hua-bai Li , et al. (106 additional authors not shown)

Abstract: We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic… ▽ More We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $μ$G, 136 $\pm$ 69 $μ$G, and 213 $\pm$ 115 $μ$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures

JCMT BISTRO survey: Magnetic Fields within the Hub-Filament Structure in IC 5146

Authors: Jia-Wei Wang, Shih-Ping Lai, Chakali Eswaraiah, Kate Pattle, James Di Francesco, Doug Johnstone, Patrick M. Koch, Tie Liu, Motohide Tamura, Ray S. Furuya, Takashi Onaka, Derek Ward-Thompson, Archana Soam, Kee-Tae Kim, Chang Won Lee, Chin-Fei Lee, Steve Mairs, Doris Arzoumanian, Gwanjeong Kim, Thiem Hoang, Jihye Hwang, Sheng-Yuan Liu, David Berry, Pierre Bastien, Tetsuo Hasegawa , et al. (108 additional authors not shown)

Abstract: We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scal… ▽ More We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfvénic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important. △ Less

Submitted 27 March, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 24 pages, 15 figures, accepted for publication in ApJ

RCW 120: A possible case of hit and run, elucidated by multi temperature dust mapping

Authors: K. A. Marsh, A. P. Whitworth

Abstract: We present resolution-enhanced images of warm dust at multiple temperatures and opacity index values in the star-forming bubble/HII region, RCW 120. The image set, representing a 4D hypercube of differential column density, was obtained using our Bayesian procedure, ppmap. The cool peripheral material ($\sim16$-22 K) exhibits ragged clumpy structure as noted previously by others. However, at highe… ▽ More We present resolution-enhanced images of warm dust at multiple temperatures and opacity index values in the star-forming bubble/HII region, RCW 120. The image set, representing a 4D hypercube of differential column density, was obtained using our Bayesian procedure, ppmap. The cool peripheral material ($\sim16$-22 K) exhibits ragged clumpy structure as noted previously by others. However, at higher temperatures ($\stackrel{>}{_\sim}26$ K) the geometry changes dramatically, showing a bubble boundary which is accurately circular in projection, except for the previously-reported opening in the north. Comparison with Spitzer 8 $μ$m data suggests that the $\sim26$-30 K dust seen by Herschel resides in the photodissociation region (PDR) surrounding the HII region. Its projected radial profile is consistent with that of a spherical shell, thus arguing against previous suggestions of cylindrical or planar geometry. The inferred geometry is, in fact, consistent with previous interpretations of the HII region as a classical Strömgren sphere, except for the fact that the ionising star (CD -38.11636; O8V) is displaced by more than half a radius from its geometric centre. None of the previously published models has satisfactorily accounted for that displacement. It could, however, be explained by proper motion of the O star at $\sim2$-4 km s$^{-1}$ since its formation, possibly due to a cloud-cloud collision. We suggest that the current spherical bubble constitutes the fossilised remnant of the initial expansion of the HII region following the formation of the star, which now continues to flee its formation site. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 8 pages, 7 figures. Accepted for publication in MNRAS

A simple approach to CO cooling in molecular clouds

Authors: A. P. Whitworth, S. E. Jaffa

Abstract: CO plays an important role in interstellar molecular clouds, both as a coolant, and as a diagnostic molecule. However, a proper evaluation of the cooling rate due to CO requires a determination of the populations of many levels, the spontaneous and stimulated radiative de-excitation rates between these levels, and the transfer of the emitted multi-line radiation; additionally, this must be done fo… ▽ More CO plays an important role in interstellar molecular clouds, both as a coolant, and as a diagnostic molecule. However, a proper evaluation of the cooling rate due to CO requires a determination of the populations of many levels, the spontaneous and stimulated radiative de-excitation rates between these levels, and the transfer of the emitted multi-line radiation; additionally, this must be done for three isotopologues. It would be useful to have a simple analytic formulation that avoided these complications; this could then be used in situations where CO plays an important role as a coolant, but the details of this role are not the main concern. We derive such a formulation here, by first considering the two asymptotic forms that obtain in the limits of (a) low volume-density and optical depth, and (b) high volume-density and optical depth. These forms are then combined in such a way as to fit the detailed numerical results from Goldsmith \& Langer (1978), which cover low temperatures, and a range of physical conditions where the interplay of thermal and sub-thermal excitation, optical-depth effects, and the contributions from rare isotopologues, are all important. The fit is obtained using the Metropolis-Hastings method, and reproduces the results of GL78 well. It is a purely local and analytic function of state --- specifically a function of the density, $ρ$, isothermal sound speed, $a$, CO abundance, $X_{_{\rm CO}}$, and velocity divergence, $\nabla\cdot{\boldsymbol\upsilon}$. As an application, we consider the cooling layer following a slow steady non-magnetic planar J-shock. We show that, if the post-shock cooling is dominated by CO and its isotopologues, the thickness of the post-shock cooling layer is very small and approximately independent of the pre-shock velocity, $\upsilon_o$, or pre-shock isothermal sound speed, $a_o$. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 8 pages, 1 figure,

Journal ref: Astronomy and Astrophysics, 2018, Vol. 611, A20

Evolution of Hubble wedges in episodic protostellar outflows

Authors: P. F. Rohde, S. Walch, D. Seifried, A. P. Whitworth, S. D. Clarke

Abstract: Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin… ▽ More Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin of these features, we have developed a new episodic outflow model for the SPH code {\sc gandalf}, which mimics the accretion and ejection behaviour of FU Ori type stars. We apply this model to simulations of star formation, invoking two types of initial conditions: spherically symmetric cores in solid-body rotation with $ρ\propto r^{-2}$, and spherically symmetric turbulent cores with density proportional to the density of a Bonnor-Ebert sphere. For a wide range of model parameters, we find that episodic outflows lead to self-regulation of the ejected mass and momentum, and we achieve acceptable results, even with relatively low resolution. Using this model, we find that recently ejected outflow bullets produce a `Hubble wedge' in the position-velocity relation. However, once such a bullet hits the leading shock front, it decelerates and aligns with older bullets to form a `Hubble-law'. Bullets can be identified as bumps in the mass-velocity relation, which can be fit with a power-law, $dM/d\upsilon_{_{\rm RAD}}\propto\upsilon_{_{\rm RAD}}^{-1.5}$. △ Less

Submitted 10 December, 2018; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: Accepted (MNRAS) - 2018 November 30. 19 pages, 10 figures

MSC Class: 85-08

Dust in the Eye of Andromeda

Authors: K. A. Marsh, A. P. Whitworth, M. W. L. Smith, O. Lomax, S. A. Eales

Abstract: We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity in… ▽ More We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity index reported by other recent studies, including the central decrease within ~ 3 kpc of the nucleus. We also investigate the detailed distribution of dust in the nuclear region, a prominent feature of which is a ~ 500 pc bar-like structure seen previously in Hα. The nature of this feature has been the subject of some debate. Our maps show it to be the site of the warmest dust, with a mean line-of-sight temperature ~ 30 K. A comparison with the stellar distribution, based on 2MASS data, provides strong evidence that it is a gravitationally induced bar. A comparison with radial velocity maps suggests the presence of an inflow towards the nucleus from opposite directions along this bar, fed presumably by the nuclear spiral with which it appears to connect. Such behaviour is common in large-scale bars in spiral galaxies, as is the phenomenon of nested bars whereby a subkiloparsec nuclear bar exists within a large-scale primary bar. We suggest that M31 represents an example of such nesting. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: 11 pages, 8 figures. Accepted for publication in MNRAS

Synthetic C$^{18}$O observations of fibrous filaments: the problems of mapping from PPV to PPP

Authors: S. D. Clarke, A. P. Whitworth, R. L. Spowage, A. Duarte-Cabral, S. T. Suri, S. E. Jaffa, S. Walch, P. C. Clark

Abstract: Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a tur… ▽ More Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a turbulent medium. Synthetic C$^{18}$O observations of the simulated filaments reveal the existence of fibres very similar to the observed ones, i.e. elongated coherent features in the resulting PPV data-cubes. Analysis of the PPP data-cubes (i.e. 3D density fields) also reveals elongated coherent features, which we term sub-filaments. Unfortunately there is very poor correspondence between the fibres and the sub-filaments in the simulations. Both fibres and sub-filaments derive from inhomogeneities in the turbulent accretion flow onto the main filament. As a consequence, fibres are often affected by line-of-sight confusion. Similarly, sub-filaments are often affected by large velocity gradients, and even velocity discontinuities. These results suggest that extreme care should be taken when using velocity coherent features to constrain the underlying substructure within a filament. △ Less

Submitted 25 June, 2018; v1 submitted 22 June, 2018; originally announced June 2018.

Comments: 15 pages + appendices, 26 figures. Accepted to MNRAS

Magnetic fields towards Ophiuchus-B derived from SCUBA-2 polarization measurements

Authors: Archana Soam, Kate Pattle, Derek Ward-Thompson, Chang Won Lee, Sarah Sadavoy, Patrick M. Koch, Gwanjeong Kim, Jungmi Kwon, Woojin Kwon, Doris Arzoumanian, David Berry, Thiem Hoang, Motohide Tamura, Sang-Sung Lee, Tie Liu, Kee-Tae Kim, Doug Johnstone, Fumitaka Nakamura, A-Ran Lyo, Takashi Onaka, Jongsoo Kim, Ray S. Furuya, Tetsuo Hasegawa, Shih-Ping Lai, Pierre Bastien , et al. (99 additional authors not shown)

Abstract: We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fi… ▽ More We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $μ$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $λ$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical. △ Less

Submitted 7 June, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: 17 pages, 10 figures, Accepted for publication in ApJ

A First Look at BISTRO Observations of The $ρ$ Oph-A core

Authors: Jungmi Kwon, Yasuo Doi, Motohide Tamura, Masafumi Matsumura, Kate Pattle, David Berry, Sarah Sadavoy, Brenda C. Matthews, Derek Ward-Thompson, Tetsuo Hasegawa, Ray S. Furuya, Andy Pon, James Di Francesco, Doris Arzoumanian, Saeko S. Hayashi, Koji S. Kawabata, Takashi Onaka, Minho Choi, Miju Kang, Thiem Hoang, Chang Won Lee, Sang-Sung Lee, Hong-Li Liu, Tie Liu, Shu-Ichiro Inutsuka , et al. (97 additional authors not shown)

Abstract: We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We iden… ▽ More We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $ρ$ Oph A core; some of them can be part of a coherent magnetic field structure in the $ρ$ Oph region. The results are consistent with previous observations of the brightest regions of $ρ$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 22 pages, 11 figures, accepted for publication in ApJ

Modelling the structure of star clusters with fractional Brownian motion

Authors: O. Lomax, M. L. Bates, A. P. Whitworth

Abstract: The degree of fractal substructure in molecular clouds can be quantified by comparing them with Fractional Brownian Motion (FBM) surfaces or volumes. These fields are self-similar over all length scales and characterised by a drift exponent $H$, which describes the structural roughness. Given that the structure of molecular clouds and the initial structure of star clusters are almost certainly lin… ▽ More The degree of fractal substructure in molecular clouds can be quantified by comparing them with Fractional Brownian Motion (FBM) surfaces or volumes. These fields are self-similar over all length scales and characterised by a drift exponent $H$, which describes the structural roughness. Given that the structure of molecular clouds and the initial structure of star clusters are almost certainly linked, it would be advantageous to also apply this analysis to clusters. Currently, the structure of star clusters is often quantified by applying $\mathcal{Q}$ analysis. $\mathcal{Q}$ values from observed targets are interpreted by comparing them with those from artificial clusters. These are typically generated using a Box-Fractal (BF) or Radial Density Profile (RDP) model. We present a single cluster model, based on FBM, as an alternative to these models. Here, the structure is parameterised by $H$, and the standard deviation of the log-surface/volume density $σ$. The FBM model is able to reproduce both centrally concentrated and substructured clusters, and is able to provide a much better match to observations than the BF model. We show that $\mathcal{Q}$ analysis is unable to estimate FBM parameters. Therefore, we develop and train a machine learning algorithm which can estimate values of $H$ and $σ$, with uncertainties. This provides us with a powerful method for quantifying the structure of star clusters in terms which relate to the structure of molecular clouds. We use the algorithm to estimate the $H$ and $σ$ for several young star clusters, some of which have no measurable BF or RDP analogue. △ Less

Submitted 17 May, 2019; v1 submitted 18 April, 2018; originally announced April 2018.

Comments: 11 Pages, accepted by MNRAS

The unexpectedly large proportion of high-mass star-forming cores in a Galactic mini-starburst

Authors: F. Motte, T. Nony, F. Louvet, K. A. Marsh, S. Bontemps, A. P. Whitworth, A. Men'shchikov, Q. Nguyen Luong, T. Csengeri, A. J. Maury, A. Gusdorf, E. Chapillon, V. Konyves, P. Schilke, A. Duarte-Cabral, P. Didelon, M. Gaudel

Abstract: Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation… ▽ More Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 20 pages, 6 Figures, 2 Tables Nature Astronomy, 2018

Journal ref: Nature Astronomy 2018

J plots: a new method for characterizing structures in the interstellar medium

Authors: S. E. Jaffa, A. P. Whitworth, S. D. Clarke, A. D. P. Howard

Abstract: Large scale surveys have brought about a revolution in astronomy. To analyse the resulting wealth of data, we need automated tools to identify, classify, and quantify the important underlying structures. We present here a method for classifying and quantifying a pixelated structure, based on its principal moments of inertia. The method enables us to automatically detect, and objectively compare, c… ▽ More Large scale surveys have brought about a revolution in astronomy. To analyse the resulting wealth of data, we need automated tools to identify, classify, and quantify the important underlying structures. We present here a method for classifying and quantifying a pixelated structure, based on its principal moments of inertia. The method enables us to automatically detect, and objectively compare, centrally condensed cores, elongated filaments and hollow rings. We illustrate the method by applying it to (i) observations of surface-density from Hi-GAL, and (ii) simulations of filament growth in a turbulent medium. We limit the discussion here to 2D data; in a future paper we will extend the method to 3D data. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 10 pages, 9 figures

Synthetic observations of protostellar multiple systems

Authors: O. Lomax, A. P. Whitworth

Abstract: Observations of protostars are often compared with synthetic observations of models in order to infer the underlying physical properties of the protostars. The majority of these models have a single protostar, attended by a disc and an envelope. However, observational and numerical evidence suggests that a large fraction of protostars form as multiple systems. This means that fitting models of sin… ▽ More Observations of protostars are often compared with synthetic observations of models in order to infer the underlying physical properties of the protostars. The majority of these models have a single protostar, attended by a disc and an envelope. However, observational and numerical evidence suggests that a large fraction of protostars form as multiple systems. This means that fitting models of single protostars to observations may be inappropriate. We produce synthetic observations of protostellar multiple systems undergoing realistic, non-continuous accretion. These systems consist of multiple protostars with episodic luminosities, embedded self-consistently in discs and envelopes. We model the gas dynamics of these systems using smoothed particle hydrodynamics and we generate synthetic observations by post-processing the snapshots using the \textsc{spamcart} Monte Carlo radiative transfer code. We present simulation results of three model protostellar multiple systems. For each of these, we generate $4\times10^4$ synthetic spectra at different points in time and from different viewing angles. We propose a Bayesian method, using similar calculations to those presented here, but in greater numbers, to infer the physical properties of protostellar multiple systems from observations. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 13 pages. Accepted by MNRAS

Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths

Authors: Richard Wünsch, Stefanie Walch, František Dinnbier, Anthony P. Whitworth

Abstract: We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this appr… ▽ More We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later. △ Less

Submitted 2 January, 2018; v1 submitted 21 August, 2017; originally announced August 2017.

Comments: 26 pages, 19 figures, 6 tables. Accepted by MNRAS

Multi-temperature mapping of dust structures throughout the Galactic Plane using the PPMAP tool with Herschel Hi-GAL data

Authors: K. A. Marsh, A. P. Whitworth, O. Lomax, S. E. Ragan, U. Becciani, L. Cambresy, A. Di Giorgio, D. Eden, D. Elia, P. Kacsuk, S. Molinari, P. Palmeirim, S. Pezzuto, N. Schneider, E. Sciacca, F. Vitello

Abstract: We describe new Hi-GAL based maps of the entire Galactic Plane, obtained using continuum data in the wavelength range 70-500 $μ$m. These maps are derived with the PPMAP procedure, and therefore represent a significant improvement over those obtained with standard analysis techniques. Specifically they have greatly improved resolution (12 arcsec) and, in addition to more accurate integrated column… ▽ More We describe new Hi-GAL based maps of the entire Galactic Plane, obtained using continuum data in the wavelength range 70-500 $μ$m. These maps are derived with the PPMAP procedure, and therefore represent a significant improvement over those obtained with standard analysis techniques. Specifically they have greatly improved resolution (12 arcsec) and, in addition to more accurate integrated column densities and mean dust temperatures, they give temperature-differential column densities, i.e., separate column density maps in twelve distinct dust temperature intervals, along with the corresponding uncertainty maps. The complete set of maps is available on-line. We briefly describe PPMAP and present some illustrative examples of the results. These include (a) multi-temperature maps of the Galactic HII region W5-E, (b) the temperature decomposition of molecular cloud column-density probability distribution functions, and (c) the global variation of mean dust temperature as a function of Galactocentric distance. Amongst our findings are: (i) a strong localised temperature gradient in W5-E in a direction orthogonal to that towards the ionising star, suggesting an alternative heating source and providing possible guidance for models of the formation of the bubble complex, and (ii) the overall radial profile of dust temperature in the Galaxy shows a monotonic decrease, broadly consistent both with models of the interstellar radiation field and with previous estimates at lower resolution. However, we also find a central temperature plateau within ~ 6 kpc of the Galactic centre, outside of which is a pronounced steepening of the radial profile. This behaviour may reflect the greater proportion of molecular (as opposed to atomic) gas in the central region of the Galaxy. △ Less

Submitted 12 July, 2017; originally announced July 2017.

Comments: 14 pages, 12 figures; accepted for publication in MNRAS

Properties of Hi-GAL clumps in the inner Galaxy]{The Hi-GAL compact source catalogue. I. The physical properties of the clumps in the inner Galaxy ($-71.0^{\circ}< \ell < 67.0^{\circ}$)

Authors: D. Elia, S. Molinari, E. Schisano, M. Pestalozzi, S. Pezzuto, M. Merello, A. Noriega-Crespo, T. J. T. Moore, D. Russeil, J. C. Mottram, R. Paladini, F. Strafella, M. Benedettini, J. P. Bernard, A. Di Giorgio, D. J. Eden, Y. Fukui, R. Plume, J. Bally, P. G. Martin, S. E. Ragan, S. E. Jaffa, F. Motte, L. Olmi, N. Schneider , et al. (61 additional authors not shown)

Abstract: Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $μ$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the port… ▽ More Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $μ$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $μ$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions. △ Less

Submitted 4 June, 2017; originally announced June 2017.

Comments: Accepted by MNRAS

First results from BISTRO -- a SCUBA-2 polarimeter survey of the Gould Belt

Authors: Derek Ward-Thompson, Kate Pattle, Pierre Bastien, Ray S. Furuya, Woojin Kwon, Shih-Ping Lai, Keping Qiu, David Berry, Minho Choi, Simon Coudé, James Di Francesco, Thiem Hoang, Erica Franzmann, Per Friberg, Sarah F. Graves, Jane S. Greaves, Martin Houde, Doug Johnstone, Jason M. Kirk, Patrick M. Koch, Jungmi Kwon, Chang Won Lee, Di Li, Brenda C. Matthews, Joseph C. Mottram , et al. (89 additional authors not shown)

Abstract: We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will a… ▽ More We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 12 pages, 8 figures, accepted for publication in ApJ

Filamentary fragmentation in a turbulent medium

Authors: S. D. Clarke, A. P. Whitworth, A. Duarte-Cabral, D. A. Hubber

Abstract: We present the results of smoothed particle hydrodynamic simulations investigating the evolution and fragmentation of filaments that are accreting from a turbulent medium. We show that the presence of turbulence, and the resulting inhomogeneities in the accretion flow, play a significant role in the fragmentation process. Filaments which experience a weakly turbulent accretion flow fragment in a t… ▽ More We present the results of smoothed particle hydrodynamic simulations investigating the evolution and fragmentation of filaments that are accreting from a turbulent medium. We show that the presence of turbulence, and the resulting inhomogeneities in the accretion flow, play a significant role in the fragmentation process. Filaments which experience a weakly turbulent accretion flow fragment in a two-tier hierarchical fashion, similar to the fragmentation pattern seen in the Orion Integral Shaped Filament. Increasing the energy in the turbulent velocity field results in more sub-structure within the filaments, and one sees a shift from gravity-dominated fragmentation to turbulence-dominated fragmentation. The sub-structure formed in the filaments is elongated and roughly parallel to the longitudinal axis of the filament, similar to the fibres seen in observations of Taurus, and suggests that the fray and fragment scenario is a possible mechanism for the production of fibres. We show that the formation of these fibre-like structures is linked to the vorticity of the velocity field inside the filament and the filament's accretion from an inhomogeneous medium. Moreover, we find that accretion is able to drive and sustain roughly sonic levels of turbulence inside the filaments, but is not able to prevent radial collapse once the filaments become supercritical. However, the supercritical filaments which contain fibre-like structures do not collapse radially, suggesting that fibrous filaments may not necessarily become radially unstable once they reach the critical line-density. △ Less

Submitted 13 March, 2017; originally announced March 2017.

Comments: (Accepted for publication in MNRAS)

Fragmentation of vertically stratified gaseous layers: monolithic or coalescence-driven collapse

Authors: František Dinnbier, Richard Wünsch, Anthony P. Whitworth, Jan Palouš

Abstract: We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by… ▽ More We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by Elmegreen & Elmegreen (1978), and consequently deviates from the dispersion relations based on the thin shell approximation (Vishniac 1983) or pressure assisted gravitational instability (Wünsch et al. 2010). In the non-linear regime, the relative importance of the confining pressure to the self-gravity is a crucial parameter controlling the qualitative course of fragmentation. When confinement of the layer is dominated by external pressure, self- gravitating condensations are delivered by a two-stage process: first the layer fragments into gravitationally bound but stable clumps, and then these clumps coalesce until they assemble enough mass to collapse. In contrast, when external pressure makes a small contribution to confinement of the layer, the layer fragments monolithically into gravitationally unstable clumps and there is no coalescence. This dichotomy persists whether the external pressure is thermal or ram. We apply these results to fragments forming in a shell swept up by an expanding H II region, and find that, unless the swept up gas is quite hot or the surrounding medium has low density, the fragments have low-mass ( ~< 3 M_Sun ), and therefore they are unlikely to spawn stars that are sufficiently massive to promote sequential self-propagating star formation. △ Less

Submitted 10 January, 2017; originally announced January 2017.

Comments: 21 pages, 12 figures, MNRAS in press

$\mathcal{Q}^{+}$: Characterising the structure of young star clusters

Authors: S. E. Jaffa, A. P. Whitworth, O. Lomax

Abstract: Many young star clusters appear to be fractal, i.e. they appear to be concentrated in a nested hierarchy of clusters within clusters. We present a new algorithm for statistically analysing the distribution of stars to quantify the level of sub-structure. We suggest that, even at the simplest level, the internal structure of a fractal cluster requires the specification of three parameters. (i) The… ▽ More Many young star clusters appear to be fractal, i.e. they appear to be concentrated in a nested hierarchy of clusters within clusters. We present a new algorithm for statistically analysing the distribution of stars to quantify the level of sub-structure. We suggest that, even at the simplest level, the internal structure of a fractal cluster requires the specification of three parameters. (i) The 3D fractal dimension, $\mathcal{D}$, measures the extent to which the clusters on one level of the nested hierarchy fill the volume of their parent cluster. (ii) The number of levels, $\mathcal{L}$, reflects the finite ratio between the linear size of the large root-cluster at the top of the hierarchy, and the smallest leaf-clusters at the bottom of the hierarchy. (iii) The volume-density scaling exponent, $\mathcal{C}=-\textrm{d}\ln[δn]/\textrm{d}\ln[L]$ measures the factor by which the excess density, $δn$, in a structure of scale $L$, exceeds that of the background formed by larger structures; it is similar, but not exactly equivalent, to the exponent in Larson's scaling relation between density and size for molecular clouds. We describe an algorithm which can be used to constrain the values of $({\cal D},{\cal L},{\cal C})$ and apply this method to artificial and observed clusters. We show that this algorithm is able to reliably describe the three dimensional structure of an artificial star cluster from the two dimensional projection, and quantify the varied structures observed in real and simulated clusters. △ Less

Submitted 7 December, 2016; v1 submitted 30 November, 2016; originally announced November 2016.

Comments: Accepted by MNRAS

A ram-pressure threshold for star formation

Authors: A. P. Whitworth

Abstract: In turbulent fragmentation, star formation occurs in condensations created by converging flows. The condensations must be sufficiently massive, dense and cool to be gravitationally unstable, so that they start to contract; {\it and} they must then radiate away thermal energy fast enough for self-gravity to remain dominant, so that they continue to contract. For the metallicities and temperatures i… ▽ More In turbulent fragmentation, star formation occurs in condensations created by converging flows. The condensations must be sufficiently massive, dense and cool to be gravitationally unstable, so that they start to contract; {\it and} they must then radiate away thermal energy fast enough for self-gravity to remain dominant, so that they continue to contract. For the metallicities and temperatures in local star forming clouds, this second requirement is only met robustly when the gas couples thermally to the dust, because this delivers the capacity to radiate across the full bandwidth of the continuum, rather than just in a few discrete spectral lines. This translates into a threshold for vigorous star formation, which can be written as a minimum ram-pressure Pcrit ~ 4 10^-11 dyne. Pcrit is independent of temperature, and corresponds to flows with molecular hydrogen number-density nH2 and velocity v satisfying nH2 v^2 > 800 cm^-3 (km/s)^2. This in turn corresponds to a minimum molecular hydrogen column-density for vigorous star formation, NH2crit ~ 4 10^21 cm^-2 (SIGMAcrit ~ 100 MSun pc^-2), and a minimum visual extinction AVcrit ~ 9. The characteristic diameter and line-density for a star-forming filament when this threshold is just exceeded -- a sweet spot for local star formation regions -- are 2Rfil ~ 0.1 pc and mufil ~ 13 MSun pc^-2. The characteristic diameter and mass for a prestellar core condensing out of such a filament are 2Rcore ~ 0.1 pc, and Mcore ~ MSun. We also show that fragmentation of a shock-compressed layer is likely to commence while the convergent flows creating the layer are still ongoing, and we stress that, under this circumstance, the phenomenology and characteristic scales for fragmentation of the layer are fundamentally different from those derived traditionally for pre-existing layers. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 20 pages, 2 fugures, published in MNRAS