Multivariate time series transformer embeddings for light curves of periodic variable stars

Context . Astronomical surveys produce time series data by observing stellar objects across multiple photometry bands. Foundational transformer-based models, such as Astromer , encode each time series as a sequence of embeddings to classify sources. However, such models operate independently on each band and therefore do not use information across wavelengths or filters. Aims . We extend the single-band Astromer framework by introducing a fusion layer that combines single photometry band observations into a unified sequence representation, thus enabling multiband analysis for downstream tasks, namely periodic variable star classification. The challenge in adapting the encoder for multiband data lies in coordinating information across bands observed at asynchronous times. Methods . We pre-trained various multiband neural network models on 600 000 high signal-to-noise light curves from the Massive Compact Halo Object ( MACHO ) survey and fine-tuned them using labeled data from the Alcock catalog (derived from MACHO ) and from the Asteroid Terrestrial-impact Last Alert System survey. Results . Our results show that both proposed multiband architectures outperform the single-band models by approximately ten percentage points in the F1-score. Jointly pre-trained multiband encoders further improve the performance compared to a collection of independently pre-trained single-band encoders, reaching an improvement of 23 percentage points in the F1-score, even when data are scarce. Furthermore, we find minimal differences in multiband performance between sampling bands asynchronously and sampling all bands on the same time steps, as most models remain within a three percentage point difference in F1-score on the full dataset. Conclusions . These results demonstrate a trade-off in the training speed and classification accuracy between single-band and multi-band encoders, with multiband models improving the F1-score performance by ten percentage points at the cost of increased pre-training time. Our results show the great potential of transformer-based multiband neural network architectures for classification tasks regarding future large-scale time-domain surveys with a greater variety of variable stars.