The distortions of absorption line profiles caused by photospheric brightness variations on the surfaces of cool, main-sequence stars can mimic or overwhelm radial velocity (RV) shifts due to the presence of exoplanets. The latest generation of precision RV spectrographs aims to detect velocity amplitudes $\lesssim 10$ cm s$^{-1}$, but requires mitigation of stellar signals. Statistical techniques are being developed to differentiate between Keplerian and activity-related velocity perturbations. Two important challenges, however, are the interpretability of the stellar activity component as RV models become more sophisticated, and ensuring the lowest-amplitude Keplerian signatures are not inadvertently accounted for in flexible models of stellar activity. For the K2V exoplanet host $\epsilon$ Eridani, we separately use ground-based photometry to constrain Gaussian processes for modeling RVs and TESS photometry with a light-curve inversion algorithm to reconstruct the stellar surface. From the reconstructions of TESS photometry, we produce an activity model, which reduces the rms scatter in RVs obtained with EXPRES from 4.72 m s$^{-1}$ to 1.98 m s$^{-1}$. We present a pilot study using the CHARA Array and MIRC-X beam combiner to directly image the starspots seen in the TESS photometry. With the limited phase coverage, our spot detections are marginal with current data but a future dedicated observing campaign should allow for imaging, as well as the stellar inclination and orientation with respect to its debris disk to be definitely determined. This work shows that stellar surface maps obtained with high cadence, time-series photometric and interferometric data can provide the constraints needed to accurately reduce RV scatter.

By controlling instrumental errors to below 10 cm s−1, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian processes (GP) have become a standard tool for modeling correlated noise in radial velocity data sets. While GPs are constrained and motivated by physical properties of the star, in some cases they are still flexible enough to absorb unresolved Keplerian signals. We apply GP regression to EXPRES radial velocity measurements of the 3.5 Gyr old chromospherically active Sun-like star, HD 101501. We obtain tight constraints on the stellar rotation period and the evolution of spot distributions using 28 seasons of ground-based photometry, as well as recent Transiting Exoplanet Survey Satellite data. Light-curve inversion was carried out on both photometry data sets to reveal the spot distribution and spot evolution timescales on the star. We find that the >5 m s−1 rms radial velocity variations in HD 101501 are well modeled with a GP stellar activity model without planets, yielding a residual rms scatter of 45 cm s−1. We carry out simulations, injecting and recovering signals with the GP framework, to demonstrate that high-cadence observations are required to use GPs most efficiently to detect low-mass planets around active stars like HD 101501. Sparse sampling prevents GPs from learning the correlated noise structure and can allow it to absorb prospective Keplerian signals. We quantify the moderate to high-cadence monitoring that provides the necessary information to disentangle photospheric features using GPs and to detect planets around active stars.

From 5.5 months of dual-band optical photometric monitoring at the 1 mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star ζ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the 1.78-d signal previously detected by Coriolis/Solar Mass Ejection Imager, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light-curve inversion algorithm, we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He II λ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II λ4686 line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in ζ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

We analyze extensive ground- and space-based ultraviolet-infrared observations of the variable white dwarf BOKS 53856. Photometry and spectroscopy were obtained from the Hubble Space Telescope, Spitzer Space Telescope, Kepler and Swift satellites, Palomar Observatory, and Kitt Peak National Observatory. Time series data spanning the entire four-year Kepler mission are used to refine the photometric rotation period of the white dwarf to 0.2557292(9) days. The associated periodic variability, with full (minimum to maximum) amplitudes of approximate to 5% and approximate to 10% in the optical and ultraviolet, respectively, is consistent with the presence of several dark spots on the white dwarf surface. A fit to the observed ultraviolet-infrared spectral energy distribution of BOKS 53856 using a grid of non-LTE synthetic spectra gives best values of effective temperature, T-WD = 31,056(19) K, surface gravity, log g = 7.913(5), mass, M-WD = 0.603(2) M-circle dot, radius, R-WD = 0.0142(4) R-circle dot, reddening, Eb-v = 0.0533(1), and distance, d = 277(10) pc. The Hubble ultraviolet spectrum does not show a detectable signature of accreted metals, and the Spitzer infrared photometry does not reveal an excess attributable to a circumstellar dust disk. Nonetheless, the ground-based optical spectra show a very weak, variable Ca II K absorption line that is phased with the Kepler light curve, being strongest at minimum brightness. We discuss BOKS 53856 and similar objects in the context of a localized enhancement of absorption from metals in dark (possibly magnetic) spots on the white dwarf surface.

Nearby active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary σ Geminorum (σ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods.

Infrared interferometry imaging of the old, magnetically active star ζ Andromedae reveals an asymmetric distribution of starspots, unlike the north–south starspot symmetry observed on the Sun, meaning the underlying dynamo mechanisms must be different. As observed on the Sun, strong magnetic fields suppress convection to create localized cooler regions that appear as sunspots or starspots. It is difficult to observe starspots on stars other than the Sun, but here Rachael Roettenbacher and colleagues present infrared interferometry images of the old, magnetically active star ζ Andromedae that reveal starspot distribution in detail. They find an asymmetric distribution, in contrast to the north–south symmetry observed on the Sun, meaning that the underlying dynamo mechanisms in the two bodies must be different. Sunspots are cool areas caused by strong surface magnetic fields that inhibit convection1,2. Moreover, strong magnetic fields can alter the average atmospheric structure3, degrading our ability to measure stellar masses and ages. Stars that are more active than the Sun have more and stronger dark spots than does the Sun, including on the rotational pole4. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal5. This leads to problems in investigating the north–south distribution of starspot active latitudes (those latitudes with more starspot activity); this distribution is a crucial constraint of dynamo theory. Polar spots, whose existence is inferred from Doppler tomography, could plausibly be observational artefacts6. Here we report imaging of the old, magnetically active star ζ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two observation epochs, whereas lower-latitude spot structures in both hemispheres do not persist between observations, revealing global starspot asymmetries. The north–south symmetry of active latitudes observed on the Sun7 is absent on ζ And, which hosts global spot patterns that cannot be produced by solar-type dynamos8.

Using the high-precision photometry from the Kepler space telescope, thousands of stars with stellar and planetary companions have been observed. The characterization of stars with companions is not always straightforward and can be contaminated by systematic and stellar influences on the light curves. Here, through a detailed analysis of starspots and eclipses, we identify KOI-1003 as a new, active RS CVn star-the first identified with data from Kepler. The Kepler light curve of this close binary system exhibits the system's primary transit, secondary eclipse, and starspot evolution of two persistent active longitudes. The near equality of the system's orbital and rotation periods indicates the orbit and primary star's rotation are nearly synchronized (P-orb = 8.360613 +/- 0.000003 days; P-rot similar to 8.23 days). By assuming the secondary star is on the main sequence, we suggest the system consists of a 1.45(-0.19)(+0.11) M-circle dot subgiant primary and a 0.59(-0.04)(+0.03) M-circle dot main-sequence companion. Our work gives a distance of 4400 +/- 600 pc and an age of t = 3.0(+2.0)(-0.5) Gyr, parameters which are discrepant with previous studies that included the star as a member of the open cluster NGC 6791.

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2 m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M-A=1.35 +/- 0.05M., MB=0.99 +/- 0.02M(circle dot), system age = 3.0. -/+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta=0.07 +/- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.

To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory.

Abstract With the precise, nearly-continuous photometry from the Kepler satellite and the sub-milliarcsecond resolving capabilities of the CHARA Array, astronomy is entering a new age for the imaging and understanding of stellar magnetic activity. We present first results from our Guest Observer Program, where 180 single-epoch surface image reconstructions of KIC 5110407 have revealed differential rotation and hints of magnetic activity cycles based on both spot and flare variations. Analysis of our larger, full dataset will establish in unprecedented detail how surface magnetic activity correlates with stellar age and spectral type. In addition to Kepler work, we have harnessed the power of the world's largest infrared interferometer to “directly” image the spotted surfaces of a few of the closest RS CVn systems, allowing a comparison of contemporaneous Doppler and light-curve inversion imaging techniques.