We study the relationship between the mid-infrared and sub-mm variability of deeply embedded protostars using the multi-epoch data from the Wide Infrared Survey Explorer ($WISE$/NEOWISE) and the ongoing James Clerk Maxwell Telescope (JCMT) transient survey. Our search for signs of stochastic (random) and/or secular (roughly monotonic in time) variability in a sample of 59 young stellar objects (YSOs) revealed that 35 are variable in at least one of the two surveys. This variability is dominated by secular changes. Of those objects with secular variability, 14 objects ($22\%$ of the sample) show correlated secular variability over mid-IR and sub-mm wavelengths. Variable accretion is the likely mechanism responsible for this type of variability. Fluxes of YSOs that vary in both wavelengths follow a relation of $\log_{10} F_{4.6}(t)=\eta \log_{10} F_{850}(t)$ between the mid-IR and sub-mm, with $\eta=5.53\pm0.29$. This relationship arises from the fact that sub-mm fluxes respond to the dust temperature in the larger envelope whereas the mid-IR emissivity is more directly proportional to the accretion luminosity. The exact scaling relation, however, depends on the structure of the envelope, the importance of viscous heating in the disc, and dust opacity laws.

A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $\mu$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC.

FU Orionis-type objects (FUors) are young low-mass stars undergoing powerful accretion outbursts. The increased accretion is often accompanied by collimated jets and energetic, large-scale molecular outflows. The extra heating during the outburst may also induce detectable geometrical, chemical, and mineralogical changes in the circumstellar material, affecting possible planet formation around these objects. V346 Nor is a southern FUor with peculiar spectral characteristics. Decades after the beginning of its outburst, it unexpectedly underwent a fading event around 2010 due to a decrease in the mass accretion rate onto the star by at least two orders of magnitude. Here we present optical and near-infrared photometry and spectroscopy obtained after the minimum. Our light curves show a gradual re-brightening of V346 Nor, with its Ks-band brightness only 1.5 mag below the outburst brightness level. Our VLT/XSHOOTER spectroscopic observations display several strong forbidden emission lines towards the source from various metals and molecular hydrogen, suggesting the launch of a new jet. Our N-band spectrum obtained with VLT/VISIR outlines a deeper silicate absorption feature than before, indicating that the geometry of the circumstellar medium has changed in the post-outburst period compared to peak brightness.

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.

We have determined the rate of large accretion events in class I and II young stellar objects (YSOs) by comparing the all-sky digitised photographic plate surveys provided by SuperCOSMOS with the latest data release from Gaia (DR2). The long mean baseline of 55 years along with a large sample of class II YSOs ($\simeq$15,000) allows us to study approximately 1 million YSO-years. We find 139 objects with $\Delta R\geq1$~mag, most of which are found at amplitudes between 1 and 3 mag. The majority of YSOs in this group show irregular variability or long-lasting fading events, which is best explained as hot spots due to accretion or by variable extinction. There is a tail of YSOs at $\Delta R\geq3$~mag and they seem to represent a different population. Surprisingly many objects in this group show high-amplitude irregular variability over timescales shorter than 10 years, in contrast with the view that high-amplitude objects always have long outbursts. However, we find 6 objects that are consistent with undergoing large, long lasting accretion events, 3 of them previously unknown. This yields an outburst recurrence timescale of 112 kyr, with a 68\% confidence interval [74 to 180] kyr. This represents the first robust determination of the outburst rate in class II YSOs and shows that YSOs in their planet-forming stage do in fact undergo large accretion events, and with timescales of $\simeq$100,000 years. In addition, we find that outbursts in the class II stage are $\simeq$10 times less frequent than during the class I stage.

We report on the source Gaia 17bpi and identify it as a new, ongoing FU Ori type outburst, associated with a young stellar object. The optical lightcurve from Gaia exhibited a 3.5 mag rise with the source appearing to plateau in mid/late 2018. Mid-infrared observations from NEOWISE also show a $>$3 mag rise that occurred in two stages, with the second one coincident with the optical brightening, and the first one preceding the optical brightening by $\sim$1.5 years. We model the outburst as having started between October and December of 2014. This wavelength-dependent aspect of young star accretion-driven outbursts has never been documented before. Both the mid-infrared and the optical colors of the object become bluer as the outburst proceeds. Optical spectroscopic characteristics in the outburst phase include: a GK-type absorption spectrum, strong wind/outflow in e.g. Mgb, NaD, H$\alpha$, KI, OI, and CaII profiles, and detection of LiI 6707 \AA. The infrared spectrum in the outburst phase is similar to that of an M-type spectrum, notably exhibiting prominent $H_2O$ and $^{12}$CO (2-0) bandhead absorption in the K-band, and likely HeI wind in the Y-band. The new FU Ori source Gaia 17bpi is associated with a little-studied dark cloud in the galactic plane, located at a distance of 1.27 kpc.

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.

We present time-series $BVI$ photometry for the Galactic globular cluster NGC 6402 (M14). The data consists of $\sim$137 images per filter, obtained using the 0.9m and 1.0m SMARTS telescopes at the Cerro Tololo Inter-American Observatory. The images were obtained during two observing runs in 2006-2007. The image-subtraction package ISIS, along with DAOPHOT II/ALLFRAME, were used to perform crowded-field photometry and search for variable stars. We identified 130 variables, 8 of which are new discoveries. The variable star population is comprised of 56 ab-type RR Lyrae stars, 54 c-type RR Lyrae, 6 type II Cepheids, 1 W UMa star, 1 detached eclipsing binary, and 12 long-period variables. We provide Fourier decomposition parameters for the RR Lyrae, and discuss the physical parameters and photometric metallicity derived therefrom. The M14 distance modulus is also discussed, based on different approaches for the calibration of the absolute magnitudes of RR Lyrae stars. The possible presence of second-overtone RR Lyrae in M14 is critically addressed, with our results arguing against this possibility. By considering all of the RR Lyrae stars as members of the cluster, we derive $\langle P_{ab}\rangle = 0.589$ d. This, together with the position of the RR Lyrae stars of both Bailey types in the period-amplitude diagram, suggests an Oosterhoff-intermediate classification for the cluster. Such an intermediate Oosterhoff type is much more commonly found in nearby extragalactic systems, and we critically discuss several other possible indications that may point to an extragalactic origin for this cluster.

We present VIRAC version 1, a near-infrared proper motion and parallax catalogue of the VISTA VVV survey for 312,587,642 unique sources averaged across all overlapping pawprint and tile images covering 560 deg$^2$ of the bulge of the Milky Way and southern disk. The catalogue includes 119 million high quality proper motion measurements, of which 47 million have statistical uncertainties below 1 mas yr$^{-1}$. In the 11$<K_s<$14 magnitude range the high quality motions have a median uncertainty of 0.67 mas yr$^{-1}$. The catalogue also includes 6,935 sources with quality-controlled 5 $\sigma$ parallaxes with a median uncertainty of 1.1 mas. The parallaxes show reasonable agreement with the TYCHO-Gaia Astrometric Solution (TGAS), though caution is advised for data with modest significance. The SQL database housing the data is made available via the web. We give example applications for studies of Galactic structure, nearby objects (low mass stars and brown dwarfs, subdwarfs, white dwarfs) and kinematic distance measurements of YSOs. Nearby objects discovered include LTT 7251 B, an L7 benchmark companion to a G dwarf with over 20 published elemental abundances, a bright L sub-dwarf, VVV 1256-6202, with extremely blue colours and nine new members of the 25 pc sample. We also demonstrate why this catalogue remains useful in the era of Gaia. Future versions will be based on profile fitting photometry, use the Gaia absolute reference frame and incorporate the longer time baseline of the VVV extended survey (VVVX).

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.

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.

The purpose of this research is to study the connection of global properties of eight young stellar clusters projected in the Vista Variables in the Via Lactea (VVV) ESO Large Public Survey disk area and their young stellar object population. The analysis in based on the combination of spectroscopic parallax-based reddening and distance determinations with main sequence and pre-main sequence ishochrone fitting to determine the basic parameters (reddening, age, distance) of the sample clusters. The lower mass limit estimations show that all clusters are low or intermediate mass (between 110 and 1800 Mo), the slope Gamma of the obtained present-day mass functions of the clusters is close to the Kroupa initial mass function. On the other hand, the young stellar objects in the surrounding cluster's fields are classified by low resolution spectra, spectral energy distribution fit with theoretical predictions, and variability, taking advantage of multi-epoch VVV observations. All spectroscopically confirmed young stellar objects (except one) are found to be massive (more than 8 Mo). Using VVV and GLIMPSE color-color cuts we have selected a large number of new young stellar object candidates, which are checked for variability and 57% are found to show at least low-amplitude variations. In few cases it was possible to distinguish between YSO and AGB classification on the basis of the light curves.

We report the discovery of an L dwarf companion to the A3V star \beta{} Circini. VVV J151721.49-585131.5, or \beta{} Cir B, was identified in a proper motion and parallax catalogue of the Vista Variables in the V\'{i}a L\'{a}ctea survey as having near infrared luminosity and colour indicative of an early L dwarf, and a proper motion and parallax consistent with that of \beta{} Cir. The projected separation of $\sim$3.6' corresponds to $6656$ au, which is unusually wide. The most recent published estimate of the age of the primary combined with our own estimate based on newer isochrones yields an age of $370-500$ Myr. The system therefore serves as a useful benchmark at an age greater than that of the Pleiades brown dwarfs and most other young L dwarf benchmarks. We have obtained a medium resolution echelle spectrum of the companion which indicates a spectral type of L1.0$\pm$0.5 and lacks the typical signatures of low surface gravity seen in younger brown dwarfs. This suggests that signs of low surface gravity disappear from the spectra of early L dwarfs by an age of $\sim370-500$ Myr, as expected from theoretical isochrones. The mass of \beta{} Cir B is estimated from the BHAC15 isochrones as $0.056\pm0.007$ M$_{\odot}$.

Context. The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is a variability survey of the Milky Way bulge and an adjacent section of the disk carried out from 2010 on ESO Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV will eventually deliver a deep near-IR atlas with photometry and positions in five passbands (ZYJHK_S) and a catalogue of 1-10 million variable point sources - mostly unknown - which require classifications. Aims. The main goal of the VVV Templates Project, that we introduce in this work, is to develop and test the machine-learning algorithms for the automated classification of the VVV light-curves. As VVV is the first massive, multi-epoch survey of stellar variability in the near-infrared, the template light-curves that are required for training the classification algorithms are not available. In the first paper of the series we describe the construction of this comprehensive database of infrared stellar variability. Methods. First we performed a systematic search in the literature and public data archives, second, we coordinated a worldwide observational campaign, and third we exploited the VVV variability database itself on (optically) well-known stars to gather high-quality infrared light-curves of several hundreds of variable stars. Results. We have now collected a significant (and still increasing) number of infrared template light-curves. This database will be used as a training-set for the machine-learning algorithms that will automatically classify the light-curves produced by VVV. The results of such an automated classification will be covered in forthcoming papers of the series.

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.

The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is an ongoing time-series, near-infrared (IR) survey of the Galactic bulge and an adjacent portion of the inner disk, covering 562 square degrees of the sky, using ESO's VISTA telescope. The survey has provided superb multi-color photometry in 5 broadband filters ($Z$, $Y$, $J$, $H$, and $K_s$), leading to the best map of the inner Milky Way ever obtained, particularly in the near-IR. The main variability part of the survey, which is focused on $K_s$-band observations, is currently underway, with bulge fields having been observed between 31 and 70 times, and disk fields between 17 and 36 times. When the survey is complete, bulge (disk) fields will have been observed up to a total of 100 (60) times, providing unprecedented depth and time coverage. Here we provide a first overview of stellar variability in the VVV data, including examples of the light curves that have been collected thus far, scientific applications, and our efforts towards the automated classification of VVV light curves.

We present BV photometry of the Galactic globular cluster NGC 6402 (M14), based on 65 V frames and 67 B frames, reaching two magnitudes below the turn-off level. This represents, to the best of our knowledge, the deepest color-magnitude diagram (CMD) of NGC 6402 available in the literature. Statistical decontamination of field stars as well as differential reddening corrections are performed in order to derive a precise ridgeline and derive physical parameters of the cluster therefrom. We discuss previous attempts to derive a reddening value for the cluster, and argue in favor of a value E(B-V) = 0.57 +/- 0.02, which is significantly higher than indicated by either the Burstein & Heiles or Schlegel et al. (corrected according to Bonifacio et al.) interstellar dust maps. Differential reddening across the face of the cluster, which we find to be present at the level of Delta E(B-V) ~ 0.17 mag, is taken into account in our analysis. We measure several metallicity indicators based on the position of the red giant branch (RGB) in the cluster CMD. These give a metallicity of [Fe/H] = -1.38 +/- 0.07 in the Zinn & West scale and [Fe/H] = -1.28 +/- 0.08 in the new Carretta et al. (UVES) scale. We also provide measurements of other important photometric parameters for this cluster, including the position of the RGB luminosity function "bump" and the horizontal branch (HB) morphology. We compare the NGC 6402 ridgeline with the one for NGC 5904 (M5) derived by Sandquist et al., and find evidence that NGC 6402 and M5 have approximately the same age, to within the uncertainties -- although the possibility that M14 may be slighlty older cannot be ruled out.