-
Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations
Authors:
Jia-Wei Wang,
Patrick M. Koch,
Seamus D. Clarke,
Gary Fuller,
Nicolas Peretto,
Ya-Wen Tang,
Hsi-Wei Yen,
Shih-Ping Lai,
Nagayoshi Ohashi,
Doris Arzoumanian,
Doug Johnstone,
Ray Furuya,
Shu-ichiro Inutsuka,
Chang Won Lee,
Derek Ward-Thompson,
Valentin J. M. Le Gouellec,
Hong-Li Liu,
Lapo Fanciullo,
Jihye Hwang,
Kate Pattle,
Frédérick Poidevin,
Mehrnoosh Tahani,
Takashi Onaka,
Mark G. Rawlings,
Eun Jung Chung
, et al. (132 additional authors not shown)
Abstract:
We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from…
▽ More
We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
GMF G214.5-1.8 as traced by CO: I -- cloud-scale CO freeze-out as a result of a low cosmic-ray ionisation rate
Authors:
S. D. Clarke,
V. A. Makeev,
Á. Sánchez-Monge,
G. M. Williams,
Y. -W. Tang,
S. Walch,
R. Higgins,
P. C. Nürnberger,
S. Suri
Abstract:
We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using IRAM 30m data of $^{12}$CO, $^{13}$CO and C$^{18}$O. We find that the $^{12}$CO (1-0) and (2-1) derived excitation temperatures are near identical and are very low, with a median of 8.2 K, showing that the gas is extremely cold across the whole cloud. Investigating the abundance of $^{13}$CO across…
▽ More
We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using IRAM 30m data of $^{12}$CO, $^{13}$CO and C$^{18}$O. We find that the $^{12}$CO (1-0) and (2-1) derived excitation temperatures are near identical and are very low, with a median of 8.2 K, showing that the gas is extremely cold across the whole cloud. Investigating the abundance of $^{13}$CO across G214.5, we find that there is a significantly lower abundance along the entire 13 pc spine of the filament, extending out to a radius of $\sim 0.8$ pc, corresponding to $A_v \gtrsim 2$ mag and $T_{dust} \lesssim 13.5$ K. Due to this, we attribute the decrease in abundance to CO freeze-out, making G214.5 the largest scale example of freeze-out yet. We construct an axisymmetric model of G214.5's $^{13}$CO volume density considering freeze-out and find that to reproduce the observed profile significant depletion is required beginning at low volume densities, $n\gtrsim2000$ cm$^{-3}$. Freeze-out at this low number density is possible only if the cosmic-ray ionisation rate is $\sim 1.9 \times 10^{-18}$ s$^{-1}$, an order of magnitude below the typical value. Using timescale arguments, we posit that such a low ionisation rate may lead to ambipolar diffusion being an important physical process along G214.5's entire spine. We suggest that if low cosmic-ray ionisation rates are more common in the outer Galaxy, and other quiescent regions, cloud-scale CO freeze-out occurring at low column and number densities may also be more prevalent, having consequences for CO observations and their interpretation.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
Unravelling the structure of magnetised molecular clouds with SILCC-Zoom: sheets, filaments and fragmentation
Authors:
S. Ganguly,
S. Walch,
D. Seifried,
S. D. Clarke,
M. Weis
Abstract:
To what extent magnetic fields affect how molecular clouds (MCs) fragment and create dense structures is an open question. We present a numerical study of cloud fragmentation using the SILCC-Zoom simulations. These simulations follow the self-consistent formation of MCs in a few hundred parsec sized region of a stratified galactic disc; and include magnetic fields, self-gravity, supernova-driven t…
▽ More
To what extent magnetic fields affect how molecular clouds (MCs) fragment and create dense structures is an open question. We present a numerical study of cloud fragmentation using the SILCC-Zoom simulations. These simulations follow the self-consistent formation of MCs in a few hundred parsec sized region of a stratified galactic disc; and include magnetic fields, self-gravity, supernova-driven turbulence, as well as a non-equilibrium chemical network. To discern the role of magnetic fields in the evolution of MCs, we study seven simulated clouds, five with magnetic fields, and two without, with a maximum resolution of 0.1 parsec. Using a dendrogram we identify hierarchical structures which form within the clouds. Overall, the magnetised clouds have more mass in a diffuse envelope with a number density between 1-100 cm$^{-3}$. We find that six out of seven clouds are sheet-like on the largest scales, as also found in recent observations, and with filamentary structures embedded within, consistent with the bubble-driven MC formation mechanism. Hydrodynamic simulations tend to produce more sheet-like structures also on smaller scales, while the presence of magnetic fields promotes filament formation. Analysing cloud energetics, we find that magnetic fields are dynamically important for less dense, mostly but not exclusively atomic structures (typically up to $\sim 100 - 1000$~cm$^{-3}$), while the denser, potentially star-forming structures are energetically dominated by self-gravity and turbulence. In addition, we compute the magnetic surface term and demonstrate that it is generally confining, and some atomic structures are even magnetically held together. In general, magnetic fields delay the cloud evolution and fragmentation by $\sim$ 1 Myr.
△ Less
Submitted 17 July, 2023;
originally announced July 2023.
-
Unveiling the formation of the massive DR21 ridge
Authors:
L. Bonne,
S. Bontemps,
N. Schneider,
R. Simon,
S. D. Clarke,
T. Csengeri,
E. Chambers,
U. Graf,
J. M. Jackson,
R. Klein,
Y. Okada,
A. G. G. M. Tielens,
M. Tiwari
Abstract:
We present new $^{13}$CO(1-0), C$^{18}$O(1-0), HCO$^{+}$(1-0) and H$^{13}$CO$^{+}$(1-0) maps from the IRAM 30m telescope, and a spectrally-resolved [CII] 158 $μ$m map observed with the SOFIA telescope towards the massive DR21 cloud. This traces the kinematics from low- to high-density gas in the cloud which allows to constrain the formation scenario of the high-mass star forming DR21 ridge. The mo…
▽ More
We present new $^{13}$CO(1-0), C$^{18}$O(1-0), HCO$^{+}$(1-0) and H$^{13}$CO$^{+}$(1-0) maps from the IRAM 30m telescope, and a spectrally-resolved [CII] 158 $μ$m map observed with the SOFIA telescope towards the massive DR21 cloud. This traces the kinematics from low- to high-density gas in the cloud which allows to constrain the formation scenario of the high-mass star forming DR21 ridge. The molecular line data reveals that the sub-filaments are systematically redshifted relative to the dense ridge. We demonstrate that [CII] unveils the surrounding CO-poor gas of the dense filaments in the DR21 cloud. We also show that this surrounding gas is organized in a flattened cloud with curved redshifted dynamics perpendicular to the ridge. The sub-filaments thus form in this curved and flattened mass reservoir. A virial analysis of the different lines indicates that self-gravity should drive the evolution of the ridge and surrounding cloud. Combining all results we propose that bending of the magnetic field, due to the interaction with a mostly atomic colliding cloud, explains the velocity field and resulting mass accretion on the ridge. This is remarkably similar to what was found for at least two nearby low-mass filaments. We tentatively propose that this scenario might be a widespread mechanism to initiate star formation in the Milky Way. However, in contrast to low-mass clouds, gravitational collapse plays a role on the pc scale of the DR21 ridge because of the higher density. This allows more effective mass collection at the centers of collapse and should facilitate massive cluster formation.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
How do supernova remnants cool? -- I. Morphology, optical emission lines, and shocks
Authors:
Ekaterina I. Makarenko,
Stefanie Walch,
Seamus D. Clarke,
Daniel Seifried,
Thorsten Naab,
Pierre C. Nürnberger,
Tim-Eric Rathjen
Abstract:
Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SN…
▽ More
Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ($[$O III$]$, $[$N II$]$, $[$S II$]$, H$α$, H$β$) are usually best observable. Our shock detection scheme shows that [S II] and [N II] emissions arise from the thin shell surrounding the SNR, while [O III], H$\rm α$, and H$\rm β$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR's complex structure and its projection onto the plane of the sky because the escaping line luminosity can be reduced by 10 -- 80\% due to absorption along the line-of-sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
A Herschel study of G214.5-1.8: a young, cold and quiescent giant molecular filament on the shell of a HI superbubble
Authors:
S. D. Clarke,
A. Sanchez-Monge,
G. M. Williams,
A. D. P. Howard,
S. Walch,
N. Schneider
Abstract:
We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using Herschel data. We find that G214.5 has a mass of $\sim$ 16,000 M$_{\odot}$, yet hosts only 15 potentially protostellar 70 $μ$m sources, making it highly quiescent compared to equally massive clouds such as Serpens and Mon R2. We show that G214.5 has a unique morphology, consisting of a narrow `Main…
▽ More
We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using Herschel data. We find that G214.5 has a mass of $\sim$ 16,000 M$_{\odot}$, yet hosts only 15 potentially protostellar 70 $μ$m sources, making it highly quiescent compared to equally massive clouds such as Serpens and Mon R2. We show that G214.5 has a unique morphology, consisting of a narrow `Main filament' running north-south and a perpendicular `Head' structure running east-west. We identify 33 distinct massive clumps from the column density maps, 8 of which are protostellar. However, the star formation activity is not evenly spread across G214.5 but rather predominantly located in the Main filament. Studying the Main filament in a manner similar to previous works, we find that G214.5 is most like a 'Bone' candidate GMF, highly elongated and massive, but it is colder and narrower than any such GMF. It also differs significantly due to its low fraction of high column density gas. Studying the radial profile, we discover that G214.5 is highly asymmetric and resembles filaments which are known to be compressed externally. Considering its environment, we find that G214.5 is co-incident, spatially and kinematically, with a HI superbubble. We discuss how a potential interaction between G214.5 and the superbubble may explain G214.5's morphology, asymmetry and, paucity of dense gas and star formation activity, highlighting the intersection of a bubble-driven interstellar medium paradigm with that of a filament paradigm for star formation.
△ Less
Submitted 4 November, 2022;
originally announced November 2022.
-
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.
-
From Bubbles and Filaments to Cores and Disks: Gas Gathering and Growth of Structure Leading to the Formation of Stellar Systems
Authors:
Jaime E. Pineda,
Doris Arzoumanian,
Philippe André,
Rachel K. Friesen,
Annie Zavagno,
Seamus D. Clarke,
Tsuyoshi Inoue,
Che-Yu Chen,
Yueh-Ning Lee,
Juan D. Soler,
Michael Kuffmeier
Abstract:
The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented…
▽ More
The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50-100 pc scales through parsec-scale molecular filaments down to $<$0.1 pc dense cores and $<$1000 au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation.
△ Less
Submitted 23 March, 2023; v1 submitted 8 May, 2022;
originally announced May 2022.
-
SILCC-Zoom: the dynamic balance in molecular cloud substructures
Authors:
S. Ganguly,
S. Walch,
S. D. Clarke,
D. Seifried
Abstract:
How molecular clouds fragment and create the dense structures which go on to form stars is an open question. We investigate the relative importance of different energy terms (kinetic, thermal, magnetic, and gravity - both self-gravity and tidal forces) for the formation and evolution of molecular clouds and their sub-structures based on the SILCC-Zoom simulations. These simulations follow the self…
▽ More
How molecular clouds fragment and create the dense structures which go on to form stars is an open question. We investigate the relative importance of different energy terms (kinetic, thermal, magnetic, and gravity - both self-gravity and tidal forces) for the formation and evolution of molecular clouds and their sub-structures based on the SILCC-Zoom simulations. These simulations follow the self-consistent formation of cold molecular clouds down to scales of 0.1 pc from the diffuse supernova-driven interstellar medium in a stratified galactic disc. We study the time evolution of seven molecular clouds (five with magnetic fields and two without) for 1.5-2 Myr. Using a dendrogram, we identify hierarchical 3D sub-structures inside the clouds with the aim to understand their dynamics and distinguish between the theories of gravo-turbulent fragmentation and global hierarchical collapse. The virial analysis shows that the dense gas is indeed dominated by the interplay of gravity and turbulence, while magnetic fields and thermal pressure are only important for fluffy, atomic structures. Over time, gravitationally bound sub-structures emerge from a marginally bound medium (viral ratio $1 \leq α_{\rm vir}^{\rm vol} <2$) as a result of large-scale supernova-driven inflows rather than global collapse. A detailed tidal analysis shows that the tidal tensor is highly anisotropic. Yet the tidal forces are generally not strong enough to disrupt either large-scale or dense sub-structures but cause their deformation. By comparing tidal and crossing time scales, we find that tidal forces do not seem to be the main driver of turbulence within the molecular clouds.
△ Less
Submitted 5 April, 2022;
originally announced April 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.
-
Synthetic CO emission and the $X_{\rm CO}$ factor of young molecular clouds: a convergence study
Authors:
Elisabeth M. A. Borchert,
Stefanie Walch,
Daniel Seifried,
Seamus D. Clarke,
Annika Franeck,
Pierre Nürnberger
Abstract:
The properties of synthetic CO emission from 3D simulations of forming molecular clouds are studied within the SILCC-Zoom project. Since the time scales of cloud evolution and molecule formation are comparable, the simulations include a live chemical network. Two sets of simulations with an increasing spatial resolution (d$x=3.9$ pc to d$x=0.06$ pc) are used to investigate the convergence of the s…
▽ More
The properties of synthetic CO emission from 3D simulations of forming molecular clouds are studied within the SILCC-Zoom project. Since the time scales of cloud evolution and molecule formation are comparable, the simulations include a live chemical network. Two sets of simulations with an increasing spatial resolution (d$x=3.9$ pc to d$x=0.06$ pc) are used to investigate the convergence of the synthetic CO emission, which is computed by post-processing the simulation data with the RADMC-3D radiative transfer code. To determine the excitation conditions, it is necessary to include atomic hydrogen and helium alongside H$_2$, which increases the resulting CO emission by ~7-26 per cent. Combining the brightness temperature of $^{12}$CO and $^{13}$CO, we compare different methods to estimate the excitation temperature, the optical depth of the CO line and hence, the CO column density. An intensity-weighted average excitation temperature results in the most accurate estimate of the total CO mass. When the pixel-based excitation temperature is used to calculate the CO mass, it is over-/underestimated at low/high CO column densities where the assumption that $^{12}$CO is optically thick while $^{13}$CO is optically thin is not valid. Further, in order to obtain a converged total CO luminosity and hence <$X_{\rm CO}$> factor, the 3D simulation must have d$x\lesssim0.1$ pc. The <$X_{\rm CO}$> evolves over time and differs for the two clouds; yet pronounced differences with numerical resolution are found. Since high column density regions with a visual extinction larger than 3~mag are not resolved for d$x\gtrsim 1$~pc, in this case the H$_2$ mass and CO luminosity both differ significantly from the higher resolution results and the local $X_{\rm CO}$ is subject to strong noise. Our calculations suggest that synthetic CO emission maps are only converged for simulations with d$x\lesssim 0.1$ pc.
△ Less
Submitted 19 November, 2021; v1 submitted 1 February, 2021;
originally announced February 2021.
-
Formation of the Musca filament: Evidence for asymmetries in the accretion flow due to a cloud-cloud collision
Authors:
L. Bonne,
S. Bontemps,
N. Schneider,
S. D. Clarke,
D. Arzoumanian,
Y. Fukui,
K. Tachihara,
T. Csengeri,
R. Guesten,
A. Ohama,
R. Okamoto,
R. Simon,
H. Yahia,
H. Yamamoto
Abstract:
Context. Dense molecular filaments are ubiquituous in the interstellar medium, yet their internal physical conditions and formation mechanism remain debated. Aims. We study the kinematics and physical conditions in the Musca filament and the Chamaeleon-Musca complex to constrain the physics of filament formation. Methods. We produced CO(2-1) isotopologue maps with the APEX telescope that cut throu…
▽ More
Context. Dense molecular filaments are ubiquituous in the interstellar medium, yet their internal physical conditions and formation mechanism remain debated. Aims. We study the kinematics and physical conditions in the Musca filament and the Chamaeleon-Musca complex to constrain the physics of filament formation. Methods. We produced CO(2-1) isotopologue maps with the APEX telescope that cut through the Musca filament. We further study a NANTEN2 $^{12}$CO(1-0) map of the Musca cloud and the HI emission of the Chamaeleon-Musca complex. Results. The Musca cloud contains multiple velocity components. Radiative transfer modelling of the CO emission indicates that the Musca filament consists of a cold ($\sim$10 K), dense (n$_{H_2}\sim$10$^4$ cm$^{-3}$) crest, which is best described with a cylindrical geometry. Connected to the crest, a separate gas component at T$\sim$15 K and n$_{H_2}\sim$10$^3$ cm$^{-3}$ is found, the so-called strands. The filament crest has a transverse velocity gradient that is linked to the kinematics of the nearby ambient cloud. Studying the large scale kinematics, we find coherence of the asymmetric kinematics from the 50 pc HI cloud down to the Musca filament. We also report a strong [C$^{18}$O]/[$^{13}$CO] abundance drop by an order of magnitude from the filament crest to the strands over a distance $<$ 0.2 pc in a weak far-ultraviolet (FUV) field. Conclusions. The dense Musca filament crest is a long-lived (several crossing times), dynamic structure that can form stars in the near future because of continuous mass accretion. This mass accretion appears to be triggered by a HI cloud-cloud collision, which bends the magnetic field around dense filaments. This bending of the magnetic field is then responsible for the observed asymmetric accretion scenario of the Musca filament, which is, for instance, seen as a V-shape in the position-velocity (PV) diagram.
△ Less
Submitted 23 October, 2020;
originally announced October 2020.
-
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.
-
Dense gas formation in the Musca filament due to the dissipation of a supersonic converging flow
Authors:
L. Bonne,
N. Schneider,
S. Bontemps,
S. D. Clarke,
A. Gusdorf,
A. Lehmann,
M. Steinke,
T. Csengeri,
S. Kabanovic,
R. Simon,
C. Buchbender,
R. Guesten
Abstract:
Observations with the Herschel Space Telescope have established that most of the star forming gas is organised in interstellar filaments, a finding that is supported by numerical simulations of the supersonic interstellar medium (ISM) where dense filamentary structures are ubiquitous. We aim to understand the formation of these dense structures by performing observations covering the $^{12}$CO(4-3…
▽ More
Observations with the Herschel Space Telescope have established that most of the star forming gas is organised in interstellar filaments, a finding that is supported by numerical simulations of the supersonic interstellar medium (ISM) where dense filamentary structures are ubiquitous. We aim to understand the formation of these dense structures by performing observations covering the $^{12}$CO(4-3), $^{12}$CO(3-2), and various CO(2-1) isotopologue lines of the Musca filament, using the APEX telescope. The observed CO intensities and line ratios cannot be explained by PDR (photodissociation region) emission because of the low ambient far-UV field that is strongly constrained by the non-detections of the [C II] line at 158 $μ$m and the [O I] line at 63 $μ$m, observed with the upGREAT receiver on SOFIA, as well as a weak [C I] 609 $μ$m line detected with APEX. We propose that the observations are consistent with a scenario in which shock excitation gives rise to warm and dense gas close to the highest column density regions in the Musca filament. Using shock models, we find that the CO observations can be consistent with excitation by J-type low-velocity shocks. A qualitative comparison of the observed CO spectra with synthetic observations of dynamic filament formation simulations shows a good agreement with the signature of a filament accretion shock that forms a cold and dense filament from a converging flow. The Musca filament is thus found to be dense molecular post-shock gas. Filament accretion shocks that dissipate the supersonic kinetic energy of converging flows in the ISM may thus play a prominent role in the evolution of cold and dense filamentary structures.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
The hierarchical fragmentation of filaments and the role of sub-filaments
Authors:
S. D. Clarke,
G. M. Williams,
S. Walch
Abstract:
Recent observations have revealed the presence of small fibres or sub-filaments within larger filaments. We present a numerical fragmentation study of fibrous filaments investigating the link between cores and sub-filaments using hydrodynamical simulations performed with the moving-mesh code Arepo. Our study suggests that cores form in two environments: (i) as isolated cores, or small chains of co…
▽ More
Recent observations have revealed the presence of small fibres or sub-filaments within larger filaments. We present a numerical fragmentation study of fibrous filaments investigating the link between cores and sub-filaments using hydrodynamical simulations performed with the moving-mesh code Arepo. Our study suggests that cores form in two environments: (i) as isolated cores, or small chains of cores, on a single sub-filament, or (ii) as an ensemble of cores located at the junction of sub-filaments. We term these isolated and hub cores respectively. We show that these core populations are statistically different from each other. Hub cores have a greater mean mass than isolated cores, and the mass distribution of hub cores is significantly wider than isolated cores. This fragmentation is reminiscent of parsec-scale hub-filament systems, showing that the combination of turbulence and gravity leads to similar fragmentation signatures on multiple scales, even within filaments. Moreover, the fact that fragmentation proceeds through sub-filaments suggests that there exists no characteristic fragmentation length-scale between cores. This is in opposition to earlier theoretical works studying fibre-less filaments which suggest a strong tendency towards the formation of quasi-periodically spaced cores, but in better agreement with observations. We also show tentative signs that global collapse of filaments preferentially form cores at both filament ends, which are more massive and dense than other cores.
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
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.
-
Determining the presence of characteristic fragmentation length-scales in filaments
Authors:
S. D. Clarke,
G. M. Williams,
J. C. Ibáñez-Mejía,
S. Walch
Abstract:
Theories suggest that filament fragmentation should occur on a characteristic fragmentation length-scale. This fragmentation length-scale can be related to filament properties, such as the width and the dynamical state of the filament. Here we present a study of a number of fragmentation analysis techniques applied to filaments, and their sensitivity to characteristic fragmentation length-scales.…
▽ More
Theories suggest that filament fragmentation should occur on a characteristic fragmentation length-scale. This fragmentation length-scale can be related to filament properties, such as the width and the dynamical state of the filament. Here we present a study of a number of fragmentation analysis techniques applied to filaments, and their sensitivity to characteristic fragmentation length-scales. We test the sensitivity to both single-tier and two-tier fragmentation, i.e. when the fragmentation can be characterised with one or two fragmentation length-scales respectively. The nearest neighbour separation, minimum spanning tree separation and two-point correlation function are all able to robustly detect characteristic fragmentation length-scales. The Fourier power spectrum and the Nth nearest neighbour technique are both poor techniques, and require very little scatter in the core spacings for the characteristic length-scale to be successfully determined. We develop a null hypothesis test to compare the results of the nearest neighbour and minimum spanning tree separation distribution with randomly placed cores. We show that a larger number of cores is necessary to successfully reject the null hypothesis if the underlying fragmentation is two-tier, N>20. Once the null is rejected we show how one may decide if the observed fragmentation is best described by single-tier or two-tier fragmentation, using either Akaike's information criterion or the Bayes factor. The analysis techniques, null hypothesis tests, and model selection approaches are all included in a new open-source Python/C library called FragMent.
△ Less
Submitted 18 January, 2019;
originally announced January 2019.
-
The CARMA-NRO Orion Survey: The filamentary structure as seen in C$^{18}$O emission
Authors:
S. T. Suri,
A. Sanchez-Monge,
P. Schilke,
S. D. Clarke,
R. J. Smith,
V. Ossenkopf-Okada,
R. Klessen,
P. Padoan,
P. Goldsmith,
H. G. Arce,
J. Bally,
J. M. Carpenter,
A. Ginsburg,
D. Johnstone,
J. Kauffmann,
S. Kong,
D. C. Lis,
S. Mairs,
T. Pillai,
J. E. Pineda,
A. Duarte-Cabral
Abstract:
We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds.…
▽ More
We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.
△ Less
Submitted 1 January, 2019;
originally announced January 2019.
-
On the resolution requirements for modelling molecular gas formation in solar neighbourhood conditions
Authors:
P. R. Joshi,
S. Walch,
D. Seifried,
S. C. O. Glover,
S. D. Clarke,
M. Weis
Abstract:
The formation of molecular hydrogen (H$_2$) and carbon monoxide (CO) is sensitive to the volume and column density distribution of the turbulent interstellar medium. In this paper, we study H$_2$ and CO formation in a large set of hydrodynamical simulations of periodic boxes with driven supersonic turbulence, as well as in colliding flows with the \textsc{Flash} code. The simulations include a non…
▽ More
The formation of molecular hydrogen (H$_2$) and carbon monoxide (CO) is sensitive to the volume and column density distribution of the turbulent interstellar medium. In this paper, we study H$_2$ and CO formation in a large set of hydrodynamical simulations of periodic boxes with driven supersonic turbulence, as well as in colliding flows with the \textsc{Flash} code. The simulations include a non-equilibrium chemistry network, gas self-gravity, and diffuse radiative transfer. We investigate the spatial resolution required to obtain a converged H$_2$ and CO mass fraction and formation history. From the numerical tests we find that H$_2$ converges at a spatial resolution of $\lesssim0.2$~pc, while the required resolution for CO convergence is $\lesssim 0.04$~pc in gas with solar metallicity which is subject to a solar neighbourhood interstellar radiation field. We derive two critical conditions from our numerical results: the simulation has to at least resolve the densities at which (1) the molecule formation time in each cell in the computational domain is equal to the dissociation time, and (2) the formation time is equal to the the typical cell crossing time. For both H$_2$ and CO, the second criterion is more restrictive. The formulae we derive can be used to check whether molecule formation is converged in any given simulation.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
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.
-
Synthetic [CII] emission maps of a simulated molecular cloud in formation
Authors:
A. Franeck,
S. Walch,
D. Seifried,
S. D. Clarke,
V. Ossenkopf-Okada,
S. C. O. Glover,
R. S. Klessen,
P. Girichidis,
T. Naab,
R. Wünsch,
P. C. Clark,
E. Pellegrini,
T. Peters
Abstract:
The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of s…
▽ More
The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{\rm CII}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{\rm CII}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{\rm CII}$ to change in environments with a lower or higher radiation field than simulated here.
△ Less
Submitted 27 September, 2018;
originally announced September 2018.
-
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.
-
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.
-
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.
-
Perturbation growth in accreting filaments
Authors:
Seamus D. Clarke,
Anthony P. Whitworth,
David A. Hubber
Abstract:
We use smoothed particle hydrodynamic simulations to investigate the growth of perturbations in infinitely long, initially sub-critical but accreting filaments. The growth of these perturbations leads to filament fragmentation and the formation of cores. Most previous work on this subject has been confined to the growth and fragmentation of equilibrium filaments and has found that there exists a p…
▽ More
We use smoothed particle hydrodynamic simulations to investigate the growth of perturbations in infinitely long, initially sub-critical but accreting filaments. The growth of these perturbations leads to filament fragmentation and the formation of cores. Most previous work on this subject has been confined to the growth and fragmentation of equilibrium filaments and has found that there exists a preferential fragmentation length scale which is roughly 4 times the filament's diameter. Our results show a more complicated dispersion relation with a series of peaks linking perturbation wavelength and growth rate. These are due to gravo-acoustic oscillations along the longitudinal axis during the sub-critical phase of growth. The positions of the peaks in growth rate have a strong dependence on both the mass accretion rate onto the filament and the temperature of the gas. When seeded with a multi-wavelength density power spectrum there exists a clear preferred core separation equal to the largest peak in the dispersion relation. Our results allow one to estimate a minimum age for a filament which is breaking up into regularly spaced fragments, as well as a maximum accretion rate. We apply the model to observations of filaments in Taurus by Tafalla & Hacar (2015) and find accretion rates consistent with those estimated by Palmeirim et al. (2013).
△ Less
Submitted 25 February, 2016; v1 submitted 20 January, 2016;
originally announced February 2016.
-
Investigating the Global Collapse of Filaments Using Smoothed Particle Hydrodynamics
Authors:
Seamus D. Clarke,
Anthony P. Whitworth
Abstract:
We use Smoothed Particle Hydrodynamic simulations of cold, uniform density, self-gravitating filaments, to investigate their longitudinal collapse timescales; these timescales are important because they determine the time available for a filament to fragment into cores. A filament is initially characterised by its line-mass, $μ$, its radius, $R$ (or equivalently its density $ρ\!=\!μ/πR^2$), and it…
▽ More
We use Smoothed Particle Hydrodynamic simulations of cold, uniform density, self-gravitating filaments, to investigate their longitudinal collapse timescales; these timescales are important because they determine the time available for a filament to fragment into cores. A filament is initially characterised by its line-mass, $μ$, its radius, $R$ (or equivalently its density $ρ\!=\!μ/πR^2$), and its aspect ratio, $A\;\,(\equiv Z/R$, where $Z$ is its half-length). The gas is only allowed to contract longitudinally, i.e. parallel to the symmetry axis of the filament (the $z$-axis). Pon et al. (2012) have considered the global dynamics of such filaments analytically. They conclude that short filaments ($A\! < \!5$) collapse along the $z$-axis more-or-less homologously, on a time-scale $t_{_{\rm HOM}} \sim 0.44\,A\,(Gρ)^{-1/2}$; in contrast, longer filaments ($A\! > \!5$) undergo end-dominated collapse, i.e. two dense clumps form at the ends of the filament and converge on the centre sweeping up mass as they go, on a time-scale $t_{_{\rm END}} \sim 0.98\,A^{1/2}\,(Gρ)^{-1/2}$. Our simulations do not corroborate these predictions. First, for all $A\! > \!2$, the collapse time satisfies a single equation \[t_{_{\rm COL}}\;\sim\;(0.49+0.26A)(Gρ)^{-1/2}\,,\] which for large $A$ is much longer than the Pon et al. prediction. Second, for all $A\! > \!2$, the collapse is end-dominated. Third, before being swept up, the gas immediately ahead of an end-clump is actually accelerated outwards by the gravitational attraction of the approaching clump, resulting in a significant ram pressure. For high aspect ratio filaments the end-clumps approach an asymptotic inward speed, due to the fact that they are doing work both accelerating and compressing the gas they sweep up. Pon et al. appear to have neglected the outward acceleration and its consequences.
△ Less
Submitted 26 February, 2015;
originally announced February 2015.