Investigating in the laboratory the process of matter accretion onto forming stars through scaled experiments is important in order to better understand star and planetary systems formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A first step has been made in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528]. in allowing such investigations. It revealed the existence of a two components stream: a hot shell surrounding a cooler inner stream. The shell was formed by matter laterally ejected upon impact and refocused by the local magnetic field. That laboratory investigation was limited tonormal incidence impacts. However, in young stellar objects, complex structure of magnetic fields causes variability of the incidentangles of the accretion columns. This led us to undertake an investigation, using laboratory plasmas, of the consequence of having a slanted accretion impacting a young star. Here we use high power laser interactions and strong magnetic field generation in the laboratory, complemented by numerical simulations, to study the asymmetry induced upon accretion structures when columns of matter impact the surface of young stars with an oblique angle. Compared to the scenario where matter accretes normal to the star surface, we observe strongly asymmetric plasma structure, strong lateral ejecta of matter, poor confinement of the accreted material and reduced heating compared to the normal incidence case. Thus, slanted accretion is a configuration that seems to be capable of inducing perturbations of the chromosphere and hence possibly influence the level of activity of the corona.

The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growing misalignments, this nozzle becomes increasingly asymmetric, disrupting jet formation. Our results thus suggest outflow/magnetic field misalignment to be a plausible key process regulating jet collimation in a variety of objects. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field.

A new experimental set-up designed to investigate the accretion dynamics in newly born stars is presented. It takes advantage of a magnetically collimated stream produced by coupling a laser-generated expanding plasma to a $2\times 10^{5}~{G}\ (20~{T})$ externally applied magnetic field. The stream is used as the accretion column and is launched onto an obstacle target that mimics the stellar surface. This setup has been used to investigate in details the accretion dynamics, as reported in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528}. Here, the characteristics of the stream are detailed and a link between the experimental plasma expansion and a 1D adiabatic expansion model is presented. Dimensionless numbers are also calculated in order to characterize the experimental flow and its closeness to the ideal MHD regime. We build a bridge between our experimental plasma dynamics and the one taking place in the Classical T Tauri Stars (CTTSs), and we find that our set-up is representative of a high plasma $\beta$ CTTS accretion case.

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.

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars.

Coronal mass ejections (CMEs), often associated with flares, are the most powerful magnetic phenomena occurring on the Sun. Stars show magnetic activity levels up to 10^4 times higher, and CME effects on stellar physics and circumstellar environments are predicted to be significant. However, stellar CMEs remain observationally unexplored. Using time-resolved high-resolution X-ray spectroscopy of a stellar flare on the active star HR 9024 observed with Chandra/HETGS, we distinctly detected Doppler shifts in S XVI, Si XIV, and Mg XII lines that indicate upward and downward motions of hot plasmas (~10-25 MK) within the flaring loop, with velocity v~100-400 km/s, in agreement with a model of flaring magnetic tube. Most notably, we also detected a later blueshift in the O VIII line which reveals an upward motion, with v=90+/-30 km/s, of cool plasma (~4 MK), that we ascribe to a CME coupled to the flare. From this evidence we were able to derive a CME mass of 1x10^21 g and a CME kinetic energy of 5x10^34 erg. These values provide clues in the extrapolation of the solar case to higher activity levels, suggesting that CMEs could indeed be a major cause of mass and angular momentum loss.

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories.

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

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

One of the aims of LSST is to perform a systematic survey of star clusters and star forming regions (SFRs) in our Galaxy. In particular, the observations obtained with LSST will make a big difference in Galactic regions that have been poorly studied in the past, such as the anticenter and the disk beyond the Galactic center, and they will have a strong impact in discovering new distant SFRs. These results can be achieved by exploiting the exquisite depth that will be attained if the wide-fast-deep (WFD) observing strategy of the main survey is also adopted for the Galactic plane, in the g, r, and i filters.

Several works have found an increase of the abundances of the s-process neutron-capture elements in the youngest Galactic stellar populations, giving important constraints to stellar and Galactic evolution. We aim to trace the abundance patterns and the time-evolution of five s-process elements in the first peak, Y and Zr, and in the second peak, Ba, La and Ce using the Gaia-ESO idr5 results. From the UVES spectra of cluster member stars, we determined the average composition of clusters with ages >0.1 Gyr. We derived statistical ages and distances of field stars, and we separated them in thin and thick disc populations. We studied the time evolution and dependence on metallicity of abundance ratios using open clusters and field stars. Using our large and homogeneous sample of open clusters, thin and thick disc stars, spanning an age range larger than 10 Gyr, we confirm an increase towards young ages of s-process abundances in the Solar neighbourhood. These trends are well defined for open clusters and stars located nearby the solar position and they may be explained by a late enrichment due to significant contribution to the production of these elements from long-living low-mass stars. At the same time, we found a strong dependence of the s-process abundance ratios with the Galactocentric distance and with the metallicity of the clusters and field stars. Our results, derived from the largest and homogeneous sample of s-process abundances in the literature, confirm the growth with decreasing stellar ages of the s-process abundances in both field and open cluster stars. At the same time, taking advantage of the abundances of open clusters located in a wide Galactocentric range, they open a new view on the dependence of the s-process evolution on the metallicity and star formation history, pointing to different behaviours at various Galactocentric distances.

Several observations of stellar jets show evidence of X-ray emitting shocks close to the launching site. In some cases, the shocked features appear to be stationary, also for YSOs at different stages of evolution. We study the case of HH 154, the jet originating from the embedded binary Class 0/I protostar IRS 5, and the case of the jet associated to DG Tau, a more evolved Class II disk-bearing source or Classical T Tauri star (CTTS), both located in the Taurus star-forming region. We aim at investigating the effect of perturbations in X-ray emitting stationary shocks in stellar jets; the stability and detectability in X-rays of these shocks; and explore the differences in jets from Class 0 to Class II sources. We performed a set of 2.5-dimensional magnetohydrodynamic numerical simulations that modelled supersonic jets ramming into a magnetized medium. The jet is formed by two components: a continously driven component that forms a quasi-stationary shock at the base of the jet; and a pulsed component constituted by blobs perturbing the shock. We explored different parameters for both components. We studied two cases: a jet less dense than the ambient medium (light jet), representing the case of HH 154; and a jet denser than the ambient (heavy jet), associated with DG Tau. We synthesized the count rate from the simulations and compared with available Chandra observations. Our model explains the formation of X-ray emitting quasi-stationary shocks observed at the base of jets in a natural way, being able to reproduce the observed jet properties at different evolutionary phases (in particular, for HH 154 and DG Tau). The jet is collimated by the magnetic field forming a quasi-stationary shock at the base which emits in X-rays even when perturbations formed by a train of blobs are present. We found similar collimation mechanisms dominating in both heavy and light jets...

We explore the structure and star formation history of the open cluster NGC 2264 (~3 Myr). We combined spectroscopic data from the Gaia-ESO Survey (GES) with multi-wavelength photometry from the Coordinated Synoptic Investigation of NGC 2264 (CSI 2264). We examined a sample of 655 cluster members, including both disk-bearing and disk-free young stars. We find a significant age spread of 4-5 Myr among cluster members. Disk-bearing objects are statistically associated with younger isochronal ages than disk-free sources. The cluster has a hierarchical structure, with two main blocks. The northern half develops around the O-type binary star S Mon; the southern half, close to the tip of the Cone Nebula, contains the most embedded regions of NGC 2264, populated mainly by objects with disks and ongoing accretion. The median ages of objects at different locations within the cluster, and the spatial distribution of disked and non-disked sources, suggest that star formation began in the north of the cluster, over 5 Myr ago, and was ignited in its southern region a few Myr later. Star formation is likely still ongoing in the most embedded regions of the cluster, while the outer regions host a widespread population of more evolved objects. We find a detectable lag between the typical age of disk-bearing objects and that of accreting objects in the inner regions of NGC 2264: the first tend to be older than the second, but younger than disk-free sources at similar locations within the cluster. This supports earlier findings that the characteristic timescales of disk accretion are shorter than those of disk dispersal, and smaller than the average age of NGC 2264 (i.e., < 3 Myr). At the same time, disks in the north of the cluster tend to be shorter-lived (~2.5 Myr) than elsewhere; this may reflect the impact of massive stars within the region (notably S Mon), that trigger rapid disk dispersal.

High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the motion of this X-ray-emitting post-shock material, we searched for a Doppler shift in the deep Chandra/HETGS observation of the CTTS TW Hya. This test should unveil the nature of this X-ray emitting plasma component in CTTS, and constrain the accretion stream geometry. We searched for a Doppler shift in the X-ray emission from TW Hya with two different methods, by measuring the position of a selected sample of emission lines, and by fitting the whole TW Hya X-ray spectrum, allowing the line-of-sight velocity to vary. We found that the plasma at T~2-4 MK has a line-of-sight velocity of 38.3+/-5.1 km/s with respect to the stellar photosphere. This result definitively confirms that this X-ray-emitting material originates in the post-shock region, at the base of the accretion stream, and not in coronal structures. The comparison of the observed velocity along the line of sight, 38.3+/-5.1 km/s, with the inferred intrinsic velocity of the post shock of TW Hya, v_post~110-120 km/s, indicates that the footpoints of the accretion streams on TW Hya are located at low latitudes on the stellar surface. Our results indicate that complex magnetic field geometries, such as that of TW Hya, permit low-latitude accretion spots. Moreover, since on TW Hya the redshift of the soft X-ray emission is very similar to that of the narrow component of the CIV resonance doublet at 1550 Ang, as found by Ardila et al. (2013), then the plasma at 2-4 MK and that at 0.1 MK likely originate in the same post-shock regions.

Accretion dynamics in the forming of young stars is still object of debate because of limitations in observations and modelling. Through scaled laboratory experiments of collimated plasma accretion onto a solid in the presence of a magnetic field, we open first window on this phenomenon by tracking, with spatial and temporal resolution, the dynamics of the system and simultaneously measuring multiband emissions. We observe in these experiments that matter, upon impact, is laterally ejected from the solid surface, then refocused by the magnetic field toward the incoming stream. Such ejected matter forms a plasma shell that envelops the shocked core, reducing escaped X-ray emission. This demonstrates one possible structure reconciling current discrepancies between mass accretion rates derived from X-ray and optical observations.

We present a spectroscopic study of the dynamics of the ionized and neutral gas throughout the Lagoon nebula (M8), using VLT/FLAMES data from the Gaia-ESO Survey. We explore the connections between the nebular gas and the stellar population of the associated star cluster NGC6530. We characterize through spectral fitting emission lines of H-alpha, [N II] and [S II] doublets, [O III], and absorption lines of sodium D doublet, using data from the FLAMES/Giraffe and UVES spectrographs, on more than 1000 sightlines towards the entire face of the Lagoon nebula. Gas temperatures are derived from line-width comparisons, densities from the [S II] doublet ratio, and ionization parameter from H-alpha/[N II] ratio. Although doubly-peaked emission profiles are rarely found, line asymmetries often imply multiple velocity components along the line of sight. This is especially true for the sodium absorption, and for the [O III] lines. Spatial maps for density and ionization are derived, and compared to other known properties of the nebula and of its massive stars 9 Sgr, Herschel 36 and HD 165052 which are confirmed to provide most of the ionizing flux. The detailed velocity fields across the nebula show several expanding shells, related to the cluster NGC6530, the O stars 9 Sgr and Herschel 36, and the massive protostar M8East-IR. The origins of kinematical expansion and ionization of the NGC6530 shell appear to be different. We are able to put constrains on the line-of-sight (relative or absolute) distances between some of these objects and the molecular cloud. The large obscuring band running through the middle of the nebula is being compressed by both sides, which might explain its enhanced density. We also find an unexplained large-scale velocity gradient across the entire nebula. At larger distances, the transition from ionized to neutral gas is studied using the sodium lines.

We present the first extensive spectroscopic study of the global population in star clusters Trumpler~16, Trumpler~14 and Collinder~232 in the Carina Nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. In addition to the standard homogeneous Survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. We find about four hundred good candidate members ranging from OB types down to slightly sub-solar masses. About one-hundred heavily-reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically-derived temperatures for nearly 300 low-mass members allows the inference of individual extinction values, and the study of the relative placement of stars along the line of sight. We find a complex spatial structure, with definite clustering of low-mass members around the most massive stars, and spatially-variable extinction. By combining the new data with existing X-ray data we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters, and of their connection to bright and dark nebulosity, and UV sources. The identification of tens of background giants enables us also to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star-formation history of the region, with Trumpler~14 stars found to be systematically younger than stars in other sub-clusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range.

The radial metallicity distribution in the Galactic thin disc represents a crucial constraint for modelling disc formation and evolution. Open clusters allow us to derive both the radial metallicity distribution and its evolution over time. In this paper we perform the first investigation of the present-day radial metallicity distribution based on [Fe/H] determinations in late type members of pre-main-sequence clusters. Because of their youth, these clusters are therefore essential for tracing the current inter-stellar medium metallicity. We used the products of the Gaia-ESO Survey analysis of 12 young regions (age<100 Myr), covering Galactocentric distances from 6.67 to 8.70 kpc. For the first time, we derived the metal content of star forming regions farther than 500 pc from the Sun. Median metallicities were determined through samples of reliable cluster members. For ten clusters the membership analysis is discussed in the present paper, while for other two clusters (Chamaeleon I and Gamma Velorum) we adopted the members identified in our previous works. All the pre-main-sequence clusters considered in this paper have close-to-solar or slightly sub-solar metallicities. The radial metallicity distribution traced by these clusters is almost flat, with the innermost star forming regions having [Fe/H] values that are 0.10-0.15 dex lower than the majority of the older clusters located at similar Galactocentric radii. This homogeneous study of the present-day radial metallicity distribution in the Galactic thin disc favours models that predict a flattening of the radial gradient over time. On the other hand, the decrease of the average [Fe/H] at young ages is not easily explained by the models. Our results reveal a complex interplay of several processes (e.g. star formation activity, initial mass function, supernova yields, gas flows) that controlled the recent evolution of the Milky Way.

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.

Context. It is generally accepted that, in Classical T Tauri Stars, the plasma from the circumstellar disc accretes onto the stellar surface with free fall velocity, and the impact generates a shock. The impact region is expected to contribute to emission in different spectral bands; many studies have confirmed that the X-rays arise from the post-shock plasma but, otherwise, there are no studies in the literature investigating the origin of the observed UV emission which is apparently correlated to accretion. Aims. We investigated the effect of radiative heating of the infalling material by the post-shock plasma at the base of the accretion stream with the aim to identify in which region a significant part of the UV emission originates. Methods. We developed a 1D hydrodynamic model describing the impact of an accretion stream onto the stellar surface; the model takes into account the gravity, the radiative cooling of an optically thin plasma, the thermal conduction, and the heating due to absorption of X-ray radiation. The latter term represents the heating of the infalling plasma due to the absorption of X-rays emitted from the post-shock region. Results. We found that the radiative heating of the pre-shock plasma plays a non-negligible role in the accretion phenomenon. In particular, the dense and cold plasma of the pre-shock accretion column is gradually heated up to few $10^5$\;K due to irradiation of X-rays arising from the shocked plasma at the impact region. This heating mechanism does not affect significantly the dynamics of the post-shock plasma. On the other hand, a region of radiatively heated gas (that we consider a precursor) forms in the unshocked accretion column and contributes significantly to UV emission. Our model naturally reproduces the luminosity of UV emission lines correlated to accretion and shows that most of the UV emission originates from the precursor.

X-ray observations of protostellar jets show evidence of strong shocks heating the plasma up to temperatures of a few million degrees. In some cases, the shocked features appear to be stationary. They are interpreted as shock diamonds. We aim at investigating the physics that guides the formation of X-ray emitting stationary shocks in protostellar jets, the role of the magnetic field in determining the location, stability, and detectability in X-rays of these shocks, and the physical properties of the shocked plasma. We performed a set of 2.5-dimensional magnetohydrodynamic numerical simulations modelling supersonic jets ramming into a magnetized medium and explored different configurations of the magnetic field. The model takes into account the most relevant physical effects, namely thermal conduction and radiative losses. We compared the model results with observations, via the emission measure and the X-ray luminosity synthesized from the simulations. Our model explains the formation of X-ray emitting stationary shocks in a natural way. The magnetic field collimates the plasma at the base of the jet and forms there a magnetic nozzle. After an initial transient, the nozzle leads to the formation of a shock diamond at its exit which is stationary over the time covered by the simulations (~ 40 - 60 yr; comparable with time scales of the observations). The shock generates a point-like X-ray source located close to the base of the jet with luminosity comparable with that inferred from X-ray observations of protostellar jets. For the range of parameters explored, the evolution of the post-shock plasma is dominated by the radiative cooling, whereas the thermal conduction slightly affects the structure of the shock.

The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to ascertain the origin of the observed spatial variations in the spectral properties of the X-ray synchrotron emission, and to predict spectral and morphological features of the resulting gamma-ray emission. We performed 3-D magnetohydrodynamic simulations modeling the evolution of SN 1006 and its interaction with the ambient cloud, and explored different model setups. By applying the REMLIGHT code on the model results, we synthesized the synchrotron X-ray emission, and compared it with actual observations, to constrain the parameters of the model. We also synthesized the leptonic and hadronic gamma-ray emission from the models, deriving constraints on the energy content of the hadrons accelerated at the southwestern limb. We found that the impact of the SN 1006 shock front with a uniform cloud with density 0.5 cm^-3 can explain the observed morphology, the azimuthal variations of the cutoff frequency of the X-ray synchrotron emission, and the shock proper motion in the interaction region. Our results show that the current upper limit for the total hadronic energy in the southwestern limb is 2.5e49 erg.

We present observations from the Gaia-ESO Survey in the lines of H$\alpha$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$\alpha$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $\eta$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $\eta$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$\alpha$ are found to arise from expanding dust reflecting the $\eta$ Car spectrum.

Runaway O- and early B-type stars passing throughout the interstellar medium at supersonic velocities and characterized by strong stellar winds may produce bow shocks that can serve as particle acceleration sites. Previous theoretical models predict the production of high energy photons by non-thermal radiative processes, but their efficiency is still debated. We aim to test and explain the possibility of emission from the bow shocks formed by runaway stars traveling through the interstellar medium by using previous theoretical models. We apply our model to AE Aurigae, the first reported star with an X-ray detected bow shock, to BD+43 3654, in which the observations failed in detecting high energy emission, and to the transition phase of a supergiant star in the late stages of its life.From our analysis, we confirm that the X-ray emission from the bow shock produced by AE Aurigae can be explained by inverse Compton processes involving the infrared photons of the heated dust. We also predict low high energy flux emission from the bow shock produced by BD+43 3654, and the possibility of high energy emission from the bow shock formed by a supergiant star during the transition phase from blue to red supergiant.Bow shock formed by different type of runaway stars are revealed as a new possible source of high energy photons in our neighbourhood.

Understanding the properties of young open clusters, such as the Initial Mass Function (IMF), star formation history and dynamic evolution, is crucial to obtain reliable theoretical predictions of the mechanisms involved in the star formation process. We want to obtain a list, as complete as possible, of confirmed members of the young open cluster Gamma Velorum, with the aim of deriving general cluster properties such as the IMF. We used all available spectroscopic membership indicators within the Gaia-ESO public archive together with literature photometry and X-ray data and, for each method, we derived the most complete list of candidate cluster members. Then, we considered photometry, gravity and radial velocities as necessary conditions to select a subsample of candidates whose membership was confirmed by using the lithium and H$\alpha$ lines and X-rays as youth indicators. We found 242 confirmed and 4 possible cluster members for which we derived masses using very recent stellar evolutionary models. The cluster IMF in the mass range investigated in this study shows a slope of $\alpha=2.6\pm0.5$ for $0.5<M/M_\odot <1.3$ and $\alpha=1.1\pm0.4$ for $0.16<M/M_\odot <0.5$ and is consistent with a standard IMF. The similarity of the IMF of the young population around $\gamma^2 $Vel to that in other star forming regions and the field suggests it may have formed through very similar processes.

The Gaia ESO Public Spectroscopic Survey (GES) is providing the astronomical community with high-precision measurements of many stellar parameters including radial velocities (RVs) of stars belonging to several young clusters and star-forming regions. One of the main goals of the young cluster observations is to study of their dynamical evolution and provide insight into their future, revealing if they will eventually disperse to populate the field, rather than evolve into bound open clusters. In this paper we report the analysis of the dynamical state of L1688 in the $\rho$~Ophiuchi molecular cloud using the dataset provided by the GES consortium. We performed the membership selection of the more than 300 objects observed. Using the presence of the lithium absorption and the location in the Hertzspung-Russell diagram, we identify 45 already known members and two new association members. We provide accurate RVs for all 47 confirmed members.A dynamical analysis, after accounting for unresolved binaries and errors, shows that the stellar surface population of L1688 has a velocity dispersion $\sigma \sim$1.14$\pm$0.35 km s$^{-1}$ that is consistent with being in virial equilibrium and is bound with a $\sim$80% probability. We also find a velocity gradient in the stellar surface population of $\sim$1.0 km s$^{-1}$pc$^{-1}$ in the northwest/southeast direction, which is consistent with that found for the pre-stellar dense cores, and we discuss the possibility of sequential and triggered star formation in L1688.

We discuss the properties of H$\alpha$ emission stars across the sample of 22035 spectra from the Gaia-ESO Survey internal data release, observed with the GIRAFFE instrument and largely belonging to stars in young open clusters. Automated fits using two independent Gaussian profiles and a third component that accounts for the nebular emission allow us to discern distinct morphological types of H$\alpha$ line profiles with the introduction of a simplified classification scheme. All in all we find 3765 stars with intrinsic emission and sort their spectra into eight distinct morphological categories: single--component emission, emission blend, sharp emission peaks, double emission, P-Cygni, inverted P-Cygni, self--absorption, and emission in absorption. We have more than one observation for 1430 stars in our sample, thus allowing a quantitative discussion of the degree of variability of H$\alpha$ emission profiles, which is expected for young, active objects. We present a catalogue of stars with properties of their H$\alpha$ emission line profiles, morphological classification, analysis of variability with time and the supplementary information from the SIMBAD, VizieR, and ADS databases. The records in SIMBAD indicate the presence of H$\alpha$ emission for roughly 25% of all stars in our catalogue, while at least 305 of them have already been more thoroughly investigated according to the references in ADS. The most frequently identified morphological categories in our sample of spectra are emission blend (23%), emission in absorption (22%), and self--absorption (16%). Objects with repeated observations demonstrate that our classification into discrete categories is generally stable through time, but categories P-Cygni and self--absorption seem less stable, which is the consequence of discrete classification rules, as well as of the fundamental change in profile shape.

This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF (JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss together the different agents affecting the structure and the evolution of disks, namely accretion, jets and winds. More information on the JEDI group and its activities can be found at \texttt{http://www.oa-roma.inaf.it/irgroup/JEDI}.

We investigate the plausibility of detecting X-ray emission from a stellar jet that impacts against a dense molecular cloud. This scenario may be usual for classical T Tauri stars with jets in dense star-forming complexes. We first model the impact of a jet against a dense cloud by 2D axisymmetric hydrodynamic simulations, exploring different configurations of the ambient environment. Then, we compare our results with XMM-Newton observations of the Herbig-Haro object HH 248, where extended X-ray emission aligned with the optical knots is detected at the edge of the nearby IC 434 cloud. Our simulations show that a jet can produce plasma with temperatures up to 10 MK, consistent with production of X-ray emission, after impacting a dense cloud. We find that jets denser than the ambient medium but less dense than the cloud produce detectable X-ray emission only at the impact onto the cloud. From the exploration of the model parameter space, we constrain the physical conditions (jet density and velocity, cloud density) that reproduce well the intrinsic luminosity and emission measure of the X-ray source possibly associated with HH 248. Thus, we suggest that the extended X-ray source close to HH 248 corresponds to the jet impacting on a dense cloud.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

The nearby (distance~350-400 pc), rich Vela OB2 association, includes $\gamma^2$ Velorum, one of the most massive binaries in the solar neighbourhood and an excellent laboratory for investigating the formation and early evolution of young clusters. Recent Gaia-ESO survey observations have led to the discovery of two kinematically distinct populations in the young (10-15 Myr) cluster immediately surrounding $\gamma^2$ Velorum. Here we analyse the results of Gaia-ESO survey observations of NGC 2547, a 35 Myr cluster located two degrees south of $\gamma^2$ Velorum. The radial velocity distribution of lithium-rich pre-main sequence stars shows a secondary population that is kinematically distinct from and younger than NGC 2547. The radial velocities, lithium absorption lines, and the positions in a colour-magnitude diagram of this secondary population are consistent with those of one of the components discovered around $\gamma^2$ Velorum. This result shows that there is a young, low-mass stellar population spread over at least several square degrees in the Vela OB2 association. This population could have originally been part of a cluster around $\gamma^2$ Velorum that expanded after gas expulsion or formed in a less dense environment that is spread over the whole Vela OB2 region.

We use the fundamental parameters delivered by the GES consortium in the first internal data release to select the members of $\gamma$ Vel and Cha I among the UVES and GIRAFFE spectroscopic observations. A total of 140 $\gamma$ Vel members and 74 Cha I members were studied. We calculated stellar luminosities through spectral energy distributions, while stellar masses were derived by comparison with evolutionary tracks. The spectral subtraction of low-activity and slowly rotating templates, which are rotationally broadened to match the $v\sin i$ of the targets, enabled us to measure the equivalent widths (EWs) and the fluxes in the H$\alpha$ and H$\beta$ lines. The H$\alpha$ line was also used for identifying accreting objects and for evaluating the mass accretion rate ($\dot M_{\rm acc}$). The distribution of $v\sin i$ for the members of $\gamma$ Vel displays a peak at about 10 km s$^{-1}$ with a tail toward faster rotators. There is also some indication of a different $v\sin i$ distribution for the members of its two kinematical populations. Only a handful of stars in $\gamma$ Vel display signatures of accretion, while many more accretors were detected in the younger Cha~I. Accreting and active stars occupy two different regions in a $T_{\rm eff}$-flux diagram and we propose a criterion for distinguishing them. We derive $\dot M_{\rm acc}$ in the ranges $10^{-11}$-$10^{-9} M_\odot$yr$^{-1}$ and $10^{-10}$-$10^{-7} M_\odot$yr$^{-1}$ for $\gamma$ Vel and Cha I accretors, respectively. We find less scatter in the $\dot M_{\rm acc}-M_\star$ relation derived through the H$\alpha$ EWs, when compared to the H$\alpha$ $10\%W$ diagnostics, in agreement with other authors.

The Einstein spontaneous rates (A-coefficients) of Fe^+ lines have been computed by several authors, with results that differ from each other up to 40%. Consequently, models for line emissivities suffer from uncertainties which in turn affect the determination of the physical conditions at the base of line excitation. We provide an empirical determination of the A-coefficient ratios of bright [Fe II] lines, which would represent both a valid benchmark for theoretical computations and a reference for the physical interpretation of the observed lines. With the ESO-VLT X-shooter instrument between 3,000 A, and 24,700 A, we obtained a spectrum of the bright Herbig-Haro object HH1. We detect around 100 [Fe II] lines, some of which with a signal-to-noise ratio > 100. Among these latter, we selected those emitted by the same level, whose de-reddened intensity ratio is a direct function of the Einstein A-coefficient ratios. From the same X-shooter spectrum, we got an accurate estimate of the extinction toward HH1 through intensity ratios of atomic species, HI, recombination lines and H_2 ro-vibrational transitions. We provide seven reliable A-ooefficient ratios between bright [Fe II] lines, which are compared with the literature determinations. In particular, the A-coefficient ratios involving the brightest near-infrared lines (12570A/16440A and 13209A/16440A) are better in agreement with the predictions by Quinet et al. (1996) Relativistic Hartree-Fock model. However, none of the theoretical models predicts A-coefficient ratios in agreement with all our determinations. We also show that literature data of near-infrared intensity ratios better agree with our determinations than with theoretical expectations.

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 study how diffuse interstellar bands (DIBs) measured toward distance-distributed target stars can be used to locate dense interstellar (IS) clouds in the Galaxy and probe a line-of-sight (LOS) kinematical structure, a potential useful tool when gaseous absorption lines are saturated or not available in the spectral range. Cool target stars are numerous enough for this purpose. We have devised automated DIB fitting methods appropriate to cool star spectra and multiple IS components. The data is fitted with a combination of a synthetic stellar spectrum, a synthetic telluric transmission, and empirical DIB profiles. In parallel, stellar distances and extinctions are estimated self-consistently by means of a 2D Bayesian method, from spectroscopically-derived stellar parameters and photometric data. We have analyzed Gaia-ESO Survey (GES) and previously recorded spectra that probe between $\sim$ 2 and 10 kpc long LOS in five different regions of the Milky Way. Depending on the observed spectral intervals, we extracted one or more of the following DIBs: $\lambda\lambda$ 6283.8, 6613.6 and 8620.4. For each field, we compared the DIB strengths with the Bayesian distances and extinctions, and the DIB Doppler velocities with the HI emission spectra. For all fields, the DIB strength and the target extinction are well correlated. In case of targets widely distributed in distance, marked steps in DIBs and extinction radial distance profiles match with each other and broadly correspond to the expected locations of spiral arms. For all fields, the DIB velocity structure agrees with HI emission spectra and all detected DIBs correspond to strong NaI lines. This illustrates how DIBs can be used to locate the Galactic interstellar gas and to study its kinematics at the kpc scale.

As the number of observed brown dwarf outflows is growing it is important to investigate how these outflows compare to the well studied jets from young stellar objects. A key point of comparison is the relationship between outflow and accretion activity and in particular the ratio between the mass outflow and accretion rates ($\dot{M}_{out}$/$\dot{M}_{acc}$). The brown dwarf candidate ISO-ChaI 217 was discovered by our group, as part of a spectro-astrometric study of brown dwarfs, to be driving an asymmetric outflow with the blue-shifted lobe having a position angle of $\sim$ 20$^{\circ}$. The aim here is to further investigate the properties of ISO-ChaI 217, the morphology and kinematics of its outflow, and to better constrain ($\dot{M}_{out}$/$\dot{M}_{acc}$). The outflow is spatially resolved in the $[SII]\lambda \lambda 6716,6731$ lines and is detected out to $\sim$ 1\farcs6 in the blue-shifted lobe and ~ 1" in the red-shifted lobe. The asymmetry between the two lobes is confirmed although the velocity asymmetry is less pronounced with respect to our previous study. Using thirteen different accretion tracers we measure log($\dot{M}_{acc}$) [M$_{sun}$/yr]= -10.6 $\pm$ 0.4. As it was not possible to measure the effect of extinction on the ISO-ChaI 217 outflow $\dot{M}_{out}$ was derived for a range of values of A$_{v}$, up to a value of A$_{v}$ = 2.5 mag estimated for the source extinction. The logarithm of the mass outflow ($\dot{M}_{out}$) was estimated in the range -11.7 to -11.1 for both jets combined. Thus $\dot{M}_{out}$/$\dot{M}_{acc}$ [\Msun/yr] lies below the maximum value predicted by magneto-centrifugal jet launching models. Finally, both model fitting of the Balmer decrements and spectro-astrometric analysis of the H$\alpha$ line show that the bulk of the H I emission comes from the accretion flow.

Context. Recent metallicity determinations in young open clusters and star-forming regions suggest that the latter may be characterized by a slightly lower metallicity than the Sun and older clusters in the solar vicinity. However, these results are based on small statistics and inhomogeneous analyses. The Gaia-ESO Survey is observing and homogeneously analyzing large samples of stars in several young clusters and star-forming regions, hence allowing us to further investigate this issue. Aims. We present a new metallicity determination of the Chamaeleon I star-forming region, based on the products distributed in the first internal release of the Gaia-ESO Survey. Methods. 48 candidate members of Chamaeleon I have been observed with the high-resolution spectrograph UVES. We use the surface gravity, lithium line equivalent width and position in the Hertzsprung-Russell diagram to confirm the cluster members and we use the iron abundance to derive the mean metallicity of the region. Results. Out of the 48 targets, we confirm 15 high probability members. Considering the metallicity measurements for 9 of them, we find that the iron abundance of Chamaeleon I is slightly subsolar with a mean value [Fe/H]=-0.08+/-0.04 dex. This result is in agreement with the metallicity determination of other nearby star-forming regions and suggests that the chemical pattern of the youngest stars in the solar neighborhood is indeed more metal-poor than the Sun. We argue that this evidence may be related to the chemical distribution of the Gould Belt that contains most of the nearby star-forming regions and young clusters.

In this paper a comprehensive analysis of VLT / X-Shooter observations of two jet systems, namely ESO-H$\alpha$ 574 a K8 classical T Tauri star and Par-Lup 3-4 a very low mass (0.13~\Msun) M5 star, is presented. Both stars are known to have near-edge on accretion disks. A summary of these first X-shooter observations of jets was given in a 2011 letter. The new results outlined here include flux tables of identified emission lines, information on the morphology, kinematics and physical conditions of both jets and, updated estimates of $\dot{M}_{out}$ / $\dot{M}_{acc}$. Asymmetries in the \eso flow are investigated while the \para jet is much more symmetric. The density, temperature, and therefore origin of the gas traced by the Balmer lines are investigated from the Balmer decrements and results suggest an origin in a jet for \eso while for \para the temperature and density are consistent with an accretion flow. $\dot{M}_{acc}$ is estimated from the luminosity of various accretion tracers. For both targets, new luminosity relationships and a re-evaluation of the effect of reddening and grey extinction (due to the edge-on disks) allows for substantial improvements on previous estimates of $\dot{M}_{acc}$. It is found that log($\dot{M}_{acc}$) = -9.15 $\pm$ 0.45~\Msun yr$^{-1}$ and -9.30 $\pm$ 0.27~\Msun yr$^{-1}$ for \eso and \para respectively. Additionally, the physical conditions in the jets (electron density, electron temperature, and ionisation) are probed using various line ratios and compared with previous determinations from iron lines. The results are combined with the luminosity of the [SII]$\lambda$6731 line to derive $\dot{M}_{out}$ through a calculation of the gas emissivity based on a 5-level atom model.

We investigate the diagnostic capabilities of the iron lines for tracing the physical conditions of the shock-excited gas in jets driven by pre-main sequence stars. We have analyzed the 300-2500 nm X-shooter spectra of two jets driven by the pre-main sequence stars ESO-Halpha 574 and Par-Lup 3-4. Both spectra are very rich in [FeII] lines over the whole spectral range; in addition, lines from [FeIII] are detected in the ESO-H\alpha 574 spectrum. NLTE codes along with codes for the ionization equilibrium are used to derive the gas excitation conditions of electron temperature and density, and fractional ionization. The iron gas-phase abundance is provided by comparing the iron lines emissivity with that of [OI] 630 nm. The [FeII] lines indicate ESO-Halpha 574 jet is, on average, colder (T_e = 9000 K), less dense (n_e = 2 10^4 cm^-3) and more ionized (x_e = 0.7) than the Par-Lup 3-4 jet (T_e = 13000 K, n_e = 6 10^4 cm^-3, x_e < 0.4), even if the existence of a higher density component (n_e = 2 10^5 cm^-3) is probed by the [FeIII] and [FeII] ultra-violet lines. Theoretical models suggest that the shock at work in ESO-Halpha 574 is faster and likely more energetic than the Par-Lup 3-4 shock. This latter feature is confirmed by the high percentage of gas-phase iron measured in ESO-Halpha 574 (50-60% of its solar abundance in comparison with less than 30% in Par-Lup 3-4), which testifies that the ESO-Halpha 574 shock is powerful enough to partially destroy the dust present inside the jet. This work demonstrates that a multiline Fe analysis can be effectively used to probe the excitation and ionization conditions of the gas in a jet without any assumption on ionic abundances. The main limitation on the diagnostics resides in the large uncertainties of the atomic data, which, however, can be overcome through a statistical approach involving many lines.

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

Protostellar jets appear at all stages of star formation when the accretion process is still at work. Jets travel at velocities of hundreds of km/s, creating strong shocks when interacting with interstellar medium. Several cases of jets have been detected in X-rays, typically showing soft emission. For the first time, we report evidence of hard X-ray emission possibly related to non-thermal processes not explained by previous models of the post-shock emission predicted in the jet/ambient interaction scenario. HH 80 is located at the south head of the jet associated to the massive protostar IRAS 18162-2048. It shows soft and hard X-ray emission in regions that are spatially separated, with the soft X-ray emission region situated behind the region of hard X-ray emission. We propose a scenario for HH 80 where soft X-ray emission is associated to thermal processes from the interaction of the jet with denser ambient matter and the hard X-ray emission is produced by synchrotron radiation at the front shock.

Attention is given to a population of 110 stars in the NGC 6611 cluster of the Eagle Nebula that have prominent near-infrared (NIR) excess and optical colours typical of pre-main sequence (PMS) stars older than 8 Myr. At least half of those for which spectroscopy exists have a Halpha emission line profile revealing active accretion. In principle, the V-I colours of all these stars would be consistent with those of young PMS objects (< 1 Myr) whose radiation is heavily obscured by a circumstellar disc seen at high inclination and in small part scattered towards the observer by the back side of the disc. However, using theoretical models it is shown here that objects of this type can only account for a few percent of this population. In fact, the spatial distribution of these objects, their X-ray luminosities, their optical brightness, their positions in the colour-magnitude diagram and the weak Li absorption lines of the stars studied spectroscopically suggest that most of them are at least 8 times older than the ~1 Myr-old PMS stars already known in this cluster and could be as old as ~30 Myr. This is the largest homogeneous sample to date of Galactic PMS stars considerably older than 8 Myr that are still actively accreting from a circumstellar disc and it allows us to set a lower limit of 7% to the disc frequency at ~16 Myr in NGC 6611. These values imply a characteristic exponential lifetime of ~6 Myr for disc dissipation.

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.

The young open cluster NGC 6611 includes among its candidate members a class of peculiar objects with interesting properties: blue stars with infrared IR excesses. These stars show excesses in IR bands, signature of the presence of a circumstellar disk, but optical colors typical of older field stars. In order to confirm their membership to the cluster, it is therefore important to use new spectroscopic observations, together with previous photometric data. We aim at confirming the membership of these objects and at investigating their physical properties to verify whether the observed colors are intrinsic or altered by the presence of the disk or by the accretion processes. We analyze the intermediate resolution spectroscopic data obtained for a subsample of blue stars in NGC 6611 with FLAMES. In particular, we focus on the study of: 1) the profile of the Ha emission line, to select stars with accretion and outflow activity; 2) the Li absorption line, used as a youth indicator; 3) the radial velocity. Using the spectroscopic analysis, it has been possible to investigate the presence of the Li absorption line as well as to discriminate between stars with inert or active disk. In particular, from the analysis of the Ha emission line we were able to infer the activity due to the accretion and outflow processes and the variability of the emission. We also investigated the binarity of the blue stars and their membership to NGC 6611.

Runaway stars produce shocks when passing through interstellar medium at supersonic velocities. Bow shocks have been detected in the mid-infrared for several high-mass runaway stars and in radio waves for one star. Theoretical models predict the production of high-energy photons by non-thermal radiative processes in a number sufficiently large to be detected in X-rays. To date, no stellar bow shock has been detected at such energies. We present the first detection of X-ray emission from a bow shock produced by a runaway star. The star is AE Aur, which was likely expelled from its birthplace by the encounter of two massive binary systems and now is passing through the dense nebula IC 405. The X-ray emission from the bow shock is detected at 30" to the northeast of the star, coinciding with an enhancement in the density of the nebula. From the analysis of the observed X-ray spectrum of the source and our theoretical emission model, we confirm that the X-ray emission is produced mainly by inverse Compton upscattering of infrared photons from dust in the shock front.

X-ray emission from about ten protostellar jets has been discovered and it appears as a feature common to the most energetic jets. Although X-ray emission seems to originate from shocks internal to jets, the mechanism forming these shocks remains controversial. One of the best studied X-ray jet is HH 154 that has been observed by Chandra over a time base of about 10 years. We analyze the Chandra observations of HH 154 by investigating the evolution of its X-ray source. We show that the X-ray emission consists of a bright stationary component and a faint elongated component. We interpret the observations by developing a hydrodynamic model describing a protostellar jet originating from a nozzle and compare the X-ray emission synthesized from the model with the X-ray observations. The model takes into account the thermal conduction and radiative losses and shows that the jet/nozzle leads to the formation of a diamond shock at the nozzle exit. The shock is stationary over the period covered by our simulations and generates an X-ray source with luminosity and spectral characteristics in excellent agreement with the observations. We conclude that the X-ray emission from HH 154 is consistent with a diamond shock originating from a nozzle through which the jet is launched into the ambient medium. We suggest that the physical origin of the nozzle could be related to the dense gas in which the HH 154 driving source is embedded and/or to the magnetic field at the jet launching/collimation region.

Protostellar jets are known to emit in a wide range of bands, from radio to IR to optical bands, and to date also about ten X-ray emitting jets have been detected, with a rate of discovery of about one per year. We aim at investigating the mechanism leading to the X-ray emission detected in protostellar jets and at constraining the physical parameters that describe the jet/ambient interaction by comparing our model predictions with observations. We perform 2D axisymmetric hydrodynamic simulations of the interaction between a supersonic jet and the ambient. The jet is described as a train of plasma blobs randomly ejected by the stellar source along the jet axis. We explore the parameter space by varying the ejection rate, the initial jet Mach number, and the initial density contrast between the ambient and the jet. We synthesized from the model the X-ray emission as it would be observed with the current X-ray telescopes. The mutual interactions among the ejected blobs and of the blobs with the ambient medium lead to complex X-ray emitting structures within the jet: irregular chains of knots; isolated knots with measurable proper motion; apparently stationary knots; reverse shocks. The predicted X-ray luminosity strongly depends on the ejection rate and on the initial density contrast between the ambient and the jet, with a weaker dependence on the jet Mach number. Our model represents the first attempt to describe the X-ray properties of all the X-ray emitting protostellar jets. The comparison between our model predictions and the observations can provide a useful diagnostic tool necessary for a proper interpretation of the observations. In particular, we suggest that the observable quantities derived from the spectral analysis of X-ray observations can be used to constrain the ejection rate, a parameter explored in our model that is not measurable by current observations.

HH objects are characterized by a complex knotty morphology detected mainly along the axis of protostellar jets in a wide range of bands. Evidence of interactions between knots formed in different epochs have been found, suggesting that jets may result from the ejection of plasma blobs from the source. We aim at investigating the physical mechanism leading to the irregular knotty structure observed in jets in different bands and the complex interactions occurring among blobs of plasma ejected from the stellar source. We perform 2D axisymmetric HD simulations of a randomly ejected pulsed jet. The jet consists of a train of blobs which ram with supersonic speed into the ambient medium. The initial random velocity of each blob follows an exponential distribution. We explore the ejection rate parameter to derive constraints on the physical properties of protostellar jets by comparison of model results with observations. Our model takes into account radiative losses and thermal conduction. We find that the mutual interactions of blobs ejected at different epochs and with different speed lead to a variety of plasma components not described by current models. The main features characterizing the random pulsed jet scenario are: single high speed knots, showing a measurable proper motion in nice agreement with observations; irregular chains of knots aligned along the jet axis and possibly interacting with each other; reverse shocks interacting with outgoing knots; oblique shocks produced by the reflection of shocks at the jet cocoon. All these structures concur to determine the morphology of the jet in different bands. We also find that the thermal conduction plays a crucial role in damping out HD instabilities that would develop within the cocoon and that contribute to the jet breaking.

Recently X-ray emission from protostellar jets has been detected with both XMM-Newton and Chandra satellites, but the physical mechanism which can give rise to this emission is still unclear. We performed an extensive exploration of a wide space of the main parameters influencing the jet/ambient interaction. Aims include: 1) to constrain the jet/ambient interaction regimes leading to the X-ray emission observed in Herbig-Haro objects in terms of the emission by a shock forming at the interaction front between a continuous supersonic jet and the surrounding medium; 2) to derive detailed predictions to be compared with optical and X-ray observations of protostellar jets; 3) to get insight into the protostellar jet's physical conditions. We performed a set of bidimensional hydrodynamic numerical simulations, in cylindrical coordinates, modeling supersonic jets ramming a uniform ambient medium. The model takes into account the most relevant physical effects, namely the thermal conduction and the radiative losses. Our model explains the observed X-ray emission from protostellar jets in a natural way. In particular we find that the case of a protostellar jet less dense than the ambient medium reproduces well the observations of the nearest Herbig-Haro object, HH154, and allows us to make detailed predictions of a possible X-ray source proper motion (vsh = 500 km/s), detectable with Chandra. Furthermore our results suggest that the simulated protostellar jets which best reproduce the X-rays observations cannot drive molecular outflows.

We have performed, in October 2005, a deep Chandra X-ray observation of HH 154. Comparison with the previous (2001) Chandra observation allows to detect proper motion down to the level predicted by models of X-ray emitting shocks in the jet. The 2005 Chandra observation of HH 154 shows unexpected morphological changes of the X-ray emission in comparison with the 2001 data. Two components are present: a stronger, point-like component with no detectable motion and a weaker component which has expanded in size by approximately 300 AU over the 4 years time base of the two observations. This expansion corresponds to approximately 500 km/s, very close to the velocity of the X-ray emitting shock in the simple theoretical models. The 2005 data show a more complex system than initially thought (and modeled), with multiple components with different properties. The observed morphology is possibly indicating a pulsed jet propagating through a non-homogeneous medium, likely with medium density decreasing with distance from the driving source. Detailed theoretical modeling and deeper X-ray observations will be needed to understand the physics of this fascinating class of sources.