We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

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.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

We present a catalogue of 618 high amplitude infrared variable stars (1 < DeltaK < 5 mag) detected by the two widely separated epochs of 2.2 um data in the UKIDSS Galactic plane survey, from searches covering 1470 deg^2. Most were discovered by a search of all fields at 30 < l < 230 deg. Sources include new dusty Mira variables, three new CV candidates, a blazar and a peculiar source that may be an interacting binary system. However, c. 60 per cent are YSOs, based on spatial association with star forming regions at distances ranging from 300 pc to over 10 kpc. This confirms our initial result in Contreras Pena et al.(Paper I) that YSOs dominate the high amplitude infrared variable sky in the Galactic disc. It is also supported by recently published VVV results at 295 < l < 350 deg. The spectral energy distributions of the YSOs indicate class I or flat spectrum systems in most cases, as in the VVV sample. A large number of variable YSOs are associated with the Cygnus X complex and other groups are associated with the North America/Pelican nebula, the Gemini OB1 molecular cloud, the Rosette complex, the Cone nebula, the W51 star forming region and the S86 and S236 HII regions. Most of the YSO variability is likely due to variable/episodic accretion on timescales of years, albeit usually less extreme than classical FUors and EXors. Luminosities at the 2010 WISE epoch range from c. 0.1 Lsun to 10^3 Lsun but only rarely exceed 10^2.5 Lsun.

High-amplitude variability in Young Stellar Objects (YSOs) is usually associated with episodic accretion events. It has not been observed so far in massive YSOs. Here, the high-amplitude variable star sample of ContrerasPe\~{n}a et al.(2016) has been used to search for highly-variable($\Delta$K$\ge$1\,mag) sources coinciding with dense clumps mapped using the 850\mum continuum emission by the ATLASGAL survey. 18 variable sources are centred on the sub-mm clump peaks, and coincide ($<$1") with a 24$\mu$m point or compact ($<$10") source. 13 of these 18 sources can be fit by YSO models. The 13 variable YSOs(VYSO) have luminosities of $\sim$10$^3$ L$_{\odot}$, an average mass of 8 M$_{\odot}$ and a range of ages up to 10$^6$ yr. 11 of these 13 VYSOs are located in the midst of infrared dark clouds. 9 of the 13 sources have $\Delta$K$>$2 mag, significantly higher compared to the mean variability of the entire VVV sample. The light curves of these objects sampled between 2010-2015 display rising, declining, or quasi-periodic behaviour but no clear periodicity. Light-curve analysis using Plavchan method show that the most prominent phased signals have periods of a few hundred days. The nature and time-scale of variations found in 6.7 Ghz methanol maser emission (MME) in massive stars are similar to that of the VYSO light curves. We argue that the origin of the observed variability is episodic accretion. We suggest that the timescale of a few hundred days may represent the frequency at which a spiralling disk feeds dense gas to the young massive star.

Context. ALMA observations of the high-mass star-forming region G35.20-0.74N have revealed the presence of a Keplerian disk in core B rotating about a massive object of 18 Msun, as computed from the velocity field. The luminosity of such a massive star would be comparable to (or higher than) the luminosity of the whole star-forming region. To solve this problem it has been proposed that core B could harbor a binary system. This could also explain the possible precession of the radio jet associated with this core, which has been suggested by its S-shaped morphology. Aims. To establish the origin of the free-free emission from core B and investigate the existence of a binary system at the center of this massive core and the possible precession of the radio jet. Methods. We carried out VLA continuum observations of G35.20-0.74N at 2 cm in the B configuration and at 1.3 cm and 7 mm in the A and B configurations. The bandwidth at 7 mm covers the CH3OH maser line at 44.069 GHz. Continuum images at 6 and 3.6 cm in the A configuration were obtained from the VLA archive. We also carried out VERA observations of the H2O maser line at 22.235 GHz. Results. The observations have revealed the presence of a binary system of UC/HC HII regions at the geometrical center of the radio jet in G35.20-0.74N. This binary system, which is associated with a Keplerian rotating disk, consists of two B-type stars of 11 and 6 Msun. The S-shaped morphology of the radio jet has been successfully explained as due to precession produced by the binary system. The analysis of the precession of the radio jet has allowed us to better interpret the IR emission in the region, which would be not tracing a wide-angle cavity open by a single outflow with a position angle of ~55 deg but two different flows: a precessing one in the NE-SW direction associated with the radio jet, and a second one in an almost E--W direction ...

The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. Mathematically the Graph- BLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of the GraphBLAS. Graphs represent connections between vertices with edges. Matrices can represent a wide range of graphs using adjacency matrices or incidence matrices. Adjacency matrices are often easier to analyze while incidence matrices are often better for representing data. Fortunately, the two are easily connected by matrix mul- tiplication. A key feature of matrix mathematics is that a very small number of matrix operations can be used to manipulate a very wide range of graphs. This composability of small number of operations is the foundation of the GraphBLAS. A standard such as the GraphBLAS can only be effective if it has low performance overhead. Performance measurements of prototype GraphBLAS implementations indicate that the overhead is low.

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

Given the rarity of young O star candidates, compact HII regions embedded in dense molecular cores continue to serve as potential sites to peer into the details of high-mass star formation. To uncover the ionizing sources of the most luminous and compact HII regions embedded in the RCW106 and RCW122 giant molecular clouds, known to be relatively nearby (2-4 kpc) and isolated, thus providing an opportunity to examine spatial scales of a few hundred to a thousand AU in size. High spatial resolution (0.3"), mid-infrared spectra (R=350), including the fine structure lines [ArIII] and [NeII], were obtained for four luminous compact HII regions, embedded inside the dense cores within the RCW106 and RCW122 molecular cloud complexes. At this resolution, these targets reveal point-like sources surrounded by nebulosity of different morphologies, uncovering details at spatial dimensions of <1000AU. The point-like sources display [ArIII] and [NeII] lines - the ratios of which are used to estimate the temperature of the embedded sources. The derived temperatures are indicative of mid-late O type objects for all the sources with [ArIII] emission. Previously known characteristics of these targets from the literature, including evidence of disk or accretion suggest that the identified sources may grow more to become early-type O stars by the end of the star formation process.

We present initial results of the first panoramic search for high-amplitude near-infrared variability in the Galactic Plane. We analyse the widely separated two-epoch K-band photometry in the 5th and 7th data releases of the UKIDSS Galactic Plane Survey. We find 45 stars with \Delta K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is not yet included in the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (12 stars). Sources in SFRs show spectral energy distributions (SEDs) that support classification as Young Stellar Objects (YSOs). This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow up of the DR5 sample shows at least four stars with clear characteristics of eruptive pre-main-sequence variables, two of which are deeply embedded. Our results support the recent concept of eruptive variability comprising a continuum of outburst events with different timescales and luminosities, but triggered by a similar physical mechanism involving unsteady accretion. Also, we find what appears to be one of the most variable classical Be stars.

We present a multiwavelength study of W40 star-forming region using IR observations in UKIRT JHK bands, Spitzer IRAC bands & Herschel PACS bands; 2.12 micron H2 narrow-band imaging; & radio observations from GMRT (610 & 1280 MHz), in a FoV of ~34'x40'. Spitzer observations along with NIR observations are used to identify 1162 Class II/III & 40 Class I sources in the FoV. The NN stellar surface density analysis shows that majority of these YSOs constitute the embedded cluster centered on the source IRS1A South. Some YSOs, predominantly younger population, are distributed along & trace the filamentary structures at lower stellar surface density. The cluster radius is obtained as 0.44pc - matching well with the extent of radio emission - with a peak density of 650pc^-2. The JHK data is used to map the extinction which is subsequently used to compute the cloud mass. It has resulted in 126 Msun & 71 Msun for the central cluster & the northern IRS5 region, respectively. H2 narrow-band imaging displays significant emission, which prominently resembles fluorescent emission arising at the borders of dense regions. Radio analysis shows this region as having blister morphology, with the radio peak coinciding with a protostellar source. Free-free emission SED analysis is used to obtain physical parameters of the overall region & the IRS5 sub-region. This multiwavelength scenario is suggestive of star formation having resulted from merging of multiple filaments to form a hub. Star formation seems to have taken place in two successive epochs, with the first epoch traced by the central cluster & the high-mass star(s) - followed by a second epoch which is spreading into the filaments as uncovered by the Class I sources & even younger protostellar sources along the filaments. The IRS5 HII region displays indications of swept-up material which has possibly led to the formation of protostars.

The formation of the highest mass stars are thought to be dominated by instabilities resulting from gravitation and radiation. Instabilities due to gravitation are commonly demonstrated by observations of fragmentation, but those due to effects of radiation have thus far not been found. Here I report on the NACO adaptive optics and mid-infrared diffraction-limited VISIR imaging data of an extemely luminous ultra-compact HII region G333.6-0.2. Two infrared sources, one bright in the near-infrared(appearing point-like) and another in the mid-infrared (resolved with an elliptical shape) are uncovered through this data, which are located at the heart of this region. These infrared sources appear to be embedded in the waist of a bipolar-shaped nebula and UCHII region, the lobes of which are separated by a dark patch. Dense filamentary features with finger/hook morphology are found; they appear to be connected to the two bright infrared sources and the sizes of these hook features are sharply limited to <5000 AU. The observed properties of this target and a large amount of previous data obtained from the literature are compared together with the results of various numerical simulations of high-mass star formation. This comparison favours the interpretation that the finger/hook-like structures likely represent radiatively driven Rayleigh-Taylor instabilities arising in the outflow cavity of a forming high-mass binary star system.

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

Context: The ESO Public Survey "VISTA Variables in the V\'ia L\'actea" (VVV) provides deep multi-epoch infrared observations for an unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. In this survey nearly 150 new open clusters and cluster candidates have been discovered. Aims: This is the second in a series of papers about young, massive open clusters observed using the VVV survey. We present the first study of six recently discovered clusters. These clusters contain at least one newly discovered Wolf-Rayet (WR) star. Methods: Following the methodology presented in the first paper of the series, wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results: We affirm that the six studied stellar groups are real young (2-7 Myr) and massive (between 0.8 and 2.2 10^3 Msol) clusters. They are highly obscured (Av ~ 5-24 mag) and compact (1-2 pc). In addition to WR stars, two of the six clusters also contain at least one red supergiant star. We claim the discovery of 8 new WR stars, and 3 stars showing WR-like emission lines which could be classified WR or OIf. Preliminary analysis provides initial masses of ~30-50 Msol for the WR stars. Finally,we discuss the spiral structure of the Galaxy using as tracers the six new clusters together with the previously studied VVV clusters.

Context: The ESO Public Survey "VISTA Variables in the V\'ia L\'actea" (VVV) provides deep multi-epoch infrared observations for an unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. Aims: The VVV observations will foster the construction of a sample of Galactic star clusters with reliable and homogeneously derived physical parameters (e.g., age, distance, and mass, etc.). In this first paper in a series, the methodology employed to establish cluster parameters for the envisioned database are elaborated upon by analyzing a subsample of 4 known young open clusters: Danks 1, Danks 2, RCW 79, and DBS 132. The analysis offers a first glimpse of the information that can be gleaned for the final cluster database from the VVV observations. Methods: Wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results: Results inferred from the deep near-infrared photometry which features mitigated uncertainties (e.g. the accuracy of the photometry is better than 0.1mag for Ks<18mag), the wide field-of-view of the VVV survey, and numerous high quality low resolution spectra (typically more than 10 per cluster), are used to establish independent cluster parameters which enable existing determinations to be assessed. An anomalous reddening law in the direction toward the Danks' clusters is found, i.e. E(J-H)/E(H-Ks)=2.20+/-0.06, which exceeds published values for the inner Galaxy. The G305 star forming complex, which includes the Danks' clusters, lies beyond the Sagittarius-Carina spiral arm and occupies the Centaurus arm. Finally, the first deep infrared color-magnitude diagram of RCW 79 is presented which reveal a sizable pre-main sequence population. A list of candidate variable stars in G305 region is reported.

We aimed to map the jets and outflows from the Serpens South star forming region and find an empirical relationship between the magnetic field and outflow orientation. Near-infrared H2 v=1-0 S(1) 2.122{\mu}m -line imaging of the \sim 30'-long filamentary shaped Serpens South star forming region was carried out. K s broadband imaging of the same region was used for continuum subraction. Candidate driving sources of the mapped jets/outflows are identified from the list of known protostars and young stars in this region, which was derived from studies using recent Spitzer and Herschel telescope observations. 14 Molecular Hydrogen emission-line objects(MHOs) are identified using our continuum-subtracted images. They are found to constitute ten individual flows. Out of these, nine flows are located in the lower-half(southern) part of the Serpens South filament, and one flow is located at the northern tip of the filament. Four flows are driven by well-identified Class 0 protostars, while the remaining six flows are driven by candidate protostars mostly in the Class I stage, based on the Spitzer and Herschel observations. The orientation of the outflows is systematically perpendicular to the direction of the near-infrared polarization vector, recently published in the literature. No significant correlation was observed between the orientation of the flows and the axis of the filamentary cloud.

We uncover the H2 flows in the Corona Australis molecular cloud and in particular identify the flows from the Coronet cluster. Near-infrared H2 v=1--0 S(1), 2.12micron-line, narrow-band imaging survey of the R CrA cloud core was carried out. We identify the best candidate-driving source for each outflow by comparing the flow properties, available proper motions, and the known/estimated properties of the driving sources. We also adopted the thumbrule of outflow power as proportional to source luminosity and inversely proportional to the source age to reach a consensus. Results: Continuum-subtracted, narrow-band images reveal several new Molecular Hydrogen emission-line Objects (MHOs). Together with previously known MHOs and Herbig-Haro objects we catalog at least 14 individual flow components of which 11 appear to be driven by the RCrA aggregate members. The flows originating in the Coronet cluster have lengths of ~0.1-0.2 pc. Eight out of nine submillimeter cores mapped in the Coronet cluster region display embedded stars driving an outflow component. Roughly 80% of the youngest objects in the Coronet are associated with outflows. The MHO flows to the west of the Coronet display lobes moving to the west and vice-versa, resulting in nondetections of the counter lobe in our deep imaging. We speculate that these counterflows may be experiencing a stunting effect in penetrating the dense central core. Conclusions:Although this work has reduced the ambiguities for many flows in the Coronet region, one of the brightest H2 feature (MHO2014) and a few fainter features in the region remain unassociated with a clear driving source. The flows from Coronet, therefore, continue to be interesting targets for future studies.

In some protostellar objects both wide angle outflows and collimated jets are seen, while in others only one is observed. Spitzer provides unprecedented sensitivity in the infrared to study both the jet and outflow features. Here, we use HiRes deconvolution to improve the visualization of spatial morphology by enhancing resolution (to sub-arcsecond levels in the IRAC bands) and removing the contaminating sidelobes from bright sources. We apply this approach to study the jet and outflow features in Cep E a young, energetic Class 0 protostar. In the reprocessed images we detect: (i) wide angle outflow seen in scattered light; (ii) morphological details on at least 29 jet driven bow shocks and jet heads or knots; (iii) three compact features in 24 micron continuum image as atomic/ionic line emission coincident with the jet heads; and, (iv) a flattened 35 arcsec size protostellar envelope seen against the interstellar background PAH emission as an absorption band across the protostar at 8 micron. By separating the protostellar photospheric scattered emission in the wide angle cavity from the jet emission we show that we can study directly the scattered light spectrum. We present the H2 emission line spectra, as observed in all IRAC bands, for 29 knots in the jets and bowshocks and use them in the IRAC color -- color space as a diagnostic of the thermal gas in the shocks driven by the jets. The data presented here will enable detailed modeling of the individual shocks retracing the history of the episodic jet activity and the associated accretion on to the protostar. The Spitzer data analysis presented here shows the richness of its archive as a resource to study the jet/outflow features in H2 and scattered light in a large homogeneous sample.

VISTA Variables in the V\'{\i}a L\'actea (VVV) is one of the six ESO Public Surveys operating on the new 4-meter Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV is scanning the Milky Way bulge and an adjacent section of the disk, where star formation activity is high. One of the principal goals of the VVV Survey is to find new star clusters of different ages. In order to trace the early epochs of star cluster formation we concentrated our search in the directions to those of known star formation regions, masers, radio, and infrared sources. The disk area covered by VVV was visually inspected using the pipeline processed and calibrated $K_{\rm S}$-band tile images for stellar overdensities. Subsequently, we examined the composite $JHK_{\rm S}$ and $ZJK_{\rm S}$ color images of each candidate. PSF photometry of $15\times15$ arcmin fields centered on the candidates was then performed on the Cambridge Astronomy Survey Unit reduced images. After statistical field-star decontamination, color-magnitude and color-color diagrams were constructed and analyzed. We report the discovery of 96 new infrared open clusters and stellar groups. Most of the new cluster candidates are faint and compact (with small angular sizes), highly reddened, and younger than 5\,Myr. For relatively well populated cluster candidates we derived their fundamental parameters such as reddening, distance, and age by fitting the solar-metallicity Padova isochrones to the color-magnitude diagrams.

BVRIJHK photometry, Spitzer-GLIMPSE photometry and HK band spectroscopy were used to study the stellar content of IRAS 19343+2026, a (proto)star/cluster candidate, located close to the Galactic plane. The data suggest that IRAS 19343+2026 is a rich cluster associated with a massive protostar of 7.6 Msol with an age of ~ 10^5 yr. Three point sources in the vicinity of the far-infrared (FIR) peak are also found to be early B type stars. The remaining (predominantly low mass) members of the cluster are best represented by a 1 - 3 Myr pre-main-sequence (PMS) population. HK band spectra of two bright and five faint point sources in the cluster confirm that the results obtained from the photometry are good representations of their young stellar object (YSO) nature. Thus, IRAS 19343+2026 is a young cluster with at least four early B-type stars classified as young (10^4 - 10^5 yr), that are surrounded by a somewhat older (1 - 3 Myr) population of low mass YSOs. Together, these results argue for a scenario in which low mass stars form prior to massive stars in a cluster forming environment. We compute the Initial Mass Function (IMF) for this cluster using the K-band luminosity function; the slope of the IMF is shallower than predicted by the Salpeter's mass function. The cluster mass, Mtot, is estimated to be in the range ~ 307 Msol (from the data completeness limit) - 585 Msol (extrapolated down to the brown dwarf limit, assuming a certain IMF).

We aim to estimate and analyse the physical properties of the infrared counterparts of HMPOs by comparing their spectral energy distributions (SED) with those predicted by radiative transfer accretion models of YSOs. The SED of 68 IRCs are extended beyond the GLIMPSE photometry to the possible limits, from the near-infrared to the millimetre wavelengths by using the 2MASS, GLIMPSE version 2.0 catalogs, MSX, IRAS and some single dish (and interferometric) (sub)mm data. An online SED fitting tool that uses 2D radiative transfer accretion models of YSOs is employed to fit the observed SED to obtain various physical parameters. The SED of IRCs were fitted by models of massive protostars with a range of masses between 5-42 Msun and ages between 10^3 and 10^6 years. The median mass and age are 10 Msun and 10^4 yrs. The envelopes are large with a mean size of ~ 0.2-0.3 pc and show a distribution that is very similar to the distribution of the sizes of 8 micron nebulae discussed in Paper I. The estimated envelope accretion rates are high with a mean value of 10^(-3) Msun/yr and show a power law dependence to mass with an exponent of 2, suggesting spherical accretion at those scales. Disks are found to exist in most of the sources with a mean mass of 10^(-1.4+-0.7) Msun. The observed infrared-millimetre SED of the infrared counterparts of HMPOs are successfully explained with an YSO accretion model. The modelled sources mostly represent proto-B stars although some of them could become O stars in future. We demonstrate that many of these results may represent a realistic picture of massive star formation, despite some of the results which may be an effect of the assumptions within the models.

We present wide-field near-IR images of Orion A. K and H2 1-0S(1) images of a contiguous 8 sqr degree region are compared to photometry from Spitzer and dust-continuum maps obtained with MAMBO and SCUBA. We also measure proper motions for H2 features in 33 outflows. We increase the number of known H2 outflows in Orion A to 116. A total of 111 H2 flows were observed with Spitzer; outflow sources are identified for at least 72 of them. The MAMBO 1200 micron maps cover 97 H2 flows; 57 of them are associated with Spitzer sources and dust cores or extended emission. The H2 jets are widely distributed and randomly orientated; the jets do not appear to be orthogonal to large-scale filaments or even to the small-scale cores. Moreover, H2 jet lengths and opening angles are not obviously correlated with indicators of outflow source age - source spectral index or (sub)millimetre core flux. We demonstrate that H2 jet sources are predominantly protostellar with flat or positive spectral indices, rather than disk-excess (or T Tauri) stars. Most protostars in molecular cores drive H2 outflows. However, not all molecular cores are associated with protostars or H2 jets. On statistical grounds, the H2 jet phase may be marginally shorter than the protostellar phase, though must be considerably shorter than the prestellar phase. In terms of their spectral index, H2 jet sources are indistinguishable from protostars. The few true protostars without H2 jets are almost certainly more evolved than their H2-jet-driving counterparts. We also find that protostars that power molecular outflows are no more (nor no less) clustered than protostars that do not. The H2 emission regions in outflows from young stars clearly weaken and fade very quickly, before the source evolves from protostar to pre-main-sequence star.

We present the results of a morphological study performed to a sample of Ultracompact (UC) HII regions with Extended Emission (EE) using Spitzer--IRAC imagery and 3.6 cm VLA conf. D radio-continuum (RC) maps. Some examples of the comparison between maps and images are presented. Usually there is an IR point source counterpart to the peak(s) of RC emission, at the position of the UC source. We find that the predominant EE morphology is the cometary, and in most cases is coincident with IR emission at 8.0 $\mu$m. Preliminary results of Spitzer--IRAC photometry of a sub-sample of 13 UC HII regions with EE based on GLIMPSE legacy data are also presented. Besides, individual IRAC photometry was performed to 19 UC sources within these 13 regions. We show that UC sources lie on specific locus, both in IRAC color-color and AM-product diagnostic diagrams. Counts of young stellar sources are presented for each region, and we conclude that a proportion of ~ 2%, ~10%, and ~88% of sources in the UC HII regions with EE are, in average, Class I, II, and III, respectively.

IRAS 21413+5442 and IRAS 21407+5441 are two massive star forming regions of high luminosity, likely associated with each other. Near-infrared photometry on these two IRAS sources was performed at UKIRT using the UFTI under excellent seeing conditions yielding an angular resolution of $\sim$ 0.5 arcsec. Our results reveal details of stellar content to a completeness limit (90%) of J = 18.5, H = 18.0, and K = 17.5 mag in the two regions. In IRAS 21413+5442, we identify a late O type star, having large (H-K) color, to be near the centre of the CO jets observed by earlier authors. The UKIRT images reveal in IRAS 21407+5441, a faint but clear compact HII region around a central high - intermediate mass star cluster. We have detected a number of sources with large (H-K) color which are not detected in J band. We also present the GMRT radio continuum map at 1.28 GHz covering the entire region surrounding the two star forming clouds. The radio continuum fluxes are used to estimate the properties of HII regions which seem to support our near-IR photometric results. Based on our radio continuum map and the archival MSX 8.2 $\mu$m image, we show that the two IRAS sources likely belong to the same parent molecular cloud and conjecture that a high mass star of large IR colors, present in between the two sources, might have triggered star formation in this region. However one can not rule out the alternative possibility that Star A could be a nearby foreground star.

The young stellar population data of the Perseus, Ophiuchus and Serpens molecular clouds are obtained from the Spitzer c2d legacy survey in order to investigate the spatial structure of embedded clusters using the nearest neighbour and minimum spanning tree method. We identify the embedded clusters in these clouds as density enhancements and analyse the clustering parameter Q with respect to source luminosity and evolutionary stage. This analysis shows that the older Class 2/3 objects are more centrally condensed than the younger Class 0/1 protostars, indicating that clusters evolve from an initial hierarchical configuration to a centrally condensed one. Only IC348 and the Serpens core, the older clusters in the sample, shows signs of mass segregation (indicated by the dependence of Q on the source magnitude), pointing to a significant effect of dynamical interactions after a few Myr. The structure of a cluster may also be linked to the turbulent energy in the natal cloud as the most centrally condensed cluster is found in the cloud with the lowest Mach number and vice versa. In general these results agree well with theoretical scenarios of star cluster formation by gravoturbulent fragmentation.

The Spitzer-GLIMPSE point source catalog and images have been used to study a sample of 381 massive protostellar candidates. IRAC-Point source photometry was used to analyse colours, magnitudes and spectral indicies of the infrared counterparts (IRCs) of massive protostellar objects and a bonafide sample of 50 point sources was obtained. Spectral energy distributions (SEDs) of these 50 sources was extended to the near-infrared and millimeter range by using 2MASS and millimeter data from the literature. An online SED fitter tool based on Monte-Carlo radiative transfer of an accretion model involving star,disk and envelope was used to fit the SEDs of the 50 sources. The IRCs to massive protostellar objects are found to successfully imitate the SEDs of evolutionary phases similar to low mass star formation. Envelope accretion, rather than disk accretion is found to be dominant in building the most massive stars. Unresolved centimeter continuum emission is associated with 27 IRCs classified as massive protostars suggesting that ionised accretion flows may play an important role along with the molecular component. The morphology of the infrared nebulae surrounding the IRCs have an unusual resemblance to the morphologies of ultra-compact HII regions suggesting that these infrared nebulae are possible precursors to the UCHII regions.

Here we present a brief overview of UC HII regions with EE (UC HII+EE) including our most recent effort aimed at searching for UC HII regions associated with extended emission.

Aims: Some theoretical scenarios suggest the formation of brown dwarfs as ejected stellar embryos in star-forming clusters. Such a formation mechanism can result in different spatial distributions of stars and substellar objects. We aim to investigate the spatial structure of stellar and substellar objects in two well sampled and nearby embedded clusters, namely IC348 and the Orion Trapezium Cluster (OTC) to test this hypothesis. Methods:Deep near-infrared K-band data complete enough to sample the substellar population in IC348 and OTC are obtained from the literature. The spatial distribution of the K-band point sources is analysed using the Minimum Spanning Tree (MST) method. The Q parameter and the spanning trees are evaluated for stellar and substellar objects as a function of cluster core radius R$_c$. Results: The stellar population in both IC348 and OTC display a clustered distribution whereas the substellar population is distributed homogeneously in space within twice the cluster core radius. Although the substellar objects do not appear to be bound by the cluster potential well, they are still within the limits of the cluster and not significantly displaced from their birth sites. Conclusions: The spatially homogeneous distribution of substellar objects is best explained by assuming higher initial velocities, distributed in a random manner and going through multiple interactions. The overall spatial coincidence of these objects with the cluster locations can be understood if these objects are nevertheless travelling slowly enough so as to feel the gravitational effect of the cluster. The observations support the formation of substellar objects as ``ejected stellar embryos''. Higher ejection velocities are necessary but net spatial displacements may not be necessary to explain the observational data.

The GLIMPSE archive was used to obtain 3.6--8.0micron, point source photometry and images for 381 massive protostellar candidates lying in the Galactic mid-plane. The colours, magnitudes and spectral indicies of sources in each of the 381 target fields were analysed and compared with the predictions of 2D radiative transfer model simulations. Although no discernable embedded clusters were found in any targets, multiple sources or associations of redenned young stellar objects were found in many sources indicating multiplicity at birth. The spectral index ($\alpha$) of these point sources in 3.6--8.0mum bands display large values of $\alpha$=2--5. A color-magnitude analog plot was used to identify 79 infrared counterparts to the HMPOs. Compact nebulae are found in 75% of the detected sources with morphologies that can be well described by core-halo, cometary, shell-like and bipolar geometries similar to those observed in ultra-compact HII regions. The IRAC band SEDs of the IR counterparts of HMPOs are best described to represent YSOs with a mass range of 8--20\msun in their Class I stages when compared with 2D radiative transfer models. They also suggest that the high $\alpha$ values represent reprocessed star/star+disk emission that is arising in the dense envelopes. Thus we are witnessing the luminous envelopes around the protostars rather than their photospheres or disks. We argue that the compact infrared nebulae likely reflect the underlying physical structure of the dense cores and are found to imitate the morphologies of known UCHII regions. Our results favour models of continuuing accretion involving both molecular and ionised accretion components to build the most massive stars rather than purely molecular rapid accretion flows.

Aims: We test the technique of spectro-astrometry as a potential method to investigate the close environment of massive young stars. Method: Archival VLT near infrared K band spectra (R=8900) of three massive young stellar objects and one Wolf-Rayet star are examined for spectro-astrometric signatures. The young stellar objects display emission lines such as Brackett gamma, CO 2-0 and CO 3-1 that are characteristic of ionised regions and molecular disks respectively. Two of the sample sources also display emission lines such as NIII and MgII that are characteristic of high temperatures. Results: Most of the emission lines show spectro-astrometric signal at various levels resulting in different positional displacements. The shapes and magnitudes of the positional displacements imply the presence of large disk/envelopes in emission and expanding shells of ionised gas. The results obtained for the source 18006-2422nr766 in particular provide larger estimates (> 300AU) on CO emitting regions indicating that in MYSOs CO may arise from inner regions of extended dense envelopes as well. Conclusions: The overall results from this study demonstrate the utility of spectro-astrometry as a potential method to constrain the sizes of various physical entities such as disks/envelopes, UCHII regions and/or ionised shells in the close environment of a massive young star.

We present wide-field near-infrared images of the DR21/W75 high-mass star forming region, obtained with the Wide Field Camera, WFCAM, on the United Kingdom Infrared Telescope. Broad-band JHK and narrow-band H2 1-0S(1) images are compared to archival mid-IR images from the Spitzer Space Telescope, and 850 micron dust-continuum maps obtained with the Submillimeter Common User Bolometer Array (SCUBA). Together these data give a complete picture of dynamic star formation across this extensive region, which includes at least four separate star forming sites in various stages of evolution. The H2 data reveal knots and bow shocks associated with more than 50 individual flows. Most are well collimated, and at least five qualify as parsec-scale flows. Most appear to be driven by embedded, low-mass protostars. The orientations of the outflows, particularly from the few higher-mass sources in the region (DR21, DR21(OH), W75N and ERO~1), show some degree of order, being preferentially orientated roughly orthogonal to the chain of dusty cores that runs north-south through DR21. Clustering may inhibit disk accretion and therefore the production of outflows; we certainly do not see enhanced outflow activity from clusters of protostars. Finally, although the low-mass protostellar outflows are abundant and widely distributed, the current generation does not provide sufficient momentum and kinetic energy to account for the observed turbulent motions in the DR21/W75 giant molecular clouds. Rather, multiple epochs of outflow activity are required over the million-year timescale for turbulent decay.

Wide field near-infrared observations and Spitzer Space Telescope IRAC observations of the DR21/W75 star formation regions are presented. The photometric data are used to analyse the extinction, stellar content and clustering in the entire region by using standard methods. A young stellar population is identified all over the observed field, which is found to be distributed in embedded clusters that are surrounded by a distributed halo population extending over a larger projected area. The Spitzer/IRAC data are used to compute a spectral index value, "alpha", for each YSO in the field. We use these data to separate pure photospheres from disk excess sources. We find a small fraction of sources with "alpha" in excess of 2 to 3 (plus a handful with "alpha"~4), which is much higher than the values found in the low mass star forming region IC348 ("alpha" < 2). The sources with high values of "alpha" spatially coincide with the densest regions of the filaments and also with the sites of massive star formation. Star formation is found to be occuring in long filaments stretching to few parsecs that are fragmented over a scale of ~ 1 pc. The spatial distribution of young stars are found to be correlated with the filamentary nebulae that are prominently revealed by 8micron and 850micron observations. Five filaments are identified that appear to converge on a center that includes the DR21/DR21(OH) regions. The morphological pattern of filaments and clustering compare well with numerical simulations of star cluster formation by Bate et al. 2003.

Aims. We present molecular line observations of the southwestern part of the IC 348 young cluster, and we use them together with NIR and mm continuum data to determine the distribution of dense gas, search for molecular outflows, and analyze the ongoing star formation activity in the region. Methods. Our molecular line data consists of C18O(1--0) and N2H+(1--0) maps obtained with the FCRAO telescope at a resolution of about 50'' and CO(2--1) data obtained with the IRAM 30m telescope at a resolution of 11''. Results. The dense gas southwest of IC 348 is concentrated in two groups of dense cores, each of them with a few solar masses of material and indications of CO depletion at high density. One of the core groups is actively forming stars, while the other seems starless. There is evidence for at least three bipolar molecular outflows in the region, two of them powered by previously identified Class 0 sources, while the other one is powered by a still not well characterized low-luminosity object. The ongoing star formation activity is producing a small stellar subgroup in the cluster. Using the observed core characteristics and the star formation rate in the cluster we propose that that similar episodes of stellar birth may have produced the subclustering observed in the halo of IC 348.

We report on the identification of 54 embedded clusters around 217 massive protostellar candidates of which 34 clusters are new detections. The embedded clusters are identified as stellar surface density enhancements in the 2 $\mu$m All Sky Survey (2MASS) data. Because the clusters are all associated with massive stars in their earliest evolutionary stage, the clusters should also be in an early stage of evolution. Thus the properties of these clusters should reflect properties associated with their formation rather than their evolution. For each cluster, we estimate the mass, the morphological type, the photometry and extinction. The clusters in our study, by their association with massive protostars and massive outflows, reinstate the notion that massive stars begin to form after the first generation of low mass stars have completed their accretion phase. Further, the observed high gas densities and accretion rates at the centers of these clusters is consistent with the hypothesis that high mass stars form by continuing accretion onto low mass stars.

We present sub-arcsecond (0.35"-0.9"), near-infrared J,H,K band photometric observations of six fields along the W51 Giant Molecular Cloud (W51 GMC). Our observations reveal four new, embedded clusters and provide a new high-resolution (0.35") view of the W51IRS2 (G49.5-0.4) region. The cluster associated with G48.9-0.3 is found to be a double cluster enclosed in a nest of near-infrared nebulosity. We construct stellar surface density maps for four major clusters in the W51 GMC. These unveil the underlying hierarchical structure. Color-color and color-magnitude diagrams for each of these clusters show clear differences in the embedded stellar populations and indicate the relative ages of these clusters. In particular, the clusters associated with the HII regions G48.9-0.3 and G49.0-0.3 are found to have a high fraction of YSOs and are therefore considered the youngest of all the near-infrared clusters in the W51 GMC. The estimated masses of the individual clusters, when summed, yield a total stellar mass of ~10^4 Msun in the W51 GMC, implying a star formation efficiency of 5-10%. These results in comparision with the CO observations of the W51 GMC, suggest for the first time, that star formation in the W51 GMC is likely triggered by a galactic spiral density wave.

We present new high-sensitivity high-resolution mm-wave observations of the Onsala 1 ultra-compact HII region that bring to light the internal structure of this massive star forming cloud. The 1.2 mm continuum map obtained with the IRAM 30-m radiotelescope (~11" resolution) shows a centrally peaked condensation of 1' size (~0.5 pc at the assumed distance of 1.8 kpc) which has been further investigated at higher resolution in the 3 mm continuum and in the emission lines of H13CO+ J=1--0 and SiO J=2--1 with the IRAM Plateau de Bure interferometer. The 3 mm data, with a resolution of ~5" X 4", displays a unresolved continuum source at the peak of the extended 1.2 mm emission and closely associated with the ultra-compact HII region. The H13CO+ map traces the central condensation in good agreement with previous NH_3 maps of Zheng et al. (1985). However, the velocity field of this central condensation, which was previously thought to arise in a rapidly rotating structure, is better explained in terms of the dense and compact component of a bipolar outflow. This interpretation is confirmed by SiO and CO observations of the full region. In fact, our new SiO data unveils the presence of multiple (at least 4) outflows in the region. In particular, there is an important center of outflow activity in the region at about 1' north of the UCHII region. Indeed the different outflows are related to different members of the Onsala 1 cluster. The data presented here support a scenario in which the phases of massive star formation begin much later in the evolution of a cluster and/or UCHII region last for much longer than 10^5 yrs.

We present sub-arcsec (FWHM ~ 0.5") J, H, K and L' images of a young stellar cluster associated with a candidate massive protostar IRAS22134+5834. The observations reveal a centrally symmetric, flattened cluster enclosing a central dark region. The central dark region is possibly a cavity within the flattened cluster. It is surrounded by a ring composed of 5 bright stars and the candidate massive protostar IRAS22134+5834. We construct JHKL' color-color and HK color-magnitude diagrams to identify the young stellar objects and estimate their spectral types. All the bright stars in the ring are found to have intrinsic infrared excess emission and are likely to be early to late B type stars. We estimate an average foreground extinction to the cluster of A_v \~ 5mag and individual extinctions to the bright stars in the range A_v ~ 20-40mag indicating possible cocoons surrounding each massive star. This ring of bright stars is devoid of any HII region. It is surrounded by an embedded cluster making this an example of a (proto)cluster that is in one of the dynamically least relaxed states. These observations are consistent with the recent non-axisymmetric calculations of Li & Nakamura, who present a star formation scenario in which a magnetically subcritical cloud fragments into multiple magnetically supercritical cores, leading to the formation of small stellar groups.

We present molecular line observations of the star-forming cloud around RNO6 along with a newly discovered nearby molecular cloud that we name RNO6NW. Both clouds display striking similarities in their cometary structures and overall kinematics. By using 13CO line observations, we estimate that these clouds have similar sizes (~4.5 pc) and masses (~200 solar masses). Both molecular clouds RNO6 and RNO6NW are active in star formation. From new high resolution near-IR narrowband images, we confirm that RNO6 hosts an embedded IR cluster that includes a Herbig Be star. A conspicuous H2 filament is found to delineate the dense cometary head of the globule. RNO6NW hosts at least two IR sources and a bipolar molecular outflow of ~0.9 pc of length and ~0.5 solar masses. We show that the cometary structure of both clouds has been created by the UV radiation from numerous OB stars lying ~1.5 degree to the north. Such OB stars are associated with the double cluster h and chi Persei, and are probably members of the PerOB1 association. Thus star formation inside these clouds has been very likely triggered by the Radiation Driven Implosion (RDI) mechanism. From comparison to RDI theoretical models, we find that the similar kinematics and morphology of both clouds is well explained if they are at a re-expansion phase. Triggered sequential star formation also explains the observed spatial distribution of the members of the near-IR cluster inside the RNO6 cloud, and the morphology of the H2 filament. We conclude that the RNO6 and RNO6NW clouds are high-mass counterparts to the cometary globules of smaller masses which have been studied up to now. Thus our observations demonstrate that the RDI mechanism can produce, not only low mass stars in small globules, but also intermediate mass stars and clusters in massive clouds.

We report the detection of the H2 v=1-0 S(1) line at 2.122um, from RNO91 in the L43 dark cloud, which is known to be a T Tauri star surrounded by a 1700AU disk structure (containing ices) and a weak outflow. The non-detection of the H2 v=2-1 S(1) line at 2.247um suggests shock excitation rather than fluorescence. The emission is extended spatially up to 9arcsec in the north-south direction. The line intensity peak (FWHM 3arcsec) corresponds to the star RNO91 which is embedded in a cocoon of gas and dust. The observed H2 emission from this cocoon may be attributed to embedded Herbig-Haro like knots. The H2 line flux in the central 2"X3"is estimated to be 7X10^{-14} ergs/sec/cm^{2}, which indicates a mass flow rate of 4X10^{-8} M_{sun}/yr. Furthermore, narrow band image taken through H2 1-0 S(1) filter is presented, which reveal a tilted disk and bipolar outflow structure that agrees with earlier observations and models. We show that this disk/outflow system is a unique case.