-
A Multi-wavelength, Multi-epoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. III. Optical Spectra
Authors:
John Wendeborn,
Catherine C. Espaillat,
Thanawuth Thanathibodee,
Connor E. Robinson,
Caeley V. Pittman,
Nuria Calvet,
James Muzerolle,
Fredrick M. Walter,
Jochen Eisloffel,
Eleonora Fiorellino,
Carlo F. Manara,
Agnes Kospal,
Peter Abraham,
Rik Claes,
Elisabetta Rigliaco,
Laura Venuti,
Justyn Campbell-White,
Pauline McGinnis,
Manuele Gangi,
Karina Mauco,
Filipe Gameiro,
Antonio Frasca,
Zhen Guo
Abstract:
Classical T Tauri Stars (CTTSs) are highly variable stars that possess gas- and dust-rich disks from which planets form. Much of their variability is driven by mass accretion from the surrounding disk, a process that is still not entirely understood. A multi-epoch optical spectral monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was conducted along with contemporaneous HST UV…
▽ More
Classical T Tauri Stars (CTTSs) are highly variable stars that possess gas- and dust-rich disks from which planets form. Much of their variability is driven by mass accretion from the surrounding disk, a process that is still not entirely understood. A multi-epoch optical spectral monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was conducted along with contemporaneous HST UV spectra and ground-based photometry in an effort to determine accretion characteristics and gauge variability in this sample. Using an accretion flow model, we find that the magnetospheric truncation radius varies between 2.5-5 R* across all of our observations. There is also significant variability in all emission lines studied, particularly Halpha, Hbeta, and Hgamma. Using previously established relationships between line luminosity and accretion, we find that, on average, most lines reproduce accretion rates consistent with accretion shock modeling of HST spectra to within 0.5 dex. Looking at individual contemporaneous observations, however, these relationships are less accurate, suggesting that variability trends differ from the trends of the population and that these empirical relationships should be used with caution in studies of variability.
△ Less
Submitted 8 August, 2024;
originally announced August 2024.
-
A Multi-wavelength, Multi-epoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. II. Photometric Light Curves
Authors:
John Wendeborn,
Catherine C. Espaillat,
Thanawuth Thanathibodee,
Connor E. Robinson,
Caeley V. Pittman,
Nuria Calvet,
Ágnes Kóspál,
Konstantin N. Grankin,
Fredrick M. Walter,
Zhen Guo,
Jochen Eislöffel
Abstract:
Classical T Tauri Stars (CTTSs) are young, low-mass stars which accrete material from their surrounding protoplanetary disk. To better understand accretion variability, we conducted a multi-epoch, multi-wavelength photometric monitoring campaign of four CTTSs: TW Hya, RU Lup, BP Tau, and GM Aur, in 2021 and 2022, contemporaneous with HST UV and optical spectra. We find that all four targets displa…
▽ More
Classical T Tauri Stars (CTTSs) are young, low-mass stars which accrete material from their surrounding protoplanetary disk. To better understand accretion variability, we conducted a multi-epoch, multi-wavelength photometric monitoring campaign of four CTTSs: TW Hya, RU Lup, BP Tau, and GM Aur, in 2021 and 2022, contemporaneous with HST UV and optical spectra. We find that all four targets display significant variability in their light curves, generally on days-long timescales (but in some cases year-to-year) often due to periodicity associated with stellar rotation and to stochastic accretion variability. Their is a strong connection between mass accretion and photometric variability in all bands, but the relationship varies per target and epoch. Thus, photometry should be used with caution as a direct measure of accretion in CTTSs.
△ Less
Submitted 31 May, 2024;
originally announced May 2024.
-
A Multi-wavelength, Multi-epoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. I. HST FUV and NUV Spectra
Authors:
John Wendeborn,
Catherine C. Espaillat,
Sophia Lopez,
Thanawuth Thanathibodee,
Connor E. Robinson,
Caeley V. Pittman,
Nuria Calvet,
Nicole Flors,
Fredrick M. Walter,
Ágnes Kóspál,
Konstantin N. Grankin,
Ignacio Mendigutía,
Hans Moritz Günther,
Jochen Eislöffel,
Zhen Guo,
Kevin France,
Eleonora Fiorellino,
William J. Fischer,
Péter Ábrahám,
Gregory J. Herczeg
Abstract:
The Classical T Tauri Star (CTTS) stage is a critical phase of the star and planet formation process. In an effort to better understand the mass accretion process, which can dictate further stellar evolution and planet formation, a multi-epoch, multi-wavelength photometric and spectroscopic monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was carried out in 2021 and 2022/2023…
▽ More
The Classical T Tauri Star (CTTS) stage is a critical phase of the star and planet formation process. In an effort to better understand the mass accretion process, which can dictate further stellar evolution and planet formation, a multi-epoch, multi-wavelength photometric and spectroscopic monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was carried out in 2021 and 2022/2023 as part of the Outflows and Disks Around Young Stars: Synergies for the Exploration of ULYSSES Spectra (ODYSSEUS) program. Here we focus on the HST UV spectra obtained by the HST Director's Discretionary Time UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. Using accretion shock modeling, we find that all targets exhibit accretion variability, varying from short increases in accretion rate by up to a factor of 3 within 48 hours, to longer decreases in accretion rate by a factor of 2.5 over the course of 1 year. This is despite the generally consistent accretion morphology within each target. Additionally, we test empirical relationships between accretion rate and UV luminosity and find stark differences, showing that these relationships should not be used to estimate the accretion rate for individual target. Our work reinforces that future multi-epoch and simultaneous multi-wavelength studies are critical in our understanding of the accretion process in low-mass star formation.
△ Less
Submitted 31 May, 2024;
originally announced May 2024.
-
Twenty-Five Years of Accretion onto the Classical T Tauri Star TW Hya
Authors:
Gregory J. Herczeg,
Yuguang Chen,
Jean-Francois Donati,
Andrea K. Dupree,
Frederick M. Walter,
Lynne A. Hillenbrand,
Christopher M. Johns-Krull,
Carlo F. Manara,
Hans Moritz Guenther,
Min Fang,
P. Christian Schneider,
Jeff A. Valenti,
Silvia H. P. Alencar,
Laura Venuti,
Juan Manuel Alcala,
Antonio Frasca,
Nicole Arulanantham,
Jeffrey L. Linsky,
Jerome Bouvier,
Nancy S. Brickhouse,
Nuria Calvet,
Catherine C. Espaillat,
Justyn Campbell-White,
John M. Carpenter,
Seok-Jun Chang
, et al. (17 additional authors not shown)
Abstract:
Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veili…
▽ More
Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of $2.51\times10^{-9}$~M$_\odot$~yr$^{-1}$ and a Gaussian distribution with a FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Ly$α$. The accretion luminosities are well correlated with He line luminosities but poorly correlated with H$α$ and H$β$ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars.
△ Less
Submitted 28 August, 2023;
originally announced August 2023.
-
Towards a comprehensive view of accretion, inner disks, and extinction in classical T Tauri stars: an ODYSSEUS study of the Orion OB1b association
Authors:
Caeley V. Pittman,
Catherine C. Espaillat,
Connor E. Robinson,
Thanawuth Thanathibodee,
Nuria Calvet,
John Wendeborn,
Jesus Hernández,
Carlo F. Manara,
Fred Walter,
Péter Ábrahám,
Juan M. Alcalá,
Sílvia H. P. Alencar,
Nicole Arulanantham,
Sylvie Cabrit,
Jochen Eislöffel,
Eleonora Fiorellino,
Kevin France,
Manuele Gangi,
Konstantin Grankin,
Gregory J. Herczeg,
Ágnes Kóspál,
Ignacio Mendigutía,
Javier Serna,
Laura Venuti
Abstract:
The coevolution of T Tauri stars and their surrounding protoplanetary disks dictates the timescales of planet formation. In this paper, we present magnetospheric accretion and inner disk wall model fits to NUV-NIR spectra of nine classical T Tauri stars in Orion OB1b as part of the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey. Using NUV-…
▽ More
The coevolution of T Tauri stars and their surrounding protoplanetary disks dictates the timescales of planet formation. In this paper, we present magnetospheric accretion and inner disk wall model fits to NUV-NIR spectra of nine classical T Tauri stars in Orion OB1b as part of the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey. Using NUV-optical spectra from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program and optical-NIR spectra from the PENELLOPE VLT Large Programme, we find that the accretion rates of these targets are relatively high for the region's intermediate age of 5.0 Myr; rates range from $0.5-17.2 \times 10^{-8}$ M$_{\odot}$/yr, with a median value of $1.2\times 10^{-8}$ M$_{\odot}$/yr. The NIR excesses can be fit with 1200-1800 K inner disk walls located at 0.05-0.10 AU from the host stars. We discuss the significance of the choice in extinction law, as the measured accretion rate depends strongly on the adopted extinction value. This analysis will be extended to the complete sample of T Tauri stars being observed through ULLYSES to characterize accretion and inner disks in star-forming regions of different ages and stellar populations.
△ Less
Submitted 29 August, 2022; v1 submitted 9 August, 2022;
originally announced August 2022.
-
Understanding Accretion Variability Through TESS Observations of Taurus
Authors:
Connor E. Robinson,
Catherine C. Espaillat,
Joseph E. Rodriguez
Abstract:
Interpreting the short-timescale variability of the accreting, young, low-mass stars known as Classical T Tauri stars remains an open task. Month-long, continuous light curves from the Transiting Exoplanet Survey Satellite (\textit{TESS}) have become available for hundreds of T Tauri stars. With this vast data set, identifying connections between the variability observed by \TESS and short-timesca…
▽ More
Interpreting the short-timescale variability of the accreting, young, low-mass stars known as Classical T Tauri stars remains an open task. Month-long, continuous light curves from the Transiting Exoplanet Survey Satellite (\textit{TESS}) have become available for hundreds of T Tauri stars. With this vast data set, identifying connections between the variability observed by \TESS and short-timescale accretion variability is valuable for characterizing the accretion process. To this end, we obtained short-cadence \TESS observations of 14 T Tauri stars in the Taurus star-formation region along with simultaneous ground-based, UBVRI-band photometry to be used as accretion diagnostics. In addition, we combine our dataset with previously published simultaneous NUV-NIR \textit{Hubble Space Telescope} spectra for one member of the sample. We find evidence that much of the short-timescale variability observed in the \TESS light curves can be attributed to changes in the accretion rate, but note significant scatter between separate nights and objects. We identify hints of time lags within our dataset that increase at shorter wavelengths which we suggest may be evidence of longitudinal density stratification of the accretion column. Our results highlight that contemporaneous, multi-wavelength observations remain critical for providing context for the observed variability of these stars.
△ Less
Submitted 30 June, 2022;
originally announced July 2022.
-
Measuring the Density Structure of an Accretion Hot Spot
Authors:
C. C. Espaillat,
C. E. Robinson,
M. M. Romanova,
T. Thanathibodee,
J. Wendeborn,
N. Calvet,
M. Reynolds,
J. Muzerolle
Abstract:
Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto the star via funnel flows which follow the stellar field lines and shock on the stellar surface leaving a hot spot with a density gradient. Previous work has inferred different densities in the hot spot, but has not been sensitive to the radial density distribution. Attempts have been made to measure this w…
▽ More
Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto the star via funnel flows which follow the stellar field lines and shock on the stellar surface leaving a hot spot with a density gradient. Previous work has inferred different densities in the hot spot, but has not been sensitive to the radial density distribution. Attempts have been made to measure this with X-ray observations, but X-ray emission only traces a fraction of the hot spot and also coronal emission. Here we report periodic ultraviolet and optical light curves of the accreting star GM Aur that display a time lag of about 1 day between their peaks. The periodicity arises as the source of the ultraviolet and optical emission moves into and out of view as it rotates along with the star. The time-lag indicates a difference in the spatial distribution of ultraviolet and optical brightness over the stellar surface. Within the framework of the magnetospheric accretion model, this indicates a radial density gradient in a hot spot on the stellar surface since different density parts of the hot spot are expected to emit radiation at different wavelengths. These results are the first observational confirmation of the magnetospheric accretion model's prediction of a density gradient in the hot spot and demonstrate the insights gained from focusing on the wavelengths where the bulk of the accretion energy can be observed.
△ Less
Submitted 1 September, 2021;
originally announced September 2021.
-
Synthetic Light Curves of Accretion Variability in T Tauri Stars
Authors:
Connor E. Robinson,
Catherine C. Espaillat,
James E Owen
Abstract:
Photometric observations of accreting, low-mass, pre-main-sequence stars (i.e., Classical T Tauri stars; CTTS) have revealed different categories of variability. Several of these classifications have been linked to changes in $\dot{M}$. To test how accretion variability conditions lead to different light-curve morphologies, we used 1D hydrodynamic simulations of accretion along a magnetic field li…
▽ More
Photometric observations of accreting, low-mass, pre-main-sequence stars (i.e., Classical T Tauri stars; CTTS) have revealed different categories of variability. Several of these classifications have been linked to changes in $\dot{M}$. To test how accretion variability conditions lead to different light-curve morphologies, we used 1D hydrodynamic simulations of accretion along a magnetic field line coupled with radiative transfer models and a simple treatment of rotation to generate synthetic light curves. We adopted previously developed metrics in order to classify observations to facilitate comparisons between observations and our models. We found that stellar mass, magnetic field geometry, corotation radius, inclination, and turbulence all play roles in producing the observed light curves and that no single parameter is entirely dominant in controlling the observed variability. While the periodic behavior of the light curve is most strongly affected by the inclination, it is also a function of the magnetic field geometry and inner disk turbulence. Objects with either pure dipole fields, strong aligned octupole components, or high turbulence in the inner disk all tend to display accretion bursts. Objects with anti-aligned octupole components or aligned, weaker octupole components tend to show light curves with slightly fewer bursts. We did not find clear monotonic trends between the stellar mass and empirical classification. This work establishes the groundwork for more detailed characterization of well-studied targets as more light curves of CTTS become available through missions such as the Transiting Exoplanet Survey Satellite (TESS).
△ Less
Submitted 12 January, 2021;
originally announced January 2021.
-
Probing CO and N$_2$ Snow Surfaces in Protoplanetary Disks with N$_2$H$^+$ Emission
Authors:
Chunhua Qi,
Karin I. Öberg,
Catherine C. Espaillat,
Connor E. Robinson,
Sean M. Andrews,
David J. Wilner,
Geoffrey A. Blake,
Edwin A. Bergin,
L. Ilsedore Cleeves
Abstract:
Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of b…
▽ More
Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$''$-0.4$''$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks.
△ Less
Submitted 24 July, 2019;
originally announced July 2019.
-
Herschel Observations of Protoplanetary Disks in Lynds 1641
Authors:
Sierra L. Grant,
Catherine C. Espaillat,
S. Thomas Megeath,
Nuria Calvet,
William J. Fischer,
Christopher J. Miller,
Kyoung Hee Kim,
Amelia M. Stutz,
Álvaro Ribas,
Connor E. Robinson
Abstract:
We analyze Herschel Space Observatory observations of 104 young stellar objects with protoplanetary disks in the ~1.5 Myr star-forming region Lynds 1641 (L1641) within the Orion A Molecular Cloud. We present spectral energy distributions from the optical to the far-infrared including new photometry from the Herschel Photodetector Array Camera and Spectrometer (PACS) at 70 microns. Our sample, take…
▽ More
We analyze Herschel Space Observatory observations of 104 young stellar objects with protoplanetary disks in the ~1.5 Myr star-forming region Lynds 1641 (L1641) within the Orion A Molecular Cloud. We present spectral energy distributions from the optical to the far-infrared including new photometry from the Herschel Photodetector Array Camera and Spectrometer (PACS) at 70 microns. Our sample, taken as part of the Herschel Orion Protostar Survey, contains 24 transitional disks, eight of which we identify for the first time in this work. We analyze the full disks with irradiated accretion disk models to infer dust settling properties. Using forward modeling to reproduce the observed nKS-[70] index for the full disk sample, we find the observed disk indices are consistent with models that have depletion of dust in the upper layers of the disk relative to the mid plane, indicating significant dust settling. We perform the same analysis on full disks in Taurus with Herschel data and find that Taurus is slightly more evolved, although both samples show signs of dust settling. These results add to the growing literature that significant dust evolution can occur in disks by ~1.5 Myr.
△ Less
Submitted 18 June, 2018;
originally announced June 2018.
-
Periodic Eclipses of the Young Star PDS 110 Discovered with WASP and KELT Photometry
Authors:
H. P. Osborn,
J. E. Rodriguez,
M. A. Kenworthy,
G. M. Kennedy,
E. E. Mamajek,
C. E. Robinson,
C. C. Espaillat,
D. J. Armstrong,
B. J. Shappee,
A. Bieryla,
D. W. Latham,
D. R. Anderson,
T. G. Beatty,
P. Berlind,
M. L. Calkins,
G. A. Esquerdo,
B. S. Gaudi,
C. Hellier,
T. W. -S. Holoien,
D. James,
C. S. Kochanek,
R. B. Kuhn,
M. B. Lund,
J. Pepper,
D. L. Pollacco
, et al. (6 additional authors not shown)
Abstract:
We report the discovery of eclipses by circumstellar disc material associated with the young star PDS 110 in the Ori OB1a association using the SuperWASP and KELT surveys. PDS 110 (HD 290380, IRAS 05209-0107) is a rare Fe/Ge-type star, a ~10 Myr-old accreting intermediate-mass star showing strong infrared excess (L$_{\rm IR}$/L$_{\rm bol}$ ~ 0.25). Two extremely similar eclipses with a depth of ~3…
▽ More
We report the discovery of eclipses by circumstellar disc material associated with the young star PDS 110 in the Ori OB1a association using the SuperWASP and KELT surveys. PDS 110 (HD 290380, IRAS 05209-0107) is a rare Fe/Ge-type star, a ~10 Myr-old accreting intermediate-mass star showing strong infrared excess (L$_{\rm IR}$/L$_{\rm bol}$ ~ 0.25). Two extremely similar eclipses with a depth of ~30\% and duration ~25 days were observed in November 2008 and January 2011. We interpret the eclipses as caused by the same structure with an orbital period of $808\pm2$ days. Shearing over a single orbit rules out diffuse dust clumps as the cause, favouring the hypothesis of a companion at ~2AU. The characteristics of the eclipses are consistent with transits by an unseen low-mass (1.8-70M$_{Jup}$) planet or brown dwarf with a circum-secondary disc of diameter ~0.3 AU. The next eclipse event is predicted to take place in September 2017 and could be monitored by amateur and professional observatories across the world.
△ Less
Submitted 29 May, 2017;
originally announced May 2017.