skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Unstable accretion in TW Hya: 3D simulations and comparisons with observations
ABSTRACT We investigate the origin of photometric variability in the classical T Tauri star TW Hya by comparing light curves obtained by Transiting Exoplanet Survey Satellite (TESS) and ground-based telescopes with light curves created using three-dimensional (3D) magnetohydrodynamic (MHD) simulations. TW Hya is modelled as a rotating star with a dipole magnetic moment, which is slightly tilted about the rotational axis. We observed that for various model parameters, matter accretes in the unstable regime and produces multiple hotspots on the star’s surface, which leads to stochastic-looking light curves similar to the observed ones. Wavelet and Fourier spectra of observed and modelled light curves show multiple quasi-periodic oscillations (QPOs) with quasi-periods from less than 0.1 to 9 d. Models show that variation in the strength and tilt of the dipole magnetosphere leads to different periodograms, where the period of the star may dominate or be hidden. The amplitude of QPOs associated with the stellar period can be smaller than that of other QPOs if the tilt of the dipole magnetosphere is small and when the unstable regime is stronger. In models with small magnetospheres, the short-period QPOs associated with rotation of the inner disc dominate and can be mistaken for a stellar period. We show that longer period (5–9 d) QPOs can be caused by waves forming beyond the corotation radius.  more » « less
Award ID(s):
2009820
PAR ID:
10575432
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
538
Issue:
1
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 480-502
Size(s):
p. 480-502
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Classical T Tauri stars (CTTSs) are young, low-mass stars that accrete material from their surrounding protoplanetary disk. To better understand accretion variability, we conducted a multiepoch, multiwavelength photometric monitoring campaign of four CTTSs, TW Hya, RU Lup, BP Tau, and GM Aur, in 2021 and 2022, contemporaneous with Hubble Space Telescope 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. There 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. 
    more » « less
  2. ABSTRACT In this paper we present the first set of 3D magnetohydrodynamic (MHD) simulations performed with the riemann geomesh code. We study the dynamics of the magnetically channeled winds of magnetic massive stars in full three dimensions using a code that is uniquely suited to spherical problems. Specifically, we perform isothermal simulations of a smooth wind on a rotating star with a tilted, initially dipolar field. We compare the mass-loss, angular momentum loss, and magnetospheric dynamics of a template star (with the properties that are reminiscent of the O4 supergiant ζ Pup) over a range of rotation rates, magnetic field strengths, and magnetic tilt angles. The simulations are run up to a quasi-steady state and the results are observed to be consistent with the existing literature, showing the episodic centrifugal breakout events of the mass outflow, confined by the magnetic field loops that form the closed magnetosphere of the star. The catalogued results provide perspective on how angular-momentum loss varies for different configurations of rotation rate, magnetic field strength, and large magnetic tilt angles. In agreement with previous 2D MHD studies, we find that high magnetic confinement reduces the overall mass-loss rate, and higher rotation increases the mass-loss rate. This and future studies will be used to estimate the angular-momentum evolution, spin-down time, and mass-loss evolution of magnetic massive stars as a function of magnetic field strength, rotation rate, and dipole tilt. 
    more » « less
  3. We present optical and infrared (IR) light curves of the enshrouded massive binary NaSt1 (WR 122) with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), and the All-Sky Automated Survey for Supernovae (ASAS-SN). The optical and IR light curves span between 2014 July and 2020 Oct., revealing periodic, sinusoidal variability from NaSt1 with a P=305.2±1.0 d period. We also present historical IR light curves taken between 1983 July and 1989 May that also indicate NaSt1 exhibits long-term IR variability on timescales of ∼decades. Fixed-period sinusoidal fits to the recent optical and IR light curves show that amplitude of NaSt1's variability is different at different wavelengths and also reveal significant phase offsets of ∼18 d between the ZTF r and PGIR J light curves.We interpret the ∼300 d period of the observed variability as the orbital period of a binary system in NaSt1. Assuming a circular orbit and adopting a range of combined stellar mass values in the range 20-100 M⊙ in NaSt1, we estimate orbital separations of ∼2-4 au. We suggest that the sinusoidal photometric variability of NaSt1 may arise from variations in the line-of-sight optical depth toward circumstellar optical/IR emitting regions throughout its orbit due to colliding-wind dust formation. We provide an interpretation on the nature of NaSt1 and speculate that the mass-transfer process may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of a Ofpe/WN9 star. Lastly, we claim that NaSt1 ceased RLOF mass transfer ≲3400 yr ago. 
    more » « less
  4. ABSTRACT Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate non-thermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anticorrelated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated with the kink growth time-scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities. 
    more » « less
  5. null (Ed.)
    ABSTRACT We report the discovery of the closest known black hole candidate as a binary companion to V723 Mon. V723 Mon is a nearby ($$d\sim 460\, \rm pc$$), bright (V ≃ 8.3 mag), evolved (Teff, giant ≃ 4440 K, and Lgiant ≃ 173 L⊙) red giant in a high mass function, f(M) = 1.72 ± 0.01 M⊙, nearly circular binary (P = 59.9 d, e ≃ 0). V723 Mon is a known variable star, previously classified as an eclipsing binary, but its All-Sky Automated Survey, Kilodegree Extremely Little Telescope, and Transiting Exoplanet Survey Satellite light curves are those of a nearly edge-on ellipsoidal variable. Detailed models of the light curves constrained by the period, radial velocities, and stellar temperature give an inclination of $$87.0^{\circ ^{+1.7^\circ }}_{-1.4^\circ }$$, a mass ratio of q ≃ 0.33 ± 0.02, a companion mass of Mcomp = 3.04 ± 0.06 M⊙, a stellar radius of Rgiant = 24.9 ± 0.7 R⊙, and a giant mass of Mgiant = 1.00 ± 0.07 M⊙. We identify a likely non-stellar, diffuse veiling component with contributions in the B and V band of $${\sim }63{{\ \rm per\ cent}}$$ and $${\sim }24{{\ \rm per\ cent}}$$, respectively. The SED and the absence of continuum eclipses imply that the companion mass must be dominated by a compact object. We do observe eclipses of the Balmer lines when the dark companion passes behind the giant, but their velocity spreads are low compared to observed accretion discs. The X-ray luminosity of the system is $$L_{\rm X}\simeq 7.6\times 10^{29}~\rm ergs~s^{-1}$$, corresponding to L/Ledd ∼ 10−9. The simplest explanation for the massive companion is a single compact object, most likely a black hole in the ‘mass gap’. 
    more » « less