Search

We use data from the Transiting Exoplanet Survey Satellite (TESS) to search for, and set limits on, optical to near-infrared photometric variability of the well-vetted, candidate James Webb Space Telescope (JWST) spectrophotometric standards. Our search of 37 of these candidate standards has revealed measurable periodic variability in 15 stars. The majority of those show variability that is less than half a percent; however, four stars are observed to vary photometrically, from minimum to maximum flux, by more than 1% (the G dwarf HD 38949 and three fainter A dwarfs). Variability of this size would likely impact the error budget in the spectrophotometric calibration of the science instruments aboard JWST. For the 22 candidate standards with no detected variability, we report upper limits on the observed changes in flux. Despite some systematic noise, all stars brighter than 12 magnitude in the TESS band show a 3 sigma upper limit on the total change in brightness of less than half a percent on time scales between an hour and multiple weeks, empirically establishing their suitability as spectrophotometric standards. We further discuss the value and limits of high-cadence, high-precision photometric monitoring with TESS as a tool to vet the suitability of stars to act as spectrophotometric standards., Comment: 17 pages, 5 figures, 3 tables, accepted for publication by AAS Journals

To investigate the asymptotic giant branch (AGB) population in the Galactic halo, we search for pulsating AGB stars at a heliocentric distance over 50kpc. Our research is based on the Catalina Southern Survey (CSS) catalogue of variables, comprising 1286 long-period variables (LPVs) with declination less than -20deg. We first focus on the 77 stars in the cap abs(b)>30deg for which spectral M-type or C-type classification can be derived from Hamburg-ESO objective prism spectra. Most of these are oxygen-rich (M-type) and very few are carbon rich. The periods are in the range 100-500 days, and CSS amplitudes are up to 3 mag. In this small sample, no halo AGB star is fainter than Ks0=12.5. This may be due to the scarcity of AGBs in the outer halo, or insufficient instrumental depth. Leaving aside spectral information, we then searched for even fainter pulsators (Ks>12.5) in the entire CSS catalogue. Gaia astrometry makes it possible to identify some contaminants. Our final result is the identification of ten candidate distant LPVs. If these ten stars obey the fundamental mode K-band period luminosity relation used for Miras and small-amplitude Miras, their distances are between 50 and 120kpc from the Sun. In a diagram showing distance versus Gaia tangential velocity, these ten stars have positions similar to those of other objects in the halo, such as globular clusters and dwarf galaxies. We also detect some underluminous AGBs that deserve further study. A detailed catalogue of the 77 high-latitude M or C stars will be made available at the CDS., 9 pages, 7 figures, accepted for publication in A&A (29 April 2021)

Heartbeat stars are a class of eccentric binary stars with short-period orbits and characteristic "heartbeat" signals in their light curves at periastron, caused primarily by tidal distortion. In many heartbeat stars, tidally excited oscillations can be observed throughout the orbit, with frequencies at exact integer multiples of the orbital frequency. Here, we characterize the tidally excited oscillations in the heartbeat stars KIC 6117415, KIC 11494130, and KIC 5790807. Using Kepler light curves and radial velocity measurements, we first model the heartbeat stars using the binary modeling software ELLC, including gravity darkening, limb darkening, Doppler boosting, and reflection. We then conduct a frequency analysis to determine the amplitudes and frequencies of the tidally excited oscillations. Finally, we apply tidal theories to stellar structure models of each system to determine whether chance resonances can be responsible for the observed tidally excited oscillations, or whether a resonance locking process is at work. We find that resonance locking is likely occurring in KIC 11494130, but not in KIC 6117415 or KIC 5790807., Comment: 20 pages, 18 figures, accepted to ApJ

PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself., accepted for publication in ApJS

The Kepler mission has provided a treasure trove of eclipsing binaries (EBs), observed at extremely high photometric precision, nearly continuously for several years. We are carrying out a survey of ~100 of these EBs to derive dynamical masses and radii with precisions of 3% or better. We use multiplexed near-infrared H band spectroscopy from the SDSS-III and -IV APOGEE instrument and optical spectroscopy from the HET High-Resolution Spectrograph to derive double-lined spectroscopic orbits and dynamical mass-ratios for the EB sample. This information is combined with Kepler photometry to derive orbital inclination, dynamical masses of the components, radii and temperatures. These measurements are directly applicable for benchmarking stellar models that are integrating the next generation of improvements, such as the magnetic suppression of convection efficiency, updated opacity tables, and fine-tuned equations of state. We selected our EB sample to include systems with low-mass ($M \lt 0.8 \; M_{\odot}$) primary or secondary components, many expected to populate the relatively sparse parameter space below $0.5 \; M_{\odot}$. In this paper, we describe our EB sample and the analysis techniques we are utilizing, and present masses and radii for two systems that inhabit particularly underpopulated regions of mass-radius-period space: KIC 2445134 and KIC 3003991. Our analysis of KIC 2445134 ($q=0.411 \pm 0.001$) yields masses and radii of $M_A = 1.29 \pm 0.03 \; M_{\odot}$, $M_B = 0.53 \pm 0.01 \; M_{\odot}$, $R_A = 1.42 \pm 0.01 \; R_{\odot}$, $R_B = 0.510 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A = 0.635 \pm 0.001$; our analysis of KIC 3003991 ($q=0.298 \pm 0.006$) yields $M_A = 0.74 \pm \; 0.04 M_{\odot}$, $M_B = 0.222 \pm \; 0.007 M_{\odot}$, $R_A = 0.84 \pm 0.01 \; R_{\odot}$, $R_B = 0.250 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A= 0.662 \pm 0.001$., 16 pages, 10 figures, Published 2019 October 1 in ApJ

The purpose of this work is to extend a sample of accurately modeled, benchmark-grade eclipsing binaries (EBs) with accurately determined masses and radii. We select four “well-behaved” Kepler binaries, KIC 2306740, KIC 4076952, KIC 5193386 and KIC 5288543, each with at least eight double-lined spectra from the Apache Point Observatory Galactic Evolution Experiment instrument that is part of the Sloan Digital Sky Surveys III and IV, and from the Hobby–Eberly High Resolution Spectrograph. We obtain masses and radii with uncertainties of 2.5% or less for all four systems. Three of these systems have orbital periods longer than 9 days, and thus populate an undersampled region of the parameter space for extremely well-characterized detached EBs. We compare the derived masses and radii against mesa mist isochrones to determine the ages of the systems. All systems were found to be coeval, showing that the results are consistent across mesa mist and phoebe.

Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on the extraction of the seismic properties of the pulsating component(s). In this paper, we present an estimation of the impact of a contaminating photometric signal with a distinct background profile on the global asteroseismic parameter $\nu_{\text{max}}$ through the analysis of synthetic red-giant binary light curves. We generate the pulsational and granulation parameters for single red giants with different masses, radii and effective temperatures from theoretical scaling relations and use them to simulate single red-giant light curves with the characteristics of Kepler long-cadence photometric data. These are subsequently blended together according to their light ratio to generate binary red-giant light curves of various configurations. We then perform a differential analysis to characterise the systematic effects of binarity on the extraction of $\nu_{\text{max}}$. We find that the $\nu_{\text{max}}$ extraction for red-giant power spectra featuring overlapping power excesses is unreliable if unconstrained priors are used. Outside of this scenario, we obtain results that are nearly identical to single-star case. We conclude that i) the photometric effects of binarity on the extraction of $\nu_{\text{max}}$ are largely negligible as long as the power excesses of the individual components do not overlap, and that ii) there is minimal advantage to using more than two super-Lorentzian components to model the granulation signal of a binary red-giant., Comment: 17 pages, 17 figures, accepted for publication in A&A

Tidal forces in eccentric binary stars known as heartbeat stars excite detectable oscillations that shed light on the processes of tidal synchronization and circularization. We examine the pulsation phases of tidally excited oscillations (TEOs) in heartbeat binary systems. The target list includes four published heartbeat binaries and four additional systems observed by {\it Kepler}. To the first order, the pulsation phases of TEOs can be explained by the geometric effect of the dominant $l=2$, $m=0$, or $\pm 2$ modes assuming pulsations are adiabatic. We found that this simple theoretical interpretation can account for more than half of the systems on the list, assuming their spin and orbit axes are aligned. We do find significant deviations from the adiabatic predictions for some other systems, especially for the misaligned binary KIC 8164262. The deviations can potentially help to probe the non-adiabaticity of pulsation modes as well as resonances in the tidal forcing., Comment: ApJ accepted

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu., Comment: 52 pages, 14 figures, aastex

We present the characterization of KIC 4142768, an eclipsing binary with two evolved A-type stars in an eccentric orbit with a period of 14 days. We measure the fundamental parameters of the two components ($M_1=2.05M_{\odot}, R_1=2.96R_{\odot}$ and $M_2=2.05M_{\odot}, R_2=2.51R_{\odot}$) by combining {\it Kepler} photometry and spectra from {\it Keck} HIRES. The measured surface rotation rates are only one-fifth of the pseudo-synchronous rate of the eccentric orbit. Fourier spectrum of the light curve reveals hybrid pulsations of $\delta$ Scuti and $\gamma$ Doradus type, with pulsation frequencies at about $15-18$ day$^{-1}$ for p modes and about $0.2-1.2$ day$^{-1}$ for low-frequency g modes. Some of the g modes are exact orbital harmonics and are likely tidally excited. Their pulsation amplitudes and phases both agree with predictions from the linear tidal theory for $l=2, m=2$ prograde modes. We examine the period spacing patterns in the free oscillating g modes and identify them mostly as prograde sectoral dipole modes. The unstable frequency range and frequency spacing of p modes and the inferred asymptotic g-mode period spacings both agree with the stellar model for the primary star evolved to a late stage of the main sequence. The inferred rotation rate of the convective core boundary is very slow, similar to the small surface rotation rate inferred from the spectroscopy. The measured surface and near-core rotation rates provide constraints for testing the mechanism of angular momentum transfer and tidal synchronization in evolved eccentric binary star systems., Comment: ApJ, accepted

Eccentric binaries known as heartbeat stars experience strong dynamical tides as the stars pass through periastron, providing a laboratory to study tidal interactions. We measure the rotation periods of 24 heartbeat systems, using the Kepler light curves to identify rotation peaks in the Fourier transform. Where possible, we compare the rotation period to the pseudosynchronization period derived by Hut 1981. Few of our heartbeat stars are pseudosynchronized with the orbital period. For four systems, we were able to identify two sets of rotation peaks, which we interpret as the rotation from both stars in the binary. The majority of the systems have a rotation period that is approximately 3/2 times the pseudosynchronization period predicted by Hut 1981, suggesting that other physical mechanisms influence the stars' rotation, or that stars typically reach tidal spin equilibrium at a rotation period slightly longer than predicted., Comment: 9 pages, 4 figures, 1 tables

Heartbeat stars (HB stars) are a class of eccentric binary stars with close periastron passages. The characteristic photometric HB signal evident in their light curves is produced by a combination of tidal distortion, heating, and Doppler boosting near orbital periastron. Many HB stars continue to oscillate after periastron and along the entire orbit, indicative of the tidal excitation of oscillation modes within one or both stars. These systems are among the most eccentric binaries known, and they constitute astrophysical laboratories for the study of tidal effects. We have undertaken a radial velocity (RV) monitoring campaign of Kepler HB stars in order to measure their orbits. We present our first results here, including a sample of 21 Kepler HB systems, where for 19 of them we obtained the Keplerian orbit and for 3 other systems we did not detect a statistically significant RV variability. Results presented here are based on 218 spectra obtained with the Keck/HIRES spectrograph during the 2015 Kepler observing season, and they have allowed us to obtain the largest sample of HB stars with orbits measured using a single instrument, which roughly doubles the number of HB stars with an RV measured orbit. The 19 systems measured here have orbital periods from 7 to 90 d and eccentricities from 0.2 to 0.9. We show that HB stars draw the upper envelope of the eccentricity - period distribution. Therefore, HB stars likely represent a population of stars currently undergoing high eccentricity migration via tidal orbital circularization, and they will allow for new tests of high eccentricity migration theories., Accepted to ApJ. Results are reported at http://www.gps.caltech.edu/~shporer/heartbeatstars/index.html

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.

Binary systems where the axis of rotation (spin) of one or both components is tilted w.r.t. the axis of revolution are called misaligned binary systems. The angle of misalignment, obliquity, has been measured for a handful of stars and extrasolar planets to date. Here we present a mathematical framework for a complete and rigorous treatment of misalignment and introduce an extension to the public PHOEBE code that implements this framework. We discuss misalignment for the Roche geometry and introduce methods for computing stellar shapes, equilibrium (generalized Lagrange) points of the potential and minimal requirements for lobe existence. Efficient parametrization of misalignment is proposed in the plane-of-sky coordinates and implementation details in PHOEBE are given alongside the proof-of-concept toy model, comparison with a known misaligned binary DI Her, and comparison with a misaligned planetary system Kepler-13. We provide important mathematical details of the model in the Appendix. This paper accompanies the release of PHOEBE 2.1, which will be available soon from its website http://phoebe-project.org., 20 pages, 10 figures, published in ApJS; release of PHOEBE 2.1 is available at http://phoebe-project.org

Heartbeat stars are eccentric (e > 0.2) ellipsoidal variables whose light curves resemble a cardiogram. We present the observations and corresponding model of KIC 3749404, a highly eccentric (e = 0.66), short period (P = 20.3 d) heartbeat star with tidally induced pulsations. A binary star model was created using PHOEBE, which we modified to include tidally induced pulsations and Doppler boosting. The morphology of the photometric periastron variation (heartbeat) depends strongly on the eccentricity, inclination and argument of periastron. We show that the inclusion of tidally induced pulsations in the model significantly changes the parameter values, specifically the inclination and those parameters dependent on it. Furthermore, we determine the rate of apsidal advance by modelling the periastron variation at the beginning and end of the 4-yr Kepler data set and dividing by the elapsed time. We compare the model with the theoretical expectations for classical and general relativistic apsidal motion and find the observed rate to be two orders of magnitude greater than the theoretical rate. We find that the observed rate cannot be explained by tidally induced pulsations alone and consequently hypothesise the presence of a third body in the system., 16 pages, 8 figures

The unparalleled photometric data obtained by NASA’s Kepler Space Telescope has led to improved understanding of red-giant stars and binary stars. We discuss the characterization of known eccentric system, containing a solar-like oscillating red-giant primary component. We also report several new binary systems that are candidates for hosting an oscillating companion. A powerful approach to study binary stars is to combine asteroseimic techniques with light curve fitting. Seismology allows us to deduce the properties of red giants. In addition, by modeling the ellipsoidal modulations we can constrain the parameters of the binary system. An valuable independent source are ground-bases, high-resolution spectrographs.

KIC 8560861 (HD 183648) is a marginally eccentric (e=0.05) eclipsing binary with an orbital period of P_orb=31.973d, exhibiting mmag amplitude pulsations on time scales of a few days. We present the results of the complex analysis of high and medium-resolution spectroscopic data and Kepler Q0 -- Q16 long cadence photometry. The iterative combination of spectral disentangling, atmospheric analysis, radial velocity and eclipse timing variation studies, separation of pulsational features of the light curve, and binary light curve analysis led to the accurate determination of the fundamental stellar parameters. We found that the binary is composed of two main sequence stars with an age of 0.9\+-0.2 Gyr, having masses, radii and temperatures of M_1=1.93+-0.12 M_sun, R_1=3.30+-0.07 R_sun, T_eff1=7650+-100 K for the primary, and M_2=1.06+-0.08 M_sun, R_2=1.11+-0.03 R_sun, T_eff2=6450+-100 K for the secondary. After subtracting the binary model, we found three independent frequencies, two of which are separated by twice the orbital frequency. We also found an enigmatic half orbital period sinusoidal variation that we attribute to an anomalous ellipsoidal effect. Both of these observations indicate that tidal effects are strongly influencing the luminosity variations of HD 183648. The analysis of the eclipse timing variations revealed both a parabolic trend, and apsidal motion with a period of (P_apse)_obs=10,400+-3,000 y, which is three times faster than what is theoretically expected. These findings might indicate the presence of a distant, unseen companion., 16 pages, 12 figures; accepted for publication in MNRAS

KIC 8569819 is an eclipsing binary star with an early F primary and G secondary in a 20.85-d eccentric orbit. The primary is a delta Sct - gamma Dor star pulsating in both p modes and g modes. Using 4 years of Kepler Mission photometric data we independently model the light curve using the traditional technique with the modelling code PHOEBE, and we study the orbital characteristics using the new frequency modulation (FM) technique. We show that both methods provide the equivalent orbital period, eccentricity and argument of periastron, thus illustrating and validating the FM technique. In the amplitude spectrum of the p-mode pulsations we also discovered an FM signal compatible with a third body in the system, a low mass M dwarf in an 861-d orbit around the primary pair. However, the eclipses show no timing variations, indicating that the FM signal is a consequence of the intrinsic change in pulsation frequency, thus providing a cautionary tale. Our analysis shows the potential of the FM technique using {\it Kepler} data, and we discuss the prospects to detect planets and brown dwarfs in Kepler data for A and F stars even in the absence of transits and with no spectroscopic radial velocity curves. This opens the possibility of finding planets orbiting hotter stars that cannot be found by traditional techniques., Comment: MNRAS in press

Context. The unparalleled photometric data obtained by NASA's Kepler Space Telescope has led to improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries that exhibit ellipsoidal modulations have been detected with Kepler. Aims. We aim to study the properties of eccentric binary systems containing a red giant star and to derive the parameters of the primary giant component. Methods. We applied asteroseismic techniques to determine the masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques were applied to extract the parameters of the system and its primary component. Stellar evolution and its effects on the evolution of the binary system were studied from theoretical models. Results. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440 days. The results of our ongoing radial velocity monitoring programme with the Hermes spectrograph reveal an eccentricity range of e = 0.2 to 0.76. As a case study we present a detailed analysis of KIC 5006817, whose rich oscillation spectrum allows for detailed seismic analysis. From seismology we constrain the rotational period of the envelope to be at least 165 d, which is roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300 ppm in the light curve. Fixing the mass and radius to the asteroseismically determined values, we find from our binary modelling a value of the gravity darkening exponent that is significantly larger than expected. Through binary modelling, we determine the mass of the secondary component to be 0.29 ± 0.03 Mmiddot; Conclusions. For KIC 5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2σ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars. © 2014 ESO., SCOPUS: ar.j, info:eu-repo/semantics/published

We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the delta Scuti and gamma Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software PHOEBE to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98+/-0.07 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do in fact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes., Comment: 19 pages, 13 figures

The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs; http://phoebe-project.org) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity., 60 pages, 15 figures, published in ApJS; accompanied by the release of PHOEBE 2.0 on http://phoebe-project.org

With the advent of high precision photometry from satellites such as Kepler and CoRoT, a whole new layer of interesting and astounding astronomical objects has been revealed: heartbeat stars are an example of such objects. Heartbeat stars are eccen- tric ellipsoidal variables that undergo strong tidal interactions when the stars are almost in contact at the time of closest approach. These interactions deform of the stars and cause a notable light curve variation in the form of a tidal pulse. A subset of these objects (∼20%) show prominent tidally induced pulsations: pulsations forced by the binary orbit. We now have a fully functional code that models binary star features (using phoebe) and stellar pulsations simultaneously, enabling a complete and accurate heartbeat star model to be determined. In this paper we show the results of our new code, which uses emcee, a variant of mcmc, to generate a full set of stellar parameters. We further highlight the interesting features of KIC 8164262, including its tidally induced pulsations and resonantly locked pulsations.

KIC 4544587 is an eclipsing binary star with clear signs of apsidal motion and indications of tidal resonance. The primary component is an early A-type δ Scuti variable, with a temperature of 8270±250 K, whilst the secondary component is an early G-type main sequence star with a temperature of 6500±310 K. The orbital period of this system is 2.18911(1) d, with the light curve demonstrating a hump after secondary minimum due to distortion and reflection. The frequency spectrum of the residual data (the original data with the binary characteristics removed) contains both pressure (p) and gravity (g) modes. Eight of the g modes are precise multiples of the orbital frequency, to an accuracy greater than 3 σ. This is a signature of resonant excitation.

We present the analysis of KIC 8164262, a heartbeat star with a high-amplitude (~1 mmag), tidally resonant pulsation (a mode in resonance with the orbit) at 229 times the orbital frequency and a plethora of tidally induced g-mode pulsations (modes excited by the orbit). The analysis combines Kepler light curves with follow-up spectroscopic data from the Keck telescope, KPNO (Kitt Peak National Observatory) 4-m Mayal telescope and the 2.7-m telescope at the McDonald observatory. We apply the binary modelling software, PHOEBE, to the Kepler light curve and radial velocity data to determine a detailed binary star model that includes the prominent pulsation and Doppler boosting, alongside the usual attributes of a binary star model (including tidal distortion and reflection). The results show that the system contains a slightly evolved F star with an M secondary companion in a highly eccentric orbit (e = 0.886). We use the results of the binary star model in a companion paper (Fuller et al., 2017) where we show that the prominent pulsation can be explained by a tidally excited oscillation mode held near resonance by a resonance locking mechanism., Comment: 13 pages, 9 figures, 4 tables