An Eclipsing 47 minute Double White Dwarf Binary at 400 pc

We present the discovery of the eclipsing double white dwarf (WD) binary WDJ 022558.21-692025.38 that has an orbital period of 47.19 min. Following identification with the Transiting Exoplanet Survey Satellite, we obtained time-series ground based spectroscopy and high-speed multi-band ULTRACAM photometry which indicate a primary DA WD of mass 0.40 +- 0.04 Msol and a 0.28 +- 0.02 Msol mass secondary WD, which is likely of type DA as well. The system becomes the third-closest eclipsing double WD binary discovered with a distance of approximately 400 pc and will be a detectable source for upcoming gravitational wave detectors in the mHz frequency range. Its orbital decay will be measurable photometrically within 10 yrs to a precision of better than 1%. The fate of the binary is to merge in approximately 41 Myr, likely forming a single, more massive WD.


INTRODUCTION
Compact white dwarf (WD) binaries are of strong astrophysical interest as descendants of mass transfer phases, possibly the progenitors of exotic merger stars (Webbink 1984;Zhang & Jeffery 2012;Cheng et al. 2019;Gvaramadze et al. 2019;Hollands et al. 2020) and are one of the likely progenitors of type Ia/.Ia supernovae (Maoz et al. 2014). Among double WD (DWD) binary systems, the subset which exit mass transfer with periods less than ∼ 10 hours are destined to return into contact within a Hubble time, primarily as a consequence of gravitational wave radiation. There are hence many DWD binaries with orbital frequencies in the mHz frequency band, allowing for detection with upcoming space-based gravitational wave detectors ★ Email: james.munday98@gmail.com like TianQuin (Luo et al. 2016) and the Laser Interferometer Space Antenna (LISA, Amaro-Seoane et al. 2017, with DWDs acting as dominant foreground sources (Nissanke et al. 2012;Korol et al. 2017;Lamberts et al. 2019).
Although the Galactic population of DWD binaries is expected to be in the hundreds of millions (e.g. Nelemans et al. 2001;Marsh 2011;Korol et al. 2022), only a small fraction of the observable population has been conclusively identified to date. This is largely the result of observational biases from the intrinsic dimness of WDs and from a detection preference with eclipsing DWD systems, which reveal the clear presence of a companion but require near edge-on inclinations. Ongoing and recent attempts to exploit all-sky variability surveys to find DWD and other compact binaries has resulted in an acceleration in the detection of such systems (e.g. Burdge Ren et al. 2023). Furthermore, the extremely low mass (ELM) WD survey (Brown et al. 2022) and the supernovae type Ia progenitor survey (Napiwotzki et al. 2020) have provided population statistics on the DWD sample from the DWDs discovered in their surveys. In particular, Napiwotzki et al. (2020) reveal a DWD fraction of 6% compared to their observed WD sample, which is consistent with the fraction expected in from population synthesis . A couple of dozen other radial-velocity variable or eclipsing DWD systems have been discovered separately to these surveys (see e.g. Kilic et al. 2010;Hallakoun et al. 2016;Kilic et al. 2021a, and references therein) which typically have orbital periods between a few hours and a day. However, the most compact of DWD binaries with an orbital period less than one hour are rare. Brown et al. (2022) find that ≈ 10% of low-mass WD binaries alone have sub-hour orbital periods in the ELM WD survey. Approximately 20 out of a total of 200 DWDs have been found to be this compact, an excessive fraction compared to the ≈ 5% predicted (Nelemans et al. 2001) that is likely consistent when considering observational biases. Eclipsing DWD binaries offer the most precise system characterisation and those with sub-hour orbital periods have the potential for their orbital decay to be measured observationally within decades. In this paper, we present the discovery of such a compact and eclipsing DWD binary with orbital period of 47.19 min. Its binarity was discovered using the Transiting Exoplanet Survey Satellite (TESS, Ricker et al. 2015), making it one of the shortest-period detached binary discovered by TESS to date. Compared to the characterised eclipsing DWD binary sample, WDJ 022558.21−692025.38 (hereafter J0225−6920) is located relatively nearby in the Galaxy at a distance of approximately 400 pc; at the time of writing, J0225−6920 becomes the third-closest eclipsing DWD discovered after NLTT 11748 (Steinfadt et al. 2010) and GALEX J1717+6757 (Vennes et al. 2011), both being of distance ≈ 180 pc.
Our observations are discussed in §2 followed by spectral modelling in §3. Then, we address light curve modelling and derived system parameters in §4. We close by discussing the projected orbital decay of J0225−6920 in §5.

Photometry
J0225−6920 was identified as a DA white dwarf by Kilkenny et al. (2015, EC02251-6933) and classified by Gentile Fusillo et al. (2019Fusillo et al. ( , 2021) as a high-probability white dwarf ( WD > 0.95). J0225−6920 was part of a 2 min cadence program with TESS in Sectors 27, 28, and 29 (program G03124, TIC 631238061). During this program, we saw clear eclipsing signatures in J0225−6920 that showed the system to be a DWD binary, although its orbital period was not immediately clear from boxed-least-squares periodogram aliases of 23.5 min and 47 min. To distinguish the orbital period, we initially obtained (white light) data from the 0.41 m PROMPT-1 telescope (Reichart et al. 2005), data from the 1.6 m Pico dos Dias Observatory (BG40 filter) and data from the South African Astronomical Observatory (SAAO) 1.0 m Lesedi telescope (white light). These revealed distinct primary and secondary eclipses for a 47 min orbital period.
Following identification, we then observed J0225−6920 with the high-speed photometer ULTRACAM (Dhillon et al. 2007) mounted on the ESO 3.5 m New Technology Telescope (NTT) during five nights in July 2021. Observations were simultaneously taken using the Super SDSS , and filters (Dhillon et al. 2021). We later observed with ULTRACAM on the NTT over three nights in September 2022 in the Super SDSS , and filters. In total, we observed more than 15 hrs on target with ULTRACAM. A full observing log is supplied in Appendix A.
All data were debiased and flat-fielded; the SAAO and PROMPT-1 data with custom scripts, and the ULTRACAM data using the HiPER-CAM reduction pipeline (Dhillon et al. 2021). Dark-current subtraction was additionally performed for the ULTRACAM data because of its hotter operational temperature. The flux of J0225−6920 was extracted through differential aperture photometry using a non-variable comparison star with Gaia DR3 source ID 4693540643322249216. We used a variable aperture size that reflected the seeing at the time of observation set to 1.8× the full-width at half-maximum of the stellar profile. All mid-exposure time-stamps were placed on a Barycentric Modified Julian Date (BMJD) Barycentric Dynamical Time (TDB) time frame.

Spectroscopy
We obtained time-series long-slit spectra from the 4.1 m SOAR telescope with the SOAR Goodman spectrograph (Clemens et al. 2004) in October 2020. We observed with a 930 1 mm −1 grating and a slit-width of 1 ′′ , resulting in a resolution of /Δ ≈ 1800 with a wavelength coverage of 3600 -5275 Å. 43 consecutive spectra were obtained with exposure time 180 s and a 7 s readout time between adjacent exposures. FeAr arc lamp spectra were taken immediately before the first and immediately after the final spectrum. To avoid any possible drift to the wavelength solution extracted from FeAr arc lamps, we placed a second star (Gaia DR3 4693540643322249216) on the slit and calibrated the wavelength solution of each science exposure using the Balmer series of this comparison, i.e. the initial wavelength solution was obtained in the reference frame of the comparison star. The systemic velocity of the binary was measured from each set of calibrated spectra (see §3.2) and the systemic velocity found through the FeAr lamp wavelength solution was used to convert from the reference frame of the comparison star to that of the barycentre, retaining the absolute wavelength calibration.
The signal-to-noise ratio (SNR) in each spectrum at the centre and wings of H is approximately 15 and 23, respectively. These reveal clear Balmer absorption features indicating that the brighter component is a DA WD (see Fig. 1). Double-lined features in the spectrum were searched for, particularly at H , but were not obviously apparent. No metals nor helium spectral lines were detected in the spectra. The observations are displayed in black, with the best-fit model overplotted in red. Each profile is vertically offset for clarity. The displayed model has a primary star of eff,1 = 25 330 K and log 1 = 6.99 dex combined with a secondary star of eff,2 = 13 750 K and log 2 = 7.60 dex. The secondary star's effective temperature that is included here was derived from a single iteration of light curve fitting using the eff,1 of a single-star fit to the co-added spectrum.

Atmospheric Parameters
The reduced spectra were co-added with a common rest wavelength. The normalised Balmer line profiles of the co-added spectrum were then fit with the DA WD models of Tremblay et al. (2011Tremblay et al. ( , 2015. Fitting the spectrum with a single DA WD model to represent the primary star 1 , we constrain the atmosphere of the brighter component 2 to have a surface gravity log 1 = 7.07 ± 0.04 dex and eff,1 = 24 250 ± 310 K. In an attempt to retrieve the surface gravity and temperature of the secondary star, and improve the accuracy of the primary star parameters, we utilised photometric flux measurements from widefield surveys and fit solutions for two unique WDs to the combined SED. At the time of writing, Gaia DR3 and SkyMAPPER (Wolf et al. 2018) survey data are available, which we utilise to model the photometric and spectroscopic data simultaneously. We dereddened these flux measurements with an extinction constant of = 0.09 (Gentile Fusillo et al. 2021) following the reddening prescription of Fitzpatrick & Massa (2007) with = 3.1. The WD models require a conversion from an Eddington flux to a flux observed at Earth, so we enforced a posterior distribution on the fitted distance based on the Gaia DR3 parallax stated in Table 1. A single DA WD model to 1 We follow the convention in this paper that the primary star is the hotter and brighter WD. 2 We present errors on the atmospheric constraints of the stars by combining in quadrature the statistical error from 2 fitting with an external error of 1.2% for eff,1 and 0.038 dex for log 1 (Liebert et al. 2005). This is primarily to account for the error induced from flux calibration of spectra. the dereddened Gaia photometric fluxes yields eff = 22 500 ± 910 K (Gentile Fusillo et al. 2021), which places a minimum temperature on the primary star since the flux contributed from the cooler companion weights a single star fit to cooler temperatures.
We then simultaneously fit all dereddened photometric flux measurements and the normalised Balmer profiles with a DA+DA model and enforced eff,2 = 13 750 K. This temperature measurement was guided by the results of our light curve fitting in §4 when using the measured eff,1 from a single star model to the co-added spectrum. The goals from fitting a DA+DA model were to improve the accuracy of the primary temperature by including the flux from a WD companion and to gain a measurement of the secondary's surface gravity. We obtain results of eff,1 = 25 330 ± 330 K, with log 1 = 6.99 ± 0.04 dex and log 2 = 7.60 ± 0.23 dex (see Fig. 1). Interpolation of these parameters with WD atmosphere models indicate a primary mass of 1 = 0.35 ± 0.01 M ⊙ for a helium core (Althaus et al. 2013) and 1 = 0.27 ± 0.01 M ⊙ for a carbon-oxygen core ) WD. A hybrid carbon-oxygen/helium core composition is also possible, and predicted masses for this composition are similar to pure-helium masses at the primary temperature (Zenati et al. 2019). Although both a carbon-oxygen and helium core WD would be possible core compositions for the primary depending on the evolutionary path of the system, carbon-oxygen-core WDs are not expected below roughly 0.33 M ⊙ (Prada Moroni & Straniero 2009), so we only consider solutions with a helium-core primary going forward. The core composition of the secondary is unclear from spectral analysis alone due to its large surface gravity error.
Lastly for atmospheric analysis, a DA+DBA WD spectral fit was performed with eff,2 again fixed at 13 750 K and we found two inconsistencies. This model significantly under-predicted the flux of the weakest Balmer transitions and we found that absorption of He I (4471 Å) would be apparent, which it is not in the observed spectrum. We consequently predict that the companion has a hydrogen rich atmosphere, perhaps of type DA also. This DA+DA WD prediction leads us to assume the spectroscopic results from the DA+DA model fit for the remainder of the paper, as listed in Table 2.

Radial Velocities
We then searched for any radial velocity variations in the normalised H absorption-line profile by modelling the full set of spectra simultaneously. The velocity of each star varies as a function of phase given by 1,2 ( ) = ± 1,2 sin 2 (1) where subscripts 1 and 2 represent the primary and the secondary stars, is the systemic velocity, is the radial velocity semiamplitude and is the orbital phase of the binary. The orbit is assumed circularised (e.g. Peters 1964) and is solved with the centre time of the exposure and our orbital ephemeris (equation 2). Furthermore, the line profiles themselves were fit with a 2-Gaussian component similar to all spectra 3 , each having a standard deviation, We measure 1 = 224.0 ± 4.4 km s −1 and = 1.6 ± 4.1 km s −1 . For clarity, we extract radial velocity measurements independently of orbital phase (as listed in Appendix D) by keeping , and , fixed from a best-fit to the stacked spectrum. Radial velocity errors are reported as the standard deviation of 1000 bootstrapping iterations. These measurements and our fit 1 and are plotted in Fig. 2. It is important to note, however, that the measured radial velocity amplitude may be underestimated due to the contribution of light from the companion in anti-phase. Even though absorption lines from the companion are not obviously apparently (see Fig. 2), its contribution can be significant (e.g. Hallakoun et al. 2016). To check if this has a significant impact in J0225−6920, we searched for a consistent radial velocity amplitude across all Balmer lines by measuring those of H − H , finding them to be within 2 of each other; H being the largest. Light curve modelling to the photometry alone agreed best with the radial velocity at H (see §4.2). As such, we report the radial velocity amplitude of H as that of the system.

Ephemeris
To determine the orbital ephemeris, we simultaneously searched for an orbital period using the photometry from all detectors and filters using the multi-band Lomb-Scargle periodogram package of Vander-Plas & Ivezić (2015). By fitting individual eclipses with our best-fit synthetic light curves for the ULTRACAM data (see §4.2), we refined the ephemeris to a higher precision. The orbital ephemeris that we obtain, centred on the primary mid-eclipse, is BMJD TDB = 59403.299199(90) + 0.03277099777 (24) (2)  47 min DWD Binary v stars is assumed. The same Galactic reddening treatment as in §3.1 was applied, which has a small effect on the light curve morphology but is most impactful to the eclipse depths (e.g. Jones et al. 2020). The synthetic light curves do not include the contribution of gravitational lensing, however, its impact on the emergent flux is suspected to be negligible for this system, amplifying the received flux by far less than 1% (Marsh 2001). Limb darkening and gravity darkening coefficients were interpolated uniquely for the two stars using tabulated values that are specific for the Super SDSS passbands. For limb darkening, we used a power-law prescription with its coefficients set according to the 3dimensional WD models of Claret et al. (2020b), or the 1-dimensional WD models of Claret et al. (2020a) when the 3-dimensional grid boundaries were exceeded. Gravity darkening coefficients were interpolated from the tables of Claret et al. (2020a). In each case, the interpolation was dependent on the effective temperature and the surface gravity of the respective star for a trial synthetic model.
We also reactivated the Doppler beaming functionality of PHOEBE 2.2 and pass beaming coefficients of Claret et al. (2020b) 5 . The synthetic 1 and 2 are solved dynamically in PHOEBE as a function of phase (which is to say through Keplerian geometry) and incorporate the impact of gravitational redshift. The effect of smearing from a finite exposure time was corrected for assuming exposures of 9 s in the -band and 3 s in the , and bands (see Table A).
To obtain a system solution, we implemented an MCMC algorithm using EMCEE (Foreman-Mackey et al. 2013), where the goodness of fit was determined by minimisation of the 2 between the synthetic and observed light curves. We fit synthetic photometry to the UL-TRACAM light curves simultaneously to improve the precision of the radii, masses and the inclination of the system. Since blackbody models are considered, the spectral energy distribution assumed for both stars in PHOEBE differs from that of a WD. The usage of blackbody models causes little issue to the , and bands as the colour between each band for WD and blackbody spectra is small, whereas the − colour deviates strongly. To overcome these issues, we model to the , and bands only and we set a Gaussian prior for the primary star temperature based on our spectroscopic solution. With a final model, we then check for consistency in the band by allowing the temperature of the secondary to be a free parameter. The secondary temperature is identical for each of the other passbands.
We started by modelling independently of 1 and found a good fit to the photometry where the synthetic 1 was within the error of that measured, giving us confidence that the measured 1 in §3.2 is reflective of the true value. Following this test, we introduced 1 into the 2 minimisation. Overall, the mass, temperature and radius of the primary and secondary star and the inclination of the system were free parameters in the MCMC. Only eff,1 included a Gaussian prior while all other parameters had no set prior.

Binary parameters
Our post-burnin results of the MCMC are presented in Table 3, with our best-fitting binary model displayed in Fig. 3. A corner plot diagram is included in Appendix B showing the covariance between ( §3.1), the spectral energy distribution for WD atmospheres would not be exact neither, making a light-curve extracted eff,2 an approximation. 5 In PHOEBE 2, this feature was disabled due to concerns of the accuracy of beaming factor computation from a non-implicit and fitted trend of the SED. In the Claret et al. (2020b) models, the beaming factor is calculated implicitly and integrated over the full SED. Therefore, the passed beaming constant is an accurate representation. free parameters. A good model fit to the light curve is obtained about the eclipses and for 1 . Furthermore, the impact of Doppler beaming in quadrature (which causes the observed flux to be larger/smaller as a star moves towards/away from the observer) is well fit in the -band; the band impacted most significantly (see e.g. Claret et al. 2020b).
From the light curve modelling, the primary star has an effective temperature of 25 550 ± 200 K and for the secondary 14 350 ± 100 K. Component masses and radii indicate that both stars likely have helium cores, falling closely on the evolutionary tracks of helium-core WD models (Althaus et al. 2013;Istrate et al. 2016), and that the system is detached. The primary star's Roche lobe is 12% filled and the secondary's is 18% filled. Interpolation of the masses and the temperatures of each WD with Istrate et al. (2016) cooling track models that feature rotational and diffusional physics indicate cooling ages for each star of 1 ≈ 7 Myr and 2 ≈ 22 Myr, such that the more massive primary star is the younger of the two. The primary mass determined from light curve fitting is consistent with the spectroscopically predicted mass of 1 = 0.35 ± 0.01 M ⊙ for a WD with a helium core (see §3.1). While the log 2 = 7.60 ± 0.23 dex for a two-star spectral fit shows a worse agreement to the light-curve inferred parameters, the measurements agree at approximately a 2 level, where such a difference could have been influenced by the systematical errors such as smearing of the secondary star to the co-added spectrum. The modelled parameters indicate that the secondary is also likely a helium core WD.

Similarities with known systems
J0225−6920 bears a resemblance to the systems ZTF J1749+0924, ZTF J2029+1534 and ZTF J0722-1839 ). These systems are also detached DWD binaries with similar effective temperatures and masses to J0225−6920, each having a shorter orbital period. All of these systems likely exited a common-envelope phase under similar conditions and follow a similar evolutionary history. In each binary, both stars are likely helium-core WDs. The characterisation of these binaries has also been used in the study of Scherbak & Fuller (2023) to investigate the common envelope efficiency parameter, , from a second mass transfer phase. Given the similarity of J0225−6920, our measurements are consistent with ≈ 0.2 − 0.4, as inferred by these authors.
A further two compact systems that are comparable to J0225−6920, having similar star masses and temperatures, are SDSS J1152+0248 (Hallakoun et al. 2016;Parsons et al. 2020) and CSS 41177 (Bours et al. 2014(Bours et al. , 2015. These eclipsing DWD binaries have an orbital period of 144 min, and 167 min, respectively. J0225−6920 appears to bridge the period gap between all of these DWD binaries, which can be utilised in the future to study the impact of tides on orbital decay (e.g. Piro 2019).
Lastly, J0225−6920 is much alike SDSS J232230.20+050942.06 in regards to the core composition, being the first double helium-core, DWD, LISA verification binary discovered (Brown et al. 2020). Systems such as these are vital gravitational wave sources for data quality verification of LISA when launched, given that the gravitational wave frequency and amplitude are solvable from the orbital parameters of the binary alone.

J0225−6920 AS A GRAVITATIONAL WAVE SOURCE
The orbit of J0225−6920 will decay due to a loss of orbital angular momentum from gravitational wave radiation and tides. If one as-  (2) and binned for clarity. PHOEBE synthetic light curves are overlaid in black. The AB magnitude is displayed on the vertical axes where all but the -band are vertically offset by = 0.23 mag, = 0.14 mag and = −0.16 mag for clarity. Top right/Bottom right: similar to the left-hand plot, but zoomed around the primary and secondary eclipses. Both right-hand-side plots are zoomed insets, where the light curves from each filter are manually offset in magnitude to enlarge the eclipses.
The orbital decay can be precisely measured with continued observations using the time of arrival of a given orbital phase, where a number of WD binaries have been characterised in this way through continued observations (Barros et al. 2007;Hermes et al. 2012;de Miguel et al. 2018;Burdge et al. 2019aBurdge et al. ,b, 2020bMunday et al. 2023). The deep eclipses of J0225−6920 make mid-eclipse timing the natural way to do so, where an earlier eclipse arrival time is explainable by an orbital decay. The deviation of the eclipse timing from a constant orbital period is dominated (to first order) by 2 /2 for general relativistic effects, with the length of time following the first eclipse. From the data presented in this paper, the arrival time difference from the full span of data is currently a mere second, however, the bright Gaia magnitude ( = 16.4 mag) and its 0.2 mag deep primary eclipses makes J0225−6920 an exemplar system for a measured that can be utilised to probe non-gravitational wave induced orbital decay in the future. This may arise in J0225−6920 due to tidal dissipation (Benacquista 2011;Piro 2011;Fuller & Lai 2012;Piro 2019), in which orbital energy transferred into rotational energy of the stars is the cause of a faster orbital decay. Additionally, measurement of the orbital decay constrains the chirp mass, thus restricting the combinations of primary and secondary star masses for system modelling (e.g. Munday et al. 2023).
We find that timing of individual J0225−6920 primary and secondary eclipses in ideal conditions with ULTRACAM is accurate to 0.5 s and 2 s, respectively. Considering this and the predicted orbital decay, we expect to obtain a period derivative precise to better than ≈5% after 5 yrs or ≈1% after 10 yrs; the precision scaling with 2 . The observed deviation of the centre of eclipse arrival time will be ≈4 s and ≈15 s after 5 yrs and 10 yrs. Furthermore, the compact binary is ideally situated in the LISA frequency band, as shown in Fig. 4  distance of 402.6 +5.9 −6.8 pc and derived orbital parameters generates a characteristic strain of (4.5 ± 0.7) × 10 −20 , with a projected SNR of 1.3±0.2 after a 4 yr LISA mission time, or 2.1±0.4 after 10 yrs (Wagg et al. 2022a,b). J0225−6920 is thus a detectable verification binary for LISA and will be the closest double helium core DWD verification binary known to date; which is a class that is expected to comprise 31% of all LISA detectable binaries (Lamberts et al. 2019). For comparison, the second-and third-closest double helium-core DWD verification binaries are SDSS J063449.92+380352.2 at 435 pc and SMSS J033816.16−813929.9 at 533 pc (Kilic et al. 2021b).
In the future, J0225−6920 is expected to merge into a single WD, likely undergoing a hot subdwarf phase during which helium is burnt to form a carbon-oxygen core (Han et al. 2002). For a circularised binary star system under purely general relativistic orbital decay, two stars will merge after a time (Peters 1964) with 0 the present-day semi-major axis, the gravitational constant and the speed of light. We hence predict J0225−6920 to merge within (0.379 ⊙ ) = 41 Myr.

CONCLUSIONS
We have discovered that J0225−6920 is an eclipsing DWD binary with an orbital period of 47.19 min. Multi-band light curve modelling indicates that the system likely consists of two helium-core WDs, where the primary has a DA spectral type, having a pure hydrogen surface composition, and the secondary likely does too. Its relatively 5 The data used to create Fig. 4 was obtained through https://gitlab.in2p3.fr/LISA/lisa-verification-binaries.
close distance and brightness makes the binary a prime candidate to measure the orbital period decay in the future, key to investigate deviations from purely general relativistic orbital decay due to tides. The binary will be a detectable source for the LISA spacecraft, and we believe it is the first such DWD verification binary discovered by the TESS mission. J0225−6920 joins a small sample of fully characterised DWD binaries which will merge within a Hubble time.
It is the closest double-helium-core DWD verification binary known to date; a class that is expected to account for 31% of all LISA detectable DWDs (Lamberts et al. 2019). Table A1. An observing log of all J0225−6920 observations acquired with ULTRACAM which were subsequently used for light-curve modelling. Nights were not noticeably impacted by clouds. The duration represents the time that the telescope was on target after acquisition. MJD mid is the MJD at the centre of the observing period. In the instrument column, UCAM is an abbreviation for ULTRACAM. Observations in the filter were 3× the cadence stated. ULTRACAM has a dead time of 24 ms between adjacent exposures.   Table D1. The extracted radial velocity measurements of J0225−6920 presented in Fig. 2. The time-stamps presented have been corrected to a barycentric reference frame (MBJD, TDB) and the velocities (V b ) have been corrected for a barycentric velocity also.