We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

Recent and upcoming stabilized spectrographs are pushing the frontier for Doppler spectroscopy to detect and characterize low-mass planets. Specifications for these instruments are so impressive that intrinsic stellar variability is expected to limit their Doppler precision for most target stars (Fischer et al. 2016). To realize their full potential, astronomers must develop new strategies for distinguishing true Doppler shifts from intrinsic stellar variability. Stellar variability due to star spots, faculae and other rotationally-linked variability are particularly concerning, as the stellar rotation period is often included in the range of potential planet orbital periods. To robustly detect and accurately characterize low-mass planets via Doppler planet surveys, the exoplanet community must develop statistical models capable of jointly modeling planetary perturbations and intrinsic stellar variability. Towards this effort, this note presents simulations of extremely high resolution, solar-like spectra created with SOAP 2.0 (arXiv:1409.3594) that includes multiple evolving star spots. We anticipate this data set will contribute to future studies developing, testing, and comparing statistical methods for measuring physical radial velocities amid contamination by stellar variability.

A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total sample of 183 Swift GRBs with afterglow plateaus and known redshifts. We added the most recent GRBs to the previous 'gold sample' (now including 45 GRBs) and obtained a relation plane with intrinsic scatter compatible within one sigma with the previous result. We compared several GRB categories, such as short with extended emission, X-ray Flashes, GRBs associated with SNe, long-duration GRBs, and the gold sample, composed only by GRBs with light curves with good data coverage and relatively flat plateaus and evaluated their relation planes. We found that they are not statistically different from the fundamental plane derived from the gold sample and that the fundamental plane still has the smallest scatter. We compared the jet opening angles tabulated in literature with the angles derived using the Eiso-Egamma relation of the method in Pescalli et al. (2015) and calculated the relation plane for a sample of long GRBs accounting for the different jet opening angles. We observed that this correction does not significantly reduce the scatter. In an extended analysis, we found that the fundamental plane is independent from several prompt and afterglow parameters.