Very low-mass stars are known to have jets and outflows, which is indicative of a scaled-down version of low-mass star formation. However, only very few outflows in very low-mass sources are well characterized. We characterize the bipolar molecular outflow of the very low-mass star Par-Lup3-4, a 0.12 M$_{\odot}$ object known to power an optical jet. We observed Par-Lup3-4 with ALMA in Bands 6 and 7, detecting both the continuum and CO molecular gas. In particular, we studied three main emission lines: CO(2-1), CO(3-2), and $^{13}$CO(3-2). Our observations reveal for the first time the base of a bipolar molecular outflow in a very low-mass star, as well as a stream of material moving perpendicular to the primary outflow of this source. The primary outflow morphology is consistent with the previously determined jet orientation and disk inclination. The outflow mass is $9.5\times10^{-7}\mathrm{M}_{\odot}$ , with an outflow rate of $4.3\times10^{-9}\mathrm{M}_{\odot}\mathrm{yr}^{-1}$ A new fitting to the spectral energy distribution suggests that Par-Lup3-4 may be a binary system. We have characterized Par-Lup3-4 in detail, and its properties are consistent with those reported in other very low-mass sources. This source provides further evidence that very low-mass sources form as a scaled-down version of low-mass stars.

Very few examples are known of red supergiant runaways, all of them descending from the more massive O-type precursors, but none from the lower mass B-type precursors, although runaway statistics among B-type stars suggest that K-type runaways must be relatively numerous. We study HD 137071, a star that has been considered so far as a normal K-type red giant. Its parallax measured by Gaia and the derived luminosity suggest that it is actually a supergiant, whereas its derived distance to the galactic plane and its spatial velocity of 54.1 km s$^{-1}$ with respect to the local standard of rest suggest that it is also a runaway star. However, intrinsic limitations in determining the trigonometric parallaxes of cool supergiants, even in the Gaia era, require accurate spectral classifications for confirmation. We reliably classify HD 137071 as a K4II star establishing its membership to the extreme Population I, which is in agreement with the luminosity derived using the Gaia DR2 parallax measurement. Kinematical data from the Gaia DR2 catalog confirm its high spatial velocity and its runaway nature. Combining the spectral classification with astrometric information, a state-of-the-art galactic potential model, and evolutionary models for high-mass stars we trace the motion of HD 137071 back to the proximities of the galactic plane and speculate on which of the two proposed mechanisms for the production of runaway stars may be responsible for the high velocity of HD 137071. The available data favor the formation of HD 137071 in a massive binary system where the more massive companion underwent a supernova explosion about 32 Myr ago.

We investigate the silicate feature of the two Class I components of HH250-IRS, a resolved binary system with a separation of $0''53$ driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13$\mu$m range using VISIR at ESO's Very Large Telescope. The silicate features of both sources are clearly different. The NW component has a broad, smooth absorption profile lacking structure. The SE component shows the silicate feature in emission, with structure longwards of 9.5$\mu$m indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differences in the components of the HH250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted in previous observations by a drop in the $L'$-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation.

Cygnus OB2 is a rich OB association in the Galaxy which has experienced intense star formation in the last 20-25 Myr. Its stellar population shows a correlation between age and Galactic longitude. Exploring the chemical composition of its stellar content we will be able to check the degree of homogeneity of the natal molecular cloud and possible effects of self-enrichment processes. Our aim is to determine silicon and oxygen abundances for a sample of eight early-type slow rotators in Cygnus OB2 in order to check possible inhomogeneities across the whole association and whether there exists a correlation of chemical composition with Galactic longitude. We have performed a spectroscopic analysis of a sample of late O and early B stars with low rotational velocity, which have been chosen so as to cover the whole association area. We have carried out an analysis based on equivalent widths of metal lines, the wings of the H Balmer lines and FASTWIND stellar atmosphere models to determine their stellar fundamental parameters as well as the silicon and oxygen surface abundances. We derive a rather homogeneous distribution of silicon and oxygen abundances across the region, with average values of 12+log(Si/H)=7.53$\pm$0.08 dex and 12+log(O/H)=8.65$\pm$0.12 dex. We find a homogeneous chemical composition in Cygnus OB2 with no clear evidence for significant chemical self-enrichment, despite indications of strong stellar winds and possible supernovae during the history of the region. Comparison with different scenarios of chemical enrichment by stellar winds and supernovae point to star forming efficiencies not significantly above 10%. The degree of homogeneity that we find is consistent with the observed Milky Way oxygen gradient based on HII regions. We also find that the oxygen scatter within Cygnus OB2 is at least of the same order than among HII regions at similar Galactocentric distance.

It has been suggested that the well-studied giant HII regions M16 and M17 may have had a common origin, being an example of large-scale triggered star formation. While some features of the distribution of the interstellar medium in the region support this interpretation, no definitive detection of an earlier population of massive stars responsible for the triggering has been made thus far. We have carried out observations looking for red supergiants in the area covered by a giant shell seen in HI and CO centered on galactic coordinates $l \sim 14^\circ 5$, $b\sim +1^\circ$ that peaks near the same radial velocity as the bulk of the emission from both giant HII regions, which are located along the shell. Red supergiants have ages in the range expected for the parent association whose most massive members could have triggered the formation of the shell and of the giant HII regions along its rim. Out of a sample of 37 bright red stars, we identify four red supergiants that confirm the existence of massive stars in the age range between $\sim 10$ and $\sim 30$~Myr in the area. At least three of them have Gaia DR2 parallaxes consistent with them being at the same distance as M16 and M17. The evidence of past massive star formation within the area of the gaseous shell lends support to the idea that it was formed by the combined action of stellar winds and ionizing radiation of the precursors of the current red supergiants. These could be the remnants of a richer population, whose most massive members have exploded already as core-collapse supernovae. The expansion of the shell against the surrounding medium, perhaps combined with the overrun of preexisting clouds, is thus a plausible trigger of the formation of a second generation of stars currently responsible for the ionization of M16 and M17.

S106 is one of the best known bipolar HII regions, thoroughly studied and modelled at infrared, submillimeter and millimeter wavelengths, and it is one of the nearest examples of the late stages of massive star formation in which the newly formed star that ionizes it is still surrounded by vast amounts of gas and dust. However, little is known about its heavily obscured central source, S106IR. The possible binarity of the central source is investigated, which is considered to be likely given the high binarity fraction among massive stars. We have carried out visible and near-infrared photometric monitoring looking for short-term variability, with special interest in that related to the presence of a close binary companion to S106IR that may produce periodic eclipses or tidal distortion of the shape of the members of the system. A periodic variability of S106IR in the J band is found with a period of 5.0 days and an amplitude of about 0.1 mag. The light curve displays a slow rise from minimum to maximum followed by a steep decrease, and can be well reproduced by a close binary system composed of two stars with different luminosity orbiting each other in an elliptical orbit of moderate eccentricity. S106IR also shows hints of short-term variability possibly related to accretion. We also report variability of four other stars previously classified as members of the S106 cluster, all of which are strong X-ray emitters. The newly discovered close binarity of S106IR adds a new element to the modeling of the nebula and to the understanding of the dynamics of the gas around the ionizing source, which suggests that the components of the binary are accreting via a circumbinary disk. Binarity also helps to explain the apparent mismatch between the spectral type of the ionizing source inferred from the nebular spectrum and its high brightness at near-infrared wavelengths.

Herbig-Haro flows are signposts of recent major accretion and outflow episodes. We aim to determine the nature and properties of the little-known outflow source HH 250-IRS, which is embedded in the Aquila clouds. We have obtained adaptive optics-assisted L-band images with the NACO instrument on the Very Large Telescope (VLT), together with N- and Q-band imaging with VISIR also on the VLT. Using the SINFONI instrument on the VLT we carried out H- and K-band integral field spectroscopy of HH 250-IRS, complemented with spectra obtained with the SpeX instrument at the InfraRed Telescope Facility (IRTF) in the JHKL bands. Finally, the SubMillimeter Array (SMA) interferometer was used to study the circumstellar environment of HH 250-IRS at 225 and 351 GHz with CO (2-1) and CO (3-2) maps and 0.9 mm and 1.3 mm continuum images. The HH 250-IRS source is resolved into a binary with 0''53 separation, corresponding to 120 AU at the adopted distance of 225 pc. The individual components show heavily veiled spectra with weak CO absorption indicative of late-type stars. Both are Class I sources, but their spectral energy distributions between 1.5 $\mu$m and 19 $\mu$m differ markedly and suggest the existence of a large cavity around one of the components. The millimeter interferometric observations indicate that the gas mainly traces a circumbinary envelope or disk, while the dust emission is dominated by one of the circumstellar envelopes. HH 250-IRS is a new addition to the handful of multiple systems where the individual stellar components, the circumstellar disks and a circumbinary disk can be studied in detail, and a rare case among those systems in which a Herbig-Haro flow is present.

The Cygnus complex is one of the most powerful star forming regions at a close distance from the Sun (~1.4 kpc). Its richest OB association Cygnus OB2 is known to harbor many tens of O-type stars and hundreds of B-type stars, providing a large homogeneous population of OB stars that can be analyzed. Many studies of its massive population have been developed in the last decades, although the total number of OB stars is still incomplete. Our aim is to increase the sample of O and B members of Cygnus OB2 and its surroundings by spectroscopically classifying 61 candidates as possible OB-type members of Cygnus OB2. We have obtained new blue intermediate-resolution spectra suitable for spectral classification of the 61 candidates in Cygnus OB2 and surroundings. We thus performed a spectral classification of the sample using He I-II and metal lines rates, as well as the Marxist Ghost Buster (MGB) software for O-type stars and the IACOB standards catalog for B-type stars. Out of the 61 candidates, we have classified 42 stars as new massive OB-type stars, earlier than B3, in Cygnus OB2 and surroundings, including 11 O-type stars. The other candidates are discarded as they display later spectral types inconsistent with membership in the association. However, the magnitude cutoff and dust extinction introduce an incompleteness. Many O and early B stars at B > 16 mag are still undiscovered in the region. Finally, we have studied the age and extinction distribution of our sample within the region, placing them in the Hertzsprung-Russell Diagram using different stellar models in order to assess age uncertainties. Massive star formation in Cygnus OB2 seems to have proceeded from lower to higher Galactic longitudes, regardless of the details of the models used. The correlation between age and Galactic longitude previously found in the region is now confirmed.

Globules and pillars, impressively revealed by the Spitzer and Herschel satellites, for example, are pervasive features found in regions of massive star formation. We studied the globule IRAS 20319+3958 in Cygnus X by means of visible and near-infrared imaging and spectroscopy, complemented with mid-infrared Spitzer/IRAC imaging, in order to obtain a census of its stellar content and the nature of its embedded sources. Our observations show that the globule contains an embedded aggregate of about 30 very young ($\lesssim 1$~Myr) stellar objects, for which we estimate a total mass of ~ 90 M_sun. The most massive members are three systems containing early B-type stars. Two of them most likely produced very compact HII regions, one of them being still highly embedded and coinciding with a peak seen in emission lines characterising the photon dominated region (PDR). Two of these three systems are resolved binaries, and one of those contains a visible Herbig Be star. An approximate derivation of the mass function of the members of the aggregate gives hints of a slope at high masses shallower than the classical Salpeter slope, and a peak of the mass distribution at a mass higher than that at which the widely adopted log-normal initial mass function peaks. The emission distribution of H$_2$ and Br gamma, tracing the PDR and the ionised gas phase, respectively, suggests that molecular gas is distributed as a shell around the embedded aggregate, filled with centrally-condensed ionised gas. Both, the morphology and the low excitation of the HII region, indicate that the sources of ionisation are the B stars of the embedded aggregate, rather than the external UV field caused by the O stars of Cygnus OB2. The youth of the embedded cluster, combined with the isolation of the globule, suggests that star formation in the globule was triggered by the passage of the ionisation front.

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 present X-Shooter/VLT observations of a sample of 36 accreting low-mass stellar and sub-stellar objects (YSOs) in the Lupus star forming region, spanning a range in mass from ~0.03 to ~1.2Msun, but mostly with 0.1Msun < Mstar < 0.5Msun. Our aim is twofold: firstly, analyse the relationship between excess-continuum and line emission accretion diagnostics, and, secondly, to investigate the accretion properties in terms of the physical properties of the central object. The accretion luminosity (Lacc), and from it the accretion rate (Macc), is derived by modelling the excess emission, from the UV to the near-IR, as the continuum emission of a slab of hydrogen. The flux and luminosity (Ll) of a large number of emission lines of H, He, CaII, etc., observed simultaneously in the range from ~330nm to 2500nm, were computed. The luminosity of all the lines is well correlated with Lacc. We provide empirical relationships between Lacc and the luminosity of 39 emission lines, which have a lower dispersion as compared to previous relationships in the literature. Our measurements extend the Pab and Brg relationships to Lacc values about two orders of magnitude lower than those reported in previous studies. We confirm that different methodologies to measure Lacc and Macc yield significantly different results: Ha line profile modelling may underestimate Macc by 0.6 to 0.8dex with respect to Macc derived from continuum-excess measures. Such differences may explain the likely spurious bi-modal relationships between Macc and other YSOs properties reported in the literature. We derive Macc in the range 2e-12 -- 4e-8 Msun/yr and conclude that Macc is proportional to Mstar^1.8(+/-0.2), with a dispersion lower by a factor of about 2 than in previous studies. A number of properties indicate that the physical conditions of the accreting gas are similar over more than 5 orders of magnitude in Macc.

We present infrared observations of the supernova remnant G21.5-0.9 with the Very Large Telescope, the Canada-France-Hawaii Telescope and the Spitzer Space Telescope. Using the VLT/ISAAC camera equipped with a narrow-band [FeII] 1.64um filter the entire pulsar wind nebula in SNR G21.5-0.9 was imaged. This led to detection of iron line-emitting material in the shape of a broken ring-like structure following the nebula's edge. The detected emission is limb-brightened. We also detect the compact nebula surrounding PSR J1833-1034, both through imaging with the CFHT/AOB-KIR instrument (K' band) and the IRAC camera (all bands) and also through polarimetric observations performed with VLT/ISAAC (Ks band). The emission from the compact nebula is highly polarised with an average value of the linear polarisation fraction $P_{L}^{avg} \simeq 0.47$, and the swing of the electric vector across the nebula can be observed. The infrared spectrum of the compact nebula can be described as a power law of index $\alpha_{IR} = 0.7 \pm 0.3$, and suggests that the spectrum flattens between the infrared and X-ray bands.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

We report on the results of an approximately 90 ks Chandra observation of a complex region that hosts multiple sites of recent and active star formation in ARA OB1a. The field is centered on the embedded cluster RCW 108-IR and includes and a large portion of the open cluster NGC 6193. We detect over 420 X-ray sources in the field and combined these data with deep near-IR, Spitzer/IRAC and MSX mid-IR data. We find about 360 of the X-ray sources have near--IR counterparts. We divide the region into 5 parts based on the X-ray point source characteristics and extended 8 micron emission. The most clearly defined regions are the central region - identified by embedded sources with high luminosities in the both the near-IR and X-ray as well as high X-ray temperatures (about 3 keV) and the eastern region - identified by low extinction and 1 keV X-ray temperatures. Other regions, identified by their directional relationship to RCW 108-IR are less uniform - representing combinations of the first two regions, independent star formation epochs, or both. Over 18% percent of the cluster members with over 100 counts exhibit flares. Overall about 50% of the stars appear to have optically thick disks when IRAC data are employed. The largest fraction of X-ray sources are best described as possessing some disk material via a more detailed extinction fitting. We estimate that the total number of pre--main sequence stars in the field is about 1600. Approximately 800 are confined to (1.1 pc) central region.

We propose the late-O, early-B star IRS2b as the ionizing source of the Flame Nebula (NGC 2024). It has been clear that such a hot, massive star must be present in this heavily obscured region, and now it has been identified. New near-infrared photometry shows that IRS2b is the most luminous and hottest star in the young star cluster embedded in the center of NGC 2024. The near-infrared observations (5' x 5') cover ~90 % of the HII region detected in radio continuum radiation, making the probability very low that the ionizing star is not present in the field. A K-band spectrum of IRS2b obtained with ISAAC on the Very Large Telescope indicates that the spectral type of IRS2b is in the range O8V - B2V. Additional arguments indicate that its spectral type is likely closer to O8 than to B2. The corresponding amount of ionizing radiation is consistent with published radio continuum and recombination line observations.

We have studied twelve very young (1-5Myr) bona fide and candidate brown dwarfs in the ChaI star forming region in terms of their kinematic properties, the occurrence of multiple systems among them as well as their rotational characteristics. Based on high-resolution spectra taken with UVES at the VLT (8.2m), radial and rotational velocities have been measured. A kinematic study of the sample showed that their radial velocity dispersion is relatively small suggesting that they are not ejected during their formation as proposed in recent formation scenarios. By means of time-resolved UVES spectra, a radial velocity survey for close companions to the targets was conducted. The radial velocities of the targets turned out to be rather constant setting upper limits for the mass Msini of possible companions to 0.1 - 2 M_Jup. These findings hint at a rather low (<10%) multiplicity fraction of the studied brown dwarfs. Furthermore, a photometric monitoring campaign of the targets yielded the determination of rotational periods for three brown dwarf candidates in the range of 2.2 to 3.4 days. These are the first rotational periods for very young brown dwarfs and among the first for brown dwarfs at all.

How frequent are brown dwarf binaries? Do brown dwarfs have planets? Are current theoretical pre-main-sequence evolutionary tracks valid down to the substellar regime? - Any detection of a companion to a brown dwarf takes us one step forward towards answering these basic questions of star formation. We report here on a search for spectroscopic and visual companions to young brown dwarfs in the Cha I star forming cloud. Based on spectra taken with UVES at the VLT, we found significant radial velocity (RV) variations for five bona-fide and candidate brown dwarfs in Cha I. They can be caused by either a (substellar or planetary) companion or stellar activity. A companion causing the detected RV variations would have about a few Jupiter masses. We are planning further UVES observations in order to explore the nature of the detected RV variations. We also found that the RV dispersion is only ~ 2km/s indicating that there is probably no run-away brown dwarf among them. Additionally a search for companions by direct imaging with the HST and SOFI (NTT) has yielded to the detection of a few companion candidates in larger orbits.