Search

Modeling based on differentiable programming holds great promise for astronomy, as it can employ techniques such as Hamiltonian Monte Carlo, gradient-based optimization, and other machine learning techniques. This new programming paradigm has motivated us to develop the first auto-differentiable spectrum model of exoplanets and brown dwarfs, ExoJAX (Kawahara et al. 2022). ExoJAX is designed to directly calculate cross-sections as functions of temperature and pressure, rather than interpolating tabulated data, to minimize errors in high-dispersion spectra modeling. However, its application was primarily proof-of-concept and limited to narrowband high-dispersion emission spectroscopy. In this paper, we have enhanced the differentiable opacity calculation using a new fast and memory-efficient algorithm, and have developed differentiable radiative transfer schemes, including emission, transmission, and reflection spectroscopy. These enhancements significantly expand the range of applications, as demonstrated through actual atmospheric retrievals: high-dispersion emission spectra of the brown dwarf GL229 B, medium-dispersion transmission spectra of the hot Saturn WASP-39 b from JWST, and high-dispersion reflection spectra of Jupiter. We obtained a C/O ratio for GL229 B consistent with its host star, constrained WASP-39 b's radial velocity from molecular fine structures at original resolution ($R \sim 3,000$), and estimated Jupiter's metallicity consistent with previous studies., Comment: 32 pages, 12 figures, submitted. ExoJAX is available at https://github.com/HajimeKawahara/exojax. ExoJAX2 is scheduled to be released after the review process

Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.

We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $\rho=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$., Comment: Published in ApJL

The distribution of small-scale magnetic fields in stellar photospheres is an important ingredient in our understanding of the magnetism of low mass stars. Their spatial distribution connects the field generated in the stellar interior with the outer corona and the large scale field, and thereby affects the space weather of planets. Unfortunately, we lack techniques that can locate them on most low-mass stars. One strategy is to localize field concentrations using the flares that occur in their vicinity. We explore a new method that adapts the spot simulation software fleck to study the modulation of flaring times as a function of active latitude. We use empirical relations to construct flare light curves similar to those available from Kepler and the Transiting Exoplanet Survey Satellite (TESS), search them for flares, and use the waiting times between flares to determine the location of active latitudes. We find that the mean and standard deviation of the waiting time distribution provide a unique diagnostic of flaring latitudes as a function of the number of active regions. Latitudes are best recovered when stars have three or less active regions that flare repeatedly, and active latitude widths below 20 deg; when either increases, the information about the active latitude location is gradually lost. We demonstrate our technique on a sample of flaring G dwarfs observed with the Kepler satellite, and furthermore suggest that combining ensemble methods for spots and flares could overcome the limitations of each individual technique for the localization of surface magnetic fields., Comment: 15 pages, 10 figures. Accepted to MNRAS. Repository with source code: https://github.com/ekaterinailin/flare-locations-ensembles-science Repository with data: https://zenodo.org/record/7996929

Fast yellow pulsating supergiants (FYPS) are a recently-discovered class of evolved massive pulsator. As candidate post-red supergiant objects, and one of the few classes of pulsating evolved massive stars, these objects have incredible potential to change our understanding of the structure and evolution of massive stars. Here we examine the lightcurves of a sample of 126 cool supergiants in the Magellanic Clouds observed by the Transiting Exoplanet Survey Satellite (\tess~) in order to identify pulsating stars. After making quality cuts and filtering out contaminant objects, we examine the distribution of pulsating stars in the Hertzprung-Russel (HR) diagram, and find that FYPS occupy a region above $\log L/L_\odot \gtrsim 5.0$. This luminosity boundary corresponds to stars with initial masses of $\sim$18-20 $M_\odot$, consistent with the most massive red supergiant progenitors of supernovae (SNe) II-P, as well as the observed properties of SNe IIb progenitors. This threshold is in agreement with the picture that FYPS are post-RSG stars. Finally, we characterize the behavior of FYPS pulsations as a function of their location in the HR diagram. We find low frequency pulsations at higher effective temperatures, higher frequency pulsations at lower temperatures, with a transition between the two behaviors at intermediate temperatures. The observed properties of FYPS make them fascinating objects for future theoretical study., Comment: Consistent with published version which contains significantly improved detection and rejection of contaminant objects. Comments welcome

Understanding magnetic activity on the surface of stars other than the Sun is important for exoplanet analyses to properly characterize an exoplanet's atmosphere and to further characterize stellar activity on a wide range of stars. Modeling stellar surface features of a variety of spectral types and rotation rates are key to understanding of the magnetic activity of these stars. Using data from Kepler, we use the starspot modeling program STarSPot (STSP) to measure the position and size of spots for KOI-340 which is an eclipsing binary consisting of a subgiant star ($\rm T_{\rm eff} = 5593 \pm 27 K; R_{*} = 1.98 \pm 0.05 R_{\odot}$) with an M-dwarf companion ($M_{*} = 0.214 \pm 0.006 M_{\odot}$). STSP uses a novel technique to measure the spot positions and radii by using the transiting secondary to study and model individual active regions on the stellar surface using high-precision photometry. We find the average size of spot features on KOI-340's primary is $\sim$10% the radius of the star, i.e. two times larger than the mean size of Solar-maximum sunspots. The spots on KOI-340 are present at every longitude and show possible signs of differential rotation. The minimum fractional spotted area of KOI-340's primary is $2\substack{+12\\ -2} \%$ while the spotted area of the Sun is at most 0.2%. One transit of KOI-340 shows a signal in the transit consistent with a plage; this plage occurs right before a dark spot indicating the plage and spot might be co-located on the surface of the star., Comment: 16 pages, 17 figures; Accepted to AJ

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., Comment: 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

White dwarf spectroscopy shows that nearly half of white dwarf atmospheres contain metals that must have been accreted from planetary material that survived the red giant phases of stellar evolution. We can use metal pollution in white dwarf atmospheres as flags, signalling recent accretion, in order to prioritize an efficient sample of white dwarfs to search for transiting material. We present a search for planetesimals orbiting six nearby white dwarfs with the CHEOPS spacecraft. The targets are relatively faint for CHEOPS, $11$ mag $< G < 12.8$ mag. We use aperture photometry data products from the CHEOPS mission as well as custom PSF photometry to search for periodic variations in flux due to transiting planetesimals. We detect no significant variations in flux that cannot be attributed to spacecraft systematics, despite reaching a photometric precision of $<2$ ppt in 60 s exposures on each target. We simulate observations to show that the small survey is sensitive primarily to Moon-sized transiting objects with periods $3$ hr $< P < 10$ hr, with radii $R \gtrsim 1000$ km., Comment: 10 pages, accepted as a Letter to the Editor in A&A

Magnetic activity in stars manifests as dark spots on their surfaces that modulate the brightness observed by telescopes. These light curves contain important information on stellar rotation. However, the accurate estimation of rotation periods is computationally expensive due to scarce ground truth information, noisy data, and large parameter spaces that lead to degenerate solutions. We harness the power of deep learning and successfully apply Convolutional Neural Networks to regress stellar rotation periods from Kepler light curves. Geometry-preserving time-series to image transformations of the light curves serve as inputs to a ResNet-18 based architecture which is trained through transfer learning. The McQuillan catalog of published rotation periods is used as ansatz to groundtruth. We benchmark the performance of our method against a random forest regressor, a 1D CNN, and the Auto-Correlation Function (ACF) - the current standard to estimate rotation periods. Despite limiting our input to fewer data points (1k), our model yields more accurate results and runs 350 times faster than ACF runs on the same number of data points and 10,000 times faster than ACF runs on 65k data points. With only minimal feature engineering our approach has impressive accuracy, motivating the application of deep learning to regress stellar parameters on an even larger scale, Comment: Accepted for the Machine Learning and the Physical Sciences Workshop, NeurIPS 2020

Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material floating in their atmospheres, in warm debris disks, or orbiting very closely, which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted. Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs., Comment: 50 pages, 12 figures, 2 tables. Published in Nature on Sept. 17, 2020. The final authenticated version is available online at: https://www.nature.com/articles/s41586-020-2713-y

Massive stars briefly pass through the yellow supergiant (YSG) phase as they evolve redward across the HR diagram and expand into red supergiants (RSGs). Higher-mass stars pass through the YSG phase again as they evolve blueward after experiencing significant RSG mass loss. These post-RSG objects offer us a tantalizing glimpse into which stars end their lives as RSGs, and why. One telltale sign of a post-RSG object may be an instability to pulsations, depending on the star's interior structure. Here we report the discovery of five YSGs with pulsation periods faster than 1 day, found in a sample of 76 cool supergiants observed by \tess at two-minute cadence. These pulsating YSGs are concentrated in a HR diagram region not previously associated with pulsations; we conclude that this is a genuine new class of pulsating star, Fast Yellow Pulsating Supergiants (FYPS). For each FYPS, we extract frequencies via iterative prewhitening and conduct a time-frequency analysis. One FYPS has an extracted frequency that is split into a triplet, and the amplitude of that peak is modulated on the same timescale as the frequency spacing of the triplet; neither rotation nor binary effects are likely culprits. We discuss the evolutionary status of FYPS and conclude that they are candidate post-RSGs. All stars in our sample also show the same stochastic low-frequency variability (SLFV) found in hot OB stars and attributed to internal gravity waves. Finally, we find four $\alpha$ Cygni variables in our sample, of which three are newly discovered., Comment: 29 pages, 13 figures, 8 tables. Accepted for publication in ApJ. Comments welcome

We present a method for detecting starspots on cool stars using the cross-correlation function (CCF) of high resolution molecular spectral templates applied to archival high-resolution spectra of G and K stars observed with HARPS/HARPS-N. We report non-detections of starspots on the Sun even when the Sun was spotted, the solar twin 18 Scorpii, and the very spotted Sun-like star HAT-P-11, suggesting that Sun-like starspot distributions will be invisible to the CCF technique, and should not produce molecular absorption signals which might be confused for signatures of exoplanet atmospheres. We detect strong TiO absorption in the T Tauri K-dwarfs LkCa 4 and AA Tau, consistent with significant coverage by cool regions. We show that despite the non-detections, the technique is sensitive to relatively small spot coverages on M dwarfs and large starspot areas on Sun-like stars., Comment: 12 pages, 8 figures, accepted to AJ

We investigate starspot distributions consistent with space-based photometry of F, G, and K stars in six stellar associations ranging in age from 10 Myr to 4 Gyr. We show that a simple light curve statistic called the "smoothed amplitude" is proportional to stellar age as $t^{-1/2}$, following a Skumanich-like spin-down relation. We marginalize over the unknown stellar inclinations by forward modeling the ensemble of light curves for direct comparison with the Kepler, K2 and TESS photometry. We sample the posterior distributions for spot coverage with Approximate Bayesian Computation. We find typical spot coverages in the range 1-10% which decrease with increasing stellar age. The spot coverage is proportional to $t^n$ where $n =-0.37 \pm 0.16$, also statistically consistent with a Skumanich-like $t^{-1/2}$ decay of starspot coverage with age. We apply two techniques to estimate the spot coverage of young exoplanet-hosting stars likely to be targeted for transmission spectroscopy with the James Webb Space Telescope, and estimate the bias in exoplanet radius measurements due to varying starspot coverage., Comment: 21 pages, 11 figures, accepted in ApJ, software in review at JOSS https://github.com/openjournals/joss-reviews/issues/2103

One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity images of the Sun taken by the HMI instrument aboard SDO to investigate the uncertainties introduced on the exoplanet radius measurements by stellar granulation and oscillations. After modeling the solar variability with a Gaussian process, we find that the amplitude of solar oscillations and granulation is of order 100 ppm -- similar to the depth of an Earth transit -- and introduces a fractional uncertainty on the depth of transit of 0.73% assuming four transits are observed over the mission duration. However, when we translate the depth measurement into a radius measurement of the planet, we find a much larger radius uncertainty of 3.6%. This is due to a degeneracy between the transit radius ratio, the limb-darkening, and the impact parameter caused by the inability to constrain the transit impact parameter in the presence of stellar variability. We find that surface brightness inhomogeneity due to photospheric granulation contributes a lower limit of only 2 ppm to the photometry in-transit. The radius uncertainty due to granulation and oscillations, combined with the degeneracy with the transit impact parameter, accounts for a significant fraction of the error budget of the PLATO mission, before detector or observational noise is introduced to the light curve. If it is possible to constrain the impact parameter or to obtain follow-up observations at longer wavelengths where limb-darkening is less significant, this may enable higher precision radius measurements., Comment: 11 pages, 8 figures, accepted in MNRAS

Star spot evolution is visible evidence of the emergence/decay of the magnetic field on stellar surface, and it is therefore important for the understanding of the underlying stellar dynamo and consequential stellar flares. In this paper, we report the temporal evolution of individual star spot area on the hot-Jupiter-hosting active solar-type star Kepler 17 whose transits occur every 1.5 days. The spot longitude and area evolution are estimated (1) from the stellar rotational modulations of Kepler data and (2) from the brightness enhancements during the exoplanet transits caused by existence of large star spots. As a result of the comparison, number of spots, spot locations, and the temporal evolution derived from the rotational modulations is largely different from those of in-transit spots. We confirm that although only two light curve minima appear per rotation, there are clearly many spots present on the star. We find that the observed differential intensity changes are sometimes consistent with the spot pattern detected by transits, but they sometimes do not match with each other. Although the temporal evolution derived from the rotational modulation differs from those of in-transit spots to a certain degree, the emergence/decay rates of in-transit spots are within an order of magnitude of those derived for sunspots as well as our previous research based only on rotational modulations. This supports a hypothesis that the emergence/decay of sunspots and extremely-large star spots on solar-type stars occur through the same underlying processes., Comment: 37 pages, 12 figures, 1 table. Accepted for publication in The Astrophysical Journal

We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future., Comment: Accepted for publication in AJ, 22 pages, 15 figures

We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective temperatures from the Gaia DR2 pipeline. We also measured their chromospheric $S$ and $\log R^\prime_\mathrm{HK}$ indices to classify the stars as active or inactive and study the connection between chromospheric activity and starspots. We then attempted to infer the starspot covering fractions on the active stars by modeling their spectra as a linear combination of hot and cool inactive stellar spectral templates. We find that it is critical to use precise colors of the stars to place stringent priors on the plausible spot covering fractions. The inferred spot covering fractions generally increase with the chromospheric activity indicator $\log R^\prime_\mathrm{HK}$, though we are largely insensitive to spot coverages smaller than $f_S \lesssim 20$%. We find a dearth of stars with small $\log R^\prime_\mathrm{HK}$ and significant spot coverages., Comment: Accepted for publication in AJ

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$., Comment: Science white paper submitted in response to the the U.S. National Academies of Science, Engineering, and Medicine's call for community input to the Astro2020 Decadal Survey; 9 pages, 3 figures

We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online., Comment: Accepted for publication in MNRAS

We present the results of an abundance analysis of three G and two K dwarfs in the Praesepe open cluster based on high-resolution, moderate signal-to-noise-ratio spectra obtained with the ARC 3.5-m Telescope at Apache Point Observatory. Using a Principle Component Analysis and the BACCHUS automated spectral analysis code, we determined stellar parameters and abundances of up to 24 elements for each of our targets, which range in temperature from 6000 to 4600 K. The average derived iron abundance for the three G stars is 0.17+/-0.07 dex, consistent with the 0.12+/-0.04 dex derived by Boesgaard et al. (2013) for their sample of 11 solar-type Praesepe members, which included these G stars. To investigate the efficacy of using automated routines to derive the abundances of cooler main-sequence stars, we compared the abundances of the K dwarfs to those of the G dwarfs. Our abundances are consistent to <=0.1~dex for 13 of the 18 elements we report for all five of the stars, providing more evidence that G and K stars in a given open cluster are chemically homogeneous. The median difference between the mean G and K stars abundances is 0.08+/-0.05 dex, despite serious challenges with the noisier data for the fainter K dwarfs. Our results are encouraging for chemical tagging, as they indicate that it may be possible to use automated abundance determination techniques to identify chemically related main-sequence stars across larger temperature ranges than are usually considered in these experiments., Comment: 11 pages (plus 29 page long table), 3 figures; accepted by ApJ

We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation consistent with a bright hemisphere followed by a relatively dark hemisphere, generating photometric variability with a semi-amplitude of 4%. We find that KIC 9652680 is darkest when its $S$-index of Ca II H & K emission is at its maximum. We interpret this anti-correlation between flux and $S$ to indicate that dark starspots in the photosphere are co-located with the bright plages in the chromosphere, as they are on the Sun. Moving to lower masses and slower rotators, we present K2 observations with simultaneous spectroscopy of EPIC 211928486 (K5V) and EPIC 211966629 (K4V), two active stars in the 650 Myr-old open cluster Praesepe. The K2 photometry reveals that both stars have rotation periods of 11.7 days; while their flux varies by 1 and 2% respectively, their Ca II H & K $S$-indices seem to hold relatively constant as a function of rotational phase. This suggests that extended chromospheric networks of plages are not concentrated into regions of emission centered on the starspots that drive rotational modulation, unlike KIC 9652680. We also note that the Ca II emission of EPIC 211928486 dipped and recovered suddenly over the duration of one rotation, suggesting that the evolution timescale of plages may be of order the rotation period., Comment: Accepted for publication in AJ

We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influences the transit depths but not the ingress/egress durations. We find that none of the planets show statistically significant evidence for self-contamination by bright or dark regions of the stellar photosphere. However, we show that small-scale magnetic activity, analogous in size to the smallest sunspots, could still be lurking in the transit photometry undetected., Comment: Accepted for publication in ApJL

We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted for. We reparameterize the transit light curve to account for "self-contamination" by photospheric inhomogeneities by splitting the parameter $R_p/R_\star$ into two parameters: one for the radius ratio -- which controls the duration of ingress and egress -- and another which measures the possibly contaminated transit depth. We show that this is equivalent to the formulation for contamination by a second star (with positive or negative flux), and that it is sensitive to time-steady inhomogeneity of the stellar photosphere. We use synthetic light curves of spotted stars at high signal-to-noise to show that the radius recovered from measurement of the ingress/egress duration can recover the true radii of planets transiting spotted stars with axisymmetric spot distributions if the limb-darkening parameters are precisely known. We fit time-averaged high signal-to-noise transit light curves from Kepler and Spitzer of ten planets to measure the planet radii and search for evidence of spot distributions. We find that this sample has a range of measured depths and ingress durations which are self-consistent, providing no strong evidence for contamination by spots. However, there is suggestive evidence for occultation of starspots on Kepler-17, and that relatively bright regions are occulted by the planets of Kepler-412 and HD 80606. Future observations with the James Webb Space Telescope may enable this technique to yield accurate planetary radii in the presence of stellar inhomogeneities., Comment: Accepted for publication in AJ

The M8V star TRAPPIST-1 hosts seven roughly Earth-sized planets and is a promising target for exoplanet characterization. Kepler/K2 Campaign 12 observations of TRAPPIST-1 in the optical show an apparent rotational modulation with a 3.3 day period, though that rotational signal is not readily detected in the Spitzer light curve at 4.5 $\mu$m. If the rotational modulation is due to starspots, persistent dark spots can be excluded from the lack of photometric variability in the Spitzer light curve. We construct a photometric model for rotational modulation due to photospheric bright spots on TRAPPIST-1 which is consistent with both the Kepler and Spitzer light curves. The maximum-likelihood model with three spots has typical spot sizes of $R_\mathrm{spot}/R_\star \approx 0.004$ at temperature $T_\mathrm{spot} \gtrsim 5300 \pm 200$ K. We also find that large flares are observed more often when the brightest spot is facing the observer, suggesting a correlation between the position of the bright spots and flare events. In addition, these flares may occur preferentially when the spots are increasing in brightness, which suggests that the 3.3 d periodicity may not be a rotational signal, but rather a characteristic timescale of active regions., Comment: Accepted by ApJ March 12, 2018

Metal pollution in white dwarf atmospheres may be the accreted remnants of planetary objects. After the discovery of disintegrating planetary objects transiting WD 1145+017 (Vanderburg et al. 2015), undergraduates in the University of Washington's Pre- Major in Astronomy Program (Pre-MAP) were inspired to collect photometry of the brightest white dwarfs to hunt for similar transiting objects around other metal-polluted white dwarfs. Prior surveys have yet to make a detection of a transiting planet orbiting a white dwarf, yet white dwarfs are still an attractive target for searches of small, rocky planetary material. Since a typical white dwarf is Earth-sized, transits of Earth-sized planets could have depths >50%, so even low S/N photometry has a chance at discovering transiting material. We identified bright, northern, metal-polluted white dwarfs in the SDSS DR10 white dwarf catalog, and observed five targets with the Astrophysical Research Consortium Small Aperture Telescope (ARCSAT) 0.5-meter telescope at Apache Point Observatory. The ARCSAT photometry had sufficient precision to detect Moon-sized objects or larger at short orbital periods, though no such planets were detected for these targets. We look forward to surveys which may find planets orbiting white dwarfs, such as NASA's TESS, ESA's PLATO, and the Evryscope., Comment: Accepted to RNAAS

Astrometry from Gaia will measure the positions of stellar photometric centroids to unprecedented precision. We show that the precision of Gaia astrometry is sufficient to detect starspot-induced centroid jitter for nearby stars in the Tycho-Gaia Astrometric Solution (TGAS) sample with magnetic activity similar to the young G-star KIC 7174505 or the active M4 dwarf GJ 1243, but is insufficient to measure centroid jitter for stars with Sun-like spot distributions. We simulate Gaia observations of stars with 10 year activity cycles to search for evidence of activity cycles, and find that Gaia astrometry alone likely can not detect activity cycles for stars in the TGAS sample, even if they have spot distributions like KIC 7174505. We review the activity of the nearby low-mass stars in the TGAS sample for which we anticipate significant detections of spot-induced jitter., Comment: Accepted for publication in MNRAS, Feb 27, 2018

HAT-P-11 is a planet-hosting K4V star in the Kepler field, with an activity cycle that bear similarities to the Sun's. The chromospheric activity of HAT-P-11 indicates that a new activity cycle is beginning. We report ground-based observations with holographic diffuser photometry to measure the starspots of HAT-P-11 in its second observed magnetic activity cycle (Cycle 1). We find the area coverage of starspots within the transit chord for UTC 2017-10-30 is 14% --- which makes this transit the most spotted HAT-P-11 transit observed to date. We suggest that we are likely observing occultations of large spot groups appearing at the beginning of Cycle 1., Comment: Submitted to RNAAS

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process., Comment: 19 pages, 8 figures, accepted for publication in ApJL

Kepler photometry of the hot Neptune host star HAT-P-11 suggests that its spot latitude distribution is comparable to the Sun's near solar maximum. We search for evidence of an activity cycle in the CaII H & K chromospheric emission $S$-index with archival Keck/HIRES spectra and observations from the echelle spectrograph on the ARC 3.5 m Telescope at APO. The chromospheric emission of HAT-P-11 is consistent with a $\gtrsim 10$ year activity cycle, which plateaued near maximum during the Kepler mission. In the cycle that we observed, the star seemed to spend more time near active maximum than minimum. We compare the $\log R^\prime_{HK}$ normalized chromospheric emission index of HAT-P-11 with other stars. HAT-P-11 has unusually strong chromospheric emission compared to planet-hosting stars of similar effective temperature and rotation period, perhaps due to tides raised by its planet., Comment: 16 pages, 8 figures; accepted to the Astrophysical Journal

We measure the starspot radii and latitude distribution on the K4 dwarf HAT-P-11 from Kepler short-cadence photometry. We take advantage of starspot occultations by its highly-misaligned planet to compare the spot size and latitude distributions to those of sunspots. We find that the spots of HAT-P-11 are distributed in latitude much like sunspots near solar activity maximum, with mean spot latitude of $16 \pm 1$ degrees. The majority of starspots of HAT-P-11 have physical sizes that closely resemble the sizes of sunspots at solar maximum. We estimate the mean spotted area coverage on HAT-P-11 is $3^{+6}_{-1}\%$, roughly two orders of magnitude greater than the typical solar spotted area., Comment: 28 pages, 16 figures, accepted for publication in the Astrophysical Journal (August 7, 2017)

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf., Comment: 42 pages, 8 figures, 2 tables

We present observations spanning 355 orbital phases of HAT-P-7 observed by Kepler from May 2009 to March 2011 (Q1-9). We find a shallower secondary eclipse depth than initially announced, consistent with a low optical albedo and detection of nearly exclusively thermal emission, without a reflected light component. We find an approximately 10 ppm perturbation to the average transit light curve near phase -0.02 that we attribute to a temperature decrease on the surface of the star, phased to the orbit of the planet. This cooler spot is consistent with planet-induced gravity darkening, slightly lagging the sub-planet position due to the finite response time of the stellar atmosphere. The brightness temperature of HAT-P-7b in the Kepler bandpass is T_B = 2733 +/- 21 K and the amplitude of the deviation in stellar surface temperature due to gravity darkening is approximately -0.18 K. The detection of the spot is not statistically unequivocal due its small amplitude, though additional Kepler observations should be able to verify the astrophysical nature of the anomaly., Comment: 5 pages, 3 figures; accepted by ApJL on Jan 18, 2013