Free Access
Issue
A&A
Volume 660, April 2022
Article Number L7
Number of page(s) 9
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202243322
Published online 11 April 2022

© ESO 2022

1. Introduction

The role of magnetic fields on different scales of star formation (from the formation of clouds to the formation of stars) is poorly understood. A major limitation in this understanding is the lack of knowledge of the three-dimensional (3D) magnetic fields in the interstellar medium (ISM). Multi-wavelength magnetic field observations are necessary to determine the 3D magnetic fields associated with star-forming regions.

Far-infrared dust polarized emission (e.g., Houde et al. 2004; Poidevin et al. 2011; Fissel et al. 2016; Pattle & Fissel 2019; Doi et al. 2020) and near-infrared and/or optical (e.g., Pereyra & Magalhães 2004; Clemens et al. 2020) starlight polarization observations have enabled us to probe the plane-of-sky orientation of magnetic fields (B) in a number of molecular clouds. Dust polarization observations revealed that magnetic field lines tend to be perpendicular to high column density (> 1021.7 cm−2) filamentary structures (Planck Collaboration Int. XXXII 2016; Planck Collaboration Int. XXXV 2016), allowing for greater mass accumulations and denser filaments (e.g., Inoue et al. 2018; Hennebelle & Inutsuka 2019).

Radio observations of Faraday rotation (Tahani et al. 2018) and Zeeman measurements (e.g., Goodman et al. 1989; Heiles 1997; Troland & Crutcher 2008) have been used to probe the line-of-sight magnetic fields (B) of molecular clouds. Tahani et al. (2018) developed a new technique based on Faraday rotation measurements for determining the strength and direction of B associated with molecular clouds1. In this technique they incorporated an approach based on relative measurements to estimate the amount of rotation measure (RM) caused by molecular clouds, using RM data from Taylor et al. (2009). To extract the magnetic field strengths from the RMs, they estimated the electron column density of the molecular cloud at the position of each RM point using a chemical evolution code (Gibson et al. 2009) and extinction maps for each cloud.

The observations of B and B are slowly paving the way for the determination of 3D magnetic fields of molecular clouds. Tahani et al. (2019) investigated the 3D magnetic field morphology of the Orion A molecular cloud using both B (Tahani et al. 2018) and B (Planck Collaboration Int. XXXV 2016) observations. The line-of-sight magnetic field observations in this region showed that the B reverses direction from one side of this filamentary-shaped cloud to the other (perpendicular to the filament axis; see Fig. 1). Tahani et al. (2019) constructed models to account for this reversal and compared their synthetic observations to the observed B and B. Using Monte Carlo simulations and chi-squared probability values (and by examining systematic biases between the two observing techniques), they concluded that an arc-shaped2 magnetic field morphology was the most likely magnetic structure for the Orion A molecular cloud.

thumbnail Fig. 1.

Observed magnetic fields of the Orion A molecular cloud. The B data were obtained by Tahani et al. (2018) and the B lines were observed by the Planck Space Observatory. The blue (red) circles show magnetic fields toward (away from) the observer, and the size of the circles represents the strength of the magnetic fields. The red and drapery lines show the B orientation. The background grayscale image is the visual extinction map obtained by Kainulainen et al. (2009). Tahani et al. (2019) suggested that an arc-shaped magnetic field morphology (as illustrated in the inset) is the most probable model that can explain the observed B reversal across the cloud.

This arc-shaped morphology is illustrated in the inset of Fig. 1. The background grayscale image in this figure depicts the visual extinction map of Orion A (in units of magnitude of visual extinction or AV; Kainulainen et al. 2009), and the red and drapery lines (made using the line integration convolution technique3, Cabral & Leedom 1993) show the B observed by the Planck Space Observatory, with blue (red) circles indicating magnetic fields toward (away from) us. The magnetic field on the Galactic north side of the cloud points toward us, while on the Galactic south side it points away from us. This B reversal was previously observed using Zeeman measurements in the study of Heiles (1997), which predicted an arc-shaped magnetic field morphology caused by recurrent shocks from nearby supernovae in the Orion-Eridanus superbubble.

An arc-shaped magnetic field morphology has also been generated in the ideal magnetohydrodynamic simulations of Inoue et al. (2018) and Inoue & Fukui (2013), and predicted in the Inutsuka et al. (2015) molecular-cloud-formation model (also see Inutsuka et al. 2016). Multiple compressions caused by expanding ISM bubbles are required for the formation of filamentary molecular clouds in the Inutsuka et al. (2015) model, where magnetic field lines with a component perpendicular to the direction of bubble expansion or propagation can be bent. Figure 2 depicts a schematic cartoon of this interaction, which occurs between a relatively dense cloud (∼103 cm−3) and a shock wave propagating in low density gas (H I) and bends the field lines around the formed filament. We refer to this cloud formation model as the shock-cloud interaction (SCI) model.

thumbnail Fig. 2.

Formation of an arc-shaped magnetic field morphology around filamentary molecular clouds, as described by Inoue et al. (2018). After the interaction between the cloud and shock-front, with the velocity v in the co-moving frame of the shock-front, a filamentary structure is formed (shown in an end-on view on the right). B shows the direction of the initial magnetic field before the collision and the morphology of the magnetic field after the collision.

While previous studies described the magnetic field morphology of the Orion A cloud as having an arc shape, its complete 3D orientation and direction (particularly when projected onto the plane of the sky) remained unknown, even in the 3D magnetic field study by Tahani et al. (2019). In this study the complete 3D magnetic field morphology of the Orion A molecular cloud, including its direction and bending orientation, is reconstructed from the data. To this end, we employ Galactic magnetic field (GMF) models, the B data of Tahani et al. (2018), and the B observed by the Planck Space Observatory (whereas in Tahani et al. 2019 only B and B observations were used). We present our approach and results in Sect. 2, discuss the role of surrounding structures, bubbles, and events in formation of this 3D field in Sect. 3, and provide a summary and conclusion in Sect. 4. Supplemental material, including the data used in this study, are provided in the appendices.

2. Results: 3D magnetic field morphology of Orion A

We recreated the 3D magnetic field shape of the Orion A molecular cloud using GMF vectors as initial magnetic fields, the orientation of B reversal, and B morphology (under the assumption of an arc-shaped magnetic morphology). To approximate the direction of the initial magnetic fields, we used the GMF model of Jansson & Farrar (2012), neglecting the isotropic random field components (caused by ISM turbulence). We refer to this structure as the “Coherent GMF” model. Figure 3 illustrates these GMF vectors projected onto the plane of the sky (red arrows) at the location of Orion A.

thumbnail Fig. 3.

Galactic magnetic field associated with the Orion A cloud. The background grayscale image shows the extinction map from Kainulainen et al. (2009). The yellow line represents the division between the tail and the head of the cloud. The region of Orion A with longitudes less than ∼210° is referred to as the “head” of Orion A, while the portion with longitudes greater than ∼210° is referred to as the “tail”. The red vectors show the Coherent GMF model projected onto the plane of the sky. The cyan line illustrates an approximate filament axis.

To best describe the GMF vectors, we employed a frame of reference in this region (see Fig. 4), with its axes pointing in the increasing directions of longitude (), latitude (), and distance () at Orion A’s plane-of-sky location (where , , and are unit vectors). The GMF direction in this region can be described as a unit vector of . This vector appears mostly parallel to Orion A when projected onto the plane of the sky, but has a large component along the line of sight (pointing toward us). This vector is consistent with previous GMF studies (e.g., Van Eck et al. 2011, see their Fig. 6), with the same inclination angle and direction (45° to 50° with respect to the plane of the sky; e.g., Oren & Wolfe 1995; Heiles 1997).

thumbnail Fig. 4.

Cartoon showing the 3D morphology of the Orion A magnetic field. The bent gray cylinder and the red vector depict the Orion A filamentary structure and the mean 3D Galactic magnetic field in the region, respectively. The blue arrow shows the 3D magnetic field morphology of the cloud and is consistent with the B and B observations shown in Fig. 1. Top panel: projection of the cloud onto the plane of the sky. Middle and lower panels: different viewing angles of the cloud. The black arrow denotes the direction in which rotations of up to 50° may be possible.

We also needed to account for the inclination angle of Orion A when reconstructing its 3D magnetic field. Großschedl et al. (2018) demonstrated that most of Orion A has a high inclination angle with respect to the plane of the sky. They divided the cloud into a head (longitudes less than ∼210°) and a tail (longitudes greater than ∼210°; see Fig. 3), with the head running approximately parallel to the plane of the sky and the tail trending into the line of sight at an approximately 70° inclination angle, as illustrated in Fig. 4.

We reconstructed the arc-shaped morphology depicted in Fig. 4 by connecting the approximate GMF vectors to the B observations via an arc that accounts for the 3D orientation of the cloud and matches the B morphology and the relative B strengths. The initial magnetic field (the GMF) is depicted in this figure as a red vector pointing toward us. The blue vector illustrates the reconstructed arc-shaped magnetic field morphology, which looks semi-convex (from our point of view) and points in the direction of decreasing latitude.

This 3D magnetic field morphology is an approximation of a large-scale configuration that ignores smaller-scale variations and fluctuations, enabling us to determine the direction and orientation of the arc, as well as the direction of magnetic fields projected onto the plane of sky. Therefore, small variations in this field could be possible as long as they remain consistent with the B and B observations, and with the GMF vectors. While the majority of these variations are insignificant, rotations up to a maximum of 50° along the black arrow shown in the middle and lower panels of Fig. 4 may be possible (also see Fig. B.1), covering possible field shapes from the front to the back of the cloud (that are consistent with the observations). The field shapes in this possible range remain semi-convex from our point of view and point in the decreasing latitude direction when projected onto the plane of the sky.

Moreover, while we did not directly consider field strengths when reconstructing the field lines, the results are consistent with observed strengths; the cloud has a stronger line-of-sight field component to its Galactic north side than to its Galactic south side. This is consistent with the B observations of Tahani et al. (2018) and Zeeman measurements of Heiles (1997). The error-weighted average strength of the B observations in both studies is twice as strong on the cloud’s Galactic north side as on the cloud’s Galactic south side. Finally, while we believe that this 3D field is the most probable and natural field morphology for the Orion A cloud, it does not completely rule out other possibilities. Future observations with a high rotation measure source density and improved B resolution will allow us to more accurately and precisely determine the 3D fields.

3. Discussion

Our reconstructed 3D magnetic morphology is consistent with the SCI cloud-formation model, in which the surrounding environment (e.g., bubbles, events, structures, or sheet-like clouds) influenced and bent the field lines. In this section we discuss the formation of this arc-shaped morphology and the environmental influences that shaped it.

3.1. Bubbles influencing the Orion A cloud

Orion-Eridanus superbubble. Orion A (located at 432 ± 2 pc from the Sun; Zucker et al. 2019) was likely shaped by energy release from the Orion OB association (Bally 2008). From the Sun the cloud is seen in the projected interior of the Orion-Eridanus superbubble. Its distance puts it on the far side of the superbubble, but its location within, along, or beyond the bubble rim remains uncertain because of the complex composite structure of the superbubble (Ochsendorf et al. 2015; Joubaud et al. 2019). The stellar content (Voss et al. 2010), stellar age gradient (Bouy & Alves 2015; Zari et al. 2017), and X-ray temperature gradient (Joubaud et al. 2019) all indicate that the bubble likely evolved in time and space from a near (150–200 pc away) to far distance (close to Orion A and B) after a series of supernova events (in the past 10–20 Myr). The approximate outline of the superbubble (based on Ochsendorf et al. 2015) is depicted in Fig. 5 as a white circle, with the background colors of light green and orange representing thermal dust (at 545 GHz) and Hα observations, respectively. Figure 6 illustrates a 3D approximation of the location of this superbubble as a gray ellipsoid, based on the 3D models of Pon et al. (2014, 2016).

thumbnail Fig. 5.

Observations of the Orion-Eridanus superbubble and its nested bubbles. The white, green dashed, and red dash-dotted circles illustrate the approximate outline of the Orion-Eridanus superbubble, dust ring, and Barnard’s loop on the plane of the sky, respectively. The light green and orange background colors represent the observations of thermal dust (at 545 GHz) obtained by the Planck Space Observatory and Hα emission (Finkbeiner 2003), respectively. The bright green region shows Orion A.

thumbnail Fig. 6.

Simplified 3D approximation of the Orion-Eridanus superbubble and its associated bubbles. Orion A is shown as an orange bent cylinder and the Sun is depicted as a yellow circle. The gray ellipsoid and red and green spheres represent the Orion-Eridanus superbubble, Barnard’s loop, and the dust ring, respectively. Due to uncertainties in the models of the Orion-Eridanus superbubble, the dust ring might also be embedded within the superbubble.

Großschedl et al. (2021), Kounkel (2020), and Pellizza et al. (2005) found coherent stellar proper motion associated with the Orion region and the superbubble. Coherent proper motions observed by Großschedl et al. (2020) and Kounkel (2020) are associated with the young stellar objects in the Orion A and B clouds, implying the influence of a feedback-driven event. Großschedl et al. (2020) refers to this event as Orion-BB (big blast), whereas Kounkel (2020) links the coherent velocities to Barnard’s loop. In general, these coherent velocities are likely indicative of feedback-driven (e.g., supernovae) impacts on the Orion A and B clouds and their young stellar objects.

Barnard’s loop. Barnard’s loop (Barnard 1894) is a complete bubble located within the Orion-Eridanus superbubble (Ochsendorf et al. 2015). The Barnard’s loop bubble, with an estimated age of 3 × 105 yr, is expanding at a velocity of 100 km s−1 (Ochsendorf et al. 2015), while the Orion-Eridanus superbubble has an expansion velocity of ∼20 km s−1 (Joubaud et al. 2019). This bubble is located between 340 pc and 400 pc from us (Großschedl et al. 2021). Using the distances of Barnard’s loop and the Orion A cloud, we find that the Barnard’s loop bubble is likely in contact with (and has interacted with) Orion A’s head, but not its tail, as depicted in Fig. 6. This could account for the tilt of Orion A’s head and the fact that the star formation rate in Orion A’s head is an order of magnitude greater than in its tail (Großschedl et al. 2020). We predict that as the Barnard’s loop bubble expands, it will interact with Orion A and trigger a new star formation sequence in the tail.

Orion dust ring.Schlafly et al. (2015) mapped the 3D dust and found a dust ring (indicating a bubble origin) between 400 pc and 550 pc (from us) in the Orion complex region. This dust ring is depicted in Fig. 5 and is seen at more positive latitudes compared to Orion A. Schlafly et al. (2015) estimate the age of the bubble to be around 10 Myr or greater because they find no evidence of Hα associated with the ring. A closer examination of Fig. 5, however, indicates that this ring may be visible in Hα observations. We roughly approximate the ring (and its original bubble) as a sphere (shown in Fig. 6) with a radius of 150 pc and a center location of ℓ = 212°, b = −11.5°, and d = 175 pc.

3.2. Gradual formation of the arc-shaped field morphology

The bubbles detailed in Sect. 3.1 have influenced the evolution (and the formation) of the Orion A cloud and its magnetic field lines. These influences can be summarized in two steps. First, recurrent interactions with old events or bubbles pushed the Orion A cloud, its ambient environment, and the field lines (toward more positive latitudes), resulting in large-scale bending of the initial Galactic magnetic field, as predicted by the simulations of Kim & Ostriker (2015). Subsequently, further interactions (likely by the dust ring) pushed the H I gas around the Orion A cloud (toward more negative latitudes), creating a sharp arc-shaped magnetic morphology associated with the Orion A cloud, as depicted in Fig. 4.

This two-step process is similar to that described in Tahani et al. (2022) for the Perseus molecular cloud (see Fig. D.2 and Appendix D). In the first step, the Galactic magnetic fields are bent approximately tangentially to expanding objects or bubbles (e.g., Orion-BB event or Orion-Eridanus superbubble). This results in magnetic field lines that are bent at Orion A’s location, point toward us on its Galactic north side, and are mostly parallel to the plane of the sky on its Galactic south side (and perpendicular to the formed filament main axis). Subsequently, the interaction with the dust ring bends these field lines even more, causing them to point away from us on the cloud’s Galactic south side. Therefore, the field lines on the Galactic north side of the cloud have a greater component along the line of sight compared to the Galactic south side of the cloud, as discussed in Sect. 2.

To examine this view further, we estimate the gas and magnetic field pressures (see Appendix E) before the formation of the Orion A cloud and in its present state. During the initial stages of the Orion A cloud formation, we find that the gas pressure is greater than the magnetic field pressure (using a Galactic magnetic field strength of ∼5 μG, particle volume density of a few to 100 cm−3, and ambient temperature of 1000 K and higher; Joubaud et al. 2019), implying that recurrent supernovae can easily bend the magnetic field lines, as seen in simulations of Kim & Ostriker (2015). Furthermore, using an average non-thermal velocity dispersion of ∼3 km s−1, we find an Alfvén Mach number of ∼3, indicating that the field lines can retain a memory of the large-scale initial field morphology (with small-scale variations; cf. Han & Zhang 2007). As a result, when the field lines interact with bubbles and events in this region (such as the Orion-Eridanus superbubble or the Orion-BB event) and bend, they remain largely coherent (Kim & Ostriker 2015), rather than becoming completely distorted and perturbed (Li & Klein 2019).

To estimate the gas and magnetic field pressures in the Orion A molecular cloud in its current state, we use the error-weighted mean magnetic fields on two sides of the cloud (for detections pointing toward us and away from us), which are 87 μG and −45 μG, respectively, based on the observations of Tahani et al. (2018). Using particle volume density and temperature values of ∼104 cm−3 (e.g., Castets et al. 1990; Dutrey et al. 1993; Johnstone & Bally 1999a,b) and 25 K (e.g., Mitchell et al. 2001; Johnstone & Bally 2006; Buckle et al. 2012; Schnee et al. 2014), we find that the magnetic field pressure of the Orion A cloud is approximately one order of magnitude greater than its gas pressure (on both sides of the cloud). Using an average non-thermal velocity dispersion of 2 km s−1 (Goicoechea et al. 2020), we find Alfvén Mach numbers of 0.4 and 0.8 for the two sides of the Orion A molecular cloud. This implies that the field lines do not deform and if their magnetic field morphology changes (due to interaction with the environment), they retain a memory of their previous field morphology.

This gradual bending of the field lines is consistent with the cloud-formation model of Inutsuka et al. (2015). We also determine the H I and CO line-of-sight velocities in this region (see Appendix C) and find no significant offset between the two, indicating that the offset velocities caused by bubble interactions have dissipated or are mostly in the plane of the sky. Given the morphology of the Orion A cloud (which is inclined along the line of sight) and that the GMF vectors are parallel to the longitude axis and point toward us along the line of sight, the interactions necessary for the formation of the arc-shaped magnetic morphology are mostly parallel to the plane of the sky. Therefore, we expect that any velocity offsets that may hint at arc-shaped magnetic field morphology formation (and perhaps are not yet completely dissipated) should be parallel to the plane of the sky.

4. Summary and conclusions

We determined the large-scale (and approximate) 3D magnetic field shape of the Orion A molecular cloud using Galactic magnetic field models as well as present line-of-sight and plane-of-sky magnetic field observations. From our perspective, this 3D field is generally semi-convex and points toward the decreasing longitude and decreasing latitude directions when projected onto the plane of the sky. To our knowledge, this is the first time that the complete 3D large-scale (∼ a few to ∼100 pc) magnetic field of the Orion A cloud has been reconstructed (including its plane-of-sky direction). We suggest that the Orion-Eridanus superbubble (or events within it, such as the Orion-BB event) and the dust ring in this region are largely responsible for the development of the arc-shaped magnetic field morphology of Orion A.

Reconstructing the 3D magnetic field morphology of Orion A, which is consistent with the Planck observations when projected onto the plane of the sky, relies mainly on the initial Galactic magnetic fields and the B observations. This 3D magnetic field is a large-scale approximation that neglects smaller-scale distortions or entanglements in the field lines. Because we are interested in relatively large-scale and approximate fields, the presence of low-density regions in the foreground of Orion A (especially those smaller than the size of the Orion A cloud and the scales at which the B reversal was observed; Rezaei et al. 2020) has no effect on our results.


2

Also referred to as bow-shaped; pronounced /bō/ as in rainbow or bow and arrow.

Acknowledgments

We would like to express our appreciation to the anonymous referee for their thoughtful comments that helped us improve the results and paper. M.T. is grateful for the helpful discussion with Pak-Shing Li. We have used LaTeX, Python and its associated libraries including astropy (Astropy Collaboration 2013) and plotly (Plotly Technologies Inc. 2015), PyCharm, Jupyter notebook, SAO Image DS9, the Starlink (Currie et al. 2014) software, the Hammurabi code, Shapr3D, and Adobe Draw for this work. For our line integration convolution plot, we used a Python function originally written by Susan Clark. QuillBot (https://quillbot.com/) was used for language editing purposes. We thank and acknowledge the Interstellar Institute and Paris-Saclay University’s Institut Pascal for hosting discussions that aided the development of this work. The Dunlap Institute is funded through an endowment established by the David Dunlap family and the University of Toronto. J.L.W. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2015-05948, and of the Canada Research Chairs program. Y.D. acknowledges the support of JSPS KAKENHI grant 18H01250. M.H. acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 772663).

References

  1. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Bally, J. 2008, in Overview of the Orion Complex, ed. B. Reipurth, 4, 459 [Google Scholar]
  3. Barnard, E. E. 1894, Popular Astron., 2, 151 [Google Scholar]
  4. Bouy, H., & Alves, J. 2015, A&A, 584, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Buckle, J. V., Davis, C. J., Di Francesco, J., et al. 2012, MNRAS, 422, 521 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cabral, B., & Leedom, L. C. 1993, Proceedings of the 20th annual conference on Computer graphics and interactive techniques [Google Scholar]
  7. Castets, A., Duvert, G., Dutrey, A., et al. 1990, A&A, 234, 469 [NASA ADS] [Google Scholar]
  8. Clemens, D. P. 1985, ApJ, 295, 422 [NASA ADS] [CrossRef] [Google Scholar]
  9. Clemens, D. P., Cashman, L. R., Cerny, C., et al. 2020, ApJS, 249, 23 [NASA ADS] [CrossRef] [Google Scholar]
  10. Currie, M. J., Berry, D. S., Jenness, T., et al. 2014, in Astronomical Data Analysis Software and Systems XXIII, eds. N. Manset,&P. Forshay, ASP Conf. Ser., 485, 391 [NASA ADS] [Google Scholar]
  11. Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 [Google Scholar]
  12. Doi, Y., Hasegawa, T., Furuya, R. S., et al. 2020, ApJ, 899, 28 [Google Scholar]
  13. Dutrey, A., Duvert, G., Castets, A., et al. 1993, A&A, 270, 468 [Google Scholar]
  14. Finkbeiner, D. P. 2003, ApJS, 146, 407 [Google Scholar]
  15. Fissel, L. M., Ade, P. A. R., Angilè, F. E., et al. 2016, ApJ, 824, 134 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gibson, D., Plume, R., Bergin, E., Ragan, S., & Evans, N. 2009, ApJ, 705, 123 [NASA ADS] [CrossRef] [Google Scholar]
  17. Goicoechea, J. R., Pabst, C. H. M., Kabanovic, S., et al. 2020, A&A, 639, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Goodman, A. A., Crutcher, R. M., Heiles, C., Myers, P. C., & Troland, T. H. 1989, ApJ, 338, L61 [NASA ADS] [CrossRef] [Google Scholar]
  19. Großschedl, J. E., Alves, J., Meingast, S., et al. 2018, A&A, 619, A106 [Google Scholar]
  20. Großschedl, J. E., Alves, J., Meingast, S., & Hasenberger, B. 2020, in Origins: From the Protosun to the First Steps of Life, eds. B. G. Elmegreen, L. V. Tóth, & M. Güdel, 345, 27 [Google Scholar]
  21. Großschedl, J. E., Alves, J., Meingast, S., & Herbst-Kiss, G. 2021, A&A, 647, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Han, J. L., & Zhang, J. S. 2007, A&A, 464, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Heiles, C. 1997, ApJS, 111, 245 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hennebelle, P., & Inutsuka, S.-I. 2019, Front. Astron. Space Sci., 6, 5 [CrossRef] [Google Scholar]
  25. HI4PI Collaboration (Ben Bekhti, N., et al.) 2016, A&A, 594, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Houde, M., Dowell, C. D., Hildebrand, R. H., et al. 2004, ApJ, 604, 717 [NASA ADS] [CrossRef] [Google Scholar]
  27. Imara, N., & Blitz, L. 2011, ApJ, 732, 78 [NASA ADS] [CrossRef] [Google Scholar]
  28. Inoue, T., & Fukui, Y. 2013, ApJ, 774, L31 [Google Scholar]
  29. Inoue, T., Hennebelle, P., Fukui, Y., et al. 2018, PASJ, 70, S53 [NASA ADS] [Google Scholar]
  30. Inutsuka, S.-I., Inoue, T., Iwasaki, K., & Hosokawa, T. 2015, A&A, 580, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Inutsuka, S. I., Inoue, T., Iwasaki, K., Hosokawa, T., & Kobayashi, M. I. N. 2016, in From Interstellar Clouds to Star-Forming Galaxies: Universal Processes?, eds. P. Jablonka, P. André, & F. van der Tak, 315, 61 [NASA ADS] [Google Scholar]
  32. Ishii, S., Nakamura, F., Shimajiri, Y., et al. 2019, PASJ, 71, S9 [Google Scholar]
  33. Jansson, R., & Farrar, G. R. 2012, ApJ, 757, 14 [NASA ADS] [CrossRef] [Google Scholar]
  34. Johnstone, D., & Bally, J. 1999a, ApJ, 510, L49 [Google Scholar]
  35. Johnstone, D., & Bally, J. 1999b, in The Physics and Chemistry of the Interstellar Medium, eds. V. Ossenkopf, J. Stutzki, & G. Winnewisser [Google Scholar]
  36. Johnstone, D., & Bally, J. 2006, ApJ, 653, 383 [NASA ADS] [CrossRef] [Google Scholar]
  37. Joubaud, T., Grenier, I. A., Ballet, J., & Soler, J. D. 2019, A&A, 631, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Kainulainen, J., Beuther, H., Henning, T., & Plume, R. 2009, A&A, 508, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Kerp, J., Winkel, B., Ben Bekhti, N., Flöer, L., & Kalberla, P. M. W. 2011, Astron. Nachr., 332, 637 [NASA ADS] [CrossRef] [Google Scholar]
  40. Kim, C.-G., & Ostriker, E. C. 2015, ApJ, 802, 99 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kong, S., Lada, C. J., Lada, E. A., et al. 2015, ApJ, 805, 58 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kounkel, M. 2020, ApJ, 902, 122 [CrossRef] [Google Scholar]
  43. Li, P. S., & Klein, R. I. 2019, MNRAS, 485, 4509 [NASA ADS] [CrossRef] [Google Scholar]
  44. Ma, Y.-H., Wang, H.-C., Li, C., & Yang, J. 2020, Res. Astron. Astrophys., 20, 060 [CrossRef] [Google Scholar]
  45. McClure-Griffiths, N. M., Pisano, D. J., Calabretta, M. R., et al. 2009, ApJS, 181, 398 [NASA ADS] [CrossRef] [Google Scholar]
  46. Mitchell, G. F., Johnstone, D., Moriarty-Schieven, G., Fich, M., & Tothill, N. F. H. 2001, ApJ, 556, 215 [Google Scholar]
  47. Ochsendorf, B. B., Brown, A. G. A., Bally, J., & Tielens, A. G. G. M. 2015, ApJ, 808, 111 [NASA ADS] [CrossRef] [Google Scholar]
  48. Oren, A. L., & Wolfe, A. M. 1995, ApJ, 445, 624 [NASA ADS] [CrossRef] [Google Scholar]
  49. Pattle, K., & Fissel, L. 2019, Front. Astron. Space Sci., 6, 15 [Google Scholar]
  50. Pellizza, L. J., Mignani, R. P., Grenier, I. A., & Mirabel, I. F. 2005, A&A, 435, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Pereyra, A., & Magalhães, A. M. 2004, ApJ, 603, 584 [CrossRef] [Google Scholar]
  52. Planck Collaboration Int. XLII. 2016, A&A, 596, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Planck Collaboration Int. XXXII. 2016, A&A, 586, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Planck Collaboration Int. XXXV. 2016, A&A, 586, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Plotly Technologies Inc. 2015, Collaborative Data Science [Google Scholar]
  56. Poidevin, F., Bastien, P., & Jones, T. J. 2011, ApJ, 741, 112 [NASA ADS] [CrossRef] [Google Scholar]
  57. Pon, A., Johnstone, D., Bally, J., & Heiles, C. 2014, MNRAS, 444, 3657 [NASA ADS] [CrossRef] [Google Scholar]
  58. Pon, A., Ochsendorf, B. B., Alves, J., et al. 2016, ApJ, 827, 42 [NASA ADS] [CrossRef] [Google Scholar]
  59. Rezaei, Kh. S., Bailer-Jones, C. A. L., Soler, J. D., & Zari, E. 2020, A&A, 643, A151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Rice, T. S., Goodman, A. A., Bergin, E. A., Beaumont, C., & Dame, T. M. 2016, ApJ, 822, 52 [Google Scholar]
  61. Schlafly, E. F., Green, G., Finkbeiner, D. P., et al. 2015, ApJ, 799, 116 [NASA ADS] [CrossRef] [Google Scholar]
  62. Schnee, S., Mason, B., Di Francesco, J., et al. 2014, MNRAS, 444, 2303 [NASA ADS] [CrossRef] [Google Scholar]
  63. Tahani, M., Plume, R., Brown, J. C., & Kainulainen, J. 2018, A&A, 614, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Tahani, M., Plume, R., Brown, J. C., Soler, J. D., & Kainulainen, J. 2019, A&A, 632, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Tahani, M., Lupypciw, W., Glover, J., et al. 2022, A&A, in press, https://doi.org/10.1051/0004-6361/202141170 [Google Scholar]
  66. Taylor, A. R., Stil, J. M., & Sunstrum, C. 2009, ApJ, 702, 1230 [Google Scholar]
  67. Troland, T. H., & Crutcher, R. M. 2008, ApJ, 680, 457 [NASA ADS] [CrossRef] [Google Scholar]
  68. Van Eck, C. L., Brown, J. C., Stil, J. M., et al. 2011, ApJ, 728, 97 [Google Scholar]
  69. Voss, R., Diehl, R., Vink, J. S., & Hartmann, D. H. 2010, A&A, 520, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Waelkens, A., Jaffe, T., Reinecke, M., Kitaura, F. S., & Enßlin, T. A. 2009, A&A, 495, 697 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Winkel, B., Kerp, J., Flöer, L., et al. 2016, A&A, 585, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Zari, E., Brown, A. G. A., de Bruijne, J., Manara, C. F., & de Zeeuw, P. T. 2017, A&A, 608, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Zucker, C., Speagle, J. S., Schlafly, E. F., et al. 2019, ApJ, 879, 125 [Google Scholar]

Appendix A: Data used in the study

In this study we used estimates of the initial magnetic field direction, line-of-sight magnetic fields, and velocities associated with each cloud. We used the Jansson & Farrar (2012) GMF model to determine the Galactic magnetic field. We employed the catalog of Tahani et al. (2018) for the B magnetic field information, and available CO and H I observations for the velocities.

A.1. Galactic magnetic field

The Jansson & Farrar (2012) model includes a two-dimensional (2D) thin-disk field component that is tightly coupled to the Galactic spiral arms, an azimuthal–toroidal halo field component, and an X-shaped vertical–out-of-plane field component. To estimate the GMF, we used the Hammurabi program4 (Waelkens et al. 2009), which is a synchrotron modeling code that has been used in different studies (e.g., Planck Collaboration Int. XLII 2016; Jansson & Farrar 2012). We find the GMF vectors within a box around the Orion A cloud (160 pc × 160 pc × 200 pc), in the longitude range of 205° to 218° and latitude range of −26° to −13°. We set a resolution of one GMF vector per 2 pc × 2 pc × 2 pc.

A.2. Line-of-sight magnetic field

We employed the B observations of Tahani et al. (2018), in which they used Faraday rotation measurements to determine B of molecular clouds. They used RM point sources along with an on-off approach based on relative measurements to decouple the molecular clouds’ contribution to RM from that caused by the rest of the Galaxy. They then calculated the strength of B employing a chemical evolution code and the Kainulainen et al. (2009) extinction maps (for more details, see Tahani et al. 2018).

A.3. Velocity information

We considered the available CO and H I velocities to explore the line-of-sight velocities of the cloud and its surrounding environment. We obtained H I velocity information from the all-sky database of the H I 4π Survey (HI4PI; HI4PI Collaboration 2016), which is based on the Effelsberg-Bonn H I Survey (EBHIS; Kerp et al. 2011; Winkel et al. 2016) and the Galactic All-Sky Survey (GASS; McClure-Griffiths et al. 2009). We used the radial velocities from the Dame et al. (2001) carbon monoxide survey to determine the cloud’s CO velocity. This catalog is a survey of the 12CO J(1-0) spectral line of the Galaxy.

Appendix B: 3D magnetic field sensitivity

As discussed in Sect. 2, the reconstructed 3D field is an approximation of the large-scale field, which allows for small-scale variations such as small rotations along the black arrow shown in the middle panel of Fig. B.1. Rotations in other directions, on the other hand, are more sensitive and limited as they do not produce the observed B or B . We note that these variations must remain consistent with the GMF, and with the B and B observations. The magnetic field that we reconstructed5 in Sect. 2 is depicted as 3D B (1) in Fig. B.1. Examining a large number of 3D field varieties, we find that rotations along the black arrow up to 50° may be possible (covering possible field morphologies from the front to the back of the cloud), resulting in the arc-shaped field depicted as 3D B(2) in Fig. B.1, which is semi-convex from our perspective and points in the direction of decreasing latitude.

thumbnail Fig. B.1.

Cartoon showing the sensitivity of the reconstructed 3D magnetic field. The Galactic magnetic field, the Orion A cloud, and the reconstructed field are shown with the red vector, bent gray cylinder, and the blue 3D B (1) arrow, respectively. Top panel: Projection onto the plane of the sky. Middle and lower panels: Different viewing angles of the cloud. The black arrow in the middle panel illustrates the direction of rotation to obtain 3D B(2), which remains consistent with the observations and is semi-convex from the observer’s point of view.

Appendix C: Velocity observations

To investigate the line-of-sight velocities of the cloud thoroughly, we first calculated the line-of-sight Galactic rotation velocities. We find a Galactic rotation velocity (LSR) of ∼ + 4 km s−1 for the Orion A cloud (l ≃ 211° and d ≃ 432 pc), using the model of Clemens (1985) with the IAU standard values of the solar distance from the Galactic center (8.5 kpc) and orbital velocity (220 km s−1).

For H I velocities, we selected locations that have a single H I peak emission and excluded those with multiple peaks, absorption, or self-absorption, as discussed in Tahani et al. (2022). Similarly, we chose points whose CO spectrum can be described primarily by a single Gaussian fit. For more accuracy, we took H I and CO data from similar coordinates. These H I and CO velocities agree with other studies at similar coordinates (e.g., Kong et al. 2015; Rice et al. 2016; Ishii et al. 2019; Ma et al. 2020). We note that the average CO and HI velocities along the main axis of the cloud are close to Galactic rotation velocities at this location (∼4 to 6 km s−1).

Subsequently, we found the molecular cloud velocities (CO) in the co-moving frame of the H I gas using

(C.1)

where vCO - HI is the cloud CO velocity in the co-moving H I frame, vCO, LSR is the CO velocity in the local standard of rest (LSR) frame, and vHI, LSR is the H I velocity of the region in the LSR frame. To account for the velocity gradients and fluctuations along the cloud, we found vCO - HI on the cloud point by point (considering the peak emission of CO and H I) and then took an average. This obtained value is 0.0 ± 0.5 km s−1 (where the uncertainty value is the standard deviation of the vCO - HI points), indicating that within the uncertainty range of the observations there is no significant line-of-sight offset between the H I and CO velocities in the Orion A cloud (also consistent with Fig. 1 of Imara & Blitz 2011).

This zero km s−1 line-of-sight average CO velocity in the co-moving H I frame emphasizes the importance of considering the plane-of-sky environmental influences and their effect on the GMF vectors. Furthermore, since the GMF vectors fall parallel to the cloud on the plane of the sky and the cloud has a large inclination angle, it is even more critical to consider plane-of-sky influences and events (e.g., coherent velocities observed by Großschedl et al. 2021) in order to fully understand the evolution of 3D magnetic fields in this region.

Appendix D: Orion bubbles

We discuss the bubbles that influence the Orion A cloud in Sect. 3.1 and Sect. 3.2. The presence of these bubbles is further demonstrated by multi-wavelength observations shown in Fig. D.1. Thermal dust, Hα, CO, and H I observations are illustrated in this multi-panel figure, with the Orion-Eridanus superbubble, dust ring, Barnard’s loop, and the λ Orionis ring outlined as white, green dashed, red dash-dotted, and blue dotted circles, respectively.

thumbnail Fig. D.1.

Multi-wavelength observations of the Orion-Eridanus superbubble. The Orion-Eridanus superbubble, dust ring, Barnard’s loop, and the λ Orionis ring are depicted as white, green dashed, red dash-dotted, and blue dotted circles, respectively. Top left panel: Thermal dust observations from the Planck Space Observatory. Top middle panel: observations (Finkbeiner 2003). Top right panel: Composite CO survey of Dame et al. (2001). Lower panel: H I observations at velocities of 2.7 km s−1 and 9.1 km s−1 from HI4PI.

In Sect. 3.2 we discuss the gradual evolution of initial magnetic field lines, resulting in the arc-shaped morphology depicted in Fig. 4. This evolution process is summarized in Fig. D.2: The first step (left panel) involves the bending of Galactic magnetic fields (by the Orion-Eridanus superbubble, or bubbles and events within, such as the Orion-BB event). At this stage field lines are bent at the Orion A location, pointing toward us on its Galactic north side, and mostly parallel to the plane of the sky on the Galactic south side of the cloud (on the plane of the sky appearing mostly perpendicular to the main axis of the filament that is formed). Following that (right panel), additional environmental influences, such as the dust ring, interact with the cloud’s surrounding H I gas, further bending these field lines, causing them to point away from us on the cloud’s Galactic south side.

thumbnail Fig. D.2.

Formation of the arc-shaped magnetic field morphology in the Orion A cloud. The large light blue circle represents a bubble within the Orion-Eridanus superbubble, or the Orion-Eridanus superbubble itself (the blue star represents the progenitor(s)). The small dark blue cylinder illustrates the Orion A filamentary molecular cloud. The 3D magnetic configurations cannot be accurately represented in 2D and a more accurate representation of the field lines is shown in Fig. 4. Left panel: Straight B represents the direction of the initial magnetic field prior to interaction with the bubble or event, while the large-scale mild bending of the magnetic field morphology after interaction with the bubble or event is represented by the curved B. Right panel: Formation of the arc-shaped magnetic field morphology around the Orion A molecular cloud due to further interactions, possibly with the dust ring.

Appendix E: Alfvén Mach number and gas and magnetic pressures

In Sect. 3.2 we calculate magnetic and gas pressures, as well as Alfvén Mach numbers, to elucidate the step-by-step evolution of the field lines that resulted in the arc-shaped morphology shown in Fig. 4. To estimate the magnetic and gas pressures, we used the following equations:

(E.1)

Here PB, B, Pgas, n, kB, and T are the magnetic pressure, total strength of magnetic field, gas pressure, particle volume density, Boltzmann constant, and temperature, respectively.

To determine the Alfvén mach number, we used the equation

(E.2)

where σv is the non-thermal velocity dispersion and vA is the Alfvén wave group velocity. The Alfvén velocity can be obtained using

(E.3)

where B is the strength of magnetic field and ρ is the volume density. Simulations by Li & Klein (2019) suggest that in the presence of strong fields (ℳA ≃ 1), magnetic field lines can coherently bend around a formed filamentary molecular cloud as a result of converging flows, whereas for ℳA ≃ 10, complete distortion of magnetic fields is expected.

All Figures

thumbnail Fig. 1.

Observed magnetic fields of the Orion A molecular cloud. The B data were obtained by Tahani et al. (2018) and the B lines were observed by the Planck Space Observatory. The blue (red) circles show magnetic fields toward (away from) the observer, and the size of the circles represents the strength of the magnetic fields. The red and drapery lines show the B orientation. The background grayscale image is the visual extinction map obtained by Kainulainen et al. (2009). Tahani et al. (2019) suggested that an arc-shaped magnetic field morphology (as illustrated in the inset) is the most probable model that can explain the observed B reversal across the cloud.

In the text
thumbnail Fig. 2.

Formation of an arc-shaped magnetic field morphology around filamentary molecular clouds, as described by Inoue et al. (2018). After the interaction between the cloud and shock-front, with the velocity v in the co-moving frame of the shock-front, a filamentary structure is formed (shown in an end-on view on the right). B shows the direction of the initial magnetic field before the collision and the morphology of the magnetic field after the collision.

In the text
thumbnail Fig. 3.

Galactic magnetic field associated with the Orion A cloud. The background grayscale image shows the extinction map from Kainulainen et al. (2009). The yellow line represents the division between the tail and the head of the cloud. The region of Orion A with longitudes less than ∼210° is referred to as the “head” of Orion A, while the portion with longitudes greater than ∼210° is referred to as the “tail”. The red vectors show the Coherent GMF model projected onto the plane of the sky. The cyan line illustrates an approximate filament axis.

In the text
thumbnail Fig. 4.

Cartoon showing the 3D morphology of the Orion A magnetic field. The bent gray cylinder and the red vector depict the Orion A filamentary structure and the mean 3D Galactic magnetic field in the region, respectively. The blue arrow shows the 3D magnetic field morphology of the cloud and is consistent with the B and B observations shown in Fig. 1. Top panel: projection of the cloud onto the plane of the sky. Middle and lower panels: different viewing angles of the cloud. The black arrow denotes the direction in which rotations of up to 50° may be possible.

In the text
thumbnail Fig. 5.

Observations of the Orion-Eridanus superbubble and its nested bubbles. The white, green dashed, and red dash-dotted circles illustrate the approximate outline of the Orion-Eridanus superbubble, dust ring, and Barnard’s loop on the plane of the sky, respectively. The light green and orange background colors represent the observations of thermal dust (at 545 GHz) obtained by the Planck Space Observatory and Hα emission (Finkbeiner 2003), respectively. The bright green region shows Orion A.

In the text
thumbnail Fig. 6.

Simplified 3D approximation of the Orion-Eridanus superbubble and its associated bubbles. Orion A is shown as an orange bent cylinder and the Sun is depicted as a yellow circle. The gray ellipsoid and red and green spheres represent the Orion-Eridanus superbubble, Barnard’s loop, and the dust ring, respectively. Due to uncertainties in the models of the Orion-Eridanus superbubble, the dust ring might also be embedded within the superbubble.

In the text
thumbnail Fig. B.1.

Cartoon showing the sensitivity of the reconstructed 3D magnetic field. The Galactic magnetic field, the Orion A cloud, and the reconstructed field are shown with the red vector, bent gray cylinder, and the blue 3D B (1) arrow, respectively. Top panel: Projection onto the plane of the sky. Middle and lower panels: Different viewing angles of the cloud. The black arrow in the middle panel illustrates the direction of rotation to obtain 3D B(2), which remains consistent with the observations and is semi-convex from the observer’s point of view.

In the text
thumbnail Fig. D.1.

Multi-wavelength observations of the Orion-Eridanus superbubble. The Orion-Eridanus superbubble, dust ring, Barnard’s loop, and the λ Orionis ring are depicted as white, green dashed, red dash-dotted, and blue dotted circles, respectively. Top left panel: Thermal dust observations from the Planck Space Observatory. Top middle panel: observations (Finkbeiner 2003). Top right panel: Composite CO survey of Dame et al. (2001). Lower panel: H I observations at velocities of 2.7 km s−1 and 9.1 km s−1 from HI4PI.

In the text
thumbnail Fig. D.2.

Formation of the arc-shaped magnetic field morphology in the Orion A cloud. The large light blue circle represents a bubble within the Orion-Eridanus superbubble, or the Orion-Eridanus superbubble itself (the blue star represents the progenitor(s)). The small dark blue cylinder illustrates the Orion A filamentary molecular cloud. The 3D magnetic configurations cannot be accurately represented in 2D and a more accurate representation of the field lines is shown in Fig. 4. Left panel: Straight B represents the direction of the initial magnetic field prior to interaction with the bubble or event, while the large-scale mild bending of the magnetic field morphology after interaction with the bubble or event is represented by the curved B. Right panel: Formation of the arc-shaped magnetic field morphology around the Orion A molecular cloud due to further interactions, possibly with the dust ring.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.