{'fraction_non_alphanumeric': 0.041058829770303965, 'fraction_numerical': 0.05087263275157817, 'mean_word_length': 4.17201646090535, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The shortest period sdB plus white dwarf binary CD-30 11223 (GALEX J1411-3053)

We report on the discovery of the shortest period binary comprising a hot subdwarf star (CD-30 11223, GALEX J1411-3053) and a massive unseen companion. Photometric data from the All Sky Automated Survey show ellipsoidal variations of the hot subdwarf primary and spectroscopic series revealed an orbital period of 70.5 minutes. The large velocity amplitude suggests the presence of a massive white dwarf in the system (M_2/M_sun>0.77) assuming a canonical mass for the hot subdwarf (0.48 M_sun), although a white dwarf mass as low as 0.75 M_sun is allowable by postulating a subdwarf mass as low as 0.44 M_sun. The amplitude of ellipsoidal variations and a high rotation velocity imposed a high-inclination to the system (i>68 deg) and, possibly, observable secondary transits (i>74 deg). At the lowest permissible inclination and assuming a subdwarf mass of ~0.48 M_sun, the total mass of the system reaches the Chandrasekhar mass limit at 1.35 M_sun and would exceed it for a subdwarf mass above 0.48 M_sun. The system should be considered, like its sibling KPD 1930+2752, a candidate progenitor for a Type Ia supernova. The system should become semi-detached and initiate mass transfer within ~30 Myr.


INTRODUCTION
A possible scenario for triggering type Ia supernovae involves a degenerate core collapse following the merger of two white dwarfs (see a recent summary in Ferrario 2012, and recent population syntheses by Toonen et al. (2012)). A search for double degenerate stars that could progress toward such an outcome remains unfulfilled with most binaries either too lean or too wide (SPY program, see Napiwotzki et al. 2001), with the possible exception of WD 2020−425 (Napiwotzki et al. 2005). A more promising variation on this scenario involves the interaction of a hot subdwarf star with a close white dwarf companion. Mass transfered from a He-star onto a CO white dwarf may ignite and even lead to detonation of the carbon core. At an early stage these objects appear as AM CVn systems (see Yungelson 2008;Solheim 2010). Hot subdwarfs, also referred to as extreme horizontal branch (EHB) stars are core helium burning stars that formed preferentially in the aftermath of close, postinteracting binaries (Han et al. 2002(Han et al. , 2003. The identification of the short-period subdwarf plus white dwarf binary KPD 1930+2752 (Maxted et al. 2000;Geier et al. 2007) lent support to the idea, and the expectation of conservative mass transfer gave KPD 1930+2752 the character of a proper type Ia candidate progenitor. A search for other massive subluminous companions to hot subdwarf stars has so far failed to provide further candidates (MUCHFUSS program, see Geier et al. 2011Geier et al. , 2012a. In this letter, we report on the discovery of a close companion to the hot subdwarf CD−30 11223. The or-bital period is the shortest known to date for such systems. The observations were carried-out as part of an intermediate-dispersion survey of ∼40 hot, ultravioletselected subdwarf stars from the catalog of Vennes et al. (2011) andNemeth et al. (2012) and aimed at constraining the fraction of close, evolved binaries in the EHB population. Kawka et al. (2010) described the orbital properties of the first two systems discovered as part of this survey.
The hot hydrogen-rich subdwarf (sdB) CD−30 11223 is relatively bright (V ≈ 12) and nearby (d = 250 − 400 pc) and was first classified by Vennes et al. (2011) in spectroscopic follow-up observations of ultraviolet sources detected in the Galaxy Evolution Explorer (GALEX). Based on low-dispersion spectra, they originally determined T eff = 26500 ± 2700 K, log g = 5.36 ± 0.26, and an abundance of helium log He/H = −1.34 ± 0.24. These parameters were refined by Nemeth et al. (2012) using higher-dispersion spectra: T eff = 30150 +290 −310 K, log g = 5.72 +0.04 −0.07 , log He/H = −1.58 +0.03 −0.05 . Nemeth et al. (2012) noted that temperature and surface gravity measurements in low dispersion spectra were systematically offset from measurements obtained at higher-dispersion and are probably less reliable. Moreover, Nemeth et al. (2012) reported the detection of nitrogen lines and measured an abundance log N/H = −4.40 +0.29 −0.67 . They found no evidence of carbon or oxygen in the atmosphere of this star. Remarkably, large radial velocity excursions between the two series of spectra were evident and prompted further investigations. Independently, photometric variations in a large ASAS data set confirmed our suspicions that this object is in a close, evolved binary. The nature of the variations (reflection effect or tidal distortion) as well as the nature of the companion may be established with joint phasing of photometric and spectroscopic data.
In the following sections we describe the spectroscopic and photometric data (Section 2) and their analysis (Section 3). Section 3.1 gives the results of a frequency analysis of the photometric data and Section 3.2 provides corroborating evidence of the short orbital period based on spectroscopic data. We also measured the mass function of the secondary and rotation velocity of the primary and analyzed the ellipsoidal variations to constrain the nature of the secondary star (Section 3.3). Finally, we summarize in Section 4.

OBSERVATIONS
The first identification spectra of the hot subdwarf CD−30 11223 were obtained with the European Southern Observatory (ESO) Faint Object Spectrograph and Camera (EFOSC2) attached to the 3.6 m New Technology Telescope (NTT) at La Silla Observatory on UT 2009 August 24 (ESO P83 program). We used grism no. 11 which has 300 lines per mm and a blaze wavelength at 4000Å. The slit width was set to 1 arcsec, which resulted in a spectral resolution of ∼14Å and a coverage from 3700 to 7250Å. We re-observed CD−30 11223 at the NTT on UT 2010 September 21 (ESO P85 program). We employed grism no. 7 (600 lines per mm) with a dispersion of 1.96Å per pixel. The slit width was set at 1 arcsec resulting in a resolution of ∼8Å . Although the 2009 and 2010 spectra only offered modest velocity accuracy (50 to 60 km s −1 ) the difference in measurements was close to 800 km s −1 suggesting the presence of a massive, close companion to the hot subdwarf. Consequently, we secured a set of intermediate dispersion spectra, first at Siding Spring Observatory (SSO) using the 2.3m telescope (April 2012) and at ESO using the NTT (August-September 2012).
We obtained a set of three spectra on UT 2012 April 29 and 30 with the Wide Field Spectrograph (WiFeS Dopita et al. 2007) attached to the 2.3m telescope at SSO. We used the B3000 and R7000 gratings which provided spectral ranges of 3200-5900Å at a resolution R = λ/∆λ = 3000, and 5300-7000Å at R = 7000, respectively. We used the RT560 dichroic beam-splitter to separate the incoming light into its red and blue components. We maximized the signal-to-noise of each observation by extracting the spectrum from the most significant ( 6) traces. Each trace was wavelength and flux calibrated prior to co-addition. The spectra were wavelength calibrated using NeAr arc spectra that were obtained following each observation. The exposure time was 600 s. The two exposures acquired on April 29 were separated by only 45 minutes but revealed a velocity shift close to 500 km s −1 .
Finally, we obtained complete orbital cycles at the NTT from UT 2012 August 30 to September 2 (ESO P89 program). We used grism no. 20 centered on Hα and with a dispersion of 0.55Å per pixel in the spectral range from 6040 to 7140Å. We set the slit width to 0.7 arcsecond resulting in a ∼ 2Å resolution. We set the exposure times to 120 s in order to minimize orbital smearing. We obtained 8 spectra on the night of August 31 and 6 spectra on September 1. The spectra revealed short (≈ 70 min) high-amplitude velocity variations and is also present at f phot /2. (Upper panel) ASAS photometric data folded on the best frequency and best fit at f phot /2. appeared to cover several orbits.
The field near CD−30 11223 was also covered with the All Sky Automated Survey (ASAS, Pojmanski 1997Pojmanski , 2001) on numerous occasions. A total of 1060 images were obtained as part of ASAS using the photometric V band. The star CD−30 11223 was not immediately identified as a variable star. We retrieved the photometric measurements using the ASAS all star catalogue page 2 .

Photometric frequency analysis
The ASAS measurements span nearly 3200 days from UT 2000 December 24 to 2009 September 13. Figure 1 shows a frequency analysis of the ASAS data. Fitting a sinusoidal function to the data set unambiguously selected a signal frequency f phot = 40.83377 ± 0.00003 d −1 with a semi-amplitude ∆m/2 = 0.047 ± 0.004 mag and mean magnitude m V = 12.342±0.003. A reflection effect on a fainter companion would display sinusoidal variations on the orbital period P orb (= 2 f −1 phot ), while tidal distortion of the primary star would induce sinusoidal variations on half the orbital period P orb /2 (= f −1 phot ) as shown in Figure 1. An analysis of radial velocity data should help elucidate the nature of these variations.

Radial velocities and secondary mass function
All spectra were obtained at high signal-to-noise ratio and revealed well defined Balmer and Hei lines. We fitted the Hα (or Hβ in the P85 data) line cores using a Gaussian function. We estimated the errors at 50-60 km s −1 in the low-dispersion data (P83 and P85) and 5-10 km s −1 in the intermediate dispersion data (P89 and SSO). Table 1 lists all radial velocity measurements. The heliocentric velocity is given as a function of the HJD time at mid-exposure. Folding the data on half the photometric frequency proved successful showing that the variations are due to tidal distortion of the bright hot subdwarf primary. The corresponding ephemeris is where T 0 is the time of passage of the hot subdwarf primary at inferior conjunction. Figure 2 shows the radial velocity data folded on the (above) ephemeris. The ASAS data preceded most of our radial velocity data by ∼ 1000 days and the accumulated error in the predicted phase is ∼ 1.7 minutes. In fact, our radial velocity measurements are offset from the predicted radial velocity curve by only 0.0246 phase (or ∼ 1.7 minutes) and firmly link the photometric variations to tidal distortion of the subdwarf primary. The residual in the more precise P89 velocity measurements is 12 km s −1 and is consistent with the estimated errors of 10 km s −1 . The velocity semi-amplitude and mean velocity of the hot subdwarf are γ 1(sdB) = 31.5 ± 1.3 km s −1 , K 1(sdB) = 386.9 ± 1.9 km s −1 , and the corresponding mass function for the unseen companion is f 2 = 0.294 ± 0.004 M ⊙ .
Finally, we generated a phase-corrected, co-added spectrum of the Hα and Heiλ6678 using all ESO P89 spectra. The exposure time limited orbital smearing to 50 km s −1 , or nearly half of the spectral resolution of the P89 data (90 km s −1 ). A rotation velocity comparable or in excess of 100 km s −1 should be readily measurable. where q = M 1 /M 2 , M 1 and R 1 are the mass and radius of the hot subdwarf primary, and M 2 is the mass of the unseen companion. The subdwarf mass is constrained independently but M 2 and i are treated as free parameters. The semi-major axis a follows Kepler's third law. The gravity darkening and limb darkening coefficients were held fixed at values suggested by Warner (1995), e.g., µ g = 1 and µ = 0.36. The subdwarf radius may be obtained from R 1 /R ⊙ = M 1 /(g/g ⊙ ). Adopting log g from Nemeth et al. (2012) and assuming M 1 = 0.48 M ⊙ , we estimated the subdwarf radius R 1 = 0.159 +0.13 −0.07 R ⊙ . We may also set an independent and more reliable constraint on the primary radius using a measurement of the rotation velocity (Fig 3). Fitting synthetic Hα profiles to co-added P89 spectra  and adopting best-fit model parameters (Nemeth et al. 2012) we measured v rot sin i = 164 ± 5 km s −1 . The models were first convolved with a Gaussian function with a full-width half-maximum set at 2Å. Also, the effect of orbital smearing was taken into account by smoothing the models with a box-car function with a full-width of 50 km s −1 . Geier & Heber (2012) found that all subdwarfs from their sample, single or in binaries with a period longer than 1.2 d are slow rotators (v rot sin i 10 km s −1 ). Clearly, CD−30 11223 is exceptional. We verified our solution by convolving with the rotation broadening function a template built using a subdwarf (GALEX J0206+1438) similar to CD−30 11223 and observed with the same set-up. We measured an identical value of v rot sin i = 164 ± 15 km s −1 but with larger error bars possibly because of a slight parameter mismatch between the two stars. Taking observed ellipsoidal variations as evidence that the subdwarf is tidally locked to the primary, the rotation and orbital velocities are related: v rot sin i K 1 = R 1 a 1 , and since a = a 1 + a 2 = a 1 (1 + q) then where q = M 1 /M 2 = a 2 /a 1 . The resulting primary radius is consistent with the gravitational radius (0.159 R ⊙ ) for all probable inclinations ( 50 • ). Figure 4 shows a set of constraints placed on the free parameters M 2 and i. The mass function defines a narrow constraint for a given subdwarf mass: The minimum secondary mass ranges from 0.75 to 0.77 M ⊙ for a subdwarf mass ranging from 0.44 to 0.48 M ⊙ . Combining mass function and tidal constraints at M 1 = 0.48 M ⊙ , the secondary mass ranges from 0.77 to 0.87 M ⊙ implying a total mass of 1.25 − 1.35 M ⊙ . The effective temperature and surface gravity of CD−30 11223 place it squarely amongst the bulk of sdB stars that are best described with the models of Dorman et al. (1993) within a narrow mass range centered on ∼ 0.48 M ⊙ . Moreover, tests showed that constraints placed by ellipsoidal variations and by the mass function are incompatible below M 1 /M ⊙ ≈ 0.44. Therefore, assuming a mass close to 0.48 M ⊙ is probably justified and the possibility that past interactions may have left an extremely low mass remnant (≈ 0.3 M ⊙ , see, e.g., Heber et al. 2003) does not need to be considered.
Uncertainties on the rotation velocity measurements are possibly affected by systematic effects such as over-or under-estimated spectral resolution or orbital smearing. Assuming a rotation velocity as low as 150 km s −1 would increase the upper mass limit to 1.05 M ⊙ and push the inclination as low as 57 • .
The possibility of shallow (∆m ≈ (R 2 /R 1 ) 2 5 mmag) secondary transits cannot be ruled out and tran- The range of primary mass is dictated by evolutionary considerations, such as the evolutionary tracks of Dorman et al. (1993) that generates remnants of 0.48 M ⊙ (Zhang et al. 2009;Fontaine et al. 2012). Next, a relation between M 2 and the inclination is derived from the measured amplitude of the ellipsoidal variations assuming M 1 = 0.48 (thin lines from bottom to top ∆m =0.086, 0.094, 0.102 magnitude and shaded area). The diagram is split between areas with detectable and non-detectable transits.
sit timing would allow us to constrain accurately the inclination and stellar radii. Transits may become observable at an inclination i such that tan (90 • − i) R 1 /a or i 74 • (where a ≈ 0.6 R ⊙ ). It is worth noting that all solutions provide a primary radius smaller than its Roche lobe radius (R L ≈ 0.2R ⊙ ) excluding the possibility of present-day interaction.

SUMMARY AND CONCLUSIONS
We found that the hot subdwarf CD−30 11223 is in a close 70.5 minutes binary with an unseen companion presumed to be a white dwarf. The orbital properties and noted photometric variations suggest the presence of a companion with a minimum mass of 0.75 M ⊙ assuming a minimum mass of 0.44 M ⊙ for the primary. However, the measured parameters of the sdB CD−30 11223 favor a canonical mass close to 0.48 M ⊙ that would correspond to a minimum mass of 0.77 M ⊙ for the secondary. Constraining the radius of the primary using a measurement of its rotation velocity and fitting the amplitude of the ellipsoidal variations set an upper limit of 0.87 M ⊙ for the secondary mass at an inclination i ≈ 68 • and a total mass of 1.35 M ⊙ for the system. However, systematic effects in the measurement of the rotation velocity and the possibility that the primary mass may exceed 0.48 M ⊙ imply that the system mass may exceed the Chandrasekhar mass limit (M Ch ). Observations of transits would help constrain the inclination of the system and the primary radius.
Following Ritter (1986), the system will initiate mass transfer in 25 to 28 Myr with a contact period of 49 minutes. While transferring mass, the system would take the appearance of an AM CVn system that could eventually trigger a type Ia supernova (Yungelson 2008;Solheim 2010). With a total life-time on the EHB that must exceed 100 Myr (Dorman et al. 1993) and the short time remaining before contact ( 30 Myr), it is likely that a substantial amount of helium may remain in the envelope of the EHB star that would then constitute a generous helium donor. Following Nelemans et al. (2001), the outcome depends on the mass of the accreting white dwarf and the amount of helium available from the donor. The minimum mass of the white dwarf ( 0.75 M ⊙ ) exceeds a threshold (0.6 M ⊙ ) below which the system would undergo nova-like, helium-flash episodes. Instead, should enough helium ( 0.3 M ⊙ ) accumulate on the massive white dwarf after the onset of stable mass transfer the system may suffer helium detonation that could light-up the C/O core triggering a type Ia supernova. On the other hand, it is possible that the EHB star could exhaust most of its helium supply before contact and merge with its massive white dwarf companion. Should the total mass of the system exceed M Ch the remnant could also trigger a type Ia supernova. Alternatively, should the system mass prove insufficient, the immediate outcome would be a hot, ultramassive white dwarf such as those sampled by Vennes et al. (1997) in the Extreme Ultraviolet Explorer survey.
We noted that Geier et al. (2012b) independently reported on this system and reached similar conclusions.
S.V. and A.K. acknowledge support from the Grant Agency of the Czech Republic (GAČR P209/10/0967, GAČR P209/12/0217). This work was also supported by the project RVO:67985815 in the Czech Republic. We thank the referee for helping us clarify the evolutionary prospects of the system. Also, we thank A. Williams (Perth Observatory) for his contribution to the project.
Note added in proof: L. Yungelson informed us that, following Fink et al. (2010, A&A, 514, 53), a supernova may be triggered by an accreted mass of helium as low as 0.1 M ⊙ for a 0.8 M ⊙ white dwarf, thereby increasing the odds of such an event taking place in CD−30 11223.