Here we report the detection of dust depletion in a misaligned inner disk around UX Tau A using JWST MIRI spectra. Midinfrared MIR continuum seesaw variability was detected in this disk by Spitzer and attributed to variable shadows cast on the outer disk by the inner disk. The JWST MIRI spectrum of UX Tau A also shows seesaw variability but with a significant decrease of emission shortward of 10 m to nearly photospheric levels. We argue that UX Tau As MIR continuum variability is due to depletion of dust in a misaligned inner disk. We speculate that this dust depletion occurs because the inner disk is misaligned from the outer disk which can disrupt the replenishment of the inner disk from the outer disk. Using contemporaneous measurements of the mass accretion rate of UX Tau A and estimating the amount of dust necessary to produce the MIR excess in the Spitzer observations we estimate a minimum dust depletion timescale of 0.1 yr. These observations show that we can indirectly detect the signatures of misaligned inner disks through MIR continuum variability and that in some cases the inner disk may be significantly depleted of dust and become optically thin.
['Classical T Tauri stars', 'Protoplanetary disks', '252', '1300', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics']
Evidence for Dust Depletion in a Misaligned Protoplanetary Disk with JWST
{'Evidence for Dust Depletion in a Misaligned Protoplanetary Disk with JWST': "C. C. Espaillat, 1, 2 T. Thanathibodee, 1, 2 Z. Zhu, 3, 4 I. Rabago, 3, 4 J. Wendeborn, 1, 2 N. Calvet, 5 L. Zamudio-Ruvalcaba, 1, 2 M. Volz, 1, 2 C. Pittman, 1, 2 M. McClure, 6 J. F. Babb, 7 R. Franco-Hern'andez, 8 E. Mac'ıas, 9 M. Reynolds, 10, 5 and P. -G. Yan 7  \n1  \nDepartment of Physics and Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Pkwy, Las Vegas, NV 89154, USA 7 Center for Astrophysics | Harvard & Smithsonian, Cambridge, MA 02138, USA  \nDepartment of Astronomy, Boston University, 725 Commonwealth Avenue, Boston, MA 02215, USA 2 Institute for Astrophysical Research, Boston University, 725 Commonwealth Avenue, Boston, MA 02215, USA 3 4 Nevada Center for Astrophysics, University of Nevada, 4505 S. Maryland Pkwy., Las Vegas, NV 89154-4002, USA 5 Department of Astronomy, University of Michigan, 1085 South University Avenue, Ann Arbor, MI 48109, USA 6 Leiden Observatory, Leiden University, PO Box 9513, NL-2300 RA Leiden, The Netherlands  \n8 Instituto de Astronom'ıa y Meteorolog'ıa, Universidad de Guadalajara, Avenida Vallarta No. 2602, Col. Arcos Vallarta Sur, CP 44130, Guadalajara, Jalisco, Mexico  \n9 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748, Garching bei Munchen, Germany 10 Department of Astronomy, Ohio State University, 140 West 18th Ave., Columbus, OH 43210, USA", 'ABSTRACT': "Here we report the detection of dust depletion in a misaligned inner disk around UX Tau A using JWST MIRI spectra. Mid-infrared (MIR) continuum 'seesaw' variability was detected in this disk by Spitzer and attributed to variable shadows cast on the outer disk by the inner disk. The JWST MIRI spectrum of UX Tau A also shows seesaw variability, but with a significant decrease of emission shortwards of 10 µ m to nearly photospheric levels. We argue that UX Tau A's MIR continuum variability is due to depletion of dust in a misaligned inner disk. We speculate that this dust depletion occurs because the inner disk is misaligned from the outer disk, which can disrupt the replenishment of the inner disk from the outer disk. Using contemporaneous measurements of the mass accretion rate of UX Tau A and estimating the amount of dust necessary to produce the MIR excess in the Spitzer observations, we estimate a minimum dust depletion timescale of ∼ 0.1 yr. These observations show that we can indirectly detect the signatures of misaligned inner disks through MIR continuum variability and that in some cases the inner disk may be significantly depleted of dust and become optically thin.  \nKeywords: accretion disks, stars: circumstellar matter, planetary systems: protoplanetary disks, stars: formation, stars: pre-main sequence", '1. INTRODUCTION': "Variability is a distinctive characteristic of low-mass ( < 2 M ⊙ ), accreting pre-main sequence stars (i.e., Classical T Tauri stars; CTTS). Extensive multi-epoch studies of CTTS in the X-ray, UV, optical, and IR wavelengths revealed significant variability at each of these wavelengths (e.g., Stauffer et al. 2015).  \nOut of the wide range of variability seen to date, at least two types are linked to dust in the inner regions of  \nCorresponding author: Catherine C. Espaillat  \ncce@bu.edu  \nthe protoplanetary disk: MIR 'seesaw' continuum variability and optical 'dipper' light curves. Seesaw variability was discovered by Spitzer , which found that the MIR continuum in many disks with large central cavities seesaws, i.e., the flux at shorter wavelengths varies inversely with the flux at longer wavelengths (Muzerolle et al. 2009; Espaillat et al. 2011; Flaherty et al. 2011). Spectral energy distribution (SED) modeling showed that this seesawing behavior could be explained by differences in the size of the shadow cast on the edge of the outer disk by the inner disk. Espaillat et al. (2011) attributed the seesaw behavior to changes in inner disk wall height resulting in differences in illumination in the outer wall, a geometrical effect. NIR imaging of dark re-  \ngions supports shadowing of the outer disk by the inner disk (e.g., Facchini et al. 2018; Zhu 2019; Nealon et al. 2020; Benisty et al. 2023). 'Dippers' were identified in Kepler surveys based on their optical light curves, which display deep dips in their brightness (Cody et al. 2014; Stauffer et al. 2015). The period of the dips is typically a few days and equals the stellar rotation period, which indicates that there is dust located at the corotation radius that obscures the star (Bouvier et al. 2007; Cody & Hillenbrand 2018).  \nHere, we study UX Tau A, which has been reported to display seesaw variability (Espaillat et al. 2011), to further explore the role of dust in the innermost disk in MIR continuum variability and the link to dippers and disk shadows. UX Tau is a quadruple system located in the Taurus-Auriga star-forming region at a distance of ∼ 147 pc (Gaia Collaboration et al. 2018). The UX Tau system is made up of components A, B, and C where A and B are separated by 5.86 '' and A and C are separated by 2.63 '' (White & Ghez 2001). UX Tau B is a close binary separated by 0.1 '' (Duchˆene et al. 2024). In this work, we focus on the A component of the UX Tau system, which has been shown to dominate the emission (White & Ghez 2001).  \nUX Tau A has a stellar mass of 1.4 M ⊙ (Zapata et al. 2020), and reported spectral types of G5-K5 (Rydgren et al. 1976; Espaillat et al. 2010; Kraus & Hillenbrand 2009; Hartigan et al. 1989) with an accretion rate of 1 × 10 -8 M ⊙ yr -1 (Espaillat et al. 2010). This object is surrounded by a pre-transitional disk (i.e., an inner disk separated from an outer disk by a large gap ∼ tens of astronomical units wide; Espaillat et al. 2014). (Sub-)mm imaging finds that the disk has a gap of ∼ 25-31 au (Andrews et al. 2011; Francis et al. 2020) as well as spiral arm structure indicating that UX Tau A has been perturbed by UX Tau C within the past 1000 yrs (M'enard et al. 2020; Zapata et al. 2020). NIR observations detect dust down to at least 23 au (Tanii et al. 2012).  \nHere we present new JWST data of UX Tau A and contemporaneous optical spectra, which we compare to archival MIR spectra and photometry as well as contemporaneous optical light curves. In Section 2, we present the new and supporting datasets. In Section 3, we analyze the optical light curves and H α profiles. We discuss possible connections between the datasets in Section 4 and end with a summary in Section 5.", '2. OBSERVATIONS & DATA REDUCTION': 'We study MIR spectra from JWST and Spitzer , MIR photometry from WISE/NEOWISE, optical light curves from TESS and ASAS-SN, and optical spectra from  \nLCOGT/NRES. In the following, we provide more details about these data.', '2.1. Mid-Infrared Data': "2.1.1. JWST  \nWe present new JWST MIRI (Rieke et al. 2015; Wright et al. 2023) MRS (Wells et al. 2015) data of UX Tau A taken on 2023 October 13 starting at 05:37:37 UT as part of GO program 1676 (PI: Espaillat). The observations were 336 s long and used the default four-point dither. Background observations were also taken with the same exposure time and setup. We follow the same reduction procedure as in Espaillat et al. (2023). In sum, we reduce the uncalibrated raw MRS data using the calibration reference file version jwst 1252.pmap and the JWST Science Calibration Pipeline v.1.15.1 (Bushouse et al. 2024). There is a known flux calibration uncertainty of ∼ 10% (Law et al in prep). The spectrum is presented in Figure 1.  \nThe Channel 1 spectra (4.9-7.5 µ m) are of UX Tau A alone while Channels 2, 3 and 4 (7.5-27.90 µ m) include the C component. At K-band, the flux ratio between A and C is 16.9 ± 1.6 and between A and B it is 3.91 ± 0.31 (White & Ghez 2001). If there was a significant contribution from the C component, there would have been a sharp increase in flux at 7.5 µ m. We conclude that the C component is not contributing significantly to the MIR.  \nUX Tau A has detections of [Ne II ] at 12.81 µ m, [Ne III ] at 15.5 µ m, and [Ar II ] at 6.98 µ m. There are also a few H 2 lines present and weak PAH features at ∼ 6.2, 7.7, 8.6, 12.0, and 12.7 µ m along with a strong PAH feature at 11.3 µ m. We follow Espaillat et al. (2023) to measure [Ne II ], [Ne III ], and [Ar II ] lines fluxes of 1.5 ± 0.1 × 10 -14 , 1.4 ± 0.3 × 10 -15 , and 1.1 ± 0.1 × 10 -14 erg cm -2 s -1 . We note that the [Ne II] line flux is consistent with the line flux measured by Szul'agyi et al. (2012).", '2.1.2. Spitzer': 'In Figure 1, we show Spitzer IRS spectra taken from the Combined Atlas of Sources with Spitzer IRS Spectra (CASSIS, Lebouteiller et al. 2015). These include lowresolution data of UX Tau A on 2008-10-01 and 200810-08. Also included are high-resolution spectra of UX Tau A on 2007-03-22 and 2008-11-05. We note that the spectrum of UX Tau A from 2004-02-08 appeared to have scaling issues with the CASSIS reduction so here we use the spectrum from Espaillat et al. (2011). The systematic absolute flux uncertainties on the IRS spectra are ∼ 5%.  \n2.1.3. WISE/NEOWISE  \nFigure 1. Top: JWST spectrum of UX Tau A with some atomic lines, H 2 lines, and PAH features identified. Bottom: MIR spectra of UX Tau A displaying MIR continuum variability over ∼ 19 yrs (see legend). The JWST MIRI MRS spectrum is shown in black, and the other spectra are from Spitzer . We include WISE/NEOWISE W1 and W2 photometry, which are shown as magenta and teal circles. The stellar photosphere (dotted line) is adapted from Espaillat et al. (2011).  \n<!-- image -->  \nIn Figure 2, we use data from from the Wide-field Infrared Survey Explorer ( WISE ) mission (Wright et al. 2010) and from the NEOWISE project (Mainzer et al. 2011). These data are in WISE bands W1 (3.4 µ m) and W2 (4.6 µ m). The WISE data are from mid-February 2010 and mid-August 2010. The NEOWISE data span 2014 February - 2023 August and are typically obtained in 2-day segments, with 10-20 observations per filter and 130-200 days between segments. The flux calibration uncertainties are less than 5%.', '2.2. Optical Data': '2.2.1. TESS and ASAS-SN  \nWe show All-Sky Automated Survey for Supernovae (ASAS-SN; Kochanek et al. 2017) and TESS light curves for UX Tau A in Figure 3. The ASAS-SN g photometry come from Sky Patrol V2.0 (Hart et al. 2023) and are contemporaneous with the JWST observations. We obtained the TESS light curves using the TESS-GAIA Light Curve (tglc) Python package (Han & Brandt 2023).  \nFigure 2. MIR light curve of UX Tau A from WISE/NEOWISE. The W1 and W2 data points are shown as light magenta and teal circles. The diamonds are the median flux taken over 2-3 days and the error bars are the standard deviation. The vertical broken line notes the time of the JWST observation.  \n<!-- image -->  \nThe TESS observations were simultaneous with JWST . TESS observed UX Tau A for ∼ 26 days in Sector 70 from 2023 September 20 to October 16 (MJD: 60207.9-60233.3) with 3.3-minute cadence. In Figure 3, we include additional TESS light curves of UX Tau A from Sectors 43, 44, and 71 with cadences ranging from 3.3-30 minutes.', '2.2.2. LCOGT': 'We observed UX Tau A using the NRES spectrograph on the Las Cumbres Global Observatory 1-meter telescope network (LCOGT). NRES is a robotic, fiberfed spectrograph providing high-resolution (R ∼ 48000) spectra from the BANZAI-NRES automatic reduction pipeline. Multiple exposures were stacked. All spectroscopic observations are contemporaneous with the JWST observations. UX Tau A was observed starting on 2023-10-11UT00:52, 2023-10-11UT23:50, 202310-14UT09:33, and 2023-10-17UT11:20. All spectra were re-normalized using a polynomial fit to the continuum and are shown in Figure 4.', '3. ANALYSIS & RESULTS': 'We measure the periods of the optical light curves and classify them following Cody et al. (2014). Then we model the H α profiles to measure accretion rates and properties of the accretion flow. These results will be discussed along with the MIR data in Section 4.', '3.1. Optical light curves': "To measure the periods in the TESS data, we subtract a linear fit from each light curve and fit a Lomb-Scargle periodogram. We measure the Q and M variability metrics following Cody & Hillenbrand (2018). Q measures the light curve's periodicity and varies between 0 and 1 where 0 is periodic and 1 is aperiodic. Q values between 0.15 and 0.85 are quasiperiodic. M measures the light curve's asymmetry and typically varies between -1 and 1. Positive values correspond to dips and negative values to bursts with the highest values corresponding to dippers and the lowest values to bursters. M values between ± 0 . 25 are symmetric.  \nThe period of UX Tau A is ∼ 3.8 days amongst the TESS sectors (43, 44, 70, and 71). We find Q values of 0.9, 0.4, 0.9, 0.9 and M values of -0.3, 0.2, 0.5, 1.2. These metrics lead to the following classifications: Sector 43 burster (B); Sector 44 quasiperiodic symmetric (QPS); Sector 70 aperiodic dipper (APD); Sector 71 aperiodic dipper (APD).", '3.2. Optical spectra': 'We model the H α profiles of UX Tau A (Figure 4) using the magnetospheric accretion flow model from Hartmann et al. (1994); Muzerolle et al. (1998, 2001). The magnetic, stellar, and disk rotation axes are aligned and the flow geometry follows a dipolar magnetic field which has the following parameters: inner radius (R i ), width  \nFigure 3. TESS (black) and ASAS-SN g-band (purple) light curves of UX Tau A. To facilitate comparison, the TESS data were scaled to match the ASAS-SN data. A vertical broken line marks the time of the JWST observation.  \n<!-- image -->  \nFigure 4. H α profiles (blue) of UX Tau A and the average accretion flow model fit (black line) along with the top 100 bestfitting models (gray lines). The horizontal and vertical dotted lines correspond to the continuum and line center, respectively.  \n<!-- image -->  \nTable 1. Results of the Magnetospheric Accretion Model  \nof the flow (W r ) at the disk plane, maximum temperature in the flow (T max ), and viewing inclination ( i ). To determine the best fits, we follow the procedure of Thanathibodee et al. (2023).  \nWe created a grid of 29,700 models with the following values: ˙ M (1 × 10 -9 -1 × 10 -7 M ⊙ yr -1 ), R i (1 . 5 -7 . 0 R ∗ ), W r (0 . 2 -1 . 8 R ∗ ), T max (8 -10 × 10 3 K), and i (20 · -60 · ). We calculate the χ 2 for each combination of the model and observed profile and we selected the models where the normalized likelihood is ≥ 0 . 5 and calculated the weighted mean of ˙ M , R i , W r , T max , and i . The best fit parameters are listed in Table 1.', '4. DISCUSSION': 'Seesaw variability is still present in the JWST spectrum of UX Tau A, but differs significantly from what has been seen previously. Strikingly, the JWST spectrum is nearly consistent with photospheric emission at the shortest wavelengths, with a very small MIR excess shortwards of 10 µ m and a broad, weak 10 µ m silicate emission feature. 1 This combination points to a small amount of ∼ micron-sized dust grains in an optically thin inner disk. Previous work has reproduced the 10 µ mfeature of UX Tau A using a grain size distribution of a -3 . 5 , where a varies between a min =0.005 µ m and a max =10 µ m (Espaillat et al. 2011).  \nWe propose that the drop to nearly photospheric levels in UX Tau A is due to dust depletion in a misaligned inner disk undergoing disk breaking. Disk breaking is a phenomenon that occurs when the inner and outer disks precess independently, possibly due to the presence of a companion (Zhu 2019). The inclination of the outer disk of UX Tau A is 37.96 · +0 . 97 -0 . 90 while the inclination of the inner disk is 73.46 · +11 . 76 -16 . 47 (Bohn et al. 2022).  \nAs the inner disk precesses, there will be times when it is aligned with the outer disk (Figure 5, left) and times when it is misaligned (Figure 5, right). Mass transfer from the outer disk to the inner disk is easiest when the disks are coplanar and is the most disrupted when the inner disk is at maximum misalignment from the outer disk. However, accretion may still occur at intermediate points. Because it only takes a small amount of dust to make the disk optically thick, any accretion can deposit sufficient dust to lead to an optically thick inner disk. Therefore, we expect that a significant MIR excess  \nwould be seen most of the time. Then when the accretion to the inner disk is significantly disrupted, most of the dust in the inner disk may accrete or drift to the star due to lack of replenishment from the outer disk. The inner disk becomes optically thin and there will be a drop at the shorter wavelengths in the JWST spectrum, while there is an increase at the longer wavelengths as there is little/no shadow cast on the outer disk. Below we estimate the dust depletion and precession timescales and consider this proposed scenario in light of the MIR light curve, optical light curves, and the MIR emission lines and features.', '4.1. Depletion and precession timescales': "We can estimate the dust depletion timescale by calculating how much mass was in the inner disk at the time of the Spitzer observations and comparing this to the accretion rate of UX Tau A. Assuming the excess MIR emission comes from a marginally optically thick inner disk at 1550 K, the lower limit of the dust mass in the inner disk is then ∼ 10 -11 M ⊙ using a dust opacity of 200 cm 2 g -1 (Birnstiel et al. 2018). We measure accretion rates of 1-2 × 10 -8 M ⊙ yr -1 contemporaneous with the JWST observations (Table 1), which would quickly deplete the inner disk of dust in 0.1 yrs. Here we adopt the typically assumed dust-to-gas ratio of 0.01, but the dust-to-gas ratio in circumstellar disks is unknown and it is not clear what it would be in the inner disk of UX Tau A.  \nWe can also consider the precession timescale of the inner disk of UX Tau A. Theoretical works show that a misaligned Jupiter-mass planet can drive misalignment and precession of the inner disk at ∼ 1000 times the planetary orbital timescale (e.g., Zhu 2019). We assume that the period measured from the TESS light curves (Figure 3) is the stellar rotation period. Then, if a misaligned Jupiter-mass planet is around the corotation radius, the inner disk's precession timescale is ∼ 500 wks or 10 yrs. If the planet is 10 Jupiter masses, the precession timescale is then 1 year.", '4.2. MIR light curve': 'There are no significant dips in the WISE/NEOWISE fluxes. Our estimated dust depletion timescale is ∼ 0.1 yrs ( ∼ 37 days) and the dust replenishment timescale would be similar. Since there are at least 130 days in between the segments of the WISE/NEOWISE data, it is possible that we missed this short-lived depleted period. If the precession timescale is 10 yrs, the last time that the inner disk was in the same position as in the JWST observation would have been in 2013, which falls in the ∼ 3.5 yrs gap in coverage between WISE and NEOWISE. It is plausible that given the sparse time cadence  \n<!-- image -->  \nFigure 5. Schematics showing that the inner and outer disks are coplanar (left) and misaligned (right). If the inner disk has nodal precession along the vertical direction, it changes from the left to right panel after half the precession cycle.  \n<!-- image -->  \nof the data and the short depletion and replenishment timescale, that another significant MIR dip was not observed.  \nThere is significant MIR emission in the last NEOWISE observation taken ∼ 43 days before the JWST observations, which implies a quick drop to nearly photospheric levels. This is still consistent with our estimated depletion timescale. During the last NEOWISE observations, the inner disk could have been marginally optically thick and then it became optically thin when observed with JWST . However, the estimated depletion timescale is a lower limit and any slightly longer timescale would not be consistent with the observations.', '4.3. Optical light curves': "One may speculate that the drop in the JWST spectrum is due to an edge-on inner disk. However, we can exclude the possibility that the inner disk of UX Tau A was close to edge-on at the time of the JWST observations since there is no evidence of extinction of the star by the inner disk in the simultaneous optical light curve (Figure 3).  \nThe dips in the TESS light curve suggest that there is dust in the inner disk along our line of sight which obscures the star a few days before and after the JWST observations. This is consistent with the JWST spectrum whose small MIR excess shows that there is still some small amount of dust grains in the inner disk dust. The exact mechanism by which this dust obscures the star in the UX Tau A system is unclear. Scenarios that have been proposed to explain dippers include dusty ac-  \ncretion flows (Nagel et al. 2024), disk warps (McGinnis et al. 2015), and dusty disk winds (Gaidos et al. 2024).  \nFinally, we see a change in UX Tau A's optical light curve behavior that could possibly be further evidence of a precessing misaligned inner disk. In 2021 (top, Figure 3), UX Tau A had a quasiperiodic symmetric light curve and two years later (bottom, Figure 3) showed a dipper light curve. This could indicate that the inner disk was at different inclinations in 2021 and 2023, changing the obscuration of the star and leading to the different light curve behavior. This may have implications on a subset of objects with 'hybrid' optical light curves which switch from quasiperiodic to dipper (McGinnis et al. 2015; Cody & Hillenbrand 2018).", '4.4. MIR atomic gas emission lines and PAH features': "The proposed scenario of a depleted inner disk also appears roughly consistent with the lines and features seen in the JWST spectrum. First, only the 11.3 µ m PAH feature was clearly seen in the Spitzer spectrum (Watson et al. 2009) whereas they are all prominent in the JWST spectrum. The PAH features were likely 'drowned out' in the Spitzer spectrum due to the significant MIR excess. In the JWST spectrum, the depletion of dust from the inner disk leads to a decrease in the continuum emission, allowing the PAHs to be conspicuous. Likewise, it is easier to see the [Ar II] and [Ne II] emission lines due to the decreased MIR continuum, as has been seen in other pre-transitional and transitional disks (Espaillat et al. 2007, 2013). In addition, the presence of the [Ar II] and [Ne II] lines confirms that there is gas in the inner disk, consistent with the measured accretion rates  \nin this work, which points to a significant gas reservoir in the inner disk despite the depletion of dust.", '5. SUMMARY & CONCLUSIONS': 'We presented a new JWST spectrum of UX Tau A which shows a striking decrease of continuum emission shortwards of 10 µ m, dropping to nearly photospheric levels. We conclude that the dust in the inner disk has been depleted to the point where the inner disk is optically thin, leading to this drop in MIR emission. This dust depletion could be due to a disruption in the replenishment of the inner disk since accretion from the outer disk may be more difficult when the inner and outer disks are misaligned.  \nMonitoring of the precession timescale via inclination measurements, MIR spectra, and/or imaging is necessary to test if this variability is periodic, which would be expected if it is due to disk precession. More detailed  \nmodeling with non-axisymmetric disks is also necessary. These approaches will help us link the variability seen in protoplanetary disks at different wavelengths.', 'ACKNOWLEDGMENTS': 'Weacknowledge support from JWST grant GO-01676. We thank the anonymous referee for a constructive report which improved the paper.  \nThe JWST data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute. The specific observations analyzed can be accessed via 10.17909/sej8-g173.  \nFacilities: ASAS-SN, JWST, LCOGT, WISE, NEOWISE, TESS', 'REFERENCES': "Kraus, A. L., & Hillenbrand, L. A. 2009, ApJ, 704, 531, doi: 10.1088/0004-637X/704/1/531  \nLebouteiller, V., Barry, D. J., Goes, C., et al. 2015, ApJS, 218, 21, doi: 10.1088/0067-0049/218/2/21  \nMainzer, A., Bauer, J., Grav, T., et al. 2011, ApJ, 731, 53, doi: 10.1088/0004-637X/731/1/53  \n- McGinnis, P. T., Alencar, S. H. P., Guimar˜aes, M. M., et al. 2015, A&A, 577, A11, doi: 10.1051/0004-6361/201425475  \nM'enard, F., Cuello, N., Ginski, C., et al. 2020, A&A, 639,  \n- L1, doi: 10.1051/0004-6361/202038356\n- Muzerolle, J., Calvet, N., & Hartmann, L. 1998, ApJ, 492, 743, doi: 10.1086/305069\n- -. 2001, ApJ, 550, 944, doi: 10.1086/319779\n- Muzerolle, J., Flaherty, K., Balog, Z., et al. 2009, ApJL, 704, L15, doi: 10.1088/0004-637X/704/1/L15\n- Nagel, E., Bouvier, J., & Duarte, A. E. 2024, A&A, 688, A61, doi: 10.1051/0004-6361/202450085\n- Nealon, R., Price, D. J., & Pinte, C. 2020, MNRAS, 493, L143, doi: 10.1093/mnrasl/slaa026  \nRieke, G. H., Wright, G. S., Boker, T., et al. 2015, PASP, 127, 584, doi: 10.1086/682252  \nRydgren, A. E., Strom, S. E., & Strom, K. M. 1976, ApJS, 30, 307, doi: 10.1086/190364  \n- Stauffer, J., Cody, A. M., McGinnis, P., et al. 2015, AJ, 149, 130, doi: 10.1088/0004-6256/149/4/130\n- Szul'agyi, J., Pascucci, I., ' Abrah'am, P., et al. 2012, ApJ, 759, 47, doi: 10.1088/0004-637X/759/1/47\n- Tanii, R., Itoh, Y., Kudo, T., et al. 2012, PASJ, 64, 124. https://arxiv.org/abs/1206.1215\n- Thanathibodee, T., Molina, B., Serna, J., et al. 2023, ApJ, 944, 90, doi: 10.3847/1538-4357/acac84\n- Watson, D. M., Leisenring, J. M., Furlan, E., et al. 2009, ApJS, 180, 84, doi: 10.1088/0067-0049/180/1/84\n- Wells, M., Pel, J. W., Glasse, A., et al. 2015, PASP, 127, 646, doi: 10.1086/682281\n- White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265, doi: 10.1086/321542\n- Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868 Wright, G. S., Rieke, G. H., Glasse, A., et al. 2023, PASP, 135, 048003, doi: 10.1088/1538-3873/acbe66\n- Zapata, L. A., Rodr'ıguez, L. F., Fern'andez-L'opez, M., et al. 2020, ApJ, 896, 132, doi: 10.3847/1538-4357/ab8fac Zhu, Z. 2019, MNRAS, 483, 4221,  \ndoi: 10.1093/mnras/sty3358"}