Models and observations indicate that the impact of matter accreting onto the surface of young stars produces regions at the base of accretion columns, in which optically thin and thick plasma components coexist. Thus an accurate description of these impacts requires to account for the effects of absorption and emission of radiation. We study the effects of radiation emerging from shock-heated plasma in impact regions on the structure of the pre-shock downfalling material. We investigate if a significant absorption of radiation occurs and if it leads to a pre-shock heating of the accreting gas. We developed a radiation hydrodynamics model describing an accretion column impacting onto the surface of a Classical T Tauri Star. The model takes into account the stellar gravity, the thermal conduction, and the effects of both radiative losses and absorption of radiation by matter in the non local thermodynamic equilibrium regime. After the impact, a hot slab of post-shock plasma develops at the base of the accretion column. Part of radiation emerging from the slab is absorbed by the pre-shock accreting material. As a result, the pre-shock accretion column gradually heats up to temperatures of $10^5$ K, forming a radiative precursor of the shock. The precursor has a thermal structure with the hottest part at $T \approx 10^5$ K, with size comparable to that of the hot slab, above the post-shock region. At larger distances the temperature gradually decreases to $T \approx 10^4$ K.

Context. Accreted material onto CTTSs is expected to form a hot quasi-periodic plasma structure that radiates in X-rays. Simulations of this phenomenon only partly match with observations. They all rely on a static model for the chromosphere model and on the assumption that radiation and matter are decoupled. Aims. We explore the effects on the structure and on the dynamics of the accretion flow of both a shock-heated chromosphere and of the coupling between radiation and hydrodynamics. Methods. We simulate accretion columns falling onto a stellar chromosphere using the 1D ALE code AstroLabE. This code solves the hydrodynamics equations along with the two first momenta equations for radiation transfer, with the help of a dedicated opacity table for the coupling between matter and radiation. We derive the total electron and ions densities from collisional-radiative NLTE ionisation equilibrium. Results. The chromospheric acoustic heating has an impact on the duration of the cycle and on the structure of the heated slab. In addition, the coupling between radiation and hydrodynamics leads to a heating of the accretion flow and the chromosphere, inducing a possible unburial of the whole column. These two last conclusions are in agreement with the computed monochromatic intensity. Both effects (acoustic heating and radiation coupling) have an influence on the amplitude and temporal variations of the net X-ray luminosity, which varies between 30 and 94% of the incoming mechanical energy flux, depending on the model considered.

Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those among them that use either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to the local thermodynamic equilibrium (LTE). This assumption may be not valid in some astrophysical cases. We present an upgraded version of the LTE radiation-hydrodynamics module implemented in the PLUTO code, originally developed by Kolb et al. (2013), which we have extended to handle non-LTE regimes. Starting from the general frequency-integrated comoving-frame equations of radiation hydrodynamics (RHD), we have justified all the assumptions made to obtain the non-LTE equations actually implemented in the module, under the FLD approximation. An operator-split method is employed, with two substeps: the hydrodynamic part is solved with an explicit method by the solvers already available in PLUTO, the non-LTE radiation diffusion and energy exchange part is solved with an implicit method. The module is implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: the radiative power loss, the Planck and Rosseland mean opacities. Our implementation has been validated through different tests, in particular radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that non-LTE approach is of paramount importance to properly model accretion shock structures. Our radiation FLD module represents a step toward the general non-LTE RHD modeling. The module is available, under request, for the community.

Classical T Tauri stars (CTTSs) are young low-mass stellar objects accreting mass from their circumstellar disks. They are characterized by high levels of coronal activity as revealed by X-ray observations. This activity may affect the disk stability and the circumstellar environment. Here we investigate if an intense coronal activity due to flares occurring close to the accretion disk may perturb the inner disk stability, disrupt the inner part of the disk and, possibly, trigger accretion phenomena with rates comparable with those observed. We model a magnetized protostar surrounded by an accretion disk through 3D magnetohydrodinamic simulations. We explore cases characterized by a dipole plus an octupole stellar magnetic field configuration and different density of the disk or by different levels of flaring activity. As a result of the simulated intense flaring activity, we observe the formation of several loops that link the star to the disk; all these loops build up a hot extended corona with an X-ray luminosity comparable with typical values observed in CTTSs. The intense flaring activity close to the disk can strongly perturb the disk stability. The flares trigger overpressure waves which travel through the disk and modify its configuration. Accretion funnels may be triggered by the flaring activity, thus contributing to the mass accretion rate of the star. Accretion rates synthesized from the simulations are in a range between 10^-10 and 10^-9M_sun yr^-1 The accretion columns can be perturbed by the flares and they can interact with each other, possibly merging together in larger streams. As a result, the accretion pattern can be rather complex: the streams are highly inhomogeneous, with a complex density structure, and clumped

Past high-energy density laboratory experiments provided insights into the physics of supernovae, supernova remnants, and the destruction of interstellar clouds. In a typical experimental setting, a laser-driven planar blast wave interacts with a compositionally-homogeneous spherical or cylindrical target. In this work we propose a new laboratory platform that accounts for curvature of the impacting shock and density stratification of the target. Both characteristics reflect the conditions expected to exist shortly after a supernova explosion in a close binary system. We provide details of a proposed experimental design (laser drive, target configuration, diagnostic system), optimized to capture the key properties of recent ejecta-companion interaction models. Good qualitative agreement found between our experimental models and their astrophysical counterparts highlights strong potential of the proposed design to probe details of the ejecta-companion interaction for broad classes of objects by means of high energy density laboratory experiments.

(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band.

Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the shock structure influencing the observability of the oscillatory behavior. Methods. We performed local 2D magneto-hydrodynamical simulations of an accretion shock impacting a chromosphere, taking optically thin radiation losses and thermal conduction into account. We investigated the effects of several perturbation types, such as clumps in the accretion stream or chromospheric fluctuations, and also explored a wide range of plasma-\beta values. Results. In the case of a weak magnetic field, the post-shock region shows chaotic motion and mixing, smoothing out the perturbations and retaining a global periodic signature. On the other hand, a strong magnetic field confines the plasma in flux tubes, which leads to the formation of fibrils that oscillate independently. Realistic values for the amplitude, length, and time scales of the perturbation are capable of bringing the fibril oscillations out of phase, suppressing the periodicity of the emission. Conclusions. The strength of a locally uniform magnetic field in YSO accretion shocks determines the structure of the post-shock region, namely, whether it will be somewhat homogeneous or if it will split up to form a collection of fibrils. In the second case, the size and shape of the fibrils is found to depend strongly on the plasma-\beta value but not on the perturbation type. Therefore, the actual value of the protostellar magnetic field is expected to play a critical role in the time dependence of the observable emission.

We present IRIS, a new generic three-dimensional (3D) spectral radiative transfer code that generates synthetic spectra, or images. It can be used as a diagnostic tool for comparison with astrophysical observations or laboratory astrophysics experiments. We have developed a 3D short-characteristic solver that works with a 3D nonuniform Cartesian grid. We have implemented a piecewise cubic, locally monotonic, interpolation technique that dramatically reduces the numerical diffusion effect. The code takes into account the velocity gradient effect resulting in gradual Doppler shifts of photon frequencies and subsequent alterations of spectral line profiles. It can also handle periodic boundary conditions. This first version of the code assumes Local Thermodynamic Equilibrium (LTE) and no scattering. The opacities and source functions are specified by the user. In the near future, the capabilities of IRIS will be extended to allow for non-LTE and scattering modeling. IRIS has been validated through a number of tests. We provide the results for the most relevant ones, in particular a searchlight beam test, a comparison with a 1D plane-parallel model, and a test of the velocity gradient effect. IRIS is a generic code to address a wide variety of astrophysical issues applied to different objects or structures, such as accretion shocks, jets in young stellar objects, stellar atmospheres, exoplanet atmospheres, accretion disks, rotating stellar winds, cosmological structures. It can also be applied to model laboratory astrophysics experiments, such as radiative shocks produced with high power lasers.

We report on the investigation of strong radiative shocks generated with the high energy, sub-nanosecond iodine laser at PALS. These shock waves are characterized by a developed radiative precursor and their dynamics is analyzed over long time scales (~50 ns), approaching a quasi-stationary limit. We present the first preliminary results on the rear side XUV spectroscopy. These studies are relevant to the understanding of the spectroscopic signatures of accretion shocks in Classical T Tauri Stars.

Stark broadened lines of the hydrogen Brackett series are computed for the conditions of stellar atmospheres and circumstellar envelopes. The computation is performed within the Model Microfield Method, which includes the ion dynamic effects and makes the bridge between the impact limit at low density and the static limit at high density and in the line wings. The computation gives the area normalized line shape, from the line core up to the static line wings.

Radiative shocks are found in various astrophysical objects and particularly at different stages of stellar evolution. Studying radiative shocks, their topology, and thermodynamical properties is therefore a starting point to understanding their physical properties. This study has become possible with the development of large laser facilities, which has provided fresh impulse to laboratory astrophysics. We present the main characteristics of radiative shocks modeled using cylindrical simulations. We focus our discussion on the importance of multi-dimensional radiative-transfer effects on the shock topology and dynamics. We present results obtained with our code HERACLES for conditions corresponding to experiments already performed on laser installations. The multi-dimensional hydrodynamic code HERACLES is specially adapted to laboratory astrophysics experiments and to astrophysical situations where radiation and hydrodynamics are coupled. The importance of the ratio of the photon mean free path to the transverse extension of the shock is emphasized. We present how it is possible to achieve the stationary limit of these shocks in the laboratory and analyze the angular distribution of the radiative flux that may emerge from the walls of the shock tube. Implications of these studies for stellar accretion shocks are presented.

Radiative shocks (also called supercritical shocks) are high Mach number shock waves that photoionize the medium ahead of the shock front and give rise to a radiative precursor. They are generated in the laboratory using high-energy or high-power lasers and are frequently present in a wide range of astronomical objects. Their modelisation in one dimension has been the subject of numerous studies, but generalization to three dimensions is not straightforward. We calculate analyticaly the absorption of radiation in a grey uniform cylinder and show how it decreases with chi R, the product of the opacity chi and of the cylinder radius R. Simple formulas, whose validity range increases when chi R diminishes, are derived for the radiation field on the axis of symmetry. Numerical calculations in three dimensions of the radiative energy density, flux and pressure created by a stationary shock wave show how the radiation decreases whith R. Finally, the bidimensional structures of both the precursor and the radiation field are calculated with time-dependent radiation hydrodynamics numerical simulations and the influence of two-dimensional effects on the electron density, the temperature, the shock velocity and the shock geometry are exhibited. These simulations show how the radiative precursor shortens, cools and slows down when R is decreased.