Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Aims: The inner regions of high-mass protostars are often invisible in the near-infrared. We aim to investigate the inner gaseous disc of IRAS11101-5829 through scattered light from the outflow cavity walls. Methods: We observed the environment of the high-mass young stellar object IRAS11101-5829 and the closest knots of its jet, HH135-136, with the VLT/SINFONI. We also retrieved archival data from the high-resolution long-slit spectrograph VLT/X-shooter. Results: We detect the first three bandheads of the $\upsilon=2-0$ CO vibrational emission for the first time in this object. It is coincident with continuum and Br$\gamma$ emission and extends up to $\sim10000$ au to the north-east and $\sim10 000$ au to the south-west. The line profiles have been modelled as a Keplerian rotating disc assuming a single ring in LTE. The model output gives a temperature of $\sim3000$ K, a CO column density of $\sim1\times10^{22}\mathrm{ cm^{-2}}$, and a projected Keplerian velocity $v_\mathrm{K}\sin i_\mathrm{disc} \sim 25\mathrm{ km s^{-1}}$, which is consistent with previous modelling in other high-mass protostars. In particular, the low value of $v_\mathrm{K}\sin i_\mathrm{disc}$ suggests that the disc is observed almost face-on, whereas the well-constrained geometry of the jet imposes that the disc must be close to edge-on. This apparent discrepancy is interpreted as the CO seen reflected in the mirror of the outflow cavity wall. Conclusions: From both jet geometry and disc modelling, we conclude that all the CO emission is seen through reflection by the cavity walls and not directly. This result implies that in the case of highly embedded objects, as for many high-mass protostars, line profile modelling alone might be deceptive and the observed emission could affect the derived physical and geometrical properties; in particular the inclination of the system can be incorrectly interpreted.

Numerous eruptive variable young stellar objects (YSOs), mostly Class I systems, were recently detected by the near-infrared Vista Variables in the Via Lactea (VVV) survey. We present an exploratory near-infrared spectroscopic variability study of 14 eruptive YSOs. The variations were sampled over 1-day and 1 to 2-year intervals and analysed in combination with VVV light curves. CO overtone absorption features are observed on 3 objects with FUor-like spectra: all show deeper absorption when they are brighter. This implies stronger emission from the circumstellar disc with a steeper vertical temperature gradient when the accretion rate is higher. This confirms the nature of fast VVV FUor-like events, in line with the accepted picture for classical FUors. The absence of Br$\gamma$ emission in a FUor-like object declining to pre-outburst brightness suggests that reconstruction of the stellar magnetic field is a slow process. Within the 1-day timescale, 60% of H$_2$-emitting YSOs show significant but modest variation, and 2/6 sources have large variations in Br$\gamma$. Over year-long timescales, H$_2$ flux variations remain modest despite up to 1.8 mag variation in $K_s$. This indicates that emission from the molecular outflow usually arises further from the protostar and is unaffected by relatively large changes in accretion rate on year-long timescales. Two objects show signs of on/off magnetospheric accretion traced by Br$\gamma$ emission. In addition, a 60% inter-night brightening of the H$_2$ outflow is detected in one YSO.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded.

It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements.

RW Aur A is a CTTS that has suddenly undergone three major dimming events since 2010. We aim to understand the dimming properties, examine accretion variability, and derive the physical properties of the inner disc traced by the CO ro-vibrational emission at NIR wavelengths (2.3 mic). We compared two epochs of X-Shooter observations, during and after the dimming. We modelled the rarely detected CO bandhead emission in both epochs to examine whether the inner disc properties had changed. The SED was used to derive the extinction properties of the dimmed spectrum and compare the infrared excess between the two epochs. Lines tracing accretion were used to derive the mass accretion rate in both states. The CO originates from a region with physical properties of T=3000 K, N$_{CO}$=1x10$^{21}$ cm$^{-2}$ and vsini=113 km/s. The extinction properties of the dimming layer were derived with the effective optical depth ranging from teff 2.5-1.5 from the UV to the NIR. The inferred mass accretion rate Macc is $1.5x 10^{-8}$ Msun/yr and $\sim 2x 10^{-8}$ Msun/yr after and during the dimming respectively. By fitting the SED, additional emission is observed in the IR during the dimming event from dust grains with temperatures of 500-700K. The physical conditions traced by the CO are similar for both epochs, indicating that the inner gaseous disc properties do not change during the dimming events. The extinction curve is flatter than that of the ISM, and large grains of a few hundred microns are thus required. When we correct for the observed extinction, Macc is constant in the two epochs, suggesting that the accretion is stable and therefore does not cause the dimming. The additional hot emission in the NIR is located at about 0.5 au from the star. The dimming events could be due to a dust-laden wind, a severe puffing-up of the inner rim, or a perturbation caused by the recent star-disc encounter.

Binarity and multiplicity appear to be a common outcome in star formation. In particular, the binary fraction of massive (OB-type) stars can be very high. In many cases, the further stellar evolution of these stars is affected by binary interactions at some stage during their lifetime. The origin of this high binarity and the binary parameters are poorly understood because observational constraints are scarce, which is predominantly due to a dearth of known young massive binary systems. We aim to identify and describe massive young binary systems in order to fill in the gaps of our knowledge of primordial binarity of massive stars, which is crucial for our understanding of massive star formation. We observed the two massive young stellar objects (MYSOs) PDS 27 and PDS 37 at the highest spatial resolution provided by VLTI/PIONIER in the H-band (1.3 mas). We applied geometrical models to fit the observed squared visibilities and closure phases. In addition, we performed a radial velocity analysis using published VLT/FORS2 spectropolarimetric and VLT/X-shooter spectroscopic observations. Our findings suggest binary companions for both objects at 12 mas (30 au) for PDS 27 and at 22-28 mas (42-54 au) for PDS 37. This means that they are among the closest MYSO binaries resolved to date. Our data spatially resolve PDS 27 and PDS 37 for the first time, revealing two of the closest and most massive ($>$8 M$_\odot$) YSO binary candidates to date. PDS 27 and PDS 37 are rare but great laboratories to quantitatively inform and test the theories on formation of such systems.

Centimeter continuum observations of protostellar jets have revealed the presence of knots of shocked gas where the flux density decreases with frequency. This spectrum is characteristic of nonthermal synchrotron radiation and implies the presence of both magnetic fields and relativistic electrons in protostellar jets. Here, we report on one of the few detections of nonthermal jet driven by a young massive star in the star-forming region G035.02$+$0.35. We made use of the NSF's Karl G. Jansky Very Large Array (VLA) to observe this region at C, Ku, and K bands with the A- and B-array configurations, and obtained sensitive radio continuum maps down to a rms of 10 $\mu$Jy beam$^{-1}$. These observations allow for a detailed spectral index analysis of the radio continuum emission in the region, which we interpret as a protostellar jet with a number of knots aligned with extended 4.5 $\mu$m emission. Two knots clearly emit nonthermal radiation and are found at similar distances, of approximately 10,000 au, each side of the central young star, from which they expand at velocities of hundreds km s$^{-1}$. We estimate both the mechanical force and the magnetic field associated with the radio jet, and infer a lower limit of $0.4\times10^{-4} $M$_{\odot}$ yr$^{-1}$ km s$^{-1}$ and values in the range $0.7-1.3 $mG, respectively.

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

The VVV survey has allowed for an unprecedented number of multi-epoch observations of the southern Galactic plane. In a recent paper,13 massive young stellar objects(MYSOs) have already been identified within the highly variable(\Delta Ks > 1 mag) YSO sample of another published work.This study aims to understand the general nature of variability in MYSOs.We present the first systematic study of variability in a large sample of candidate MYSOs.We examined the data for variability of the putative driving sources of all known Spitzer EGOs and bright 24 \mu m sources coinciding with the peak of 870 \mu m detected ATLASGAL clumps, a total of 718 targets. Of these, 190 point sources (139 EGOs and 51 non-EGOs) displayed variability (IQR > 0.05, \Delta Ks > 0.15 mag). Light-curves(LCs) have been sub-classified into eruptive, dipper, fader, short-term-variable and long-period variable-YSO categories.Lomb-Scargle periodogram analysis of periodic LCs was carried out. 1 - 870 \mu m spectral energy distributions of the variable sources were fitted with YSO models to obtain representative properties. 41% of the variable sources are represented by > 4Msun objects, and only 6% by > 8Msun objects.The highest-mass objects are mostly non-EGOs,deeply embedded.By placing them on the HR diagram we show that most lower mass,EGO type objects are concentrated on the putative birth-line position, while the luminous non-EGO type objects group around the ZAMS track.Some of the most luminous far infrared sources in the massive clumps and infrared quiet driving sources of EGOs have been missed out by this study owing to an uniform sample selection method.A high rate of detectable variability in EGO targets (139 out of 153 searched) implies that near-infrared variability in MYSOs is closely linked to the accretion phenomenon and outflow activity.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

This short review examines the most recent observational results obtained for jets driven by high-mass young stellar objects (HMYSOs). Reviewed data range from near-infrared to radio wavelengths, with particular emphasis on the IR regime, and are presented along with perspectives for future observations achievable with newly available or forthcoming facilities.

Context: Protostellar jets in high-mass young stellar objects (HMYSOs) play a key role in the understanding of star formation and provide us with an excellent tool to study fundamental properties of HMYSOs. Aims: We aim at studying the physical and kinematic properties of the near-IR (NIR) jet of IRAS\,13481-6124 from au to parsec scales. Methods: Our study includes NIR data from the Very Large Telescope instruments SINFONI, CRIRES, and ISAAC. Information about the source and its immediate environment is retrieved with SINFONI. The technique of spectro-astrometry is performed with CRIRES to study the jet on au scales. The parsec-scale jet and its kinematic and dynamic properties are investigated using ISAAC. Results: The SINFONI spectra in the $H$ and $K$ band are rich in emission lines that are mainly associated with ejection and accretion processes. Spectro-astrometry is applied to the Br$\gamma$ line, and for the first time, to the Br$\alpha$ line, revealing their jet origin with milliarcsecond-scale photocentre displacements ($11-15$\,au). This allows us to constrain the kinematics of the au-scale jet and to derive its position angle ($\sim216^{\circ}$). ISAAC spectroscopy reveals H$_2$ emission along the parsec-scale jet, which allows us to infer kinematic and dynamic properties of the NIR parsec-scale jet. The mass-loss rate inferred for the NIR jet is $\dot{M}_\mathrm{ejec}\sim10^{-4}\mathrm{\,M_\odot\,yr^{-1}}$ and the thrust is $\dot{P}\sim10^{-2}\mathrm{\,M_\odot\,yr^{-1}\,km\,s^{-1}}$ , which is roughly constant for the formation history of the young star. A tentative estimate of the ionisation fraction is derived for the massive jet by comparing the radio and NIR mass-loss rates. An ionisation fraction $\lesssim8\%$ is obtained, which means that the bulk of the ejecta is traced by the NIR jet and that the radio jet only delineates a small portion of it.

HD50138 is a Herbig B[e] star with a circumstellar disc detected at IR and mm wavelength. Its brightness makes it a good candidate for NIR interferometry observations. We aim to resolve, spatially and spectrally, the continuum and hydrogen emission lines in the 2.12-2.47 micron region, to shed light on the immediate circumstellar environment of the star. VLTI/AMBER K-band observations provide spectra, visibilities, differential phases, and closure phases along three long baselines for the continuum, and HI emission in Br$\gamma$ and five high-n Pfund lines. By computing the pure-line visibilities, we derive the angular size of the different line-emitting regions. A simple LTE model was created to constrain the physical conditions of HI emitting region. The continuum region cannot be reproduced by a geometrical 2D elongated Gaussian fitting model. We estimate the size of the region to be 1 au. We find the Br$\gamma$ and Pfund lines come from a more compact region of size 0.4 au. The Br$\gamma$ line exhibits an S-shaped differential phase, indicative of rotation. The continuum and Br$\gamma$ line closure phase show offsets of $\sim$-25$\pm$5 $^o$ and 20$\pm$10$^o$, respectively. This is evidence of an asymmetry in their origin, but with opposing directions. We find that we cannot converge on constraints for the HI physical parameters without a more detailed model. Our analysis reveals that HD50138 hosts a complex circumstellar environment. Its continuum emission cannot be reproduced by a simple disc brightness distribution. Similarly, several components must be evoked to reproduce the interferometric observables within the Br$\gamma$, line. Combining the spectroscopic and interferometric data of the Br$\gamma$ and Pfund lines favours an origin in a wind region with a large opening angle. Finally, our results point to an evolved source.

Context. Recent observations of the massive young stellar object S255 NIRS 3 have revealed a large increase in both methanol maser flux density and IR emission, which have been interpreted as the result of an accretion outburst, possibly due to instabilities in a circumstellar disk. This indicates that this type of accretion event could be common in young/forming early-type stars and in their lower mass siblings, and supports the idea that accretion onto the star may occur in a non-continuous way. Aims. As accretion and ejection are believed to be tightly associated phenomena, we wanted to confirm the accretion interpretation of the outburst in S255 NIRS 3 by detecting the corresponding burst of the associated thermal jet. Methods. We monitored the radio continuum emission from S255 NIRS 3 at four bands using the Karl G. Jansky Very Large Array. The millimetre continuum emission was also observed with both the Northern Extended Millimeter Array of IRAM and the Atacama Large Millimeter/submillimeter Array. Results. We have detected an exponential increase in the radio flux density from 6 to 45 GHz starting right after July 10, 2016, namely about 13 months after the estimated onset of the IR outburst. This is the first ever detection of a radio burst associated with an IR accretion outburst from a young stellar object. The flux density at all observed centimetre bands can be reproduced with a simple expanding jet model. At millimetre wavelengths we infer a marginal flux increase with respect to the literature values and we show this is due to free-free emission from the radio jet. Abridged.

There is growing evidence for disk-mediated accretion being the dominant mode of star formation across nearly the whole stellar mass spectrum. The stochastic nature of this process has been realized which implies an inherent source variability. It can be traced more easily for low-mass YSOs (LMYSOs) since high-mass YSOs (HMYSOs) are still embedded even when reaching the ZAMS. While variable reflection nebulae around LMYSOs were among the earliest signs of star formation, little is known on the variability of scattered light from embedded clusters, the birthplaces of HMYSOs. Since the few most massive stars dominate this emission, their variability is literally reflected in scattered light. Moreover, because of their high luminosity, for a given ambient dust density and source distance, the associated nebulosities are much larger than those of LMYSOs. In this case, the light travel time becomes substantial. So the apparent brightness distribution constitutes a light echo, shaped by both the HMYSO variability history and the spatial distribution of the scattering medium. We report on early results of a NIR variability study of HMYSOs associated with Class II methanol masers which aims at revealing a possible correlation between maser flux density and infrared brightness. Additionally, relevant findings for the eruptive HMYSO S255IR-NIRS3 are presented.

Methanol and water masers indicate young stellar objects. They often exhibit flares, and a fraction shows periodic activity. Several mechanisms might explain this behavior but the lack of concurrent infrared (IR) data complicates to identify the cause. Recently, 6.7 GHz methanol maser flares were observed, triggered by accretion bursts of high-mass YSOs which confirmed the IR-pumping of these masers. This suggests that regular IR changes might lead to maser periodicity. Hence, we scrutinized space-based IR imaging of YSOs associated with periodic methanol masers. We succeeded to extract the IR light curve from NEOWISE data for the intermediate mass YSO G107.298+5.639. Thus, for the first time a relationship between the maser and IR variability could be established. While the IR light curve shows the same period of ~34.6 days as the masers, its shape is distinct from that of the maser flares. Possible reasons for the IR periodicity are discussed.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.

Aims: We want to study the physical properties of the ionized jet emission in the vicinity of an O-type young stellar object (YSO), and estimate how efficient is the transfer of energy and momentum from small- to large-scale outflows. Methods: We conducted Karl G. Jansky Very Large Array (VLA) observations, at both 22 and 45 GHz, of the compact and faint radio continuum emission in the high-mass star-forming region G023.01-00.41, with an angular resolution between 0.3" and 0.1", and a thermal rms of the order of 10 uJy/beam. Results: We discovered a collimated thermal (bremsstrahlung) jet emission, with a radio luminosity (L_rad) of 24 mJy kpc^2 at 45 GHz, in the inner 1000 AU from an O-type YSO. The radio thermal jet has an opening angle of 44 degrees and brings a momentum rate of 8 10^-3 M_sun yr^-1 km/s. By combining the new data with previous observations of the molecular outflow and water maser shocks, we can trace the outflow emission from its driving source through the molecular clump, across more than two order of magnitude in length (500 AU-0.2 pc). We find that the momentum-transfer efficiency, between the inner jet emission and the extended outflow of entrained ambient gas, is near unity. This result suggests that the large-scale flow is swept-up by the mechanical force of the radio jet emission, which originates in the inner 1000 AU from the high-mass YSO.

We present the discovery of 816 high amplitude infrared variable stars ($\Delta K_{\rm s} >$ 1 mag) in 119 deg$^2$ of the Galactic midplane covered by the Vista Variables in the Via Lactea (VVV) survey. Almost all are new discoveries and about 50$\%$ are YSOs. This provides further evidence that YSOs are the commonest high amplitude infrared variable stars in the Galactic plane. In the 2010-2014 time series of likely YSOs we find that the amplitude of variability increases towards younger evolutionary classes (class I and flat-spectrum sources) except on short timescales ($<$25 days) where this trend is reversed. Dividing the likely YSOs by light curve morphology, we find 106 with eruptive light curves, 45 dippers, 39 faders, 24 eclipsing binaries, 65 long-term periodic variables (P$>$100 days) and 162 short-term variables. Eruptive YSOs and faders tend to have the highest amplitudes and eruptive systems have the reddest SEDs. Follow up spectroscopy in a companion paper verifies high accretion rates in the eruptive systems. Variable extinction is disfavoured by the 2 epochs of colour data. These discoveries increase the number of eruptive variable YSOs by a factor of at least 5, most being at earlier stages of evolution than the known FUor and EXor types. We find that eruptive variability is at least an order of magnitude more common in class I YSOs than class II YSOs. Typical outburst durations are 1 to 4 years, between those of EXors and FUors. They occur in 3 to 6\% of class I YSOs over a 4 year time span.

In a companion work (Paper I) we detected a large population of highly variable Young Stellar Objects (YSOs) in the Vista Variables in the Via Lactea (VVV) survey, typically with class I or flat spectrum spectral energy distributions and diverse light curve types. Here we present infrared spectra (0.9--2.5 $\mu$m) of 37 of these variables, many of them observed in a bright state. The spectra confirm that 15/18 sources with eruptive light curves have signatures of a high accretion rate, either showing EXor-like emission features ($\Delta$v=2 CO, Br$\gamma$) and/or FUor-like features ($\Delta$v=2 CO and H$_{2}$O strongly in absorption). Similar features were seen in some long term periodic YSOs and faders but not in dippers or short-term variables. The sample includes some dusty Mira variables (typically distinguished by smooth Mira-like light curves), 2 cataclysmic variables and a carbon star. In total we have added 19 new objects to the broad class of eruptive variable YSOs with episodic accretion. Eruptive variable YSOs in our sample that were observed at bright states show higher accretion luminosities than the rest of the sample. Most of the eruptive variables differ from the established FUor and EXor subclasses, showing intermediate outburst durations and a mixture of their spectroscopic characteristics. This is in line with a small number of other recent discoveries. Since these previously atypical objects are now the majority amongst embedded members of the class, we propose a new classification for them as MNors. This term (pronounced emnor) follows V1647 Ori, the illuminating star of McNeil's Nebula.

During their formation phase stars gain most of their mass in violent episodic accretion events, such as observed in FU Orionis (FUor) and EXor stars. V346 Normae is a well-studied FUor that underwent a strong outburst beginning in $\sim1980$. Here, we report photometric and spectroscopic observations which show that the visual/near-infrared brightness has decreased dramatically between the 1990s and 2010 (${\Delta}R\approx10.9^{\rm m}$, ${\Delta}J\approx7.8^{\rm m}$, ${\Delta}K\approx5.8^{\rm m}$). The spectral properties of this fading event cannot be explained with variable extinction alone, but indicate a drop in accretion rate by 2-3 orders of magnitude, marking the first time that a member of the FUor class has been observed to switch to a very low accretion phase. Remarkably, in the last few years (2011-2015) V346 Nor has brightened again at all near-infrared wavelengths, indicating the onset of a new outburst event. The observed behaviour might be consistent with the clustered luminosity bursts that have been predicted by recent gravitational instability and fragmentation models for the early stages of protostellar evolution. Given V346 Nor's unique characteristics (concerning outburst duration, repetition frequency, and spectroscopic diagnostics), our results also highlight the need for revisiting the FUor/EXor classification scheme.

To probe the circumstellar environment of IRAS 13481-6124, a 20 M_sun high-mass young stellar object (HMYSO) with a parsec-scale jet and accretion disc, we investigate the origin of its Br\gamma-emission line through NIR interferometry. We present the first AMBER/VLTI observations of the Br\gamma-emitting region in an HMYSO at R~1500. Our AMBER/VLTI observations reveal a spatially and spectrally resolved Br\gamma-line in emission with a strong P Cygni profile, indicating outflowing matter with a terminal velocity of ~500 km/s. Visibilities, differential phases, and closure phases are detected in our observations within the spectral line and in the adjacent continuum. Both total visibilities (continuum plus line emitting region) and pure-line visibilities indicate that the Br\gamma-emitting region is more compact (2-4 mas in diameter or ~6-13 au at 3.2 kpc) than the continuum-emitting region (~5.4 mas or ~17 au). The absorption feature is also spatially resolved at the longest baselines (81 and 85 m) and has a visibility that is slightly smaller than the continuum-emitting region. The differential phases at the four longest baselines display an \u2018S\u2019-shaped structure across the line, peaking in the blue- and red-shifted high-velocity components. The calibrated photocentre shifts are aligned with the known jet axis, i.e they are probably tracing an ionised jet. The high-velocity components (v_r~100-500 km/s) are located far from the source, whereas the low-velocity components (0-100 km/s) are observed to be closer, indicating a strong acceleration of the gas flow in the inner 10 au. Finally, a non-zero closure phase along the continuum is detected. By comparing our observations with the synthetic images of the continuum around 2.16 um, we confirm that this feature originates from the asymmetric brightness distribution of the continuum owing to the inclination of the inner disc.

The origin of the near-infrared (NIR) HI emission lines in young stellar objects are not yet understood. To probe it, we present multi-epoch LBT-LUCIFER spectroscopic observations of the Pa{\delta}, Pa{\beta}, and Br{\gamma} lines observed in the Herbig star VVSer, along with VLTI-AMBER Br{\gamma} spectro-interferometric observations at medium resolution. Our spectroscopic observations show line profile variability in all the HI lines. The strongest variability is observed in the redshifted part of the line profiles. The Br{\gamma} spectro-interferometric observations indicate that the Br{\gamma} line emitting region is smaller than the continuum emitting region. To interpret our results, we employed radiative transfer models with three different flow configurations: magnetospheric accretion, a magneto-centrifugally driven disc wind, and a schematic bipolar outflow. Our models suggest that the HI line emission in VVSer is dominated by the contribution of an extended wind, perhaps a bipolar outflow. Although the exact physical process for producing such outflow is not known, this model is capable of reproducing the averaged single-peaked line profiles of the HI lines. Additionally, the observed visibilities, differential and closure phases are best reproduced when a wind is considered. Nevertheless, the complex line profiles and variability could be explained by changes in the relative contribution of the magnetosphere and/or winds to the line emission. This might indicate that the NIR HI lines are formed in a complex inner disc region where inflow and outflow components might coexist. Furthermore, the contribution of each of these mechanisms to the line appears time variable, suggesting a non-steady accretion/ejection flow.

We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object (YSO), at very high spatial resolution (10-100 AU). We performed milli-arcsecond polarimetric observations of the strong CH3OH maser emission observed in the vicinity of an O-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH3OH masing gas previously measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2 degree/AU). Overall, the velocity field vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30 degrees. We interpret this finding as the contribution of a turbulent velocity field of about 3.5 km/s, responsible for braking up the alignment between the velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk plane, with an average speed of 7 km/s. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 500-600 AU, and inward for shorter radii.

We analyse the main physical parameters and the circumstellar environment of the young Herbig Be star HD 98922. We present AMBER/VLTI high spectral resolution (R =12000) interferometric observations across the Br$\gamma$ line, accompanied by UVES high-resolution spectroscopy and SINFONI-AO assisted near-infrared integral field spectroscopic data. To interpret our observations, we develop a magneto-centrifugally driven disc-wind model. Our analysis of the UVES spectrum shows that HD 98922 is a young (~5x10^5 yr) Herbig Be star (SpT=B9V), located at a distance of 440(+60-50) pc, with a mass accretion rate of ~9+/-3x10^(-7) M_sun yr^(-1). SINFONI K-band AO-assisted imaging shows a spatially resolved circumstellar disc-like region (~140 AU in diameter) with asymmetric brightness distribution. Our AMBER/VLTI UT observations indicate that the Br$\gamma$ emitting region (radius ~0.31+/-0.04 AU) is smaller than the continuum emitting region (inner dust radius ~0.7+/-0.2 AU), showing significant non-zero V-shaped differential phases (i.e. non S-shaped, as expected for a rotating disc). The value of the continuum-corrected pure Br$\gamma$ line visibility at the longest baseline (89 m) is ~0.8+/-0.1, i.e. the Br$\gamma$ emitting region is partially resolved. Our modelling suggests that the observed Br$\gamma$ line-emitting region mainly originates from a disc wind with a half opening angle of 30deg, and with a mass-loss rate of ~2x10(-7) M_sun yr^(-1). The observed V-shaped differential phases are reliably reproduced by combining a simple asymmetric continuum disc model with our Br$\gamma$ disc-wind model. The Br$\gamma$ emission of HD 98922 can be modelled with a disc wind that is able to approximately reproduce all interferometric observations if we assume that the intensity distribution of the dust continuum disc is asymmetric.

FUor outbursts in young stellar objects (YSOs) are the most dramatic events among episodic accretion phenomena. The origin of these bursts is not clear: disk instabilities and/or disk perturbations by an external body being the most viable hypotheses. Here, we report our VLT/SINFONI high angular resolution AO-assisted observations of 2MASS J06593158-0405277, which is undergoing a recently discovered FUor outburst. Our observations reveal the presence of an extended disc-like structure around the FUor, a very low-mass companion (2MASS J06593158-0405277B) at ~100 au in projection, and, possibly, a third closer companion at ~11 au. These sources appear to be young, displaying accretion signatures. Assuming the components are physically linked, 2MASS J06593158-0405277 would then be one of the very few triple systems observed in FUors.

Accretion and ejection are tightly connected and represent the fundamental mechanisms regulating star formation. However, the exact physical processes involved are not yet fully understood. We present high angular and spectral resolution observations of the Br Gamma emitting region in the Herbig Ae star HD163296 (MWC275) in order to probe the origin of this line and constrain the physical processes taking place at sub-AU scales in the circumstellar region. By means of VLTI-AMBER observations at high spectral resolution (R~12000), we studied interferometric visibilities, wavelength-differential phases, and closure phases across the Br Gamma line of HD163296. To constrain the physical origin of the Br Gamma line in Herbig Ae stars, all the interferometric observables were compared with the predictions of a line radiative transfer disc wind model. The measured visibilities clearly increase within the Br Gamma line, indicating that the Br Gamma emitting region is more compact than the continuum. By fitting a geometric Gaussian model to the continuum-corrected Br Gamma visibilities, we derived a compact radius of the Br Gamma emitting region of ~0.07+/-0.02AU (Gaussian half width at half maximum; or a ring-fit radius of ~0.08+/-0.02AU). To interpret the observations, we developed a magneto-centrifugally driven disc wind model. Our best disc wind model is able to reproduce, within the errors, all the interferometric observables and it predicts a launching region with an outer radius of ~0.04AU. However, the intensity distribution of the entire disc wind emitting region extends up to ~0.16AU. Our observations, along with a detailed modelling of the Br Gamma emitting region, suggest that most of the Br Gamma emission in HD163296 originates from a disc wind with a launching region that is over five times more compact than previous estimates of the continuum dust rim radius.

We aim at deriving the main physical properties of massive jets from near-IR observations, comparing them to those of a large sample of jets from low-mass YSOs, and relating them to the main features of their driving sources. We present a NIR imaging (H2 and Ks) and low-resolution spectroscopic (0.95-2.50 um) survey of 18 massive jets towards GLIMPSE extended green objects, driven by intermediate- and high-mass YSOs, which have Lbol between 4x10^2 and 10^5 Lsun. As in low-mass jets, H2 is the primary NIR coolant, detected in all the analysed flows, whereas the most important ionic tracer is [FeII], detected in half of the sampled jets. Our analysis indicates that the emission lines originate from shocks at high temperatures and densities. No fluorescent emission is detected along the flows, regardless of the source Lbol. On average, the physical parameters of these massive jets (i.e. Av, temperature, column density, mass, and luminosity) have higher values than those measured in their low-mass counterparts. The morphology of the H2 flows is varied, mostly depending on the complex, dynamic, and inhomogeneous environment in which these massive jets form and propagate. All flows and jets in our sample are collimated, showing large precession angles. Additionally, the presence of both knots and jets suggests that the ejection process is continuous with burst episodes, as in low-mass YSOs. We compare the flow H2 luminosity with the source Lbol confirming the tight correlation between these two quantities. Five sources, however, display a lower L(H2)/Lbol efficiency, which might be related to YSO evolution. Most important, the inferred L(H2) vs Lbol relationship agrees well with the correlation between the momentum flux of the CO outflows and the bolometric luminosities of high-mass YSOs indicating that outflows from high-mass YSOs are momentum driven, as are their low-mass counterparts.

As part of the POISSON project (Protostellar Optical-Infrared Spectral Survey on NTT), we present the results of the analysis of low-resolution NIR spectra 0.9-2.4 um) of two samples of YSOs in Lupus and Serpens (52 and 17 objects), with masses 0.1-2.0 Msun and ages from 10^5 to a few 10^7 yr. After determining the accretion parameters of the Lup and Ser targets by analysing their HI near-IR emission features, we added the results to those from previous regions (investigated in POISSON with the same methodology). We obtained a final catalogue (143 objects) of mass accretion rates (Macc) derived in a homogeneous fashion and analysed how Macc correlates with M* and how it evolves in time. We derived the accretion luminosity (Lacc) and Macc for Lup and Ser objects from the Br_gamma line by using relevant empirical relationships from the literature that connect HI line luminosity and Lacc. To minimise the biases and also for self-consistency, we re-derived mass and age for each source using the same set of evolutionary tracks. We observe a correlation MaccM*^2.2, similarly to what has previously been observed in several star-forming clouds. The time variation of Macc is roughly consistent with the expected evolution in viscous disks, with an asymptotic decay that behaves as t^-1.6. However, Macc values are characterised by a large scatter at similar ages and are on average higher than the predictions of viscous models. Although part of the scattering may be related to the employed empirical relationship and to uncertainties on the single measurements, the general distribution and decay trend of the Macc points are real. These findings might be indicative of a large variation in the initial mass of the disks, of fairly different viscous laws among disks, of varying accretion regimes, and of other mechanisms that add to the dissipation of the disks, such as photo-evaporation.

(Abridged) METHODS: We performed SMA observations at 1.3 mm with both the most extended and compact array configurations, providing sub-arcsecond and high sensitivity maps of various molecular lines, including both hot-core and outflow tracers. We also reconstruct the spectral energy distribution of the region from millimeter to near infrared wavelengths, using the Herschel/Hi-GAL maps, as well as archival data. RESULTS: From the spectral energy distribution, we derive a bolometric luminosity of about 4x10^4 Lsun. Our interferometric observations reveal that the distribution of dense gas and dust in the HMC is significantly flattened and extends up to a radius of 8000 AU from the center of radio continuum and maser emission in the region. The equatorial plane of this HMC is strictly perpendicular to the elongation of the collimated bipolar outflow, as imaged on scales of about 0.1-0.5 pc in the main CO isotopomers as well as in the SiO(5-4) line. In the innermost HMC regions (ca. 1000 AU), the velocity field traced by the CH3CN(12_K-11_K) line emission shows that molecular gas is both expanding along the outflow direction following a Hubble-law, and rotating about the outflow axis, in agreement with the (3-D) velocity field traced by methanol masers. The velocity field associated with rotation indicates a dynamical mass of 19 Msun at the center of the core. The latter is likely to be concentrated in a single O9.5 ZAMS star, consistent with the estimated bolometric luminosity of G023.01-00.41. The physical properties of the CO(2-1) outflow emission, such as its momentum rate 6x10^-3 Msun km/s /yr and its outflow rate 2x10^-4 Msun/yr, support our estimates of the luminosity (and mass) of the embedded young stellar object.

Context. HD 85567 is an enigmatic object exhibiting the B[e] phenomenon, i.e. an infrared excess and forbidden emission lines in the optical. The object's evolutionary status is uncertain and there are conflicting claims that it is either a young stellar object or an evolved, interacting binary. Aims. To elucidate the reason for the B[e] behaviour of HD 85567, we have observed it with the VLTI and AMBER. Methods. Our observations were conducted in the K-band with moderate spectral resolution (R~1500, i.e. 200 km/s). The spectrum of HD 85567 exhibits Br gamma and CO overtone bandhead emission. The interferometric data obtained consist of spectrally dispersed visibilities, closure phases and differential phases across these spectral features and the K-band continuum. Results. The closure phase observations do not reveal evidence of asymmetry. The apparent size of HD 85567 in the K-band was determined by fitting the visibilities with a ring model. The best fitting radius, 0.8 +/- 0.3 AU, is relatively small making HD 85567 undersized in comparison to the size-luminosity relationship based on YSOs of low and intermediate luminosity. This has previously been found to be the case for luminous YSOs, and it has been proposed that this is due to the presence of an optically thick gaseous disc. We demonstrate that the differential phase observations over the CO bandhead emission are indeed consistent with the presence of a compact (~1 AU) gaseous disc interior to the dust sublimation radius. Conclusions. The observations reveal no sign of binarity. However, the data do indicate the presence of a gaseous disc interior to the dust sublimation radius. We conclude that the data are consistent with the hypothesis that HD 85567 is a YSO with an optically thick gaseous disc within a larger dust disc that is being photo-evaporated from the outer edge.

Molecular outflows from very low-mass stars (VLMSs) and brown dwarfs (BDs) have been studied very little, and only a few objects have been directly imaged. Using VLT SINFONI K-band observations, we spatially resolved, for the first time, the H2 emission around IRS54, a ~0.1-0.2 Msun Class I source. The molecular emission shows a complex structure delineating a large outflow cavity and an asymmetric molecular jet. In addition, new [FeII] VLT ISAAC observations at 1.644um allowed us to discover the atomic jet counterpart which extends down to the central source. The outflow structure is similar to those found in low-mass Class I young stellar objects (YSOs) and Classical TTauri stars (CTTSs). However, its Lacc/Lbol ratio is very high (~80%), and the derived mass accretion rate is about one order of magnitude higher than in objects with similar mass, pointing to the young nature of the investigated source.

We present a detailed spectroscopic investigation of the young eruptive star PV Cep, to improve our understanding of its nature and characterise its circumstellar environment after its last outburst in 2004. The analysis of our medium-resolution spectroscopy in the near-IR (0.9-2.35 um), collected in 2012 at the Large Binocular Telescope with the IR spectrograph LUCIFER, allows us to infer the main stellar parameters (visual extinction, accretion luminosity, mass accretion and ejection rates), and model the inner disc, jet, and wind. The NIR spectrum displays several strong emission lines associated with accretion/ejection activity and circumstellar environment. Our analysis shows that the brightness of PV Cep is fading, as well as the mass accretion rate (2x10^-7 Msun/yr^-1 in 2012 vs ~5x10^-6 Msun/yr^-1 in 2004), which is more than one order of magnitude lower than in the outburst phase. Among the several emission lines, only the [FeII] intensity increased after the outburst. The observed [FeII] emission delineates blue- and red-shifted lobes, both with high- and low-velocity components, which trace an asymmetric jet and wind, respectively. The observed emission in the jet has a dynamical age of ~8 years, indicating that it was produced during the last outburst. The mass ejection rate in both lobes is 1.5x10^-7 Msun/yr^-1, approximately matching the high accretion rate observed during and immediately after the outburst . The observed jet/outflow asymmetries are consistent with an inhomogeneous medium. Our modelling of the CO emission hints at a small-scale gaseous disc ring, extending from ~0.2-0.4 AU to ~3 AU from the source, with an inner temperature of ~3000 K. Our HI lines modelling indicates that most of the observed emission comes from an expanding disc wind at Te=10000 K. The line profiles are strongly affected by scattering, disc screening, and outflow self-absorption.

Molecular outflows from very low-mass stars (VLMSs) and brown dwarfs have been studied very little. So far, only a few CO outflows have been observed, allowing us to map the immediate circumstellar environment. We present the first spatially resolved H2 emission around IRS54 (YLW52), a ~0.1-0.2 Msun Class I source. By means of VLT SINFONI K-band observations, we probed the H2 emission down to the first ~50 AU from the source. The molecular emission shows a complex structure delineating a large outflow cavity and an asymmetric molecular jet. Thanks to the detection of several H2 transitions, we are able to estimate average values along the jet-like structure (from source position to knot D) of Av~28 mag, T~2000-3000 K, and H2 column density N(H2)~1.7x10^17 cm^-2. This allows us to estimate a mass loss rate of ~2x10^-10 Msun/yr for the warm H2 component . In addition, from the total flux of the Br Gamma line, we infer an accretion luminosity and mass accretion rate of 0.64 Lsun and ~3x10^-7 Msun/yr, respectively. The outflow structure is similar to those found in low-mass Class I and CTTS. However, the Lacc/Lbol ratio is very high (~80%), and the mass accretion rate is about one order of magnitude higher when compared to objects of roughly the same mass, pointing to the young nature of the investigated source.

We present results of AO-assisted K-band IFU spectroscopy of the massive young star IRAS13481-6124 performed with ESO's VLT/SINFONI instrument. Our spectro-astrometric analysis of the Br{\gamma} line revealed a photo-center shift with respect to the adjacent continuum of \sim1 AU at a distance of 3.1 kpc. The position angle of this shift matches with that of the outflow which confirms that the massive star is indeed the driving source. Furthermore, a velocity gradient along the major disk axis was found which hints at the rotational sense of the ionized region, and thus of the disk as well. The gradient is not consistent with Keplerian motion but points to rigid rotation of the innermost disk. Notably, emission of H_2 is absent from source while both shocked and fluorescent H_2 emission are observed in its immediate surroundings.

Characterising stellar and circumstellar properties of embedded young stellar objects (YSOs) is mandatory for understanding the early stages of the stellar evolution. This task requires the combination of both spectroscopy and photometry, covering the widest possible wavelength range, to disentangle the various protostellar components and activities. As part of the POISSON project, we present a multi-wavelength spectroscopic and photometric investigation of embedded YSOs in L1641, aimed to derive the stellar parameters and evolutionary stages and to infer their accretion properties. Our database includes low-resolution optical-IR spectra from the NTT and Spitzer (0.6-40 um) and photometric data covering a spectral range from 0.4 to 1100 um, which allow us to construct the YSOs spectral energy distributions (SEDs) and to infer the main stellar parameters. The SED analysis allows us to group our 27 YSOs into nine Class I, eleven Flat, and seven Class II objects. However, on the basis of the derived stellar properties, only six Class I YSOs have an age of ~10^5 yr, while the others are older 5x10^5-10^6 yr), and, among the Flat sources, three out of eleven are more evolved objects (5x10^6-10^7 yr), indicating that geometrical effects can significantly modify the SED shapes. Inferred mass accretion rates (Macc) show a wide range of values (3.6x10^-9 to 1.2x10^-5 M_sun yr^-1), which reflects the age spread observed in our sample. Average values of mass accretion rates, extinction, and spectral indices decrease with the YSO class. The youngest YSOs have the highest Macc, whereas the oldest YSOs do not show any detectable jet activity in either images and spectra. We also observe a clear correlation among the YSO Macc, M*, and age, consistent with mass accretion evolution in viscous disc models.

Most of our knowledge about star formation is based on studies of low-mass stars, whereas very little is known about the properties of the circumstellar material around young and embedded intermidiate-mass T Tauri stars (IMTTSs). We present an analysis of the excitation and accretion properties of the young IMTTS DK Cha. The nearly face-on configuration of this source allows us to have direct access to the star-disk system through the excavated envelope and outflow cavity. Based on LR optical and IR spectroscopy obtained with SofI/EFOSC2 on the NTT we derive the spectrum of DK Cha from ~0.6 to 2.5 \mu m. From the detected lines we probe the conditions of the gas that emits the HI IR emission lines. In addition, we derive the mass accretion rate (Macc) from the relationships that connect the luminosity of the Br\gamma\ and Pa\beta\ lines with the accretion luminosity (Lacc). The observed optical/IR spectrum is extremely rich in forbidden and permitted atomic and molecular emission lines, which makes this source similar to very active low-mass T Tauri stars. Some of the permitted emission lines are identified as being excited by fluorescence. We derive Brackett decrements and compare them with different excitation mechanisms. The Pa\beta/Br\gamma\ ratio is consistent with optically thick emission in LTE at a temperature of ~3500 K, originated from a compact region of ~5 Rsun in size: but the line opacity decreases in the Br lines for high quantum numbers n_{up}. A good fit to the data is obtained assuming an expanding gas in LTE, with an electron density at the wind base of ~10^13 cm-3. In addition, we find that the observed Brackett ratios are very similar to those reported in previous studies of low-mass CTTSs and Class I sources, indicating that these ratios are not dependent on masses and ages. Finally, Lacc~9 Lsun and Macc~3x10-7 Msun/yr values were found.

We present the results of the analysis of LR optical-NIR spectra (0.6-2.4 um) of a sample 47 YSOs in the ChaI and II star-forming clouds. These data are part of the POISSON project (Protostellar Optical-Infrared Spectral Survey on NTT). The aim is to determine the accretion luminosity (Lacc) and mass accretion rate (Macc) of the sources through the analysis of the detected emission features. We also aim at verifying the reliability and consistency of the existing empirical relationships connecting emission line luminosity and Lacc. We employ five tracers (OI-6300A, Ha, CaII-8542A, Pab, and Brg) to derive the accretion luminosity. The tracers provide Lacc values showing different scatters when plotted as a function of L*. The Brg seems to be the most reliable, because it gives the minimum Lacc dispersion over the entire range of L*, whereas the other tracers provide much more scattered Lacc values, which are not expected for our homogeneous sample. The comparison between Lacc(Brg) and Lacc obtained from the other tracers also shows systematic differences among the empirical relationships. These may probably be ascribed to different excitation mechanisms contributing to the line emission, which may vary between our sample and those where the relationships were calibrated. Adopting the Lacc derived from Brg, we find Lacc=0.1L*-1L* for all sources, and Macc of the order of 10^-7-10^-9 Msun/yr. The Macc derived in ChaI are proportional to M*^2, as found in other low-mass star-forming regions. The discrepancies observed in the case of Lacc(Brg) and Lacc(Pab) can be related to different intrinsic Pab/Brg, ratios. The derived ratios show the existence of two different emission modalities, one that agrees with predictions of both wind and accretion models, the other suggesting optically thick emission from relatively small regions (10^21-10^22 cm^-3) with gas at low temperatures (<4000K).

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 deg^2 field including the Corona Australis (CrA) star-forming region. At a distance of 130 pc, CrA is one of the closest regions known to be actively forming stars, particularly within its embedded association, the Coronet. Using the Spitzer data, we identify 51 young stellar objects (YSOs) in CrA which include sources in the well-studied Coronet cluster as well as distributed throughout the molecular cloud. Twelve of the YSOs discussed are new candidates, one of which is located in the Coronet. Known YSOs retrieved from the literature are also added to the list, and a total of 116 candidate YSOs in CrA are compiled. Based on these YSO candidates, the star formation rate is computed to be 12 M_o Myr^-1, similar to that of the Lupus clouds. A clustering analysis was also performed, finding that the main cluster core, consisting of 68 members, is elongated (having an aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density of 150 pc^-2. In addition, we analyze outflows and jets in CrA by means of new CO and H_2 data. We present 1.3 mm interferometric continuum observations made with the Submillimeter Array (SMA) covering R CrA, IRS 5, IRS 7, and IRAS 18595-3712 (IRAS 32). We also present multi-epoch H_2 maps and detect jets and outflows, study their proper motions, and identify exciting sources. The Spitzer and ISAAC/VLT observations of IRAS 32 show a bipolar precessing jet, which drives a CO (2-1) outflow detected in the SMA observations. There is also clear evidence for a parsec-scale precessing outflow, E-W oriented, and originating in the SMA 2 region, likely driven by SMA 2 or IRS 7A.

Strong outbursts in very young and embedded protostars are rare and not yet fully understood. They are believed to originate from an increase of the mass accretion rate onto the source. We report the discovery of a strong outburst in a low-mass embedded young stellar object (YSO), namely 2MASS-J05424848-0816347 or [CTF93]216-2, as well as its photometric and spectroscopic follow-up. Using near- to mid-IR photometry and NIR low-resolution spectroscopy, we monitor the outburst, deriving its magnitude, duration, as well as the enhanced accretion luminosity and mass accretion rate. [CTF93]216-2 increased in brightness by ~4.6, 4.0, 3.8, and 1.9 mag in the J, H, Ks bands and at 24 um, respectively, corresponding to an L_bol increase of ~20 L_sun. Its early spectrum, probably taken soon after the outburst, displays a steep almost featureless continuum, with strong CO band heads and H_2O broad-band absorption features, and Br gamma line in emission. A later spectrum reveals more absorption features, allowing us to estimate T_eff~3200 K, M~0.25 M_sun, and mass accretion rate~1.2x10^{-6} M_sun yr^{-1}. This makes it one of the lowest mass YSOs with a strong outburst so far discovered.

IRAS04325+2402 (herafter IRAS04325) is a complex protostellar system hosting two young stellar objects (AB and C in the following) at a separation of 1250AU. Here we present new deep Gemini imaging and spectroscopy for the system covering the wavelength regime from 1-12mu as well as Sub-Millimeter Array interferometry at 870mu, in combination with Spitzer and literature data. Object AB is a low-mass star with a disk/envelope system and an outflow cavity, which is prominently seen in infrared images. Object C, previously suspected to be a brown dwarf, is likely a very low mass star, with an effective temperature of ~3400K. It features an edge-on disk and an elongated envelope, and shows strong indications for accretion and ejection activity. Both objects are likely to drive parsec-scale molecular outflows. The two objects are embedded in an isolated, dense molecular cloud core. High extinction, lack of X-ray emission, and relatively high bolometric luminosity argue for a very young age below 1Myr. The disk/outflow systems of AB and C are misaligned by ~60deg against each other and by 80 and 40deg against the orbital plane of the binary. The system might be a good case for primordial misalignment, as opposed to misalignment caused by dynamical interactions, because the outflow direction is constant and the realignment timescale is likely larger than the system age. This favours turbulent fragmentation, rather than rotational fragmentation, as the formation scenario. We show that the spectral energy distributions and images for the two objects can be reproduced with radiative transfer models for disk/envelope systems. Our analysis provides reassurance in the established paradigm for the structure and early evolution of YSOs, but stresses the importance of developing 3D models with sophisticated dust chemistry (abridged).

We present a multi-epoch (20 years baseline) kinematical investigation of HH52, 53, and 54 at optical and near-IR wavelengths, along with medium and high- resolution spectroscopic analyses, probing the kinematical and physical time variability conditions of the gas along the flows. By means of multi-epoch and multi-wavelength narrow-band images, we derived proper motions, tangential velocities, velocity and flux variability of the knots. Radial velocities and physical parameters of the gas were derived from spectroscopy. Finally, spatial velocities and inclination of the flows were obtained by combining both imaging and spectroscopy. The P.M. analysis reveals three distinct, partially overlapping outflows. In 20 years, about 60% of the knots show some degree of flux variability. Our set of observations apparently indicates acceleration and deceleration in a variety of knots along the jets. For about 20% of the knots, mostly coincident with working surfaces or interacting knots along the flows, a relevant variability in both flux and velocity is observed. We argue that both variabilities are related and that all or part of the kinetic energy lost by the interacting knots is successively radiated. The analysis indicates the presence of very light, ionised, and hot flows, impacting a denser medium. Several knots are deflected. At least for a couple of them (HH54 G and G0), the deflection originates from the collision of the two. For the more massive parts of the flow, the deflection is likely the result of the flow collision with a dense cloud or with clumps.

We present here an extensive analysis of the protostellar jet driven by IRAS 20126+4104, deriving the kinematical, dynamical, and physical conditions of the H2 gas along the flow. The jet has been investigated by means of near-IR H2 and [FeII] narrow-band imaging, high resolution spectroscopy of the 1-0S(1) line (2.12 um), NIR (0.9-2.5 um) low resolution spectroscopy, along with ISO-SWS and LWS spectra (from 2.4 to 200 um). The flow shows a complex morphology. In addition to the large-scale jet precession presented in previous studies, we detect a small-scale wiggling close to the source, that may indicate the presence of a multiple system. The peak radial velocities of the H2 knots range from -42 to -14 km s^-1 in the blue lobe, and from -8 to 47 km s^-1 in the red lobe. The low resolution spectra are rich in H_2 emission, and relatively faint [FeII] (NIR), [OI] and [CII] (FIR) emission is observed in the region close to the source. A warm H2 gas component has an average excitation temperature that ranges between 2000 K and 2500 K. Additionally, the ISO-SWS spectrum reveals the presence of a cold component (520 K), that strongly contributes to the radiative cooling of the flow and plays a major role in the dynamics of the flow. The estimated L(H2) of the jet is 8.2+/-0.7 L_sun, suggesting that IRAS20126+4104 has an accretion rate significantly increased compared to low-mass YSOs. This is also supported by the derived mass flux rate from the H2 lines (Mflux(H2)~7.5x10^-4 M_sun yr^-1). The comparison between the H2 and the outflow parameters strongly indicates that the jet is driving, at least partially, the outflow. As already found for low-mass protostellar jets, the measured H2 outflow luminosity is tightly related to the source bolometric luminosity.

Aims. We present a detailed study of the infrared line emission (1-200 micron) in the Herbig-Haro object HH54. Our database comprises: high- (R~9000) and low- (R~600) resolution spectroscopic data in the near-infrared band (1-2.5 micron); mid-infrared spectrophotometric images (5-12 micron); and, far-IR (45-200 micron, R~200) spectra acquired with the ISO satellite. As a result, we provide the detection of and the absolute fluxes for more than 60 molecular features (mainly from H2 in the near- and mid-infrared and from H2O, CO and OH in the far-infrared) and 23 ionic lines. Methods. The H2 lines, coming from levels from v=0 to v=4 have been interpreted in the context of a state-of-art shock code, whose output parameters are adopted as input to a Large Velocity Gradient computation in order to interpret the FIR emission of CO, H2O and OH. Results. The H2 emission can be interpreted as originating in either steady-state J-type shocks or in quasi-steady J-type shocks with magnetic precursor. However, our multi-species analysis shows that only a model of a J-type shock with magnetic precursor (v_shock=18 km/s, n_H=10^4 cm^-3, B=100 microG, age=400 yr) can account for both the observed H2 emission and the CO and H2O lines. Such a model predicts a H2O abundance of ~ 7 10^-5, in agreement with estimations from other shock models for outflows associated with low mass protostars. We can exclude the possibility that the observed atomic lines arise in the same shock as the molecular lines, and give arguments in favour of the presence of a further high-velocity, fully dissociative shock component in the region. Finally, in view of the forthcoming spectroscopic facilities on board of the Herschel satellite, we provide predictions for H2O lines considered to be the most suitable for diagnostic purposes.

We present a NIR analysis of a sample of H2 outflows from young embedded sources to compare the physical properties and cooling mechanisms of the different flows. The sample comprises 23 outflows driven by Class 0 and I sources having low-intermediate luminosity. We have obtained narrow band images in H2 2.12um and FeII 1.64um and spectroscopic observations in the range 1-2.5um. H2 line ratios have been used to estimate the visual extinction and average temperature of the molecular gas. Av values range from 2 to 15mag; average temperatures range between 2000 and 4000K. In several knots a stratification of temperatures is found with maximum values up to 5000K. Such a stratification is more commonly observed in those knots which also show FeII emission, while a thermalized gas at a single temperature is generally found in knots emitting only in molecular lines. Combining narrow band imaging with the parameters derived from the spectral analysis, we are able to measure the total luminosity of the H2 and FeII shocked regions in each flow. H2 is the major IR coolant with an average L(H2)/L(FeII) ratio of 10^{-2}. 83% of the sources have a L(H2)/L(bol) ratio of 0.04, irrespective of the Class of the driving source, while a smaller group of sources have L(H2)/L(bol) an order of magnitude smaller.Such a separation reveals the non-homogeneous behaviour of ClassI, where sources with very different outflow activity can be found. This is consistent with other studies showing that among ClassI one can find objects with different accretion properties and it demonstrates that the H2 power in the jet can be a powerful tool to identify the most active sources of this class.

We present near infrared (IR) spectra (0.98--2.5 mu m) from the group of Herbig-Haro (HH) objects comprising HH99: a series of knots forming a bow, HH99B, and a separate knot, HH99A. Observations of H2, [Fe II] and [C I] are used to constrain shock model parameters and determine the origin of the emission. Previous work has shown that it is likely that the atomic and ionic emission arises in regions of higher ionization than the molecular emission. On the basis of observations, it has often been suggested that the [Fe II] and [C I] emission could arise, for example, at the apex of a bow shock, with the H2 emission produced in the wings. Accordingly, we have combined models of C-component and J-type shocks in order to reproduce the observed H2, [Fe II] and [C I] line spectra. We can account for the H2 emission towards the HH99B complex by means of J-type shocks with magnetic precursors. We derive shock velocities of about 30 km s-1 and ages of order 100 yr; the pre-shock gas has a density of approximately 10^4 cm-3 and the magnetic field is approximately 100 mu G. The J-type shocks required to reproduce the intensities of the observed [Fe II] and [C I] lines have velocities of 50 km s-1. It is necessary to assume that Fe has been previously eroded from grains, probably by the earlier passage of a C-type shock wave. Thus, our analysis supports the view that molecular outflows are episodic phenomena whose observed emission arises in shock waves.

The complete near infrared (0.9-2.5 micron) spectra in three different star forming regions (HH24-26, HH72 and BHR71) are presented and analyzed in the framework of shock excitation models. The spectra are dominated by H_2 rovibrational emission (vibrational state nu up to 5, excitation energy T_ex up to 35000 K), while emission from ionized material, recognizable from [Fe II] and [S II] lines, is significantly fainter. The analysis of the H_2 excitation diagrams points to the existence of two different excitation regimes: whilst condensations observed only in the infrared appear to have temperatures rarely exceeding 3000 K and can be modelled in the framework of steady-state C-shock models, the infrared counterparts of Herbig Haro (HH) objects exhibit a temperature stratification with components up to more than 5000 K. The H_2 emission from representative HH objects (HH26A, HH72A and HH320A) has been successfully modelled by planar J-shocks with magnetic precursors, for which the main parameters (pre-shock density, speed) are derived. However, these same models are unable to reproduce the observed atomic and ionic emission, which probably arises from a distinct and perhaps more embedded region with respect to that traced by the H_2. Some of the physical parameters of such regions (fractional ionization, density) have been estimated in HH72, on the basis of the observed ionic lines.

On the 13th of May 2002, supernova 2002cv was discovered using a near-infrared camera working at the AZT-24 1.1m telescope at Campo Imperatore (AQ-Italy). After the infrared detection a simultaneous photometric follow-up was started at optical wavelengths. The preliminary results confirm a heavily obscured object with a V-K color not lower than 6 magnitudes, making SN 2002cv the most reddened supernova ever observed. This finding, along with the recent discovery of another obscured supernova, suggests a critical revision of the rates known to date. The estimate of the visual extinction and the light curves are provided here. These latter indicate that our SN 2002cv observations are the earliest available for a type-Ia supernova at IR wavelengths.