X-ray emission from dense plasma in CTTSs: Hydrodynamic modeling of the accretion shock

TL;DR

This study uses hydrodynamic simulations to analyze the stability and properties of accretion shock-heated plasma in classical T Tauri stars, explaining observed X-ray emissions and their variability.

Contribution

It provides the first detailed 1-D hydrodynamic model including chromospheric effects to study the dynamics and stability of accretion shocks in CTTSs.

Findings

Shocked region forms a hot slab with specific density and temperature.

The slab exhibits quasi-periodic instabilities and cyclic formation.

Simulated X-ray luminosity and spectrum match observations of MP Mus.

Abstract

High spectral resolution X-ray observations of CTTSs demonstrate the presence of plasma at T~2-3X10^6 K and n_e~10^11-10^13 cm^-3, unobserved in non-accreting stars. Stationary models suggest that this emission is due to shock-heated accreting material, but they do not allow to analyze the stability of such material and its position in the stellar atmosphere. We investigate the dynamics and the stability of shock-heated accreting material in CTTSs and the role of the stellar chromosphere in determining the position and the thickness of the shocked region. We perform 1-D HD simulations of the impact of the accretion flow onto chromosphere of a CTTS, including the effects of gravity, radiative losses from optically thin plasma, thermal conduction and a well tested detailed model of the stellar chromosphere. Here we present the results of a simulation based on the parameters of the CTTS MP Mus. We find that the accretion shock generates an hot slab of material above the chromosphere with a maximum thickness of 1.8X10^9 cm, density n_e~10^11-10^2 cm^-3, temperature T~3X10^6 K and uniform pressure equal to the ram pressure of the accretion flow (~450 dyn cm^-2). The base of the shocked region penetrates the chromosphere and stays where the ram pressure is equal to the thermal pressure. The system evolves with quasi-periodic instabilities of the material in the slab leading to cyclic disappearance and re-formation of the slab. For an accretion rate of ~10^-10 M_sun yr^-1, the shocked region emits a time-averaged X-ray luminosity L_X~7X10^29 erg s^-1, which is comparable to the X-ray luminosity observed in CTTSs of the same mass. Furthermore, the X-ray spectrum synthesized from the simulation matches in detail all the main features of the O VIII and O VII lines of the star MP Mus.