Showing 1–50 of 100 results for author: Davenport, J R A

A Multiwavelength Survey of Nearby M dwarfs: Optical and Near-Ultraviolet Flares and Activity with Contemporaneous TESS, Kepler/K2, \textit{Swift}, and HST Observations

Authors: Rishi R. Paudel, Thomas Barclay, Allison Youngblood, Elisa V. Quintana, Joshua E. Schlieder, Laura D. Vega, Emily A. Gilbert, Rachel A. Osten, Sarah Peacock, Isaiah I. Tristan, Dax L. Feliz, Patricia T. Boyd, James R. A. Davenport, Daniel Huber, Adam F. Kowalski, Teresa A. Monsue, Michele L. Silverstein

Abstract: We present a comprehensive multiwavelength investigation into flares and activity in nearby M~dwarf stars. We leverage the most extensive contemporaneous dataset obtained through the Transiting Exoplanet Sky Survey (TESS), Kepler/K2, the Neil Gehrels Swift Observatory (\textit{Swift}), and the Hubble Space Telescope (HST), spanning the optical and near-ultraviolet (NUV) regimes. In total, we obser… ▽ More We present a comprehensive multiwavelength investigation into flares and activity in nearby M~dwarf stars. We leverage the most extensive contemporaneous dataset obtained through the Transiting Exoplanet Sky Survey (TESS), Kepler/K2, the Neil Gehrels Swift Observatory (\textit{Swift}), and the Hubble Space Telescope (HST), spanning the optical and near-ultraviolet (NUV) regimes. In total, we observed 213 NUV flares on 24 nearby M dwarfs, with $\sim$27\% of them having detected optical counterparts, and found that all optical flares had NUV counterparts. We explore NUV/optical energy fractionation in M dwarf flares. Our findings reveal a slight decrease in the ratio of optical to NUV energies with increasing NUV energies, a trend in agreement with prior investigations on G-K stars' flares at higher energies. Our analysis yields an average NUV fraction of flaring time for M0-M3 dwarfs of 2.1\%, while for M4-M6 dwarfs, it is 5\%. We present an empirical relationship between NUV and optical flare energies and compare to predictions from radiative-hydrodynamic and blackbody models. We conducted a comparison of the flare frequency distribution (FFDs) of NUV and optical flares, revealing the FFDs of both NUV and optical flares exhibit comparable slopes across all spectral subtypes. NUV flares on stars affect the atmospheric chemistry, the radiation environment, and the overall potential to sustain life on any exoplanets they host. We find that early and mid-M dwarfs (M0-M5) have the potential to generate NUV flares capable of initiating abiogenesis. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 41 pages, 22 figures, Accepted for publication in the Astrophysical Journal

Searching the SN 1987A SETI Ellipsoid with TESS

Authors: Bárbara Cabrales, James R. A. Davenport, Sofia Z. Sheikh, Steve Croft, Andrew P. V. Siemion, Daniel Giles, Ann Marie Cody

Abstract: The SETI Ellipsoid is a strategy for technosignature candidate selection which assumes that extraterrestrial civilizations who have observed a galactic-scale event -- such as supernova 1987A -- may use it as a Schelling point to broadcast synchronized signals indicating their presence. Continuous wide-field surveys of the sky offer a powerful new opportunity to look for these signals, compensating… ▽ More The SETI Ellipsoid is a strategy for technosignature candidate selection which assumes that extraterrestrial civilizations who have observed a galactic-scale event -- such as supernova 1987A -- may use it as a Schelling point to broadcast synchronized signals indicating their presence. Continuous wide-field surveys of the sky offer a powerful new opportunity to look for these signals, compensating for the uncertainty in their estimated time of arrival. We explore sources in the TESS continuous viewing zone, which corresponds to 5% of all TESS data, observed during the first three years of the mission. Using improved 3D locations for stars from Gaia Early Data Release 3, we identified 32 SN 1987A SETI Ellipsoid targets in the TESS continuous viewing zone with uncertainties better than 0.5 ly. We examined the TESS light curves of these stars during the Ellipsoid crossing event and found no anomalous signatures. We discuss ways to expand this methodology to other surveys, more targets, and different potential signal types. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 17 pages, 6 figures, AJ published

Every Datapoint Counts: Stellar Flares as a Case Study of Atmosphere Aided Studies of Transients in the LSST Era

Authors: Riley W. Clarke, James R. A. Davenport, John Gizis, Melissa L. Graham, Xiaolong Li, Willow Fortino, Ian Sullivan, Yusra Alsayyad, James Bosch, Robert A. Knop, Federica Bianco

Abstract: Due to their short timescale, stellar flares are a challenging target for the most modern synoptic sky surveys. The upcoming Vera C. Rubin Legacy Survey of Space and Time (LSST), a project designed to collect more data than any precursor survey, is unlikely to detect flares with more than one data point in its main survey. We developed a methodology to enable LSST studies of stellar flares, with a… ▽ More Due to their short timescale, stellar flares are a challenging target for the most modern synoptic sky surveys. The upcoming Vera C. Rubin Legacy Survey of Space and Time (LSST), a project designed to collect more data than any precursor survey, is unlikely to detect flares with more than one data point in its main survey. We developed a methodology to enable LSST studies of stellar flares, with a focus on flare temperature and temperature evolution, which remain poorly constrained compared to flare morphology. By leveraging the sensitivity expected from the Rubin system, Differential Chromatic Refraction can be used to constrain flare temperature from a single-epoch detection, which will enable statistical studies of flare temperatures and constrain models of the physical processes behind flare emission using the unprecedentedly high volume of data produced by Rubin over the 10-year LSST. We model the refraction effect as a function of the atmospheric column density, photometric filter, and temperature of the flare, and show that flare temperatures at or above ~4,000K can be constrained by a single g-band observation at airmass X > 1.2, given the minimum specified requirement on single-visit relative astrometric accuracy of LSST, and that a surprisingly large number of LSST observations is in fact likely be conducted at X > 1.2, in spite of image quality requirements pushing the survey to preferentially low X. Having failed to measure flare DCR in LSST precursor surveys, we make recommendations on survey design and data products that enable these studies in LSST and other future surveys. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 18 pages, 16 figures, 1 table. Submitted to the Astrophysical Journal Supplement Series

The weird and the wonderful in our Solar System: Searching for serendipity in the Legacy Survey of Space and Time

Authors: Brian Rogers, Chris J. Lintott, Steve Croft, Megan E. Schwamb, James R. A. Davenport

Abstract: We present a novel method for anomaly detection in Solar System object data, in preparation for the Legacy Survey of Space and Time. We train a deep autoencoder for anomaly detection and use the learned latent space to search for other interesting objects. We demonstrate the efficacy of the autoencoder approach by finding interesting examples, such as interstellar objects, and show that using the… ▽ More We present a novel method for anomaly detection in Solar System object data, in preparation for the Legacy Survey of Space and Time. We train a deep autoencoder for anomaly detection and use the learned latent space to search for other interesting objects. We demonstrate the efficacy of the autoencoder approach by finding interesting examples, such as interstellar objects, and show that using the autoencoder, further examples of interesting classes can be found. We also investigate the limits of classic unsupervised approaches to anomaly detection through the generation of synthetic anomalies and evaluate the feasibility of using a supervised learning approach. Future work should consider expanding the feature space to increase the variety of anomalies that can be uncovered during the survey using an autoencoder. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted by AJ

The Sonified Hertzsprung-Russell Diagram

Authors: Daniela Huppenkothen, Juan Pampin, James R. A. Davenport, James Wenlock

Abstract: Understanding the physical properties of stars, and putting these properties into the context of stellar evolution, is a core challenge in astronomical research. A key visualization in studying stellar evolution is the Hertzsprung-Russell diagram (HRD), organizing data about stellar luminosity and colour into a form that is informative about stellar structure and evolution. However, connecting the… ▽ More Understanding the physical properties of stars, and putting these properties into the context of stellar evolution, is a core challenge in astronomical research. A key visualization in studying stellar evolution is the Hertzsprung-Russell diagram (HRD), organizing data about stellar luminosity and colour into a form that is informative about stellar structure and evolution. However, connecting the HRD with other sources of information, including stellar time series, is an outstanding challenge. Here we present a new method to turn stellar time series into sound. This method encodes physically meaningful features such that auditory comparisons between sonifications of different stars preserve astrophysical differences between them. We present an interactive multimedia version of the HRD that combines both visual and auditory components and that allows exploration of different types of stars both on and off the main sequence through both visual and auditory media. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Comments: 8 pages, 5 figures; accepted for publication in the proceedings of "The 28th International Conference on Auditory Display (ICAD 2023) - Special Session on Astronomical Data Sonification"

A Pilot Search for Gravitational Self-Lensing Binaries with the Zwicky Transient Facility

Authors: Allison Crossland, Eric C. Bellm, Courtney Klein, James R. A. Davenport, Thomas Kupfer, Steven L. Groom, Russ R. Laher, Reed Riddle

Abstract: Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such "self-lensing" signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility… ▽ More Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such "self-lensing" signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility (ZTF). We present a systematic search for self-lensing signals in Galactic binaries from a subset of high-cadence ZTF data taken in 2018. We identify 19 plausible candidates from the search, although because each candidate is observed to only brighten once, other origins such as stellar flares are more likely. We discuss prospects for more comprehensive future searches of the ZTF data. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 12 pages. Submitted to the Open Journal of Astrophysics

APO & SMARTS flare star campaign observations I. Blue wing asymmetries in chromospheric lines during mid M dwarf flares from simultaneous spectroscopic and photometric observation data

Authors: Yuta Notsu, Adam F. Kowalski, Hiroyuki Maehara, Kosuke Namekata, Kenji Hamaguchi, Teruaki Enoto, Isaiah I. Tristan, Suzanne L. Hawley, James R. A. Davenport, Satoshi Honda, Kai Ikuta, Shun Inoue, Keiichi Namizaki, Daisaku Nogami, Kazunari Shibata

Abstract: We conducted the time-resolved simultaneous optical spectroscopic and photometric observations of mid M dwarf flare stars YZ CMi, EV Lac, and AD Leo. Spectroscopic observations were obtained using Apache Point Observatory 3.5m and Small \& Moderate Aperture Research Telescope System 1.5m telescopes during 31 nights. Among the 41 detected flares, seven flares showed clear blue wing asymmetries in t… ▽ More We conducted the time-resolved simultaneous optical spectroscopic and photometric observations of mid M dwarf flare stars YZ CMi, EV Lac, and AD Leo. Spectroscopic observations were obtained using Apache Point Observatory 3.5m and Small \& Moderate Aperture Research Telescope System 1.5m telescopes during 31 nights. Among the 41 detected flares, seven flares showed clear blue wing asymmetries in the H$αあるふぁ$ line, with various correspondences in flare properties. The duration of the blue wing asymmetries range from 20 min to 2.5 hours, including a flare showing the shift from blue to red wing asymmetry. Blue wing asymmetries can be observed during both white-light and candidate non white-light flares. All of the seven flares showed blue wing asymmetries also in the H$βべーた$ line, but there are large varieties on which other chromospheric lines showed blue wing asymmetries. One among the 7 flares was also observed with soft X-ray spectroscopy, which enabled us to estimate the flare magnetic field and length of the flare loop. The line-of-sight velocities of the blue-shifted components range from -73 to -122 km s$^{-1}$. Assuming that the blue-shifts were caused by prominence eruptions, the mass of upward moving plasma was estimated to be 10$^{15}$ -- 10$^{19}$ g, which are roughly on the relation between flare energy and erupting mass expected from solar coronal mass ejections (CMEs). Although further investigations are necessary for understanding the observed various properties, these possible prominence eruptions on M-dwarfs could evolve into CMEs, assuming the similar acceleration mechanism with solar eruptions. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 123 pages, 102 figures, 7 tables. Accepted for publication in The Astrophysical Journal

Signal Synchronization Strategies and Time Domain SETI with Gaia DR3

Authors: Andy Nilipour, James R. A. Davenport, Steve Croft, Andrew P. V. Siemion

Abstract: Spatiotemporal techniques for signal coordination with actively transmitting extraterrestrial civilizations, without the need for prior communication, can constrain technosignature searches to a significantly smaller coordinate space. With the variable star catalog from Gaia Data Release 3, we explore two related signaling strategies: the SETI Ellipsoid, and that proposed by Seto, which are both b… ▽ More Spatiotemporal techniques for signal coordination with actively transmitting extraterrestrial civilizations, without the need for prior communication, can constrain technosignature searches to a significantly smaller coordinate space. With the variable star catalog from Gaia Data Release 3, we explore two related signaling strategies: the SETI Ellipsoid, and that proposed by Seto, which are both based on the synchronization of transmissions with a conspicuous astrophysical event. This dataset contains more than 10 million variable star candidates with light curves from the first three years of Gaia's operational phase, between 2014 and 2017. Using four different historical supernovae as source events, we find that less than 0.01% of stars in the sample have crossing times, the times at which we would expect to receive synchronized signals on Earth, within the date range of available Gaia observations. For these stars, we present a framework for technosignature analysis that searches for modulations in the variability parameters by splitting the stellar light curve at the crossing time. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: 22 pages, 13 figures, published in AJ

Journal ref: 2023 AJ 166 79

Catalog of Integrated-Light Star Cluster Light Curves in TESS

Authors: Tobin M. Wainer, Gail Zasowski, Joshua Pepper, Tom Wagg, Christina L. Hedges, Vijith Jacob Poovelil, Tara Fetherolf, James R. A. Davenport, P. Marios Christodoulou, Jack T. Dinsmore, Avi Patel, Kameron Goold, Benjamin J. Gibson

Abstract: We present the first integrated light, TESS-based light curves for star clusters in the Milky Way, Small Magellanic Cloud, and Large Magellanic Cloud. We explore the information encoded in these light curves, with particular emphasis on variability. We describe our publicly available package ELK, which is designed to extract the light curves by applying principal component analysis to perform back… ▽ More We present the first integrated light, TESS-based light curves for star clusters in the Milky Way, Small Magellanic Cloud, and Large Magellanic Cloud. We explore the information encoded in these light curves, with particular emphasis on variability. We describe our publicly available package ELK, which is designed to extract the light curves by applying principal component analysis to perform background light correction, and incorporating corrections for TESS systematics, allowing us to detect variability on time scales shorter than ~10 days. We perform a series of checks to ensure the quality of our light curves, removing observations where systematics are identified as dominant features, and deliver light curves for 348 previously-cataloged open and globular clusters. Where TESS has observed a cluster in more than one observing sectors, we provide separate light curves for each sector (for a total of 2204 light curves). We explore in detail the light curves of star clusters known to contain high-amplitude Cepheid and RR Lyrae variable stars, and confirm that the variability of these known variables is still detectable when summed together with the light from thousands of other stars. We also demonstrate that even some low-amplitude stellar variability is preserved when integrating over a stellar population. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 17 pages, 9 Figures, 3 Tables, Accepted to AAS Journals

Journal ref: AJ 166 106 (2023)

Enabling Stellar Flare Science in the Roman Galactic Bulge Survey: Cadence, Filters, and the Read-Out Strategy Matter

Authors: Guadalupe Tovar Mendoza, Robert F. Wilson, Allison Youngblood, Laura D. Vega, Thomas Barclay, James R. A. Davenport, Jordan Ealy

Abstract: As was discovered with other wide field, precise imagers, the stable photometry necessary for the microlensing surveys is well-suited to general stellar astrophysics, including stellar flares, which are important for understanding stellar magnetic activity and even the space weather environments of exoplanets. Large stellar flare surveys have never been performed before in the Roman spectral range… ▽ More As was discovered with other wide field, precise imagers, the stable photometry necessary for the microlensing surveys is well-suited to general stellar astrophysics, including stellar flares, which are important for understanding stellar magnetic activity and even the space weather environments of exoplanets. Large stellar flare surveys have never been performed before in the Roman spectral range, and Roman may reveal new information about flare emission mechanisms (blackbody, recombination continuum, chromospheric emission lines) and how flare rates change with stellar age and metallicity. For instance, the Galactic Bulge stars will be much older than the typical studied flare stars, and Roman's wide field and exquisite imaging may provide sufficient statistics to probe the flare behavior and properties of such an old stellar population. However, the information yield will likely depend on sky location, cadence, read-out strategy, and filter choices. Stellar flare timescales range from seconds to hours, so Roman may only be able to resolve the longest, most energetic events. However, because a single exposure in the Galactic Bulge Time Domain Survey will consist of several non-destructive reads, short duration events can be modeled from the flux variations within a single exposure. Here we provide a proof of concept, showing that flare morphologies can be significantly better constrained if sub-exposure data are analyzed. As a result, we advocate that such data be made publicly available for Roman flare studies, with minimal on-board processing. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: Roman Core Community Survey White Paper, 6 pages, 2 figures

Discovery of Gaia17bpp, a Giant Star with the Deepest and Longest Known Dimming Event

Authors: Anastasios Tzanidakis, James R. A. Davenport, Eric C. Bellm, Yuankun Wang

Abstract: We report the serendipitous discovery of Gaia17bpp/2MASS J19372316+1759029, a star with a deep single large-amplitude dimming event of $\sim$4.5 magnitudes that lasted over 6.5 years. Using the optical to IR spectral energy distribution (SED), we constrain the primary star to be a cool giant M0-III star with effective temperature $T_{\text{eff}}$=3,850 K and radius R=58 R$_{\odot}$. Based on the S… ▽ More We report the serendipitous discovery of Gaia17bpp/2MASS J19372316+1759029, a star with a deep single large-amplitude dimming event of $\sim$4.5 magnitudes that lasted over 6.5 years. Using the optical to IR spectral energy distribution (SED), we constrain the primary star to be a cool giant M0-III star with effective temperature $T_{\text{eff}}$=3,850 K and radius R=58 R$_{\odot}$. Based on the SED fitting, we obtained a bimodal posterior distribution of primary stellar masses at 1.5 M${\odot}$ and 3.7 M${\odot}$. Within the last 66 years of photometric coverage, no other significant dimming events of this depth and duration were identified in the optical light curves. Using a Gaussian Process, we fit a high-order Gaussian model to the optical and IR light curves and conclude the dimming event exhibits moderate asymmetries from optical to IR. At the minimum of the dimming event, the (W$_{1}$-W$_{2}$) color was bluer by $\sim$0.2 mag relative to the primary star outside the dimming event. The ingress and egress colors show a shallow reddening profile. We suggest that the main culprit of the dimming event is likely due to the presence of a large, optically thick disk transiting the primary giant star. By fitting a monochromatic transit model of an oblate disk transiting a star, we found good agreement with a slow-moving, 0.005 km sec$^{-1}$, disk with a $\sim$1.4 AUえーゆー radius. We propose that Gaia17bpp belongs to a rare binary star population similar to the Epsilon Aurigae system, which consists of a secondary star enshrouded by an optically thick debris disk. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 16 pages, 12 figures, submitted to ApJ, comments welcome!

Real-Time Technosignature Strategies with SN 2023ixf

Authors: James R. A. Davenport, Sofia Z. Sheikh, Wael Farah, Andy Nilipour, Bárbara Cabrales, Steve Croft, Alexander W. Pollak, Andrew P. V. Siemion

Abstract: Several technosignature techniques focus on historic events such as SN 1987A as the basis to search for coordinated signal broadcasts from extraterrestrial agents. The recently discovered SN 2023ixf in the spiral galaxy M101 is the nearest Type II supernova in over a decade, and will serve as an important benchmark event. Here we review the potential for SN 2023ixf to advance ongoing techonsignatu… ▽ More Several technosignature techniques focus on historic events such as SN 1987A as the basis to search for coordinated signal broadcasts from extraterrestrial agents. The recently discovered SN 2023ixf in the spiral galaxy M101 is the nearest Type II supernova in over a decade, and will serve as an important benchmark event. Here we review the potential for SN 2023ixf to advance ongoing techonsignature searches, particularly signal-synchronization techniques such as the "SETI Ellipsoid". We find that more than 100 stars within 100 pc are already close to intersecting this SETI Ellipsoid, providing numerous targets for real-time monitoring within ~3$^\circ$ of SN 2023ixf. We are commencing a radio technosignature monitoring campaign of these targets with the Allen Telescope Array and the Green Bank Telescope. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 3 Pages, 1 Figure, Submitted to RNAAS. Collaboration welcomed

Deep Drilling in the Time Domain with DECam: Survey Characterization

Authors: Melissa L. Graham, Robert A. Knop, Thomas Kennedy, Peter E. Nugent, Eric Bellm, Márcio Catelan, Avi Patel, Hayden Smotherman, Monika Soraisam, Steven Stetzler, Lauren N. Aldoroty, Autumn Awbrey, Karina Baeza-Villagra, Pedro H. Bernardinelli, Federica Bianco, Dillon Brout, Riley Clarke, William I. Clarkson, Thomas Collett, James R. A. Davenport, Shenming Fu, John E. Gizis, Ari Heinze, Lei Hu, Saurabh W. Jha , et al. (19 additional authors not shown)

Abstract: This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5… ▽ More This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 22 pages, 17 figures, 2 tables. Accepted to MNRAS

Rubin Observatory LSST Transients and Variable Stars Roadmap

Authors: Kelly M. Hambleton, Federica B. Bianco, Rachel Street, Keaton Bell, David Buckley, Melissa Graham, Nina Hernitschek, Michael B. Lund, Elena Mason, Joshua Pepper, Andrej Prsa, Markus Rabus, Claudia M. Raiteri, Robert Szabo, Paula Szkody, Igor Andreoni, Simone Antoniucci, Barbara Balmaverde, Eric Bellm, Rosaria Bonito, Giuseppe Bono, Maria Teresa Botticella, Enzo Brocato, Katja Bucar Bricman, Enrico Cappellaro , et al. (57 additional authors not shown)

Abstract: The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the T… ▽ More The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 202 pages (in book format) 34 figures plus chapter heading figures (13)

From Data to Software to Science with the Rubin Observatory LSST

Authors: Katelyn Breivik, Andrew J. Connolly, K. E. Saavik Ford, Mario Jurić, Rachel Mandelbaum, Adam A. Miller, Dara Norman, Knut Olsen, William O'Mullane, Adrian Price-Whelan, Timothy Sacco, J. L. Sokoloski, Ashley Villar, Viviana Acquaviva, Tomas Ahumada, Yusra AlSayyad, Catarina S. Alves, Igor Andreoni, Timo Anguita, Henry J. Best, Federica B. Bianco, Rosaria Bonito, Andrew Bradshaw, Colin J. Burke, Andresa Rodrigues de Campos , et al. (75 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

The Properties of Fast Yellow Pulsating Supergiants: FYPS Point the Way to Missing Red Supergiants

Authors: Trevor Z. Dorn-Wallenstein, Emily M. Levesque, James R. A. Davenport, Kathryn F. Neugent, Brett M. Morris, K. Azalee Bostroem

Abstract: 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 M… ▽ More 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. △ Less

Submitted 16 November, 2022; v1 submitted 23 June, 2022; originally announced June 2022.

Comments: Consistent with published version which contains significantly improved detection and rejection of contaminant objects. Comments welcome

Searching the SETI Ellipsoid with Gaia

Authors: James R. A. Davenport, Barbara Cabrales, Sofia Sheikh, Steve Croft, Andrew P. V. Siemion, Daniel Giles, Ann Marie Cody

Abstract: The SETI Ellipsoid is a geometric method for prioritizing technosignature observations based on the strategy of receiving signals synchronized to conspicuous astronomical events. Precise distances to nearby stars from Gaia makes constraining Ellipsoid crossing times possible. Here we explore the utility of using the Gaia Catalog of Nearby Stars to select targets on the SN 1987A SETI Ellipsoid, as… ▽ More The SETI Ellipsoid is a geometric method for prioritizing technosignature observations based on the strategy of receiving signals synchronized to conspicuous astronomical events. Precise distances to nearby stars from Gaia makes constraining Ellipsoid crossing times possible. Here we explore the utility of using the Gaia Catalog of Nearby Stars to select targets on the SN 1987A SETI Ellipsoid, as well the Ellipsoids defined by 278 classical novae. Less than 8% of stars within the 100 pc sample are inside the SN 1987A SETI Ellipsoid, meaning the vast majority of nearby stars are still viable targets for monitoring over time. We find an average of 734 stars per year within the 100 pc volume will intersect the Ellipsoid from SN 1987A, with ~10% of those having distance uncertainties from Gaia better than 0.1 lyr. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 11 pages, 7 figures, Submitted to the AAS Journals

370 New Eclipsing Binary Candidates from TESS Sectors 1-26

Authors: Erin L. Howard, James R. A. Davenport, Kevin R. Covey

Abstract: We present 370 candidate eclipsing binaries (EBs), identified from ~510,000 short cadence TESS light curves. Our statistical criteria identify 5,105 light curves with features consistent with eclipses (~1% of the initial sample). After visual confirmation of the light curves, we have a final sample of 2,288 EB candidates. Among these, we find 370 sources that were not included in the catalog recen… ▽ More We present 370 candidate eclipsing binaries (EBs), identified from ~510,000 short cadence TESS light curves. Our statistical criteria identify 5,105 light curves with features consistent with eclipses (~1% of the initial sample). After visual confirmation of the light curves, we have a final sample of 2,288 EB candidates. Among these, we find 370 sources that were not included in the catalog recently published by Prsa et al. We publish our full sample of 370 new EB candidates, and statistical features used for their identification, reported per observation sector. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: Accepted by RNAAS. Table 1 is available here: https://doi.org/10.3847/2515-5172/ac6e42

Llamaradas Estelares: Modeling the Morphology of White-Light Flares

Authors: Guadalupe Tovar Mendoza, James R. A. Davenport, Eric Agol, James A. G. Jackman, Suzanne L. Hawley

Abstract: Stellar variability is a limiting factor for planet detection and characterization, particularly around active M-type stars. Here we revisit one of the most active stars from the Kepler mission, the M4 star GJ 1243, and use a sample of 414 flare events from 11 months of 1-minute cadence light curves to study the empirical morphology of white-light stellar flares. We use a Gaussian process detrendi… ▽ More Stellar variability is a limiting factor for planet detection and characterization, particularly around active M-type stars. Here we revisit one of the most active stars from the Kepler mission, the M4 star GJ 1243, and use a sample of 414 flare events from 11 months of 1-minute cadence light curves to study the empirical morphology of white-light stellar flares. We use a Gaussian process detrending technique to account for the underlying starspots. We present an improved analytic, continuous flare template that is generated by stacking the flares onto a scaled time and amplitude and uses a Markov Chain Monte Carlo analysis to fit the model. Our model is defined using classical flare events, but can also be used to model complex, multi-peaked flare events. We demonstrate the utility of our model using TESS data at the 10-minute, 2-minute and 20-second cadence modes. Our new flare model code is made publicly available on GitHub. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: 20 pages, 10 figures, accepted for publication in AJ. New flare model is available online https://github.com/lupitatovar/Llamaradas-Estelares

The Effect of Stellar Contamination on Low-resolution Transmission Spectroscopy: Needs Identified by NASA's Exoplanet Exploration Program Study Analysis Group 21

Authors: Benjamin V. Rackham, Néstor Espinoza, Svetlana V. Berdyugina, Heidi Korhonen, Ryan J. MacDonald, Benjamin T. Montet, Brett M. Morris, Mahmoudreza Oshagh, Alexander I. Shapiro, Yvonne C. Unruh, Elisa V. Quintana, Robert T. Zellem, Dániel Apai, Thomas Barclay, Joanna K. Barstow, Giovanni Bruno, Ludmila Carone, Sarah L. Casewell, Heather M. Cegla, Serena Criscuoli, Catherine Fischer, Damien Fournier, Mark S. Giampapa, Helen Giles, Aishwarya Iyer , et al. (36 additional authors not shown)

Abstract: 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 spectru… ▽ More 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. △ Less

Submitted 17 March, 2023; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 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

Simultaneous Multiwavelength Flare Observations of EV Lacertae

Authors: Rishi R. Paudel, Thomas Barclay, Joshua E. Schlieder, Elisa V. Quintana, Emily A. Gilbert, Laura D. Vega, Allison Youngblood, Michele Silverstein, Rachel A. Osten, Michael A. Tucker, Daniel Huber, Aaron Do, Kenji Hamaguchi, D. J. Mullan, John E. Gizis, Teresa A. Monsue, Knicole D. Colón, Patricia T. Boyd, James R. A. Davenport, Lucianne Walkowicz

Abstract: We present the first results of our ongoing project conducting simultaneous multiwavelength observations of flares on nearby active M dwarfs. We acquired data of the nearby dM3.5e star EV Lac using 5 different observatories: NASA's Transiting Exoplanet Survey Satellite (TESS), NASA's Neil Gehrels Swift Observatory (\textit{Swift}), NASA's Neutron Interior Composition Explorer (NICER), the Universi… ▽ More We present the first results of our ongoing project conducting simultaneous multiwavelength observations of flares on nearby active M dwarfs. We acquired data of the nearby dM3.5e star EV Lac using 5 different observatories: NASA's Transiting Exoplanet Survey Satellite (TESS), NASA's Neil Gehrels Swift Observatory (\textit{Swift}), NASA's Neutron Interior Composition Explorer (NICER), the University of Hawaii 2.2-m telescope (UH88) and the Las Cumbres Observatory Global Telescope (LCOGT) Network. During the $\sim$25 days of TESS observations, we acquired three simultaneous UV/X-ray observations using \textit{Swift} that total $\sim$18 ks, 21 simultaneous epochs totaling $\sim$98 ks of X-ray data using NICER, one observation ($\sim$ 3 hours) with UH88, and one observation ($\sim$ 3 hours) with LCOGT. We identified 56 flares in the TESS light curve with estimated energies in the range log $E_{\rm T}$ (erg) = (30.5 - 33.2), nine flares in the \textit{Swift} UVM2 light curve with estimated energies in the range log $E_{UV}$ (erg) = (29.3 - 31.1), 14 flares in the NICER light curve with estimated minimum energies in the range log $E_{N}$ (erg) = (30.5 - 32.3), and 1 flare in the LCOGT light curve with log $E_{L}$ (erg) = 31.6. We find that the flare frequency distributions (FFDs) of TESS and NICER flares have comparable slopes, $βべーた_{T}$ = -0.67$\pm$0.09 and $βべーた_{N}$ = -0.65$\pm$0.19, and the FFD of UVOT flares has a shallower slope ($βべーた_{U}$ = -0.38$\pm$0.13). Furthermore, we do not find conclusive evidence for either the first ionization potential (FIP) or the inverse FIP effect during coronal flares on EV Lac. △ Less

Submitted 10 August, 2021; originally announced August 2021.

Comments: 28 pages, 16 Figures, Accepted for publication in the Astrophysical Journal,

Giant white-light flares on fully convective stars occur at high latitudes

Authors: Ekaterina Ilin, Katja Poppenhaeger, Sarah J. Schmidt, Silva P. Järvinen, Elisabeth R. Newton, Julián D. Alvarado-Gómez, J. Sebastian Pineda, James R. A. Davenport, Mahmoudreza Oshagh, Ilya Ilyin

Abstract: White-light flares are magnetically driven localized brightenings on the surfaces of stars. Their temporal, spectral, and statistical properties present a treasury of physical information about stellar magnetic fields. The spatial distributions of magnetic spots and associated flaring regions help constrain dynamo theories. Moreover, flares are thought to crucially affect the habitability of exopl… ▽ More White-light flares are magnetically driven localized brightenings on the surfaces of stars. Their temporal, spectral, and statistical properties present a treasury of physical information about stellar magnetic fields. The spatial distributions of magnetic spots and associated flaring regions help constrain dynamo theories. Moreover, flares are thought to crucially affect the habitability of exoplanets that orbit these stars. Measuring the location of flares on stars other than the Sun is challenging due to the lack of spatial resolution. Here we present four fully convective stars observed with the Transiting Exoplanet Survey Satellite (TESS) that displayed large, long-duration flares in white-light which were modulated in brightness by the stars' fast rotation. This allowed us to determine the loci of these flares directly from the light curves. All four flares occurred at latitudes between 55 deg and 81 deg, far higher than typical solar flare latitudes. Our findings are evidence that strong magnetic fields tend to emerge close to the stellar rotational poles for fully convective stars, and suggest that the impact of flares on the habitability of exoplanets around small stars could be weaker than previously thought. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 23 pages, 18 figures, accepted to MNRAS, see conference poster https://doi.org/10.5281/zenodo.4558791 for short summary

The Rise and Fall of the Eclipsing Binary HS Hydrae

Authors: James R. A. Davenport, Diana Windemuth, Karen Warmbein, Erin L. Howard, Courtney Klein, Jessica Birky

Abstract: HS Hydrae is a short period eclipsing binary (P_orb=1.57 day) that belongs to a rare group of systems observed to have rapidly changing inclinations. This evolution is due to a third star on an intermediate orbit, and results in significant differences in eclipse depths and timings year-to-year. Zasche & Paschke (2012) revealed that HS Hydrae's eclipses were rapidly fading from view, predicting th… ▽ More HS Hydrae is a short period eclipsing binary (P_orb=1.57 day) that belongs to a rare group of systems observed to have rapidly changing inclinations. This evolution is due to a third star on an intermediate orbit, and results in significant differences in eclipse depths and timings year-to-year. Zasche & Paschke (2012) revealed that HS Hydrae's eclipses were rapidly fading from view, predicting they would cease around 2022. Using 25 days of photometric data from Sector 009 of the Transiting Exoplanet Survey Satellite (TESS), we find that the primary eclipses for HS Hydrae were only 0.00173+/-0.00007 mag in depth in March 2019. This data from TESS likely represents the last eclipses detected from HS Hydrae. We also searched the Digitization of the Harvard Astronomical Plate Collection (DASCH) archive for historic data from the system. With a total baseline of over 125 years, this unique combination of data sets - from photographic plates to precision space-based photometry - allows us to trace the emergence and decay of eclipses from HS Hydrae, and further constrain its evolution. Recent TESS observations from Sector 035 confirm that eclipses have ceased for HS Hya, and we estimate they will begin again in 2195. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 13 pages, 5 figures, AJ Accepted

How to organize an online conference -- Lessons learned from Cool Stars 20.5 (virtually cool)

Authors: Hans Moritz Günther, James R. A. Davenport, Scott Wolk, Shaun Gallagher

Abstract: The virtual meeting was a success. Several people told us that this was "the best virtual meeting they had seen so far", which, a year into the pandemic and without a commercial provider in the back, is a great success. The biggest point of criticism was the timing: We had programming from UTC 17:00-22:00 (evening and night in central Europe, afternoon on the US East Coast, during the day in South… ▽ More The virtual meeting was a success. Several people told us that this was "the best virtual meeting they had seen so far", which, a year into the pandemic and without a commercial provider in the back, is a great success. The biggest point of criticism was the timing: We had programming from UTC 17:00-22:00 (evening and night in central Europe, afternoon on the US East Coast, during the day in South America and on the US West coast, but in the middle of the night for Asia and Australia). There is no good solution, but at least some variation in session time might go a long way to make it easier for all to attend at least some sessions. Feedback also indicates that the schedule was too compressed. Poster sessions and social contacts with the tool Gathertown worked out really well for all that used it. Our way of combining several services (Zoom for plenary and break-out rooms, Zenodo for uploading and viewing posters and proceedings, Google forms for registration and abstract submission, gathertown) allowed for a very low-cost meeting with little overhead (total cost: 600 $ for gathertown, zoom was provided through an institutional subscription, just 4 people on the LOC). △ Less

Submitted 18 May, 2021; originally announced May 2021.

Photometric Classifications of Evolved Massive Stars: Preparing for the Era of Webb and Roman with Machine Learning

Authors: Trevor Z. Dorn-Wallenstein, James R. A. Davenport, Daniela Huppenkothen, Emily M. Levesque

Abstract: In the coming years, next-generation space-based infrared observatories will significantly increase our samples of rare massive stars, representing a tremendous opportunity to leverage modern statistical tools and methods to test massive stellar evolution in entirely new environments. Such work is only possible if the observed objects can be reliably classified. Spectroscopic observations are infe… ▽ More In the coming years, next-generation space-based infrared observatories will significantly increase our samples of rare massive stars, representing a tremendous opportunity to leverage modern statistical tools and methods to test massive stellar evolution in entirely new environments. Such work is only possible if the observed objects can be reliably classified. Spectroscopic observations are infeasible with more distant targets, and so we wish to determine whether machine learning methods can classify massive stars using broadband infrared photometry. We find that a Support Vector Machine classifier is capable of coarsely classifying massive stars with labels corresponding to hot, cool, and emission line stars with high accuracy, while rejecting contaminating low mass giants. Remarkably, 76\% of emission line stars can be recovered without the need for narrowband or spectroscopic observations. We classify a sample of ${\sim}2500$ objects with no existing labels, and identify fourteen candidate emission line objects. Unfortunately, despite the high precision of the photometry in our sample, the heterogeneous origins of the labels for the stars in our sample severely inhibits our classifier from distinguishing classes of stars with more granularity. Ultimately, no large and homogeneously labeled sample of massive stars currently exists. Without significant efforts to robustly classify evolved massive stars -- which is feasible given existing data from large all-sky spectroscopic surveys -- shortcomings in the labeling of existing data sets will hinder efforts to leverage the next-generation of space observatories. △ Less

Submitted 22 March, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 23 pages, 16 figures, 6 tables. Accepted for publication in ApJ. Comments welcome

Stellar Rotation in the K2 Sample: Evidence for Modified Spindown

Authors: Tyler A. Gordon, James R. A. Davenport, Ruth Angus, Daniel Foreman-Mackey, Eric Agol, Kevin R. Covey, Marcel Agüeros, David Kipping

Abstract: We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap… ▽ More We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap in the rotation period-color diagram first reported by \cite{McQuillan2013}. While the gap was tentatively detected in \cite{Reinhold2020}, this work represents the first robust detection of the gap in K2 data for field stars. This is significant because K2 observed along many lines of sight at wide angular separation, in contrast to Kepler's single line of sight. Together with recent results for rotation in open clusters, we interpret this gap as evidence for a departure from the $t^{-1/2}$ Skumanich spin down law, rather than an indication of a bimodal star formation history. We provide maximum likelihood estimates and uncertainties for all parameters of the quasi-periodic light curve model for each of the 284,834 stars in our sample. △ Less

Submitted 1 April, 2021; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: 17 pages, 14 figures, submitted to ApJ, revised April 1, 2020

Flares in Open Clusters with K2. II. Pleiades, Hyades, Praesepe, Ruprecht 147, and M67

Authors: Ekaterina Ilin, Sarah J. Schmidt, Katja Poppenhäger, James R. A. Davenport, Martti H. Kristiansen, Mark Omohundro

Abstract: Flares, energetic eruptions on the surfaces of stars, are an unmistakable manifestation of magnetically driven emission. Their occurrence rates and energy distributions trace stellar characteristics such as mass and age. But before flares can be used to constrain stellar properties, the flaring-age-mass relation requires proper calibration. This work sets out to quantify flaring activity of indepe… ▽ More Flares, energetic eruptions on the surfaces of stars, are an unmistakable manifestation of magnetically driven emission. Their occurrence rates and energy distributions trace stellar characteristics such as mass and age. But before flares can be used to constrain stellar properties, the flaring-age-mass relation requires proper calibration. This work sets out to quantify flaring activity of independently age-dated main sequence stars for a broad range of spectral types using optical light curves obtained by the Kepler satellite. Drawing from the complete K2 archive, we searched 3435 $\sim 80$ day long light curves of 2111 open cluster members for flares using the open-source software packages K2SC to remove instrumental and astrophysical variability from K2 light curves, and AltaiPony to search and characterize the flare candidates. We confirmed a total of 3844 flares on high probability open cluster members with ages from zero age main sequence (Pleiades) to 3.6 Gyr (M67). We extended the mass range probed in the first study of this series to span from Sun-like stars to mid-M dwarfs. We added the Hyades (690 Myr) to the sample as a comparison cluster to Praesepe (750 Myr), the 2.6 Gyr old Ruprecht 147, and several hundred light curves from the late K2 Campaigns in the remaining clusters. The flare energy distribution was similar in the entire parameter space, following a power law relation with exponent $αあるふぁ\approx 1.84-2.39$. The flaring rates declined with age, and declined faster for higher mass stars. We found evidence that a rapid decline in flaring activity occurred in M1-M2 dwarfs around Hyades/Praesepe age, when these stars spun down to rotation periods of about 10 days, while higher mass stars had already transitioned to lower flaring rates, and lower mass stars still resided in the saturated activity regime. (abridged) △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 27 pages, 13 figures. Accepted to A&A

Journal ref: A&A 645, A42 (2021)

Short Term Variability of Evolved Massive Stars with TESS II: A New Class of Cool, Pulsating Supergiants

Authors: Trevor Z. Dorn-Wallenstein, Emily M. Levesque, Kathryn F. Neugent, James R. A. Davenport, Brett M. Morris, Keyan Gootkin

Abstract: 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… ▽ More 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 $αあるふぁ$ Cygni variables in our sample, of which three are newly discovered. △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: 29 pages, 13 figures, 8 tables. Accepted for publication in ApJ. Comments welcome

Deconfusing the Confusogram: Getting New Insights from Zeeman Doppler Imaging

Authors: J. Sebastian Pineda, James R. A. Davenport

Abstract: Zeeman Doppler Imaging is a powerful tool for characterizing the strength and topology of stellar magnetic fields. In this research note, we present a new way to visualize the typical results from ZDI for an ensemble of stars, addressing some of the concerns with the standard `confusogram' approach to illustrating the data. Our publically available plotting methods further enable an accessible mea… ▽ More Zeeman Doppler Imaging is a powerful tool for characterizing the strength and topology of stellar magnetic fields. In this research note, we present a new way to visualize the typical results from ZDI for an ensemble of stars, addressing some of the concerns with the standard `confusogram' approach to illustrating the data. Our publically available plotting methods further enable an accessible means to consider variability in the inferred magnetic field topologies from repeated observations, as we demonstrate with the literature ZDI data on M dwarfs. △ Less

Submitted 5 June, 2020; originally announced June 2020.

Comments: Published in RNAAS; Code publicly available

10 Years of Stellar Activity for GJ 1243

Authors: James. R. A. Davenport, Guadalupe Tovar Mendoza, Suzanne L. Hawley

Abstract: The flaring M4 dwarf GJ 1243 has become a benchmark for studying stellar flare and starspot activity thanks to the exceptional photometric monitoring archive from the Kepler mission. New light curves from the TESS mission for this star allow precise stellar activity characterization over more than a decade timescale. We have carried out the first flare and starspot analysis of GJ 1243 from over 50… ▽ More The flaring M4 dwarf GJ 1243 has become a benchmark for studying stellar flare and starspot activity thanks to the exceptional photometric monitoring archive from the Kepler mission. New light curves from the TESS mission for this star allow precise stellar activity characterization over more than a decade timescale. We have carried out the first flare and starspot analysis of GJ 1243 from over 50 days of data from TESS Sectors 14 and 15. Using 133 flare events detected in the 2-minute cadence TESS data, we compare the cumulative flare frequency distributions, and find the flare activity for GJ 1243 is unchanged between the Kepler and TESS epochs. Two distinct starspot groups are found in the TESS data, with the primary spot having the same rotational period and phase as seen in Kepler. The phase of the secondary spot feature is consistent with the predicted location of the secondary starspot and measurement of weak differential rotation, suggesting this secondary spot may be long-lived and stable in both latitude and longitude. As expected for this highly active star, the constant spot and flare activity reveal no sign of solar-like activity cycles over 10 years. However, we highlight the unique ability for Kepler and TESS to use flare rates to detect activity cycles. △ Less

Submitted 20 May, 2020; originally announced May 2020.

Comments: 14 pages, 5 figures, AJ accepted

Temporal Evolution of Spatially-Resolved Individual Star Spots on a Planet-Hosting Solar-type Star: Kepler 17

Authors: Kosuke Namekata, James R. A. Davenport, Brett M. Morris, Suzanne L. Hawley, Hiroyuki Maehara, Yuta Notsu, Shin Toriumi, Kai Ikuta, Shota Notsu, Satoshi Honda, Daisaku Nogami, Kazunari Shibata

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 37 pages, 12 figures, 1 table. Accepted for publication in The Astrophysical Journal

A blueprint of state-of-the-art techniques for detecting quasi-periodic pulsations in solar and stellar flares

Authors: Anne-Marie Broomhall, James R. A. Davenport, Laura A. Hayes, Andrew R. Inglis, Dmitrii Y. Kolotkov, James A. McLaughlin, Tishtrya Mehta, Valery M. Nakariakov, Yuta Notsu, David J. Pascoe, Chloe E. Pugh, Tom Van Doorsselaere

Abstract: Quasi-periodic pulsations (QPPs) appear to be a common feature observed in the light curves of both solar and stellar flares. However, their quasi-periodic nature, along with the fact that they can be small in amplitude and short-lived, makes QPPs difficult to unequivocally detect. In this paper, we test the strengths and limitations of state-of-the-art methods for detecting QPPs using a series of… ▽ More Quasi-periodic pulsations (QPPs) appear to be a common feature observed in the light curves of both solar and stellar flares. However, their quasi-periodic nature, along with the fact that they can be small in amplitude and short-lived, makes QPPs difficult to unequivocally detect. In this paper, we test the strengths and limitations of state-of-the-art methods for detecting QPPs using a series of hare-and-hounds exercises. The hare simulated a set of flares, both with and without QPPs of a variety of forms, while the hounds attempted to detect QPPs in blind tests. We use the results of these exercises to create a blueprint for anyone who wishes to detect QPPs in real solar and stellar data. We present eight clear recommendations to be kept in mind for future QPP detections, with the plethora of solar and stellar flare data from new and future satellites. These recommendations address the key pitfalls in QPP detection, including detrending, trimming data, accounting for colored noise, detecting stationary-period QPPs, detecting QPPs with nonstationary periods, and ensuring that detections are robust and false detections are minimized. We find that QPPs can be detected reliably and robustly by a variety of methods, which are clearly identified and described, if the appropriate care and due diligence are taken. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: Accepted for publication in ApJSS, 46 pages, 29 figures, 7 tables

High Fidelity Imaging of the Inner AUえーゆー Mic Debris Disk: Evidence of Differential Wind Sculpting?

Authors: John P. Wisniewski, Adam F. Kowalski, James R. A. Davenport, Glenn Schneider, Carol A. Grady, Leslie Hebb, Kellen D. Lawson, Jean-Charles Augereau, Anthony Boccaletti, Alexander Brown, John H. Debes, Andras Gaspar, Thomas K. Henning, Dean C. Hines, Marc J. Kuchner, Anne-Marie Lagrange, Julien Milli, Elie Sezestre, Christopher C. Stark, Christian Thalmann

Abstract: We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AUえーゆー Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morpholog… ▽ More We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AUえーゆー Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morphology at the time of our observations. The loop has a projected width of 1.5 au and rises 2.3 au above the disk midplane. We also explored TESS photometric observations of AUえーゆー Mic that are consistent with evidence of two starspot complexes in the system. The likely co-alignment of the stellar and disk rotational axes breaks degeneracies in detailed spot modeling, indicating that AUえーゆー Mic's projected magnetic field axis is offset from its rotational axis. We speculate that small grains in AUえーゆー Mic's disk could be sculpted by a time-dependent wind that is influenced by this offset magnetic field axis, analogous to co-rotating Solar interaction regions that sculpt and influence the inner and outer regions of our own Heliosphere. Alternatively, if the observed spot modulation is indicative of a significant mis-alignment of the stellar and disk rotational axes, we suggest the disk could still be sculpted by the differential equatorial versus polar wind that it sees with every stellar rotation. △ Less

Submitted 2 September, 2019; v1 submitted 23 July, 2019; originally announced July 2019.

Comments: Accepted to ApJL

SETI in the Spatio-Temporal Survey Domain

Authors: James. R. A. Davenport

Abstract: Traditional searches for extraterrestrial intelligence (SETI) or "technosignatures" focus on dedicated observations of single stars or regions in the sky to detect excess or transient emission from intelligent sources. The newest generation of synoptic time domain surveys enable an entirely new approach: spatio-temporal SETI, where technosignatures may be discovered from spatially resolved sources… ▽ More Traditional searches for extraterrestrial intelligence (SETI) or "technosignatures" focus on dedicated observations of single stars or regions in the sky to detect excess or transient emission from intelligent sources. The newest generation of synoptic time domain surveys enable an entirely new approach: spatio-temporal SETI, where technosignatures may be discovered from spatially resolved sources or multiple stars over time. Current optical time domain surveys such as ZTF and the Evryscope can probe 10-100 times more of the "Cosmic Haystack" parameter space volume than many radio SETI investigations. Small-aperture, high cadence surveys like Evryscope can be comparable in their Haystack volume completeness to deeper surveys including LSST. Investigations with these surveys can also be conducted at a fraction of the cost of dedicated SETI surveys, since they make use of data already being gathered. However, SETI methodology has not widely utilized such surveys, and the field is in need of new search algorithms that can account for signals in both the spatial and temporal domains. Here I describe the broad potential for modern wide-field time domain optical surveys to revolutionize our search for technosignatures, and illustrate some example SETI approaches using transiting exoplanets to form a distributed beacon. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 8 pages, 3 figures, AJ submitted

Do Kepler superflare stars really include slowly-rotating Sun-like stars ? - Results using APO 3.5m telescope spectroscopic observations and Gaia-DR2 data -

Authors: Yuta Notsu, Hiroyuki Maehara, Satoshi Honda, Suzanne L. Hawley, James R. A. Davenport, Kosuke Namekata, Shota Notsu, Kai Ikuta, Daisaku Nogami, Kazunari Shibata

Abstract: We report the latest view of Kepler solar-type (G-type main-sequence) superflare stars, including recent updates with Apache Point Observatory (APO) 3.5m telescope spectroscopic observations and Gaia-DR2 data. First, we newly conducted APO3.5m spectroscopic observations of 18 superflare stars found from Kepler 1-min time cadence data. More than half (43 stars) are confirmed to be "single" stars, a… ▽ More We report the latest view of Kepler solar-type (G-type main-sequence) superflare stars, including recent updates with Apache Point Observatory (APO) 3.5m telescope spectroscopic observations and Gaia-DR2 data. First, we newly conducted APO3.5m spectroscopic observations of 18 superflare stars found from Kepler 1-min time cadence data. More than half (43 stars) are confirmed to be "single" stars, among 64 superflare stars in total that have been spectroscopically investigated so far in this APO3.5m and our previous Subaru/HDS observations. The measurements of $v\sin i$ (projected rotational velocity) and chromospheric lines (Ca II H\&K and Ca II 8542Å) support the brightness variation of superflare stars is caused by the rotation of a star with large starspots. We then investigated the statistical properties of Kepler solar-type superflare stars by incorporating Gaia-DR2 stellar radius estimates. As a result, the maximum superflare energy continuously decreases as the rotation period $P_{\mathrm{rot}}$ increases. Superflares with energies $\lesssim 5\times10^{34}$ erg occur on old, slowly-rotating Sun-like stars ($P_{\mathrm{rot}}\sim$25 days) approximately once every 2000--3000 years, while young rapidly-rotating stars with $P_{\mathrm{rot}}\sim$ a few days have superflares up to $10^{36}$ erg. The maximum starspot area does not depend on the rotation period when the star is young, but as the rotation slows down, it starts to steeply decrease at $P_{\mathrm{rot}}\gtrsim$12 days for Sun-like stars. These two decreasing trends are consistent since the magnetic energy stored around starspots explains the flare energy, but other factors like spot magnetic structure should also be considered. △ Less

Submitted 29 March, 2019; originally announced April 2019.

Comments: 71 pages, 31 figures, 10 tables. Accepted for publication in The Astrophysical Journal (on March 29, 2019)

Rotation Period Evolution in Low-Mass Binary Stars: The Impact of Tidal Torques and Magnetic Braking

Authors: David P. Fleming, Rory Barnes, James R. A. Davenport, Rodrigo Luger

Abstract: We examine how tides, stellar evolution, and magnetic braking shape the rotation period (P$_{rot}$) evolution of low-mass stellar binaries up to orbital periods (P$_{orb}$) of 100 d across a wide range tidal dissipation parameters using two common equilibrium tidal models. We find that many binaries with P$_{orb} < 20$ d tidally lock, and most with $P_{orb} < 4$ d tidally lock into synchronous rot… ▽ More We examine how tides, stellar evolution, and magnetic braking shape the rotation period (P$_{rot}$) evolution of low-mass stellar binaries up to orbital periods (P$_{orb}$) of 100 d across a wide range tidal dissipation parameters using two common equilibrium tidal models. We find that many binaries with P$_{orb} < 20$ d tidally lock, and most with $P_{orb} < 4$ d tidally lock into synchronous rotation on circularized orbits. At short P$_{orb}$, tidal torques produce a population of fast rotators that single-star only models of magnetic braking fail to produce. In many cases, we show that the competition between magnetic braking and tides produces a population of subsynchronous rotators that persists for Gyrs, even in short P$_{orb}$ binaries, qualitatively reproducing the subsynchronous eclipsing binaries (EBs) discovered in the Kepler field by Lurie et al. (2017). Both equilibrium tidal models predict that binaries can tidally-interact out to P$_{orb} \approx 80$ d, while the Constant Phase Lag tidal model predicts that binaries can tidally lock out to P$_{orb} \approx 100$ d. Tidal torques often force the P$_{rot}$ evolution of stellar binaries to depart from the long-term magnetic braking-driven spin down experienced by single stars, revealing that P$_{rot}$ is not be a valid proxy for age in all cases, i.e. gyrochronology can underpredict ages by up to $300\%$ unless one accounts for binarity. We suggest that accurate determinations of orbital eccentricties and P$_{rot}$ can be used to discriminate between which equilibrium tidal models best describes tidal interactions in low-mass binary stars. △ Less

Submitted 1 July, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: Accepted, ApJ

Short Term Variability of Evolved Massive Stars with TESS

Authors: Trevor Z. Dorn-Wallenstein, Emily M. Levesque, James R. A. Davenport

Abstract: We present the first results from a study of TESS Sector 1 and 2 light curves for eight evolved massive stars in the LMC: six yellow supergiants (YSGs) and two luminous blue variables (LBVs), including S Doradus. We use an iterative prewhitening procedure to characterize the short-timescale variability in all eight stars. The periodogram of one of the YSGs, HD 269953, displays multiple strong peak… ▽ More We present the first results from a study of TESS Sector 1 and 2 light curves for eight evolved massive stars in the LMC: six yellow supergiants (YSGs) and two luminous blue variables (LBVs), including S Doradus. We use an iterative prewhitening procedure to characterize the short-timescale variability in all eight stars. The periodogram of one of the YSGs, HD 269953, displays multiple strong peaks at higher frequencies than its fellows. While the field surrounding HD 269953 is quite crowded, it is the brightest star in the region, and has infrared colors indicating it is dusty. We suggest HD 269953 may be in a post-red supergiant evolutionary phase. We find a signal with a period of $\sim5$ days for the LBV HD 269582. The periodogram of S Doradus shows a complicated structure, with peaks below frequencies of 1.5 cycles per day. We fit the shape of the background noise of all eight light curves, and find a red noise component in all of them. However, the power law slope of the red noise and the timescale over which coherent structures arise changes from star to star. Our results highlight the potential for studying evolved massive stars with TESS. △ Less

Submitted 20 May, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 24 pages, 11 figures, 12 Tables. Accepted, ApJ. Comments welcome

The Solar Benchmark: Rotational Modulation of the Sun Reconstructed from Archival Sunspot Records

Authors: Brett M. Morris, James R. A. Davenport, Helen A. C. Giles, Leslie Hebb, Suzanne L. Hawley, Ruth Angus, Peter A. Gilman, Eric Agol

Abstract: 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… ▽ More 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. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

The Evolution of Flare Activity with Stellar Age

Authors: James R. A. Davenport, Kevin R. Covey, Riley W. Clarke, Austin C. Boeck, Jonathan Cornet, Suzanne L. Hawley

Abstract: Using a recent census of flare stars from the Kepler survey, we have explored how flare activity evolves across stellar main sequence lifetimes. We utilize a sample of 347 stars with robust flare activity detections, and which have rotation periods measured via starspot modulations in their Kepler light curves. We consider three separate methods for quantifying flare activity from optical light cu… ▽ More Using a recent census of flare stars from the Kepler survey, we have explored how flare activity evolves across stellar main sequence lifetimes. We utilize a sample of 347 stars with robust flare activity detections, and which have rotation periods measured via starspot modulations in their Kepler light curves. We consider three separate methods for quantifying flare activity from optical light curves, and compare their utility for comparing flare activity between stars of differing ages and luminosities. These metrics include: the fractional luminosity emitted in flares, the specific rate of flares emitted at a given energy, and a model for the entire flare frequency distribution. With all three approaches we find that flare activity decreases for all low-mass stars as they spin-down, and thus with age. Most striking is the evolution of the flare occurrence frequency distributions, which show no significant change in the power law slope with age. Since our sample is preferentially constructed of younger, more active stars, our model over-predicts the super-flare rate previously estimated for the Sun. Finally, we parameterize our best-fit model of the flare frequency distribution for ease in predicting the rates of flares and their associated impacts on planet habitability and detection. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 21 pages, 11 figures, ApJ Accepted 2018-Dec-31

Flares in Open Clusters with K2. I. M45 (Pleiades), M44 (Praesepe) and M67

Authors: Ekaterina Ilin, Sarah J. Schmidt, James R. A. Davenport, Klaus G. Strassmeier

Abstract: The presence and strength of a stellar magnetic field and activity is rooted in a star's fundamental parameters such as mass and age. Can flares serve as an accurate stellar "clock"? To explore if we can quantify an activity-age relation in the form of a flaring-age relation, we measured trends in the flaring rates and energies for stars with different masses and ages. We investigated the time… ▽ More The presence and strength of a stellar magnetic field and activity is rooted in a star's fundamental parameters such as mass and age. Can flares serve as an accurate stellar "clock"? To explore if we can quantify an activity-age relation in the form of a flaring-age relation, we measured trends in the flaring rates and energies for stars with different masses and ages. We investigated the time-domain photometry provided by Kepler's follow-up mission K2 and searched for flares in three solar metallicity open clusters with well-known ages, M45 (0.125 Gyr), M44 (0.63 Gyr), and M67 (4.3 Gyr). We updated and employed the automated flare finding and analysis pipeline Appaloosa, originally designed for Kepler. We introduced a synthetic flare injection and recovery subroutine to ascribe detection and energy recovery rates for flares in a broad energy range for each light curve. We collected a sample of 1 761 stars, mostly late-K to mid-M dwarfs and found 751 flare candidates with energies ranging from $4\cdot10^{32}$ erg to $6\cdot10^{34}$ erg, of which 596 belong to M45, 155 to M44, and none to M67. We find that flaring activity depends both on $T_\mathrm{eff}$, and age. But all flare frequency distributions have similar slopes with $αあるふぁ\approx2.0-2.4$, supporting a universal flare generation process. We discuss implications for the physical conditions under which flares occur, and how the sample's metallicity and multiplicity affect our results. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 17 pages, 11 figures, appendix. Accepted to A&A

Journal ref: A&A 622, A133 (2019)

A Significant Over-Luminosity in the Transiting Brown Dwarf CWW 89Ab

Authors: Thomas G. Beatty, Caroline V. Morley, Jason L. Curtis, Adam Burrows, James R. A. Davenport, Benjamin T. Montet

Abstract: We observed eclipses of the transiting brown dwarf CWW 89Ab at 3.6um and 4.5um using Spitzer/IRAC. The CWW 89 binary system is a member of the $3.0\pm0.25$ Gyr-old open cluster Ruprecht 147, and is composed of a Sun-like primary and an early M-dwarf secondary separated by a projected distance of 25 AUえーゆー. CWW 89Ab has a radius of $0.937\pm0.042$ RJ and a mass of $36.5\pm0.1$ MJ, and is on a 5.3 day o… ▽ More We observed eclipses of the transiting brown dwarf CWW 89Ab at 3.6um and 4.5um using Spitzer/IRAC. The CWW 89 binary system is a member of the $3.0\pm0.25$ Gyr-old open cluster Ruprecht 147, and is composed of a Sun-like primary and an early M-dwarf secondary separated by a projected distance of 25 AUえーゆー. CWW 89Ab has a radius of $0.937\pm0.042$ RJ and a mass of $36.5\pm0.1$ MJ, and is on a 5.3 day orbit about CWW 89A with a non-zero eccentricity of $e=0.19$ (Curtis et al. 2016). We strongly detect the eclipses of CWW 89Ab in both Spitzer channels as $δでるた_{3.6}=1147\pm213$ ppm and $δでるた_{4.5}=1097\pm225$ ppm after correcting for the dilution from CWW 89B. After accounting for the irradiation that CWW 89Ab receives from its host star, these measurements imply that the brown dwarf has an internal luminosity of $\log(\mathrm{L_{bol}/\mathrm{L}_\odot})=-4.19\pm0.14$. This is 16 times, or $9.3\,σしぐま$, higher than model predictions given the known mass, radius, and age of CWW 89Ab. As we discuss, this over-luminosity is not explainable by an inaccurate age determination, additional stellar heating, nor tidal heating. Instead, we suggest that the anomalous luminosity of CWW 89Ab is caused by a dayside temperature inversion -- though a significant error in the evolutionary models is also a possibility. Importantly, a temperature inversion would require a super-stellar C/O ratio in CWW 89Ab's atmosphere. If this is indeed the case, it implies that CWW 89Ab is a 36.5 MJ object that formed via core accretion processes. Finally, we use our measurement of CWW 89Ab's orbital eccentricity, improved via these observations, to constrain the tidal quality factors of the brown dwarf and the host star CWW 89A to be $Q_{BD}>10^{4.15}$ and $Q_*>10^{9}$, respectively. △ Less

Submitted 24 September, 2018; v1 submitted 30 July, 2018; originally announced July 2018.

Comments: 19 pages, 12 figures. AJ in press. Updated based on comments received

Rotating Stars from Kepler Observed with Gaia DR2

Authors: James R. A. Davenport, Kevin R. Covey

Abstract: We have matched the astrometric data from Gaia Data Release 2 to the sample of stars with measured rotation periods from Kepler. Using 30,305 stars with good distance estimates, we select 16,248 as being likely main sequence single stars centered within a 0.5 mag region about a 1 Gyr isochrone, removing many sub-giants and unresolved binary stars from the sample. The rotation period bimodality, or… ▽ More We have matched the astrometric data from Gaia Data Release 2 to the sample of stars with measured rotation periods from Kepler. Using 30,305 stars with good distance estimates, we select 16,248 as being likely main sequence single stars centered within a 0.5 mag region about a 1 Gyr isochrone, removing many sub-giants and unresolved binary stars from the sample. The rotation period bimodality, originally discovered by McQuillan et al. (2013), is clearly recovered for stars out to 525pc, but is not detectable at further distances. This bimodality correlates with Galactic height as well, dropping strongly for stars above Z>90 pc. We also find a significant width in the stellar main sequence of $ΔでるたM_G$ ~0.25 mag, as well as a coherent gradient of increasing rotation periods orthogonal to the main sequence. We interpret this as a signature of stellar angular momentum loss over time, implying a corresponding diagonal age gradient across the main sequence. Stellar evolution models predict changes in color and luminosity that are consistent in amplitude, but not in direction, with those required to produce the gradient we have detected. This rotation gradient suggests that main sequence evolution produces offsets in color-magnitude space that are significantly more orthogonal to the zero-age main sequence than models currently predict, and may provide new tests for both stellar evolution and gyrochronology models. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: ApJ Submitted

Possible Bright Starspots on TRAPPIST-1

Authors: Brett M. Morris, Eric Agol, James R. A. Davenport, Suzanne L. Hawley

Abstract: 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 $μみゅー$m. If the rotational modulation is due to starspots, persiste… ▽ More 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 $μみゅー$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. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: Accepted by ApJ March 12, 2018

Spotting Stellar Activity Cycles in Gaia Astrometry

Authors: Brett M. Morris, Eric Agol, James R. A. Davenport, Suzanne L. Hawley

Abstract: 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… ▽ More 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. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS, Feb 27, 2018

The First Post-Kepler Brightness Dips of KIC 8462852

Authors: Tabetha S. Boyajian, Roi Alonso, Alex Ammerman, David Armstrong, A. Asensio Ramos, K. Barkaoui, Thomas G. Beatty, Z. Benkhaldoun, Paul Benni, Rory Bentley, Andrei Berdyugin, Svetlana Berdyugina, Serge Bergeron, Allyson Bieryla, Michaela G. Blain, Alicia Capetillo Blanco, Eva H. L. Bodman, Anne Boucher, Mark Bradley, Stephen M. Brincat, Thomas G. Brink, John Briol, David J. A. Brown, J. Budaj, A. Burdanov , et al. (181 additional authors not shown)

Abstract: 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 "Els… ▽ More 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. △ Less

Submitted 2 January, 2018; originally announced January 2018.

Comments: 19 pages, 8 figures, accepted for publication in ApJL

The GALEX View of "Boyajian's Star" (KIC 8462852)

Authors: James. R. A. Davenport, Kevin R. Covey, Riley W. Clarke, Zachery Laycock, Scott W. Fleming, Tabetha S. Boyajian, Benjamin T. Montet, Bernie Shiao, Chase C. Million, David J. Wilson, Manuel Olmedo, Eric E. Mamajek, Daniel Olmedo, Miguel Chavez, Emanuele Bertone

Abstract: The enigmatic star KIC 8462852, informally known as "Boyajian's Star", has exhibited unexplained variability from both short timescale (days) dimming events, and years-long fading in the Kepler mission. No single physical mechanism has successfully explained these observations to date. Here we investigate the ultraviolet variability of KIC 8462852 on a range of timescales using data from the GALEX… ▽ More The enigmatic star KIC 8462852, informally known as "Boyajian's Star", has exhibited unexplained variability from both short timescale (days) dimming events, and years-long fading in the Kepler mission. No single physical mechanism has successfully explained these observations to date. Here we investigate the ultraviolet variability of KIC 8462852 on a range of timescales using data from the GALEX mission that occurred contemporaneously with the Kepler mission. The wide wavelength baseline between the Kepler and GALEX data provides a unique constraint on the nature of the variability. Using 1600 seconds of photon-counting data from four GALEX visits spread over 70 days in 2011, we find no coherent NUV variability in the system on 10-100 second or months timescales. Comparing the integrated flux from these 2011 visits to the 2012 NUV flux published in the GALEX-CAUSE Kepler survey, we find a 3% decrease in brightness for KIC 8462852. We find this level of variability is significant, but not necessarily unusual for stars of similar spectral type in the GALEX data. This decrease coincides with the secular optical fading reported by Montet & Simon (2016). We find the multi-wavelength variability is somewhat inconsistent with typical interstellar dust absorption, but instead favors a R$_V$ = 5.0 $\pm$ 0.9 reddening law potentially from circumstellar dust. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: 8 pages, 4 figures, ApJ Accepted

Flare Activity of Wide Binary Stars with Kepler

Authors: Riley W. Clarke, James R. A. Davenport, Kevin R. Covey, Christoph Baranec

Abstract: We present an analysis of flare activity in wide binary stars using a combination of value-added data sets from the NASA Kepler mission. The target list contains a set of previously discovered wide binary star systems identified by proper motions in the Kepler field. We cross-matched these systems with estimates of flare activity for $\sim$200,000 stars in the Kepler field, allowing us to compare… ▽ More We present an analysis of flare activity in wide binary stars using a combination of value-added data sets from the NASA Kepler mission. The target list contains a set of previously discovered wide binary star systems identified by proper motions in the Kepler field. We cross-matched these systems with estimates of flare activity for $\sim$200,000 stars in the Kepler field, allowing us to compare relative flare luminosity between stars in coeval binaries. From a sample of 184 previously known wide binaries in the Kepler field, we find 58 with detectable flare activity in at least 1 component, 33 of which are similar in mass (q $>$ 0.8). Of these 33 equal-mass binaries, the majority display similar ($\pm 1$ dex) flare luminosity between both stars, as expected for stars of equal mass and age. However, we find two equal-mass pairs where the secondary (lower mass) star is more active than its counterpart, and two equal-mass pairs where the primary star is more active. The stellar rotation periods are also anomalously fast for stars with elevated flare activity. Pairs with discrepant rotation and activity qualitatively seem to have lower mass ratios. These outliers may be due to tidal spin-up, indicating these wide binaries could be hierarchical triple systems. We additionally present high resolution adaptive optics images for two wide binary systems to test this hypothesis. The demographics of stellar rotation and magnetic activity between stars in wide binaries may be useful indicators for discerning formation scenarios of these systems. △ Less

Submitted 12 December, 2017; originally announced December 2017.

Comments: 15 pages, 10 figures, accepted by ApJ

Infrared Flares from M Dwarfs: a Hinderance to Future Transiting Exoplanet Studies

Authors: James R. A Davenport

Abstract: Many current and future exoplanet missions are pushing to infrared (IR) wavelengths where the flux contrast between the planet and star is more favorable (Deming et al. 2009), and the impact of stellar magnetic activity is decreased. Indeed, a recent analysis of starspots and faculae found these forms of stellar activity do not substantially impact the transit signatures or science potential for F… ▽ More Many current and future exoplanet missions are pushing to infrared (IR) wavelengths where the flux contrast between the planet and star is more favorable (Deming et al. 2009), and the impact of stellar magnetic activity is decreased. Indeed, a recent analysis of starspots and faculae found these forms of stellar activity do not substantially impact the transit signatures or science potential for FGKM stars with JWST (Zellem et al. 2017). However, this is not true in the case of flares, which I demonstrate can be a hinderance to transit studies in this note. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: 2 pages, 1 figure. Accepted to Research Notes of the American Astronomical Society

Tidal Synchronization and Differential Rotation of Kepler Eclipsing Binaries

Authors: John C. Lurie, Karl Vyhmeister, Suzanne L. Hawley, Jamel Adilia, Andrea Chen, James R. A. Davenport, Mario Juric, Michael Puig-Holzman, Kolby L. Weisenburger

Abstract: Few observational constraints exist for the tidal synchronization rate of late-type stars, despite its fundamental role in binary evolution. We visually inspected the light curves of 2278 eclipsing binaries (EBs) from the Kepler Eclipsing Binary Catalog to identify those with starspot modulations, as well as other types of out-of-eclipse variability. We report rotation periods for 816 EBs with sta… ▽ More Few observational constraints exist for the tidal synchronization rate of late-type stars, despite its fundamental role in binary evolution. We visually inspected the light curves of 2278 eclipsing binaries (EBs) from the Kepler Eclipsing Binary Catalog to identify those with starspot modulations, as well as other types of out-of-eclipse variability. We report rotation periods for 816 EBs with starspot modulations, and find that 79% of EBs with orbital periods less than ten days are synchronized. However, a population of short period EBs exists with rotation periods typically 13% slower than synchronous, which we attribute to the differential rotation of high latitude starspots. At 10 days, there is a transition from predominantly circular, synchronized EBs to predominantly eccentric, pseudosynchronized EBs. This transition period is in good agreement with the predicted and observed circularization period for Milky Way field binaries. At orbital periods greater than about 30 days, the amount of tidal synchronization decreases. We also report 12 previously unidentified candidate $δでるた$ Scuti and $γがんま$ Doradus pulsators, as well as a candidate RS CVn system with an evolved primary that exhibits starspot occultations. For short period contact binaries, we observe a period-color relation, and compare it to previous studies. As a whole, these results represent the largest homogeneous study of tidal synchronization of late-type stars. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: Accepted for publication in the Astronomical Journal. EB rotation periods and classifications available at https://github.com/jlurie/decatur/blob/master/decatur/data/final_catalog.csv

Chromospheric Activity of HAT-P-11: an Unusually Active Planet-Hosting K Star

Authors: Brett M. Morris, Suzanne L. Hawley, Leslie Hebb, Charli Sakari, James. R. A. Davenport, Howard Isaacson, Andrew W. Howard, Benjamin T. Montet, Eric Agol

Abstract: 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… ▽ More 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. △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: 16 pages, 8 figures; accepted to the Astrophysical Journal