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

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

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

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.

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 have discovered a dwarf nova (DN) of type SU UMa in Kepler data which is 7.0 arcsec from the G-type exoplanet survey target KIC 4378554. The DN appears as a background source in the pixel aperture of the foreground G star. We extracted only the pixels where the DN is present and observed the source to undergo five outbursts -- one a superoutburst -- over a timespan of 22 months. The superoutburst was triggered by a normal outburst, a feature that has been seen in all DNe superoutburst observed by Kepler. Superhumps during the super outburst had a period of 1.842+/-0.004 h and we see a transition from disc-dominated superhump signal to a mix of disc and accretion stream impact. Predictions of the number of DNe present in Kepler data based on previously published space densities vary from 0.3 to 258. An investigation of the background pixels targets would lead to firmer constraints on the space density of DN.

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

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to explore the mechanism by which it was ejected from the open cluster of its birth. Here we present recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic binary and physical parameters of the O star primary. We also present the first XMM-Newton observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place limits on the emission from a proposed compact companion in the system, and we rule out a quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved secondary companion is also not detected in our optical spectra of the binary, and it is difficult to conclude that a gravitational interaction could have ejected this runaway binary with a low mass optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent black hole with conditions unfavorable for accretion at the time of our observations., Accepted to AJ

Non-radial pulsations (NRPs) are a proposed mechanism for the formation of decretion disks around Be stars and are important tools to study the internal structure of stars. NGC 3766 has an unusually large fraction of transient Be stars, so it is an excellent location to study the formation mechanism of Be star disks. High resolution spectroscopy can reveal line profile variations from NRPs, allowing measurements of both the degree, l, and azimuthal order, m. However, spectroscopic studies require large amounts of time with large telescopes to achieve the necessary high S/N and time domain coverage. On the other hand, multi-color photometry can be performed more easily with small telescopes to measure l only. Here, we present representative light curves of Be stars and non-emitting B stars in NGC 3766 from the CTIO 0.9m telescope in an effort to study NRPs in this cluster., 4 pages, to appear in the proceedings of IAU Symposium 266: Star Clusters

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.

A classical nova occurs when material accreting onto the surface of a white dwarf in a close binary system ignites in a thermonuclear runaway. Complex structures observed in the ejecta at late stages could result from interactions with the companion during the common envelope phase. Alternatively, the explosion could be intrinsically bipolar, resulting from a localized ignition on the surface of the white dwarf or as a consequence of rotational distortion. Studying the structure of novae during the earliest phases is challenging because of the high spatial resolution needed to measure their small sizes. Here we report near-infrared interferometric measurements of the angular size of Nova Delphini 2013, starting from one day after the explosion and continuing with extensive time coverage during the first 43 days. Changes in the apparent expansion rate can be explained by an explosion model consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of the ejected material changes as it expands. We detect an ellipticity in the light distribution, suggesting a prolate or bipolar structure that develops as early as the second day. Combining the angular expansion rate with radial velocity measurements, we derive a geometric distance to the nova of 4.54 +/- 0.59 kpc from the Sun., Comment: Published in Nature. 32 pages. Final version available at http://www.nature.com/nature/journal/v515/n7526/full/nature13834.html

We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: {\theta} Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753+/-0.009 mas for {\theta} Cyg, 0.539+/-0.007 mas for 16 Cyg A and 0.490+/-0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii ({\theta} Cyg: 1.48+/-0.02 Rsun, 16 Cyg A: 1.22+/-0.02 Rsun, 16 Cyg B: 1.12+/-0.02 Rsun), effective temperatures ({\theta} Cyg: 6749+/-44 K, 16 Cyg A: 5839+/-42 K, 16 Cyg B: 5809+/-39 K), and masses ({\theta} Cyg: 1.37+/-0.04 Msun, 16 Cyg A: 1.07+/-0.05 Msun, 16 Cyg B: 1.05+/-0.04 Msun) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts., Comment: 9 pages, 4 figures, and 5 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

The Kepler target KIC 5110407, a K-type star, shows strong quasi-periodic light curve fluctuations likely arising from the formation and decay of spots on the stellar surface rotating with a period of 3.4693 days. Using an established light-curve inversion algorithm, we study the evolution of the surface features based on Kepler space telescope light curves over a period of two years (with a gap of .25 years). At virtually all epochs, we detect at least one large spot group on the surface causing a 1-10% flux modulation in the Kepler passband. By identifying and tracking spot groups over a range of inferred latitudes, we measured the surface differential rotation to be much smaller than that found for the Sun. We also searched for a correlation between the seventeen stellar flares that occurred during our observations and the orientation of the dominant surface spot at the time of each flare. No statistically-significant correlation was found except perhaps for the very brightest flares, suggesting most flares are associated with regions devoid of spots or spots too small to be clearly discerned using our reconstruction technique. While we may see hints of long-term changes in the spot characteristics and flare statistics within our current dataset, a longer baseline of observation will be needed to detect the existence of a magnetic cycle in KIC 5110407., 32 pages, 15 figures, accepted to ApJ

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

We present the results of a study of differential rotation on the K2 IV primary of the RS CVn binary II Pegasi (HD 224085) performed by inverting light curves to produce images of the dark starspots on its surface. The data were obtained in the standard Johnson B and V filter passbands via the Tennessee State University T3 0.4-m Automated Photometric Telescope from JD 2447115.8086 to 2454136.6221 (1987 November 16 to 2007 February 5). The observations were subdivided into 68 data sets consisting of pairs of B and V light curves, which were then inverted using a constrained non-linear inversion algorithm that makes no a priori assumptions regarding the number of spots or their shapes. The resulting surface images were then assigned to 21 groups corresponding to time intervals over which we could observe the evolution of a given group of spots (except for three groups consisting of single data sets). Of these 21 groups, six showed convincing evidence of differential rotation over time intervals of several months. For the others, the spot configuration was such that differential rotation was neither exhibited nor contraindicated. The differential rotation we infer is in the same sense as that on the Sun: lower latitudes have shorter rotation periods. From plots of the range in longitude spanned by the spotted regions vs. time, we obtain estimates of the differential rotation coefficient k defined in earlier work by Henry et al., and show that our results for its value are consistent with the value obtained therein., Submitted to The Astronomical Journal

LS I +61 303 and LS 5039 are two of only a handful of known high mass X-ray binaries (HMXBs) that exhibit very high energy emission in the MeV-TeV range, and these "gamma-ray binaries" are of renewed interest due to the recent launch of the Fermi Gamma-ray Space Telescope. Here we present new radial velocities of both systems based on recent red and blue optical spectra. Both systems have somewhat discrepant orbital solutions available in the literature, and our new measurements result in improved orbital elements and resolve the disagreements. The improved geometry of each orbit will aid in studies of the high energy emission region near each source., Accepted to ApJ, 13 pages, preprint2 style

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

Massive O- and B-type stars are “cosmic engines” in the Universe and can be the dominant source of luminosity in a galaxy. The class of Be stars are rapidly rotating B-type stars that lose mass in an equatorial, circumstellar disk (Porter & Rivinius 2003) and cause Balmer and other line emission. Currently, we are unsure as to why these stars rotate so quickly but three scenarios are possible: they may have been born as rapid rotators, spun up by binary mass transfer, or spun up during the main-sequence evolution of B stars. In order to investigate these scenarios for this population of massive stars, we have been spectroscopically observing a set of 115 field Be stars with the Kitt Peak Coudè Feed telescope in both the Hα and Hγ wavelength regimes since 2004. This time baseline allows for examination of variability properties of the circumstellar disks as well as determine candidates for closer examination for binarity. We find that 90% of the observed stars show some variability with 8% showing significant variability over the 5-year baseline. Such values may be compared with the significant variability seen in some clusters such as NGC 3766 (McSwain 2008). Also, while ~20% of the sample consists of known binaries, we find that another 15–30% of the sample shows indications of binarity.

Nonradial pulsations (NRPs) are a possible formation mechanism for the equatorial disks surrounding Be stars. The open cluster NGC 3766 has a high fraction of transient Be stars, Be stars that have been observed with both emission due to a circumstellar disk and a non-emitting B-type spectrum. Because of the large fraction of transient Be stars, this cluster is a prime location for studying the formation mechanisms of Be star disks. We observed NGC 3766 for more than 25 nights over three years to generate Strömgren uvby light curves of the Be population. We present the results of a period search to investigate the presence of NRPs.