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.

The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated. Although efforts have been made in the past to reconcile this dispute using asteroseismology, results were inconclusive. In an attempt to resolve this controversy, we study four more evolved planet-hosting stars in this paper using asteroseismology, and we revisit previous results to make an informed study of the whole ensemble in a self-consistent way. For the four new stars, we measure their masses by locating their characteristic oscillation frequency, $\mathrm{\nu}_{\mathrm{max}}$, from their radial velocity time series observed by SONG. For two stars, we are also able to measure the large frequency separation, $\mathrm{\Delta\nu}$, helped by extended SONG single-site and dual-site observations and new TESS observations. We establish the robustness of the $\mathrm{\nu}_{\mathrm{max}}$-only-based results by determining the stellar mass from $\mathrm{\Delta\nu}$, and from both $\mathrm{\Delta\nu}$ and $\mathrm{\nu}_{\mathrm{max}}$. We then compare the seismic masses of the full ensemble of 16 stars with the spectroscopic masses from three different literature sources. We find an offset between the seismic and spectroscopic mass scales that is mass-dependent, suggesting that the previously claimed overestimation of spectroscopic masses only affects stars more massive than about 1.6 M$_\mathrm{\odot}$.

Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young.

Measuring phase modulation in pulsating stars has proved to be a highly successful way of finding binary systems. The class of pulsating main-sequence A and F variables known as delta Scuti stars are particularly good targets for this, and the \textit{Kepler} sample of these has been almost fully exploited. However, some \textit{Kepler} $\delta$ Scuti stars have incorrect temperatures in stellar properties catalogues, and were missed in previous analyses. We used an automated pulsation classification algorithm to find 93 new $\delta$ Scuti pulsators among tens of thousands of F-type stars, which we then searched for phase modulation attributable to binarity. We discovered 10 new binary systems and calculated their orbital parameters, which we compared with those of binaries previously discovered in the same way. The results suggest that some of the new companions may be white dwarfs.

Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the {\delta} Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for {\delta} Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to under-sampling near periapsis. We introduce a new approach that simultaneously forward-models each sample in the light curve and show that this method improves upon current sensitivity limits - especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around {\delta} Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the SNR of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler {\delta} Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.

We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial-velocity observations with the Minerva-Australis telescope array. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of $M_P=0.134^{+0.023}_{-0.022}$$\rm{M_J}$ ($42.6^{+7.3}_{-7.0}$$\rm{M_{\oplus}}$), a radius of $R_P=0.626^{+0.013}_{-0.012}$$\rm{R_J}$ ($7.02^{+0.15}_{-0.13}$$\rm{R_{\oplus}}$), and an orbit with eccentricity $0.242^{+0.040}_{-0.065}$ and period $18.38827\pm0.00072$$\rm{days}$. TOI-257b orbits a bright ($\mathrm{V}=7.570$mag) somewhat evolved late F-type star with $M_*=1.390\pm0.046$$\rm{M_{\odot}}$, $R_*=1.888\pm0.033$$\rm{R_{\odot}}$, $T_{\rm eff}=6075\pm90$$\rm{K}$, and $v\sin{i}=11.3\pm0.5$km/s. Additionally, we statistically validate a second non-transiting sub-Saturn mass planet on a $\sim71$day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems.

Detached eclipsing binaries allow stellar masses and radii to be measured with unrivalled accuracy. While inspecting light curves obtained with the Transiting Exoplanet Survey Satellite (TESS), we noticed that the A0 III star $\alpha$ Dra shows clear and well-separated primary and secondary eclipses. This star is known to be a single-lined spectroscopic binary, with a period of 51.5 d and an eccentricity of 0.43. The currently available TESS observations cover two 27-d sectors and the light curve shows a primary eclipse (depth 9%) and a secondary eclipse (depth 2%), separated in time by 38.5 days. Additional TESS observations of $\alpha$ Dra will come from TESS Sectors 16, 21 and 22, and we predict that an eclipse will be visible in each of these. With a $V$ magnitude of 3.68, $\alpha$ Dra is one of the brightest known detached eclipsing binaries.

We perform a search for rapidly oscillating Ap stars in the Kepler long-cadence data, where true oscillations above the Nyquist limit of 283.21 {\mu}Hz can be reliably distinguished from aliases as a consequence of the barycentric time corrections applied to the Kepler data. We find evidence for rapid oscillations in six stars: KIC 6631188, KIC 7018170, KIC 10685175, KIC 11031749, KIC 11296437 and KIC 11409673, and identify each star as chemically peculiar through either pre-existing classifications or spectroscopic measurements. For each star, we identify the principal pulsation mode, and are able to observe several additional pulsation modes in KIC 7018170. We find that KIC 7018170 and KIC 11409673 both oscillate above their theoretical acoustic cutoff frequency, whilst KIC 11031749 oscillates at the cutoff frequency within uncertainty. All but KIC 11031749 exhibit strong amplitude modulation consistent with the oblique pulsator model, confirming their mode geometry and periods of rotation.

We present the first results from the Transiting Exoplanet Survey Satellite (TESS) on the rotational and pulsational variability of magnetic chemically peculiar A-type stars. We analyse TESS 2-min cadence data from sectors 1 and 2 on a sample of 83 stars. Five new rapidly oscillating Ap (roAp) stars are announced. One of these pulsates with periods around 4.7 min, making it the shortest period roAp star known to date. Four out of the five new roAp stars are multiperiodic. Three of these, and the singly-periodic one show the presence of rotational mode splitting. Individual frequencies are provided in all cases. In addition, seven previously known roAp stars are analysed. Additional modes of oscillation are found in some stars, while in others we are able to distinguish the true pulsations from possible aliases present in the ground-based data. We find that the pulsation amplitude in the TESS filter is typically a factor 6 smaller than that in the $B$ filter which is usually used for ground-based observations. For four roAp stars we set constraints on the inclination angle and magnetic obliquity, through the application of the oblique pulsator model. We also confirm the absence of roAp-type pulsations down to amplitude limits of 6 and 13 micromag, respectively, in two of the best characterised non-oscillating Ap (noAp) stars. We announce 27 new rotational variables along with their rotation periods, and provide different rotation periods for seven other stars. Finally, we discuss how these results challenge state-of-the-art pulsation models for roAp stars.

We have examined the period-luminosity (P-L) relation for delta Scuti stars using Gaia DR2 parallaxes. We included 228 stars from the catalogue of Rodriguez et al. (2000), as well as 1124 stars observed in the four-year Kepler mission. For each star we considered the dominant pulsation period, and used DR2 parallaxes and extinction corrections to determine absolute V magnitudes. Many stars fall along a sequence in the P-L relation coinciding with fundamental-mode pulsation, while others pulsate in shorter-period overtones. The latter stars tend to have higher effective temperatures, consistent with theoretical calculations. Surprisingly, we find an excess of stars lying on a ridge with periods half that of the fundamental. We suggest this may be due to a 2:1 resonance between the third or fourth overtone and the fundamental mode.

We study the fraction of stars in and around the {\delta} Scuti instability strip that are pulsating, using {\em Gaia} DR2 parallaxes to derive precise luminosities. We classify a sample of over 15,000 {\em Kepler} A and F stars into {\delta} Sct and non-{\delta} Sct stars, paying close attention to variability that could have other origins. We find that 18 per cent of the {\delta} Sct stars have their dominant frequency above the Kepler long-cadence Nyquist frequency (periods < 1 hr), and 30 per cent have some super-Nyquist variability. We analyse the pulsator fraction as a function of effective temperature and luminosity, finding that many stars in the {\delta} Sct instability strip do not pulsate. The pulsator fraction peaks at just over 70 per cent in the middle of the instability strip. The results are insensitive to the amplitude threshold used to identify the pulsators. We define a new empirical instability strip based on the observed pulsator fraction that is systematically hotter than theoretical strips currently in use. The stellar temperatures, luminosities, and pulsation classifications are provided in an online catalogue.