Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight Into Time Series Analysis project. As a result, we present a complete sample having 44, 998, 752 variable star candidates (VVV-CVSC), which include accurate individual coordinates, near-IR magnitudes (ZYJHKs), extinctions A(Ks), variability indices, periods, amplitudes, among other parameters to assess the science. Unfortunately, a side effect of having a highly complete sample, is also having a high level of contamination by non-variable (contamination ratio of non-variables to variables is slightly over 10:1). To deal with this, we also provide some flags and parameters that can be used by the community to de-crease the number of variable candidates without heavily decreasing the completeness of the sample. In particular, we cross-identified 339,601 of our sources with Simbad and AAVSO databases, which provide us with information for these objects at other wavelegths. This sub-sample constitutes a unique resource to study the corresponding near-IR variability of known sources as well as to assess the IR variability related with X-ray and Gamma-Ray sources. On the other hand, the other 99.5% sources in our sample constitutes a number of potentially new objects with variability information for the heavily crowded and reddened regions of the Galactic Plane and Bulge. The present results also provide an important queryable resource to perform variability analysis and to characterize ongoing and future surveys like TESS and LSST.

We present a comprehensive review of all observations of the eclipsing recurrent Nova LMC 1968 in the Large Magellanic Cloud which was previously observed in eruption in 1968, 1990, 2002, 2010, and most recently in 2016. We derive a probable recurrence time of $6.2 \pm 1.2$ years and provide the ephemerides of the eclipse. In the ultraviolet-optical-IR photometry the light curve shows high variability right from the first observation around two days after eruption. Therefore no colour changes can be substantiated. Outburst spectra from 2016 and 1990 are very similar and are dominated by H and He lines longward of 2000 Angstrom. Interstellar reddening is found to be E(B-V) = $0.07\pm0.01$. The super soft X-ray luminosity is lower than the Eddington luminosity and the X-ray spectra suggest the mass of the WD is larger than 1.3 M$_\odot$. Eclipses in the light curve suggest that the system is at high orbital inclination. On day four after the eruption a recombination wave was observed in Fe II ultraviolet absorption lines. Narrow line components are seen after day 6 and explained as being due to reionisation of ejecta from a previous eruption. The UV spectrum varies with orbital phase, in particular a component of the He II 1640 Angstrom emission line, which leads us to propose that early-on the inner WD Roche lobe might be filled with a bound opaque medium prior to the re-formation of an accretion disk. Both this medium and the ejecta can cause the delay in the appearance of the soft X-ray source.

We present high resolution spectroscopy of the yellow symbiotic star AG Draconis with ESPaDOnS at the {\it Canada-France-Hawaii Telescope}. Our analysis is focused on the profiles of Raman scattered \ion{O}{VI} features centered at 6825 \AA\ and 7082 \AA, which are formed through Raman scattering of \ion{O}{VI}$\lambda\lambda$1032 and 1038 with atomic hydrogen. These features are found to exhibit double component profiles with conspicuously enhanced red parts. Assuming that the \ion{O}{vi} emission region constitutes a part of the accretion flow around the white dwarf, Monte Carlo simulations for \ion{O}{VI} line radiative transfer are performed to find that the overall profiles are well fit with the accretion flow azimuthally asymmetric with more matter on the entering side than on the opposite side. As the mass loss rate of the giant component is increased, we find that the flux ratio $F(6825)/F(7082)$ of Raman 6825 and 7082 features decreases and that our observational data are consistent with a mass loss rate $\dot M\sim 2 \times 10^{-7} {\rm\ M_{\odot}\ yr^{-1}}$. We also find that additional bipolar components moving away with a speed $\sim 70{\rm\ km\ s^{-1}}$ provide considerably improved fit to the red wing parts of Raman features. The possibility that the two Raman profiles differ is briefly discussed in relation to the local variation of the \ion{O}{VI} doublet flux ratio.

Symbiotic stars (SySts) are long-period interacting binaries composed of a hot compact star, an evolved giant star, and a tangled network of gas and dust nebulae. They represent unique laboratories for studying a variety of important astrophysical problems, and have also been proposed as possible progenitors of SNIa. Presently, we know 257 SySts in the Milky Way and 69 in external galaxies. However, these numbers are still in striking contrast with the predicted population of SySts in our Galaxy. Because of other astrophysical sources that mimic SySt colors, no photometric diagnostic tool has so far demonstrated the power to unambiguously identify a SySt, thus making the recourse to costly spectroscopic follow-up still inescapable. In this paper we present the concept, commissioning, and science verification phases, as well as the first scientific results, of RAMSES II - a Gemini Observatory Instrument Upgrade Project that has provided each GMOS instrument at both Gemini telescopes with a set of narrow-band filters centered on the Raman OVI 6830 A band. Continuum-subtracted images using these new filters clearly revealed known SySts with a range of Raman OVI line strengths, even in crowded fields. RAMSES II observations also produced the first detection of Raman OVI emission from the SySt LMC 1 and confirmed Hen 3-1768 as a new SySt - the first photometric confirmation of a SySt. Via Raman OVI narrow-band imaging, RAMSES II provides the astronomical community with the first purely photometric tool for hunting SySts in the local Universe.

We report the discovery of VVV-WIT-07, an unique and intriguing variable source presenting a sequence of recurrent dips with a likely deep eclipse in July 2012. The object was found serendipitously in the near-IR data obtained by the VISTA Variables in the V\'ia L\'actea (VVV) ESO Public Survey. Our analysis is based on VVV variability, multicolor, and proper motion (PM) data. Complementary data from the VVV eXtended survey (VVVX) as well as archive data and spectroscopic follow-up observations aided in the analysis and interpretation of VVV-WIT-07. A search for periodicity in the VVV Ks-band light curve of VVV-WIT-07 results in two tentative periods at P~322 days and P~170 days. Colors and PM are consistent either with a reddened MS star or a pre-MS star in the foreground disk. The near-IR spectra of VVV-WIT-07 appear featureless, having no prominent lines in emission or absorption. Features found in the light curve of VVV-WIT-07 are similar to those seen in J1407 (Mamajek's object), a pre-MS K5 dwarf with a ring system eclipsing the star or, alternatively, to KIC 8462852 (Boyajian's star), an F3 IV/V star showing irregular and aperiodic dips in its light curve. Alternative scenarios, none of which is fully consistent with the available data, are also briefly discussed, including a young stellar object, a T Tauri star surrounded by clumpy dust structure, a main sequence star eclipsed by a nearby extended object, a self-eclipsing R CrB variable star, and even a long-period, high-inclination X-ray binary.

V735 Sgr was known as an enigmatic star with rapid brightness variations. Long-term OGLE photometry, brightness measurements in infrared bands, and recently obtained moderate resolution spectrum from the 6.5-m Magellan telescope show that this star is an active young stellar object of Herbig Ae/Be type.

The split main sequences (MSs) and extended MS turnoffs (eMSTOs) detected in a few young clusters have demonstrated that these stellar systems host multiple populations differing in a number of properties such as rotation and, possibly, age.We analyze Hubble Space Telescope photometry for thirteen clusters with ages between ~40 and ~1000 Myrs and of different masses. Our goal is to investigate for the first time the occurrence of multiple populations in a large sample of young clusters. We find that all the clusters exhibit the eMSTO phenomenon and that MS stars more massive than ~1.6 solar masses define a blue and red MS, with the latter hosting the majority of MS stars. The comparison between the observations and isochrones suggests that the blue MSs are made of slow-rotating stars, while the red MSs host stars with rotational velocities close to the breakup value. About half of the bright MS stars in the youngest clusters are H-alpha emitters. These Be stars populate the red MS and the reddest part of the eMSTO thus supporting the idea that the red MS is made of fast rotators. We conclude that the split MS and the eMSTO are a common feature of young clusters in both Magellanic Clouds. The phenomena of a split MS and an eMSTO occur for stars that are more massive than a specific threshold which is independent of the host-cluster mass. As a by-product, we report the serendipitous discovery of a young SMC cluster, GALFOR1.

We investigate the escape of Ly$\beta$ from emission nebulae with a significant population of excited hydrogen atoms in the level $n=2$, rendering them optically thick in H$\alpha$. The transfer of Ly$\beta$ line photons in these optically thick regions is complicated by the presence of another scattering channel leading to re-emission of H$\alpha$, alternating their identities between Ly$\beta$ and H$\alpha$. In this work, we develop a Monte Carlo code to simulate the transfer of Ly$\beta$ line photons incorporating the scattering channel into H$\alpha$. Both H$\alpha$ and Ly$\beta$ lines are formed through diffusion in frequency space, where a line photon enters the wing regime after a fairly large number of resonance scatterings with hydrogen atoms. Various line profiles of H$\alpha$ and Ly$\beta$ emergent from our model nebulae are presented. It is argued that the electron temperature is a critical parameter which controls the flux ratio of emergent Ly$\beta$ and H$\alpha$. Specifically for $T=3 \times 10^4{\rm\ K}$ and H$\alpha$ line center optical depth $\tau_\alpha=10$, the number flux ratio of emergent Ly$\beta$ and H$\alpha$ is $\sim 49$ percent, which is quite significant. We propose that the leaking Ly$\beta$ can be an interesting source for the formation of H$\alpha$ wings observed in many symbiotic stars and active galactic nuclei. Similar broad H$\alpha$ wings are also expected in Ly$\alpha$ emitting halos found in the early universe, which can be potentially probed by the {\it James Webb Telescope} in the future.

Regular intrinsic brightness variations observed in many stars are caused by pulsations. These pulsations provide information on the global and structural parameters of the star. The pulsation periods range from seconds to years, depending on the compactness of the star and properties of the matter that forms its outer layers. Here, we report the discovery of more than a dozen of previously unknown short-period variable stars: blue large-amplitude pulsators. These objects show very regular brightness variations with periods in the range of 20-40 min and amplitudes of 0.2-0.4 mag in the optical passbands. The phased light curves have a characteristic sawtooth shape, similar to the shape of classical Cepheids and RR Lyrae-type stars pulsating in the fundamental mode. The objects are significantly bluer than main sequence stars observed in the same fields, which indicates that all of them are hot stars. Follow-up spectroscopy confirms a high surface temperature of about 30,000 K. Temperature and colour changes over the cycle prove the pulsational nature of the variables. However, large-amplitude pulsations at such short periods are not observed in any known type of stars, including hot objects. Long-term photometric observations show that the variable stars are very stable over time. Derived rates of period change are of the order of 10^-7 per year and, in most cases, they are positive. According to pulsation theory, such large-amplitude oscillations may occur in evolved low-mass stars that have inflated helium-enriched envelopes. The evolutionary path that could lead to such stellar configurations remains unknown.

We present a detailed modeling of the two broad bands observed at 6825 \AA\ and 7082 \AA\ in Sanduleak's star, a controversial object in the Large Magellanic Cloud. These bands are known to originate from Raman-scattering of O VI $\lambda\lambda$ 1032 and 1038 photons with atomic hydrogen and are only observed in bona fide symbiotic stars. Our high-resolution spectrum obtained with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph at the Magellan-Clay Telescope reveals, quite surprisingly, that the profiles of the two bands look very different: while the Raman 6825 \AA\ band shows a single broad profile with a redward extended bump, the Raman 7082 \AA\ band exhibits a distinct triple-peak profile. Our model suggests that the O VI emission nebula can be decomposed into a red, blue and central emission regions from an accretion disk, a bipolar outflow and a further compact, optically thick region. We also perform Monte Carlo simulations with the aim of fitting the observed flux ratio $F(6825)/F(7082) \sim 4.5$, which indicate that the neutral region in Sanduleak's star is characterized by the column density $N_{HI} \sim 1 \times 10^{23} {\rm\ cm^{-2}}$.

Variable stars of RR Lyrae type are a prime tool to obtain distances to old stellar populations in the Milky Way, and one of the main aims of the Vista Variables in the Via Lactea (VVV) near-infrared survey is to use them to map the structure of the Galactic Bulge. Due to the large number of expected sources, this requires an automated mechanism for selecting RR Lyrae,and particularly those of the more easily recognized type ab (i.e., fundamental-mode pulsators), from the 10^6-10^7 variables expected in the VVV survey area. In this work we describe a supervised machine-learned classifier constructed for assigning a score to a K_s-band VVV light curve that indicates its likelihood of being ab-type RR Lyrae. We describe the key steps in the construction of the classifier, which were the choice of features, training set, selection of aperture and family of classifiers. We find that the AdaBoost family of classifiers give consistently the best performance for our problem, and obtain a classifier based on the AdaBoost algorithm that achieves a harmonic mean between false positives and false negatives of ~7% for typical VVV light curve sets. This performance is estimated using cross-validation and through the comparison to two independent datasets that were classified by human experts.

A dynamical characterization of dark matter dominated dwarf galaxies from their observed effects on galactic disks (i.e. Galactoseismology) has remained an elusive goal. Here, we present preliminary results from spectroscopic observations of three clustered Cepheid candidates identified from $K$-band light curves towards Norma. The average heliocentric radial velocity of these stars is $\sim$ 156 km/s, which is large and distinct from that of the Galaxy's stellar disk. These objects at $l \sim 333 ^\circ$ and $b \sim -1 ^\circ$ are therefore halo stars; using the $3.6~\micron$ period-luminosity relation of Type I Cepheids, they are at $\sim$ 73 kpc. Our ongoing $I$-band photometry indicates variability on the same time scale as the period determined from the $K_{s}$-band light curve. Distances determined from the $K$-band period-luminosity relation and the 3.6 $\micron$ period-luminosity relation are comparable. The observed radial velocity of these stars agrees roughly with predictions from dynamical models. If these stars are indeed members of the predicted dwarf galaxy that perturbed the outer HI disk of the Milky Way, this would mark the first application of Galactoseismology.

Powered by a supermassive black hole with an accretion disk, the spectra of active galactic nuclei (AGNs) are characterized by prominent emission lines including Balmer lines. The unification schemes of AGNs require the existence of a thick molecular torus that may hide the broad emission line region from the view of observers near the equatorial direction. In this configuration, one may expect that the far UV radiation from the central engine can be Raman scattered by neutral hydrogen to reappear around Balmer and Paschen emission lines which can be identified with broad wings. We produce H$\alpha$, H$\beta$ and Pa$\alpha$ wings using a Monte Carlo technique to investigate their properties. The neutral scattering region is assumed to be a cylindrical torus specified by the inner and outer radii and the height. While the covering factor of the scattering region affects the overall strengths of the wings, the wing widths are primarily dependent on the neutral hydrogen column density $N_{\rm HI}$ being roughly proportional to $N_{\rm HI}^{1/2}$. In particular, with $N_{\rm HI}=10^{23}{\rm\ cm^{-2}}$ the H$\alpha$ wings typically show a width $\sim 2\times 10^4{\rm\ km\ s^{-1}}$. We also find that H$\alpha$ and Pa$\alpha$ wing profiles are asymmetric with the red part stronger than the blue part and an opposite behavior is seen for H$\beta$ wings.

We present the results of a spectroscopic follow-up of various puzzling variable objects detected in the OGLE-III Galactic disk and bulge fields. The sample includes mainly short-period multi-mode pulsating stars that could not have been unambiguously classified as either delta Sct or beta Cep type stars based on photometric data only, also stars with irregular fluctuations mimicking cataclysmic variables and stars with dusty shells, and periodic variables displaying brightenings in their light curves that last for more than half of the period. The obtained low-resolution spectra show that all observed short-period pulsators are of delta Sct type, the stars with irregular fluctuations are young stellar objects, and the objects with regular brightenings are A type stars or very likely Ap stars with strong magnetic field responsible for the presence of bright caps around magnetic poles on their surface. We also took spectra of objects designated OGLE-GD-DSCT-0058 and OGLE-GD-CEP-0013. An estimated effective temperature of 33,000 K in OGLE-GD-DSCT-0058 indicates that it cannot be a delta Sct type variable. This very short-period (0.01962 d) high-amplitude (0.24 mag in the I-band) object remains a mystery. It may represent a new class of variable stars. The spectrum of OGLE-GD-CEP-0013 confirms that this is a classical Cepheid despite a peculiar shape of its light curve. The presented results will help in proper classification of variable objects in the OGLE Galactic Variability Survey.

[ABRIDGED] $\omega$ Centauri (NGC 5139) contains large numbers of variable stars of different types and, in particular, more than a hundred RR Lyrae stars. We have conducted a variability survey of $\omega$ Cen in the NIR, using ESO's 4.1m Visible and Infrared Survey Telescope for Astronomy (VISTA). This is the first paper of a series describing our results. $\omega$ Cen was observed using VIRCAM mounted on VISTA. A total of 42 and 100 epochs in $J$ and $K_{\rm S}$, respectively, were obtained, distributed over a total timespan of 352 days. PSF photometry was performed, and periods of the known variable stars were improved when necessary using an ANOVA analysis. An unprecedented homogeneous and complete NIR catalogue of RR Lyrae stars in the field of $\omega$ Cen was collected, allowing us to study, for the first time, all the RR Lyrae stars associated to the cluster, except 4 located far away from the cluster center. Membership status, subclassifications between RRab and RRc subtypes, periods, amplitudes, and mean magnitudes were derived for all the stars in our sample. Additionally, 4 new RR Lyrae stars were discovered, 2 of them with high probability of being cluster members. The distribution of $\omega$ Cen stars in the Bailey (period-amplitude) diagram is also discussed. Reference lines in this plane, for both Oosterhoff type I (OoI) and II (OoII) components, are provided. In the present paper, we clarify the status of many (candidate) RR Lyrae stars that had been unclear in previous studies. This includes stars with anomalous positions in the color-magnitude diagram, uncertain periods or/and variability types, and possible field interlopers. We conclude that $\omega$ Cen hosts a total of 88 RRab and 101 RRc stars, for a grand total of 189 likely members. We confirm that most RRab stars in the cluster belong to an OoII component, as previously found using visual data.

Stellar variability in the near-infrared (NIR) remains largely unexplored. The exploitation of public science archives with data-mining methods offers a perspective for the time-domain exploration of the NIR sky. We perform a comprehensive search for stellar variability using the optical-NIR multi-band photometric data in the public Calibration Database of the WFCAM Science Archive (WSA), with the aim of contributing to the general census of variable stars, and to extend the current scarce inventory of accurate NIR light curves for a number of variable star classes. We introduce new variability indices designed for multi-band data with correlated sampling, and apply them for pre-selecting variable star candidates, i.e., light curves that are dominated by correlated variations, from noise-dominated ones. Pre-selection criteria are established by robust numerical tests for evaluating the response of variability indices to colored noise characteristic to the data. We find 275 periodic variable stars and an additional 44 objects with suspected variability with uncertain periods or apparently aperiodic variation. Only 44 of these objects had been previously known, including 11 RR~Lyrae stars in the outskirts of the globular cluster M3 (NGC~5272). We provide a preliminary classification of the new variable stars that have well-measured light curves, but the variability types of a large number of objects remain ambiguous. We classify most of the new variables as contact binary stars, but we also find several pulsating stars, among which 34 are probably new field RR~Lyrae and 3 are likely Cepheids. We also identify 32 highly reddened variable objects close to previously known dark nebulae, suggesting that these are embedded young stellar objects. We publish our results and all light-curve data as the WFCAM Variable Star Catalog.

The VISTA Variables in the V\'ia L\'actea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

NGC1851 is surrounded by a stellar component that extends more than ten times beyond the tidal radius. Although the nature of this stellar structure is not known, it has been suggested to be a sparse halo of stars or associated with a stellar stream. We analyse the nature of this intriguing stellar component surrounding NGC1851 by investigating its radial velocities and chemical composition, in particular in comparison with those of the central cluster analysed in a homogeneous manner. In total we observed 23 stars in the halo with radial velocities consistent with NGC1851, and for 15 of them we infer [Fe/H] abundances. Our results show that: (i) stars dynamically linked to NGC1851 are present at least up to ~2.5 tidal radii, supporting the presence of a halo of stars surrounding the cluster; (ii) apart from the NGC1851 radial velocity-like stars, our observed velocity distribution agrees with that expected from Galactic models, suggesting that no other sub-structure (such as a stream) at different radial velocities is present in our field; (iii) the chemical abundances for the s-process elements Sr and Ba are consistent with the s-normal stars observed in NGC1851; (iv) all halo stars have metallicities, and abundances for the other studied elements Ca, Mg and Cr, consistent with those exhibited by the cluster. The complexity of the whole NGC1851 cluster+halo system may agree with the scenario of a tidally-disrupted dwarf galaxy in which NGC1851 was originally embedded.

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.

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 an abundance analysis of 96 horizontal branch (HB) stars in NGC2808, a globular cluster exhibiting a complex multiple stellar population pattern. These stars are distributed in different portions of the HB and cover a wide range of temperature. By studying the chemical abundances of this sample, we explore the connection between HB morphology and the chemical enrichment history of multiple stellar populations. For stars lying on the red HB, we use GIRAFFE and UVES spectra to determine Na, Mg, Si, Ca, Sc, Ti, Cr, Mn, Fe, Ni, Zn, Y, Ba, and Nd abundances. For colder, blue HB stars, we derive abundances for Na, primarily from GIRAFFE spectra. We were also able to measure direct NLTE He abundances for a subset of these blue HB stars with temperature higher than ~9000 K. Our results show that: (i) HB stars in NGC2808 show different content in Na depending on their position in the color-magnitude diagram, with blue HB stars having higher Na than red HB stars; (ii) the red HB is not consistent with an uniform chemical abundance, with slightly warmer stars exhibiting a statistically significant higher Na content; and (iii) our subsample of blue HB stars with He abundances shows evidence of enhancement with respect to the predicted primordial He content by Delta(Y)=+0.09+-0.01. Our results strongly support theoretical models that predict He enhancement among second generation(s) stars in globular clusters and provide observational constraints on the second-parameter governing HB morphology.

Symbiotic stars are long-orbital-period interacting-binaries characterized by extended emission over the whole electromagnetic range and by complex photometric and spectroscopic variability. In this paper, the first of a series, we present OGLE light curves of all the confirmed symbiotic stars in the Large Magellanic Cloud, with one exception. By careful visual inspection and combined time-series analysis techniques, we investigate for the first time in a systematic way the photometric properties of these astrophysical objects, trying in particular to distinguish the nature of the cool component (e.g., Semi-Regular Variable vs. OGLE Small-Amplitude Red Giant), to provide its first-order pulsational ephemerides, and to link all this information with the physical parameters of the binary system as a whole. Among the most interesting results, there is the discovery of a 20-year-long steady fading of Sanduleak's star, a peculiar symbiotic star known to produce the largest stellar jet ever discovered. We discuss by means of direct examples the crucial need for long-term multi-band observations to get a real understanding of symbiotic and other interacting binary stars. We eventually introduce BOMBOLO, a multi-band simultaneous imager for the SOAR 4m Telescope, whose design and construction we are currently leading.

Near-IR data of Classical Novae contain useful information about the ejected gas mass and the thermal emission by dust formed during eruption, and provide independent methods to classify the objects according to the colour of their progenitors, and the fading rate and features seen after eruption. The VISTA Variables in the Via Lactea survey (VVV) is a near-IR ESO Public Survey mapping the MW bulge and southern plane. Data taken during 2010-2011 covered the entire area in the JHKs bands plus some epochs in Ks-band of the ongoing VVV variability campaign. We used the novae list provided by VSX/AAVSO catalogue to search for all objects within the VVV area. We used the VVV data to create a near-IR catalogue of the known Galactic novae in the 562 sq.deg. area covered by VVV. The VVV near-IR catalogue of novae contains JHKs photometry of 93 objects completed as of December 2012. We also checked in the ongoing VVV variability campaign for the light-curves of novae that erupted in the last years. VVV images can also be used to discover and study novae by searching for the expanding shell. Since objects are seen at different distances and reddening levels, the colour-magnitude and colour-colour diagrams show the novae spread in magnitude as well as in colour. Dereddened colours and reddening-free indices were used with caution and cannot be a good approach in all cases since the distance and spectral features prevent more conclusive results for some extreme objects. Light-curves for some recent novae are presented. Thanks to its high spatial resolution in the near-IR, and large Ks-range, the VVV survey can be a major contributor for the search and study of novae in the most crowded and high-extinction regions of the Milky Way. The VVV survey area contains ~35 of all known novae in the Galaxy.

We present photometric and spectroscopic analysis of AE For -- a detached eclipsing binary composed of two late K dwarfs. The masses of the components are found to be 0.6314 +- 0.0035 and 0.6197 +- 0.0034 Msun and the radii to be 0.67 +- 0.03 and 0.63$ +- 0.03 Rsun for primary and secondary component, respectively. Both components are significantly oversized compared to theoretical models, which we attribute to their high activity. They show Halpha, Hbeta, Hgamma, Ca H and Ca K lines in emission, and are heavily spotted, causing large variations of the light curve.

We report on the discovery of large-amplitude flickering from V648 Car (= SS73-17), a poorly studied object listed amongst the very few hard X-ray emitting symbiotic stars. We performed milli-magnitude precision optical photometry with the Swope Telescope at the Las Campanas Observatory, Chile, and found that V648 Car shows large U-band variability over time scales of minutes. To our knowledge, it is amongst the largest flickering of a symbiotic star ever reported. Our finding supports the hypothesis that symbiotic WDs producing hard X-rays are predominantly powered by accretion, rather than quasi-steady nuclear burning, and have masses close to the Chandrasekhar limit. No significant periodicity is evident from the flickering light curve. The ASAS long-term V light curve suggests the presence of a tidally distorted giant accreting via Roche Lobe overflow, and a binary period of about 520 days. On the basis of the outstanding physical properties of V648 Car as hinted by its fast and long-term optical variability, as well as by its nature as hard X-ray emitter, we therefore call for simultaneous follow-up observations in different bands, ideally combined with time-resolved optical spectroscopy.

The ESO Public Survey VISTA Variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for ~5 years. In this paper we describe the progress of the survey observations in the first observing season, the observing strategy and quality of the data obtained. The observations are carried out on the 4-m VISTA telescope in the ZYJHKs filters. In addition to the multi-band imaging the variability monitoring campaign in the Ks filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHKs filters taken in the 2010 observing season. The typical image quality is ~0.9-1.0". The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHKs images in the disk area and 90% of the JHKs images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28x10^8 stellar sources in the bulge and 1.68x10^8 in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63x10^8 stellar sources. About 10% of these are double detections due to overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHKs bands extend typically ~4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for Ks=15-18 mag has rms ~35-175 mas.

Highly-collimated gas ejections are among the most dramatic structures in the Universe, observed to emerge from very different astrophysical systems - from active galactic nuclei down to young brown dwarf stars. Even with the huge span in spatial scales, there is convincing evidence that the physics at the origin of the phenomenon, namely the acceleration and collimation mechanisms, is the same in all classes of jets. Here we report on the discovery of a giant, highly-collimated jet from Sanduleak's star in the Large Magellanic Cloud (LMC). With a physical extent of 14 parsecs at the distance of the LMC, it represents the largest stellar jet ever discovered, and the first resolved stellar jet beyond the Milky Way. The kinematics and extreme chemical composition of the ejecta from Sanduleak's star bear strong resemblance with the low-velocity remnants of SN1987A and with the outer filaments of the most famous supernova progenitor candidate, i.e., eta Carinae. Moreover, the precise knowledge of the jet's distance implies that it will be possible to derive accurate estimates of most of its physical properties. Sanduleak's bipolar outflow will thus become a crucial test-bed for future theoretical modeling of astrophysical jets.

Vista Variables in the Via Lactea (VVV) is an ESO Public Survey that is performing a variability survey of the Galactic bulge and part of the inner disk using ESO's Visible and Infrared Survey Telescope for Astronomy (VISTA). The survey covers 520 deg^2 of sky area in the ZYJHK_S filters, for a total observing time of 1929 hours, including ~ 10^9 point sources and an estimated ~ 10^6 variable stars. Here we describe the current status of the VVV Survey, in addition to a variety of new results based on VVV data, including light curves for variable stars, newly discovered globular clusters, open clusters, and associations. A set of reddening-free indices based on the ZYJHK_S system is also introduced. Finally, we provide an overview of the VVV Templates Project, whose main goal is to derive well-defined light curve templates in the near-IR, for the automated classification of VVV light curves.

We are performing a search for symbiotic stars using IPHAS, the INT Halpha survey of the northern Galactic plane, and follow-up observations. Candidate symbiotic stars are selected on the basis of their IPHAS and near-IR colours, and spectroscopy and photometry are obtained to determine their nature. We present here observations of the symbiotic star candidate IPHAS J205836.43+503307.2. The optical spectrum shows the combination of a number of emission lines, among which are the high-excitation species of [OIII], HeII, [Ca V], and [Fe VII], and a red continuum with the features of a star at the cool end of the carbon star sequence. The nebular component is spatially resolved: the analysis of the spatial profile of the [NII]6583 line in the spectrum indicates a linear size of ~2.5 arcsec along the east-west direction. Its velocity structure suggests an aspherical morphology. The near-infrared excess of the source, which was especially strong in 1999, indicated that a thick circumstellar dust shell was also present in the system. The carbon star has brightened in the last decade by two to four magnitudes at red and near-infrared wavelengths. Photometric monitoring during a period of 60 days from November 2010 to January 2011 reveals a slow luminosity decrease of 0.2 magnitudes. From the observed spectrophotometric properties and variability, we conclude that the source is a new Galactic symbiotic star of the D-type, of the rare kind that contains a carbon star, likely a carbon Mira. Only two other systems of this type are known in the Galaxy.

The modeling of UV and optical spectra emitted from the symbiotic system AG Draconis, adopting collision of the winds, predicts soft X-ray bremsstrahlung from nebulae downstream of the reverse shock with velocities > 150 km/s and intensities comparable to those of the white dwarf black body flux. At outbursts, the envelop of debris, which corresponds to the nebula downstream of the high velocity shocks (700-1000 km/s) accompanying the blast wave, absorbs the black body soft X-ray flux from the white dwarf, explains the broad component of the H and He lines, and leads to low optical-UV-X-ray continuum fluxes. The high optical-UV flux observed at the outbursts is explained by bremsstrahlung downstream of the reverse shock between the stars. The depletion of C, N, O, and Mg relative to H indicates that they are trapped into dust grains and/or into diatomic molecules, suggesting that the collision of the wind from the white dwarf with the dusty shells, ejected from the red giant with about 1 year periodicity, leads to the U-band fluctuations during the major bursts.

We have collected continuum data of a sample of D-type symbiotic stars. By modelling their spectral energy distribution in a colliding-wind theoretical scenario we have found the common characteristics to all the systems: 1) at least two dust shells are clearly present, one at \sim 1000 K and the other at \sim 400 K; they dominate the emission in the IR; 2) the radio data are explained by thermal self-absorbed emission from the reverse shock between the stars; while 3) the data in the long wavelength tail come from the expanding shock outwards the system; 4) in some symbiotic stars, the contribution from the WD in the UV is directly seen. Finally, 5) for some objects soft X-ray emitted by bremsstrahlung downstream of the reverse-shock between the stars are predicted. The results thus confirm the validity of the colliding wind model and the important role of the shocks. The comparison of the fluxes calculated at the nebula with those observed at Earth reveals the distribution throughout the system of the different components, in particular the nebulae and the dust shells. The correlation of shell radii with the orbital period shows that larger radii are found at larger periods. Moreover, the temperatures of the dust shells regarding the sample are found at 1000 K and <=400 K, while, in the case of late giants, they spread more uniformly throughout the same range.

MASYS is the AKARI spectroscopic survey of Symbiotic Stars in the Magellanic Clouds, and one of the European Open Time Observing Programmes approved for the AKARI (Post-Helium) Phase-3. It is providing the first ever near-IR spectra of extragalactic symbiotic stars. The observations are scheduled to be completed in July 2009.

We investigate the symbiotic star BI Crucis through a comprehensive and self-consistent analysis of the spectra emitted in three different epochs: 60's, 70's, and late 80's. In particular, we would like to find out the physical conditions in the shocked nebula and in the dust shells, as well as their location within the symbiotic system, by exploiting both photometric and spectroscopic data from radio to UV. We suggest a model which, on the basis of optical imaging, emission line ratios and spectral energy distribution profile, is able to account for collision of the winds, formation of lobes and jets by accretion onto the WD, as well as for the interaction of the blast wave from a past, unrecorded outburst with the ISM. We have found that the spectra observed throughout the years show the marks of the different processes at work within BI Cru, perhaps signatures of a post-outburst evolution. We then call for new infrared and millimeter observations, potentially able to resolve the inner structure of the symbiotic nebula.

Context. We analyse the line and continuum spectra of the symbiotic system CH Cygni. Aims. To show that the colliding-wind model is valid to explain this symbiotic star at different phases. Methods. Peculiar observed features such as flickering, radio variation, X-ray emission, as well as the distribution of the nebulae and shells throughout the system are investigated by modelling the spectra at different epochs. The models account consistently for shock and photoionization and are constrained by absolute fluxes. Results. We find that the reverse shock between the stars leads to the broad lines observed during the active phases, as well as to radio and hard X-ray emission, while the expanding shock is invoked to explain the data during the transition phases.

Context. In 1985, at the end of the active phase 1977-1986, a broad (4000 km/s) Ly alpha line appeared in the symbiotic system CH Cygni that had never been observed previously. Aims. In this work we investigate the origin of this anomalous broad Ly alpha line. Methods. We suggest a new interpretation of the broad Ly alpha based on the theory of charge transfer reactions between ambient hydrogen atoms and post-shock protons at a strong shock front. Results. We have found that the broad Ly alpha line originated from the blast wave created by the outburst, while the contemporary optical and UV lines arose from the nebula downstream of the expanding shock in the colliding wind scenario.

In this paper we analyse the spectra of D-type SS H1-36 within a colliding-wind scenario. We aim to analyse the properties of this object taking into account the observational data along the whole electromagnetic spectrum, in order to derive a self-consistent picture able to interpret the nature of the system as a whole. After constraining the relative physical conditions by modelling more than 40 emission lines from radio to UV, we are able to explain the continuum spectral energy distribution by taking into account all the emitting contributions arising from both the stars, the dust shells and the gaseous nebulae. A comprehensive model of the radio spectra allows to reproduce the different slopes of the radio profile and the turnover frequency, as well as the different size of the observed shocked envelope at different frequencies in the light of the different contributions from the expanding and reverse nebulae. The IR continuum unveils the presence of two dust shells with different radii and temperatures, which might be a distinctive feature of D-type symbiotic systems as a class of objects. The broad profiles of IR lines direct us to investigate whether an X-ray jet may be present. This insight leads us to indicate H1-36 as a promising X-ray target and to encourage observations and studies which consistently take into account the complex nature of symbiotic stars throughout the whole electromagnetic spectrum.

We present a comprehensive and self-consistent modelling of the D' type symbiotic star (SS) HD330036 from radio to UV. Within a colliding-wind scenario, we analyse the continuum, line and dust spectra by means of SUMA, a code that simulates the physical conditions of an emitting gaseous cloud under the coupled effect of ionization from an external radiation source and shocks. We find that the UV lines are emitted from high density gas between the stars downstream of the reverse shock, while the optical lines are emitted downstream of the shock propagating outwards the system. As regards with the continuum SED, three shells are identified in the IR, at 850K, 320 K and 200 K with radii r = 2.8 10^13 cm, 4 10^14$ cm, and 10^15 cm, respectively, adopting a distance to Earth d=2.3 kpc: interestingly, all these shells appear to be circumbinary. The analysis of the unexploited ISO-SWS spectrum reveals that both PAHs and crystalline silicates coexist in HD330036, with PAHs associated to the internal shell at 850 K, and crystalline silicates stored into the cool shells at 320 K and 200 K. Strong evidence that crystalline silicates are shaped in a disk-like structure is derived on the basis of the relative band strengths. Finally, we suggest that shocks can be a reliable mechanism in activating the annealing and the consequent crystallization processes. We show that a consistent interpretation of gas and dust spectra emitted by SS can be obtained by models which accounts for the coupled effect of the photoionizing flux and of shocks. The VLTI/MIDI proposal recently accepted by ESO aims to verify and better constrain some of our results by means of IR interferometric observations.

We investigate the IR spectral features of a sample of D-type symbiotic stars. Analyzing unexploited ISO-SWS data, deriving the basic observational parameters of dust bands and comparing them with respect to those observed in other astronomical sources, we try to highlight the effect of environment on grain chemistry and physic. We find strong amorphous silicate emission bands at 10 micron and 18 micron in a large fraction of the sample. The analysis of the 10 micron band, along with a direct comparison with several astronomical sources, reveals that silicate dust in symbiotic stars shows features between the characteristic circumstellar environments and the interstellar medium. This indicates an increasing reprocessing of grains in relation to specific symbiotic behavior of the objects. A correlation between the central wavelength of the 10 and 18 micron dust bands is found. By the modeling of IR spectral lines we investigate also dust grains conditions within the shocked nebulae. Both the unusual depletion values and the high sputtering efficiency might be explained by the formation of SiO moleculae, which are known to be a very reliable shock tracer. We conclude that the signature of dust chemical disturbance due to symbiotic activity should be looked for in the outer, circumbinary, expanding shells where the environmental conditions for grain processing might be achieved. Symbiotic stars are thus attractive targets for new mid-infrared and mm observations.