Research
I am primarily interested in binary interaction: how the exchange of mass, energy and angular momentum between two stars changes their evolution and final fates. My work combines targeted observations of small samples of binary interaction products with data-driven methods applied to large-scale surveys. A full list of papers can be found on my publications page.


Intermediate-mass helium stars in the Magellanic Clouds
Drout et al. (2023) and Götberg et al. (2023) discovered a population of hot, helium-rich stars with masses from 2 to 8 M⊙ in the Magellanic Clouds, which are thought to have lost their hydrogen envelopes via interaction with a binary companion. Using multi-epoch Magellan/MagE spectroscopy, I measured radial velocities for eight of these helium stars. Five are in binaries with orbital periods from 16 hours to 300 days and companions that are invisible in the optical spectra, while three show no binary motion. Based on their orbits, the optically dark companions could be low-mass stars near the main sequence, low-mass stripped stars, or compact objects. The binaries are hydrogen-poor and the apparently-single stars are nearly hydrogen-free, which suggests that the two groups formed through different evolutionary pathways. I also contributed to the Stripped-Star Ultraviolet Magellanic Cloud Survey (SUMS) led by Bethany Ludwig, which identified hundreds of stripped star candidates from their UV excess.
Laroche et al. (2026), arXiv:2609.30384; Ludwig, Drout, Götberg, Lang & Laroche (2026), ApJ 999, 73

Data-driven models for Gaia XP spectra
Data-driven models of stellar spectra are often trained to map stellar labels to spectra, which means their performance can be limited by the stellar label systematics. With Josh Speagle, I developed a variational auto-encoder that learns to generate Gaia BP/RP (XP) spectra and their intrinsic scatter without relying on stellar labels. The model reconstructs XP spectra as well as label-dependent models. When applied to giant stars with APOGEE abundances it recovers the [α/M] bimodality, robustly showing that XP spectra contain [α/M] information. Because the model needs only spectra to train, it can be applied to any and all of the ~220 million stars with XP spectra in Gaia DR3, and to the new XP spectra coming in Gaia DR4 (Dec. 2026!). I have also applied a similar autoencoder framework to identify unresolved binaries in photometric surveys with Tobias Géron, beginning with Rubin observations of 47 Tucanae.
Laroche & Speagle (2025), ApJ 979, 5; Laroche & Speagle (2023), ICML ML4Astro workshop; Géron, Laroche et al. (2026), arXiv:2609.05616

Ultra-light dark matter and strong lensing
With Jo Bovy and Daniel Gilman, I used the flux ratios of eleven quadruply-imaged quasars to constrain ultra-light dark matter (ULDM), including the wave-interference fluctuations ULDM should produce in the host halo density profile, calibrated against simulations. These quantum fluctuations can masquerade as dark matter halos, perturbing the flux ratios in much the same way, and including them substantially changes the inferred particle mass. Even so, the data disfavor particle masses below 10−21.5 eV.
Laroche et al. (2022), MNRAS 517, 1867
Talks
Radial velocity survey of stars stripped in binaries in the Magellanic Clouds, TASTY talk, Department of Astronomy & Astrophysics, University of Toronto.
