Most time-series models assume that the data come from observations that are equally spaced in time. However, this assumption does not hold in many diverse scientific fields, such as astronomy, finance, and climatology, among others. There are some techniques that fit unequally spaced time series, such as the continuous-time autoregressive moving average (CARMA) processes. These models are defined as the solution of a stochastic differential equation. It is not uncommon in astronomical time series, that the time gaps between observations are large. Therefore, an alternative suitable approach to modeling astronomical time series with large gaps between observations should be based on the solution of a difference equation of a discrete process. In this work we propose a novel model to fit irregular time series called the complex irregular autoregressive (CIAR) model that is represented directly as a discrete-time process. We show that the model is weakly stationary and that it can be represented as a state-space system, allowing efficient maximum likelihood estimation based on the Kalman recursions. Furthermore, we show via Monte Carlo simulations that the finite sample performance of the parameter estimation is accurate. The proposed methodology is applied to light curves from periodic variable stars, illustrating how the model can be implemented to detect poor adjustment of the harmonic model. This can occur when the period has not been accurately estimated or when the variable stars are multiperiodic. Last, we show how the CIAR model, through its state space representation, allows unobserved measurements to be forecast.

Time series observations are ubiquitous in astronomy, and are generated to distinguish between different types of supernovae, to detect and characterize extrasolar planets and to classify variable stars. These time series are usually modeled using a parametric and/or physical model that assumes independent and homoscedastic errors, but in many cases these assumptions are not accurate and there remains a temporal dependency structure on the errors. This can occur, for example, when the proposed model cannot explain all the variability of the data or when the parameters of the model are not properly estimated. In this work we define an autoregressive model for irregular discrete-time series, based on the discrete time representation of the continuous autoregressive model of order 1. We show that the model is ergodic and stationary. We further propose a maximum likelihood estimation procedure and assess the finite sample performance by Monte Carlo simulations. We implement the model on real and simulated data from Gaussian as well as other distributions, showing that the model can flexibly adapt to different data distributions. We apply the irregular autoregressive model to the residuals of a transit of an extrasolar planet to illustrate errors that remain with temporal structure. We also apply this model to residuals of an harmonic fit of light-curves from variable stars to illustrate how the model can be used to detect incorrect parameter estimation.

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

RR Lyrae stars (RRLs) are tracers of the Milky Way's fossil record, holding valuable information on its formation and early evolution. Owing to the high interstellar extinction endemic to the Galactic plane, distant RRLs lying at low Galactic latitudes have been elusive. We attained a census of 1892 high-confidence RRLs by exploiting the near-infrared photometric database of the VVV survey's disk footprint spanning $\sim$70$^\circ$ of Galactic longitude, using a machine-learned classifier. Novel data-driven methods were employed to accurately characterize their spatial distribution using sparsely sampled multi-band photometry. The RRL metallicity distribution function (MDF) was derived from their $K_s$-band light curve parameters using machine-learning methods. The MDF shows remarkable structural similarities to both the spectroscopic MDF of red clump giants and the MDF of bulge RRLs. We model the MDF with a multi-component density distribution and find that the number density of stars associated with the different model components systematically changes with both the Galactocentric radius and vertical distance from the Galactic plane, equivalent to weak metallicity gradients. Based on the consistency with results from the ARGOS survey, three MDF modes are attributed to the old disk populations, while the most metal-poor RRLs are probably halo interlopers. We propose that the dominant [Fe/H] component with a mean of $-1$ dex might correspond to the outskirts of an ancient Galactic spheroid or classical bulge component residing in the central Milky Way. The physical origins of the RRLs in this study need to be verified by kinematical information.

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.

The VISTA Variables in the V\'ia L\'actea (VVV) is a near-IR time-domain survey of the Galactic bulge and southern plane. One of the main goals of this survey is to reveal the 3D structure of the Milky Way through their variable stars. Particularly the RR Lyrae stars have been massively discovered in the inner regions of the bulge ($-8^\circ \lesssim b \lesssim -1^\circ$) by optical surveys such as OGLE and MACHO but leaving an unexplored window of more than $\sim 47$ sq deg ($-10.0^\circ \lesssim \ell \lesssim +10.7^\circ$ and $-10.3^\circ \lesssim b \lesssim -8.0^\circ$) observed by the VVV Survey. Our goal is to characterize the RR Lyrae stars in the outer bulge in terms of their periods, amplitudes, Fourier coefficients, and distances, in order to evaluate the 3D structure of the bulge in this area. The distance distribution of RR Lyrae stars will be compared to the one of red clump stars that is known to trace a X-shaped structure in order to determine if these two different stellar populations share the same Galactic distribution. We report the detection of more than 1000 RR Lyrae ab-type stars in the VVV Survey located in the outskirts of the Galactic bulge. Some of them are possibly associated with the Sagittarius Dwarf Spheroidal Galaxy. We calculated colors, reddening, extinction, and distances of the detected RR Lyrae stars in order to determine the outer bulge 3D structure. Our main result is that, at the low galactic latitudes mapped here, the RR Lyrae stars trace a centrally concentrated spheroidal distribution. This is a noticeably different spatial distribution to the one traced by red clump stars known to follow a bar and X-shape structure. We estimate the completeness of our RRab sample in $80\%$ for $K_{\rm s}\lesssim15$ mag.

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 a ground based optical transmission spectrum of the inflated sub-Jupiter mass planet WASP-6b. The spectrum was measured in twenty spectral channels from 480 nm to 860nm using a series of 91 spectra over a complete transit event. The observations were carried out using multi-object differential spectrophotometry with the IMACS spectrograph on the Baade telescope at Las Campanas Observatory. We model systematic effects on the observed light curves using principal component analysis on the comparison stars, and allow for the presence of short and long memory correlation structure in our Monte Carlo Markov Chain analysis of the transit light curves for WASP-6. The measured transmission spectrum presents a general trend of decreasing apparent planetary size with wavelength and lacks evidence for broad spectral features of Na and K predicted by clear atmosphere models. The spectrum is consistent with that expected for scattering that is more efficient in the blue, as could be caused by hazes or condensates in the atmosphere of WASP-6b. WASP-6b therefore appears to be yet another massive exoplanet with evidence for a mostly featureless transmission spectrum, underscoring the importance that hazes and condensates can have in determining the transmission spectra of exoplanets.

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.