COBALT-BLUE: a fast and accurate model for gravitational self-lensing binary systems

We present COBALT-BLUE, a fast public model for the optical light curve of a luminous star in a non-accreting binary with a compact object. The model includes gravitational self-lensing, Doppler beaming, ellipsoidal variations, eccentric orbital geometry, and wavelength-dependent limb darkening. In sufficiently edge-on systems, it predicts sharp periodic self-lensing flares together with phase-locked orbital modulation outside the flare. We improve upon earlier models by including non-linear limb darkening laws based on stellar atmosphere simulations to calculate the self-lensing flare profile. This enables light curves to be self-consistently calculated in multiple wavebands from one set of physical parameters, thus offering the prospect of constraining binary parameters such as compact object mass from photometry alone if multiple self-lensing flares are observed. We also include the option to utilise simpler limb darkening prescriptions to save computational expense. As a proof of principle, we fit our model to ZTF data of 14 candidate self-lensing sources. The model describes the flare morphology well for most candidates, but we find that the best fitting parameters are not physically plausible. We conclude that the observed flares are most likely stellar flares and not self-lensing events. In future, we plan to use the model to systematically search optical survey data for signatures of binary motion.