GW Orionis: a Pre-Main-Sequence Triple with a Warped Disk and a Torn-Apart Ring As Benchmark for Disk Hydrodynamics
Total Page:16
File Type:pdf, Size:1020Kb
main-sequence (PMS) multiple systems provide us with a unique laboratory (for general reviews on PMS binaries GW Orionis: and multiples see for instance Dˆuchene & Kraus 2013 and A pre-main-sequence triple with a Reipurth et al. 2014). For these systems we are able to warped disk and a torn-apart ring as directly measure the 3-dimensional orbits and dynamical benchmark for disk hydrodynamics masses of the perturbing bodies and can image how the disk responds to the perturbation. Stefan Kraus, University of Exeter (email: [email protected]) One system that has the potential to serve as such a “rosetta stone” for hydrodynamic studies, is the PMS triple GW Orionis. This system is one of the brightest and best- studied T Tauri multiple systems. It is located in the λ Orionis star-forming region (388 pc; Kounkel et al. 2017) and has an age of ∼ 1 million years (Calvet et al. 2004). With orbital periods of ∼ 9 months and 11 years, the or- bital periods and separations are just in the right range to enable a full orbit characterisation, while expecting at the same time strong interactions between the disk and the stars. The triple star system GW Ori is long known to be a single-lined spectroscopic binary with a 242 day period (Mathieu et al. 1991). Prato et al. (2018) reported the detection of lines associated with Pre-main-sequence multiples as benchmark the secondary. Observations with the IOTA infrared inter- for disk hydrodynamics ferometer resolved the inner binary and discovered a third component (Berger et al. 2011). Building on the work Understanding how bodies interact with each other and by Mathieu and coworkers, Czekala et al. presented in with disk material holds the key to understanding the 2017 an impressive set of spectra that were obtained over architecture of stellar systems and of planetary systems. 35years with the Fred L. Whipple Observatory and Oak While the interactions between point sources can be de- Ridge Observatory and that provide radial velocities for scribed by simple gravity, interactions with disk material the primary and secondary. Using these radial velocities require further knowledge about viscosity and gas+dust and the IOTA astrometry, Czekala and colleagues derived microphysics that need to be included when simulating first orbit solutions for all 3 stars. These orbit solutions in- disk-body interactions. As a result of our limited knowl- dicated a 11.5 year orbit period for the tertiary and hinted edge in these areas it is, for instance, still difficult to esti- at a signficant misalignment between the stellar orbits and mate the mass of a gap-opening planet from the morphol- the disk, although the small orbital arc covered by the ogy of a gap observed in a protoplanetary disk, or, to de- IOTA astrometry resulted in degeneracies in the orbit fits. arXiv:2012.06578v1 [astro-ph.SR] 5 Oct 2020 rive with certainty whether a gap is opened by a planet in- Between 2008 and 2019, the VLT Interferometer and the stead of by other processes. Furthermore, numerical sim- CHARA array were used to monitor the astrometric orbit ulations continue to unveil new dynamical processes that of the inner binary and the tertiary (Kraus et al. 2020a). might shape protoplanetary disk structures and affect the The resulting orbits are shown in Figure 1 and correspond planet populations forming from these disks. One example to a near-circular (e =0.069±0.009), 241.62±0.05day or- is disk tearing that might occur in the disks around mul- bit for the inner binary and a 4216.8±4.6day orbit for the tiple stars whose orbital angular momentum vectors are tertiary with significant eccentricity (e = 0.379 ± 0.003), misaligned with respect to the rotation axis of the disk. where the mutual inclination between the orbital planes is ◦ Based on computer simulations it has been proposed in 13.9 ± 1.1 . 2012 that the resulting gravitational torques could tear The precise masses of the components in the GWOri sys- the disk apart and cause rings to separate from the disk tem has long been a matter of debate. Mathieu et al. and to precess independently around the central objects (1995) estimated the mass using evolutionary tracks, yield- (Nixon et al. 2012, 2013, Facchini et al. 2013, 2018). In ing 2.5 M⊙ and 0.5 M⊙ for the primary and secondary, order to test such theories and to calibrate the fundamen- respectively. Other workers estimated the mass from the tal parameters involved in hydrodynamic simulations, pre- H-band flux ratio and derived more equal mass ratios (e.g. 1 20 20 2010.0 C D 2009.9 2011.0 4 A B 2009.0 10 10 2005.9 2011.9 2003.9 2009.9 2019.7 [km/s] [km/s] 2011.9 0 0 A,A-B 2 A,AB-C RV 2010.0 C RV 2012.8 -10 -10 2019.7 A A+B 0 -20 -20 2015.0 B E F 2018.1 2015.0 30 30 dDEC[mas] 2003.9 20 20 2011.0 -2 10 10 2005.9 [km/s] [km/s] 0 0 2005.0 2018.1 B,A-B 2009.0 2012.8 2005.0 -10 -10 B,AB-C RV RV -20 -20 -4 GWOriC GWOriB GWOriB -30 -30 4 2 0 -2 -4 30 20 10 0 -10 -20 -30 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 Phase A-B Phase AB-C dRA [mas] dRA [mas] Figure 1: Top: GW Ori orbit fit to the astrometry of the inner binary (A-B; panel A) and the tertiary (AB-C; panel B), where the primary is in the origin of the coordinate sys- tem (Kraus et al. 2020a). Panels C-F show the fit to the radial velocity data presented by Czekala et al. 2017. Left: Aperture-synthesis image of GW Ori, obtained with the MIRC-X instrument and CHARA on 2019 August 27. The white dot marks the location of the center-of-mass of the system, while the blue/brown/white curves give the best-fit orbits derived for the component A/B/C. From Kraus et al. (2020). 3.2 M⊙ and 2.7 M⊙; Prato et al. 2018). These discrepan- To fit the mid-infrared SED, Fang et al. (2014, 2017) de- cies can likely be explained by the non-negligle & time- rived a dust-depleted gap at ∼ 45 au. Further evidence variable contributions from circumbinary/circumtertiary for a truncated or gapped disk structure comes from the dust emission biasing the near-infrared flux ratio (Kraus line profile of the CO fundamental lines, where Najita et al. 2020a). The dynamical masses resulting from the et al. (2003) noted that the line profile exhibits a nar- orbit solution are MA = 2.47 ± 0.33, MB = 1.43 ± 0.18, row+broad emission component and that the line width and MC =1.36 ± 0.28 M⊙ (Kraus et al. 2020a). increases towards the more energetic transitions. The sys- tem also shows signposts of active accretion, in particular Brγ-emission (Folha & Emerson 2001) and strong UV veil- A highly dynamical environment in the in- −7 −1 ing (M˙ acc =3 × 10 M⊙ yr , Calvet et al. 2004). ner few au There has been a long debate on the viewing geometry The spectral energy distribution (SED) of GWOri fea- of the system: Based on measurements of the stellar rota- P . v i tures strong excess emission from mid-infrared to millime- tion period ( =3 25 days), the rotation velocity ( sin = . ± . −1 ter wavelengths, indicating the presence of circumstellar 43 0 2 5kms ) and estimates of the stellar radius, Bou- vier & Bertout 1989 estimated the inclination of the sys- dust. SED modeling suggested the presence of a circumbi- ∼ ◦ nary disk extending from around 1.2 au (Fang et al. 2017), tem to 15 , i.e. close to face-on. Using a similar method, but other observational data, Mathieu et al. 1995 obtained 2.1au (Artymowicz & Ludow 1994) or 3.3au (Mathieu et ◦ al. 1995) outwards. This circumbinary disk was also re- an inclination of 30 . However, photometric observations solved spatially, where the VLTI and CHARA visibilities also reported Algol-like eclipses (e.g. Shevchenko et al. associate 16 ± 2% of the H-band flux with circumbinary 1992, Lamzin et al. 1998, Czekala et al. 2017) that have +2.5 been interpreted as evidence for a nearly edge-on disk ori- material located 2−0.8 au from the inner binary (Kraus et al. 2020a). entation. Variability on longer time scales has been reported as well, 2 ABThermal!dust!emission!1.3mm Thermal!dust!emission!1.3mm R3asym DC DAB R3 R1 R2 N 100!au 20!au E C 12CO! moment! 0! D 12CO! moment! 1! Bi! et! al.! 2020 +! continuum! (contours) +! continuum! (contours) 12CO! mom.! 1 C2 C1 Figure 2: ALMA observations of GW Ori at 0.12” (right, Bi et al. 2020a) and 0.02” resolution (left, panels A-D, Kraus et al. 2020a). The top row shows continuum images, while the bottom row shows 12CO moment maps. including dramatic changes in the near-infrared SED on and seen under intermediate inclination 37 ± 1◦ (Fig. 2A), timescales of ∼ 20 yrs (Fang et al. 2014). This variabil- representing the angular moment vector of the cloud that ity might be linked with a 0.2 mag-amplitude sinusoidal feeds the disk. The Eastern side of the disk is tilted to- variation in the V-band light curve (Czekala et al.