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.

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

We investigate the properties of X-ray emission from accretion shocks in classical T Tauri stars (CTTSs), generated where the infalling material impacts the stellar surface. Both observations and models of the accretion process reveal several aspects that are still unclear: the observed X-ray luminosity in accretion shocks is below the predicted value, and the density versus temperature structure of the shocked plasma, with increasing densities at higher temperature, deduced from the observations, is at odds with that proposed in the current picture of accretion shocks. To address these open issues we investigate whether a correct treatment of the local absorption by the surrounding medium is crucial to explain the observations. To this end, we describe the impact of an accretion stream on a CTTS by considering a magnetohydrodynamic model. From the model results we synthesize the X-ray emission from the accretion shock by producing maps and spectra. We perform density and temperature diagnostics on the synthetic spectra, and we directly compare the results with the observations. Our model shows that the X-ray fluxes inferred from the emerging spectra are lower than expected because of the complex local absorption by the optically thick material of the chromosphere and of the unperturbed stream. Moreover, our model including the effects of local absorption explains in a natural way the apparently puzzling pattern of density versus temperature observed in the X-ray emission from accretion shocks.

We wish to test the origins proposed for the extended hot H2 at 2000K around the atomic jet from the T Tauri star DGTau, in order to constrain the wide-angle wind structure and the possible presence of an MHD disk wind. We present flux calibrated IFS observations in H2 1-0 S(1) obtained with SINFONI/VLT. Thanks to spatial deconvolution by the PSF and to accurate correction for uneven slit illumination, we performed a thorough analysis and modeled the morphology, kinematics, and surface brightness. We also compared our results with studies in [FeII], [OI], and FUV-pumped H2. The limb-brightened H2 emission in the blue lobe is strikingly similar to FUV-pumped H2 imaged 6yr later, confirming that they trace the same hot gas and setting an upper limit of 12km/s on any expansion proper motion. The wide-angle H2 rims are at lower blueshifts than probed by narrow long-slit spectra. We confirm that they extend to larger angle and to lower speed the onion-like velocity structure observed in optical atomic lines. The latter is shown to be steady over more/equal than 4yr but undetected in [FeII] by SINFONI, probably due to strong iron depletion. The H2 rim thickness less/equal than 14AU rules out excitation by C-shocks, and J-shock speeds are constrained to 10km/s. We find that explaining the H2 wide-angle emission with a shocked layer requires either a recent outburst (15yr) into a pre-existing ambient outflow or an excessive wind mass flux. A slow photoevaporative wind from the dense irradiated disk surface and an MHD disk wind heated by ambipolar diffusion seem to be more promising and need to be modeled in more detail.

(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.

The presence of a close, low-mass companion is thought to play a substantial and perhaps necessary role in shaping post-Asymptotic Giant Branch and Planetary Nebula outflows. During post-main-sequence evolution, radial expansion of the primary star, accompanied by intense winds, can significantly alter the binary orbit via tidal dissipation and mass loss. To investigate this, we couple stellar evolution models (from the zero-age main-sequence through the end of the post-main sequence) to a tidal evolution code. The binary's fate is determined by the initial masses of the primary and the companion, the initial orbit (taken to be circular), and the Reimers mass-loss parameter. For a range of these parameters, we determine whether the orbit expands due to mass loss or decays due to tidal torques. Where a common envelope (CE) phase ensues, we estimate the final orbital separation based on the energy required to unbind the envelope. These calculations predict period gaps for planetary and brown dwarf companions to white dwarfs. The upper end of the gap is the shortest period at which a CE phase is avoided. The lower end is the longest period at which companions survive their CE phase. For binary systems with 1 $M_\odot$ progenitors, we predict no Jupiter-mass companions with periods $\lesssim$270 days. Once engulfed, Jupiter-mass companions do not survive a CE phase. For binary systems consisting of a 1 $M_\odot$ progenitor with a companion 10 times the mass of Jupiter, we predict a period gap between $\sim$0.1 and $\sim$380 days. These results are consistent with both the detection of a $\sim$50 $M_{\rm J}$ brown dwarf in a $\sim$0.003 AU ($\sim$0.08 day) orbit around the white dwarf WD 0137-349 and the tentative detection of a $\sim$2 $M_{\rm J}$ planet in a $\gtrsim$2.7 AU ($\gtrsim$4 year) orbit around the white dwarf GD66.

In order to explain the inflated radii of some transiting extrasolar giant planets, we investigate a tidal heating scenario for the inflated planets WASP-4b, WASP-6b, WASP-12b, WASP-15b, and TrES-4. To do so, we assume that they retain a nonzero eccentricity, possibly by dint of continuing interaction with a third body. We calculate the amount of extra heating in the envelope that is then required to fit the radius of each planet, and we explore how this additional power depends on the planetary atmospheric opacity and on the mass of a heavy-element central core. There is a degeneracy between the core mass $M_{\rm core}$ and the heating $\dot{E}_{\rm heating}$. Therefore, in the case of tidal heating, there is for each planet a range of the couple $\{M_{\rm core},e^2/Q'_p\}$ that can lead to the same radius, where $Q'_p$ is the tidal dissipation factor and $e$ is the eccentricity. With this in mind, we also investigate the case of the non-inflated planet HAT-P-12b, which can admit solutions combining a heavy-element core and tidal heating. A substantial improvement of the measured eccentricities of such planetary systems could simplify this degeneracy by linking the two unknown parameters $\{M_{\rm core},Q'_p\}$. Further independent constraints on either of these parameters would, through our calculations, constrain the other.

There are now many known exoplanets with Msin(i) within a factor of two of Neptune's, including the transiting planets GJ436b and HAT-P-11b. Planets in this mass-range are different from their more massive cousins in several ways that are relevant to their radiative properties and thermal structures. By analogy with Neptune and Uranus, they are likely to have metal abundances that are an order of magnitude or more greater than those of larger, more massive planets. This increases their opacity, decreases Rayleigh scattering, and changes their equation of state. Furthermore, their smaller radii mean that fluxes from these planets are roughly an order of magnitude lower than those of otherwise identical gas giant planets. Here, we compute a range of plausible radiative equilibrium models of GJ436b and HAT-P-11b. In addition, we explore the dependence of generic Neptune-mass planets on a range of physical properties, including their distance from their host stars, their metallicity, the spectral type of their stars, the redistribution of heat in their atmospheres, and the possible presence of additional optical opacity in their upper atmospheres.

The radii of some transiting extrasolar giant planets are larger than would be expected by the standard theory. We address this puzzle with the model of coupled radius-orbit tidal evolution developed by \citet{Ibgui_and_Burrows_2009}. The planetary radius is evolved self-consistently with orbital parameters, under the influence of tidal torques and tidal dissipation in the interior of the planet. A general feature of this model, which we have previously demonstrated in the generic case, is that a possible transient inflation of the planetary radius can temporarily interrupt its standard monotonic shrinking and can lead to the inflated radii that we observe. In particular, a bloated planet with even a circular orbit may still be inflated due to an earlier episode of tidal heating. We have modified our model to include an orbital period dependence of the tidal dissipation factor in the star, $Q'_{\ast} \propto P^{\gamma}$, $-1 \leqslant \gamma \leqslant 1$. With this model, we search, for a tidally heated planet, orbital and radius evolutionary tracks that fall within the observational limits of the radius, the semimajor axis, and the eccentricity of the planet in its current estimated age range. We find that, for some inflated planets (WASP-6b and WASP-15b), there are such tracks; for another (TrES-4), there are none; and for still others (WASP-4b and WASP-12b), there are such tracks, but our model might imply that we are observing the planets at a special time. Finally, we stress that there is a two to three order-of-magnitude timescale uncertainty of the inspiraling phase of the planet into its host star, arising from uncertainties in the tidal dissipation factor in the star $Q'_{\ast}$.

Some transiting extrasolar giant planets have measured radii larger than predicted by the standard theory. In this paper, we explore the possibility that an earlier episode of tidal heating can explain such radius anomalies and apply the formalism we develop to HD 209458b as an example. We find that for strong enough tides the planet's radius can undergo a transient phase of inflation that temporarily interrupts canonical, monotonic shrinking due to radiative losses. Importantly, an earlier episode of tidal heating can result in a planet with an inflated radius, even though its orbit has nearly circularized. Moreover, we confirm that at late times, and under some circumstances, by raising tides on the star itself a planet can spiral into its host. We note that a 3$\times$ to 10$\times$solar planet atmospheric opacity with no tidal heating is sufficient to explain the observed radius of HD 209458b. However, our model demonstrates that with an earlier phase of episodic tidal heating we can fit the observed radius of HD 209458b even with lower (solar) atmospheric opacities. This work demonstrates that, if a planet is left with an appreciable eccentricity after early inward migration and/or dynamical interaction, coupling radius and orbit evolution in a consistent fashion that includes tidal heating, stellar irradiation, and detailed model atmospheres might offer a generic solution to the inflated radius puzzle for transiting extrasolar giant planets such as WASP-12b, TrES-4, and WASP-6b.

To explain their observed radii, we present theoretical radius-age trajectories for the extrasolar giant planets (EGPs) TrES-4, XO-3b, and HAT-P-1b. We factor in variations in atmospheric opacity, the presence of an inner heavy-element core, and possible heating due to orbital tidal dissipation. A small, yet non-zero, degree of core heating is needed to explain the observed radius of TrES-4, unless its atmospheric opacity is significantly larger than a value equivalent to that at 10$\times$solar metallicity with equilibrium molecular abundances. This heating rate is reasonable, and corresponds for an energy dissipation parameter ($Q_p$) of $\sim10^{3.8}$ to an eccentricity of $\sim$0.01, assuming 3$\times$solar atmospheric opacity and a heavy-element core of $M_c = 30$ $M_{\oplus}$. For XO-3b, which has an observed orbital eccentricity of 0.26, we show that tidal heating needs to be taken into account to explain its observed radius. Furthermore, we reexamine the core mass needed for HAT-P-1b in light of new measurements and find that it now generally follows the correlation between stellar metallicity and core mass suggested recently. Given various core heating rates, theoretical grids and fitting formulae for a giant planet's equilibrium radius and equilibration timescale are provided for planet masses $M_p=$ 0.5, 1.0, and 1.5 $M_J$ with $a =$ 0.02-0.06 AU, orbiting a G2V star. When the equilibration timescale is much shorter than that of tidal heating variation, the ``effective age'' of the planet is shortened, resulting in evolutionary trajectories more like those of younger EGPs. Motivated by the work of \citet{jackson08a,jackson08b}, we suggest that this effect could indeed be important in better explaining some observed transit radii.

We calculate a new suite of albedo models for close-in extrasolar giant planets and compare with the recent stringent upper limit for HD 209458b of Rowe et al. using MOST. We find that all models without scattering clouds are consistent with this optical limit. We explore the dependence on wavelength and waveband, metallicity, the degree of heat redistribution, and the possible presence of thermal inversions and find a rich diversity of behaviors. Measurements of transiting extrasolar giant planets (EGPs) at short wavelengths by MOST, Kepler, and CoRoT, as well as by proposed dedicated multi-band missions, can complement measurements in the near- and mid-IR using {\it Spitzer} and JWST. Collectively, such measurements can help determine metallicity, compositions, atmospheric temperatures, and the cause of thermal inversions (when they arise) for EGPs with a broad range of radii, masses, degrees of stellar insolation, and ages. With this paper, we reappraise and highlight the diagnostic potential of albedo measurements of hot EGPs shortward of $\sim$1.3 $\mu$m.