Showing 1–50 of 61 results for author: Falle, S

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.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

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)

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Shocking interactions of supernova remnants with atomic and molecular clouds -- the interplay between shocks, thermal instability and gravity in the large cloud regime

Authors: M. M. Kupilas, J. M. Pittard, C. J. Wareing, S. A. E. G. Falle

Abstract: Using the adaptive mesh refinement code MG, we perform 3D hydrodynamic simulations of a supernova-cloud interaction in the "large cloud regime". The cloud is initially atomic and evolving due to the thermal instability (TI) and gravity. We study interactions in a "pre-TI" and "post-TI" stage when cold and dense clumps are present, and compare these results to idealised shock-cloud scenarios in the… ▽ More Using the adaptive mesh refinement code MG, we perform 3D hydrodynamic simulations of a supernova-cloud interaction in the "large cloud regime". The cloud is initially atomic and evolving due to the thermal instability (TI) and gravity. We study interactions in a "pre-TI" and "post-TI" stage when cold and dense clumps are present, and compare these results to idealised shock-cloud scenarios in the "small cloud regime", and a scenario without shocks. On aggregate, the supernova disruption is significantly weaker than that from an idealised shock due to the supernova impact being instantaneous, and not continuous. In both supernova-cloud interactions, we observe two shocks impact the cloud, followed by the development of a weak 10 km s$^{-1}$ upstream flow on the cloud interface, and a global ambient pressure drop. When the cloud is still atomic, it expands due to this drop. Additionally, the TI is triggered at the front of the cloud, causing the formation of a cap-like structure with clumps embedded inside. The upstream flow converges in this region, resulting in a lobe-like cloud morphology. When the cloud is molecular, the transmitted shock disrupts the inter-clump material and causes the clumps' outer envelopes to expand slightly and form tail-like morphologies. These effects are less pronounced than those in our shock-cloud scenarios, and more pronounced that those in our un-shocked scenario. After 3.5 Myrs, the effects from the supernova decay and the cloud returns to an almost indistinguishable state from an un-shocked cloud, in spite of the global ambient pressure drop. In neither supernova-cloud scenario do we see any local gravitational collapse. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 15 pages, 9 figures

Journal ref: MNRAS 513, 2022, 3345

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filament… ▽ More We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

Revealing the diverse magnetic field morphologies in Taurus dense cores with sensitive sub-millimeter polarimetry

Authors: Chakali Eswaraiah, Di Li, Ray S. Furuya, Tetsuo Hasegawa, Derek Ward-Thompson, Keping Qiu, Nagayoshi Ohashi, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, David Berry, Yasuo Doi, Tao-Chung Ching, Shih-Ping Lai, Jia-Wei Wang, Patrick M. Koch, Jungmi Kwon, Woojin Kwon, Pierre Bastien, Doris Arzoumanian, Simon Coudé, Archana Soam, Lapo Fanciullo, Hsi-Wei Yen, Junhao Liu , et al. (120 additional authors not shown)

Abstract: We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protost… ▽ More We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $μ$G, 44$\pm$16 $μ$G, and 12$\pm$5 $μ$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 20 pages (10 are main), 5 figures (3 are main), and 3 tables (2 are main); Accepted for its publication in ApJL

Dust polarized emission observations of NGC 6334; BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network

Authors: D. Arzoumanian, R. Furuya, T. Hasegawa, M. Tahani, S. Sadavoy, C. L. H. Hull, D. Johnstone, P. M. Koch, S. -i. Inutsuka, Y. Doi, T. Hoang, T. Onaka, K. Iwasaki, Y. Shimajiri, T. Inoue, N. Peretto, P. André, P. Bastien, D. Berry, H. -R. V. Chen, J. Di Francesco, C. Eswaraiah, L. Fanciullo, L. M. Fissel, J. Hwang , et al. (123 additional authors not shown)

Abstract: [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We presen… ▽ More [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$μ$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($θ_{B}$). We derive the power spectrum of the intensity and $θ_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $θ_{B}$ varies coherently along the filaments. The observed power spectrum of $θ_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $θ_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $θ_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A78 (2021)

Interactions of a shock with a molecular cloud at various stages of its evolution due to thermal instability and gravity

Authors: M. M. Kupilas, C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: Using the adaptive mesh refinement code MG, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of thermal in… ▽ More Using the adaptive mesh refinement code MG, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of thermal instability and gravity. The second case introduces a shock whilst the cloud is still in the warm atomic phase, and the third scenario introduces a shock once the molecular cloud has formed. The shocks accelerate the global collapse of the clouds with both experiencing local gravitational collapse prior to this. When the cloud is still atomic, the evolution is shock dominated and structures form due to dynamical instabilities within a radiatively cooled shell. While the transmitted shock can potentially trigger the thermal instability, this is prevented as material is shocked multiple times on the order of a cloud crushing time-scale. When the cloud is molecular, the post-shock flow is directed via the pre-existing structure through low-density regions in the inter-clump medium. The clumps are accelerated and deformed as the flow induces clump-clump collisions and mergers that collapse under gravity. For a limited period, both shocked cases show a mixture of Kolmogorov and Burgers turbulence-like velocity and logarithmic density power spectra, and strongly varying density spectra. The clouds presented in this work provide realistic conditions that will be used in future feedback studies. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Observations of magnetic fields surrounding LkH$α$ 101 taken by the BISTRO survey with JCMT-POL-2

Authors: Nguyen Bich Ngoc, Pham Ngoc Diep, Harriet Parsons, Kate Pattle, Thiem Hoang, Derek Ward-Thompson, Le Ngoc Tram, Charles L. H. Hull, Mehrnoosh Tahani, Ray Furuya, Pierre Bastien, Keping Qiu, Tetsuo Hasegawa, Woojin Kwon, Yasuo Doi, Shih-Ping Lai, Simon Coude, David Berry, Tao-Chung Ching, Jihye Hwang, Archana Soam, Jia-Wei Wang, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun , et al. (124 additional authors not shown)

Abstract: We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to… ▽ More We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrmμ$G. The mass-to-magnetic-flux ratio in critical value units, $λ\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$α$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $δB/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $α=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$α$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

The JCMT BISTRO survey: alignment between outflows and magnetic fields in dense cores/clumps

Authors: Hsi-Wei Yen, Patrick M. Koch, Charles L. H. Hull, Derek Ward-Thompson, Pierre Bastien, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Tao-Chung Ching, Eun Jung Chung, Simon Coude, James Di Francesco, Pham Ngoc Diep, Yasuo Doi, Chakali Eswaraiah, Sam Falle, Gary Fuller, Ray S. Furuya, Ilseung Han, Jennifer Hatchell, Martin Houde, Shu-ichiro Inutsuka, Doug Johnstone, Ji-hyun Kang , et al. (21 additional authors not shown)

Abstract: We compare the directions of molecular outflows of 62 low-mass Class 0 and I protostars in nearby (<450 pc) star-forming regions with the mean orientations of the magnetic fields on 0.05-0.5 pc scales in the dense cores/clumps where they are embedded. The magnetic field orientations were measured using the JCMT POL-2 data taken by the BISTRO-1 survey and from the archive. The outflow directions we… ▽ More We compare the directions of molecular outflows of 62 low-mass Class 0 and I protostars in nearby (<450 pc) star-forming regions with the mean orientations of the magnetic fields on 0.05-0.5 pc scales in the dense cores/clumps where they are embedded. The magnetic field orientations were measured using the JCMT POL-2 data taken by the BISTRO-1 survey and from the archive. The outflow directions were observed with interferometers in the literature. The observed distribution of the angles between the outflows and the magnetic fields peaks between 15 and 35 degrees. After considering projection effects, our results could suggest that the outflows tend to be misaligned with the magnetic fields by 50+/-15 degrees in three-dimensional space and are less likely (but not ruled out) randomly oriented with respect to the magnetic fields. There is no correlation between the misalignment and the bolometric temperatures in our sample. In several sources, the small-scale (1000-3000 au) magnetic fields is more misaligned with the outflows than their large-scale magnetic fields, suggesting that the small-scale magnetic field has been twisted by the dynamics. In comparison with turbulent MHD simulations of core formation, our observational results are more consistent with models in which the energy densities in the magnetic field and the turbulence of the gas are comparable. Our results also suggest that the misalignment alone cannot sufficiently reduce the efficiency of magnetic braking to enable formation of the observed number of large Keplerian disks with sizes larger than 30-50 au. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 27 pages, 14 figures, Accepted by ApJ

Striations, integrals, hourglasses and collapse -- thermal instability driven magnetic simulations of molecular clouds

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: The MHD version of the adaptive mesh refinement (AMR) code, MG, has been employed to study the interaction of thermal instability, magnetic fields and gravity through 3D simulations of the formation of collapsing cold clumps on the scale of a few parsecs, inside a larger molecular cloud. The diffuse atomic initial condition consists of a stationary, thermally unstable, spherical cloud in pressure… ▽ More The MHD version of the adaptive mesh refinement (AMR) code, MG, has been employed to study the interaction of thermal instability, magnetic fields and gravity through 3D simulations of the formation of collapsing cold clumps on the scale of a few parsecs, inside a larger molecular cloud. The diffuse atomic initial condition consists of a stationary, thermally unstable, spherical cloud in pressure equilibrium with lower density surroundings and threaded by a uniform magnetic field. This cloud was seeded with 10% density perturbations at the finest initial grid level around n=1.1 cm^{-3} and evolved with self-gravity included from the outset. Several cloud diameters were considered (100 pc, 200 pc and 400 pc) equating to several cloud masses (17,000 Msun, 136,000 Msun and 1.1x10^6 Msun). Low-density magnetic-field-aligned striations were observed as the clouds collapse along the field lines into disc-like structures. The induced flow along field lines leads to oscillations of the sheet about the gravitational minimum and an integral-shaped appearance. When magnetically supercritical, the clouds then collapse and generate hourglass magnetic field configurations with strongly intensified magnetic fields, reproducing observational behaviour. Resimulation of a region of the highest mass cloud at higher resolution forms gravitationally-bound collapsing clumps within the sheet that contain clump-frame supersonic (M~5) and super-Alfvenic (M_A~4) velocities. Observationally realistic density and velocity power spectra of the cloud and densest clump are obtained. Future work will use these realistic initial conditions to study individual star and cluster feedback. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 20 pages, 11 figures + appendix and figure. Raw data: doi.org/10.5518/897

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

Thermal instability revisited

Authors: Samuel Falle, Christopher Wareing, Julian Pittard

Abstract: Field's linear analysis of thermal instability is repeated using methods related to Whitham's theory of wave hierarchies, which brings out the physically relevant parameters in a much clearer way than in the original analysis. It is also used for the stability of non-equilibrium states and we show that for gas cooling behind a shock, the usual analysis is only quantitatively valid for shocks that… ▽ More Field's linear analysis of thermal instability is repeated using methods related to Whitham's theory of wave hierarchies, which brings out the physically relevant parameters in a much clearer way than in the original analysis. It is also used for the stability of non-equilibrium states and we show that for gas cooling behind a shock, the usual analysis is only quantitatively valid for shocks that are just able to trigger a transition to the cold phase. A magnetic field can readily be included and we show that this does not change the stability criteria. By considering steady shock solutions, we show that almost all plausible initial conditions lead to a magnetically dominated state on the unstable part of the equilibrium curve. These results are used to analyse numerical calculations of perturbed steady shock solutions and of shocks interacting with a warm cloud. △ Less

Submitted 4 January, 2020; originally announced January 2020.

Comments: 17 pages, 15 figures

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

Sheets, filaments and clumps - high resolution simulations of how the thermal instability can form molecular clouds

Authors: C. J. Wareing, S. A. E. G. Falle, J. M. Pittard

Abstract: This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.… ▽ More This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.1cm^{-3} and evolved with self-gravity included. No magnetic field was imposed. Resimulations at a higher resolution of a region extracted from this simulation (down to 0.039pc), show that the thermal instability forms sheets, then filaments and finally clumps. The width of the filaments increases over time, in one particular case from 0.26 to 0.56pc. Thereafter clumps with sizes of around 5pc grow at the intersections of filaments. 21 distinct clumps, with properties similar to those observed in molecular clouds, are found by using the FellWalker algorithm to find minima in the gravitational potential. Not all of these are gravitationally bound, but the convergent nature of the flow and increasing central density suggest they are likely to form stars. Further simulation of the most massive clump shows the gravitational collapse to a density >10^6 cm^{-3}. These results provide realistic initial conditions that can be used to study feedback in individual clumps, interacting clumps and the entire molecular cloud complex. △ Less

Submitted 12 March, 2019; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: 19 pages, 10 figures. MNRAS accepted. Moderate revision and extension

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

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

A new mechanical stellar wind feedback model for the Rosette Nebula

Authors: C. J. Wareing, J. M. Pittard, N. J. Wright, S. A. E. G. Falle

Abstract: The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data r… ▽ More The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD46223, from NGC 2244 occurred 1.73 (+0.34,-0.25)Myr (1$σ$ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD46150, the age is set for the famous ionised region at more than ten times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40M$_{\odot}$ star in a thin sheet-like molecular cloud. We form the 135,000M$_{\odot}$ cloud from thermally-unstable diffuse interstellar medium under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5Myr, including cavity size, stellar age, magnetic field and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: This is the Version of Record, published RAS Open Access CC BY online in MNRAS, 2018 Feb 13. 15 Pages, 9 figures

Journal ref: MNRAS, Volume 475, Issue 3, Pages 3598-3612, 2018

Hydrodynamic simulations of mechanical stellar feedback in a molecular cloud formed by thermal instability

Authors: Christopher J. Wareing, Julian M. Pittard, Samuel A. E. G. Falle

Abstract: We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with… ▽ More We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with randomly distributed high density condensations. The stellar winds are introduced using appropriate non-rotating Geneva stellar evolution models. In the 15 Msun star case, the wind has very little effect, spreading around a few neighbouring clumps before becoming overwhelmed by the cloud collapse. In contrast, in the 40 Msun, 60 Msun and 120 Msun star cases, the more powerful stellar winds create large cavities and carve channels through the cloud, breaking out into the surrounding tenuous medium during the wind phase and considerably altering the cloud structure. After 4.97 Myrs, 3.97 Myrs and 3.01 Myrs respectively, the massive stars explode as supernovae (SNe). The wind-sculpted surroundings considerably affect the evolution of these SN events as they both escape the cloud along wind-carved channels and sweep up remaining clumps of cloud/wind material. The `cloud' as a coherent structure does not survive the SN from any of these stars, but only in the 120 Msun case is the cold molecular material completely destabilised and returned to the unstable thermal phase. In the 40 Msun and 60 Msun cases, coherent clumps of cold material are ejected from the cloud by the SN, potentially capable of further star formation. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 33 pages, 26 figures, accepted to MNRAS 2017 June 6. Raw data available from http://doi.org/10.5518/201

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

The interaction of hydrodynamic shocks with self-gravitating clouds

Authors: S. A. E. G. Falle, B. Vaidya, T. W. Hartquist

Abstract: We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and no… ▽ More We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and not in others and we derive a simple estimate for the shock strength required for this to occur. These results are relevant to dense cores and Bok globules in star forming regions subjected to shocks produced by stellar feedback. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 9 pages 9 figures

Ambipolar diffusion regulated collapse of filaments threaded by perpendicular magnetic fields

Authors: C. A. Burge, S. Van Loo, S. A. E. G. Falle, T. W. Hartquist

Abstract: We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic… ▽ More We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic flux is redistributed within the filament from the centre towards the envelope. The rate of the collapse is inversely proportional to the fractional ionisation and two gravitationally-driven ambipolar diffusion regimes for the collapse are observed as predicted in Mouschovias & Morton (1991). For high values of the ionisation coefficient, that is $X \geq 10^{-7}$, the gas is strongly coupled to the magnetic field and the Jeans length is larger than the ambipolar diffusion length scale. Then the collapse is governed by magnetically-regulated ambipolar diffusion. For $X \lesssim 10^{-8}$, the gas is weakly coupled to the magnetic field and the magnetic support is removed by gravitationally-dominated ambipolar diffusion. Here, neutrals and ions only collide sporadically, that is the ambipolar diffusion length scale is larger than the Jeans length. When decaying turbulence is included, additional support is provided to the filament. This slows down the collapse of the filament even in the absence of a magnetic field. When a magnetic field is present, the collapse rate increases by a ratio smaller than for the non-magnetic case. This is because of a speed-up of the ambipolar diffusion due to larger magnetic field gradients generated by the turbulence and because the ambipolar diffusion aids the dissipation of turbulence below the ambipolar diffusion length scale. The highest increase in the rate is observed for the lowest ionisation coefficient and the highest turbulent intensity. △ Less

Submitted 22 September, 2016; originally announced September 2016.

Comments: 10 pages, 11 figures, accepted for publication by A&A

Journal ref: A&A 596, A28 (2016)

Magnetohydrodynamic simulations of mechanical stellar feedback in a sheet-like molecular cloud

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morpho… ▽ More We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morphology that in projection appears filamentary. The stellar winds are introduced using appropriate Geneva stellar evolution models. In the 15 solar mass star case, the wind forms a narrow bipolar cavity with minimal effect on the parent cloud. In the 40 solar mass star case, the more powerful stellar wind creates a large cylindrical cavity through the centre of the cloud. After 12.5 Myrs and 4.97 Myrs respectively, the massive stars explode as supernovae (SNe). In the 15 solar mass star case, the SN material and energy is primarily deposited into the molecular cloud surroundings over ~10^5 years before the SN remnant escapes the cloud. In the 40 solar mass star case, a significant fraction of the SN material and energy rapidly escapes the molecular cloud along the wind cavity in a few tens of kiloyears. Both SN events compress the molecular cloud material around them to higher densities (so may trigger further star formation), and strengthen the magnetic field, typically by factors of 2-3 but up to a factor of 10. Our simulations are relevant to observations of bubbles in flattened ring-like molecular clouds and bipolar HII regions. △ Less

Submitted 16 November, 2016; v1 submitted 16 May, 2016; originally announced May 2016.

Comments: 29 pages, 20 figures, MNRAS accepted. Additions after review include an extended analysis in a new Section 4, clarification of numerous minor points throughout and quantitative data figures in a new Appendix

MHD simulation of the formation of clumps and filaments in quiescent diffuse medium by thermal instability

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle, S. Van Loo

Abstract: We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range o… ▽ More We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure. △ Less

Submitted 8 March, 2016; v1 submitted 17 January, 2016; originally announced January 2016.

Comments: 17 pages, 10 figures. Accepted to MNRAS, minor changes and additions to text, primarily in Section 6

Magnetic Fields and Galactic Star Formation Rates

Authors: Sven Van Loo, Jonathan C. Tan, Sam A. E. G. Falle

Abstract: The regulation of galactic-scale star formation rates (SFRs) is a basic problem for theories of galaxy formation and evolution: which processes are responsible for making observed star formation rates so inefficient compared to maximal rates of gas content divided by dynamical timescale? Here we study the effect of magnetic fields of different strengths on the evolution of giant molecular clouds (… ▽ More The regulation of galactic-scale star formation rates (SFRs) is a basic problem for theories of galaxy formation and evolution: which processes are responsible for making observed star formation rates so inefficient compared to maximal rates of gas content divided by dynamical timescale? Here we study the effect of magnetic fields of different strengths on the evolution of giant molecular clouds (GMCs) within a kiloparsec patch of a disk galaxy and resolving scales down to $\simeq0.5\:{\rm{pc}}$. Including an empirically motivated prescription for star formation from dense gas ($n_{\rm{H}}>10^5\:{\rm{cm}^{-3}}$) at an efficiency of 2\% per local free-fall time, we derive the amount of suppression of star formation by magnetic fields compared to the nonmagnetized case. We find GMC fragmentation, dense clump formation and SFR can be significantly affected by the inclusion of magnetic fields, especially in our strongest investigated $B$-field case of $80\:{\rmμ}$G. However, our chosen kpc-scale region, extracted from a global galaxy simulation, happens to contain a starbursting cloud complex that is only modestly affected by these magnetic fields and likely requires internal star formation feedback to regulate its SFR. △ Less

Submitted 25 January, 2015; v1 submitted 27 November, 2014; originally announced November 2014.

Comments: 7 pages, 5 figures. Accepted by ApJ Letters

Numerical Simulations of a Shock Interacting with Multiple Magnetized Clouds

Authors: R. Alūzas, J. M. Pittard, S. A. E. G. Falle, T. W. Hartquist

Abstract: We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour o… ▽ More We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $β$ in cloud material can vary rapidly as clouds collide with one another, but low values of $β$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds. △ Less

Submitted 24 July, 2014; originally announced July 2014.

Comments: 24 pages, 30 figures. Accepted for publication in MNRAS

Ram pressure stripping of the hot gaseous haloes of galaxies using the k-epsilon sub-grid turbulence model

Authors: J. L. Close, J. M. Pittard, T. W. Hartquist, S. A. E. G. Falle

Abstract: We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stri… ▽ More We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: 10 pages, 8 figures, 3 tables; accepted for publication in MNRAS

The Responses of Magnetically Sub-Critical Cores to Shocks

Authors: B. Vaidya, T. W. Hartquist, S. A. E. G. Falle

Abstract: An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even fo… ▽ More An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even for relatively weak shocks. Self-gravity plays only a minor role in determining the highest density that is reached, but it does play a role in the subsequent evolution. The densities at comparable times differ by a factor of a few for shocks initially propagating perpendicularly or obliquely to the magnetic field in the ambient medium. △ Less

Submitted 8 May, 2013; originally announced May 2013.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS

Convergence of AMR and SPH simulations - I. Hydrodynamical resolution and convergence tests

Authors: D. A. Hubber, S. A. E. G. Falle, S. P. Goodwin

Abstract: We compare the results for a set of hydrodynamical tests performed with the AMR finite volume code, MG and the SPH code, SEREN. The test suite includes shock tube tests, with and without cooling, the non-linear thin-shell instability and the Kelvin-Helmholtz instability. The main conclusions are : (i) the two methods converge in the limit of high resolution and accuracy in most cases. All tests sh… ▽ More We compare the results for a set of hydrodynamical tests performed with the AMR finite volume code, MG and the SPH code, SEREN. The test suite includes shock tube tests, with and without cooling, the non-linear thin-shell instability and the Kelvin-Helmholtz instability. The main conclusions are : (i) the two methods converge in the limit of high resolution and accuracy in most cases. All tests show good agreement when numerical effects (e.g. discontinuities in SPH) are properly treated. (ii) Both methods can capture adiabatic shocks and well-resolved cooling shocks perfectly well with standard prescriptions. However, they both have problems when dealing with under-resolved cooling shocks, or strictly isothermal shocks, at high Mach numbers. The finite volume code only works well at 1st order and even then requires some additional artificial viscosity. SPH requires either a larger value of the artificial viscosity parameter, alpha_AV, or a modified form of the standard artificial viscosity term using the harmonic mean of the density, rather than the arithmetic mean. (iii) Some SPH simulations require larger kernels to increase neighbour number and reduce particle noise in order to achieve agreement with finite volume simulations. However, this is partly due to the need to reduce noise that can corrupt the growth of small-scale perturbations. In contrast, instabilities seeded from large-scale perturbations do not require more neighbours and hence work well with standard SPH formulations and converge with the finite volume simulations. (iv) For purely hydrodynamical problems, SPH simulations take an order of magnitude longer to run than finite volume simulations when running at equivalent resolutions, i.e. when they both resolve the underlying physics to the same degree. This requires about 2-3 times as many particles as the number of cells. △ Less

Submitted 11 April, 2013; originally announced April 2013.

Comments: 19 pages, 9 figures, 1 table, accepted for publication in MNRAS

Kiloparsec-Scale Simulations of Magnetised Molecular Clouds in Disk Galaxies

Authors: Sven Van Loo, Michael J. Butler, Jonathan C. Tan, Sam A. E. G. Falle

Abstract: We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt r… ▽ More We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. The evolution of these simulated GMCs and clumps is determined by self-gravity balanced by turbulent pressure support and the large scale galactic shear. While the cloud structures and densities significantly change during their evolution, they remain roughly in virial equilibrium for time scales exceeding the free-fall time of GMCs, indicating that energy from the galactic shear continuously cascades down. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but this yields global star formation rates that are more than ~two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to the over-production of dense clump gas. To explain this discrepancy, we anticipate magnetic fields to provide additional support. Low-resolution simulations indeed show that the magnetic field reduces the star formation rate. △ Less

Submitted 28 November, 2012; v1 submitted 2 November, 2012; originally announced November 2012.

Comments: For proceedings of ASTRONUM 2012, 6 pages, 5 figures

Sputtering in oblique C-type shocks

Authors: S. Van Loo, I. Ashmore, P. Caselli, S. A. E. G. Falle, T. W. Hartquist

Abstract: We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering comme… ▽ More We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced. △ Less

Submitted 18 September, 2012; originally announced September 2012.

Comments: 8 pages, 4 figures and 3 tables. Accepted by MNRAS

Magnetic Fields and Star Formation

Authors: S. Van Loo, T. W. Hartquist, S. A. E. G. Falle

Abstract: Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problem… ▽ More Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problems that have received attention in the last few years. The coverage is not comprehensive, and the authors have drawn on their own results more and touched more briefly on those of others than they would in a normal review. Theirs is a personal contribution to the issue honouring Leon, one of the truly great gentlemen, wits, and most insightful of astrophysicists. △ Less

Submitted 31 July, 2012; originally announced July 2012.

Comments: 13 pages, no figures, to be published in Astronomy and Geophysics

Numerical simulations of shocks encountering clumpy regions

Authors: R. Alūzas, J. M. Pittard, T. W. Hartquist, S. A. E. G. Falle, R. Langton

Abstract: We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud m… ▽ More We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution. △ Less

Submitted 26 June, 2012; originally announced June 2012.

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

Journal ref: Monthly Notices of the Royal Astronomical Society 425 (2012) 2212-2227

Cosmic Rays and Radiative Instabilities

Authors: T. W. Hartquist, A. Y. Wagner, S. A. E. G. Falle, J. M. Pittard, S. Van Loo

Abstract: In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the th… ▽ More In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent. △ Less

Submitted 5 August, 2011; originally announced August 2011.

Comments: 5 pages, to appear in the Proceedings of Cosmic Rays & Their Interstellar Medium Environment CRISM-2011 (Montpellier, 26th June - 1st July 2011), to be published in an issue of MEMORIE della Società Astronomica Italiana, eds. A. Marcowith, A.Bykov, K.Ferriere, Y.Gallant, P.Hennebelle, T.Montmerle

Shock-triggered formation of magnetically-dominated clouds. II. Weak shock-cloud interaction in three dimensions

Authors: S. Van Loo, S. A. E. G. Falle, T. W. Hartquist

Abstract: To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the sh… ▽ More To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the shocked cloud with a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 3D, we investigate the effect of different magnetic field orientations including parallel, perpendicular and oblique to the shock normal. We find that the angle between the shock normal and the magnetic field must be small to produce clouds with properties similar to observed molecular clouds. △ Less

Submitted 30 March, 2010; originally announced March 2010.

Comments: 10 pages, 12 figures, accepted for publication in MNRAS

The Turbulent Destruction of Clouds - II. Mach Number Dependence, Mass-loss Rates, and Tail Formation

Authors: J. M. Pittard, T. W. Hartquist, S. A. E. G. Falle

Abstract: The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagn… ▽ More The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3. △ Less

Submitted 15 February, 2010; v1 submitted 10 February, 2010; originally announced February 2010.

Comments: 19 pages, 19 figures, accepted by MNRAS (extra figures after end of manuscript removed)

Transient evolution of C-type shocks in dusty regions of varying density

Authors: I. Ashmore, S. Van Loo, P. Caselli, S. A. E. G. Falle, T. W. Hartquist

Abstract: Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-… ▽ More Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-type shocks interacting with density perturbations. We include a self-consistent calculation of the thermal and ionisation balances and a fluid treatment of grains. We identify features that develop when a multifluid shock encounters a density inhomogeneity to investigate whether any part of the precursor region ever behaves in a quasi-steady fashion. If it does the shock may be modelled approximately without solving the time-dependent hydromagnetic equations. Simulations were made for initially steady oblique C-type shocks encountering density inhomogeneities. For a semi-finite inhomogeneity with a density larger than the surrounding medium, a transmitted shock evolves from being J-type to a steady C-type shock on a timescale comparable to the ion-flow time through it. A sufficiently upstream part of the precursor of an evolving J-type shock is quasi-steady. The ion-flow timescale is also relevant for the evolution of a shock moving into a region of decreasing density. The models for shocks propagating into regions in which the density increases and then decreases to its initial value cannot be entirely described in terms of the results obtained for monotonically increasing and decreasing densities. For the latter model, the long-term evolution to a C-type shock cannot be approximated by quasi-steady models. △ Less

Submitted 25 November, 2009; originally announced November 2009.

Comments: 11 pages, 9 figures

Time-dependent simulations of steady C-type shocks

Authors: S. Van Loo, I. Ashmore, P. Caselli, S. A. E. G. Falle, T. W. Hartquist

Abstract: Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging and studied the effect of single-sized and multiple sized grains on the shock structure. Our models… ▽ More Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging and studied the effect of single-sized and multiple sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionisation structures including appropriate microphysics. At low densities the grains do not play any significant role in the shock dynamics. At high densities, the ionisation fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allows a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena. △ Less

Submitted 16 January, 2009; originally announced January 2009.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

How Interstellar Chemistry (and Astrochemistry More Generally) Became Useful

Authors: T. W. Hartquist, S. Van Loo, S. A. E. G. Falle

Abstract: In 1986 Alex Dalgarno published a paper entitled "Is Interstellar Chemistry Useful?" By the middle 1970s, and perhaps even earlier, Alex had hoped that astronomical molecules would prove to: possess significant diagnostic utility; control many of the environments in which they exist; stimulate a wide variety of physicists and chemists who are at least as fascinated by the mechanisms forming and… ▽ More In 1986 Alex Dalgarno published a paper entitled "Is Interstellar Chemistry Useful?" By the middle 1970s, and perhaps even earlier, Alex had hoped that astronomical molecules would prove to: possess significant diagnostic utility; control many of the environments in which they exist; stimulate a wide variety of physicists and chemists who are at least as fascinated by the mechanisms forming and removing the molecules as by astronomy. His own research efforts have contributed greatly to the realization of that hope. This paper contains a few examples of: how molecules are used to diagnose large-scale dynamics in astronomical sources including star forming regions and supernovae; the ways in which molecular processes control the evolution of astronomical objects such as dense cores destined to become stars and very evolved giant stars; theoretical and laboratory investigations that elucidate the processes producing and removing astronomical molecules and allow their detection. △ Less

Submitted 17 December, 2008; originally announced December 2008.

Comments: 15 pages, 1 figure, to be published in Proceedings of the Dalgarno Celebratory Symposium:Contributions to Atomic, Molecular, and Optical Physics, Astrophysics and Atmospheric Physics, Eds. James Babb, Kate Kirby and Hossein Sadeghpour

A cosmic-ray precursor model for a Balmer-dominated shock in Tycho's supernova remnant

Authors: A. Y. Wagner, J. -J. Lee, J. C. Raymond, T. W. Hartquist, S. A. E. G. Falle

Abstract: We present a time-dependent cosmic-ray modified shock model for which the calculated H-alpha emissivity profile agrees well with the H-alpha flux increase ahead of the Balmer-dominated shock at knot g in Tycho's supernova remnant, observed by Lee et al (2007). The backreaction of the cosmic ray component on the thermal component is treated in the two-fluid approximation, and we include thermal p… ▽ More We present a time-dependent cosmic-ray modified shock model for which the calculated H-alpha emissivity profile agrees well with the H-alpha flux increase ahead of the Balmer-dominated shock at knot g in Tycho's supernova remnant, observed by Lee et al (2007). The backreaction of the cosmic ray component on the thermal component is treated in the two-fluid approximation, and we include thermal particle injection and energy transfer due to the acoustic instability in the precursor. The transient state of our model that describes the current state of the shock at knot g, occurs during the evolution from a thermal gas dominated shock to a smooth cosmic-ray dominated shock. Assuming a distance of 2.3 kpc to Tycho's remnant we obtain values for the cosmic ray diffusion coefficient, the injection parameter, and the time scale for the energy transfer of 10^{24} cm^{2} s^{-1}, 4.2x10^{-3}, and 426 y, respectively. We have also studied the parameter space for fast (300 km s^{-1} - 3000 km s^{-1}), time-asymptotically steady shocks and have identified a branch of solutions, for which the temperature in the cosmic ray precursor typically reaches 2-6x10^{4} K and the bulk acceleration of the flow through the precursor is less than 10 km s^{-1}. These solutions fall into the low cosmic ray acceleration efficiency regime and are relatively insensitive to shock parameters. This low cosmic ray acceleration efficiency branch of solutions may provide a natural explanation for the line broadening of the H-alpha narrow component observed in non-radiative shocks in many supernova remnants. △ Less

Submitted 15 September, 2008; originally announced September 2008.

Comments: 13 pages, 5 figures, to be published in ApJ

Journal ref: Astrophys.J.690:1412-1423,2009

The turbulent destruction of clouds - I. A k-epsilon treatment of turbulence in 2D models of adiabatic shock-cloud interactions

Authors: J. M. Pittard, S. A. E. G. Falle, T. W. Hartquist, J. E. Dyson

Abstract: The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous… ▽ More The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous single-cloud studies. In reality, the flow behind the shock is also likely to be turbulent, with non-uniform density, pressure and velocity structure created as the shock sweeps over inhomogenities upstream of the cloud. To address these twin issues we use a sub-grid compressible k-epsilon turbulence model to estimate the properties of the turbulence generated in shock-cloud interactions and the resulting increase in the transport coefficients that the turbulence brings. A detailed comparison with the output from an inviscid hydrodynamical code puts these new results into context. We find that cloud destruction in inviscid and k-epsilon models occurs at roughly the same speed when the post-shock flow is smooth and when the density contrast between the cloud and inter-cloud medium is less than 100. However, there are increasing and significant differences as this contrast increases. Clouds subjected to strong ``buffeting'' by a highly turbulent post-shock environment are destroyed significantly quicker. Additional calculations with an inviscid code where the post-shock flow is given random, grid-scale, motions confirms the more rapid destruction of the cloud. Our results clearly show that turbulence plays an important role in shock-cloud interactions, and that environmental turbulence adds a new dimension to the parameter space which has hitherto been studied (abridged). △ Less

Submitted 28 July, 2008; originally announced July 2008.

Comments: 31 pages, 22 figures, accepted for publication in MNRAS

Dust in Interstellar Clouds, Evolved Stars and Supernovae

Authors: T. W. Hartquist, S. Van Loo, S. A. E. G. Falle, P. Caselli, I. Ashmore

Abstract: Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to… ▽ More Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest. △ Less

Submitted 9 July, 2008; originally announced July 2008.

Comments: 4 pages, to be published in the proceedings of Fifth International Conference on Physics of Dusty Plasmas

Journal ref: AIP Conf.Proc.1041:101-104,2008

The effect of ambipolar resistivity on the formation of dense cores

Authors: S. Van Loo, S. A. E. G. Falle, T. W. Hartquist, A. J. Barker

Abstract: We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar diffusion then plays an important role in the formation process. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D nume… ▽ More We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar diffusion then plays an important role in the formation process. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D numerical simulations, we follow the evolution of the fast-mode wave. The simulations are done with a multifluid, adaptive mesh refinement MHD code. Initial perturbations with wavelengths that are 2 orders of magnitude larger than the dissipation length are strongly affected by the ion-neutral drift. Only in situations where there are large variations in the magnetic field corresponding to a highly turbulent gas can fast-mode waves generate dense cores. This means that, in most cores, no substructure can be produced. However, Core D of TMC-1 is an exception to this case. Due to its atypically high ionisation fraction, waves with wavelengths up to 3 orders of magnitude greater than the dissipation length can be present. Such waves are only weakly affected by ambipolar diffusion and can produce dense substructure without large wave-amplitudes. Our results also explain the observed transition from Alfvenic turbulent motion at large scales to subsonic motions at the level of dense cores. △ Less

Submitted 31 March, 2008; originally announced March 2008.

Comments: 7 pages; 7 figures; accepted by A&A

The generation of optical emission-line filaments in galaxy clusters

Authors: Edward Pope, Julian Pittard, Thomas Hartquist, Sam Falle

Abstract: Recent data support the idea that the filaments observed in H_alpha emission near the centres of some galaxy clusters were shaped by bulk flows within their intracluster media. We present numerical simulations of evaporated clump material interacting with impinging winds to investigate this possibility. In each simulation, a clump falls due to gravity while the drag of a wind retards the fall of… ▽ More Recent data support the idea that the filaments observed in H_alpha emission near the centres of some galaxy clusters were shaped by bulk flows within their intracluster media. We present numerical simulations of evaporated clump material interacting with impinging winds to investigate this possibility. In each simulation, a clump falls due to gravity while the drag of a wind retards the fall of evaporated material leading to elongation of the tail. However, we find that long filaments can only form if the outflowing wind velocity is sufficiently large, 10^8 cm s^-1. Otherwise, the tail material sinks almost as quickly as the cloud. For reasonable values of parameters, the morphological structure of a tail is qualitatively similar to those observed in clusters. Under certain conditions, the kinematics of the tail resemble those reported in Hatch et al.(2006). A comparison of the observations with the numerical results indicates that the filaments are likely to be a few tens of Myrs old. We also present arguments which suggest that the momentum transfer, from an outflowing wind, in the formation of these filaments is probably significant. As a result, tail formation could play a role in dissipating some of the energy injected by a central AGN close to the cluster centre where it is needed most. The trapping of energy by the cold gas may provide an additional feedback mechanism that helps to regulate the heating of the central regions of galaxy clusters and couple the AGN to the ICM. △ Less

Submitted 12 February, 2008; v1 submitted 14 January, 2008; originally announced January 2008.

Comments: 15 pages, 18 figures. Altered density values in table 1, which were a factor of 1000 too small. Also extended caption for table 1 to show the mass outflow rate. Accepted for publication in MNRAS

Shock-triggered formation of magnetically-dominated clouds

Authors: S. Van Loo, S. A. E. G. Falle, T. W. Hartquist, T. J. T. Moore

Abstract: To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud, we follow the dynamical evolution of the cloud with a time-dependent axisymmetric magnetohydrodynamic code. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium interacts with the cloud. As a fast-mode… ▽ More To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud, we follow the dynamical evolution of the cloud with a time-dependent axisymmetric magnetohydrodynamic code. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium interacts with the cloud. As a fast-mode shock propagates through the cloud, the gas behind it becomes thermally unstable. The $β$ value of the gas also becomes much smaller than the initial value of order unity. These conditions are ideal for magnetohydrodynamic waves to produce high-density clumps embedded in a rarefied warm medium. A slow-mode shock follows the fast-mode shock. Behind this shock a dense shell forms, which subsequently fragments. This is a primary region for the formation of massive stars. Our simulations show that only weak and moderate-strength shocks can form cold clouds which have properties typical of giant molecular clouds. △ Less

Submitted 4 June, 2007; originally announced June 2007.

Comments: 7 pages, 6 figures, accepted by Astronomy and Astrophysics