Showing 1–50 of 53 results for author: Pope, B

Differentiable Optics with dLux II: Optical Design Maximising Fisher Information

Authors: Louis Desdoigts, Benjamin Pope, Michael Gully-Santiago, Peter Tuthill

Abstract: The design of astronomical hardware operating at the diffraction limit requires optimization of physical optical simulations of the instrument with respect to desired figures of merit, such as throughput or astrometric accuracy. These systems can be high dimensional, with highly nonlinear relationships between outputs and the adjustable parameters of the hardware. In this series of papers we prese… ▽ More The design of astronomical hardware operating at the diffraction limit requires optimization of physical optical simulations of the instrument with respect to desired figures of merit, such as throughput or astrometric accuracy. These systems can be high dimensional, with highly nonlinear relationships between outputs and the adjustable parameters of the hardware. In this series of papers we present and apply dLux, an open-source end-to-end differentiable optical modelling framework. Automatic differentiation enables not just efficient high-dimensional optimization of astronomical hardware designs, but also Bayesian experimental design directly targeting the precision of experimental outcomes. Automatic second derivatives enable the exact and numerically stable calculation of parameter covariance forecasts, and higher derivatives of these enable direct optimization of these forecasts. We validate this method against analytic theory and illustrate its utility in evaluating the astrometric precision of a parametrized telescope model, and the design of a diffractive pupil to achieve optimal astrometric performance for exoplanet searches. The source code and tutorial software are open source and publicly available, targeting researchers who may wish to harness dLux for their own optical simulation problems. △ Less

Submitted 16 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Comments: 12 pages, 6 figures

Differentiable Optics with dLux I: Deep calibration of Flat Field and Phase Retrieval with Automatic Differentiation

Authors: Louis Desdoigts, Benjamin Pope, Jordan Dennis, Peter Tuthill

Abstract: The sensitivity limits of space telescopes are imposed by uncalibrated errors in the point spread function, photon-noise, background light, and detector sensitivity. These are typically calibrated with specialized wavefront sensor hardware and with flat fields obtained on the ground or with calibration sources, but these leave vulnerabilities to residual time-varying or non-common path aberrations… ▽ More The sensitivity limits of space telescopes are imposed by uncalibrated errors in the point spread function, photon-noise, background light, and detector sensitivity. These are typically calibrated with specialized wavefront sensor hardware and with flat fields obtained on the ground or with calibration sources, but these leave vulnerabilities to residual time-varying or non-common path aberrations and variations in the detector conditions. It is therefore desirable to infer these from science data alone, facing the prohibitively high dimensional problems of phase retrieval and pixel-level calibration. We introduce a new Python package for physical optics simulation, dLux, which uses the machine learning framework JAX to achieve GPU acceleration and automatic differentiation (autodiff), and apply this to simulating astronomical imaging. In this first of a series of papers, we show that gradient descent enabled by autodiff can be used to simultaneously perform phase retrieval and calibration of detector sensitivity, scaling efficiently to inferring millions of parameters. This new framework enables high dimensional optimization and inference in data analysis and hardware design in astronomy and beyond, which we explore in subsequent papers in this series. △ Less

Submitted 16 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Comments: 13 pages, 6 figures

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems, Vol. 9, Issue 2, 028007 (June 2023)

The James Webb Interferometer: Space-based interferometric detections of PDS 70 b and c at 4.8 $μ$m

Authors: Dori Blakely, Doug Johnstone, Gabriele Cugno, Anand Sivaramakrishnan, Peter Tuthill, Ruobing Dong, Benjamin J. S. Pope, Loïc Albert, Max Charles, Rachel A. Cooper, Matthew De Furio, Louis Desdoigts, René Doyon, Logan Francis, Alexandra Z. Greenbaum, David Lafrenière, James P. Lloyd, Michael R. Meyer, Laurent Pueyo, Shrishmoy Ray, Joel Sánchez-Bermúdez, Anthony Soulain, Deepashri Thatte, Thomas Vandal

Abstract: We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of… ▽ More We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of both PDS 70 b and c provide evidence for circumplanetary disk emission through fitting SED models to these new measurements and those found in the literature. We also newly detect emission within the disk gap at an SNR of $\sim$4, at a position angle of $207^{+11}_{-10}$ degrees, and an unconstrained separation within $\sim$200 mas. Follow-up observations will be needed to determine the nature of this emission. We place a 5$σ$ upper limit of $Δ$mag = 7.56 on the contrast of the candidate PDS 70 d at 4.8 $μ$m, which indicates that if the previously observed emission at shorter wavelengths is due to a planet, this putative planet has a different atmospheric composition than PDS 70 b or c. Finally, we place upper limits on emission from any additional planets in the disk gap. We find an azimuthally averaged 5$σ$ upper limit of $Δ$mag $\approx$ 7.5 at separations greater than 125 mas. These are the deepest limits to date within $\sim$250 mas at 4.8 $μ$m and the first space-based interferometric observations of this system. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

The Sydney Radio Star Catalogue: properties of radio stars at megahertz to gigahertz frequencies

Authors: Laura N. Driessen, Joshua Pritchard, Tara Murphy, George Heald, Jan Robrade, Barnali Das, Stefan Duchesne, David L. Kaplan, Emil Lenc, Christene R. Lynch, Benjamin J. S. Pope, Kovi Rose, Beate Stelzer, Yuanming Wang, Andrew Zic

Abstract: We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and pro… ▽ More We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor instruments; however we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We find that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the SRG spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 21 pages, 8 Figures, 7 Tables, Submitted to PASA

Coronagraphic Data Post-processing Using Projections on Instrumental Modes

Authors: Yinzi Xin, Laurent Pueyo, Romain Laugier, Leonid Pogorelyuk, Ewan S. Douglas, Benjamin J. S. Pope, Kerri L. Cahoy

Abstract: Directly observing exoplanets with coronagraphs is impeded by the presence of speckles from aberrations in the optical path, which can be mitigated in hardware with wavefront control as well as in post-processing. This work explores using an instrument model in post-processing to separate astrophysical signals from residual aberrations in coronagraphic data. The effect of wavefront error (WFE) on… ▽ More Directly observing exoplanets with coronagraphs is impeded by the presence of speckles from aberrations in the optical path, which can be mitigated in hardware with wavefront control as well as in post-processing. This work explores using an instrument model in post-processing to separate astrophysical signals from residual aberrations in coronagraphic data. The effect of wavefront error (WFE) on the coronagraphic intensity consists of a linear contribution and a quadratic contribution. When either of the terms is much larger than the other, the instrument response can be approximated by a transfer matrix mapping WFE to detector plane intensity. From this transfer matrix, a useful projection onto instrumental modes that removes the dominant error modes can be derived. We apply this projection to synthetically generated Roman Space Telescope hybrid Lyot coronagraph (HLC) data to extract "robust observables," which can be used instead of raw data for applications such as detection testing. The projection improves planet flux ratio detection limits by about 28% in the linear regime and by over a factor of 2 in the quadratic regime, illustrating that robust observables can increase sensitivity to astrophysical signals and improve the scientific yield from coronagraphic data. While this approach does not require additional information such as observations of reference stars or modulations of a deformable mirror, it can and should be combined with these other techniques, acting as a model-informed prior in an overall post-processing strategy. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 19 pages, 10 figures

Phenomenology and periodicity of radio emission from the stellar system AU Microscopii

Authors: Sanne Bloot, Joseph R. Callingham, Harish K. Vedantham, Robert D. Kavanagh, Benjamin J. S. Pope, Juan B. Climent, José Carlos Guirado, Luis Peña-Moñino, Miguel Pérez-Torres

Abstract: Stellar radio emission can measure a star's magnetic field strength and structure, plasma density and dynamics, and the stellar wind pressure impinging on exoplanet atmospheres. However, properly interpreting the radio data often requires temporal baselines that cover the rotation of the stars, orbits of their planets and any longer-term stellar activity cycles. Here we present our monitoring camp… ▽ More Stellar radio emission can measure a star's magnetic field strength and structure, plasma density and dynamics, and the stellar wind pressure impinging on exoplanet atmospheres. However, properly interpreting the radio data often requires temporal baselines that cover the rotation of the stars, orbits of their planets and any longer-term stellar activity cycles. Here we present our monitoring campaign on the young, active M dwarf AU Microscopii with the Australia Telescope Compact Array between 1.1 and 3.1 GHz. With over 250 hours of observations, these data represent the longest radio monitoring campaign on a single main-sequence star to date. We find that AU Mic produces a wide variety of radio emission, for which we introduce a phenomenological classification scheme predicated on the polarisation fraction and time-frequency structure of the emission. The radio emission detected on AU Mic can be broadly categorised into five distinct types of bursts, and broadband quiescent emission. The radio bursts are highly circularly polarised and periodic with the rotation period of the star, implying that the emission is beamed. It is therefore most likely produced by the electron cyclotron maser instability. We present a model to show that the observed emission can be explained with auroral rings on the magnetic poles. The total intensity of the broadband emission is stochastic, but we show that its circular polarisation fraction is also periodic with the rotation of the star. We present a qualitative model to describe the periodicity in the polarisation fraction of the broadband emission using gyromagnetic emission, and infer a magnetic obliquity of at least 20 degrees from the observed variation in polarisation fraction. Finally, we show that the radio emission might be evolving on long timescales, hinting at a potential stellar magnetic activity cycle. △ Less

Submitted 19 December, 2023; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A, typos corrected

Astrometry and Precise Radial Velocities Yield a Complete Orbital Solution for the Nearby Eccentric Brown Dwarf LHS 1610 b

Authors: Evan Fitzmaurice, Gudmundur Stefánsson, Robert D. Kavanagh, Suvrath Mahadevan, Caleb I. Cañas, Joshua N. Winn, Paul Robertson, Joe P. Ninan, Simon Albrecht, J. R. Callingham, William D. Cochran, Megan Delamer, Shubham Kanodia, Andrea S. J. Lin, Marcus L. Marcussen, Benjamin J. S. Pope, Lawrence W. Ramsey, Arpita Roy, Harish Vedantham, Jason T. Wright

Abstract: We characterize the LHS 1610 system, a nearby ($d=9.7$ pc) M5 dwarf hosting a brown dwarf in a $10.6$ day, eccentric ($e \sim 0.37$) orbit. A joint fit of the available Gaia two-body solution, discovery radial velocities (RVs) from TRES, and new RVs obtained with the Habitable-zone Planet Finder, yields an orbital inclination of $117.2\pm0.9^\circ$ and a mass constraint of $50.9\pm0.9$ M$_J$. This… ▽ More We characterize the LHS 1610 system, a nearby ($d=9.7$ pc) M5 dwarf hosting a brown dwarf in a $10.6$ day, eccentric ($e \sim 0.37$) orbit. A joint fit of the available Gaia two-body solution, discovery radial velocities (RVs) from TRES, and new RVs obtained with the Habitable-zone Planet Finder, yields an orbital inclination of $117.2\pm0.9^\circ$ and a mass constraint of $50.9\pm0.9$ M$_J$. This gives LHS 1610 b the second most precise mass of brown dwarfs orbiting M stars within 25pc. We highlight a discrepancy between the Gaia two-body solution eccentricity ($e=0.52 \pm 0.03$) and that from the RVs ($e=0.3702\pm0.0003$), which requires the astrometric time-series release (Gaia DR4) for further diagnostics. With a flare rate of $0.28\pm 0.07$ flares/day from TESS photometry, and a rotation period of $84 \pm 8$ days, LHS 1610 joins other mid M stars -- including Proxima Centauri and YZ Ceti -- as nearby mid M dwarfs with flare rates on the higher end for their long rotation periods. These stars are promising candidates for searching for sub-Alfvénic star-companion interactions, raising the question whether LHS 1610 b could be driving the flares on its host star. However, the available TESS photometry is insufficient to confirm or rule out any orbital phase-dependence of the flares. We show that the LHS 1610 system, as a nearby mid M star with a large, short-period companion, is a promising target to look for evidence of star-companion interactions or aural emission from the brown dwarf at radio wavelengths. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 24 pages, 7 figures, 3 tables. Submitted to AAS Journals on Oct 11, 2023

JWST/NIRCam discovery of the first Y+Y brown dwarf binary: WISE J033605.05$-$014350.4

Authors: Per Calissendorff, Matthew De Furio, Michael Meyer, Loïc Albert, Christian Aganze, Mohamad Ali-Dib, Daniella C. Bardalez Gagliuffi, Frederique Baron, Charles A. Beichman, Adam J. Burgasser, Michael C. Cushing, Jacqueline Kelly Faherty, Clémence Fontanive, Christopher R. Gelino, John E. Gizis, Alexandra Z. Greenbaum, J. Davy Kirkpatrick, Sandy K. Leggett, Frantz Martinache, David Mary, Mamadou N'Diaye, Benjamin J. S. Pope, Thomas L Roellig, Johannes Sahlmann, Anand Sivaramakrishnan , et al. (3 additional authors not shown)

Abstract: We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in… ▽ More We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in F150W and F480M, respectively. At a distance of 10$\,$pc based on its Spitzer parallax, and assuming a random inclination distribution, the physical separation is approximately 1$\,$au. Evolutionary models predict for that an age of 1-5 Gyr, the companion mass is about 4-12.5 Jupiter masses around the 7.5-20 Jupiter mass primary, corresponding to a companion-to-host mass fraction of $q=0.61\pm0.05$. Under the assumption of a Keplerian orbit the period for this extreme binary is in the range of 5-9 years. The system joins a small but growing sample of ultracool dwarf binaries with effective temperatures of a few hundreds of Kelvin. Brown dwarf binaries lie at the nexus of importance for understanding the formation mechanisms of these elusive objects, as they allow us to investigate whether the companions formed as stars or as planets in a disk around the primary. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 8 pages, 3 figures, 1 table. Accepted for publication in Astrophysical Journal Letters

V-LoTSS: The Circularly-Polarised LOFAR Two-metre Sky Survey

Authors: J. R. Callingham, T. W. Shimwell, H. K. Vedantham, C. G. Bassa, S. P. O'Sullivan, T. W. H. Yiu, S. Bloot, P. N. Best, M. J. Hardcastle, M. Haverkorn, R. D. Kavanagh, L. Lamy, B. J. S. Pope, H. J. A. Röttgering, D. J. Schwarz, C. Tasse, R. J. van Weeren, G. J. White, P. Zarka, D. J. Bomans, A. Bonafede, M. Bonato, A. Botteon, M. Bruggen, K. T. Chyży , et al. (22 additional authors not shown)

Abstract: We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $μ$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polar… ▽ More We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $μ$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polarisation is measured to be $\approx$0.06%, and our survey is complete at flux densities $\geq1$ mJy. A detection is considered reliable when the circularly-polarised fraction exceeds 1%. We find the population of circularly-polarised sources is composed of four distinct classes: stellar systems, pulsars, active galactic nuclei, and sources unidentified in the literature. The stellar systems can be further separated into chromospherically-active stars, M dwarfs, and brown dwarfs. Based on the circularly-polarised fraction and lack of an optical counterpart, we show it is possible to infer whether the unidentified sources are likely unknown pulsars or brown dwarfs. By the completion of this survey of the northern sky, we expect to detect 300$\pm$100 circularly-polarised sources. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A, 15 pages, 8 figures, 1 table. The catalogue will be publicly available at http://lofar-surveys.org/ and via Vizier shortly

Journal ref: A&A 670, A124 (2023)

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- IV. Aperture Masking Interferometry

Authors: Anand Sivaramakrishnan, Peter Tuthill, James P. Lloyd, Alexandra Z. Greenbaum, Deepashri Thatte, Rachel A. Cooper, Thomas Vandal, Jens Kammerer, Joel Sanchez-Bermudez, Benjamin J. S. Pope, Dori Blakely, Loïc Albert, Neil J. Cook, Doug Johnstone, André R. Martel, Kevin Volk, Anthony Soulain, Étienne Artigau, David Lafrenière, Chris J. Willott, Sébastien Parmentier, K. E. Saavik Ford, Barry McKernan, M. Begoña Vila, Neil Rowlands , et al. (14 additional authors not shown)

Abstract: The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, pres… ▽ More The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, present some sample science cases, explain its operational observing strategies, indicate how AMI proposals can be developed with data simulations, and how AMI data can be analyzed. We also present key results from commissioning AMI. Since the allied Kernel Phase Imaging (KPI) technique benefits from AMI operational strategies, we also cover NIRISS KPI methods and analysis techniques, including a new user-friendly KPI pipeline. The NIRISS KPI bright limit is $\simeq 8$ W2 magnitudes. AMI (and KPI) achieve an inner working angle of $\sim 70$ mas that is well inside the $\sim 400$ mas NIRCam inner working angle for its circular occulter coronagraphs at comparable wavelengths. △ Less

Submitted 7 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 30 pages, 10 figures

Modelling Cosmic Radiation Events in the Tree-ring Radiocarbon Record

Authors: Qingyuan Zhang, Utkarsh Sharma, Jordan A. Dennis, Andrea Scifo, Margot Kuitems, Ulf Buentgen, Mathew J. Owens, Michael W. Dee, Benjamin J. S. Pope

Abstract: Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked… ▽ More Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked in the annual growth rings of trees. To interpret high-resolution tree-ring radiocarbon measurements therefore necessitates modelling the entire global carbon cycle. Here, we introduce 'ticktack', the first open-source Python package that connects box models of the carbon cycle with modern Bayesian inference tools. We use this to analyse all public annual 14C tree data, and infer posterior parameters for all six known Miyake events. They do not show a consistent relationship to the solar cycle, and several display extended durations that challenge either astrophysical or geophysical models. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: Accepted Proceedings of the Royal Society A. 19 pages 6 figures body, 12 pages appendices which are supplementary material in the published version

Solar Energetic Particle Ground-Level Enhancements and the Solar Cycle

Authors: Mathew Owens, Luke Barnard, Benjamin Pope, Mike Lockwood, Ilya Usoskin, Eleanna Asvestari

Abstract: Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C,… ▽ More Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C, $^{10}$Be and $^{36}$Cl) records that are attributed to significantly larger solar energetic particle (SEP) events than have been observed during the space age. Despite both storms and SEPs being driven by solar eruptive phenomena, the event-by-event correspondence between extreme storms and extreme SEPs is low. Thus it should not be assumed a priori that the solar cycle patterns found for storms also hold for SEPs and the cosmogenic-isotope events. In this study we investigate the solar cycle trends in the timing and magnitude of the 67 SEP ground-level enhancements (GLEs) recorded by neutron monitors since the mid 1950s. Using a number of models of GLE occurrence probability, we show that GLEs are around a factor four more likely around solar maximum than around solar minimum, and that they preferentially occur earlier in even-numbered solar cycles than in odd-numbered cycles. There are insufficient data to conclusively determine whether larger solar cycles produce more GLEs. Implications for putative space-weather events in the cosmogenic-isotope records are discussed. We find that GLEs tend to cluster within a few tens of days, likely due to particularly productive individual active regions, and with approximately 11-year separations, owing to the solar cycle ordering. But these timescales do not explain cosmogenic-isotope spikes which require multiple extreme SEP events over consecutive years. △ Less

Submitted 27 July, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 22 pages, 9 figures

The TESS View of LOFAR Radio-Emitting Stars

Authors: Benjamin J. S. Pope, Joseph R. Callingham, Adina D. Feinstein, Maximilian N. Günther, Harish K. Vedantham, Megan Ansdell, Timothy W. Shimwell

Abstract: The recent detection of the M dwarf GJ 1151 at 144 MHz low radio frequencies using LOFAR has been interpreted as evidence of an exoplanet magnetically interacting with its host star. This would be the first exoplanet detected around a main sequence star by a radio telescope. Radial velocity confirmation of such a planet has proven inconclusive, and it remains possible that the radio emission could… ▽ More The recent detection of the M dwarf GJ 1151 at 144 MHz low radio frequencies using LOFAR has been interpreted as evidence of an exoplanet magnetically interacting with its host star. This would be the first exoplanet detected around a main sequence star by a radio telescope. Radial velocity confirmation of such a planet has proven inconclusive, and it remains possible that the radio emission could be generated by a stellar coronal process. Using data from TESS, we shed light on this question by probing the stellar activity and flares of GJ 1151 as well as fourteen other M dwarfs detected by LOFAR. GJ 1151 and three other star-planet interaction candidates are found to be inactive, with no rotational modulation and few, if any, flares. The remainder of the LOFAR detected M dwarfs flare frequently. We consider it unlikely that stellar activity is responsible for the bright, circularly-polarized emission from GJ 1151 and its three analogs, and support the star-planet magnetic interaction interpretation. △ Less

Submitted 14 October, 2021; v1 submitted 10 October, 2021; originally announced October 2021.

Comments: ApJL: https://iopscience.iop.org/article/10.3847/2041-8213/ac230c

The population of M dwarfs observed at low radio frequencies

Authors: J. R. Callingham, H. K. Vedantham, T. W. Shimwell, B. J. S. Pope, I. E. Davis, P. N. Best, M. J. Hardcastle, H. J. A. Rottgering, J. Sabater, C. Tasse, R. J. van Weeren, W. L. Williams, P. Zarka, F. de Gasperin, A. Drabent

Abstract: Coherent low-frequency ($\lesssim 200$ MHz) radio emission from stars encodes the conditions of the outer corona, mass-ejection events, and space weather. Previous low-frequency searches for radio emitting stellar systems have lacked the sensitivity to detect the general population, instead largely focusing on targeted studies of anomalously active stars. Here we present 19 detections of coherent… ▽ More Coherent low-frequency ($\lesssim 200$ MHz) radio emission from stars encodes the conditions of the outer corona, mass-ejection events, and space weather. Previous low-frequency searches for radio emitting stellar systems have lacked the sensitivity to detect the general population, instead largely focusing on targeted studies of anomalously active stars. Here we present 19 detections of coherent radio emission associated with known M~dwarfs from a blind flux-limited low-frequency survey. Our detections show that coherent radio emission is ubiquitous across the M~dwarf main sequence, and that the radio luminosity is independent of known coronal and chromospheric activity indicators. While plasma emission can generate the low-frequency emission from the most chromospherically active stars of our sample, the origin of the radio emission from the most quiescent sources is yet to be ascertained. Large-scale analogues of the magnetospheric processes seen in gas-giant planets likely drive the radio emission associated with these quiescent stars. The slowest-rotating stars of this sample are candidate systems to search for star-planet interaction signatures. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: Accepted to Nature Astronomy; 31 pages, 9 figures, and 2 tables

A search for transits among the δ Scuti variables in Kepler

Authors: Daniel R. Hey, Benjamin T. Montet, Benjamin J. S. Pope, Simon J. Murphy, Timothy R. Bedding

Abstract: We search for transits around all known pulsating δ Sct variables (6500 K < Teff < 10 000 K) in the long-cadence Kepler data after subtracting the pulsation signal through an automated routine. To achieve this, we devise a simple and computationally inexpensive method for distinguishing between low-frequency pulsations and transits in light curves. We find 3 new candidate transit events that were… ▽ More We search for transits around all known pulsating δ Sct variables (6500 K < Teff < 10 000 K) in the long-cadence Kepler data after subtracting the pulsation signal through an automated routine. To achieve this, we devise a simple and computationally inexpensive method for distinguishing between low-frequency pulsations and transits in light curves. We find 3 new candidate transit events that were previously hidden behind the pulsations, but caution that they are likely to be false positive events. We also examined the Kepler Objects of Interest catalog and identify 13 additional host stars which show δ Sct pulsations. For each star in our sample, we use the non-detection of pulsation timing variations for a planet that is known to be transiting a δ Sct variable to obtain both an upper limit on the mass of the planet and the expected radial velocity semi-amplitude of the host star. Simple injection tests of our pipeline imply 100% recovery for planets of 0.5 RJup or greater. Extrapolating our number of Kepler δ Sct stars, we expect 12 detectable planets above 0.5 RJup in TESS. Our sample contains some of the hottest known transiting planets around evolved stars, and is the first complete sample of transits around δ Sct variables. We make available our code and pulsation-subtracted light curves to facilitate further analysis. △ Less

Submitted 8 August, 2021; originally announced August 2021.

Comments: 18 pages, accepted to AAS journals

Asteroseismology of iota Draconis and Discovery of an Additional Long-Period Companion

Authors: Michelle L. Hill, Stephen R. Kane, Tiago L. Campante, Zhexing Li, Paul A. Dalba, Timothy D. Brandt, Timothy R. White, Benjamin J. S. Pope, Keivan G. Stassun, Benjamin J. Fulton, Enrico Corsaro, Tanda Li, J. M. Joel Ong, Timothy R. Bedding, Diego Bossini, Derek L. Buzasi, William J. Chaplin, Margarida S. Cunha, Rafael A. Garcia, Sylvain N. Breton, Marc Hon, Daniel Huber, Chen Jiang, Cenk Kayhan, James S. Kuszlewicz , et al. (3 additional authors not shown)

Abstract: Giant stars as known exoplanet hosts are relatively rare due to the potential challenges in acquiring precision radial velocities and the small predicted transit depths. However, these giant host stars are also some of the brightest in the sky and so enable high signal-to-noise follow-up measurements. Here we report on new observations of the bright (V ~ 3.3) giant star $ι$ Draconis ($ι$ Dra), kno… ▽ More Giant stars as known exoplanet hosts are relatively rare due to the potential challenges in acquiring precision radial velocities and the small predicted transit depths. However, these giant host stars are also some of the brightest in the sky and so enable high signal-to-noise follow-up measurements. Here we report on new observations of the bright (V ~ 3.3) giant star $ι$ Draconis ($ι$ Dra), known to host a planet in a highly eccentric ~511 day period orbit. TESS observations of the star over 137 days reveal asteroseismic signatures, allowing us to constrain the stellar radius, mass, and age to ~2%, ~6%, and ~28%, respectively. We present the results of continued radial velocity monitoring of the star using the Automated Planet Finder over several orbits of the planet. We provide more precise planet parameters of the known planet and, through the combination of our radial velocity measurements with Hipparcos and Gaia astrometry, we discover an additional long-period companion with an orbital period of ~$68^{+60}_{-36}$ years. Mass predictions from our analysis place this sub-stellar companion on the border of the planet and brown dwarf regimes. The bright nature of the star combined with the revised orbital architecture of the system provides an opportunity to study planetary orbital dynamics that evolve as the star moves into the giant phase of its evolution. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 20 pages, 11 Figures, Accepted for publication in the Astronomical Journal

Phase Retrieval and Design with Automatic Differentiation

Authors: Alison Wong, Benjamin Pope, Louis Desdoigts, Peter Tuthill, Barnaby Norris, Chris Betters

Abstract: The principal limitation in many areas of astronomy, especially for directly imaging exoplanets, arises from instability in the point spread function (PSF) delivered by the telescope and instrument. To understand the transfer function, it is often necessary to infer a set of optical aberrations given only the intensity distribution on the sensor - the problem of phase retrieval. This can be import… ▽ More The principal limitation in many areas of astronomy, especially for directly imaging exoplanets, arises from instability in the point spread function (PSF) delivered by the telescope and instrument. To understand the transfer function, it is often necessary to infer a set of optical aberrations given only the intensity distribution on the sensor - the problem of phase retrieval. This can be important for post-processing of existing data, or for the design of optical phase masks to engineer PSFs optimized to achieve high contrast, angular resolution, or astrometric stability. By exploiting newly efficient and flexible technology for automatic differentiation, which in recent years has undergone rapid development driven by machine learning, we can perform both phase retrieval and design in a way that is systematic, user-friendly, fast, and effective. By using modern gradient descent techniques, this converges efficiently and is easily extended to incorporate constraints and regularization. We illustrate the wide-ranging potential for this approach using our new package, Morphine. Challenging applications performed with this code include precise phase retrieval for both discrete and continuous phase distributions, even where information has been censored such as heavily-saturated sensor data. We also show that the same algorithms can optimize continuous or binary phase masks that are competitive with existing best solutions for two example problems: an Apodizing Phase Plate (APP) coronagraph for exoplanet direct imaging, and a diffractive pupil for narrow-angle astrometry. The Morphine source code and examples are available open-source, with a similar interface to the popular physical optics package Poppy. △ Less

Submitted 2 July, 2021; originally announced July 2021.

TESS Data for Asteroseismology: Photometry

Authors: Rasmus Handberg, Mikkel N. Lund, Timothy R. White, Oliver J. Hall, Derek L. Buzasi, Benjamin J. S. Pope, Jonas S. Hansen, Carolina von Essen, Lindsey Carboneau, Daniel Huber, Roland K. Vanderspek, Michael M. Fausnaug, Peter Tenenbaum, Jon M. Jenkins, the T'DA Collaboration

Abstract: Over the last two decades, asteroseismology has increasingly proven to be the observational tool of choice for the study of stellar physics, aided by the high quality of data available from space-based missions such as CoRoT, Kepler, K2 and TESS. TESS in particular will produce more than an order of magnitude more such data than has ever been available before. While the standard TESS mission pro… ▽ More Over the last two decades, asteroseismology has increasingly proven to be the observational tool of choice for the study of stellar physics, aided by the high quality of data available from space-based missions such as CoRoT, Kepler, K2 and TESS. TESS in particular will produce more than an order of magnitude more such data than has ever been available before. While the standard TESS mission products include light curves from 120-sec observations suitable for both exoplanet and asteroseismic studies, they do not include light curves for the vastly larger number of targets observed by the mission at a longer 1800-sec cadence in Full Frame Images (FFIs). To address this lack, the TESS Data for Asteroseismology (T'DA) group under the TESS Asteroseismic Science Consortium (TASC), has constructed an open-source pipeline focused on producing light curves for all stars observed by TESS at all cadences, currently including stars down to a TESS magnitude of 15. The pipeline includes target identification, background estimation and removal, correction of FFI timestamps, and a range of potential photometric extraction methodologies, though aperture photometry is currently the default approach. For the brightest targets, we transparently apply a halo photometry algorithm to construct a calibrated light curve from unsaturated pixels in the image. In this paper, we describe in detail the algorithms, functionality, and products of this pipeline, and summarize the noise metrics for the light curves. Companion papers will address the removal of systematic noise sources from our light curves, and a stellar variability classification from these. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 27 pages, 14 figures, Accepted for publication in The Astronomical Journal

Low-frequency monitoring of flare star binary CR Draconis: Long-term electron-cyclotron maser emission

Authors: J. R. Callingham, B. J. S. Pope, A. D. Feinstein, H. K. Vedantham, T. W. Shimwell, P. Zarka, C. Tasse, L. Lamy, K. Veken, S. Toet, J. Sabater, P. N. Best, R. J. van Weeren, H. J. A. Röttgering, T. P. Ray

Abstract: Recently detected coherent low-frequency radio emission from M dwarf systems shares phenomenological similarities with emission produced by magnetospheric processes from the gas giant planets of our Solar System. Such beamed electron-cyclotron maser emission can be driven by a star-planet interaction or a breakdown in co-rotation between a rotating plasma disk and a stellar magnetosphere. Both mod… ▽ More Recently detected coherent low-frequency radio emission from M dwarf systems shares phenomenological similarities with emission produced by magnetospheric processes from the gas giant planets of our Solar System. Such beamed electron-cyclotron maser emission can be driven by a star-planet interaction or a breakdown in co-rotation between a rotating plasma disk and a stellar magnetosphere. Both models suggest that the radio emission could be periodic. Here we present the longest low-frequency interferometric monitoring campaign of an M dwarf system, composed of twenty-one $\approx$8 hour epochs taken in two series of observing blocks separated by a year. We achieved a total on-source time of 6.5 days. We show that the M dwarf binary CR Draconis has a low-frequency 3$σ$ detection rate of 90$^{+5}_{-8}$% when a noise floor of $\approx$0.1 mJy is reached, with a median flux density of 0.92 mJy, consistent circularly polarised handedness, and a median circularly polarised fraction of 66%. We resolve three bright radio bursts in dynamic spectra, revealing the brightest is elliptically polarised, confined to 4 MHz of bandwidth centred on 170 MHz, and reaches a flux density of 205 mJy. The burst structure is mottled, indicating it consists of unresolved sub-bursts. Such a structure shares a striking resemblance with the low-frequency emission from Jupiter. We suggest the near-constant detection of high brightness temperature, highly-circularly-polarised radiation that has a consistent circular polarisation handedness implies the emission is produced via the electron-cyclotron maser instability. Optical photometric data reveal the system has a rotation period of 1.984$\pm$0.003 days. We observe no periodicity in the radio data, but the sampling of our radio observations produces a window function that would hide the near two-day signal. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: Accepted for publication in A&A, 16 pages, 7 figures, 2 tables

Journal ref: A&A 648, A13 (2021)

Kernel Phase and Coronagraphy with Automatic Differentiation

Authors: Benjamin J. S. Pope, Laurent Pueyo, Yinzi Xin, Peter G. Tuthill

Abstract: The accumulation of aberrations along the optical path in a telescope produces distortions and speckles in the resulting images, limiting the performance of cameras at high angular resolution. It is important to achieve the highest possible sensitivity to faint sources such as planets, using both hardware and data analysis software. While analytic methods are efficient, real systems are better-mod… ▽ More The accumulation of aberrations along the optical path in a telescope produces distortions and speckles in the resulting images, limiting the performance of cameras at high angular resolution. It is important to achieve the highest possible sensitivity to faint sources such as planets, using both hardware and data analysis software. While analytic methods are efficient, real systems are better-modelled numerically, but such models with many parameters can be hard to understand, optimize and apply. Automatic differentiation software developed for machine learning now makes calculating derivatives with respect to aberrations straightforward for arbitrary optical systems. We apply this powerful new tool to enhance high-angular-resolution astronomical imaging. Self-calibrating observables such as the 'closure phase' or 'bispectrum' have been widely used in optical and radio astronomy to mitigate optical aberrations and achieve high-fidelity imagery. Kernel phases are a generalization of closure phases in the limit of small phase errors. Using automatic differentiation, we reproduce existing kernel phase theory within this framework and demonstrate an extension to the Lyot coronagraph, finding self-calibrating combinations of speckles which are resistant to phase noise, but only in the very high-wavefront-quality regime. As an illustrative example, we reanalyze Palomar adaptive optics observations of the binary alpha Ophiuchi, finding consistency between the new pipeline and the existing standard. We present a new Python package 'morphine' that incorporates these ideas, with an interface similar to the popular package poppy, for optical simulation with automatic differentiation. These methods may be useful for designing improved astronomical optical systems by gradient descent. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: Accepted ApJ

The extreme colliding-wind system Apep: resolved imagery of the central binary and dust plume in the infrared

Authors: Y. Han, P. G. Tuthill, R. M. Lau, A. Soulain, J. R. Callingham, P. M. Williams, P. A. Crowther, B. J. S. Pope, B. Marcote

Abstract: The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark… ▽ More The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf-Rayet binary: a system that has so far remained spatially unresolved. Here we present an updated proper motion study deriving the expansion speed of Apep's dust plume over a two-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high-angular-resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf-Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf-Rayet binary and lending further support to the anisotropic wind model. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society. 17 pages, 9 figures, 6 tables

EarthFinder Probe Mission Concept Study: Characterizing nearby stellar exoplanet systems with Earth-mass analogs for future direct imaging

Authors: Peter Plavchan, Gautam Vasisht, Chas Beichman, Heather Cegla, Xavier Dumusque, Sharon Wang, Peter Gao, Courtney Dressing, Fabienne Bastien, Sarbani Basu, Thomas Beatty, Andrew Bechter, Eric Bechter, Cullen Blake, Vincent Bourrier, Bryson Cale, David Ciardi, Jonathan Crass, Justin Crepp, Katherine de Kleer, Scott Diddams, Jason Eastman, Debra Fischer, Jonathan Gagné, Scott Gaudi , et al. (30 additional authors not shown)

Abstract: EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to th… ▽ More EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Comments: NASA Probe Mission concept white paper for 2020 Astrophysics National Academies Decadal Survey

Two Wolf-Rayet stars at the heart of colliding-wind binary Apep

Authors: J. R. Callingham, P. A. Crowther, P. M. Williams, P. G. Tuthill, Y. Han, B. J. S. Pope, B. Marcote

Abstract: Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the diff… ▽ More Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS, 10 pages, 9 figures, 1 table

Forward modeling the orbits of companions to pulsating stars from their light travel time variations

Authors: Daniel R. Hey, Simon J. Murphy, Daniel Foreman-Mackey, Timothy R. Bedding, Benjamin J. S. Pope, David W. Hogg

Abstract: Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve… ▽ More Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to under-sampling near periapsis. We introduce a new approach that simultaneously forward-models each sample in the light curve and show that this method improves upon current sensitivity limits - especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the SNR of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Comments: 14 pages, accepted for publication in AJ

Coherent radio emission from a quiescent red dwarf indicative of star-planet interaction

Authors: H. K. Vedantham, J. R. Callingham, T. W. Shimwell, C. Tasse, B. J. S. Pope, M. Bedell, I. Snellen, P. Best, M. J. Hardcastle, M. Haverkorn, A. Mechev, S. P. O'Sullivan, H. J. A. Röttgering, G. J. White

Abstract: Low frequency ($ν\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active f… ▽ More Low frequency ($ν\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfvénic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: Accepted version, Nat Astron (2020)

A Mystery in Chamaeleon: Serendipitous Discovery of a Galactic Symbiotic Nova

Authors: Lachlan Lancaster, Jenny Greene, Yuan-Sen Ting, Sergey E. Koposov, Benjamin J. S. Pope, Rachael L. Beaton

Abstract: We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past… ▽ More We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past, preceding an increase in optical brightness by $\sim\! 8\,$magnitudes and a persistence at this luminosity for about 3 years, followed by a period of $\sim\! 1.4\,{\rm mag}\,{\rm yr}^{-1}$ dimming. The object's current optical luminosity seems to be dominated by H$α$ emission, which also exhibits blue-shifted absorption (a P-Cygni-like profile). After consideration of a number of theories to explain these myriad observations, we determine that V$^*$ CN Cha is most likely a symbiotic (an evolved star-white dwarf binary) system which has undergone a long-duration, low luminosity, nova. Interpreted in this way, the outburst in V$^*$ CN Cha is among the lowest luminosity novae ever observed. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: Submitted to AJ, 19 pages, 9 figure

No Massive Companion to the Coherent Radio-Emitting M Dwarf GJ 1151

Authors: Benjamin J. S. Pope, Megan Bedell, Joseph R. Callingham, Harish K. Vedantham, Ignas A. G. Snellen, Adrian M. Price-Whelan, Timothy W. Shimwell

Abstract: The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the… ▽ More The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the mass of any close companion to GJ 1151 at Msini < 5.6 M earth. With no stellar, brown dwarf, or giant planet companion likely in a close orbit, our data are consistent with detected radio emission emerging from a magnetic interaction between a short-period terrestrial-mass planet and GJ 1151. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Journal ref: Published in ApJL: https://iopscience.iop.org/article/10.3847/2041-8213/ab5b99

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Authors: Daniel Lecoanet, Matteo Cantiello, Eliot Quataert, Louis-Alexandre Couston, Keaton J. Burns, Benjamin J. S. Pope, Adam S. Jermyn, Benjamin Favier, Michael Le Bars

Abstract: Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, det… ▽ More Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars. △ Less

Submitted 12 November, 2019; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: Accepted to ApJL. Supplementary information available at: http://www.princeton.edu/~lecoanet/data/SupplementaryInfo_Lecoanet2019.pdf

Extended Aperture Photometry of K2 RR Lyrae stars

Authors: Emese Plachy, László Molnár, Attila Bódi, Marek Skarka, Pál Szabó, Róbert Szabó, Péter Klagyivik, Ádám Sódor, Benjamin J. S. Pope

Abstract: The \textit{Kepler} space telescope observed thousands of RR Lyrae stars in the K2 mission. In this paper we present our photometric solutions using extended apertures in order to conserve the flux of the stars to the highest possible extent. With this method we are able to avoid most of the problems that RR Lyrae light curves produced by other pipelines suffer from. For post-processing we apply t… ▽ More The \textit{Kepler} space telescope observed thousands of RR Lyrae stars in the K2 mission. In this paper we present our photometric solutions using extended apertures in order to conserve the flux of the stars to the highest possible extent. With this method we are able to avoid most of the problems that RR Lyrae light curves produced by other pipelines suffer from. For post-processing we apply the K2SC pipeline to our light curves. We provide the EAP (Extended Aperture Photometry) of 432 RR Lyrae stars observed in campaigns 3, 4, 5, and 6. We also provide subclass classifications based on Fourier parameters. We investigated in particular the presence of the Blazhko effect in the stars, and found it to be 44.7\% among the RRab stars, in agreement with results from independent samples. We found that the amplitude and phase modulation in the Blazhko stars may behave rather differently, at least over the length of a K2 Campaign. We also identified four anomalous Cepheid candidates in the sample one of which is potentially the first Blazhko-modulated member of its class. △ Less

Submitted 1 September, 2019; originally announced September 2019.

Comments: 20 pages, 13 figures, accepted in ApJS, light curves available at https://konkoly.hu/KIK/data_en.html

The K2 Bright Star Survey I: Methodology and Data Release

Authors: Benjamin J. S. Pope, Timothy R. White, Will M. Farr, Jie Yu, Michael Greklek-McKeon, Daniel Huber, Conny Aerts, Suzanne Aigrain, Timothy R. Bedding, Tabetha Boyajian, Orlagh L. Creevey, David W. Hogg

Abstract: While the Kepler Mission was designed to look at tens of thousands of faint stars (V > 12), brighter stars that saturated the detector are important because they can be and have been observed very accurately by other instruments. By analyzing the unsaturated scattered-light `halo' around these stars, we have retrieved precise light curves of most of the brightest stars in K2 fields from Campaign~4… ▽ More While the Kepler Mission was designed to look at tens of thousands of faint stars (V > 12), brighter stars that saturated the detector are important because they can be and have been observed very accurately by other instruments. By analyzing the unsaturated scattered-light `halo' around these stars, we have retrieved precise light curves of most of the brightest stars in K2 fields from Campaign~4 onwards. The halo method does not depend on the detailed cause and form of systematics, and we show that it is effective at extracting light curves from both normal and saturated stars. The key methodology is to optimize the weights of a linear combination of pixel time series with respect to an objective function. We test a range of such objective functions, finding that lagged Total Variation, a generalization of Total Variation, performs well on both saturated and unsaturated K2 targets. Applying this to the bright stars across the K2 Campaigns reveals stellar variability ubiquitously, including effects of stellar pulsation, rotation, and binarity. We describe our pipeline and present a catalogue of the 161 bright stars, with classifications of their variability, asteroseismic parameters for red giants with well-measured solar-like oscillations, and remarks on interesting objects. These light curves are publicly available as a High Level Science Product from the Mikulski Archive for Space Telescopes (MAST). △ Less

Submitted 26 August, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: Accepted ApJS. Paper documents data release of K2 Halo light curves, available from http://archive.stsci.edu/doi/resolve/resolve.html?doi=10.17909/t9-6wj4-eb32 or https://github.com/benjaminpope/k2halo. v2: corrected typo

The Kepler Smear Campaign: Light curves for 102 Very Bright Stars

Authors: Benjamin J. S. Pope, Guy R. Davies, Keith Hawkins, Timothy R. White, Amalie Stokholm, Allyson Bieryla, David W. Latham, Madeline Lucey, Conny Aerts, Suzanne Aigrain, Victoria Antoci, Timothy R. Bedding, Dominic M. Bowman, Douglas A. Caldwell, Ashley Chontos, Gilbert A. Esquerdo, Daniel Huber, Paula Jofre, Simon J. Murphy, Timothy van Reeth, Victor Silva Aguirre, Jie Yu

Abstract: We present the first data release of the Kepler Smear Campaign, using collateral 'smear' data obtained in the Kepler four-year mission to reconstruct light curves of 102 stars too bright to have been otherwise targeted. We describe the pipeline developed to extract and calibrate these light curves, and show that we attain photometric precision comparable to stars analyzed by the standard pipeline… ▽ More We present the first data release of the Kepler Smear Campaign, using collateral 'smear' data obtained in the Kepler four-year mission to reconstruct light curves of 102 stars too bright to have been otherwise targeted. We describe the pipeline developed to extract and calibrate these light curves, and show that we attain photometric precision comparable to stars analyzed by the standard pipeline in the nominal Kepler mission. In this paper, aside from publishing the light curves of these stars, we focus on 66 red giants for which we detect solar-like oscillations, characterizing 33 of these in detail with spectroscopic chemical abundances and asteroseismic masses as benchmark stars. We also classify the whole sample, finding nearly all to be variable, with classical pulsations and binary effects. All source code, light curves, TRES spectra, and asteroseismic and stellar parameters are publicly available as a Kepler legacy sample. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 35 pages, accepted ApJS

Low-frequency gravity waves in blue supergiants revealed by high-precision space photometry

Authors: Dominic M. Bowman, Siemen Burssens, May G. Pedersen, Cole Johnston, Conny Aerts, Bram Buysschaert, Mathias Michielsen, Andrew Tkachenko, Tamara M. Rogers, Philipp V. F. Edelmann, Rathish P. Ratnasingam, Sergio Simón-Díaz, Norberto Castro, Ehsan Moravveji, Benjamin J. S. Pope, Timothy R. White, Peter De Cat

Abstract: Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pe… ▽ More Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: Publication date: 6 May 2019 at 16:00h (GMT)

Journal ref: Nature Astronomy, 2019

Asteroseismology of the Hyades red giant and planet host epsilon Tauri

Authors: Torben Arentoft, Frank Grundahl, Timothy R. White, Ditte Slumstrup, Rasmus Handberg, Mikkel N. Lund, Karsten Brogaard, Mads F. Andersen, Victor Silva Aguirre, Chunguang Zhang, Xiaodian Chen, Zhengzhou Yan, Benjamin J. S. Pope, Daniel Huber, Hans Kjeldsen, Jørgen Christensen-Dalsgaard, Jens Jessen-Hansen, Victoria Antoci, Søren Frandsen, Timothy R. Bedding, Pere L. Palle, Rafael A. Garcia, Licai Deng, Marc Hon, Dennis Stello , et al. (1 additional authors not shown)

Abstract: Asteroseismic analysis of solar-like stars allows us to determine physical parameters such as stellar mass, with a higher precision compared to most other methods. Even in a well-studied cluster such as the Hyades, the masses of the red giant stars are not well known, and previous mass estimates are based on model calculations (isochrones). The four known red giants in the Hyades are assumed to be… ▽ More Asteroseismic analysis of solar-like stars allows us to determine physical parameters such as stellar mass, with a higher precision compared to most other methods. Even in a well-studied cluster such as the Hyades, the masses of the red giant stars are not well known, and previous mass estimates are based on model calculations (isochrones). The four known red giants in the Hyades are assumed to be clump (core-helium-burning) stars based on their positions in colour-magnitude diagrams, however asteroseismology offers an opportunity to test this assumption. Using asteroseismic techniques combined with other methods, we aim to derive physical parameters and the evolutionary stage for the planet hosting star epsilon Tau, which is one of the four red giants located in the Hyades. We analysed time-series data from both ground and space to perform the asteroseismic analysis. By combining high signal-to-noise (S/N) radial-velocity data from the ground-based SONG network with continuous space-based data from the revised Kepler mission K2, we derive and characterize 27 individual oscillation modes for epsilon Tau, along with global oscillation parameters such as the large frequency separation and the ratio between the amplitude of the oscillations measured in radial velocity and intensity as a function of frequency. The latter has been measured previously for only two stars, the Sun and Procyon. Combining the seismic analysis with interferometric and spectroscopic measurements, we derive physical parameters for epsilon Tau, and discuss its evolutionary status. △ Less

Submitted 22 February, 2019; v1 submitted 18 January, 2019; originally announced January 2019.

Comments: 13 pages, 13 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

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

Anisotropic winds in Wolf-Rayet binary identify potential gamma-ray burst progenitor

Authors: J. R. Callingham, P. G. Tuthill, B. J. S. Pope, P. M. Williams, P. A. Crowther, M. Edwards, B. Norris, L. Kedziora-Chudczer

Abstract: The massive evolved Wolf-Rayet stars sometimes occur in colliding-wind binary systems in which dust plumes are formed as a result of the collision of stellar winds. These structures are known to encode the parameters of the binary orbit and winds. Here, we report observations of a previously undiscovered Wolf-Rayet system, 2XMM J160050.7-514245, with a spectroscopically determined wind speed of… ▽ More The massive evolved Wolf-Rayet stars sometimes occur in colliding-wind binary systems in which dust plumes are formed as a result of the collision of stellar winds. These structures are known to encode the parameters of the binary orbit and winds. Here, we report observations of a previously undiscovered Wolf-Rayet system, 2XMM J160050.7-514245, with a spectroscopically determined wind speed of $\approx$3400 km s$^{-1}$. In the thermal infrared, the system is adorned with a prominent $\approx$12$''$ spiral dust plume, revealed by proper motion studies to be expanding at only $\approx$570 km s$^{-1}$. As the dust and gas appear coeval, these observations are inconsistent with existing models of the dynamics of such colliding wind systems. We propose that this contradiction can be resolved if the system is capable of launching extremely anisotropic winds. Near-critical stellar rotation is known to drive such winds, suggesting this Wolf-Rayet system as a potential Galactic progenitor system for long-duration gamma-ray bursts. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: Accepted to Nature Astronomy on 2018 Oct 1. 27 pages, 9 figures, and 4 tables

Exoplanet Transits with Next-Generation Radio Telescopes

Authors: Benjamin J. S. Pope, Paul Withers, Joseph R. Callingham, Marissa F. Vogt

Abstract: Nearly everything we know about extrasolar planets to date comes from optical astronomy. While exoplanetary aurorae are predicted to be bright at low radio frequencies (< 1 GHz), we consider the effect of an exoplanet transit on radio emission from the host star. As radio emission from solar-like stars is concentrated in active regions, a planet occulting a starspot can cause a disproportionately… ▽ More Nearly everything we know about extrasolar planets to date comes from optical astronomy. While exoplanetary aurorae are predicted to be bright at low radio frequencies (< 1 GHz), we consider the effect of an exoplanet transit on radio emission from the host star. As radio emission from solar-like stars is concentrated in active regions, a planet occulting a starspot can cause a disproportionately deep transit which should be detectable with major radio arrays currently under development, such as the Square Kilometre Array (SKA). We calculate the radiometric sensitivity of the SKA stages and components, finding that SKA2-Mid can expect to detect transits around the very nearest solar-like stars and many cool dwarfs. The shape of this radiometric light curve will be affected by scintillation and lensing from the planet's magnetosphere and thereby encode magnetospheric parameters. Furthermore, these transits will also probe the distribution of stellar activity across a star's surface, and will help scrub out contamination from stellar activity on exoplanet transmission spectra and radial velocity spectra. This radio window on exoplanets and their host stars is therefore a valuable complement to existing optical tools. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: 11 pages, 2 figures. Under review at MNRAS

EarthFinder: A Precise Radial Velocity Probe Mission Concept For the Detection of Earth-Mass Planets Orbiting Sun-like Stars

Authors: Peter Plavchan, Bryson Cale, Patrick Newman, Bahaa Hamze, Natasha Latouf, William Matzko, Chas Beichman, David Ciardi, Bill Purcell, Paul Lightsey, Heather Cegla, Xavier Dumusque, Vincent Bourrier, Courtney Dressing, Peter Gao, Gautam Vasisht, Stephanie Leifer, Sharon Wang, Jonathan Gagne, Samantha Thompson, Jonathan Crass, Andrew Bechter, Eric Bechter, Cullen Blake, Sam Halverson , et al. (18 additional authors not shown)

Abstract: EarthFinder is a Probe Mission concept selected for study by NASA for input to the 2020 astronomy decadal survey. This study is currently active and a final white paper report is due to NASA at the end of calendar 2018. We are tasked with evaluating the scientific rationale for obtaining precise radial velocity (PRV) measurements in space, which is a two-part inquiry: What can be gained from going… ▽ More EarthFinder is a Probe Mission concept selected for study by NASA for input to the 2020 astronomy decadal survey. This study is currently active and a final white paper report is due to NASA at the end of calendar 2018. We are tasked with evaluating the scientific rationale for obtaining precise radial velocity (PRV) measurements in space, which is a two-part inquiry: What can be gained from going to space? What can't be done form the ground? These two questions flow down to these specific tasks for our study - Identify the velocity limit, if any, introduced from micro- and macro-telluric absorption in the Earth's atmosphere; Evaluate the unique advantages that a space-based platform provides to emable the identification and mitigation of stellar acitivity for multi-planet signal recovery. △ Less

Submitted 11 March, 2018; originally announced March 2018.

Comments: Submitted to the National Academies Committee on Exoplanet Science Strategy

Aldebaran b's temperate past uncovered in planet search data

Authors: Will M. Farr, Benjamin J. S. Pope, Guy R. Davies, Thomas S. H. North, Timothy R. White, Jim W. Barrett, Andrea Miglio, Mikkel N. Lund, Victoria Antoci, Mads Fredslund Andersen, Frank Grundahl, Daniel Huber

Abstract: The nearby red giant Aldebaran is known to host a gas giant planetary companion from decades of ground-based spectroscopic radial velocity measurements. Using Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) models, we show that these historic data also contain evidence of acoustic oscillations in the star itself, and verify this result with further dedicated ground-based s… ▽ More The nearby red giant Aldebaran is known to host a gas giant planetary companion from decades of ground-based spectroscopic radial velocity measurements. Using Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) models, we show that these historic data also contain evidence of acoustic oscillations in the star itself, and verify this result with further dedicated ground-based spectroscopy and space-based photometry with the Kepler Space Telescope. From the frequency of these oscillations we determine the mass of Aldebaran to be $1.16 \pm 0.07 \, M_\odot$, and note that this implies its planet will have been subject to insolation comparable to the Earth for some of the star's main sequence lifetime. Our approach to sparse, irregularly sampled time series astronomical observations has the potential to unlock asteroseismic measurements for thousands of stars in archival data, and push to lower-mass planets around red giant stars. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: 24 pages, 7 figures (including appendices); submitted to ApJL; paper text, figures, data, and code at https://github.com/farr/Aldebaran

Beyond the Kepler/K2 bright limit: variability in the seven brightest members of the Pleiades

Authors: T. R. White, B. J. S. Pope, V. Antoci, P. I. Pápics, C. Aerts, D. R. Gies, K. Gordon, D. Huber, G. H. Schaefer, S. Aigrain, S. Albrecht, T. Barclay, G. Barentsen, P. G. Beck, T. R. Bedding, M. Fredslund Andersen, F. Grundahl, S. B. Howell, M. J. Ireland, S. J. Murphy, M. B. Nielsen, V. Silva Aguirre, P. G. Tuthill

Abstract: The most powerful tests of stellar models come from the brightest stars in the sky, for which complementary techniques, such as astrometry, asteroseismology, spectroscopy, and interferometry can be combined. The K2 Mission is providing a unique opportunity to obtain high-precision photometric time series for bright stars along the ecliptic. However, bright targets require a large number of pixels… ▽ More The most powerful tests of stellar models come from the brightest stars in the sky, for which complementary techniques, such as astrometry, asteroseismology, spectroscopy, and interferometry can be combined. The K2 Mission is providing a unique opportunity to obtain high-precision photometric time series for bright stars along the ecliptic. However, bright targets require a large number of pixels to capture the entirety of the stellar flux, and bandwidth restrictions limit the number and brightness of stars that can be observed. To overcome this, we have developed a new photometric technique, that we call halo photometry, to observe very bright stars using a limited number of pixels. Halo photometry is simple, fast and does not require extensive pixel allocation, and will allow us to use K2 and other photometric missions, such as TESS, to observe very bright stars for asteroseismology and to search for transiting exoplanets. We apply this method to the seven brightest stars in the Pleiades open cluster. Each star exhibits variability; six of the stars show what are most-likely slowly pulsating B-star (SPB) pulsations, with amplitudes ranging from 20 to 2000 ppm. For the star Maia, we demonstrate the utility of combining K2 photometry with spectroscopy and interferometry to show that it is not a 'Maia variable', and to establish that its variability is caused by rotational modulation of a large chemical spot on a 10 d time scale. △ Less

Submitted 24 August, 2017; originally announced August 2017.

Comments: Published in Monthly Notices of the Royal Astronomical Society. 20 pages, 10 figures and 10 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, 471, p.2882-2901 (2017)

A Low-Mass Exoplanet Candidate Detected By ${\it K2}$ Transiting the Praesepe M Dwarf JS 183

Authors: Joshua Pepper, Ed Gillen, Hannu Parviainen, Lynne A. Hillenbrand, Ann Marie Cody, Suzanne Aigrain, John Stauffer, Frederick J. Vrba, Trevor David, Jorge Lillo-Box, Keivan G. Stassun, Kyle E. Conroy, Benjamin J. S. Pope, David Barrado

Abstract: We report the discovery of a repeating photometric signal from a low-mass member of the Praesepe open cluster that we interpret as a Neptune-sized transiting planet. The star is JS 183 (HSHJ 163, EPIC 211916756) with $T_{\rm eff} = 3325\pm100$ K, $M_{*} = 0.44\pm0.04$ $M_{\odot}$, $R_{*} = 0.44\pm0.03$ $R_{\odot}$, and $\log{g_*} = 4.82\pm0.06$. The planet has an orbital period of 10.134588 days a… ▽ More We report the discovery of a repeating photometric signal from a low-mass member of the Praesepe open cluster that we interpret as a Neptune-sized transiting planet. The star is JS 183 (HSHJ 163, EPIC 211916756) with $T_{\rm eff} = 3325\pm100$ K, $M_{*} = 0.44\pm0.04$ $M_{\odot}$, $R_{*} = 0.44\pm0.03$ $R_{\odot}$, and $\log{g_*} = 4.82\pm0.06$. The planet has an orbital period of 10.134588 days and a radius of $R_{P}= 0.32\pm0.02$ $R_J$. Since the star is faint at $V=16.5$ and $J=13.3$, we are unable to obtain a measured radial-velocity orbit, but we can constrain the companion mass to below about 1.7 $M_J$, and thus well below the planetary boundary. JS 183b (since designated as K2-95b) is the second transiting planet found with ${\it K2}$ that resides in a several hundred Myr open cluster; both planets orbit mid-M dwarf stars and are approximately Neptune-sized. With a well-determined stellar density from the planetary transit, and with an independently known metallicity from its cluster membership, JS 183 provides a particularly valuable test of stellar models at the fully convective boundary. We find that JS 183 is the lowest-density transit host known at the fully convective boundary, and that its very low density is consistent with current models of stars just above the fully convective boundary but in tension with the models just below the fully convective boundary. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Comments: 30 pages, 11 figures, Accepted for publication in the Astronomical Journal

Journal ref: The Astronomical Journal, Volume 153, Number 4, 2017

Testing stellar evolution models with the retired A star HD 185351

Authors: Jakob G. Hjørringgaard, Victor Silva Aguirre, Tim R. White, Daniel Huber, Benjamin J. S. Pope, Luca Casagrande, Anders B. Justesen, Jørgen Christensen-Dalsgaard

Abstract: The physical parameters of the retired A star HD 185351 were analysed in great detail by Johnson et al. (2014) using interferometry, spectroscopy and asteroseismology. Results from all independent methods are consistent with HD 185351 having a mass in excess of $1.5\mathrm{M}_{\odot}$. However, the study also showed that not all observational constraints could be reconciled in stellar evolutionary… ▽ More The physical parameters of the retired A star HD 185351 were analysed in great detail by Johnson et al. (2014) using interferometry, spectroscopy and asteroseismology. Results from all independent methods are consistent with HD 185351 having a mass in excess of $1.5\mathrm{M}_{\odot}$. However, the study also showed that not all observational constraints could be reconciled in stellar evolutionary models, leading to mass estimates ranging from $\sim 1.6-1.9\mathrm{M}_{\odot}$ and casting doubts on the accuracy of stellar properties determined from asteroseismology. Here we solve this discrepancy and construct a theoretical model in agreement with all observational constraints on the physical parameters of HD 185351. The effects of varying input physics are examined as well as considering the additional constraint of the observed g-mode period spacing. This quantity is found to be sensitive to the inclusion of additional mixing from the convective core during the main sequence, and can be used to calibrate the overshooting efficiency using low-luminosity red giant stars. A theoretical model with metallicity $\left[\mathrm{Fe/H}\right]=0.16$dex, mixing-length parameter $α_{\mathrm{MLT}}=2.00$, and convective overshooting efficiency parameter $f=0.030$ is found to be in complete agreement with all observational constraints for a stellar mass of $M\simeq1.60\mathrm{M}_{\odot}$. △ Less

Submitted 19 October, 2016; originally announced October 2016.

Comments: 7 pages, 7 figures. To be published in mnras

Kernel Phase and Kernel Amplitude in Fizeau Imaging

Authors: Benjamin J. S. Pope

Abstract: Kernel phase interferometry is an approach to high angular resolution imaging which enhances the performance of speckle imaging with adaptive optics. Kernel phases are self-calibrating observables that generalize the idea of closure phases from non-redundant arrays to telescopes with arbitrarily shaped pupils, by considering a matrix-based approximation to the diffraction problem. In this paper I… ▽ More Kernel phase interferometry is an approach to high angular resolution imaging which enhances the performance of speckle imaging with adaptive optics. Kernel phases are self-calibrating observables that generalize the idea of closure phases from non-redundant arrays to telescopes with arbitrarily shaped pupils, by considering a matrix-based approximation to the diffraction problem. In this paper I discuss the recent history of kernel phase, in particular in the matrix-based study of sparse arrays, and propose an analogous generalization of the closure amplitude to kernel amplitudes. This new approach can self-calibrate throughput and scintillation errors in optical imaging, which extends the power of kernel phase-like methods to symmetric targets where amplitude and not phase calibration can be a significant limitation, and will enable further developments in high angular resolution astronomy. △ Less

Submitted 1 September, 2016; originally announced September 2016.

Comments: Accepted MNRAS, 10 pages

Supernovae and Single-Year Anomalies in the Atmospheric Radiocarbon Record

Authors: Michael Dee, Benjamin Pope, Daniel Miles, Sturt Manning, Fusa Miyake

Abstract: Single-year spikes in radiocarbon production are caused by intense bursts of radiation from space. Supernovae emit both high-energy particle and electromagnetic radiation, but it is the latter that is most likely to strike the atmosphere all at once and cause a surge in 14C production. In the 1990s, it was claimed that the supernova in 1006 CE produced exactly this effect. With the 14C spikes in t… ▽ More Single-year spikes in radiocarbon production are caused by intense bursts of radiation from space. Supernovae emit both high-energy particle and electromagnetic radiation, but it is the latter that is most likely to strike the atmosphere all at once and cause a surge in 14C production. In the 1990s, it was claimed that the supernova in 1006 CE produced exactly this effect. With the 14C spikes in the years 775 and 994 CE now attributed to extreme solar events, attention has returned to the question of whether historical supernovae are indeed detectable using annual 14C measurements. Here, we combine new and existing measurements over six documented and putative supernovae, and conclude that no such astrophysical event has yet left a distinct imprint on the past atmospheric 14C record. △ Less

Submitted 7 August, 2016; originally announced August 2016.

Comments: 10 Pages; Radiocarbon, accepted

Journal ref: Radiocarbon / FirstView Article / August 2016, pp 1 - 10

Planet Formation Imager (PFI): science vision and key requirements

Authors: Stefan Kraus, John D. Monnier, Michael J. Ireland, Gaspard Duchene, Catherine Espaillat, Sebastian Hoenig, Attila Juhasz, Chris Mordasini, Johan Olofsson, Claudia Paladini, Keivan Stassun, Neal Turner, Gautam Vasisht, Tim J. Harries, Matthew R. Bate, Jean-Francois Gonzalez, Alexis Matter, Zhaohuan Zhu, Olja Panic, Zsolt Regaly, Alessandro Morbidelli, Farzana Meru, Sebastian Wolf, John Ilee, Jean-Philippe Berger , et al. (53 additional authors not shown)

Abstract: The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the prot… ▽ More The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs. △ Less

Submitted 16 August, 2016; v1 submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 5 figures, Proceedings of SPIE 2016

Transiting exoplanet candidates from K2 Campaigns 5 and 6

Authors: Benjamin J. S. Pope, Hannu Parviainen, Suzanne Aigrain

Abstract: We introduce a new transit search and vetting pipeline for observations from the K2 mission, and present the candidate transiting planets identified by this pipeline out of the targets in Campaigns 5 and 6. Our pipeline uses the Gaussian Process-based K2SC code to correct for the K2 pointing systematics and simultaneously model stellar variability. The systematics-corrected, variability-detrended… ▽ More We introduce a new transit search and vetting pipeline for observations from the K2 mission, and present the candidate transiting planets identified by this pipeline out of the targets in Campaigns 5 and 6. Our pipeline uses the Gaussian Process-based K2SC code to correct for the K2 pointing systematics and simultaneously model stellar variability. The systematics-corrected, variability-detrended light curves are searched for transits with the Box Least Squares method, and a period-dependent detection threshold is used to generate a preliminary candidate list. Two or three individuals vet each candidate manually to produce the final candidate list, using a set of automatically-generated transit fits and assorted diagnostic tests to inform the vetting. We detect 147 single-planet system candidates and 5 multi-planet systems, independently recovering the previously-published hot~Jupiters EPIC 212110888b, WASP-55b (EPIC 212300977b) and Qatar-2b (EPIC 212756297b). We also report the outcome of reconnaissance spectroscopy carried out for all candidates with Kepler magnitude Kp < 13, identifying 12 targets as likely false positives. We compare our results to those of other K2 transit search pipelines, noting that ours performs particularly well for variable and/or active stars, but that the results are very similar overall. All the light curves and code used in the transit search and vetting process are publicly available, as are the follow-up spectra. △ Less

Submitted 3 June, 2016; originally announced June 2016.

Comments: Accepted to MNRAS

K2SC: Flexible systematics correction and detrending of K2 light curves using Gaussian Process regression

Authors: Suzanne Aigrain, Hannu Parviainen, Benjamin Pope

Abstract: We present K2SC (K2 Systematics Correction), a Python pipeline to model instrumental systematics and astrophysical variability in light curves from the K2 mission. K2SC uses Gaussian process regression to model position-dependent systematics and time-dependent variability simultaneously, enabling the user to remove both (e.g., for transit searches) or to remove systematics while preserving variabi… ▽ More We present K2SC (K2 Systematics Correction), a Python pipeline to model instrumental systematics and astrophysical variability in light curves from the K2 mission. K2SC uses Gaussian process regression to model position-dependent systematics and time-dependent variability simultaneously, enabling the user to remove both (e.g., for transit searches) or to remove systematics while preserving variability (for variability studies). For periodic variables, K2SC automatically computes estimates of the period, amplitude and evolution timescale of the variability. We apply K2SC to publicly available K2 data from campaigns 3--5, showing that we obtain photometric precision approaching that of the original Kepler mission. We compare our results to other publicly available K2 pipelines, showing that we obtain similar or better results, on average. We use transit injection and recovery tests to evaluate the impact of K2SC on planetary transit searches in K2 PDC (Pre-search Data Conditioning) data, for planet-to-star radius ratios down Rp/Rstar = 0.01 and periods up to P = 40 d, and show that K2SC significantly improves the ability to distinguish between correct and false detections, particularly for small planets. K2SC can be run automatically on many light curves, or manually tailored for specific objects such as pulsating stars or large amplitude eclipsing binaries. It can be run on ASCII and FITS light curve files, regardless of their origin. Both the code and the processed light curves are publicly available, and we provide instructions for downloading and using them. The methodology used by K2SC will be applicable to future transit search missions such as TESS and PLATO. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS. The processed data for campaigns 3 to 6 will be available shortly at https://archive.stsci.edu/prepds/k2sc/ and the K2SC code is available at https://github.com/OxES/k2sc

The Palomar Kernel Phase Experiment: Testing Kernel Phase Interferometry for Ground-based Astronomical Observations

Authors: Benjamin Pope, Peter Tuthill, Sasha Hinkley, Michael J. Ireland, Alexandra Greenbaum, Alexey Latyshev, John D. Monnier, Frantz Martinache

Abstract: At present, the principal limitation on the resolution and contrast of astronomical imaging instruments comes from aberrations in the optical path, which may be imposed by the Earth's turbulent atmosphere or by variations in the alignment and shape of the telescope optics. These errors can be corrected physically, with active and adaptive optics, and in post-processing of the resulting image. A re… ▽ More At present, the principal limitation on the resolution and contrast of astronomical imaging instruments comes from aberrations in the optical path, which may be imposed by the Earth's turbulent atmosphere or by variations in the alignment and shape of the telescope optics. These errors can be corrected physically, with active and adaptive optics, and in post-processing of the resulting image. A recently-developed adaptive optics post-processing technique, called kernel phase interferometry, uses linear combinations of phases that are self-calibrating with respect to small errors, with the goal of constructing observables that are robust against the residual optical aberrations in otherwise well-corrected imaging systems. Here we present a direct comparison between kernel phase and the more established competing techniques, aperture masking interferometry, point spread function (PSF) fitting and bispectral analysis. We resolve the alpha Ophiuchi binary system near periastron, using the Palomar 200-Inch Telescope. This is the first case in which kernel phase has been used with a full aperture to resolve a system close to the diffraction limit with ground-based extreme adaptive optics observations. Excellent agreement in astrometric quantities is found between kernel phase and masking, and kernel phase significantly outperforms PSF fitting and bispectral analysis, demonstrating its viability as an alternative to conventional non-redundant masking under appropriate conditions. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: Accepted to MNRAS

HII 2407: A Low-Mass Eclipsing Binary Revealed by K2 Observations of the Pleiades

Authors: Trevor J. David, John Stauffer, Lynne A. Hillenbrand, Ann Marie Cody, Kyle Conroy, Keivan G. Stassun, Benjamin Pope, Suzanne Aigrain, Ed Gillen, Andrew Collier Cameron, David Barrado, L. M. Rebull, Howard Isaacson, Geoffrey W. Marcy, Celia Zhang, Reed L. Riddle, Carl Ziegler, Nicholas M. Law, Christoph Baranec

Abstract: The star HII 2407 is a member of the relatively young Pleiades star cluster and was previously discovered to be a single-lined spectroscopic binary. It is newly identified here within $Kepler$/$K2$ photometric time series data as an eclipsing binary system. Mutual fitting of the radial velocity and photometric data leads to an orbital solution and constraints on fundamental stellar parameters. Whi… ▽ More The star HII 2407 is a member of the relatively young Pleiades star cluster and was previously discovered to be a single-lined spectroscopic binary. It is newly identified here within $Kepler$/$K2$ photometric time series data as an eclipsing binary system. Mutual fitting of the radial velocity and photometric data leads to an orbital solution and constraints on fundamental stellar parameters. While the primary has arrived on the main sequence, the secondary is still pre-main-sequence and we compare our results for the $M/M_\odot$ and $R/R_\odot$ values with stellar evolutionary models. We also demonstrate that the system is likely to be tidally synchronized. Follow-up infrared spectroscopy is likely to reveal the lines of the secondary, allowing for dynamically measured masses and elevating the system to benchmark eclipsing binary status. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: 9 pages, 4 figures, accepted for publication in ApJ

Photometry of Very Bright Stars with Kepler and K2 Smear Data

Authors: Benjamin Pope, Timothy White, Daniel Huber, Simon Murphy, Tim Bedding, Douglas Caldwell, Aleksa Sarai, Suzanne Aigrain, Thomas Barclay

Abstract: High-precision time series photometry with the Kepler satellite has been crucial to our understanding both of exoplanets, and via asteroseismology, of stellar physics. After the failure of two reaction wheels, the Kepler satellite has been repurposed as Kepler-2 (K2), observing fields close to the ecliptic plane. As these fields contain many more bright stars than the original Kepler field, K2 pro… ▽ More High-precision time series photometry with the Kepler satellite has been crucial to our understanding both of exoplanets, and via asteroseismology, of stellar physics. After the failure of two reaction wheels, the Kepler satellite has been repurposed as Kepler-2 (K2), observing fields close to the ecliptic plane. As these fields contain many more bright stars than the original Kepler field, K2 provides an unprecedented opportunity to study nearby objects amenable to detailed follow-up with ground-based instruments. Due to bandwidth constraints, only a small fraction of pixels can be downloaded, with the result that most bright stars which saturate the detector are not observed. We show that engineering data acquired for photometric calibration, consisting of collateral `smear' measurements, can be used to reconstruct light curves for bright targets not otherwise observable with Kepler/K2. Here we present some examples from Kepler Quarter 6 and K2 Campaign 3, including the delta Scuti variables HD 178875 and 70 Aqr, and the red giant HR 8500 displaying solar-like oscillations. We compare aperture and smear photometry where possible, and also study targets not previously observed. These encouraging results suggest this new method can be applied to most Kepler and K2 fields. △ Less

Submitted 30 September, 2015; originally announced October 2015.

Comments: Accepted to MNRAS Letters, 5 pages, 5 figures

Wavefront sensing from the image domain with the Oxford-SWIFT integral field spectrograph

Authors: Benjamin Pope, Niranjan Thatte, Rick Burruss, Matthias Tecza, Fraser Clarke, Garret Cotter

Abstract: The limits for adaptive-optics (AO) imaging at high contrast and high resolution are determined by residual phase errors from non-common-path aberrations not sensed by the wavefront sensor, especially for integral field spectrographs, where phase diversity techniques are complicated by the image slicer. We present the first application of kernel phase-based wavefront sensing to ground-based AO, wh… ▽ More The limits for adaptive-optics (AO) imaging at high contrast and high resolution are determined by residual phase errors from non-common-path aberrations not sensed by the wavefront sensor, especially for integral field spectrographs, where phase diversity techniques are complicated by the image slicer. We present the first application of kernel phase-based wavefront sensing to ground-based AO, where an asymmetric pupil mask and a single image are sufficient to map aberrations up to high order. We push toward internally diffraction-limited performance with the Oxford-SWIFT integral field spectrograph coupled with the PALM-3000 extreme AO system on the Palomar 200-inch telescope. This represents the first observation in which the PALM-3000 + SWIFT internal point-spread-function has closely approached the Airy pattern. While this can only be used on SWIFT with an internal stimulus source, as at short wavelengths the uncorrected atmospheric wavefront errors are still > 1 radian, this nevertheless demonstrates the feasibility of detecting non-common-path errors with this method as an active optics paradigm for near-infrared, AO-corrected instruments at Palomar. We note that this is a particularly promising approach for correcting integral field spectrographs, as the diversity of many narrowband images provides strong constraints on the wavefront error estimate, and the average of reconstructions from many narrow bands can be used to improve overall reconstruction quality. △ Less

Submitted 2 July, 2014; originally announced July 2014.

Comments: 8 Pages, 3 Figures, submitted to SPIE Astronomical Telescopes and Instrumentation 2014

A Demonstration of Wavefront Sensing and Mirror Phasing from the Image Domain

Authors: Benjamin Pope, Nick Cvetojevic, Anthony Cheetham, Frantz Martinache, Barnaby Norris, Peter Tuthill

Abstract: In astronomy and microscopy, distortions in the wavefront affect the dynamic range of a high contrast imaging system. These aberrations are either imposed by a turbulent medium such as the atmosphere, by static or thermal aberrations in the optical path, or by imperfectly phased subapertures in a segmented mirror. Active and adaptive optics (AO), consisting of a wavefront sensor and a deformable m… ▽ More In astronomy and microscopy, distortions in the wavefront affect the dynamic range of a high contrast imaging system. These aberrations are either imposed by a turbulent medium such as the atmosphere, by static or thermal aberrations in the optical path, or by imperfectly phased subapertures in a segmented mirror. Active and adaptive optics (AO), consisting of a wavefront sensor and a deformable mirror, are employed to address this problem. Nevertheless, the non-common-path between the wavefront sensor and the science camera leads to persistent quasi-static speckles that are difficult to calibrate and which impose a floor on the image contrast. In this paper we present the first experimental demonstration of a novel wavefront sensor requiring only a minor asymmetric obscuration of the pupil, using the science camera itself to detect high order wavefront errors from the speckle pattern produced. We apply this to correct errors imposed on a deformable microelectromechanical (MEMS) segmented mirror in a closed loop, restoring a high quality point spread function (PSF) and residual wavefront errors of order $\sim 10$ nm using 1600 nm light, from a starting point of $\sim 300$ nm in piston and $\sim 0.3$ mrad in tip-tilt. We recommend this as a method for measuring the non-common-path error in AO-equipped ground based telescopes, as well as as an approach to phasing difficult segmented mirrors such as on the \emph{James Webb Space Telescope} primary and as a future direction for extreme adaptive optics. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: 9 pages, 6 figures