Search

A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P, Comment: 12 pages, 4 figures, published in Nature Astronomy

Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature., Invited review in press at RASTI. Based on the ExoPAG SAG21 report (arXiv:2201.09905v1) and refined via feedback from three reviewers. 75 pages, 30 figures, 5 tables

We present H-band interferometric observations of the red supergiant (RSG) AZ Cyg made with the Michigan Infra-Red Combiner (MIRC) at the six-telescope Center for High Angular Resolution Astronomy (CHARA) Array. The observations span 5 years (2011-2016), offering insight into the short and long-term evolution of surface features on RSGs. Using a spectrum of AZ Cyg obtained with SpeX on the NASA InfraRed Telescope Facility (IRTF) and synthetic spectra calculated from spherical MARCS, spherical PHOENIX, and SAtlas model atmospheres, we derive $T_{\text{eff}}$ is between $3972 K$ and $4000 K$ and $\log~g$ between $-0.50$ and $0.00$, depending on the stellar model used. Using fits to the squared visibility and Gaia parallaxes we measure its average radius $R=911^{+57}_{-50}~R_{\odot}$. Reconstructions of the stellar surface using our model-independent imaging codes SQUEEZE and OITOOLS.jl show a complex surface with small bright features that appear to vary on a timescale of less than one year and larger features that persist for more than one year. 1D power spectra of these images suggest a characteristic size of $0.52-0.69~R_{\star}$ for the larger, long lived features. This is close to the values of $0.51-0.53~R_{\star}$ derived from 3D RHD models of stellar surfaces. We conclude that interferometric imaging of this star is in line with predictions of 3D RHD models but that short-term imaging is needed to more stringently test predictions of convection in RSGs., 31 pages, 8 tables, 14 figures, accepted to ApJ

Aims: We aim to study the spot evolution and differential rotation in the magnetically active cool K-type giant star sigma Gem from broadband photometry and continuous spectroscopic observations that span 150 nights. Methods: We use high-resolution, high signal-to-noise ratio spectra obtained with the Hertzsprung SONG telescope to reconstruct surface (photospheric) temperature maps with Doppler imaging techniques. The 303 observations span 150 nights and allow for a detailed analysis of the spot evolution and surface differential rotation. The Doppler imaging results are compared to simultaneous broadband photometry from the Tennessee State University T3 0.4 m Automated Photometric Telescope. The activity from the stellar chromosphere, which is higher in the stellar atmosphere, is also studied using SONG observations of Balmer H alpha line profiles and correlated with the photospheric activity. Results: The temperature maps obtained during eight consecutive stellar rotations show mainly high-latitude or polar spots, with the main spot concentrations above latitude 45 deg. The spots concentrate around phase 0.25 near the beginning of our observations and around phase 0.75 towards the end. The photometric observations confirm a small jump in spot phases that occurred in February 2016. The cross-correlation of the temperature maps reveals rather strong solar-like differential rotation, giving a relative surface differential rotation coefficient of $\alpha$ = 0.10 +/- 0.02. There is a weak correlation between the locations of starspots and enhanced emission in the chromosphere at some epochs., Comment: 19 pages, 10 figures, accepted for publication in Astronomy & Astrophysics, in the published version the Tables B1 and B2 will only be available via CDS

Long-baseline optical interferometry has provided a method with which stars are not only resolved, but directly imaged. Resolving features on stellar surfaces is becoming routine for observers using the Center for High Angular Resolution Astronomy (CHARA) Array or the Very Large Telescope Interferometer (VLTI), in particular. The CHARA Array presently has six 1-m telescopes at permanent locations with baselines ranging from 34-331 m.1 Currently for imaging, the six CHARA Array telescopes are combined with the Michigan InfraRed Combiner-eXeter (MIRC-X; formerly the Michigan Infrared Combiner, MIRC).2 Working in H- and J -bands, MIRC-X provides an angular resolution of θH ∼ 0.5 mas and θJ ∼ 0.4 mas, respectively. Four telescopes at a time can be combined at VLTI, either the 8.2-m Unit Telescopes (UTs) or the 1.8-m Auxiliary Telescopes (ATs). The ATs can be moved to provide multiple configurations to increase the possible baselines (up to 130 m). The Precision Integrated-Optics Near-Infrared Imaging ExpeRiment (PIONIER),3 which operates in H-band (θH ∼ 1.3 mas), is the current beam combiner frequently used for imaging. In the following sections, the types of stars imaged by long-baseline optical interferometry will be discussed with significant astrophysical results highlighted.

By controlling instrumental errors to below 10 cm/s, 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 datasets. 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 $TESS$ data. Light curve inversion was carried out on both photometry datasets to reveal the spot distribution and spot evolution timescales on the star. We find that the $> 5$ m/s 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. 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., 25 pages, 16 figures, accepted to AJ

Exoplanet discoveries have reached into the realm of terrestrial planets that are becoming the subject of atmospheric studies. One such discovery is LHS 3844b, a 1.3 Earth radius planet in a 0.46 day orbit around an M4.5-5 dwarf star. Follow-up observations indicate that the planet is largely devoid of substantial atmosphere. This lack of significant atmosphere places astrophysical and geophysical constraints on LHS 3844b, primarily the degree of volatile outgassing and the rate of atmosphere erosion. We estimate the age of the host star as $7.8\pm1.6$ Gyrs and find evidence of an active past comparable to Proxima Centauri. We use geodynamical models of volcanic outgassing and atmospheric erosion to show that the apparent lack of atmosphere is consistent with a volatile-poor mantle for LHS 3844b. We show the core is unlikely to host enough C to produce a sufficiently volatile-poor mantle, unless the bulk planet is volatile-poor relative to Earth. While we cannot rule out a giant impact stripping LHS 3844b's atmosphere, we show this mechanism would require significant mantle stripping, potentially leaving LHS 3844b as an Fe-rich "super-Mercury". Atmospheric erosion by smaller impacts is possible, but only if the planet has already begun degassing and is bombarded by $10^3$ impactors of radius 500-1000 km traveling at escape velocity. We discuss formation and migration scenarios that could account for a volatile poor origin, including the potential for an unobserved massive companion planet. A relatively volatile-poor composition of LHS 3844b suggests that the planet formed interior to the system snow-line., Comment: 23 pages, 14 figures, accepted for publication in the Planetary Science Journal

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., Comment: 24 pages, 12 figures, 5 tables, accepted for publication in AJ

Alpha Ophiuchi (Rasalhague) is a nearby rapidly rotating A5IV star which has been imaged by infrared interferometry. $\alpha$ Oph is also part of a known binary system, with a companion semi-major axis of $\sim$430 milli-arcseconds and high eccentricity of 0.92. The binary companion provides the unique opportunity to measure the dynamical mass to compare with the results of rapid rotator evolution models. The lack of data near periastron passage limited the precision of mass measurements in previous work. We add new interferometric data from the MIRC combiner at the CHARA Array as well as new Keck adaptive optics imaging data with NIRC2, including epochs taken near periastron passage. We also obtained new radial velocities of both components at Fairborn Observatory. Our updated combined orbit for the system drastically reduces the errors of the orbital elements, and allows for precise measurement of the primary star mass at the few percent level. Our resulting primary star mass of $2.20\pm0.06$ M$_{\odot}$ agrees well with predictions from imaging results, and matches evolution models with rotation when plotting on an HR diagram. However, to truly distinguish between non-rotating and rotating evolution models for this system we need $\sim$1\% errors on mass, which might be achieved once the distance is known to higher precision in future Gaia releases. We find that the secondary mass of $0.824\pm0.023$ M$_{\odot}$ is slightly under-luminous when compared to stellar evolution models. We show that $\alpha$ Oph is a useful reference source for programs that need $\pm$1 milli-arcsecond astrometry., Comment: 17 pages, 7 figures, accepted for publication in ApJ

We present temperature maps of RS CVn star λ Andromedae, reconstructed from interferometric data acquired in 2010 and 2011 by the Michigan Infra-Red Combiner instrument at the Center for High Angular Resolution Astronomy Array. To constrain the stellar parameters required for this imaging task, we first modeled the star using our GPU-accelerated code SImulation and Modeling Tool for Optical Interferometry. The stellar surface was then imaged using our open source interferometric imaging code ROTational Image Reconstruction, in the process further refining the estimation of stellar parameters. We report that the measured angular diameter is 2.742 ± 0.010 mas with a limb-darkening coefficient of 0.231 ± 0.024. While our images are consistent with those of prior works, we provide updated physical parameters for λ Andromedae (R ⋆ = 7.78 ± 0.05 R ⊙, M ⋆ = 1.24 ± 0.72 M ⊙, and log L/L ⊙ = 1.46 ± 0.04).

The EXPRES Stellar-Signals Project is providing sets of high-fidelity, spectroscopic and photometric observations to enable direct comparisons of various approaches for disentangling stellar signals and true radial velocities (RVs). We will provide all EXPRES RVs, meta data, and activity indicators as well as high-precision photometric data from the Fairborn Automatic Photoelectric Telescopes (APTs) for HD 101501, HD 34411, HD 217014, and HD 10700. Intrinsic stellar variability and the resulting apparent RVs are widely believed to dominate the error budget for extremely precise radial-velocity (EPRV) measurements. Several new methods to disentangle photospheric velocities from Keplerian velocities are being developed throughout the EPRV community. In addition to releasing data sets for testing these methods, the EXPRES Stellar-Signals Project will publish a summary of the current state of the field circa 2020 to guide next steps toward mitigating photospheric velocities in EPRV data. More information can be found on http://exoplanets.astro.yale.edu/science/activity.php.

We aim at detecting and characterizing the main-sequence companions of a sample of known and suspected Galactic binary Cepheids. We used the multi-telescope interferometric combiners MIRC and PIONIER to detect and measure the astrometric positions of the high-contrast companions orbiting 16 bright Galactic Cepheids. We made use of the CANDID algorithm to search for the companions and set detection limits. We also present new high-precision RVs which were used to fit radial pulsation and orbital velocities. We report the detection of the companions orbiting the Cepheids U Aql, BP Cir, and S Mus, possible detections for FF Aql, Y Car, BG Cru, X Sgr, V350 Sgr, and V636 Sco, while no component is detected around U Car, YZ Car, T Mon, R Mus, S Nor, W Sgr and AH Vel. For U Aql and S Mus, we performed a preliminary orbital fit combining astrometry with new high-precision RVs, providing the full set of orbital and pulsation parameters. Assuming the distance, we estimated preliminary masses of M(U Aql) = 4.97+/-0.62Msun and M(S Mus) = 4.63+/-0.99Msun. For YZ Car, W Sgr, V350 Sgr, and V636 Sco, we revised the spectroscopic orbits using new high-precision RVs, while we updated the pulsation parameters for BP Cir, BG Cru, S Nor and AH Vel. Our interferometric observations also provide measurements of the angular diameters, that can be used in a Baade-Wesselink type analysis. We have now several astrometric detections of Cepheid companions. When RVs of the companions will be available, such systems will provide accurate and independent masses and distances. Orbital parallaxes with an accuracy, 29 pages, 11 figures, 23 tables. Accepted for publication in A&A

The physical nature of time variable objects is often inferred from photometric light-curves and spectroscopic variations. Long-baseline optical interferometry has the power to resolve the spatial structure of time variable sources directly in order to measure their physical properties and test the physics of the underlying models. Recent interferometric studies of variable objects include measuring the angular expansion and spatial structure during the early stages of novae outbursts, studying the transits and tidal distortions of the components in eclipsing and interacting binaries, measuring the radial pulsations in Cepheid variables, monitoring changes in the circumstellar discs around rapidly rotating massive stars, and imaging starspots. Future applications include measuring the image size and centroid displacements in gravitational microlensing events, and imaging the transits of exoplanets. Ongoing and upcoming photometric surveys will dramatically increase the number of time-variable objects detected each year, providing many potential targets to observe interferometrically. For short-lived transient events, it is critical for interferometric arrays to have the flexibility to respond rapidly to targets of opportunity and optimize the selection of baselines and beam combiners to provide the necessary resolution and sensitivity to resolve the source as its brightness and size change. We discuss the science opportunities made possible by resolving variable sources using long baseline optical interferometry.

Cepheid stars play a considerable role as extragalactic distances indicators, thanks to the simple empirical relation between their pulsation period and their luminosity. They overlap with that of secondary distance indicators, such as Type Ia supernovae, whose distance scale is tied to Cepheid luminosities. However, the Period-Luminosity (P-L) relation still lacks a calibration to better than 5%. Using an original combination of interferometric astrometry with optical and ultraviolet spectroscopy, we measured the geometrical distance d = 720.35+/-7.84 pc of the 3.33 d period Cepheid V1334 Cyg with an unprecedented accuracy of +/-1 %, providing the most accurate distance for a Cepheid. Placing this star in the P-L diagram provides an independent test of existing period-luminosity relations. We show that the secondary star has a significant impact on the integrated magnitude, particularly at visible wavelengths. Binarity in future high precision calibrations of the P-L relations is not negligible, at least in the short-period regime. Subtracting the companion flux leaves V1334 Cyg in marginal agreement with existing photometric-based P-L relations, indicating either an overall calibration bias or a significant intrinsic dispersion at a few percent level. Our work also enabled us to determine the dynamical masses of both components, M1 = 4.288 +/- 0.133 Msun (Cepheid) and M2 = 4.040 +/- 0.048 Msun (companion), providing the most accurate masses for a Galactic binary Cepheid system., Comment: Accepted for publication in ApJ

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007., Accepted for publication on ApJ; 32 pages, 16 figures, 10 tables

Starspots and flares are indicators of stellar magnetic activity that can both be studied in greater detail by utilizing the long-term, space-based archive of the Kepler satellite. Here, we aim to investigate a subset of the Kepler archive to reveal a connection between the starspots and the stellar flares, in order to provide insight into the overall stellar magnetic field. We use the flare-finding algorithm FLATW'RM in conjunction with a new suite of algorithms that aim to locate the local minima caused by starspot groups. We compare the phase difference between the time of maximum flux of a flare and the time of minimum stellar flux due to a starspot group. The strongest flares do not appear to be correlated to the largest starspot group present, but are also not uniformly distributed in phase with respect to the starspot group. The weaker flares, however, do show an increased occurrence close to the starspot groups., Accepted to ApJ; 11 pages, 6 figures

Combining visual and spectroscopic orbits of binary stars leads to a determination of the full 3D orbit, individual masses, and distance to the system. We present a full analysis of the evolved binary system $\delta$ Delphini using astrometric data from the MIRC and PAVO instruments on the CHARA long-baseline interferometer, 97 new spectra from the Fairborn Observatory, and 87 unpublished spectra from Lick Observatory. We determine the full set of orbital elements for $\delta$ Del, along with masses of $1.78 \pm 0.07$ $M_{\odot}$ and $1.62 \pm 0.07$ $M_{\odot}$ for each component, and a distance of $63.61 \pm 0.89$ pc. These results are important in two contexts: for testing stellar evolution models and defining the detection capabilities for future planet searches. We find that the evolutionary state of this system is puzzling, as our measured flux ratios, radii, and masses imply a $\sim$ 200 Myr age difference between the components using standard stellar evolution models. Possible explanations for this age discrepancy include mass transfer scenarios with a now ejected tertiary companion. For individual measurements taken over a span of 2 years we achieve $2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry., Comment: 28 pages, 14 figures, accepted for publication in ApJ

Archives of long photometric surveys, like the Kepler database, are a gold mine for studying flares. However, identifying them is a complex task; while in the case of single-target observations it can be easily done manually by visual inspection, this is nearly impossible for years-long time series for several thousand targets. Although there exist automated methods for this task, several problems are difficult (or impossible) to overcome with traditional fitting and analysis approaches. We introduce a code for identifying and analyzing flares based on machine learning methods, which are intrinsically adept at handling such data sets. We used the RANSAC (RANdom SAmple Consensus) algorithm to model light curves, as it yields robust fits even in case of several outliers, like flares. The light curve is divided into search windows, approximately in the order of the stellar rotation period. This search window is shifted over the data set, and a voting system is used to keep false positives to a minimum: only those flare candidate points are kept that were identified in several windows as a flare. The code was tested on the K2 observations of the TRAPPIST-1, and on the long cadence data of KIC 1722506. The detected flare events and flare energies are consistent with earlier results from manual inspections., Comment: 6 pages, 5 figures, accepted to A&A

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 $\sigma$ Geminorum ($\sigma$ 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., Comment: Accepted to ApJ, 19 pages, 16 figures, 5 tables

We report new CHARA/MIRC interferometric observations of the Cepheid archetype $\delta$ Cep, which aimed at detecting the newly discovered spectroscopic companion. We reached a maximum dynamic range $\Delta H $ = 6.4, 5.8, and 5.2 mag, respectively within the relative distance to the Cepheid $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. Our observations did not show strong evidence of a companion. We have a marginal detection at $3\sigma$ with a flux ratio of 0.21\%, but nothing convincing as we found other possible probable locations. We ruled out the presence of companion with a spectral type earlier than F0V, A1V and B9V, respectively for the previously cited ranges $r$. From our estimated sensitivity limits and the Cepheid light curve, we derived lower-limit magnitudes in the $H$ band for this possible companion to be $H_\mathrm{comp} > 9.15, 8.31$ and 7.77 mag, respectively for $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. We also found that to be consistent with the predicted orbital period, the companion has to be located at a projected separation $< 24$ mas with a spectral type later than a F0V star., Comment: Accepted for publication in MNRAS

Magnetic fields are essential ingredients of many physical processes in the interiors and envelopes of cool stars. Yet their direct detection and characterisation is notoriously difficult, requiring high-quality observations and advanced analysis techniques. Significant progress has been recently achieved by several types of direct magnetic field studies on the surfaces of cool active stars. In particular, complementary techniques of the field topology mapping with polarisation data and total magnetic flux measurements from intensity spectra have been systematically applied to different classes of active stars leading to interesting and occasionally controversial results. In this paper we summarise the current status of direct magnetic field studies of cool stars, and investigations of surface inhomogeneities caused by the field, based on the material presented at the Cool Stars 19 splinter session., Summary of the splinter session "Surface Magnetism of Cool Stars" at the Cool Stars 19 conference; to be published in Astronomische Nachrichten

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_\mathrm{orb} = 8.360613\pm0.000003$ days; $P_\mathrm{rot} \sim 8.23$ days). By assuming the secondary star is on the main sequence, we suggest the system consists of a $1.45^{+0.11}_{-0.19} \ M_\odot$ subgiant primary and a $0.59^{+0.03}_{-0.04} \ M_\odot$ main-sequence companion. Our work gives a distance of $4400 \pm 600$ pc and an age of $t = 3.0^{-0.5}_{+2.0}$ Gyr, parameters which are discrepant with previous studies that included the star as a member of the open cluster NGC 6791., Comment: 21 pages, 19 figures, accepted to ApJ

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, Spi ...

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.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, 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., Comment: 14 pages, 13 figures, Accepted to ApJ

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively., Comment: Accepted for publication in Astronomy and Astrophysics

The signatures of planets hosted by M dwarfs are more readily detected with transit photometry and radial velocity methods than those of planets around larger stars. Recently, transit photometry was used to discover seven planets orbiting the late-M dwarf TRAPPIST-1. Three of TRAPPIST-1's planets fall in the Habitable Zone, a region where liquid water could exist on the planetary surface given appropriate planetary conditions. We aim to investigate the habitability of the TRAPPIST-1 planets by studying the star's activity and its effect on the planets. We analyze previously-published space- and ground-based light curves and show the photometrically-determined rotation period of TRAPPIST-1 appears to vary over time due to complicated, evolving surface activity. The dramatic changes of the surface of TRAPPIST-1 suggest that rotation periods determined photometrically may not be reliable for this and similarly active stars. While the activity of the star is low, we use the premise of the "cosmic shoreline" to provide evidence that the TRAPPIST-1 environment has potentially led to the erosion of possible planetary atmospheres by extreme ultraviolet stellar emission., Accepted to ApJ, 8 pages, 5 figures, 1 table

We thank Robert Wilson (University of Florida) for providing a custom version of his code to compute images of spotted stellar surfaces and for his help with using it. This work is based upon observations obtained with the Georgia State University (GSU) Center for High Angular Resolution Astronomy (CHARA) array at Mount Wilson Observatory. The CHARA array is supported by the National Science Foundation under grant numbers AST-1211929 and AST-1411654. Institutional support has been provided by the GSU College of Arts and Sciences and the GSU Office of the Vice President for Research and Economic Development. The MIRC instrument at the CHARA array was funded by the University of Michigan. F.B., R.R., and J.D.M. acknowledge support from NSF-AST 1210972 and 1108963. G.T. acknowledges partial support from NSF grant AST-1509375. S.K. acknowledges support from an STFC Rutherford Fellowship (ST/J004030/1) and ERC Starting Grant (grant agreement no. 639889). This work is also based on observations made with the Nordic Optical Telescope (NOT), operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France. This research has made use of the Jean-Marie Mariotti Center SearchCal service13 codeveloped by FIZEAU and LAOG/IPAG and of the CDS astronomical databases SIMBAD and VIZIER.14 This research has made use of the Washington Double Star Catalog, maintained at the U.S. Naval Observatory. We thank Nicholas Elias II for discussions. We thank Dimitri Pourbaix for maintaining and providing access to the SB9 database of RV measurements of spectroscopic binaries.

Sunspots are cool areas caused by strong surface magnetic fields inhibiting convection. Moreover, strong magnetic fields can alter the average atmospheric structure, degrading our ability to measure stellar masses and ages. Stars more active than the Sun have more and stronger dark spots than in the solar case, including on the rotational pole itself. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars than the Sun, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal. This leads to problems in investigating the north-south distribution of starspot active latitudes (those latitudes with more spot activity), which are crucial constraints of dynamo theory. Polar spots, inferred only from Doppler tomography, could plausibly be observational artifacts, casting some doubt on their very existence. Here we report imaging of the old, magnetically-active star $\zeta$ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two epochs, while 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 Sun is absent on $\zeta$ And, which hosts global spot patterns that cannot be produced by solar-type dynamos., Comment: Published in Nature, 33 pages, 8 figures, 5 tables

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., Accepted to ApJ, 11 pages, 6 figures, 8 tables