Abstract
Some accreting binary systems containing a white dwarf (such as classical novae or persistent supersoft sources) are seen to emit low-energy X-rays with temperatures of ∼ 106 K and luminosities exceeding 1035 erg s−1. These X-rays are thought to originate from nuclear burning on the white dwarf surface, either caused by a thermonuclear runaway (classical novae) or a high mass-accretion rate that sustains steady nuclear burning (persistent sources). The discovery of transient supersoft X-rays from ASASSN-16oh challenged these ideas, as no clear signatures of mass ejection indicative of a classical nova eruption were detected, and the origin of these X-rays remains controversial. It was unclear whether this star was one of a kind or representative of a larger, as yet undiscovered, group. Here, we present the discovery of 29 stars located in the direction of the Magellanic Clouds exhibiting long-duration, symmetrical optical outbursts similar to that seen in ASASSN-16oh. We observed one of these objects during an optical outburst and found it to be emitting transient supersoft X-rays, while no signatures of mass ejection (indicative of a classical nova eruption) were detected. We therefore propose that these objects form a homogeneous group of transient supersoft X-ray sources, which we dub "millinovae" because their optical luminosities are approximately a 1000 times fainter than those of ordinary classical novae.
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
A classical nova eruption occurs on the surface of a white dwarf that is accreting matter from its nondegenerate companion via an accretion disk. A layer of fresh, usually hydrogen-rich gas accumulates on the white dwarf surface and becomes degenerate. The gas is heated, and when hydrogen ultimately ignites, a thermonuclear runaway ejects about 10−5–10−4 M⊙ of matter from the system with velocities from hundreds to thousands of km s−1 (M. F. Bode & A. Evans 2008). As the ejecta expands and becomes optically thin, the exposed white dwarf surface is observed to emit supersoft X-rays due to stable hydrogen burning of the leftover gas (J. P. Osborne 2015; K. L. Page & A. W. Shaw 2022). The maximum X-ray luminosities of classical novae are typically 1036–1038 erg s−1 (e.g., M. Henze et al. 2014). In the optical passbands, they normally reach absolute magnitudes from −6 to −10 (M. F. Bode & A. Evans 2008).
The accretion disk surrounding the white dwarf is subject to a thermal instability that can lead to dwarf nova outbursts. They occur when gas in the disk reaches a critical temperature, and its viscosity abruptly changes, increasing the mass accretion rate onto the white dwarf (J.-P. Lasota 2001). If the instability develops in the outer regions of the accretion disk (outside-in outbursts), the rise to maximum brightness is rapid. However, the light curves of inside-out outbursts (which start in the inner disk regions) are more symmetric, with a relatively slow rise to the peak (J.-P. Lasota 2001). Dwarf novae in outburst reach optical absolute magnitudes from +6 to 0, depending on their orbital period (longer orbital periods mean larger accretion disks, and hence, they become brighter; J. Patterson 2011). The X-ray luminosities of dwarf novae in outbursts are typically 1029–1032 erg s−1, about 6 orders of magnitude lower than those of classical novae (A. D. Schwope et al. 2024; A. C. Rodriguez et al. 2024).
Therefore, the discovery of supersoft X-ray emission from the optical transient ASASSN-16oh (T. J. Maccarone et al. 2019) was a surprise. The X-ray emission was consistent with that of a ∼ 900,000 K blackbody and a luminosity of about 6.7 × 1036 erg s−1, similar to that seen from persistent supersoft X-ray sources and classical novae. However, the optical properties of the transient (light-curve shape, its amplitude, peak absolute magnitude, and optical spectra) did not match those of classical novae at all.
ASASSN-16oh was first detected in 2016 December by the All-Sky Automated Survey for Supernovae (ASASSN; B. J. Shappee et al. 2014) as a V = 16.9 transient in the Small Magellanic Cloud (SMC). Its archival light curve (upper panel of Figure 1) collected by the Optical Gravitational Lensing Experiment (OGLE; A. Udalski et al. 2015) revealed that the optical outburst commenced at least 4 months before the ASASSN discovery. The object exhibited irregular variability before (and after) the outburst, with a mean brightness of about I = 20.3 and V = 21.1. It peaked at I ≈ 16.6, which corresponds to MI = −2.5 ± 0.1, assuming the SMC distance modulus of μ = 18.977 ± 0.016 (D. Graczyk et al. 2020) and reddening of (D. M. Skowron et al. 2021). The full optical light curve does not match those of classical novae, which normally rise to maximum within hours to a few days (not months). Additionally, the optical spectra revealed narrow (FWHM = 164 km s−1) Balmer and He ii emission lines (T. J. Maccarone et al. 2019), unlike the broad features (several hundred to a few thousand km s−1) usually observed in classical novae.
The unusual outburst light curve, narrow optical emission lines, and X-ray spectrum led T. J. Maccarone et al. (2019) to propose that the X-rays come from a spreading layer—a belt around the white dwarf's equator near the inner disk edge in which a large fraction of the total accretion energy is emitted. However, this interpretation remains controversial, with alternative models being based on a nonejecting nova outburst (Y. Hillman et al. 2019; M. Kato et al. 2020).
Y. Hillman et al. (2019) proposed that both the optical flux and X-ray emission in ASASSN-16oh originated from a hot white dwarf that underwent a non-mass-ejecting thermonuclear runaway. They argued that such nonejecting nova events are possible if the mass accretion rate onto the white dwarf is sufficiently high (at least a few times 10−7 M⊙ yr−1) and claimed that the light curve of ASASSN-16oh is best modeled by a 1.1 M⊙ white dwarf accreting at 3.5–5 × 10−7 M⊙ yr−1. Nonetheless, their model predicted an asymmetric outburst (with a rapid rise to the peak followed by a slower decline), in stark contrast to observations. Furthermore M. Kato et al. (2020) later criticized this model, pointing out that Hillman et al. had overestimated the optical brightness by 6 mag.
On the other hand, M. Kato et al. (2020) proposed that the outburst of ASASSN-16oh was triggered by a sudden mass accretion caused by disk instability. In this model, the optical flux originates from the irradiated accretion disk, and X-rays are coming from the hot hydrogen-burning white dwarf surface. They found that the X-ray light curve of ASASSN-16oh is best explained by a massive (1.32 M⊙) white dwarf model.
Given the unusual and unexplained optical and X-ray properties of ASASSN-16oh, we decided to investigate whether this star is one of a kind or representative of a larger, as yet undiscovered, group. We thus searched for objects with outburst properties similar to those seen in ASASSN-16oh among the 20 yr-long light curves of 76 million stars observed toward the Magellanic Clouds by the OGLE survey (A. Udalski et al. 2015).
2. OGLE Photometric Data
The OGLE photometric data have been collected using the 1.3 m Warsaw Telescope located at Las Campanas Observatory, Chile, during the OGLE-III (2001–2009; A. Udalski 2003) and OGLE-IV (2010–2024; A. Udalski et al. 2015) phases of the project. Different instruments were used in different phases of the survey. The OGLE-III camera consisted of eight 2048 × 4096 CCD detectors, with a plate scale of 0.26 arcsec pixel−1. The OGLE-IV camera had the same pixel scale but comprised thirty-two 2048 × 4102 detectors, providing a field of view of 1.4 deg2. The typical exposure time was 150 s. The photometric data were extracted using a custom version of the difference image analysis algorithm (A. B. Tomaney & A. P. S. Crotts 1996; C. Alard & R. H. Lupton 1998) implemented by P. R. Woźniak (2000). The OGLE-III and OGLE-IV data have been reduced using the same reference image (P. Mróz et al. 2024), enabling us to obtain homogeneous, long-term light curves of all analyzed objects. Photometry was calibrated to the standard I- and V-band systems. See A. Udalski (2003) and A. Udalski et al. (2015) for a detailed description of the observing setup and data calibration. All data presented in this Letter are available to the astronomical community from https://ftp.astrouw.edu.pl/ogle/ogle4/millinovae/.
3. Results
We started with a list of 72,303 objects (54,996 in the Large Magellanic Cloud, LMC, and 17,307 in the SMC) whose light curves contain at least five consecutive data points magnified with respect to the remaining light curve (see P. Mróz et al. 2024 for more details about the selection of outbursting events). We required that the light curve contain at least one outburst with an amplitude larger than 1 mag and duration between 10 and 600 days. Subsequently, we cross-matched selected objects with the Gaia Data Release 3 (DR3) archive (Gaia Collaboration et al. 2023) and removed artifacts due to high-proper-motion stars (those with proper motions > 10 mas yr−1), which were the largest source of contamination. The light curves of the remaining objects (7412 in the LMC and 2689 in the SMC) were visually inspected, and we selected the initial sample of objects with triangle-shaped, symmetrical outbursts similar to those of ASASSN-16oh.
We vetted all selected objects using additional data, including multicolor photometry, Gaia parallaxes and proper motions, and sky images. We removed objects that could be classified as classical novae, supernovae, and active galactic nuclei, leaving us with a list of 29 objects (22 in the LMC and seven in the SMC). Example light curves are presented in Figure 1 (see also Figures A1 and A2 in Appendix A). All objects show long-duration (weeks to months), symmetrical (triangle-shaped) outbursts with amplitudes in I ranging from 1.0 to 3.7 mag. Some outbursts are deceivingly symmetric, which led to them being classified as gravitational microlensing events (OGLE-mNOVA-12 = MACHO-LMC-7 and OGLE-mNOVA-13 = MACHO-LMC-23, C. Alcock et al. 2000; OGLE-mNOVA-25 = OGLE-LMC-02, Ł. Wyrzykowski et al. 2009) although such events should be achromatic and nonrepeating if they really are due to microlensing. The list of all our selected objects is presented in Table 1.
Table 1. List of Detected Objects
Name | R.A. | Decl. | Iq | (V − I)q | ΔI | ΔT | τd | E(V − I) | Comments | ||
---|---|---|---|---|---|---|---|---|---|---|---|
(hh:mm:ss) | (deg:arcmin:arcsec) | (mag) | (mag) | (mag) | (mag) | (mag) | (day) | (day mag−1) | (mag) | ||
OGLE-mNOVA-01 | 20.273 ± 0.016 | 0.784 ± 0.056 | 16.595 ± 0.012 | 0.198 ± 0.008 | 3.7 | > 200 | 33 | 0.059 | ASASSN-16oh | ||
OGLE-mNOVA-02 | 20.981 ± 0.017 | 1.088 ± 0.058 | 19.216 ± 0.050 | ... | 1.8 | > 70 | ... | 0.042 | ⋯ | ||
OGLE-mNOVA-03 | 21.063 ± 0.040 | 0.502 ± 0.186 | 19.956 ± 0.055 | 0.501 ± 0.055 | 1.1 | 55 | 33 | 0.028 | ⋯ | ||
OGLE-mNOVA-04 | 21.008 ± 0.025 | 0.611 ± 0.239 | 19.706 ± 0.064 | ... | 1.3 | 83 | 26 | 0.046 | ⋯ | ||
OGLE-mNOVA-05 | 20.745 ± 0.017 | 0.894 ± 0.068 | 18.885 ± 0.022 | 0.420 ± 0.017 | 1.9 | 74 | 39 | 0.016 | candidate | ||
OGLE-mNOVA-06 | 20.916 ± 0.017 | 1.051 ± 0.073 | 19.268 ± 0.029 | 0.354 ± 0.033 | 1.6 | 102 | 20 | 0.011 | ⋯ | ||
OGLE-mNOVA-07 | 18.276 ± 0.010 | 0.574 ± 0.015 | 17.301 ± 0.012 | 0.545 ± 0.008 | 1.0 | > 235 | 213 | 0.068 | ⋯ | ||
OGLE-mNOVA-08 | 20.197 ± 0.013 | 0.987 ± 0.019 | 18.255 ± 0.024 | ... | 1.9 | 78 | 27 | 0.128 | ⋯ | ||
OGLE-mNOVA-09 | 20.057 ± 0.015 | 1.428 ± 0.064 | 19.015 ± 0.046 | 1.267 ± 0.028 | 1.0 | 115 | 105 | 0.122 | candidate | ||
OGLE-mNOVA-10 | 20.673 ± 0.026 | 0.789 ± 0.110 | 19.586 ± 0.091 | 0.771 ± 0.128 | 1.1 | 80 | 41 | 0.119 | ⋯ | ||
OGLE-mNOVA-11 | 20.884 ± 0.039 | 1.244 ± 0.129 | 18.151 ± 0.021 | 0.420 ± 0.023 | 2.7 | 142 | 29 | 0.100 | ⋯ | ||
OGLE-mNOVA-12 | 20.814 ± 0.033 | 0.776 ± 0.171 | 19.407 ± 0.029 | 0.570 ± 0.014 | 1.4 | 95 | 65 | 0.083 | MACHO-LMC-7 | ||
OGLE-mNOVA-13 | 20.096 ± 0.011 | 0.847 ± 0.018 | 18.909 ± 0.048 | ... | 1.2 | 131 | 30 | 0.124 | MACHO-LMC-23 | ||
OGLE-mNOVA-14 | 20.423 ± 0.018 | 0.372 ± 0.061 | 18.847 ± 0.058 | 0.156 ± 0.091 | 1.6 | 95 | 17 | 0.092 | ⋯ | ||
OGLE-mNOVA-15 | 20.790 ± 0.022 | 0.886 ± 0.078 | 18.540 ± 0.018 | 0.478 ± 0.015 | 2.2 | > 120 | 34 | 0.183 | candidate | ||
OGLE-mNOVA-16 | 20.545 ± 0.019 | 0.327 ± 0.031 | 19.202 ± 0.045 | 0.223 ± 0.061 | 1.3 | 42 | 24 | 0.072 | ⋯ | ||
OGLE-mNOVA-17 | 20.461 ± 0.015 | 0.893 ± 0.087 | 19.312 ± 0.065 | 0.810 ± 0.027 | 1.1 | 120 | 49 | 0.105 | ⋯ | ||
OGLE-mNOVA-18 | 20.366 ± 0.011 | 0.904 ± 0.018 | 18.570 ± 0.037 | 0.756 ± 0.078 | 1.8 | 105 | 58 | 0.094 | ⋯ | ||
OGLE-mNOVA-19 | 21.179 ± 0.031 | 0.838 ± 0.084 | 19.775 ± 0.037 | 0.541 ± 0.045 | 1.4 | 48 | 21 | 0.102 | ⋯ | ||
OGLE-mNOVA-20 | 19.648 ± 0.011 | 0.118 ± 0.018 | 18.259 ± 0.020 | 0.146 ± 0.024 | 1.4 | 108 | 35 | 0.078 | ⋯ | ||
OGLE-mNOVA-21 | 19.886 ± 0.011 | 0.656 ± 0.017 | 18.492 ± 0.025 | 0.498 ± 0.020 | 1.4 | > 110 | 39 | 0.078 | ⋯ | ||
OGLE-mNOVA-22 | 18.574 ± 0.010 | 0.580 ± 0.014 | 17.653 ± 0.015 | ... | 0.9 | > 800 | ... | 0.139 | ⋯ | ||
OGLE-mNOVA-23 | 21.167 ± 0.047 | 0.700 ± 0.302 | 20.006 ± 0.051 | 0.451 ± 0.029 | 1.2 | 33 | 14 | 0.099 | ⋯ | ||
OGLE-mNOVA-24 | 21.032 ± 0.048 | 0.671 ± 0.094 | 18.661 ± 0.036 | 0.487 ± 0.027 | 2.4 | > 280 | 48 | 0.082 | ⋯ | ||
OGLE-mNOVA-25 | 20.404 ± 0.014 | 0.428 ± 0.022 | 19.228 ± 0.080 | 0.520 ± 0.059 | 1.2 | 81 | 39 | 0.060 | OGLE-LMC-02 | ||
OGLE-mNOVA-26 | 20.764 ± 0.027 | 0.097 ± 0.116 | 19.366 ± 0.057 | −0.023 ± 0.025 | 1.4 | 43 | 148 | 0.124 | ⋯ | ||
OGLE-mNOVA-27 | 20.763 ± 0.028 | 1.010 ± 0.025 | 19.273 ± 0.069 | ... | 1.5 | 40 | 14 | 0.263 | ⋯ | ||
OGLE-mNOVA-28 | 20.818 ± 0.023 | 1.079 ± 0.227 | 19.439 ± 0.034 | 0.537 ± 0.055 | 1.4 | 85 | 26 | 0.141 | ⋯ | ||
OGLE-mNOVA-29 | 20.747 ± 0.019 | 0.410 ± 0.027 | 19.758 ± 0.078 | 0.655 ± 0.088 | 1.0 | > 120 | ... | 0.105 | ⋯ |
Note. The table provides equatorial coordinates (for the epoch J2000), mean I-band brightness and V − I color in quiescence, mean I-band brightness and V − I color in outburst, amplitude of the outburst in I-band ΔI, duration of the highest-amplitude outburst ΔT, mean decline rate τd, and color excess E(V − I) (D. M. Skowron et al. 2021) toward the detected objects.
Download table as: ASCIITypeset image
Given the unique spectral and X-ray characteristics of ASASSN-16oh, we also decided to start the near-real-time monitoring of these selected objects in 2023 September, with the goal of triggering follow-up observations should any enter an outburst state.
4. Outburst of OGLE-mNOVA-11
The outburst of OGLE-mNOVA-11 started shortly afterwards, between 2023 October 15.3 and 26.3 and reached a maximum of I = 18.15 ± 0.02 and V = 18.57 ± 0.01 on 2023 December 6.1 (Figure 2), an amplitude of almost 3 mag. This corresponds to MI = −0.5 ± 0.1 and MV = −0.2 ± 0.1, assuming the LMC distance modulus of μ = 18.477 ± 0.004 (G. Pietrzyński et al. 2019) and reddening (D. M. Skowron et al. 2021). The star returned to quiescence (in I) in late 2024 February, about 120 days after the outburst onset. The outburst duration was shorter in V.
Download figure:
Standard image High-resolution imageWe obtained a set of low-resolution spectra of OGLE-mNOVA-11 from 2023 November 22 to December 10 (Figure 3) with the Robert Stobie Spectrograph (H. A. Kobulnicky et al. 2003; E. B. Burgh et al. 2003) mounted on the Southern African Large Telescope (SALT), as part of the SALT Large Science Programme on Transients (Appendix B). The spectra cover the range 3500–6700 Å and reveal narrow (FWHM = 247 ± 28 km s−1) emission lines (Balmer lines, He ii 4686 Å). The Bowen blend (a complex of C iii and N iii lines in the range 4640–4650 Å) is also clearly detected from 2023 December 2 to 10. The emission lines are redshifted with a mean radial velocity of 278.3 ± 4.5 km s−1, close to that of the LMC systemic value (262.2 ± 3.4 km s−1; R. P. van der Marel et al. 2002). Both the light-curve shape and optical spectra are very different from those of classical novae.
Download figure:
Standard image High-resolution imageOGLE-mNOVA-11 was also observed by the Neil Gehrels Swift Observatory (N. Gehrels et al. 2004) five times from 2023 December 6 to 30 (Appendix C). Swift's X-Ray Telescope (XRT; D. N. Burrows et al. 2005) detected a faint X-ray source at the transient's position, with a 0.3–10 keV count rate declining from count s−1 (December 6) to count s−1 (December 30); see Figure 2. All the data were combined to create a single X-ray spectrum, which can be fitted with a blackbody of K absorbed by a column density cm−2. If we fix NH to the LMC value (1.3 × 1021 cm−2), then the temperature is slightly better constrained to K. Assuming this blackbody model, the observed LX (0.3–10 keV) ≈ 0.9 × 1035 erg s−1, or 3.6 × 1035 erg s−1 when corrected for absorption, about an order of magnitude smaller than that of ASASSN-16oh, but still substantially higher than dwarf nova X-ray luminosities.
5. Discussion and Conclusions
Overall, the optical (light-curve shape, narrow emission lines, strong He II emission) and X-ray properties of OGLE-mNOVA-11 match well with those of ASASSN-16oh. We thus propose that they, and probably the remaining objects described in this study, form a homogeneous group of transient supersoft X-ray sources. We dub them "millinovae" because their optical luminosities are roughly a 1000 times fainter than those of ordinary classical novae.11
Using ASASSN-16oh and OGLE-mNOVA-11 as prototypes, millinovae should fulfill the following criteria:
- 1.They exhibit symmetrical, triangle-shaped outbursts in the optical bands (V, I).
- 2.Outbursts last from a month to several months, i.e., substantially longer than the typical outbursts of short-orbital-period dwarf novae.
- 3.The peak absolute magnitude in the optical bands is between that of classical novae and short-orbital-period dwarf novae.
- 4.The optical spectra show narrow (FWHM smaller than a few hundred km s−1) emission lines of H and He ii (that is, no signatures of mass ejection are seen in the spectra).
- 5.Soft X-ray emission (declining with decreasing optical luminosity) is seen during the outburst.
Only two systems (ASASSN-16oh and OGLE-mNOVA-11) meet all these five criteria. The remaining objects, given the lack of dedicated spectroscopic and X-ray observations during outbursts, remain millinova candidates.
Their sky location (upper panels of Figure 4) and the fact that most of the selected stars occupy a relatively narrow region in the color–magnitude diagram (middle left panel of Figure 4) indicate that millinovae are all in the Magellanic Clouds. We used the long-term astrometric time series of OGLE observations to measure proper motions of selected objects (Appendix D), which are consistent (with one exception) with those of stars located in the Magellanic Clouds (Gaia Collaboration et al. 2018). The radial velocities of ASASSN-16oh (T. J. Maccarone et al. 2019) and OGLE-mNOVA-11 (Appendix B) are consistent with the systemic velocities of the SMC and LMC, respectively. However, we cannot rule out that some objects are located in the foreground Milky Way disk or in background galaxies although we consider this very unlikely.
Download figure:
Standard image High-resolution imageThe middle left panel of Figure 4 shows that most objects in quiescence occupy the same region of the color–magnitude diagram, 0.5 ≤ (V − I)0 ≤ 1.0, 1.5 ≤ MI ≤ 2.5, which indicates that the accretion disk dominates their quiescent luminosity. Such a range is consistent with the absolute magnitudes of dwarf novae in the Magellanic Clouds (M. M. Shara et al. 2003), indicating that they are essentially nova-like disks. In outburst (middle right panel of Figure 4), objects move toward bluer colors (0.0 ≤ (V − I)0 ≤ 0.5) and higher luminosities (−0.5 < MI < 1.0), as expected given the strong irradiation of the disk by the supersoft component.
The mean absolute magnitudes of millinovae in quiescence are MI = 1.72 ± 0.70 and MV = 2.37 ± 0.80. In outburst, the absolute magnitudes are MI = 0.18 ± 0.83 and MV = 0.62 ± 0.97 and are about a 1000 times fainter than classical novae (lower right panel of Figure 4). Millinovae are also much brighter than most known dwarf novae (lower panels of Figure 4).
The spectroscopic properties of ASASSN-16oh and OGLE-mNOVA-11 and the (V − I) color evolution of OGLE-mNOVA-11 indicate that their optical outbursts likely have a disk origin. Because dwarf novae follow a well-known relation between the orbital period and absolute magnitude (both in quiescence and outburst; J. Patterson 2011), relatively high absolute magnitudes of millinovae indicate that their orbital periods should be on the order of a few days. Indeed, one of the analyzed objects (OGLE-mNOVA-08) shows clear eclipsing variability with a 4.830893(47) day period (Figure 5). (We searched for possible periodic brightness variations in quiescence and during outbursts in all 29 objects but detected a statistically significant signal only in the case of OGLE-mNOVA-08.)
Download figure:
Standard image High-resolution imageFurthermore, the rate of decline from the outburst of dwarf novae, τd, is also known to be proportional to Porb (the "Bailey relation"; J. Bailey 1975; see Figure 6), as expected, given that longer periods imply larger disks. We find that most millinovae have decline rates in the range 20 < τd < 70 day mag−1, which corresponds to Porb from 3 to 15 days, assuming that the Bailey relation can be extrapolated to such long periods.
Download figure:
Standard image High-resolution imageThe "spreading layer" (R. Kippenhahn & H. C. Thomas 1978; N. A. Inogamov & R. A. Sunyaev 1999; A. L. Piro & L. Bildsten 2004; T. J. Maccarone et al. 2019) model can naturally produce the soft X-ray temperatures that have been seen in many dwarf nova outbursts, albeit in shorter-period systems such as U Gem and SS Cyg (see, e.g., P. J. Wheatley et al. 2003). However, the long orbital periods of millinovae imply the existence of large disks and hence likely large mass-transfer rates (e.g., J. Echevarrìa 1994) capable of producing higher X-ray luminosities. Nevertheless, there are difficulties with this model that have been noted (e.g., M. Kato et al. 2020) in that there are high mass-transfer-rate recurrent novae in quiescence (e.g., U Sco, RS Oph) where such a soft component would have been expected but has never been seen.
The most obvious solution to explaining the properties of millinovae is then to find a way to trigger thermonuclear burning on the white dwarf surface without producing a nova flash. One such attempt to do this is the nonejecting thermonuclear runaway event of Y. Hillman et al. (2019), but the resulting light curves of their ASASSN-16oh model are completely different to the smooth rise and fall that are observed, not to mention the huge color variation they predict, which is definitely excluded by the data.
An alternative route to such a trigger has been sought by M. Kato et al. (2020), who build on the already known, ongoing high mass-transfer rate in these systems (their brightness and variability properties are consistent with nova-like in the Magellanic Clouds; M. M. Shara et al. 2003), thereby keeping the disk (and white dwarf) in a hot, active state. Furthermore, they require a massive white dwarf, close to 1.2 M⊙ in order to initiate steady burning without a strong wind or shell being ejected. While we do not yet know how many of these long symmetric outburst systems produce detectable supersoft X-rays, if they were all to do so, then this would have important implications for the theory of nonejecting supersoft X-ray sources, as it would constrain the parameters of systems capable of behaving in this way. Furthermore, it would provide a route for allowing the white dwarf to continue to grow in mass and hence become a potential Type Ia supernova progenitor.
We recognize, however, that an added complication to the M. Kato et al. (2020) models is that they produce close to Eddington-limited events, whereas some of our millinovae clearly do not reach such a level in X-rays (Appendix E). There may be severe geometrical constraints on the fluxes observed, or the nuclear burning might not cover the entire white dwarf surface. Moreover, some of these models (see also M. Kato et al. 2022) predict the highest, close-to-Eddington flux to be in the far-UV range (rather than soft X-rays), which can be tested with dedicated follow-up observations of future millinova outbursts. Nevertheless, we believe that this new group of millinovae, all likely representing long-period, high mass-transfer-rate cataclysmic variables, opens an important new route for study, with the added benefit of being a well-constrained population in the Magellanic Clouds.
Acknowledgments
We thank all the OGLE observers for their contribution to the collection of the photometric data over the decades. This research was funded in part by National Science Centre, Poland, grant OPUS 2021/41/B/ST9/00252 awarded to P.M. K.L.P. acknowledges funding from the UK Space Agency. A part of this work is based on observations made with the Southern African Large Telescope (SALT), with the Large Science Programme on transients 2021-2-LSP-001 (PI: D.A.H.B.). Polish participation in SALT was funded by grant No. MEiN 2021/WK/01.
Facility: OGLE, Swift, SALT, eROSITA. -
Appendix A: Light Curves
Figures A1 and A2 show the light curves of the remaining stars exhibiting long-duration symmetrical outbursts, which we have classified as "millinovae."
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageAppendix B: SALT Observations of OGLE-mNOVA-11
Optical spectra of OGLE-mNOVA-11 were obtained with the Robert Stobie Spectrograph (H. A. Kobulnicky et al. 2003; E. B. Burgh et al. 2003) mounted on SALT (D. A. H. Buckley et al. 2006) in long-slit mode. Two 1200 s exposure spectra were taken each night on 2023 November 22 and December 2, 3, 6, 8, and 10, using the PG0700 grating at an angle of 3∘ with a slit width of 15, giving a resolving power of 460 < R < 800. The data were prereduced using the PySALT package (S. M. Crawford et al. 2010), which includes bias subtraction, gain and cross-talk corrections, and mosaicking. An Argon lamp was utilized for wavelength calibration, and the spectrophotometric standard star, EG 21, was used for flux calibration, both performed with IRAF. The spectra were stacked so as to obtain one average spectrum per night.
The software package, Fityk12 (M. Wojdyr 2010), was used for spectral analysis. The SALT spectra were normalized before the emission lines were modeled using the method of A. W. Shafter (1983) in order to determine the central wavelengths of the He ii (4686 Å), Hβ, and Hα emission lines.. This entailed fitting a Gaussian function through 95% of the line wings. In Fityk, the errors of the fitted parameters were obtained for a confidence level of 1σ (68.27%). The measured central wavelengths were then utilized to calculate the barycentric-corrected radial velocities.
Table B1 presents the barycentric radial velocities of He ii (4686 Å), Hβ, and Hα. Because of the larger scatter in our Hβ velocities (due to its lower-significance detections), we also calculated the weighted mean radial velocity from He ii and Hα, obtaining values that ranged from 232.1 to 311.3 km s−1, with a mean of 278.3 ± 4.5 km s−1, which is very close to the established LMC line-of-sight velocity (262.2 ± 3.4 km s−1; R. P. van der Marel et al. 2002). The number of radial velocity measurements prevents us from determining the orbital period unequivocally, but it is likely to be shorter than 4 days. If we restrict orbital periods to be longer than 1 day and assume a circular orbit, the best-fit periods are 2.42, 1.68, 2.08, and 1.89 days.
Table B1. SALT Radial Velocities of OGLE-mNOVA-11
Time | BJD | He ii | Hβ | Hα | Mean RV |
---|---|---|---|---|---|
(km s−1) | (km s−1) | (km s−1) | (km s−1) | ||
2023 Nov 22.93 | 2460271.43 | 313.7 ± 25.7 | 199.9 ± 132.0 | 310.9 ± 11.7 | 311.3 ± 10.6 |
2023 Dec 2.88 | 2460281.38 | 308.5 ± 8.6 | 329.4 ± 25.6 | 306.7 ± 18.2 | 308.2 ± 7.8 |
2023 Dec 3.84 | 2460282.34 | 203.0 ± 15.3 | 203.0 ± 99.5 | 279.8 ± 18.4 | 234.3 ± 11.8 |
2023 Dec 6.84 | 2460285.34 | 224.1 ± 19.1 | 475.5 ± 95.5 | 252.8 ± 30.9 | 232.1 ± 16.3 |
2023 Dec 8.83 | 2460287.33 | 284.2 ± 22.3 | 143.1 ± 59.5 | 241.8 ± 22.7 | 263.4 ± 15.9 |
2023 Dec 10.83 | 2460289.33 | 245.2 ± 13.4 | 417.6 ± 124.2 | 257.4 ± 19.7 | 249.0 ± 11.1 |
Download table as: ASCIITypeset image
Appendix C: Swift Observations of OGLE-mNOVA-11
C.1. Swift X-Ray Data
The Neil Gehrels Swift Observatory (N. Gehrels et al. 2004) performed five Target of Opportunity observations of OGLE-mNOVA-11 (Target ID 16408) between 2023 December 6 and December 30, collecting ≈ 1.5 ks during each pointing, with the exception of December 17, when only ≈ 0.5 ks was obtained. The XRT (D. N. Burrows et al. 2005) was operated in photon counting mode, and the data were processed using HEASoft 6.32, together with the most up-to-date calibration files available. Given the faintness of the source, a small extraction region of 10 pixels (236) radius was used for the source, and a larger, 60 pixel (1414) radius circle, offset from, but close to, the source was used to estimate the background contribution.
In observations where the source was detected at the 3σ level, the 0.3–10 keV count rate (corrected for all point-spread function losses) is given in Table C1. In the shorter exposure on December 17, the source was not detected, and a 3σ upper limit is given instead.
Table C1. Swift XRT Observations of OGLE-mNOVA-11
Obs. ID | Start Time | Exp. Time | Count Rate |
---|---|---|---|
(UT) | (s) | (10−3 s−1) | |
00016408001 | 2023-12-06 16:04:57 | 1594 | |
00016408002 | 2023-12-13 09:21:57 | 1499 | 8 ± 2 |
00016408003 | 2023-12-17 07:10:56 | 479 | < 20 |
00016408004 | 2023-12-20 04:27:56 | 1605 | 11 ± 2 |
00016408005 | 2023-12-30 16:46:56 | 1529 |
Download table as: ASCIITypeset image
There is no significant change in the hardness ratio during these observations though it should be noted that we are dealing with a limited number of source counts (<50 in total across all five observations, the vast majority below 1 keV). Therefore, all the data were combined to create a single spectrum. Cash statistics (W. Cash 1979) were used to fit this spectrum with a combination of an optically thick blackbody for the soft emission and an optically thin thermal plasma to account for the few higher-energy photons, both absorbed by the same hydrogen column. While the optically thin plasma can only be constrained to have a temperature of >4.2 keV (>4.9 × 107 K), the best-fitting blackbody temperature is 36 eV (418,000 K), absorbed by a column density of cm−2. This fit results in a C-stat value of 19.9 for 25 degrees of freedom.
If NH is fixed to the LMC value of 1.3 × 1021 cm−2 (HI4PI Collaboration et al. 2016), the blackbody temperature becomes slightly better constrained, with eV (T = 607,000 K), for C-stat/degrees of freedom = 21.2/26. This fit leads to an LX (0.3–10 keV) of 9.0 × 1034 erg s−1, or 3.6 × 1035 erg s−1 unabsorbed (assuming a distance of 49.59 kpc; G. Pietrzyński et al. 2019).
C.2. Swift UV Data
OGLE-mNOVA-11 was also observed by the Ultraviolet and Optical Telescope (UVOT; P. W. A. Roming et al. 2005) on board the Neil Gehrels Swift Observatory. Observations took place during a period when one of the three onboard gyroscopes malfunctioned, resulting in increased noise and affecting the image quality of UVOT data (B. Cenko 2023). Moreover, the nearest source detected in the UVOT images was slightly offset from the target position by 36–48 and was likely a blend of the target and a nearby (54) bright (V = 16.996, V − I = −0.097) constant star. We measured the brightness of the blend using the UVOTSOURCE tool, which performs aperture photometry. We used an aperture with a radius of 15'' centered on the position of OGLE-mNOVA-11.
We then estimated the brightness of the neighbor from the analysis of the UV–optical color–magnitude diagrams, assuming that the neighbor's color is consistent with colors of other main-sequence stars of similar magnitude. We, therefore, created UV–optical color–magnitude diagrams by cross-matching the UVOT catalogs with the OGLE photometric map. We then calculated the mean UV − V color of main-sequence stars in the brightness range ∣V − 16.996∣ < Δ, where Δ = 0.15 mag. We found
(We also tested Δ = 0.1 and Δ = 0.2 and obtained similar results.) That enabled us to estimate the brightness of the neighbor in the UV filters:
and estimate the UV brightness of the target by subtracting the brightness of the neighbor star. Results are reported in Table C2. All magnitudes are in the Vega system.
Table C2. Swift UVOT Observations of OGLE-mNOVA-11
Obs. ID | uvw1 | uvm2 | uvw2 |
---|---|---|---|
00016408001 | ... | 15.838 ± 0.046 | ... |
00016408002 | 16.448 ± 0.094 | 16.042 ± 0.078 | 16.019 ± 0.068 |
00016408003 | ... | 15.843 ± 0.068 | ... |
00016408004 | ... | 16.166 ± 0.059 | ... |
00016408005 | ... | 16.197 ± 0.074 | ... |
Download table as: ASCIITypeset image
Appendix D: OGLE Proper Motions
We calculated the proper motions of detected stars using the OGLE time-series astrometric data (OGLE-Uranus; A. Udalski et al. 2024, in preparation). In short, we measured each star's position (x, y) in individual OGLE images collected between 2010 and 2024 (in some cases, 2010–2020 if later data were unavailable). We then cross-matched the positions of bright stars with those calculated using the Gaia DR3 data (Gaia Collaboration et al. 2016, 2023), taking into account their proper motion and heliocentric parallax. That allowed us to find the transformation between the positions (x, y) in individual OGLE images to the equatorial coordinates in the Gaia reference frame (L. Lindegren et al. 2021). The proper motions in R.A. μα and decl. μδ for all detected stars are reported in Table D1. Because most of our targets are very faint (I ≈ 20–21), the proper motion uncertainties are large (the median error bar is 2 mas yr−1). Nonetheless, in all cases (but one), the proper motions are consistent within the quoted errors with the Magellanic Clouds' proper motions (Gaia Collaboration et al. 2018). The only exception is OGLE-mNOVA-15 with a total proper motion of 23.5 ± 3.8 mas yr−1, making it most likely to be a foreground object, and we therefore treat it as a candidate millinova.
Table D1. Millinovae Proper Motions
Object | μα (OGLE) | μδ (OGLE) | μα (Gaia) | μδ (Gaia) |
---|---|---|---|---|
(mas yr−1) | (mas yr−1) | (mas yr−1) | (mas yr−1) | |
OGLE-mNOVA-01 | 1.66 ± 2.03 | 0.14 ± 2.39 | ⋯ | ⋯ |
OGLE-mNOVA-02 | 1.31 ± 2.35 | −2.50 ± 2.31 | ⋯ | ⋯ |
OGLE-mNOVA-03 | 1.60 ± 3.22 | −4.03 ± 2.68 | ⋯ | ⋯ |
OGLE-mNOVA-04 | 1.25 ± 2.78 | −3.15 ± 2.53 | ⋯ | ⋯ |
OGLE-mNOVA-05 | −0.37 ± 1.44 | −0.44 ± 1.49 | 0.55 ± 0.51 | −1.19 ± 0.65 |
OGLE-mNOVA-06 | −0.67 ± 2.04 | 0.29 ± 2.16 | ⋯ | ⋯ |
OGLE-mNOVA-07 | 0.74 ± 0.09 | −1.39 ± 0.10 | 0.80 ± 0.17 | −1.18 ± 0.18 |
OGLE-mNOVA-08 | 0.85 ± 0.91 | 0.84 ± 0.92 | ⋯ | ⋯ |
OGLE-mNOVA-09 | −1.26 ± 1.16 | 0.08 ± 1.16 | ⋯ | ⋯ |
OGLE-mNOVA-10 | 2.13 ± 1.26 | −0.13 ± 1.51 | ⋯ | ⋯ |
OGLE-mNOVA-11 | 4.35 ± 2.18 | 4.08 ± 2.37 | ⋯ | ⋯ |
OGLE-mNOVA-12 | −3.00 ± 2.51 | −7.29 ± 2.64 | ⋯ | ⋯ |
OGLE-mNOVA-13 | 2.11 ± 0.74 | −0.25 ± 0.66 | 2.53 ± 2.28a | 10.09 ± 2.60a |
OGLE-mNOVA-14 | 1.96 ± 2.56 | −2.29 ± 1.99 | ⋯ | ⋯ |
OGLE-mNOVA-15 | 12.35 ± 3.04 | −20.00 ± 4.11 | ⋯ | ⋯ |
OGLE-mNOVA-16 | 3.32 ± 1.37 | −9.88 ± 3.42 | ⋯ | ⋯ |
OGLE-mNOVA-17 | 1.44 ± 1.32 | 2.71 ± 1.87 | ⋯ | ⋯ |
OGLE-mNOVA-18 | 0.91 ± 1.40 | 0.38 ± 1.67 | ⋯ | ⋯ |
OGLE-mNOVA-19 | −4.35 ± 3.39 | −1.48 ± 2.37 | ⋯ | ⋯ |
OGLE-mNOVA-20 | −0.34 ± 1.67 | 1.88 ± 0.51 | ⋯ | ⋯ |
OGLE-mNOVA-21 | 0.19 ± 0.61 | 0.02 ± 0.68 | ⋯ | ⋯ |
OGLE-mNOVA-22 | 1.82 ± 0.09 | 0.71 ± 0.10 | 2.24 ± 0.19 | 1.16 ± 0.24 |
OGLE-mNOVA-23 | −1.93 ± 4.22 | 3.85 ± 6.32 | ⋯ | ⋯ |
OGLE-mNOVA-24 | 6.43 ± 4.11 | 1.45 ± 2.32 | ⋯ | ⋯ |
OGLE-mNOVA-25 | −0.27 ± 2.01 | −0.13 ± 2.81 | ⋯ | ⋯ |
OGLE-mNOVA-26 | 4.13 ± 1.71 | 1.49 ± 0.97 | ⋯ | ⋯ |
OGLE-mNOVA-27 | −0.33 ± 1.44 | 3.21 ± 2.41 | ⋯ | ⋯ |
OGLE-mNOVA-28 | −1.22 ± 2.05 | 0.05 ± 1.94 | ⋯ | ⋯ |
OGLE-mNOVA-29 | −1.85 ± 2.75 | −0.54 ± 2.33 | ⋯ | ⋯ |
Note. aThe Gaia DR3 (Gaia Collaboration et al. 2016, 2023) astrometric solution is unreliable because the renormalized unit weight error is larger than 1.4 (L. Lindegren et al. 2021).
Download table as: ASCIITypeset image
Appendix E: eROSITA Upper Limits
The first catalog of X-ray sources detected by the eROSITA telescope array (P. Predehl et al. 2021) on board the Spectrum-Roentgen-Gamma (SRG; R. Sunyaev et al. 2021) satellite was recently published by A. Merloni et al. (2024). The number of exposures, and so the depth of the survey, depends on sky position. Thanks to the SRG scanning law, objects located near the ecliptic poles are observed most frequently. Coincidentally, the LMC is located near the south ecliptic pole.
We cross-matched our list with the eROSITA catalog (A. Merloni et al. 2024), but none of the objects from Table 1 were detected by eROSITA. The catalog of A. Merloni et al. (2024) contains X-ray sources located in the western Galactic hemisphere detected during the first 6 months of eROSITA operations (from 2019 December 12 to 2020 June 11, that is, 2,458,829 < JD < 2,459,011). The overlap with OGLE observations is relatively small because the OGLE-IV operations were suspended on 2020 March 18 (JD = 2,458,927) due to the COVID-19 pandemic. Nonetheless, we found that four objects (OGLE-mNOVA-06, OGLE-mNOVA-10, OGLE-mNOVA-19, and OGLE-mNOVA-24) were in outburst during the first eROSITA All-Sky Survey. All four objects were relatively faint and reached a peak I ∼ 19.55–19.98 during that time. Assuming that the X-ray flux scales proportionally with peak I-band flux, and taking OGLE-mNOVA-11 as a benchmark, the maximal expected X-ray flux is ∼ (6–8) × 10−14 erg cm−2 s−1 in the 0.3–10 keV range. The expected flux in the eROSITA soft and medium energy bands (0.2–2.3 keV) is about 6% larger at (6–9) × 10−14 erg cm−2 s−1. However, each object must have been scanned by eROSITA multiple times, and the time-averaged flux is likely to be smaller.
We queried the eROSITA upper limits service (M. Krumpe et al. 2024) to get upper flux limits in the 0.2–2.3 keV range for all four objects discussed above. These limits are calculated at a 99.87% (3σ) one-sided confidence level, assuming a power-law spectrum with a photon index Γ = 2.0 and a column density NH = 3 × 1020 cm−2. We obtained limits of (5.1, 1.9, 2.0, 3.8) × 10−14 erg cm−2 s−1 at the positions of OGLE-mNOVA-(06, 10, 19, 24), respectively. We followed a procedure outlined in M. Krumpe et al. (2024) to estimate upper limits for the best-fitting absorbed blackbody model presented in Appendix C. Using the appropriate eROSITA calibration files and XSPEC version 12.14.0 (K. A. Arnaud 1996), we estimated the energy-to-count conversion factor of 1.037 × 1012 cm2 erg−1, which is only 3.4% smaller than that calculated using the fiducial model in M. Krumpe et al. (2024). Therefore, the upper limits for our best-fitting model are only 3.4% larger than default ones. In the 0.2–5 keV range, the corresponding limits are (7.5, 3.7, 2.5, 5.4) × 10−14 erg cm−2 s−1. We thus conclude that the expected time-averaged flux is close to or below the current eROSITA limits. It is therefore possible that the selected stars may be detected in future eROSITA data releases once more data are processed or may be detected in data from individual scans.
Footnotes
- 11
They should not be confused with "micronovae," a much lower-luminosity and much faster event seen in a number of Galactic magnetic cataclysmic variables (S. Scaringi et al. 2022).
- 12