
Elorrieta, F., Eyheramendy, S., & Palma, W. (2019). Discretetime autoregressive model for unequally spaced timeseries observations. Astron. Astrophys., 627, 11 pp.
Abstract: Most timeseries 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 continuoustime 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 discretetime process. We show that the model is weakly stationary and that it can be represented as a statespace 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.



Elorrieta, F., Eyheramendy, S., Palma, W., & Ojeda, C. (2021). A novel bivariate autoregressive model for predicting and forecasting irregularly observed time series. Mon. Not. Roy. Astron. Soc., 505(1), 1105–1116.
Abstract: In several disciplines, it is common to find time series measured at irregular observational times. In particular, in astronomy there are a large number of surveys that gather information over irregular time gaps and in more than one passband. Some examples are PanSTARRS, ZTF, and also the LSST. However, current commonly used time series models that estimate the time dependence in astronomical light curves consider the information of each band separately (e.g, CIAR, IAR, and CARMA models) disregarding the dependence that might exist between different passbands. In this paper, we propose a novel bivariate model for irregularly sampled time series, called the Bivariate Irregular Autoregressive (BIAR) model. The BIAR model assumes an autoregressive structure on each time series; it is stationary, and it allows to estimate the autocorrelation, the crosscorrelation and the contemporary correlation between two unequally spaced time series. We implemented the BIAR model on light curves, in the g and r bands, obtained from the ZTF alerts processed by the ALeRCE broker. We show that if the light curves of the two bands are highly correlated, the model has more accurate forecast and prediction using the bivariate model than a similar method that uses only univariate information. Further, the estimated parameters of the BIAR are useful to characterize longperiod variable stars and to distinguish between classes of stochastic objects, providing promising features that can be used for classification purposes.



Mancini, L., Sarkis, P., Henning, T., Bakos, G. A., Bayliss, D., Bento, J., et al. (2020). The highly inflated giant planet WASP174b. Astron. Astrophys., 633, 12 pp.
Abstract: Context. The transiting exoplanetary system WASP174 was reported to be composed by a mainsequence F star (V = 11.8 mag) and a giant planet, WASP174b (orbital period Porb = 4.23 days). However only an upper limit was placed on the planet mass (<1.3 MJup), and a highly uncertain planetary radius (0.71.7 RJup) was determined.Aims. We aim to better characterise both the star and the planet and precisely measure their orbital and physical parameters.Methods. In order to constrain the mass of the planet, we obtained new measurements of the radial velocity of the star and joined them with those from the discovery paper. Photometric data from the HATSouth survey and new multiband, highquality (precision reached up to 0.37 mmag) photometric followup observations of transit events were acquired and analysed for getting accurate photometric parameters. We fit the model to all the observations, including data from the TESS space telescope, in two different modes: incorporating the stellar isochrones into the fit, and using an empirical method to get the stellar parameters. The two modes resulted to be consistent with each other to within 2<sigma>.Results. We confirm the grazing nature of the WASP174b transits with a confidence level greater than 5 sigma, which is also corroborated by simultaneously observing the transit through four optical bands and noting how the transit depth changes due to the limbdarkening effect. We estimate that approximate to 76% of the disk of the planet actually eclipses the parent star at midtransit of its transit events. We find that WASP174b is a highlyinflated hot giant planet with a mass of Mp = 0.330 +/ 0.091 MJup and a radius of Rp = 1.435 +/ 0.050 RJup, and is therefore a good target for transmissionspectroscopy observations. With a density of rho (p) = 0.135 +/ 0.042 g cm(3), it is amongst the lowestdensity planets ever discovered with precisely measured mass and radius.



Sandford, E., Espinoza, N., Brahm, R., & Jordan, A. (2019). Estimation of singly transiting K2 planet periods with Gaia parallaxes. Mon. Not. Roy. Astron. Soc., 489(3), 3149–3161.
Abstract: When a planet is only observed to transit once, direct measurement of its period is impossible. It is possible, however, to constrain the periods of single transiters, and this is desirable as they are likely to represent the cold and far extremes of the planet population observed by any particular survey. Improving the accuracy with which the period of single transiters can be constrained is therefore critical to enhance the longperiod planet yield of surveys. Here, we combine Gaia parallaxes with stellar models and broadband photometry to estimate the stellar densities of K2 planet host stars, then use that stellar density information to model individual planet transits and infer the posterior period distribution. We show that the densities we infer are reliable by comparing with densities derived through asteroseismology, and apply our method to 27 validation planets of known (directly measured) period, treating each transit as if it were the only one, as well as to 12 true single transiters. When we treat eccentricity as a free parameter, we achieve a fractional period uncertainty over the true single transits of 94(58)(+87) per cent, and when we fix e = 0, we achieve fractional period uncertainty 15(6)(+30) per cent, a roughly threefold improvement over typical period uncertainties of previous studies.



Sedaghati, E., MacDonald, R. J., CasasayasBarris, N., Hoeijmakers, H. J., Boffin, H. M. J., Rodler, F., et al. (2021). A spectral survey of WASP19b with ESPRESSO. Mon. Not. Roy. Astron. Soc., 505(1), 435–458.
Abstract: Highresolution precision spectroscopy provides a multitude of robust techniques for probing exoplanetary atmospheres. We present multiple VLT/ESPRESSO transit observations of the hotJupiter exoplanet WASP19b with previously published but disputed atmospheric features from low resolution studies. Through spectral synthesis and modelling of the RossiterMcLaughlin (RM) effect we calculate stellar, orbital and physical parameters for the system. From narrowband spectroscopy we do not detect any of Hi, Fei, Mgi, Cai, Nai, and Ki neutral species, placing upper limits on their line contrasts. Through crosscorrelation analyses with atmospheric models, we do not detect Fei and place a 3 sigma upper limit of on its mass fraction, from injection and retrieval. We show the inability to detect the presence of H2O for known abundances, owing to lack of strong absorption bands, as well as relatively low S/N ratio. We detect a barely significant peak (3.02 +/ 0.15 sigma) in the crosscorrelation map for TiO, consistent with the subsolar abundance previously reported. This is merely a hint for the presence of TiO and does not constitute a confirmation. However, we do confirm the presence of previously observed enhanced scattering towards blue wavelengths, through chromatic RM measurements, pointing to a hazy atmosphere. We finally present a reanalysis of lowresolution transmission spectra of this exoplanet, concluding that unocculted starspots alone cannot explain previously detected features. Our reanalysis of the FORS2 spectra of WASP19b finds a similar to 100x subsolar TiO abundance, precisely constrained to , consistent with the TiO hint from ESPRESSO. We present plausible paths to reconciliation with other seemingly contradicting results.

