Issue |
A&A
Volume 634, February 2020
|
|
---|---|---|
Article Number | L11 | |
Number of page(s) | 7 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/201937154 | |
Published online | 10 February 2020 |
Letter to the Editor
Spiral arms and instability within the AFGL 4176 mm1 disc
1
School of Physics & Astronomy, E.C. Stoner Building, The University of Leeds, Leeds LS2 9JT, UK
e-mail: k.g.johnston@leeds.ac.uk
2
Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
3
Institute of Astronomy and Astrophysics, Eberhard Karls University Tübingen, Auf der Morgenstelle 10, 72076 Tübingen, Germany
4
Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
5
Moscow Institute of Physics and Technology, 9 Institutskiy Per., Dolgoprudny 141701, Russia
6
Ural Federal University, 19 Mira St., Ekaterinburg 620075, Russia
7
ESO Headquarters, Karl-Schwarzschild-Str. 2, 85748 Garching bei München, Germany
8
Aperio Software, Headingley Enterprise and Arts Centre, Bennett Road, Headingley, Leeds LS6 3HN, UK
Received:
20
November
2019
Accepted:
21
January
2020
We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ∼1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH3CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH3CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.
Key words: accretion / accretion disks / circumstellar matter / stars: formation / stars: massive / techniques: interferometric
© ESO 2020
1. Introduction
Recently, the Atacama Large Millimeter/submillimeter Array (ALMA) has detected spiral structures in the midplanes of both protostellar (e.g. Tobin et al. 2016; Lee et al. 2019) and protoplanetary discs (e.g. Pérez et al. 2016; Huang et al. 2018; Kurtovic et al. 2018) around low-mass (< 2 M⊙) stars. In this Letter, we present observations which show spiral structure within the disc of the forming massive star AFGL 4176 mm1, which lies at a distance of 4.2 kpc (Green & McClure-Griffiths 2011). This result accompanies recent observations of spiral structure in the discs or inner envelopes around other massive stars (e.g. Csengeri et al. 2018; Maud et al. 2019).
We previously presented ALMA observations of the AFGL 4176 region, uncovering one of the best examples of a near-Keplerian disc around an O-type star (Johnston et al. 2015, hereafter J15). Radiative transfer modelling by J15 of these data determined that the best-fitting disc mass and radius were 12 M⊙ and 2000 au, respectively, and the inclination was close to face-on (i = 30°).
Spiral structure in discs around young high-mass stars has been anticipated by a host of simulations (e.g. Krumholz et al. 2007; Kuiper et al. 2011; Harries et al. 2017; Meyer et al. 2018; Ahmadi et al. 2019), which predict near-Keplerian discs with spiral arms, sizes of the order of hundreds to thousands of astronomical units, and masses of several ∼10% of their host stars. In some cases, disc fragments are predicted (e.g. Meyer et al. 2018; Ahmadi et al. 2019; Rosen et al. 2019), which may have also been observed recently for the first time (Ilee et al. 2018).
The spiral structure seen in these simulations stems mainly from their high disc masses, producing discs that are gravitationally unstable. The Toomre parameter Q (Toomre 1964) can be used as a measure of the level of gravitational instability within a disc:
where cs is the sound speed, κ is the epicyclic frequency, G is the gravitational constant, and Σ is the disc surface density. A range of simulations have shown that values of Q < 1.5−1.7 are required for gravitational instability to occur (Durisen et al. 2007, 2008, and references therein).
The importance of gravitational instabilities in high-mass star formation are far-reaching as they likely play a role in observed variability due to accretion bursts (e.g. Caratti o Garatti et al. 2017; Hunter et al. 2017). Additionally, they have the potential to explain the high binary and multiple fraction for high-mass stars (Duchêne & Kraus 2013).
ALMA spectral setup.
2. Observations
We observed AFGL 4176 with the 12 m antenna array of ALMA Cycle 5 under program 2017.1.00700.S (PI: Johnston). The observations were taken in dual polarisation mode in Band 6 at a frequency of ∼250 GHz (1.2 mm), with one pointing centred on 13h43m01.689s −62°08′51.25″ (ICRS). The observations were taken in three different configurations on 3 October 2017, 4 November 2017, and 25 November 2017, with baseline lengths ranging between 118.1 m−15.0 km, 113.0 m−13.9 km, and 92.1 m−8.5 km, and measurements of precipitable water vapour of 1.1, 0.5, and 0.6 mm on each date, respectively. Therefore, the angular resolution for the continuum is ∼30 mas or ∼130 au (d = 4.2 kpc, robust = 0.5 weighting) and the largest angular scale that the observations are sensitive to is ∼1.6″ (or 6700 au, corresponding to a baseline of 92.1 m). The number of antennas for each observation was 42, 43, and 44 (with 30, 41, and 43 having useful data, respectively). The majority of the 12 antennas removed from the first observation were flagged during the pipeline reduction due to poor phase stability because of high wind. The primary beam sizes range between 20.5 and 22.0″. The spectral setup is described in Table 1. The bandpass and absolute flux calibrator was J1427−4206, the phase calibrator was J1337−6509, the check source was J1308−6707, and the pointing calibrator was J1424−6807.
The data reduction was carried out using Common Astronomy Software Applications (CASA, McMullin et al. 2007) version 5.1.1-5 r40000 via the pipeline version 40896 (Pipeline-CASA51-P2-B), and imaging was carried out in CASA version 5.3.0. In addition to the pipeline calibration, we performed phase-only self-calibration to the continuum and applied these solutions to the line data before imaging. The continuum was imaged using multiscale clean and a Briggs weighting parameter of 0.5. The central frequency of the combined continuum image made from the line-free channels is 248.057 GHz (1.21 mm). The line data presented in this Letter (K-ladder transitions of CH3CN) were imaged using multiscale clean, a robust parameter of 1.0, and a channel width of 1 km s−1 to improve sensitivity.
3. Results and discussion
3.1. 1.2 mm continuum
Figure 1 shows the 1.21 mm continuum emission from the AFGL 4176 region. The new observations show that the source AFGL 4176 mm1 presented in J15 now splits into four sources: mm1-main, mm1b, mm1c, and mm1d. The source mm2 is still single and shows an elongation, which extends in the direction away from mm1.
Fig. 1. 1.21 mm continuum emission of the AFGL 4176 region in greyscale and contours (σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 50, 100). The beam is shown in the bottom left corner (0.034″ × 0.028″, PA = −33.5°), and a scalebar is shown in the bottom right corner. The millimetre sources are labelled. |
Most notably, the new observations show the disc source mm1-main has a significant amount of sub-structure. The west side of the disc contains a spiral arm which emanates from the north-west of the central region of the disc, whereas the east side resembles an arc-like structure which curves like a closing parenthesis away from the centre of the disc. Meyer et al. (2018) present simulated ALMA 1.2 mm continuum observations of a radiation-hydrodynamical simulation of a system with similar properties as AFGL 4176 (their Fig. 16). Compared to our Fig. 1, there are several similarities. If the eastern arc in our observations is interpreted as another spiral arm, which has a fragment forming in the southern part of the arc, both discs have an asymmetric two-armed morphology. In addition, the end of the north-eastern spiral arm in the simulated observation resembles the morphology of the eastern arc-like structure, with a peak or possible fragment in the south-east end of the arc. Thus, the morphology of the dust in our observations is consistent with that seen in simulations of gravitationally unstable discs around forming massive stars.
The peak flux of mm1-main is 10.72 ± 0.02 mJy beam−1 at a position of 13:43:01.681 −62:08:51.32 (ICRS), and the integrated flux above 1σ is 42.67 ± 0.03 mJy. The integrated flux for mm1 is slightly smaller than that determined by J15 with a resolution of ∼0.3″ (50 ± 4 mJy), but it is still within 2σ. Using the same assumptions as J15, including T = 190 K, κν = 0.24 cm2 g−1, and a total-mass-to-dust ratio of 154, we find that the disc mass is 6.6 M⊙, which should be compared to the mass of ∼8 M⊙ determined by J15. We note however that full radiative-transfer modelling in J15 determined a disc mass of 12 M⊙.
There is a central bright compact source at the centre of the mm1 disc. We fitted this central source with a 2D Gaussian using imfit in CASA and subtracted it, thus finding a peak flux of 10.4 mJy beam−1, an integrated flux of 13.9 mJy, as well as convolved and deconvolved sizes of 35.6 × 35.5 mas, PA = 55° and 21.7 × 11.6 mas (corresponding to 91 × 49 au), PA = 56°, respectively. Within the positional accuracy of the Australia Telescope Compact Array (ATCA; ∼0.3″), there is a coincident compact continuum source at 1.23 cm measured with ATCA with an integrated flux of 1.33 mJy (Johnston et al. 2019). Therefore, the spectral index between the centimetre and millimetre compact continuum sources is ∼0.9. If the compact sources at both wavelengths are entirely due to ionised gas, then this would indicate there is compact, partially-optically-thick ionised gas associated with the centre of the disc, which may be a jet or a hypercompact HII region. In addition, removal of the compact source from the 1.21 mm emission would reduce the total disc mass from 6.6 to 4.4 M⊙. However, we note that a fraction of the compact millimetre source is likely also due to the rising surface density expected at the centre of the disc (as seen in the simulated observation of Meyer et al. 2018).
By using an ellipse with PA = 60°, we measured the semi-major axis of mm1-main including the spiral arms at 5σ to be ∼0.25″, giving a disc radius of ∼1000 au. We also used the ratio of the semi-major-to-minor axes of the disc to determine an estimate of the disc inclination. The smallest semi-minor axis that is consistent with the shape of the disc is 0.15″ (this is illustrated in Fig. A.1), providing an upper limit to the inclination of 51°. This is consistent with the inclination of 30° found for the best-fitting model of J15.
In Fig. 2 we show the residual image of mm1-main. The residual image has the central compact source removed and is deprojected, assuming a position angle of 60° and an inclination angle of 30°, which we assume in the remainder of this analysis. Figure 3 shows the same data as a function of radius from mm1-main and position angle around the disc. After masking Fig. 3 so that only the western (now right) spiral arm remains, we find the peak pixel and therefore radius for each position angle, which is plotted as a black line in Fig. 3.
Fig. 2. Deprojected 1.21 mm continuum image of AFGL 4176 mm1-main in greyscale and contours, with the central compact source subtracted (contours: σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 30). The rotated and deprojected beam is shown in the bottom left corner. The best fitting logarithmic and Archimedean spiral models (cf. Fig. 3) are plotted as blue dashed and red solid lines on the western (right) spiral, respectively. The same models rotated 180° are also shown. |
Fig. 3. Flux density as a function of radius and position angle. The peak radius for each position angle is plotted as a black line and points for the western spiral arm, and the green dot-dashed and blue dashed lines show the best fits for logarithmic and Archimedean spiral models, respectively. A positional error bar is shown in the top left corner. |
Similar to the analysis that was carried out for the Disk Substructures at High Angular Resolution Project (DSHARP; e.g. Kurtovic et al. 2018; Huang et al. 2018), we fitted the peak positions of the western spiral arm with both logarithmic and Archimedean spirals using the curve_fit function in scipy with the Levenberg-Marquardt algorithm. We assume that the error in determining the peak position in radius is half the major axis of the beam divided by the lowest signal-to-noise ratio of the detected emission in the arm, that is, five (Taylor et al. 1999). We also weighted each datapoint by the spacing in beams between the sampled position angles at each radius to account for correlation within the scale of a beam. The equations for the models, as a function of radius r and position angle θ around the disc, and best fit parameters are given in Table 2. With a lower χ2 of 8.9, the Archimedean spiral model more closely reproduces the data.
Spiral models and associated fit parameters.
Symmetric, tightly-wound logarithmic spirals with a constant pitch angle are expected from simulations of gravitational instabilities (GI) in discs (e.g. Cossins et al. 2009), whereas asymmetric spirals with a varying pitch angle are expected in the case of companion-induced spiral arms (e.g. Rafikov 2002; Forgan et al. 2018). However, when these models are compared using simulated observations, spiral arms in discs produced by a companion can look similar to those predicted for GI (e.g. Meru et al. 2017). Nevertheless, the fact that there is not a clear symmetrical spiral pattern on the other side of the disc suggests that the spiral structure is generated to some degree by the existence of one or more nearby companions. In fact, the other continuum sources in the field seen in Fig. 1 may act as these. These observed companions may have been produced by previous episodes of gravitational instability within the disc and later ejected, or formed from the surrounding cloud.
3.2. CH3CN line emission
Figure 4 presents the first moment map of the CH3CN J = 13−12 K = 3 line emission from AFGL 4176, which is overlaid with contours of 1.21 mm continuum emission. The velocity gradient across the disc, which was originally seen by J15, is apparent in the data. However, further details can now also be seen, such as a velocity gradient across the inner part of the redshifted spiral arm. On the larger scales, the CH3CN emission shows a bridge of gas between mm1-main and mm1c, which extends in a north-west direction from the blueshifted side of the disc. Therefore, the gas emission from the AFGL 4176 region shows evidence of tidal interactions between mm1-main and possible companions, such as mm1c. The direction of this extended north-west arm or bridge mirrors the redshifted spiral on the west side of the disc, indicating mm1c may be exciting the spiral structure within mm1-main. An alternative explanation is that the process of ejection of mm1c from the disc could produce a stretched gas link between the secondary fragment or protostar and the primary.
Fig. 4. CH3CN J = 13−12 K = 3 first moment map in colourscale and 1.21 mm continuum in contours (same as in Fig. 1). The beam is shown in the bottom left corner (0.043″ × 0.038″, PA = −32.6°), and a scalebar is shown in the bottom right corner. |
Similar to the analysis carried out by J15, we created excitation temperature, velocity, linewidth, and column density maps across AFGL 4176 by fitting the CH3CN and CHCN J = 13−12 K-ladder emission associated with each pixel with the Markov chain Monte Carlo χ2 minimisation algorithm within the Centre d’Analyse Scientifique de Spectres Instrumentaux et Synthétiques tool (CASSIS)1. We used the same setup as described in J15 but with slightly updated initial parameter ranges (Tex = 50−350 K, vLSR = −58 – −48 km s−1, vFWHM = 0.5−15 km s−1, and Nmol = 1 × 1015–1 × 1018 cm−2). We determined the excitation temperature Tex and velocity across the disc, as shown in Fig. 5. As we assume local thermodynamic equilibrium (LTE) in our analysis, Tex therefore provides an estimate of the kinetic temperature. Maps of the remaining parameters (linewidth and column density) are presented in Figs. B.1 and B.2. The maps in Fig. 5 are deprojected and are clipped to the 4σ level in the continuum maps (shown as contours).
Fig. 5. Results of CASSIS modelling. Panel a: excitation temperature. Panel b: observed velocity. Both maps have been deprojected similarly to Fig. 2 (we note that the velocities are not yet corrected for inclination). The 1.21 mm continuum emission as shown in Fig. 2 is overplotted as contours in both panels. The maps are masked to values above 4σ in the continuum map. The deprojected beam is shown in the bottom left corner of each panel. |
Figure 5a shows that the temperature is lower towards the arc of continuum emission on the left, but higher in the inner leading section of the western spiral arm on the right, and directly to the upper left of the central source. We note that increased temperatures within spiral arms are seen in many simulations of self-gravitating discs (e.g. Boley & Durisen 2008). Figure 5b shows the velocity gradient across the entire disc as well as a velocity gradient across the inner part of the disc within several hundred astronomical units.
3.3. Resolved Toomre Q map
Similar to the analyses conducted by Ahmadi et al. (2018) and Maud et al. (2019), we used the results from previous sections to construct a map of the Toomre parameter Q across the disc, given in Eq. (1). The sound speed cs can be calculated from the temperature, using the equation
where k is the Boltzmann constant, T is the temperature found from CASSIS, mp is the mass of a proton, and 2.37 is the mean molecular weight per free particle (Kauffmann et al. 2008). The angular velocity Ω was found by subtracting the systematic velocity (vLSR = −52.5 km s−1) from the velocity map determined by CASSIS, deprojecting the velocities, and then dividing by the radius r at each point in the disc to convert from tangential to angular velocities. Here we assume that the disc is close to Keplerian so that κ = Ω; however, we have calculated κ and the resulting map of Q without this assumption in Appendix C. Finally, the surface density Σ can be calculated using the equation
where g is the total-mass-to-dust ratio (154 ± 15%, derived from the ratio of the total mass to mass of hydrogen of 1.402 and the ratio of the mass of dust to mass of hydrogen of 0.0091 from Tables 1.4 and 23.1 of Draine 2011), Speak is the peak flux of the ALMA continuum with the central compact source removed, Θ is the ALMA continuum beam area, B(ν, T) is the black body function, and κν is the opacity (0.24 cm2 g−1 at 1.21 mm for RV = 5.5, Draine 2003a,b). Little is known about the properties of dust around massive forming stars, but as they form quickly, we assume that the dust is typical of star forming clouds, represented by the Draine Milky-Way RV = 5.5 model.
Figure 6 shows the resulting map of Q across the disc. The purple areas in the map show parts of the disc which are Toomre unstable (Q < 1.5−1.7, Durisen et al. 2007). The low Q values in the centre of the disc are due to an artefact caused by the velocity quickly shifting from blue to redshifted. However, further out in the disc we can see that the left (eastern) arc is clearly Toomre unstable, which may already contain a fragment at its lower end, as well as the densest parts of the western spiral on the right. Thus, our analysis shows that this disc has the right conditions to undergo fragmentation that may produce further companions.
Fig. 6. Deprojected map of Q, similar to that shown for Fig. 5. The low Q values in the centre of the disc are an artefact caused by the velocity quickly shifting from blue to redshifted. |
We note that the CH3CN does not trace the midplane, but hotter material higher in the disc atmosphere. In addition, although the disc is not globally optically thick, it may be locally optically thick, which would mean the surface density in the densest parts of the disc is underestimated. These optically thick parts of the disc would have to be unresolved, as the < 1 mJy beam−1 emission from the outer disc has a brightness temperature of < 20 K. Hence by comparison with Fig. 5a, the outer disc is not optically thick at our resolution. From Eqs. (1)–(3) we can see that lowering the temperature and raising the surface density act to lower Q. Therefore, in this case, our analysis provides an upper limit on Q and thus a lower limit on the level of instability of the disc.
In the case the assumed dust opacity κν is a factor of two higher, this would mean that only parts of the disc in Fig. 6 with Q < 0.75 − 0.85 would be unstable. However, these values still occur in the left (eastern) side of the disc.
4. Conclusions
We present high-resolution (∼30 mas) 1.2 mm ALMA observations of the AFGL 4176 disc. We find that the disc is ∼1000 au in radius and contains a spiral arm, which we fitted with logarithmic and Archimedean spiral models. The spiral structure is asymmetric and the observed CH3CN gas emission shows evidence of gas arcs, which link AFGL 4176 mm1-main to other millimetre sources in the field, providing evidence of interactions with companions. We produce a map of the Toomre parameter across the disc and find that the disc and spiral arms have already or are likely to fragment in future, possibly forming companions, thus showing that both current companions and gravitational instabilities play a combined role in creating structure in the AFGL 4176 mm1-main disc. These results provide evidence that instability and fragmentation within discs is a possible pathway to explain the large number of multiple systems observed around high-mass stars.
CASSIS is currently being developed by IRAP-UPS/CNRS (http://cassis.irap.omp.eu).
Acknowledgments
We thank the referee for their comments which improved this manuscript, and John Ilee for useful discussions. AA and HB acknowledge support from the ERC under Horizon 2020 via ERC Consolidator Grant CSF-648505. H. B. also acknowledges support from the DFG in the Collaborative Research Center (SFB 881) “The Milky Way System” (subproject B1). RK acknowledges support via his Emmy Noether Research Group funded by the DFG under grant no. KU 2849/3-1 and KU 2849/3-2. PB acknowledges support by the Russian Science Foundation (grant 18-72-10132).
References
- Ahmadi, A., Beuther, H., Mottram, J. C., et al. 2018, A&A, 618, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ahmadi, A., Kuiper, R., & Beuther, H. 2019, A&A, 632, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Boley, A. C., & Durisen, R. H. 2008, ApJ, 685, 1193 [NASA ADS] [CrossRef] [Google Scholar]
- Caratti o Garatti, A., Stecklum, B., Garcia Lopez, R., et al. 2017, Nat. Phys., 13, 276 [NASA ADS] [CrossRef] [Google Scholar]
- Cossins, P., Lodato, G., & Clarke, C. J. 2009, MNRAS, 393, 1157 [Google Scholar]
- Csengeri, T., Bontemps, S., Wyrowski, F., et al. 2018, A&A, 617, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Draine, B. T. 2003a, ARA&A, 41, 241 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T. 2003b, ApJ, 598, 1017 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium (Princeton: Princeton University Press) [Google Scholar]
- Duchêne, G., & Kraus, A. 2013, ARA&A, 51, 269 [Google Scholar]
- Durisen, R. H., Boss, A. P., Mayer, L., et al. 2007, Protostars and Planets V (Tucson: University of Arizona Press), 607 [Google Scholar]
- Durisen, R. H., Hartquist, T. W., & Pickett, M. K. 2008, Ap&SS, 317, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Forgan, D. H., Ilee, J. D., & Meru, F. 2018, ApJ, 860, L5 [NASA ADS] [CrossRef] [Google Scholar]
- Green, J. A., & McClure-Griffiths, N. M. 2011, MNRAS, 417, 2500 [NASA ADS] [CrossRef] [Google Scholar]
- Harries, T. J., Douglas, T. A., & Ali, A. 2017, MNRAS, 471, 4111 [NASA ADS] [CrossRef] [Google Scholar]
- Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018, ApJ, 869, L43 [NASA ADS] [CrossRef] [Google Scholar]
- Hunter, T. R., Brogan, C. L., MacLeod, G., et al. 2017, ApJ, 837, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Ilee, J. D., Cyganowski, C. J., Brogan, C. L., et al. 2018, ApJ, 869, L24 [NASA ADS] [CrossRef] [Google Scholar]
- Johnston, K. G., Robitaille, T. P., Beuther, H., et al. 2015, ApJ, 813, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Johnston, K. G., Hoare, M. G., Beuther, H., et al. 2019, ApJ, submitted [Google Scholar]
- Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, II, N. J., & Lee, C. W. 2008, A&A, 487, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Krumholz, M. R., Klein, R. I., & McKee, C. F. 2007, ApJ, 665, 478 [NASA ADS] [CrossRef] [Google Scholar]
- Kuiper, R., Klahr, H., Beuther, H., & Henning, T. 2011, ApJ, 732, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Kurtovic, N. T., Pérez, L. M., Benisty, M., et al. 2018, ApJ, 869, L44 [NASA ADS] [CrossRef] [Google Scholar]
- Lee, C. F., Li, Z. Y., & Turner, N. J. 2019, Nat. Astron., 466, 24 [NASA ADS] [CrossRef] [Google Scholar]
- Maud, L. T., Cesaroni, R., Kumar, M. S. N., et al. 2019, A&A, 627, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
- Meru, F., Juhász, A., Ilee, J. D., et al. 2017, ApJ, 839, L24 [NASA ADS] [CrossRef] [Google Scholar]
- Meyer, D. M.-A., Kuiper, R., Kley, W., Johnston, K. G., & Vorobyov, E. 2018, MNRAS, 473, 3615 [NASA ADS] [CrossRef] [Google Scholar]
- Pérez, L. M., Carpenter, J. M., Andrews, S. M., et al. 2016, Science, 353, 1519 [NASA ADS] [CrossRef] [Google Scholar]
- Rafikov, R. R. 2002, ApJ, 569, 997 [NASA ADS] [CrossRef] [Google Scholar]
- Rosen, A. L., Li, P. S., Zhang, Q., & Burkhart, B. 2019, ApJ, 887, 108 [NASA ADS] [CrossRef] [Google Scholar]
- Taylor, G. B., Carilli, C. L., & Perley, R. A. 1999, ASP Conf. Ser., 180, 127 [Google Scholar]
- Tobin, J. J., Kratter, K. M., Persson, M. V., et al. 2016, Nature, 538, 483 [NASA ADS] [CrossRef] [Google Scholar]
- Toomre, A. 1964, ApJ, 139, 1217 [Google Scholar]
Appendix A: Measuring the geometry of the mm1-main disc from continuum emission
Figure A.1 shows a zoom-in of the 1.21 mm continuum emission from mm1-main with ellipses used to measure its inclination over plotted (cf. Sect. 3.1). The red ellipse shows the smallest reasonable semi-minor axis of 0.15″, which corresponds to i = 51°, whereas the blue dashed ellipse shows the geometry expected for i = 30°.
Fig. A.1. 1.21 mm continuum emission of AFGL 4176 mm1-main in greyscale and contours (σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 50, 100). The beam is shown in the bottom left corner (0.034″ × 0.028″, PA = −33.5°), and a scalebar is in the right corner. Two ellipses with semi-major axes of 1000 au or ∼0.24″ are shown, one with a semi-minor axis of 0.15″ (solid red, corresponding to i = 51°) and the other with a semi-minor axis of 0.22″ (dashed blue, corresponding to i = 30°). |
Appendix B: Maps of column density and linewidth from CASSIS modelling
Figures B.1 (column density) and B.2 (linewidth) show the remaining maps resulting from the CASSIS LTE modelling. The column density shown in Fig. B.1 is highest in the leading section of the right or western spiral arm, which is offset from the spiral peak seen in the continuum, as well as in the lower end of the left or eastern arc, towards the continuum peak within the arc. The linewidth shown in Fig. B.2 shows a ring of high values around the position of the subtracted continuum source, which peaks to the lower right.
Appendix C: Calculation of the Toomre parameter without assuming Keplerian rotation
Here we determine κ, the epicyclic frequency without the assumption of Keplerian rotation and use this to determine a new map of Q. The equation relating κ to Ω in this case is
where r is the radius within the disc. Inserting Ω = v/r we get
To determine κ across the disc, we used the velocity map shown in Fig. 5 after subtracting the systematic velocity and deprojection of the velocities. The differential in Eq. (C.2) was determined by finding the radial gradient of rv across the map. The resulting value for κ was inserted into Eq. (1) to determine the Toomre parameter Q. A map of Q for the above assumptions is shown in Fig. C.1.
Fig. C.1. Deprojected map of Q without the assumption that κ = Ω. Undefined areas of the figure are enclosed by a purple contour; these regions are gravitationally unstable. |
In some regions of the disc, the derivative in Eq. (C.2) is negative and therefore κ is undefined, leading to gaps in the map. The physical meaning of a negative κ2 is that radially moving a parcel of gas does not lead to oscillation, but instability. As κ2 tends to zero, velocity shear has a reduced balancing effect against gravity, and as κ2 goes negative, the velocity shear contributes to rather than suppresses instability. Thus, regions where Q cannot be computed are unstable. To indicate this, we have drawn a purple contour around the area in which κ2 was calculated. Our map of Q (Fig. C.1), which was computed using κ, indeed shows in the eastern region that values of Q start to show purple areas indicating instability, before being undefined, indicating strong instability. Our results are therefore consistent with the map derived from the assumption that κ = Ω, which indicates that the eastern (left) region was unstable, although in this Q map the western (right) spiral is stable.
All Tables
All Figures
Fig. 1. 1.21 mm continuum emission of the AFGL 4176 region in greyscale and contours (σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 50, 100). The beam is shown in the bottom left corner (0.034″ × 0.028″, PA = −33.5°), and a scalebar is shown in the bottom right corner. The millimetre sources are labelled. |
|
In the text |
Fig. 2. Deprojected 1.21 mm continuum image of AFGL 4176 mm1-main in greyscale and contours, with the central compact source subtracted (contours: σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 30). The rotated and deprojected beam is shown in the bottom left corner. The best fitting logarithmic and Archimedean spiral models (cf. Fig. 3) are plotted as blue dashed and red solid lines on the western (right) spiral, respectively. The same models rotated 180° are also shown. |
|
In the text |
Fig. 3. Flux density as a function of radius and position angle. The peak radius for each position angle is plotted as a black line and points for the western spiral arm, and the green dot-dashed and blue dashed lines show the best fits for logarithmic and Archimedean spiral models, respectively. A positional error bar is shown in the top left corner. |
|
In the text |
Fig. 4. CH3CN J = 13−12 K = 3 first moment map in colourscale and 1.21 mm continuum in contours (same as in Fig. 1). The beam is shown in the bottom left corner (0.043″ × 0.038″, PA = −32.6°), and a scalebar is shown in the bottom right corner. |
|
In the text |
Fig. 5. Results of CASSIS modelling. Panel a: excitation temperature. Panel b: observed velocity. Both maps have been deprojected similarly to Fig. 2 (we note that the velocities are not yet corrected for inclination). The 1.21 mm continuum emission as shown in Fig. 2 is overplotted as contours in both panels. The maps are masked to values above 4σ in the continuum map. The deprojected beam is shown in the bottom left corner of each panel. |
|
In the text |
Fig. 6. Deprojected map of Q, similar to that shown for Fig. 5. The low Q values in the centre of the disc are an artefact caused by the velocity quickly shifting from blue to redshifted. |
|
In the text |
Fig. A.1. 1.21 mm continuum emission of AFGL 4176 mm1-main in greyscale and contours (σ = 23 μJy beam−1 × −5, 5, 7, 10, 12, 16, 20, 25, 50, 100). The beam is shown in the bottom left corner (0.034″ × 0.028″, PA = −33.5°), and a scalebar is in the right corner. Two ellipses with semi-major axes of 1000 au or ∼0.24″ are shown, one with a semi-minor axis of 0.15″ (solid red, corresponding to i = 51°) and the other with a semi-minor axis of 0.22″ (dashed blue, corresponding to i = 30°). |
|
In the text |
Fig. B.1. Results of CASSIS LTE modelling similar to Fig. 5 but for column density. |
|
In the text |
Fig. B.2. Results of CASSIS LTE modelling similar to Fig. 5 but for linewidth. |
|
In the text |
Fig. C.1. Deprojected map of Q without the assumption that κ = Ω. Undefined areas of the figure are enclosed by a purple contour; these regions are gravitationally unstable. |
|
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.