NY Vir Research Articles

Testing the planetary hypothesis of NY Virginis: anticipated change in the eclipse timing trend within the next five years

ABSTRACT Regarding the recent debate about the potential circumbinary exoplanets around NY Virginis (NY Vir), we analysed mid-eclipse timings of NY Vir using archival photometric data as well as our own observations and the ones from TESS. For this purpose, we first modelled the available eclipse light curves of the binary system to determine the masses of the stars. Then, we measured mid-eclipse times by fitting the light-curve model cycle-to-cycle to the light curves from TESS and our observations. By fitting a Newtonian eclipse timing model to the data, which takes both the light-time effect and potential mutual gravitational interactions into account, we derived orbital parameters and masses of the potential circumbinary planets assuming both eccentric and circular orbits. The models without a quadratic term that can model any possible secular trend, converged to comparable results. Dynamical stability tests show that our Newtonian timing solution corresponds to stable orbital configurations for two circumbinary planets with masses ∼2.3 and ∼4.0 MJup in orbits with very low eccentricity. Our analyses show that the addition of quadratic term for modelling the eclipse timing variation may induce the planetary orbits to be eccentric, hence more likely to be unstable. According to our findings, an upward trend in the eclipse timings followed by a downward one within the next five years is expected due to binary motion induced by circumbinary planets.

Eclipse timing variations in post-common envelope binaries: Are they a reliable indicator of circumbinary companions?

ABSTRACT Post-common envelope binary systems evolve when matter is transferred from the primary star at a rate that cannot be accommodated by its secondary companion. A common envelope forms, which is subsequently ejected resulting in a system with a binary period frequently between 2 and 3 h. Where circumbinary companions are predicted, it remains unclear whether they form before or after the common envelope ejection. From observations of eclipse timing variations (ETVs), exoplanet data bases e.g. NASA Exoplanet Archive, list typically a dozen systems with confirmed circumbinary planets. Here, we examine seven of these systems, discuss other possible causes, and consider whether, for these dynamic systems, the ETV methodology is a reliable indicator of planetary companions. The systems selected were those where we could determine precise eclipse timings, free from significant extraneous effects such as pulsations, and present 163 new times of minima permitting us to test existing models. Over 30 circumbinary models have been proposed for these seven systems and note that all, other than the latest model for NY Vir, which remains to be fully tested, fail within a year to accurately predict eclipse times. In examining alternative mechanisms, we find that magnetic effects could contribute significantly in two of the seven systems studied. We conclude that the structure of these dynamic systems, with the extreme temperature differences and small binary separations, is not fully understood and that many factors may contribute to the observed ETVs.

Filling the gap between synchronized and non-synchronized sdBs in short-period sdBV+dM binaries with TESS: TIC 137608661, a new system with a well-defined rotational splitting

ABSTRACT TIC 137608661/TYC 4544-2658-1/FBS 0938+788 is a new sdBV+dM reflection-effect binary discovered by the TESS space mission with an orbital period of 7.21 h. In addition to the orbital frequency and its harmonics, the Fourier transform of TIC 137608661 shows many g-mode pulsation frequencies from the subdwarf B (sdB) star. The amplitude spectrum is particularly simple to interpret as we immediately see several rotational triplets of equally spaced frequencies. The central frequencies of these triplets are equally spaced in period with a mean period spacing of 270.12 s, corresponding to consecutive l = 1 modes. From the mean frequency spacing of 1.25 μHz we derive a rotation period of 4.6 d in the deep layers of the sdB star, significantly longer than the orbital period. Among the handful of sdB+dM binaries for which the sdB rotation was measured through asteroseismology, TIC 137608661 is the non-synchronized system with both the shortest orbital period and the shortest core rotation period. Only NY Vir has a shorter orbital period but it is synchronized. From a spectroscopic follow-up of TIC 137608661 we measure the radial velocities of the sdB star, determine its atmospheric parameters, and estimate the rotation rate at the surface of the star. This measurement allows us to exclude synchronized rotation also in the outer layers and suggests a differential rotation, with the surface rotating faster than the core, as found in few other similar systems. Furthermore, an analysis of the spectral energy distribution of TIC 137608661, together with a comparison between sdB pulsation properties and asteroseismic models, gives us further elements to constrain the system.

New observations of the eclipsing binary system NY Vir and its candidate circumbinary planets

ABSTRACT As a result of various studies, it has been determined that several post-common envelope eclipsing binaries have variations in their orbital periods. These variations are thought to be caused by the existence of additional bodies in the system (hypothetical stars or planets) and/or other physical effects (such as angular momentum loss, magnetic activity) of the binary system. It is also known that the sdB+M eclipsing system NY Vir has shown such variations in the last decade, indicating additional objects and/or other physical effects. In this work, we present 51 new eclipse times for this system, which extend the time-span of it is O − C diagram by about 3 yr, obtained between 2015 and 2021 using two different telescopes in Turkey. The data obtained in the last 3 yr show a new trend in the O − C diagram differently from the predictions of the previous studies. Our model is consistent with the new O − C diagram, which is statistically well fitted with the quadratic term and the additional two planets with masses of M3 = 2.74 MJup and M4 = 5.59 MJup. However, the orbital period variation can also be related to magnetic activity. In order to better understand the mechanism causing the changes in the orbital period, new observation data are needed that will show at least one full cycle of the change in the O − C diagram.

Convection physics and tidal synchronization of the subdwarf binary NY Virginis

Asteroseismological analysis of NY Vir suggests that at least the outer 55 per cent of the star (in radius) rotates as a solid body and is tidally synchronized to the orbit. Detailed calculation of tidal dissipation rates in NY Vir fails to account for this synchronization. Recent observations of He core burning stars suggest that the extent of the convective core may be substantially larger than that predicted with theoretical models. We conduct a parametric investigation of sdB models generated with the Cambridge STARS code to artificially extend the radial extent of the convective core. These models with extended cores still fail to account for the synchronization. Tidal synchronization may be achievable with a non-MLT treatment of convection.

The quest for stable circumbinary companions to post-common envelope sdB eclipsing binaries

Context. Period variations have been detected in a number of eclipsing close compact binary subdwarf B stars (sdBs) and these have often been interpreted as being caused by circumbinary massive planets or brown dwarfs. According to canonical binary models, the majority of sdB systems are produced from low mass stars with degenerate cores where helium is ignited in flashes. Various evolutionary scenarios have been proposed for these stars, but a definite mechanism remains to be established. Equally puzzling is the formation of these putative circumbinary objects which must have formed from the remaining post-common envelope circumbinary disk or survived its evolution.Aim. In this paper we review the eclipse time variations (ETVs) exhibited by seven such systems (EC 10246-2707, HS 0705+6700, HS 2231+2441, J08205+0008, NSVS 07826147, NSVS 14256825, and NY Vir) and explore whether there is conclusive evidence that the ETVs observed over the last two decades can reliably predict the presence of one or more circumbinary bodies.Methods. We report 246 new observations of the seven sdB systems made between 2013 September and 2017 July using a worldwide network of telescopes. We combined our new data with previously published measurements to analyse the ETVs of these systems.Results. Our data show that period variations cannot be modelled simply on the basis of circumbinary objects. This implies that more complex processes may be taking place in these systems. These difficulties are compounded by the secondary star not being spectroscopically visible. From ETVs, it has historically been suggested that five of the seven binary systems reported here had circumbinary objects. Based on our recent observations and analysis, only three systems remain serious contenders. We find agreement with other observers that at least a decade of observations is required to establish reliable ephemerides. With longer observational baselines it is quite conceivable that the data will support the circumbinary object hypothesis of these binary systems. Also, we generally agree with other observers that higher values of (O–C) residuals are found with secondary companions of spectral type M5/6 (or possibly earlier as a result of an Applegate type mechanism).

Orbital Period Variations in the NY Vir System, Revisited in the Light of New Data

Abstract NY Virginis is an eclipsing binary system with a subdwarf B primary and an M type dwarf secondary. Recent studies (Qian et al. 2012; Lee et al. 2014) suggested the presence of two circumbinary planets with a few Jovian masses within the system. Lee et al. (2014) examined the orbital stabilities of the suggested planets, using the best-fit parameters derived from their eclipse timing variation analysis. They found that the outer companion should be ejected from the system in about 800 000 years. An observational report from Pulley et al. (2016) pointed out that the recent mideclipse times of the binary deviate significantly from the models suggested by Lee et al. (2014). In fact, variations in the orbital period of the system had already been recognized by many authors, but the parameters of these variations vary significantly as new data accumulate. Here, we analyze the eclipse timing variations of the NY Vir system, using new mid-eclipse times that we have obtained together with earlier published measurements in order to understand the nature of the system and constrain its parameters.

The pulsating sdB+M eclipsing system NY Virginis and its circumbinary planets

We searched for circumbinary planets orbiting NY Vir in historical eclipse times including our long-term CCD data. Sixty-eight times of minimum light with accuracies better than 10 s were used for the ephemeris computations. The best fit to those timings indicated that the orbital period of NY Vir has varied due to a combination of two sinusoids with periods of $P_3$=8.2 yr and $P_4$=27.0 yr and semi-amplitudes of $K_3$=6.9 s and $K_4$=27.3 s, respectively. The periodic variations most likely arise from a pair of light-time effects due to the presence of third and fourth bodies that are gravitationally bound to the eclipsing pair. We have derived the orbital parameters and the minimum masses, $M_3 \sin i_3$ = 2.8 M$\rm_{Jup}$ and $M_4 \sin i_4$ = 4.5 M$\rm_{Jup}$, of both objects. A dynamical analysis suggests that the outer companion is less likely to orbit the binary on a circular orbit. Instead we show that future timing data might push its eccentricity to moderate values for which the system exhibits long-term stability. The results demonstrate that NY Vir is probably a star-planet system, which consists of a very close binary star and two giant planets. The period ratio $P_{3}/P_{4}$ suggests that a long-term gravitational interaction between them would result in capture into a nearly 3:10 mean motion resonance. When the presence of the circumbinary planets is verified and understood more comprehensively, the formation and evolution of this planetary system should be advanced greatly.

Period and amplitude variations in post-common-envelope eclipsing binaries observed with SuperWASP

Period or amplitude variations in eclipsing binaries may reveal the presence of additional massive bodies in the system, such as circumbinary planets. Here, we have studied twelve previously-known eclipsing post-common-envelope binaries for evidence of such light curve variations, on the basis of multi-year observations in the SuperWASP archive. The results for HW Vir provided strong evidence for period changes consistent with those measured by previous studies, and help support a two-planet model for the system. ASAS J102322-3737.0 exhibited plausible evidence for a period increase not previously suggested; while NY Vir, QS Vir and NSVS 14256825 afforded less significant support for period change, providing some confirmation to earlier claims. In other cases, period change was not convincingly observed; for AA Dor and NSVS 07826147, previous findings of constant period were confirmed. This study allows us to present hundreds of new primary eclipse timings for these systems, and further demonstrates the value of wide-field high-cadence surveys like SuperWASP for the investigation of variable stars.

CIRCUMBINARY PLANETS ORBITING THE RAPIDLY PULSATING SUBDWARF B-TYPE BINARY NY Vir

We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.

Photometric observations and orbital period variations of HS 0705 + 6700 and NY Vir

We present photometric observations of two post-common-envelope stars, NY Vir (=PG 1336-018) and HS 0705 + 6700. The V band CCD observation of NY Vir was performed by a 40 cm telescope at Ege University Observatory and the R band observations of HS 0705 + 6700 were performed by 100 cm telescope at TÜBİTAK National Observatory. The new light curves were analyzed by the WD code and the physical parameters of stars were determined. We obtained new mid-eclipse timings for HS 0705 + 6700 and combined them with those previously published data. The analysis of the O–C residuals yields a period of about 8.06 ± 0.28 yr and an amplitude of 98.5 s for the system HS 0705 + 6700, which is attributed to the third star physically bounded to the evolved eclipsing pair. A mass function of 1.2 × 10 −4 M ⊙ for the third star is obtained. The existence of a third star is also confirmed by the light curve analysis, indicating light contribution of about 0.043 at phase 0.25 in R-bandpass of the eclipsing pair. Using mass-luminosity relationship of the low mass stars we estimate a mass of 0.12 M ⊙ with an orbital inclination of about 20°. The O–C residuals obtained for the system NY Vir were represented by a downward parabola which indicates orbital period decrease in the system. Using the coefficient of quadratic term we calculate a rate of orbital period decrease of about dP/ dt = −4.09 × 10 −8days yr −1. The period decrease we have measured in NY Vir may be explained by angular momentum loss from the binary system.

The orbital periods of AA Dor and NY Vir

New timings of eclipses made between 2000 and 2010 are presented for two binary systems with hot subdwarf primary stars. In the case of AA Dor, an sdOB star with a very cool secondary, the period is found to be constant at a level of about 10 -14 d per orbit. In the case of NY Vir, a rapidly pulsating sdBV r with a cool companion, the period is discovered to be decreasing at a rate of -11.2 x 10 -13 d per orbit.

The binary properties of the pulsating subdwarf B eclipsing binary PG 1336-018 (NY Virginis)

We present an unbiased orbit solution and mass determination of the components of the eclipsing binary NY Vir as a critical test for the formation scenarios of subdwarf B stars. We obtained high-resolution time series VLT/UVES spectra and high-speed multicolour VLT/ULTRACAM photometric observations of NY Vir, a rapidly pulsating subdwarf B star in a short period eclipsing binary. Combining the radial velocity curve obtained from the VLT/UVES spectra with the VLT/ULTRACAM multicolour lightcurves, we determined numerical orbital solutions for this eclipsing binary. Due to the large number of free parameters and their strong correlations, no unique solution could be found, only families of solutions. We present three solutions of equal statistical significance, two of which are compatible with the primary having gone through a core He-flash and a common-envelope phase described by the $\alpha$-formalism. These two models have an sdB primary of 0.466 \msol and 0.389 \msol, respectively. Finally, we report the detection of the Rossiter-McLaughlin effect for NY Vir.

A Whole Earth Telescope campaign on the pulsating subdwarf B binary system PG 1336−018 (NY Vir)

We present results from a multisite (‘Whole Earth Telescope’) photometric campaign on PG 1336−018, the close eclipsing binary system containing a pulsating subdwarf B (sdB) star. The main part of the campaign (1999 April) resulted in ∼172 h of observations, representing