Showing 1–50 of 58 results for author: Howard, S

Giant planet interiors and atmospheres

Authors: Ravit Helled, Saburo Howard

Abstract: Studying the interiors of the outer planets is crucial for a comprehensive understanding of our planetary system, and provides key knowledge on the origin of the solar system, the behavior of materials at extreme conditions, the relation between the planetary interior and atmosphere, and exoplanet characterization. In this review, we summarize the current understanding of the interiors and atmosph… ▽ More Studying the interiors of the outer planets is crucial for a comprehensive understanding of our planetary system, and provides key knowledge on the origin of the solar system, the behavior of materials at extreme conditions, the relation between the planetary interior and atmosphere, and exoplanet characterization. In this review, we summarize the current understanding of the interiors and atmospheres of the giant planets in the solar system: Jupiter, Saturn, Uranus and Neptune. We describe the principles of interior structure and evolution models and discuss the link between the planetary atmospheres and deep interiors. We summarize the current understanding of the bulk compositions and internal structures of the outer planets and the remaining future challenges. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Preprint of a chapter for the 'Encyclopedia of Astrophysics' (Editor-in-Chief Ilya Mandel, Section Editor Dimitri Veras) to be published by Elsevier as a Reference Module

High-energy spectra of LTT 1445A and GJ 486 reveal flares and activity

Authors: H. Diamond-Lowe, G. W. King, A. Youngblood, A. Brown, W. S. Howard, J. G. Winters, D. J. Wilson, K. France, J. M. Mendonça, L. A. Buchhave, L. Corrales, L. Kreidberg, A. A. Medina, J. L. Bean, Z. K. Berta-Thompson, T. M. Evans-Soma, C. Froning, G. M. Duvvuri, E. M. -R. Kempton, Y. Miguel, J. S. Pineda, C. Schneider

Abstract: The high-energy radiative output, from the X-ray to the ultraviolet, of exoplanet host stars drives photochemical reactions and mass loss in the upper regions of planetary atmospheres. In order to place constraints on the atmospheric properties of the three closest terrestrial exoplanets transiting M dwarfs, we observe the high-energy spectra of the host stars LTT1445A and GJ486 in the X-ray with… ▽ More The high-energy radiative output, from the X-ray to the ultraviolet, of exoplanet host stars drives photochemical reactions and mass loss in the upper regions of planetary atmospheres. In order to place constraints on the atmospheric properties of the three closest terrestrial exoplanets transiting M dwarfs, we observe the high-energy spectra of the host stars LTT1445A and GJ486 in the X-ray with XMM-Newton and Chandra and in the ultraviolet with HST/COS and STIS. We combine these observations with estimates of extreme ultraviolet flux, reconstructions of the Ly-a lines, and stellar models at optical and infrared wavelengths to produce panchromatic spectra from 1A--20um for each star. While LTT1445Ab, LTT1445Ac, and GJ486b do not possess primordial hydrogen-dominated atmospheres, we calculate that they are able to retain pure CO2 atmospheres if starting with 10, 15, and 50% of Earth's total CO2 budget, respectively, in the presence of their host stars' stellar wind. We use age-activity relationships to place lower limits of 2.2 and 6.6 Gyr on the ages of the host stars LTT1445A and GJ486. Despite both LTT1445A and GJ486 appearing inactive at optical wavelengths, we detect flares at ultraviolet and X-ray wavelengths for both stars. In particular, GJ486 exhibits two flares with absolute energies of 10^29.5 and 10^30.1 erg (equivalent durations of 4357+/-96 and 19724+/-169 s) occurring three hours apart, captured with HST/COS G130M. Based on the timing of the observations, we suggest that these high-energy flares are related and indicative of heightened flaring activity that lasts for a period of days, but our interpretations are limited by sparse time-sampling. Consistent high-energy monitoring is needed to determine the duration and extent of high-energy activity on individual M dwarfs, as well as the population as a whole. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 22 pages, accepted to A&A

NeuralCMS: A deep learning approach to study Jupiter's interior

Authors: Maayan Ziv, Eli Galanti, Amir Sheffer, Saburo Howard, Tristan Guillot, Yohai Kaspi

Abstract: NASA's Juno mission provided exquisite measurements of Jupiter's gravity field that together with the Galileo entry probe atmospheric measurements constrains the interior structure of the giant planet. Inferring its interior structure range remains a challenging inverse problem requiring a computationally intensive search of combinations of various planetary properties, such as the cloud-level tem… ▽ More NASA's Juno mission provided exquisite measurements of Jupiter's gravity field that together with the Galileo entry probe atmospheric measurements constrains the interior structure of the giant planet. Inferring its interior structure range remains a challenging inverse problem requiring a computationally intensive search of combinations of various planetary properties, such as the cloud-level temperature, composition, and core features, requiring the computation of ~10^9 interior models. We propose an efficient deep neural network (DNN) model to generate high-precision wide-ranged interior models based on the very accurate but computationally demanding concentric MacLaurin spheroid (CMS) method. We trained a sharing-based DNN with a large set of CMS results for a four-layer interior model of Jupiter, including a dilute core, to accurately predict the gravity moments and mass, given a combination of interior features. We evaluated the performance of the trained DNN (NeuralCMS) to inspect its predictive limitations. NeuralCMS shows very good performance in predicting the gravity moments, with errors comparable with the uncertainty due to differential rotation, and a very accurate mass prediction. This allowed us to perform a broad parameter space search by computing only ~10^4 actual CMS interior models, resulting in a large sample of plausible interior structures, and reducing the computation time by a factor of 10^5. Moreover, we used a DNN explainability algorithm to analyze the impact of the parameters setting the interior model on the predicted observables, providing information on their nonlinear relation. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 8 pages, 6 figures, 4 tables, accepted for publication in A&A

Spiral Galaxies' Disk Dominance and Testing MOND

Authors: Gene Byrd, Sethanne Howard

Abstract: Many Hubble type S galaxies have a flat rotation curve, V extending into the outer disk. The surface brightness there is too low to create V from stars of a reasonable mass-to-light ratio. To maintain a stable disk, a massive dark matter halo is usually assumed to dominate in creating V. Gravitational arm amplification is used here to estimate whether V is created by a disk mixture, or by a 100 pe… ▽ More Many Hubble type S galaxies have a flat rotation curve, V extending into the outer disk. The surface brightness there is too low to create V from stars of a reasonable mass-to-light ratio. To maintain a stable disk, a massive dark matter halo is usually assumed to dominate in creating V. Gravitational arm amplification is used here to estimate whether V is created by a disk mixture, or by a 100 percent halo. The disk surface density is not estimated from the disk surface brightness and an assumed mass to light ratio. The observed arm pitch angle is presumed to be long term maximally amplified. This permits calculation of the disk surface mass density with radius relative to a 100 percent disk which would fully create V. NGC 7217's tightly wound 4.8 degree arm pitch angle indicates a disk component contributing 43 percent of the observed V for the range over the spiral arms (important but not dominant). It has too low a disk density to create V with the rest presumably arising from a halo. V is proportional to the square root of surface density. By contrast a loose arm galaxy disk can dominate in creating V. For example NGC 3198's arm pitch angle of 30 degrees implies a disk surface density sufficient to account for 88 percent of V which is a dominant fraction. Little dark matter halo is needed for galaxies like this. Observationally, the dark outer disk may be molecular hydrogen. NGC 3198's dominant disk also provides a counter example to the need for any universal (MOND) force law modification over the outer disk spiral arms. △ Less

Submitted 1 April, 2024; originally announced April 2024.

The ArgusSpec Prototype: Autonomous Spectroscopic Follow-up of Flares Detected by Large Array Telescopes

Authors: Nathan W. Galliher, Thomas Procter, Nicholas M. Law, Hank Corbett, Ward S. Howard, Alan Vasquez Soto, Ramses Gonzalez, Lawrence Machia, Jonathan Carney, William J. Marshall

Abstract: ArgusSpec is a prototype autonomous spectroscopic follow-up instrument designed to characterize flares detected by the Argus Pathfinder telescope array by taking short exposure (30 s) broadband spectra (370 - 750 nm) at low resolutions (R~150 at 500 nm). The instrument is built from consumer off-the-shelf astronomical equipment, assembled inside a shipping container, and deployed alongside the Arg… ▽ More ArgusSpec is a prototype autonomous spectroscopic follow-up instrument designed to characterize flares detected by the Argus Pathfinder telescope array by taking short exposure (30 s) broadband spectra (370 - 750 nm) at low resolutions (R~150 at 500 nm). The instrument is built from consumer off-the-shelf astronomical equipment, assembled inside a shipping container, and deployed alongside the Argus Pathfinder at a dark sky observing site in Western North Carolina. The \$35k prototype ArgusSpec was designed, built, and deployed in under a year, largely from existing parts, and has been operating on-sky since March 2023. With current hardware and software, the system is capable of receiving an observation, slewing, performing autonomous slit acquisition, and beginning data acquisition within an average of 32 s. With Argus Pathfinder's 1-second-cadence survey reporting alerts of rising sources within 2 s of onset, ArgusSpec can reach new targets well within a minute of the start of the event. As built, ArgusSpec can observe targets down to a 20$σ$ limiting magnitude of $m_V$~13 at 30 s cadence with an optical resolution of R~150 (at 500 nm). With automated rapid acquisition demonstrated, later hardware upgrades will significantly improve the limiting magnitude, and potentially enable deep spectroscopy by the coaddition of data from an array of ArgusSpec systems. ArgusSpec's primary science driver is the characterization of the blackbody evolution of flares from nearby M-dwarfs. Large flares emitted by these stars could have significant impacts on the potential habitability of any orbiting exoplanets, but our current understanding of these events is in large part built on observations from a handful of active stars. ArgusSpec will characterize large numbers of flares, building a spectroscopic library of the most extreme events from a wide variety of stellar masses and ages. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 19 pages, 12 figures, 2 tables, Accepted to PASP

Three-dimensional atmospheric dynamics of Jupiter from ground-based Doppler imaging spectroscopy in the visible

Authors: François-Xavier Schmider, Patrick Gaulme, Raúl Morales-Juberías, Jason Jackiewicz, Ivan Gonçalves, Tristan Guillot, Amy A. Simon, Michael H. Wong, Thomas Underwood, David Voelz, Cristo Sanchez, Riley DeColibus, Sarah A. Kovac, Sean Sellers, Doug Gilliam, Patrick Boumier, Thierry Appourchaux, Julien Dejonghe, Jean Pierre Rivet, Steve Markham, Saburo Howard, Lyu Abe, Djamel Mekarnia, Masahiro Ikoma, Hidekazu Hanayama , et al. (3 additional authors not shown)

Abstract: We present three-dimensional (3D) maps of Jupiter's atmospheric circulation at cloud-top level from Doppler-imaging data obtained in the visible domain with JIVE, the second node of the JOVIAL network, which is mounted on the Dunn Solar Telescope at Sunspot, New Mexico. We report on 12 nights of observations between May 4 and May 30, 2018, representing a total of about 80 hours. Firstly, the avera… ▽ More We present three-dimensional (3D) maps of Jupiter's atmospheric circulation at cloud-top level from Doppler-imaging data obtained in the visible domain with JIVE, the second node of the JOVIAL network, which is mounted on the Dunn Solar Telescope at Sunspot, New Mexico. We report on 12 nights of observations between May 4 and May 30, 2018, representing a total of about 80 hours. Firstly, the average zonal wind profile derived from our data is compatible with that derived from cloud-tracking measurements performed on Hubble Space Telescope images obtained in April 2018 from the Outer Planet Atmospheres Legacy (OPAL) program. Secondly, we present the first ever two-dimensional maps of Jupiter's atmospheric circulation from Doppler measurements. The zonal velocity map highlights well-known atmospheric features, such as the equatorial hot spots and the Great Red Spot (GRS). In addition to zonal winds, we derive meridional and vertical velocity fields from the Doppler data. The motions attributed to vertical flows are mainly located at the boundary between the equatorial belts and tropical zones, which could indicate active motion in theses regions. Qualitatively, these results compare well to recent Juno data that have unveiled the three-dimensional structure of Jupiter's wind field. To the contrary, the motions attributed to meridional circulation are very different from what is obtained by cloud tracking, except at the GRS. Because of limitations with data resolution and processing techniques, we acknowledge that our measurement of vertical or meridional flows of Jupiter are still to be confirmed. △ Less

Submitted 16 April, 2024; v1 submitted 28 December, 2023; originally announced December 2023.

Journal ref: Planet. Sci. J. 5 100 (2024)

On the hypothesis of an inverted Z-gradient inside Jupiter

Authors: Saburo Howard, Tristan Guillot, Steve Markham, Ravit Helled, Simon Müller, David J. Stevenson, Jonathan I. Lunine, Yamila Miguel, Nadine Nettelmann

Abstract: Models of Jupiter s interior struggle to agree with measurements of the atmospheric composition. Interior models favour a subsolar or solar abundance of heavy elements Z while atmospheric measurements suggest a supersolar abundance. One potential solution may be the presence of an inverted Z-gradient, namely an inward decrease of Z, which implies a larger heavy element abundance in the atmosphere… ▽ More Models of Jupiter s interior struggle to agree with measurements of the atmospheric composition. Interior models favour a subsolar or solar abundance of heavy elements Z while atmospheric measurements suggest a supersolar abundance. One potential solution may be the presence of an inverted Z-gradient, namely an inward decrease of Z, which implies a larger heavy element abundance in the atmosphere than in the outer envelope. We investigate two scenarios in which the inverted Z gradient is located either where helium rain occurs (Mbar level) or at upper levels (kbar level) where a radiative region could exist. We aim to assess how plausible these scenarios are. We calculate interior and evolution models of Jupiter with such inverted Z-gradient and use constraints on the stability and the formation of an inverted Z-gradient. We find that an inverted Z-gradient at the location of helium rain cannot work as it requires a late accretion and of too much material. We find interior models with an inverted Z-gradient at upper levels, due to a radiative zone preventing downward mixing, that could satisfy the present gravity field of the planet. However, our evolution models suggest that this second scenario might not be in place. An inverted Z-gradient in Jupiter could be stable. Yet, its presence either at the Mbar level or kbar level is rather unlikely. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

Characterizing the Near-infrared Spectra of Flares from TRAPPIST-1 During JWST Transit Spectroscopy Observations

Authors: Ward S. Howard, Adam F. Kowalski, Laura Flagg, Meredith A. MacGregor, Olivia Lim, Michael Radica, Caroline Piaulet, Pierre-Alexis Roy, David Lafrenière, Björn Benneke, Alexander Brown, Néstor Espinoza, René Doyon, Louis-Philippe Coulombe, Doug Johnstone, Nicolas B. Cowan, Ray Jayawardhana, Jake D. Turner, Lisa Dang

Abstract: We present the first analysis of JWST near-infrared spectroscopy of stellar flares from TRAPPIST-1 during transits of rocky exoplanets. Four flares were observed from 0.6--2.8 $μ$m with NIRISS and 0.6--3.5 $μ$m with NIRSpec during transits of TRAPPIST-1b, f, and g. We discover P$α$ and Br$β$ line emission and characterize flare continuum at wavelengths from 1--3.5 $μ$m for the first time. Observed… ▽ More We present the first analysis of JWST near-infrared spectroscopy of stellar flares from TRAPPIST-1 during transits of rocky exoplanets. Four flares were observed from 0.6--2.8 $μ$m with NIRISS and 0.6--3.5 $μ$m with NIRSpec during transits of TRAPPIST-1b, f, and g. We discover P$α$ and Br$β$ line emission and characterize flare continuum at wavelengths from 1--3.5 $μ$m for the first time. Observed lines include H$α$, P$α$-P$ε$, Br$β$, He I $λ$0.7062$μ$m, two Ca II infrared triplet (IRT) lines, and the He I IRT. We observe a reversed Paschen decrement from P$α$-P$γ$ alongside changes in the light curve shapes of these lines. The continuum of all four flares is well-described by blackbody emission with an effective temperature below 5300 K, lower than temperatures typically observed at optical wavelengths. The 0.6--1 $μ$m spectra were convolved with the TESS response, enabling us to measure the flare rate of TRAPPIST-1 in the TESS bandpass. We find flares of 10$^{30}$ erg large enough to impact transit spectra occur at a rate of 3.6$\substack{+2.1 \\ -1.3}$ flare d$^{-1}$, $\sim$10$\times$ higher than previous predictions from K2. We measure the amount of flare contamination at 2 $μ$m for the TRAPPIST-1b and f transits to be 500$\pm$450 and 2100$\pm$400 ppm, respectively. We find up to 80% of flare contamination can be removed, with mitigation most effective from 1.0--2.4 $μ$m. These results suggest transits affected by flares may still be useful for atmospheric characterization efforts. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 29 pages, 17 figures, 3 tables, accepted to The Astrophysical Journal

TOI-199 b: A well-characterized 100-day transiting warm giant planet with TTVs seen from Antarctica

Authors: Melissa J. Hobson, Trifon Trifonov, Thomas Henning, Andrés Jordán, Felipe Rojas, Nestor Espinoza, Rafael Brahm, Jan Eberhardt, Matías I. Jones, Djamel Mekarnia, Diana Kossakowski, Martin Schlecker, Marcelo Tala Pinto, Pascal José Torres Miranda, Lyu Abe, Khalid Barkaoui, Philippe Bendjoya, François Bouchy, Marco Buttu, Ilaria Carleo, Karen A. Collins, Knicole D. Colón, Nicolas Crouzet, Diana Dragomir, Georgina Dransfield , et al. (27 additional authors not shown)

Abstract: We present the spectroscopic confirmation and precise mass measurement of the warm giant planet TOI-199 b. This planet was first identified in TESS photometry and confirmed using ground-based photometry from ASTEP in Antarctica including a full 6.5$\,$h long transit, PEST, Hazelwood, and LCO; space photometry from NEOSSat; and radial velocities (RVs) from FEROS, HARPS, CORALIE, and CHIRON. Orbitin… ▽ More We present the spectroscopic confirmation and precise mass measurement of the warm giant planet TOI-199 b. This planet was first identified in TESS photometry and confirmed using ground-based photometry from ASTEP in Antarctica including a full 6.5$\,$h long transit, PEST, Hazelwood, and LCO; space photometry from NEOSSat; and radial velocities (RVs) from FEROS, HARPS, CORALIE, and CHIRON. Orbiting a late G-type star, TOI-199\,b has a $\mathrm{104.854_{-0.002}^{+0.001} \, d}$ period, a mass of $\mathrm{0.17\pm0.02 \, M_J}$, and a radius of $\mathrm{0.810\pm0.005 \, R_J}$. It is the first warm exo-Saturn with a precisely determined mass and radius. The TESS and ASTEP transits show strong transit timing variations, pointing to the existence of a second planet in the system. The joint analysis of the RVs and TTVs provides a unique solution for the non-transiting companion TOI-199 c, which has a period of $\mathrm{273.69_{-0.22}^{+0.26} \, d}$ and an estimated mass of $\mathrm{0.28_{-0.01}^{+0.02} \, M_J}$. This period places it within the conservative Habitable Zone. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 33 pages, 23 figures. Accepted for publication in AJ

Exoplanet Interior Retrievals: core masses and metallicities from atmospheric abundances

Authors: Sanne Bloot, Yamila Miguel, Michaël Bazot, Saburo Howard

Abstract: The mass and distribution of metals in the interiors of exoplanets are essential for constraining their formation and evolution processes. Nevertheless, with only masses and radii measured, the determination of exoplanet interior structures is degenerate, and so far simplified assumptions have mostly been used to derive planetary metallicities. In this work, we present a method based on a state-of… ▽ More The mass and distribution of metals in the interiors of exoplanets are essential for constraining their formation and evolution processes. Nevertheless, with only masses and radii measured, the determination of exoplanet interior structures is degenerate, and so far simplified assumptions have mostly been used to derive planetary metallicities. In this work, we present a method based on a state-of-the-art interior code, recently used for Jupiter, and a Bayesian framework, to explore the possibility of retrieving the interior structure of exoplanets. We use masses, radii, equilibrium temperatures, and measured atmospheric metallicities to retrieve planetary bulk metallicities and core masses. Following results on the giant planets in the solar system and recent development in planet formation, we implement two interior structure models: one with a homogeneous envelope and one with an inhomogeneous one. Our method is first evaluated using a test planet and then applied to a sample of 37 giant exoplanets with observed atmospheric metallicities from the pre-JWST era. Although neither internal structure model is preferred with the current data, it is possible to obtain information on the interior properties of the planets, such as the core mass, through atmospheric measurements in both cases. We present updated metal mass fractions, in agreement with recent results on giant planets in the solar system. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

TOI-1130: A photodynamical analysis of a hot Jupiter in resonance with an inner low-mass planet

Authors: J. Korth, D. Gandolfi, J. Šubjak, S. Howard, S. Ataiee, K. A. Collins, S. N. Quinn, A. J. Mustill, T. Guillot, N. Lodieu, A. M. S. Smith, M. Esposito, F. Rodler, A. Muresan, L. Abe, S. H. Albrecht, A. Alqasim, K. Barkaoui, P. G. Beck, C. J. Burke, R. P. Butler, D. M. Conti, K. I. Collins, J. D. Crane, F. Dai , et al. (37 additional authors not shown)

Abstract: The TOI-1130 is a known planetary system around a K-dwarf consisting of a gas giant planet, TOI-1130 c, on an 8.4-day orbit, accompanied by an inner Neptune-sized planet, TOI-1130 b, with an orbital period of 4.1 days. We collected precise radial velocity (RV) measurements of TOI-1130 with the HARPS and PFS spectrographs as part of our ongoing RV follow-up program. We perform a photodynamical mode… ▽ More The TOI-1130 is a known planetary system around a K-dwarf consisting of a gas giant planet, TOI-1130 c, on an 8.4-day orbit, accompanied by an inner Neptune-sized planet, TOI-1130 b, with an orbital period of 4.1 days. We collected precise radial velocity (RV) measurements of TOI-1130 with the HARPS and PFS spectrographs as part of our ongoing RV follow-up program. We perform a photodynamical modeling of the HARPS and PFS RVs, and transit photometry from the Transiting Exoplanet Survey Satellite (TESS) and the TESS Follow-up Observing Program. We determine the planet masses and radii of TOI-1130 b and TOI-1130 c to be Mb = 19.28 $\pm$ 0.97 M$_\oplus$ and Rb = 3.56 $\pm$ 0.13 R$_\oplus$, and Mc = 325.59 $\pm$ 5.59 M$_\oplus$ and Rc = 13.32+1.55-1.41 R$_\oplus$, respectively. We spectroscopically confirm TOI-1130 b that was previously only validated. We find that the two planets orbit with small eccentricities in a 2:1 resonant configuration. This is the first known system with a hot Jupiter and an inner lower mass planet locked in a mean-motion resonance. TOI-1130 belongs to the small yet increasing population of hot Jupiters with an inner low-mass planet that challenges the pathway for hot Jupiter formation. We also detect a linear RV trend possibly due to the presence of an outer massive companion. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 19 pages, Accepted to A&A

Journal ref: A&A 675, A115 (2023)

The apparent absence of forward scattering in the HD 53143 debris disk

Authors: Christopher C. Stark, Bin Ren, Meredith A. MacGregor, Ward S. Howard, Spencer A. Hurt, Alycia J. Weinberger, Glenn Schneider, Elodie Choquet

Abstract: HD 53143 is a mature Sun-like star and host to a broad disk of dusty debris, including a cold outer ring of planetesimals near 90 AU. Unlike most other inclined debris disks imaged at visible wavelengths, the cold disk around HD 53143 appears as disconnected "arcs" of material, with no forward scattering side detected to date. We present new, deeper Hubble Space Telescope (HST) Space Telescope Ima… ▽ More HD 53143 is a mature Sun-like star and host to a broad disk of dusty debris, including a cold outer ring of planetesimals near 90 AU. Unlike most other inclined debris disks imaged at visible wavelengths, the cold disk around HD 53143 appears as disconnected "arcs" of material, with no forward scattering side detected to date. We present new, deeper Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) coronagraphic observations of the HD 53143 debris disk and show that the forward scattering side of the disk remains undetected. By fitting our KLIP-reduced observations via forward modeling with an optically thin disk model, we show that fitting the visible wavelength images with an azimuthally symmetric disk with unconstrained orientation results in an unphysical edge-on orientation that is at odds with recent ALMA observations, while constraining the orientation to that observed by ALMA results in nearly isotropically scattering dust. We show that the HD53143 host star exhibits significant stellar variations due to spot rotation and revisit age estimates for this system. △ Less

Submitted 14 April, 2023; originally announced April 2023.

Comments: 13 pages, 8 figures, 2 tables Accepted for publication in ApJ

Journal ref: ApJ 945 131 (2023)

Three Saturn-mass planets transiting F-type stars revealed with TESS and HARPS

Authors: Angelica Psaridi, François Bouchy, Monika Lendl, Babatunde Akinsanmi, Keivan G. Stassun, Barry Smalley, David J. Armstrong, Saburo Howard, Solène Ulmer-Moll, Nolan Grieves, Khalid Barkaoui, Joseph E. Rodriguez, Edward M. Bryant, Olga Suárez, Tristan Guillot, Phil Evans, Omar Attia, Robert A. Wittenmyer, Samuel W. Yee, Karen A. Collins, George Zhou, Franck Galland, Léna Parc, Stéphane Udry, Pedro Figueira , et al. (40 additional authors not shown)

Abstract: While the sample of confirmed exoplanets continues to increase, the population of transiting exoplanets around early-type stars is still limited. These planets allow us to investigate the planet properties and formation pathways over a wide range of stellar masses and study the impact of high irradiation on hot Jupiters orbiting such stars. We report the discovery of TOI-615b, TOI-622b, and TOI-26… ▽ More While the sample of confirmed exoplanets continues to increase, the population of transiting exoplanets around early-type stars is still limited. These planets allow us to investigate the planet properties and formation pathways over a wide range of stellar masses and study the impact of high irradiation on hot Jupiters orbiting such stars. We report the discovery of TOI-615b, TOI-622b, and TOI-2641b, three Saturn-mass planets transiting main sequence, F-type stars. The planets were identified by the Transiting Exoplanet Survey Satellite (TESS) and confirmed with complementary ground-based and radial velocity observations. TOI-615b is a highly irradiated ($\sim$1277 $F_{\oplus}$) and bloated Saturn-mass planet (1.69$^{+0.05}_{-0.06}$$R_{Jup}$ and 0.43$^{+0.09}_{-0.08}$$M_{Jup}$) in a 4.66 day orbit transiting a 6850 K star. TOI-622b has a radius of 0.82$^{+0.03}_{-0.03}$$R_{Jup}$ and a mass of 0.30$^{+0.07}_{-0.08}$~$M_{Jup}$ in a 6.40 day orbit. Despite its high insolation flux ($\sim$600 $F_{\oplus}$), TOI-622b does not show any evidence of radius inflation. TOI-2641b is a 0.39$^{+0.02}_{-0.04}$$M_{Jup}$ planet in a 4.88 day orbit with a grazing transit (b = 1.04$^{+0.05}_{-0.06 }$) that results in a poorly constrained radius of 1.61$^{+0.46}_{-0.64}$$R_{Jup}$. Additionally, TOI-615b is considered attractive for atmospheric studies via transmission spectroscopy with ground-based spectrographs and $\textit{JWST}$. Future atmospheric and spin-orbit alignment observations are essential since they can provide information on the atmospheric composition, formation and migration of exoplanets across various stellar types. △ Less

Submitted 11 May, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 16 pages, 17 figures, submitted to A&A

Journal ref: A&A 675, A39 (2023)

The Evryscope Fast Transient Engine: Real-Time Detection for Rapidly Evolving Transients

Authors: Hank Corbett, Jonathan Carney, Ramses Gonzalez, Octavi Fors, Nathan Galliher, Amy Glazier, Ward S. Howard, Nicholas M. Law, Robert Quimby, Jeffrey K. Ratzloff, Alan Vasquez Soto

Abstract: Astrophysical transients with rapid development on sub-hour timescales are intrinsically rare. Due to their short durations, events like stellar superflares, optical flashes from gamma-ray bursts, and shock breakouts from young supernovae are difficult to identify on timescales that enable spectroscopic followup. This paper presents the Evryscope Fast Transient Engine (EFTE), a new data reduction… ▽ More Astrophysical transients with rapid development on sub-hour timescales are intrinsically rare. Due to their short durations, events like stellar superflares, optical flashes from gamma-ray bursts, and shock breakouts from young supernovae are difficult to identify on timescales that enable spectroscopic followup. This paper presents the Evryscope Fast Transient Engine (EFTE), a new data reduction pipeline designed to provide low-latency transient alerts from the Evryscopes, a North-South pair of ultra-wide-field telescopes with an instantaneous footprint covering 38% of the entire sky, and tools for building long-term light curves from Evryscope data. EFTE leverages the optical stability of the Evryscopes by using a simple direct image subtraction routine suited to continuously monitoring the transient sky at minute cadence. Candidates are produced within the base Evryscope two-minute cadence for 98.5% of images, and internally filtered using VetNet, a convolutional neural network real-bogus classifier. EFTE provides an extensible, robust architecture for transient surveys probing similar timescales, and serves as the software testbed for the real-time analysis pipelines and public data distribution systems for the Argus Array, a next generation all-sky observatory with a data rate 62x higher than Evryscope. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 28 pages, 18 figures. in press, ApJS

Jupiter's interior from Juno: Equation-of-state uncertainties and dilute core extent

Authors: S. Howard, T. Guillot, M. Bazot, Y. Miguel, D. J. Stevenson, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, N. Nettelmann, B. Idini, S. Bolton

Abstract: The Juno mission has provided measurements of Jupiter s gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. In this study, we investigated the differences betwe… ▽ More The Juno mission has provided measurements of Jupiter s gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. In this study, we investigated the differences between the state-of-the-art equations of state and their impact on the properties of interior models. Accounting for uncertainty on the hydrogen and helium equation of state, we assessed the span of the interior features of Jupiter. We carried out an extensive exploration of the parameter space and studied a wide range of interior models using Markov chain Monte Carlo (MCMC) simulations. To consider the uncertainty on the equation of state, we allowed for modifications of the equation of state in our calculations. Our models harbour a dilute core and indicate that Jupiter s internal entropy is higher than what is usually assumed from the Galileo probe measurements. We obtain solutions with extended dilute cores, but contrary to other recent interior models of Jupiter, we also obtain models with small dilute cores. The dilute cores in such solutions extend to 20% of Jupiter s mass, leading to better agreement with formation evolution models. We conclude that the equations of state used in Jupiter models have a crucial effect on the inferred structure and composition. Further explorations of the behaviour of hydrogen helium mixtures at the pressure and temperature conditions in Jupiter will help to constrain the interior of the planet, and therefore its origin. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Journal ref: A&A 672, A33 (2023)

Accounting for non-ideal mixing effects in the hydrogen-helium equation of state

Authors: Saburo Howard, Tristan Guillot

Abstract: The equation of state for hydrogen and helium is fundamental for studying stars and giant planets. It has been shown that because of interactions at atomic and molecular levels, the behaviour of a mixture of hydrogen and helium cannot be accurately represented by considering these elements separately. This paper aims at providing a simple method to account for interactions between hydrogen and hel… ▽ More The equation of state for hydrogen and helium is fundamental for studying stars and giant planets. It has been shown that because of interactions at atomic and molecular levels, the behaviour of a mixture of hydrogen and helium cannot be accurately represented by considering these elements separately. This paper aims at providing a simple method to account for interactions between hydrogen and helium in interior and evolution models of giant planets. Using on the one hand ab initio simulations that involve a system of interacting hydrogen and helium particles and pure equations of state for hydrogen and helium on the other, we derived the contributions in density and entropy of the interactions between hydrogen and helium particles. We show that relative variations of up to 15% in density and entropy arise when non-ideal mixing is accounted for. These non-ideal mixing effects must be considered in interior models of giant planets based on accurate gravity field measurements, particularly in the context of variations in the helium-to-hydrogen ratio. They also affect the mass-radius relation of exoplanets. We provide a table that contains the volume and entropy of mixing as a function of pressure and temperature. This table is to be combined with pure hydrogen and pure helium equations of state to obtain an equation of state that self-consistently includes mixing effects for any hydrogen and helium mixing ratio and may be used to model the interior structure and evolution of giant planets to brown dwarfs. Non-linear mixing must be included in accurate calculations of the equations of state of hydrogen and helium. Uncertainties on the equation of state still exist, however. Ab initio calculations of the behaviour of the hydrogen-helium mixture in the megabar regime for various compositions should be performed in order to gain accuracy. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Journal ref: A&A 672, L1 (2023)

TOI-2525 b and c: A pair of massive warm giant planets with a strong transit timing variations revealed by TESS

Authors: Trifon Trifonov, Rafael Brahm, Andres Jordan, Christian Hartogh, Thomas Henning, Melissa J. Hobson, Martin Schlecker, Saburo Howard, Finja Reichardt, Nestor Espinoza, Man Hoi Lee, David Nesvorny, Felipe I. Rojas, Khalid Barkaoui, Diana Kossakowski, Gavin Boyle, Stefan Dreizler, Martin Kuerster, Rene Heller, Tristan Guillot, Amaury H. M. J. Triaud, Lyu Abe, Abdelkrim Agabi, Philippe Bendjoya, Nicolas Crouzet , et al. (22 additional authors not shown)

Abstract: TOI-2525 is a K-type star with an estimated mass of M = 0.849$_{-0.033}^{+0.024}$ M$_\odot$ and radius of R = 0.785$_{-0.007}^{+0.007}$ R$_\odot$ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semi-amplitude of a $\sim$6 hours. We performed TTV… ▽ More TOI-2525 is a K-type star with an estimated mass of M = 0.849$_{-0.033}^{+0.024}$ M$_\odot$ and radius of R = 0.785$_{-0.007}^{+0.007}$ R$_\odot$ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semi-amplitude of a $\sim$6 hours. We performed TTV dynamical, and photo-dynamical light curve analysis of the TESS data, combined with radial velocity (RV) measurements from FEROS and PFS, and we confirmed the planetary nature of these companions. The TOI-2525 system consists of a transiting pair of planets comparable to Neptune and Jupiter with estimated dynamical masses of $m_{\rm b}$ = 0.088$_{-0.004}^{+0.005}$ M$_{\rm Jup.}$, and $m_{\rm c}$ = 0.709$_{-0.033}^{+0.034}$ M$_{\rm Jup.}$, radius of $r_b$ = 0.88$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$ and $r_c$ = 0.98$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$, and with orbital periods of $P_{\rm b}$ = 23.288$_{-0.002}^{+0.001}$ days and $P_{\rm c}$ = 49.260$_{-0.001}^{+0.001}$ days for the inner and the outer planet, respectively. The period ratio is close to the 2:1 period commensurability, but the dynamical simulations of the system suggest that it is outside the mean motion resonance (MMR) dynamical configuration. TOI-2525 b is among the lowest density Neptune-mass planets known to date, with an estimated median density of $ρ_{\rm b}$ = 0.174$_{-0.015}^{+0.016}$ g\,cm$^{-3}$. The TOI-2525 system is very similar to the other K-dwarf systems discovered by TESS, TOI-2202 and TOI-216, which are composed of almost identical K-dwarf primary and two warm giant planets near the 2:1 MMR. △ Less

Submitted 11 February, 2023; originally announced February 2023.

Comments: Accepted for publication in AJ

HIP 33609 b: An Eccentric Brown Dwarf Transiting a V=7.3 Rapidly Rotating B-Star

Authors: Noah Vowell, Joseph E. Rodriguez, Samuel N. Quinn, George Zhou, Andrew Vanderburg, Andrew W. Mann, Matthew J. Hooton, Keivan G. Stassun, Saburo Howard, Allyson Bieryla, David W. Latham, Steve B. Howell, Tristan Guillot, Carl Ziegler, Karen A. Collins, Theron W. Carmichael, Jon M. Jenkins, Avi Shporer, Lyu ABE, Philippe Bendjoya, Jonathan L. Bush, Marco Buttu, Kevin I. Collins, Jason D. Eastman, Matthew J. Fields , et al. (19 additional authors not shown)

Abstract: We present the discovery and characterization of HIP 33609 b, a transiting warm brown dwarf orbiting a late B star, discovered by NASA's Transiting Exoplanet Survey Satellite TESS as TOI-588 b. HIP 33609 b is a large (R$_{b}$ = 1.580$_{-0.070}^{+0.074}$ R$_{J}$) brown dwarf on a highly eccentric (e = 0.560$_{-0.031}^{+0.029}$) orbit with a 39-day period. The host star is a bright (V = 7.3 mag), T… ▽ More We present the discovery and characterization of HIP 33609 b, a transiting warm brown dwarf orbiting a late B star, discovered by NASA's Transiting Exoplanet Survey Satellite TESS as TOI-588 b. HIP 33609 b is a large (R$_{b}$ = 1.580$_{-0.070}^{+0.074}$ R$_{J}$) brown dwarf on a highly eccentric (e = 0.560$_{-0.031}^{+0.029}$) orbit with a 39-day period. The host star is a bright (V = 7.3 mag), T$_{eff}$ = 10,400$_{-660}^{+800}$ K star with a mass of M$_{*}$ = 2.383$_{-0.095}^{+0.10}$ M$_{\odot}$ and radius of R$_{*}$ = 1.863$_{-0.082}^{+0.087}$ R$_{\odot}$, making it the hottest transiting brown dwarf host star discovered to date. We obtained radial velocity measurements from the CHIRON spectrograph confirming the companion's mass of M$_{b}$ = 68.0$_{-7.1}^{+7.4}$ M$_{J}$ as well as the host star's rotation rate ($vsini_{*} = 55.6 \pm 1.8$ km/s). We also present the discovery of a new comoving group of stars, designated as MELANGE-6, and determine that HIP 33609 is a member. We use a combination of rotation periods and isochrone models fit to the cluster members to estimate an age of 150 $\pm$ 25 Myr. With a measured mass, radius, and age, HIP 33609 b becomes a benchmark for substellar evolutionary models. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 16 pages, 5 figures, 4 tables, Submitted to AAS Journals

TESS spots a mini-neptune interior to a hot saturn in the TOI-2000 system

Authors: Lizhou Sha, Andrew M. Vanderburg, Chelsea X. Huang, David J. Armstrong, Rafael Brahm, Steven Giacalone, Mackenna L. Wood, Karen A. Collins, Louise D. Nielsen, Melissa J. Hobson, Carl Ziegler, Steve B. Howell, Pascal Torres-Miranda, Andrew W. Mann, George Zhou, Elisa Delgado-Mena, Felipe I. Rojas, Lyu Abe, Trifon Trifonov, Vardan Adibekyan, Sérgio G. Sousa, Sergio B. Fajardo-Acosta, Tristan Guillot, Saburo Howard, Colin Littlefield , et al. (30 additional authors not shown)

Abstract: Hot jupiters (P < 10 d, M > 60 $\mathrm{M}_\oplus$) are almost always found alone around their stars, but four out of hundreds known have inner companion planets. These rare companions allow us to constrain the hot jupiter's formation history by ruling out high-eccentricity tidal migration. Less is known about inner companions to hot Saturn-mass planets. We report here the discovery of the TOI-200… ▽ More Hot jupiters (P < 10 d, M > 60 $\mathrm{M}_\oplus$) are almost always found alone around their stars, but four out of hundreds known have inner companion planets. These rare companions allow us to constrain the hot jupiter's formation history by ruling out high-eccentricity tidal migration. Less is known about inner companions to hot Saturn-mass planets. We report here the discovery of the TOI-2000 system, which features a hot Saturn-mass planet with a smaller inner companion. The mini-neptune TOI-2000 b ($2.70 \pm 0.15 \,\mathrm{R}_\oplus$, $11.0 \pm 2.4 \,\mathrm{M}_\oplus$) is in a 3.10-day orbit, and the hot saturn TOI-2000 c ($8.14^{+0.31}_{-0.30} \,\mathrm{R}_\oplus$, $81.7^{+4.7}_{-4.6} \,\mathrm{M}_\oplus$) is in a 9.13-day orbit. Both planets transit their host star TOI-2000 (TIC 371188886, V = 10.98, TESS magnitude = 10.36), a metal-rich ([Fe/H] = $0.439^{+0.041}_{-0.043}$) G dwarf 174 pc away. TESS observed the two planets in sectors 9-11 and 36-38, and we followed up with ground-based photometry, spectroscopy, and speckle imaging. Radial velocities from CHIRON, FEROS, and HARPS allowed us to confirm both planets by direct mass measurement. In addition, we demonstrate constraining planetary and stellar parameters with MIST stellar evolutionary tracks through Hamiltonian Monte Carlo under the PyMC framework, achieving higher sampling efficiency and shorter run time compared to traditional Markov chain Monte Carlo. Having the brightest host star in the V band among similar systems, TOI-2000 b and c are superb candidates for atmospheric characterization by the JWST, which can potentially distinguish whether they formed together or TOI-2000 c swept along material during migration to form TOI-2000 b. △ Less

Submitted 31 May, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: v3 adds RV frequency analysis; 25 pages, 11 figures, 14 tables; revision submitted to MNRAS; machine-readable tables available as ancillary files; posterior samples available from Zenodo at https://doi.org/10.5281/zenodo.7683293 and source code at https://doi.org/10.5281/zenodo.7988268

The Mouse that Squeaked: A small flare from Proxima Cen observed in the millimeter, optical, and soft X-ray with Chandra and ALMA

Authors: Ward S. Howard, Meredith A. MacGregor, Rachel Osten, Jan Forbrich, Steven R. Cranmer, Isaiah Tristan, Alycia J. Weinberger, Allison Youngblood, Thomas Barclay, R. O. Parke Loyd, Evgenya L. Shkolnik, Andrew Zic, David J. Wilner

Abstract: We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that… ▽ More We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that had gone undetected until recently, making it difficult to interpret these events within the current multi-wavelength picture of the flaring process. The May 6 event is the smallest stellar millimeter flare detected to date. We compare the relationship between the soft X-ray and millimeter emission to that observed in solar flares. The X-ray and optical flare energies of 10$^{30.3\pm0.2}$ and 10$^{28.9\pm0.1}$ erg, respectively, the coronal temperature of T=11.0$\pm$2.1 MK, and the emission measure of 9.5$\pm$2.2 X 10$^{49}$ cm$^{-3}$ are consistent with M-X class solar flares. We find the soft X-ray and millimeter emission during quiescence are consistent with the Gudel-Benz Relation, but not during the flare. The millimeter luminosity is >100X higher than that of an equivalent X1 solar flare and lasts only seconds instead of minutes as seen for solar flares. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 12 pages, 5 figures. Accepted to ApJ

ALMA Images the Eccentric HD 53143 Debris Disk

Authors: Meredith A. MacGregor, Spencer A. Hurt, Christopher C. Stark, Ward S. Howard, Alycia J. Weinberger, Bin Ren, Glenn Schneider, Elodie Choquet, Dmitri Mawet

Abstract: We present ALMA 1.3 mm observations of the HD~53143 debris disk - the first infrared or millimeter image produced of this ~1 Gyr-old solar-analogue. Previous HST STIS coronagraphic imaging did not detect flux along the minor axis of the disk which could suggest a face-on geometry with two 'clumps' of dust. These ALMA observations reveal a disk with a strikingly different structure. In order to fit… ▽ More We present ALMA 1.3 mm observations of the HD~53143 debris disk - the first infrared or millimeter image produced of this ~1 Gyr-old solar-analogue. Previous HST STIS coronagraphic imaging did not detect flux along the minor axis of the disk which could suggest a face-on geometry with two 'clumps' of dust. These ALMA observations reveal a disk with a strikingly different structure. In order to fit models to the millimeter visibilities and constrain the uncertainties on the disk parameters, we adopt an MCMC approach. This is the most eccentric debris disk observed to date with a forced eccentricity of $0.21\pm0.02$, nearly twice that of the Fomalhaut debris disk, and also displays apocenter glow. Although this eccentric model fits the outer debris disk well, there are significant interior residuals remaining that may suggest a possible edge-on inner disk, which remains unresolved in these observations. Combined with the observed structure difference between HST and ALMA, these results suggest a potential previous scattering event or dynamical instability in this system. We also note that the stellar flux changes considerably over the course of our observations, suggesting flaring at millimeter wavelengths. Using simultaneous TESS observations, we determine the stellar rotation period to be $9.6\pm0.1$ days. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 9 pages, 3 figures, 2 tables, published in ApJ Letters

The Flaring TESS Objects of Interest: Flare Rates for all Two Minute Cadence TESS Planet Candidates

Authors: Ward S. Howard

Abstract: Although more than 5000 TESS Objects of Interest have been cataloged, no comprehensive survey of the flare rates of their host stars exists. We perform the first flare survey of all 2250 non-retired TOIs with 2 min cadence light curves to measure or place upper limits on their flare rates. We find 93 candidates orbit flare stars and measure their flare frequency distributions. Across the sample, T… ▽ More Although more than 5000 TESS Objects of Interest have been cataloged, no comprehensive survey of the flare rates of their host stars exists. We perform the first flare survey of all 2250 non-retired TOIs with 2 min cadence light curves to measure or place upper limits on their flare rates. We find 93 candidates orbit flare stars and measure their flare frequency distributions. Across the sample, TOIs of <1.5R_Earth orbit flare stars more frequently than do TOIs of 1.5<R<2.75R_Earth, 2.75<R<4R_Earth, or R<4R_Earth. We sort all TOI host stars by their flare rate/upper limit, stellar mass, and distance to create a flare ranking metric (FRM) to determine suitability for follow-up. The FRM of each TOI is then checked against the expected signal-to-noise of atmospheric features in transmission spectroscopy to locate the most promising targets. We find 1/4 of terrestrial M-dwarf planets amenable to transmission spectroscopy orbit flare stars. However, none of the M-dwarf hosts to terrestrial planets are currently flaring at sufficient levels for >99.9% atmospheric ozone depletion. We give the first upper limits on the flare rate of the host star to TOI 700 d and explore the flare rates incident on young planets such as DS Tuc Ab. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 6 pages, 3 figures, 1 table. Accepted for publication in MNRAS Letters. Ancillary machine-readable table included

Jupiter's inhomogeneous envelope

Authors: Y. Miguel, M. Bazot, T. Guillot, S. Howard, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, S. K. Atreya, J. E. P. Connerney, D. Durante, L. Kulowski, J. I. Lunine, D. Stevenson, S. Bolton

Abstract: While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this… ▽ More While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature 5-15 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 Earth masses, independently on the structure model (with or without dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogenous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our solar system and beyond. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 662, A18 (2022)

No Such Thing as a Simple Flare: Substructure and QPPs Observed in a Statistical Sample of 20 Second Cadence TESS Flares

Authors: Ward S. Howard, Meredith A. MacGregor

Abstract: A 20 second cadence TESS monitoring campaign of 226 low-mass flare stars during Cycle 3 recorded 3792 stellar flares of >10^32 erg. We explore the time-resolved emission and substructure in 440 of the largest flares observed at high S/N, 97% of which released energies of >10^33 erg. We discover degeneracy present at 2 minute cadence between sharply-peaked and weakly-peaked flares is common, althou… ▽ More A 20 second cadence TESS monitoring campaign of 226 low-mass flare stars during Cycle 3 recorded 3792 stellar flares of >10^32 erg. We explore the time-resolved emission and substructure in 440 of the largest flares observed at high S/N, 97% of which released energies of >10^33 erg. We discover degeneracy present at 2 minute cadence between sharply-peaked and weakly-peaked flares is common, although 20 second cadence breaks these degeneracies. We better resolve the rise phases and find 46% of large flares exhibit substructure during the rise phase. We observe 49 candidate quasi-periodic pulsations (QPP) and confirm 17 at 3+ sigma. Most of our QPPs have periods less than 10 minutes, suggesting short period optical QPPs are common. We find QPPs in both the rise and decay phases of flares, including a rise-phase QPP in a large flare from Proxima Cen. We confirm the Davenport et al. (2014) template provides a good fit to most classical flares observed at high cadence, although 9% favor Gaussian peaks instead. We characterize the properties of complex flares, finding 17% of complex flares exhibit "peak-bump" morphologies composed of a large, highly impulsive peak followed by a second more gradual Gaussian peak. We also estimate the UVC surface fluences of temperate planets at flare peak and find 1/3 of 10^34 erg flares reach the D90 dose of D. Radiodurans in just 20 seconds in the absence of an atmosphere. △ Less

Submitted 8 December, 2021; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 16 pages, 13 figures, 2 tables (machine-readable ancillary files included). Accepted for publication in ApJ after minor revisions

EvryFlare IV: Detection of periodicity in flare occurrence from cool stars with TESS

Authors: Ward S. Howard, Nicholas M. Law

Abstract: Phased flaring, or the periodic occurrence of stellar flares, may probe electromagnetic star-planet interaction (SPI), binary interaction, or magnetic conditions in spots. For the first time, we explore flare periodograms for a large sample of flare stars to identify periodicity due to magnetic interactions with orbiting companions, magnetic reservoirs, or rotational phase. Previous large surveys… ▽ More Phased flaring, or the periodic occurrence of stellar flares, may probe electromagnetic star-planet interaction (SPI), binary interaction, or magnetic conditions in spots. For the first time, we explore flare periodograms for a large sample of flare stars to identify periodicity due to magnetic interactions with orbiting companions, magnetic reservoirs, or rotational phase. Previous large surveys have explored periodicity at the stellar rotation period, but we do not assume periods must correspond with rotation in this work. Two min TESS light curves of 284 cool stars are searched for periods from 1-10 d using two newly-developed periodograms. Because flares are discrete events in noisy and incomplete data, typical periodograms are not well-suited to detect phased flaring. We construct and test a new Bayesian likelihood periodogram and a modified Lomb-Scargle periodogram. We find 6 candidates with a false-alarm probability below 1%. Three targets are >3-sigma detections of flare periodicity; the others are plausible candidates which cannot be individually confirmed. Periods range from 1.35 to 6.7 d and some, but not all, correlate with the stellar rotation period or its 1/2 alias. Periodicity from 2 targets may persist from TESS Cycle 1 into Cycle 3. The periodicity does not appear to persist for the others. Long-term changes in periodicity may result from the spot evolution observed from each candidate, which suggests magnetic conditions play an important role in sustaining periodicity. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: Accepted to the Astrophysical Journal. 14 pages, 9 figures

Low-Cost Access to the Deep, High-Cadence Sky: the Argus Optical Array

Authors: Nicholas M. Law, Hank Corbett, Nathan W. Galliher, Ramses Gonzalez, Alan Vasquez, Glenn Walters, Lawrence Machia, Jeff Ratzloff, Kendall Ackley, Chris Bizon, Christopher Clemens, Steven Cox, Steven Eikenberry, Ward S. Howard, Amy Glazier, Andrew W. Mann, Robert Quimby, Daniel Reichart, David Trilling

Abstract: New mass-produced, wide-field, small-aperture telescopes have the potential to revolutionize ground-based astronomy by greatly reducing the cost of collecting area. In this paper, we introduce a new class of large telescope based on these advances: an all-sky, arcsecond-resolution, 1000-telescope array which builds a simultaneously high-cadence and deep survey by observing the entire sky all night… ▽ More New mass-produced, wide-field, small-aperture telescopes have the potential to revolutionize ground-based astronomy by greatly reducing the cost of collecting area. In this paper, we introduce a new class of large telescope based on these advances: an all-sky, arcsecond-resolution, 1000-telescope array which builds a simultaneously high-cadence and deep survey by observing the entire sky all night. As a concrete example, we describe the Argus Array, a 5m-class telescope with an all-sky field of view and the ability to reach extremely high cadences using low-noise CMOS detectors. Each 55 GPix Argus exposure covers 20% of the entire sky to g=19.6 each minute and g=21.9 each hour; a high-speed mode will allow sub-second survey cadences for short times. Deep coadds will reach g=23.6 every five nights over 47% of the sky; a larger-aperture array telescope, with an étendue close to the Rubin Observatory, could reach g=24.3 in five nights. These arrays can build two-color, million-epoch movies of the sky, enabling sensitive and rapid searches for high-speed transients, fast-radio-burst counterparts, gravitational-wave counterparts, exoplanet microlensing events, occultations by distant solar system bodies, and myriad other phenomena. An array of O(1,000) telescopes, however, would be one of the most complex astronomical instruments yet built. Standard arrays with hundreds of tracking mounts entail thousands of moving parts and exposed optics, and maintenance costs would rapidly outpace the mass-produced-hardware cost savings compared to a monolithic large telescope. We discuss how to greatly reduce operations costs by placing all optics in a thermally controlled, sealed dome with a single moving part. Coupled with careful software scope control and use of existing pipelines, we show that the Argus Array could become the deepest and fastest Northern sky survey, with total costs below $20M. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: 17 pages, 5 figures, 2 tables

Rotation periods of TESS Objects of Interest from the Magellan-TESS Survey with multiband photometry from Evryscope and TESS

Authors: Ward S. Howard, Johanna Teske, Hank Corbett, Nicholas M. Law, Sharon Xuesong Wang, Jeffrey K. Ratzloff, Nathan W. Galliher, Ramses Gonzalez, Alan Vasquez Soto, Amy L. Glazier, Joshua Haislip

Abstract: Stellar RV jitter due to surface activity may bias the RV semi-amplitude and mass of rocky planets. The amplitude of the jitter may be estimated from the uncertainty in the rotation period, allowing the mass to be more accurately obtained. We find candidate rotation periods for 17 out of 35 TESS Objects of Interest (TOI) hosting <3 R_Earth planets as part of the Magellan-TESS Survey, which is the… ▽ More Stellar RV jitter due to surface activity may bias the RV semi-amplitude and mass of rocky planets. The amplitude of the jitter may be estimated from the uncertainty in the rotation period, allowing the mass to be more accurately obtained. We find candidate rotation periods for 17 out of 35 TESS Objects of Interest (TOI) hosting <3 R_Earth planets as part of the Magellan-TESS Survey, which is the first-ever statistically robust study of exoplanet masses and radii across the photo-evaporation gap. Seven periods are 3+ sigma detections, two are 1.5+ sigma, and 8 show plausible variability but the periods remain unconfirmed. The other 18 TOIs are non-detections. Candidate rotators include the host stars of the confirmed planets L 168-9 b, the HD 21749 system, LTT 1445 A b, TOI 1062 b, and the L 98-59 system. 13 candidates have no counterpart in the 1000 TOI rotation catalog of Canto Martins et al. (2020). We find periods for G3-M3 dwarfs using combined light curves from TESS and the Evryscope all-sky array of small telescopes, sometimes with longer periods than would be possible with TESS alone. Secure periods range from 1.4 to 26 d with Evryscope-measured photometric amplitudes as small as 2.1 mmag in g'. We also apply Monte Carlo sampling and a Gaussian Process stellar activity model from the code exoplanet to the TESS light curves of 6 TOIs to confirm the Evryscope periods. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted to the Astronomical Journal. 22 pages, 11 figures. Machine-readable table included

TIC 172900988: A Transiting Circumbinary Planet Detected in One Sector of TESS Data

Authors: Veselin B. Kostov, Brian P. Powell, Jerome A. Orosz, William F. Welsh, William Cochran, Karen A. Collins, Michael Endl, Coel Hellier, David W. Latham, Phillip MacQueen, Joshua Pepper, Billy Quarles, Lalitha Sairam, Guillermo Torres, Robert F. Wilson, Serge Bergeron, Pat Boyce, Allyson Bieryla, Robert Buchheim, Caleb Ben Christiansen, David R. Ciardi, Kevin I. Collins, Dennis M. Conti, Scott Dixon, Pere Guerra , et al. (64 additional authors not shown)

Abstract: We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a… ▽ More We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a prominent apsidal motion of the binary orbit, caused by the dynamical interactions between the binary and the planet. A comprehensive photodynamical analysis of the TESS, archival and follow-up data yields stellar masses and radii of M1 = 1.2384 +/- 0.0007 MSun and R1 = 1.3827 +/- 0.0016 RSun for the primary and M2 = 1.2019 +/- 0.0007 MSun and R2 = 1.3124 +/- 0.0012 RSun for the secondary. The radius of the planet is R3 = 11.25 +/- 0.44 REarth (1.004 +/- 0.039 RJup). The planet's mass and orbital properties are not uniquely determined - there are six solutions with nearly equal likelihood. Specifically, we find that the planet's mass is in the range of 824 < M3 < 981 MEarth (2.65 < M3 < 3.09 MJup), its orbital period could be 188.8, 190.4, 194.0, 199.0, 200.4, or 204.1 days, and the eccentricity is between 0.02 and 0.09. At a V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g. Rossiter-McLaughlin effect and transit spectroscopy. △ Less

Submitted 27 August, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 57 pages, 30 figures, 25 tables; Accepted AJ

Discovery of an Extremely Short Duration Flare from Proxima Centauri Using Millimeter through FUV Observations

Authors: Meredith A. MacGregor, Alycia J. Weinberger, R. O. Parke Loyd, Evgenya Shkolnik, Thomas Barclay, Ward S. Howard, Andrew Zic, Rachel A. Osten, Steven R. Cranmer, Adam F. Kowalski, Emil Lenc, Allison Youngblood, Anna Estes, David J. Wilner, Jan Forbrich, Anna Hughes, Nicholas M. Law, Tara Murphy, Aaron Boley, Jaymie Matthews

Abstract: We present the discovery of an extreme flaring event from Proxima Cen by ASKAP, ALMA, HST, TESS, and the du Pont Telescope that occurred on 2019 May 1. In the millimeter and FUV, this flare is the brightest ever detected, brightening by a factor of >1000 and >14000 as seen by ALMA and HST, respectively. The millimeter and FUV continuum emission trace each other closely during the flare, suggesting… ▽ More We present the discovery of an extreme flaring event from Proxima Cen by ASKAP, ALMA, HST, TESS, and the du Pont Telescope that occurred on 2019 May 1. In the millimeter and FUV, this flare is the brightest ever detected, brightening by a factor of >1000 and >14000 as seen by ALMA and HST, respectively. The millimeter and FUV continuum emission trace each other closely during the flare, suggesting that millimeter emission could serve as a proxy for FUV emission from stellar flares and become a powerful new tool to constrain the high-energy radiation environment of exoplanets. Surprisingly, optical emission associated with the event peaks at a much lower level with a time delay. The initial burst has an extremely short duration, lasting for <10 sec. Taken together with the growing sample of millimeter M dwarf flares, this event suggests that millimeter emission is actually common during stellar flares and often originates from short burst-like events. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 11 pages, 4 figures, 1 appendix, published in ApJ Letters

Revisiting the HD 21749 Planetary System with Stellar Activity Modeling

Authors: Tianjun Gan, Sharon Xuesong Wang, Johanna K. Teske, Shude Mao, Ward S. Howard, Nicholas M. Law, Natasha E. Batalha, Andrew Vanderburg, Diana Dragomir, Chelsea X. Huang, Fabo Feng, R. Paul Butler, Jeffrey D. Crane, Stephen A. Shectman, Yuri Beletsky, Avi Shporer, Benjamin T. Montet, Jennifer A. Burt, Adina D. Feinstein, Erin Flowers, Sangeetha Nandakumar, Mauro Barbieri, Hank Corbett, Jeffrey K. Ratzloff, Nathan Galliher , et al. (4 additional authors not shown)

Abstract: HD 21749 is a bright ($V=8.1$ mag) K dwarf at 16 pc known to host an inner terrestrial planet HD 21749c as well as an outer sub-Neptune HD 21749b, both delivered by TESS. Follow-up spectroscopic observations measured the mass of HD 21749b to be $22.7\pm2.2\ M_{\oplus}$ with a density of $7.0^{+1.6}_{-1.3}$ g~cm$^{-3}$, making it one of the densest sub-Neptunes. However, the mass measurement was su… ▽ More HD 21749 is a bright ($V=8.1$ mag) K dwarf at 16 pc known to host an inner terrestrial planet HD 21749c as well as an outer sub-Neptune HD 21749b, both delivered by TESS. Follow-up spectroscopic observations measured the mass of HD 21749b to be $22.7\pm2.2\ M_{\oplus}$ with a density of $7.0^{+1.6}_{-1.3}$ g~cm$^{-3}$, making it one of the densest sub-Neptunes. However, the mass measurement was suspected to be influenced by stellar rotation. Here we present new high-cadence PFS RV data to disentangle the stellar activity signal from the planetary signal. We find that HD 21749 has a similar rotational timescale as the planet's orbital period, and the amplitude of the planetary orbital RV signal is estimated to be similar to that of the stellar activity signal. We perform Gaussian Process (GP) regression on the photometry and RVs from HARPS and PFS to model the stellar activity signal. Our new models reveal that HD 21749b has a radius of $2.86\pm0.20\ R_{\oplus}$, an orbital period of $35.6133\pm0.0005$ d with a mass of $M_{b}=20.0\pm2.7\ M_{\oplus}$ and a density of $4.8^{+2.0}_{-1.4}$ g~cm$^{-3}$ on an eccentric orbit with $e=0.16\pm0.06$, which is consistent with the most recent values published for this system. HD 21749c has an orbital period of $7.7902\pm0.0006$ d, a radius of $1.13\pm0.10\ R_{\oplus}$, and a 3$σ$ mass upper limit of $3.5\ M_{\oplus}$. Our Monte Carlo simulations confirm that without properly taking stellar activity signals into account, the mass measurement of HD 21749b is likely to arrive at a significantly underestimated error bar. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Comments: 20 pages, 16 figures, accepted to MNRAS

Orbital Foregrounds for Ultra-Short Duration Transients

Authors: Hank Corbett, Nicholas M. Law, Alan Vasquez Soto, Ward S. Howard, Amy Glazier, Ramses Gonzalez, Jeffrey K. Ratzloff, Nathan Galliher, Octavi Fors, Robert Quimby

Abstract: Reflections from objects in Earth orbit can produce sub-second, star-like optical flashes similar to astrophysical transients. Reflections have historically caused false alarms for transient surveys, but the population has not been systematically studied. We report event rates for these orbital flashes using the Evryscope Fast Transient Engine, a low-latency transient detection pipeline for the Ev… ▽ More Reflections from objects in Earth orbit can produce sub-second, star-like optical flashes similar to astrophysical transients. Reflections have historically caused false alarms for transient surveys, but the population has not been systematically studied. We report event rates for these orbital flashes using the Evryscope Fast Transient Engine, a low-latency transient detection pipeline for the Evryscopes. We select single-epoch detections likely caused by Earth satellites and model the event rate as a function of both magnitude and sky position. We measure a rate of $1800^{+600}_{-280}$ sky$^{-1}$ hour$^{-1}$, peaking at $m_g = 13.0$, for flashes morphologically degenerate with real astrophysical signals in surveys like the Evryscopes. Of these, $340^{+150}_{-85}$ sky$^{-1}$ hour$^{-1}$ are bright enough to be visible to the naked eye in typical suburban skies with a visual limiting magnitude of $V\approx4$. These measurements place the event rate of orbital flashes orders of magnitude higher than the combined rate of public alerts from all active all-sky fast-timescale transient searches, including neutrino, gravitational-wave, gamma-ray, and radio observatories. Short-timescale orbital flashes form a dominating foreground for un-triggered searches for fast transients in low-resolution, wide-angle surveys. However, events like fast radio bursts (FRBs) with arcminute-scale localization have a low probability ($\sim10^{-5}$) of coincidence with an orbital flash, allowing optical surveys to place constraints on their potential optical counterparts in single images. Upcoming satellite internet constellations, like SpaceX Starlink, are unlikely to contribute significantly to the population of orbital flashes in normal operations. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: 8 pages, 4 figures

Journal ref: ApJL, 903, L27 (2020)

EvryFlare III: Temperature Evolution and Habitability Impacts of Dozens of Superflares Observed Simultaneously by Evryscope and TESS

Authors: Ward S. Howard, Hank Corbett, Nicholas M. Law, Jeffrey K. Ratzloff, Nathan Galliher, Amy L. Glazier, Ramses Gonzalez, Alan Vasquez Soto, Octavi Fors, Daniel del Ser, Joshua Haislip

Abstract: Superflares may provide the dominant source of biologically relevant UV radiation to rocky habitable zone M-dwarf planets (M-Earths), altering planetary atmospheres and conditions for surface life. The combined line and continuum flare emission has usually been approximated by a 9000 K blackbody. If superflares are hotter, then the UV emission may be 10X higher than predicted from the optical. How… ▽ More Superflares may provide the dominant source of biologically relevant UV radiation to rocky habitable zone M-dwarf planets (M-Earths), altering planetary atmospheres and conditions for surface life. The combined line and continuum flare emission has usually been approximated by a 9000 K blackbody. If superflares are hotter, then the UV emission may be 10X higher than predicted from the optical. However, it is unknown for how long M-dwarf superflares reach temperatures above 9000 K. Only a handful of M-dwarf superflares have been recorded with multi-wavelength high-cadence observations. We double the total number of events in the literature using simultaneous Evryscope and TESS observations to provide the first systematic exploration of the temperature evolution of M-dwarf superflares. We also increase the number of superflaring M-dwarfs with published time-resolved blackbody evolution by ~10X. We measure temperatures at 2 min cadence for 42 superflares from 27 K5-M5 dwarfs. We find superflare peak temperatures (defined as the mean of temperatures corresponding to flare FWHM) increase with flare energy and impulse. We find the amount of time flares emit at temperatures above 14,000 K depends on energy. We discover 43% of the flares emit above 14,000 K, 23% emit above 20,000 K and 5% emit above 30,000 K. The largest and hottest flare briefly reached 42,000 K. Some do not reach 14,000 K. During superflares, we estimate M-Earths orbiting <200 Myr stars typically receive a top-of-atmosphere UV-C flux of ~120 W m^-2 and up to 10^3 W m^-2, 100-1000X the time-averaged XUV flux from Proxima Cen. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 14 pages, 9 figures, 2 tables. Accepted for publication to the Astrophysical Journal. In press

EVR-CB-004: An Inflated Hot Subdwarf O star + Unseen WD Companion in a Compact Binary Discovered with the Evryscope

Authors: Jeffrey K. Ratzloff, Thomas Kupfer, Brad N. Barlow, David Schneider, Thomas R. Marsh, Ulrich Heber, Kyle A. Corcoran, Evan Bauer, Steven Hammerich, Henry T. Corbett, Amy Glazier, Ward S. Howard, Nicholas M. Law

Abstract: We present the discovery of EVR-CB-004, a close binary with a remnant stellar core and an unseen white dwarf companion. The analysis in this work reveals the primary is potentially an inflated hot subdwarf (sdO) and more likely is a rarer post-blue horizontal branch (post-BHB) star. Post-BHBs are the short-lived shell-burning final stage of a blue horizontal star or hot subdwarf before transitioni… ▽ More We present the discovery of EVR-CB-004, a close binary with a remnant stellar core and an unseen white dwarf companion. The analysis in this work reveals the primary is potentially an inflated hot subdwarf (sdO) and more likely is a rarer post-blue horizontal branch (post-BHB) star. Post-BHBs are the short-lived shell-burning final stage of a blue horizontal star or hot subdwarf before transitioning to a WD. This object was discovered using Evryscope photometric data in a southern-all-sky hot subdwarf variability survey. The photometric light curve for EVR-CB-004 shows multi-component variability from ellipsoidal deformation of the primary and from Doppler boosting as well as gravitational limb darkening. EVR-CB-004 is one of just a handful of known systems, and has a long period (6.08426 hours) and large amplitude ellipsoidal modulation (16.0 $\%$ change in brightness from maximum to minimum) for these extremely close binary systems, while the properties of the primary make it a truly unique system. EVR-CB-004 also shows a peculiar low-amplitude (less than $1\%$) sinusoidal light curve variation with a period that is a 1/3 resonance of the binary period. We tentatively identify this additional variation source as a tidally-induced resonant pulsation, and we suggest followup observations that could verify this interpretation. From the evolutionary state of the system, its components, and its mass fraction, EVR-CB-004 is a strong merger candidate to form a single high-mass ($\approx1.2M_{\odot}$) WD. EVR-CB-004 offers a glimpse into a brief phase of a remnant core evolution and secondary variation, not seen before in a compact binary. △ Less

Submitted 13 September, 2020; originally announced September 2020.

Comments: Accepted to ApJ September 7, 2020

Evryscope-South Survey of Upper- and Pre-main Sequence Solar Neighborhood Stars

Authors: Nathan W. Galliher, Jeffrey K. Ratzloff, Henry Corbett, Nicholas M. Law, Ward S. Howard, Amy L. Glazier, Alan Vasquez Soto, Ramses Gonzalez

Abstract: Using photometric data collected by Evryscope-South, we search for nearby young variable systems on the upper-main sequence (UMS) and pre-main sequence (PMS). The Evryscopes are all-sky high-cadence telescope arrays operating in the Northern and Southern hemispheres. We base our search on a Gaia-selected catalog of young neighborhood upper- and pre-main sequence stars which were chosen through bot… ▽ More Using photometric data collected by Evryscope-South, we search for nearby young variable systems on the upper-main sequence (UMS) and pre-main sequence (PMS). The Evryscopes are all-sky high-cadence telescope arrays operating in the Northern and Southern hemispheres. We base our search on a Gaia-selected catalog of young neighborhood upper- and pre-main sequence stars which were chosen through both astrometric and photometric criteria. We analyze 44,971 Evryscope-South light curves in search of variability. We recover 615 variables, with 378 previously known, and 237 new discoveries including 84 young eclipsing binary (EB) candidates. We discover a new highly eccentric binary system and recover a further four previously known systems, with periods ranging from 299 to 674 hr. We find 158 long-period (>50 hr) candidate EB systems, 9 from the PMS and 149 from the UMS, which will allow constraints on the mass-radius-age relation. These long-period EBs include a 179.3 hr PMS system and a 867.8 hr system from the UMS. For PMS variable candidates we estimate system ages, which range from 1 to 23 Myr for non-EBs and from 2 to 17 Myr for EBs. Other non-EB discoveries that show intrinsic variability will allow relationships between stellar rotation rates, ages, activity, and mass to be characterized. △ Less

Submitted 28 April, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: Final version published in PASP, typos corrected, 16 pages, 15 figures, 4 tables, 2 supplementary tables

Evryscope and K2 Constraints on TRAPPIST-1 Superflare Occurrence and Planetary Habitability

Authors: Amy L. Glazier, Ward S. Howard, Hank Corbett, Nicholas M. Law, Jeffrey K. Ratzloff, Octavi Fors, Daniel del Ser

Abstract: The nearby ultracool dwarf TRAPPIST-1 possesses several Earth-sized terrestrial planets, three of which have equilibrium temperatures that may support liquid surface water, making it a compelling target for exoplanet characterization. TRAPPIST-1 is an active star with frequent flaring, with implications for the habitability of its planets. Superflares (stellar flares whose energy exceeds 10^33 erg… ▽ More The nearby ultracool dwarf TRAPPIST-1 possesses several Earth-sized terrestrial planets, three of which have equilibrium temperatures that may support liquid surface water, making it a compelling target for exoplanet characterization. TRAPPIST-1 is an active star with frequent flaring, with implications for the habitability of its planets. Superflares (stellar flares whose energy exceeds 10^33 erg) can completely destroy the atmospheres of a cool star's planets, allowing ultraviolet radiation and high-energy particles to bombard their surfaces. However, ultracool dwarfs emit little ultraviolet flux when quiescent, raising the possibility of frequent flares being necessary for prebiotic chemistry that requires ultraviolet light. We combine Evryscope and Kepler observations to characterize the high-energy flare rate of TRAPPIST-1. The Evryscope is an array of 22 small telescopes imaging the entire Southern sky in g' every two minutes. Evryscope observations, spanning 170 nights over 2 years, complement the 80-day continuous short-cadence K2 observations by sampling TRAPPIST-1's long-term flare activity. We update TRAPPIST-1's superflare rate, finding a cumulative rate of 4.2 (+1.9 -0.2) superflares per year. We calculate the flare rate necessary to deplete ozone in the habitable-zone planets' atmospheres, and find that TRAPPIST-1's flare rate is insufficient to deplete ozone if present on its planets. In addition, we calculate the flare rate needed to provide enough ultraviolet flux to power prebiotic chemistry. We find TRAPPIST-1's flare rate is likely insufficient to catalyze some of the Earthlike chemical pathways thought to lead to RNA synthesis, and flux due to flares in the biologically relevant UV-B band is orders of magnitude less for any TRAPPIST-1 planet than has been experienced by Earth at any time in its history. △ Less

Submitted 17 August, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: 12 pages, 9 figures. Accepted to The Astrophysical Journal, in press

Multi-wavelength Photometry and Progenitor Analysis of the Nova V906 Car

Authors: Jerrick Wee, Nadejda Blagorodnova, Bryan Edward Penprase, Jett Pierce Facey, Taiga Morioka, Hank Corbett, Brad N. Barlow, Thomas Kupfer, Nicholas M. Law, Jeffrey K. Ratzloff, Ward S. Howard, Ramses Gonzalez Chavez, Amy Glazier, Alan Vasquez Soto, Takashi Horiuchi

Abstract: We present optical and infrared photometry of the classical nova V906 Car, also known as Nova Car 2018 and ASASSN-18fv, discovered by ASASS-SN survey on 16.32 March 2018 UT (MJD 58193.0). The nova reached its maximum on MJD 58222.56 at $V_{\rm{max}} = 5.84 \pm 0.09$ mag and had decline times of $t_{2,V} = 26.2 $ d and $t_{3,V} = 33.0 $ d. The data from Evryscope shows that the nova had already bri… ▽ More We present optical and infrared photometry of the classical nova V906 Car, also known as Nova Car 2018 and ASASSN-18fv, discovered by ASASS-SN survey on 16.32 March 2018 UT (MJD 58193.0). The nova reached its maximum on MJD 58222.56 at $V_{\rm{max}} = 5.84 \pm 0.09$ mag and had decline times of $t_{2,V} = 26.2 $ d and $t_{3,V} = 33.0 $ d. The data from Evryscope shows that the nova had already brightened to $g'\simeq 13$\,mag five days before discovery, as compared to its quiescent magnitude of $g=$20.13$\pm$0.03. The extinction towards the nova, as derived from high resolution spectroscopy, shows an estimate consistent with foreground extinction to the Carina Nebula of $A_V = 1.11_{-0.39}^{+0.54}$. The light curve resembles a rare C (cusp) class nova with a steep decline slope of $α=-3.94$ post cusp flare. From the lightcurve decline rate, we estimate the mass of white dwarf to be $M_{WD}$ = $ < 0.8$M\textsubscript{\(\odot\)}, consistent with $M_{WD}=0.71^{+0.23}_{-0.19}$ derived from modelling the accretion disk of the system in quiescence. The donor star is likely a K-M dwarf of 0.23-0.43\,\Msun, which is being heated by its companion. △ Less

Submitted 25 June, 2020; originally announced June 2020.

Comments: 13 pages, 10 figures, 7 tables. Accepted to ApJ

Hot Subdwarf All Southern Sky Fast Transit Survey with the Evryscope

Authors: Jeffrey K. Ratzloff, Brad N. Barlow, Peter Nemeth, Henry T. Corbett, Stephen Walser, Nathan W. Galliher, Amy Glazier, Ward S. Howard, Nicholas M. Law

Abstract: We have conducted a survey of candidate hot subdwarf stars in the southern sky searching for fast transits, eclipses, and sinusoidal like variability in the Evryscope light curves. The survey aims to detect transit signals from Neptune size planets to gas-giants, and eclipses from M-dwarfs and brown dwarfs. The other variability signals are primarily expected to be from compact binaries and reflec… ▽ More We have conducted a survey of candidate hot subdwarf stars in the southern sky searching for fast transits, eclipses, and sinusoidal like variability in the Evryscope light curves. The survey aims to detect transit signals from Neptune size planets to gas-giants, and eclipses from M-dwarfs and brown dwarfs. The other variability signals are primarily expected to be from compact binaries and reflection effect binaries. Due to the small size of hot subdwarfs, transit and eclipse signals are expected to last only twenty minutes, but with large signal depths (up to completely eclipsing if the orientation is edge on). With its 2-minute cadence and continuous observing Evryscope is well placed to recover these fast transits and eclipses. The very large field of view (8150 sq. deg.) is critical to obtain enough hot subdwarf targets, despite their rarity. We identified 11,000 potential hot subdwarfs from the 9.3M Evryscope light curves for sources brighter than mg = 15. With our machine learning spectral classifier, we flagged high-confidence targets and estimate the total hot subdwarfs in the survey to be 1400. The light curve search detected three planet transit candidates, shown to have stellar companions from followup analysis. We discovered several new compact binaries (including two with unseen degenerate companions, and several others with potentially rare secondaries), two eclipsing binaries with M-dwarf companions, as well as new reflection effect binaries and others with sinusoidal like variability. The hot subdwarf discoveries identified here are spectroscopically confirmed and we verified the Evryscope discovery light curve with TESS light curves when available. Four of the discoveries are in the process of being published in separate followup papers, and we discuss the followup potential of several of the other discoveries. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: Accepted to ApJ, December 2019

EVR-CB-001: An evolving, progenitor, white dwarf compact binary discovered with the Evryscope

Authors: Jeffrey K. Ratzloff, Brad N. Barlow, Thomas Kupfer, Kyle A. Corcoran, Stephan Geier, Evan Bauer, Henry T. Corbett, Ward S. Howard, Amy Glazier, Nicholas M. Law

Abstract: We present EVR-CB-001, the discovery of a compact binary with an extremely low mass ($.21 \pm 0.05 M_{\odot}$) helium core white dwarf progenitor (pre-He WD) and an unseen low mass ($.32 \pm 0.06 M_{\odot}$) helium white dwarf (He WD) companion. He WDs are thought to evolve from the remnant helium-rich core of a main-sequence star stripped during the giant phase by a close companion. Low mass He W… ▽ More We present EVR-CB-001, the discovery of a compact binary with an extremely low mass ($.21 \pm 0.05 M_{\odot}$) helium core white dwarf progenitor (pre-He WD) and an unseen low mass ($.32 \pm 0.06 M_{\odot}$) helium white dwarf (He WD) companion. He WDs are thought to evolve from the remnant helium-rich core of a main-sequence star stripped during the giant phase by a close companion. Low mass He WDs are exotic objects (only about .2$\%$ of WDs are thought to be less than .3 $M_{\odot}$), and are expected to be found in compact binaries. Pre-He WDs are even rarer, and occupy the intermediate phase after the core is stripped, but before the star becomes a fully degenerate WD and with a larger radius ($\approx .2 R_{\odot}$) than a typical WD. The primary component of EVR-CB-001 (the pre-He WD) was originally thought to be a hot subdwarf (sdB) star from its blue color and under-luminous magnitude, characteristic of sdBs. The mass, temperature ($T_{\rm eff}=18,500 \pm 500 K$), and surface gravity ($\log(g)=4.96 \pm 0.04$) solutions from this work are lower than values for typical hot subdwarfs. The primary is likely to be a post-RGB, pre-He WD contracting into a He WD, and at a stage that places it nearest to sdBs on color-magnitude and $T_{\rm eff}$-$\log(g)$ diagrams. EVR-CB-001 is expected to evolve into a fully double degenerate, compact system that should spin down and potentially evolve into a single hot subdwarf star. Single hot subdwarfs are observed, but progenitor systems have been elusive. △ Less

Submitted 3 October, 2019; v1 submitted 4 September, 2019; originally announced September 2019.

Comments: 14 pages, 11 figures. Published in The Astrophysical Journal

EvryFlare II: Rotation Periods of the Cool Flare Stars in TESS Across Half the Southern Sky

Authors: Ward S. Howard, Hank Corbett, Nicholas M. Law, Jeffrey K. Ratzloff, Nathan Galliher, Amy Glazier, Octavi Fors, Daniel del Ser, Joshua Haislip

Abstract: We measure rotation periods and sinusoidal amplitudes in Evryscope light curves for 122 two-minute K5-M4 TESS targets selected for strong flaring. The Evryscope array of telescopes has observed all bright nearby stars in the South, producing two-minute cadence light curves since 2016. Long-term, high-cadence observations of rotating flare stars probe the complex relationship between stellar rotati… ▽ More We measure rotation periods and sinusoidal amplitudes in Evryscope light curves for 122 two-minute K5-M4 TESS targets selected for strong flaring. The Evryscope array of telescopes has observed all bright nearby stars in the South, producing two-minute cadence light curves since 2016. Long-term, high-cadence observations of rotating flare stars probe the complex relationship between stellar rotation, starspots, and superflares. We detect periods from 0.3487 to 104 d, and observe amplitudes from 0.008 to 0.216 g' mag. We find the Evryscope amplitudes are larger than those in TESS with the effect correlated to stellar mass (p-value=0.01). We compute the Rossby number (Ro), and find our sample selected for flaring has twice as many intermediate rotators (0.04<Ro<0.4) as fast (Ro<0.04) or slow (Ro>0.44) rotators; this may be astrophysical or a result of period-detection sensitivity. We discover 30 fast, 59 intermediate, and 33 slow rotators. We measure a median starspot coverage of 13% of the stellar hemisphere and constrain the minimum magnetic field strength consistent with our flare energies and spot coverage to be 500 G, with later-type stars exhibiting lower values than earlier-types. We observe a possible change in superflare rates at intermediate periods. However, we do not conclusively confirm the increased activity of intermediate rotators seen in previous studies. We split all rotators at Ro~0.2 into Prot<10 d and Prot>10 d bins to confirm short-period rotators exhibit higher superflare rates, larger flare energies, and higher starspot coverage than do long-period rotators, at p-values of 3.2 X 10^-5, 1.0 X 10^-5, and 0.01, respectively. △ Less

Submitted 26 May, 2020; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: 16 pages, 8 figures, 3 tables. Ancillary machine-readable files included. Accepted for publication in ApJ (proofs submitted). Includes significant new material, including starspot color that depends on stellar mass, more rotation periods, potential changes in activity during spin-down, and examples of binary rotators

Variables in the Southern Polar Region Evryscope 2016 Dataset

Authors: Jeffrey K. Ratzloff, Henry T. Corbett, Nicholas M. Law, Brad N. Barlow, Amy Glazier, Ward S. Howard, Octavi Fors, Daniel del Ser, Trifon Trifonov

Abstract: The regions around the celestial poles offer the ability to find and characterize long-term variables from ground-based observatories. We used multi-year Evryscope data to search for high-amplitude (~5% or greater) variable objects among 160,000 bright stars (Mv < 14.5) near the South Celestial Pole. We developed a machine learning based spectral classifier to identify eclipse and transit candidat… ▽ More The regions around the celestial poles offer the ability to find and characterize long-term variables from ground-based observatories. We used multi-year Evryscope data to search for high-amplitude (~5% or greater) variable objects among 160,000 bright stars (Mv < 14.5) near the South Celestial Pole. We developed a machine learning based spectral classifier to identify eclipse and transit candidates with M-dwarf or K-dwarf host stars - and potential low-mass secondary stars or gas giant planets. The large amplitude transit signals from low-mass companions of smaller dwarf host stars lessens the photometric precision and systematics removal requirements necessary for detection, and increases the discoveries from long-term observations with modest light curve precision. The Evryscope is a robotic telescope array that observes the Southern sky continuously at 2-minute cadence, searching for stellar variability, transients, transits around exotic stars and other observationally challenging astrophysical variables. In this study, covering all stars 9 < Mv < 14.5, in declinations -75 to -90 deg, we recover 346 known variables and discover 303 new variables, including 168 eclipsing binaries. We characterize the discoveries and provide the amplitudes, periods, and variability type. A 1.7 Jupiter radius planet candidate with a late K-dwarf primary was found and the transit signal was verified with the PROMPT telescope network. Further followup revealed this object to be a likely grazing eclipsing binary system with nearly identical primary and secondary K5 stars. Radial velocity measurements from the Goodman Spectrograph on the 4.1 meter SOAR telescope of the likely-lowest-mass targets reveal that six of the eclipsing binary discoveries are low-mass (.06 - .37 solar mass) secondaries with K-dwarf primaries, strong candidates for precision mass-radius measurements. △ Less

Submitted 7 May, 2019; originally announced May 2019.

Comments: 32 pages, 17 figures, accepted to PASP

Building the Evryscope: Hardware Design and Performance

Authors: Jeffrey K. Ratzloff, Nicholas M. Law, Octavi Fors, Henry T. Corbett, Ward S. Howard, Daniel del Ser, Joshua Haislip

Abstract: The Evryscope is a telescope array designed to open a new parameter space in optical astronomy, detecting short timescale events across extremely large sky areas simultaneously. The system consists of a 780 MPix 22-camera array with an 8150 sq. deg. field of view, 13" per pixel sampling, and the ability to detect objects down to Mg=16 in each 2 minute dark-sky exposure. The Evryscope, covering 18,… ▽ More The Evryscope is a telescope array designed to open a new parameter space in optical astronomy, detecting short timescale events across extremely large sky areas simultaneously. The system consists of a 780 MPix 22-camera array with an 8150 sq. deg. field of view, 13" per pixel sampling, and the ability to detect objects down to Mg=16 in each 2 minute dark-sky exposure. The Evryscope, covering 18,400 sq.deg. with hours of high-cadence exposure time each night, is designed to find the rare events that require all-sky monitoring, including transiting exoplanets around exotic stars like white dwarfs and hot subdwarfs, stellar activity of all types within our galaxy, nearby supernovae, and other transient events such as gamma ray bursts and gravitational-wave electromagnetic counterparts. The system averages 5000 images per night with ~300,000 sources per image, and to date has taken over 3.0M images, totaling 250TB of raw data. The resulting light curve database has light curves for 9.3M targets, averaging 32,600 epochs per target through 2018. This paper summarizes the hardware and performance of the Evryscope, including the lessons learned during telescope design, electronics design, a procedure for the precision polar alignment of mounts for Evryscope-like systems, robotic control and operations, and safety and performance-optimization systems. We measure the on-sky performance of the Evryscope, discuss its data-analysis pipelines, and present some example variable star and eclipsing binary discoveries from the telescope. We also discuss new discoveries of very rare objects including 2 hot subdwarf eclipsing binaries with late M-dwarf secondaries (HW Vir systems), 2 white dwarf / hot subdwarf short-period binaries, and 4 hot subdwarf reflection binaries. We conclude with the status of our transit surveys, M-dwarf flare survey, and transient detection. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 24 pages, 24 figures, accepted PASP

EvryFlare I: Long-term Evryscope Monitoring of Flares from the Cool Stars Across Half the Southern Sky

Authors: Ward S. Howard, Hank Corbett, Nicholas M. Law, Jeffrey K. Ratzloff, Amy L. Glazier, Octavi Fors, Daniel del Ser, Joshua Haislip

Abstract: We search for superflares from 4,068 cool stars in 2+ years of Evryscope photometry, focusing on those with high-cadence data from both Evryscope and TESS. The Evryscope array of small telescopes observed 575 flares from 284 stars, with a median energy of 10^34.0 erg. Since 2016, Evryscope has enabled the detection of rare events from all stars observed by TESS through multi-year, high-cadence con… ▽ More We search for superflares from 4,068 cool stars in 2+ years of Evryscope photometry, focusing on those with high-cadence data from both Evryscope and TESS. The Evryscope array of small telescopes observed 575 flares from 284 stars, with a median energy of 10^34.0 erg. Since 2016, Evryscope has enabled the detection of rare events from all stars observed by TESS through multi-year, high-cadence continuous observing. We report ~2X the previous largest number of 10^34 erg high-cadence flares from nearby cool stars. We find 8 flares with amplitudes of 3+ g' magnitudes, with the largest reaching 5.6 magnitudes and releasing 10^36.2 erg. We observe a 10^34 erg superflare from TOI-455 (LTT 1445), a mid-M with a rocky planet candidate. We measure the superflare rate per flare-star and quantify the average flaring of active stars as a function of spectral type, including superflare rates, FFDs, and typical flare amplitudes in g'. We confirm superflare morphology is broadly consistent with magnetic re-connection. We estimate starspot coverage necessary to produce superflares, and hypothesize maximum-allowed superflare energies and waiting-times between flares corresponding to 100% coverage of the stellar hemisphere. We observe decreased flaring at high galactic latitudes. We explore the effects of superflares on ozone loss to planetary atmospheres: we observe 1 superflare with sufficient energy to photo-dissociate all ozone in an Earth-like atmosphere in one event. We find 17 stars that may deplete an Earth-like atmosphere via repeated flaring. Of the 1822 stars around which TESS may discover temperate rocky planets, we observe 14.6% +/- 2% emit large flares. △ Less

Submitted 23 July, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 18 pages, 14 figures, 5 tables. Ancillary machine-readable files included. Accepted for publication in ApJ. Includes new sections on starspot coverage and superflare emission from starspots

Bright Opportunities for Atmospheric Characterization of Small Planets: Masses and Radii of K2-3 b, c, d and GJ3470 b from Radial Velocity Measurements and Spitzer Transits

Authors: Molly R. Kosiarek, Ian J. M. Crossfield, Kevin K. Hardegree-Ullman, John H. Livingston, Bjorn Benneke, Sarah Blunt, Gregory W. Henry, Ward S. Howard, David Berardo, Benjamin J. Fulton, Lea A. Hirsch, Andrew W. Howard, Howard Isaacson, Erik A. Petigura, Evan Sinukoff, Lauren Weiss, X. Bonfils, Courtney D. Dressing, Heather A. Knutson, Joshua E. Schlieder, Michael Werner, Varoujan Gorjian, Jessica Krick, Farisa Y. Morales, Nicola Astudillo-Defru , et al. (14 additional authors not shown)

Abstract: We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairb… ▽ More We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairborn Observatory to determine the characteristic stellar activity timescales for our Gaussian Process fit, including the stellar rotation period and activity region decay timescale. The stellar rotation signals for both stars are evident in the radial velocity data and are included in our fit using a Gaussian process trained on the photometry. We find the masses of K2-3 b, K2-3 c and GJ3470 b to be 6.48$^{+0.99}_{-0.93}$, 2.14$^{+1.08}_{-1.04}$, and 12.58$^{+1.31}_{-1.28}$ M$_\oplus$ respectively. K2-3 d was not significantly detected and has a 3-$σ$ upper limit of 2.80 M$_\oplus$. These two systems are training cases for future TESS systems; due to the low planet densities ($ρ$ $<$ 3.7 g cm$^{-3}$) and bright host stars (K $<$ 9 mag), they are among the best candidates for transmission spectroscopy in order to characterize the atmospheric compositions of small planets. △ Less

Submitted 14 February, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 21 pages, 11 figures, accepted to AJ

Journal ref: AJ 157, 97 (2019)

On the Origin of Multiple Populations During Massive Star Cluster Formation

Authors: Corey S. Howard, Ralph E. Pudritz, Alison Sills, William E. Harris

Abstract: We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accret… ▽ More We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accretion and rapid merging of protoclusters. Chemical enrichment is a dynamic process happening as the young cluster assembles, so that the original (1P) and enriched (2P) subpopulations of stars form almost simultaneously. Here we test two simple and opposite extremes for the injection of enriched material into the intracluster gas: we assume either continuous injection in a way that tracks the star formation rate; or sudden injection by a single instantaneous event. Using helium abundance as a proxy for the enrichment, we find that realistic multiple population features can be reproduced by injecting a total helium mass amounting to a few percent of the cluster's total mass. The differences in individual growth histories can lead to widely differing 1P/2P outcomes. These models suggest that dual or multiple populations can emerge rapidly in massive star clusters undergoing the typical mode of star cluster formation. △ Less

Submitted 5 April, 2019; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: Accepted to MNRAS. 11 pages, 5 figures

A universal route for the formation of massive star clusters in giant molecular clouds

Authors: Corey S. Howard, Ralph E. Pudritz, William E. Harris

Abstract: Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show tha… ▽ More Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show that GC-sized clusters are naturally produced in radiation-hydrodynamic simulations of isolated 10$^7$ M$_{\odot}$ Giant Molecular Clouds (GMCs) with properties typical of the local universe, even under the influence of radiative feedback. In all cases, these massive clusters grow to GC-level masses within 5 Myr via a roughly equal combination of filamentary gas accretion and mergers with several less massive clusters. Lowering the heavy-element abundance of the GMC by a factor of 10 reduces the opacity of the gas to radiation and better represents the high-redshift formation conditions of GCs. This results in higher gas accretion leading to a mass increase of the largest cluster by a factor of ~4. When combined with simulations of less massive GMCs (10$^{4-6}$ M$_{\odot}$), a clear relation emerges between the maximum YMC mass and the mass of the host GMC. Our results demonstrate that YMCs, and potentially GCs, are a simple extension of local cluster formation to more massive clouds and do not require suggested exotic formation scenarios. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: Submitted version to Nature Astronomy. Final published paper, and online supplementary material, available at: https://www.nature.com/articles/s41550-018-0506-0

The First Naked-Eye Superflare Detected from Proxima Centauri

Authors: Ward S. Howard, Matt A. Tilley, Hank Corbett, Allison Youngblood, R. O. Parke Loyd, Jeffrey K. Ratzloff, Nicholas M. Law, Octavi Fors, Daniel del Ser, Evgenya L. Shkolnik, Carl Ziegler, Erin E. Goeke, Aaron D. Pietraallo, Joshua Haislip

Abstract: Proxima b is a terrestrial-mass planet in the habitable-zone of Proxima Centauri. Proxima Centauri's high stellar activity however casts doubt on the habitability of Proxima b: sufficiently bright and frequent flares and any associated proton events may destroy the planet's ozone layer, allowing lethal levels of UV flux to reach its surface. In March 2016, the Evryscope observed the first naked-ey… ▽ More Proxima b is a terrestrial-mass planet in the habitable-zone of Proxima Centauri. Proxima Centauri's high stellar activity however casts doubt on the habitability of Proxima b: sufficiently bright and frequent flares and any associated proton events may destroy the planet's ozone layer, allowing lethal levels of UV flux to reach its surface. In March 2016, the Evryscope observed the first naked-eye-brightness superflare detected from Proxima Centauri. Proxima increased in optical flux by a factor of ~68 during the superflare and released a bolometric energy of 10^33.5 erg, ~10X larger than any previously-detected flare from Proxima. Over the last two years the Evryscope has recorded 23 other large Proxima flares ranging in bolometric energy from 10^30.6 erg to 10^32.4 erg; coupling those rates with the single superflare detection, we predict at least five superflares occur each year. Simultaneous high-resolution HARPS spectroscopy during the Evryscope superflare constrains the superflare's UV spectrum and any associated coronal mass ejections. We use these results and the Evryscope flare rates to model the photochemical effects of NOx atmospheric species generated by particle events from this extreme stellar activity, and show that the repeated flaring may be sufficient to reduce the ozone of an Earth-like atmosphere by 90% within five years; complete depletion may occur within several hundred kyr. The UV light produced by the Evryscope superflare would therefore have reached the surface with ~100X the intensity required to kill simple UV-hardy microorganisms, suggesting that life would have to undergo extreme adaptations to survive in the surface areas of Proxima b exposed to these flares. △ Less

Submitted 7 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 7 pages, 5 figures, accepted in ApJ Letters. Includes minor changes, extra habitability calculations

Cataloging the Visible Universe through Bayesian Inference at Petascale

Authors: Jeffrey Regier, Kiran Pamnany, Keno Fischer, Andreas Noack, Maximilian Lam, Jarrett Revels, Steve Howard, Ryan Giordano, David Schlegel, Jon McAuliffe, Rollin Thomas, Prabhat

Abstract: Astronomical catalogs derived from wide-field imaging surveys are an important tool for understanding the Universe. We construct an astronomical catalog from 55 TB of imaging data using Celeste, a Bayesian variational inference code written entirely in the high-productivity programming language Julia. Using over 1.3 million threads on 650,000 Intel Xeon Phi cores of the Cori Phase II supercomputer… ▽ More Astronomical catalogs derived from wide-field imaging surveys are an important tool for understanding the Universe. We construct an astronomical catalog from 55 TB of imaging data using Celeste, a Bayesian variational inference code written entirely in the high-productivity programming language Julia. Using over 1.3 million threads on 650,000 Intel Xeon Phi cores of the Cori Phase II supercomputer, Celeste achieves a peak rate of 1.54 DP PFLOP/s. Celeste is able to jointly optimize parameters for 188M stars and galaxies, loading and processing 178 TB across 8192 nodes in 14.6 minutes. To achieve this, Celeste exploits parallelism at multiple levels (cluster, node, and thread) and accelerates I/O through Cori's Burst Buffer. Julia's native performance enables Celeste to employ high-level constructs without resorting to hand-written or generated low-level code (C/C++/Fortran), and yet achieve petascale performance. △ Less

Submitted 30 January, 2018; originally announced January 2018.

Comments: accepted to IPDPS 2018

MSC Class: 85A35; 68W10; 62P35 ACM Class: J.2; D.1.3; G.3; I.2; D.2

Laser-only adaptive optics achieves significant image quality gains compared to seeing-limited observations over the entire sky

Authors: Ward S. Howard, Nicholas M. Law, Carl A. Ziegler, Christoph Baranec, Reed Riddle

Abstract: Adaptive optics laser guide star systems perform atmospheric correction of stellar wavefronts in two parts: stellar tip-tilt and high-spatial-order laser-correction. The requirement of a sufficiently bright guide star in the field-of-view to correct tip-tilt limits sky coverage. Here we show an improvement to effective seeing without the need for nearby bright stars, enabling full sky coverage by… ▽ More Adaptive optics laser guide star systems perform atmospheric correction of stellar wavefronts in two parts: stellar tip-tilt and high-spatial-order laser-correction. The requirement of a sufficiently bright guide star in the field-of-view to correct tip-tilt limits sky coverage. Here we show an improvement to effective seeing without the need for nearby bright stars, enabling full sky coverage by performing only laser-assisted wavefront correction. We used Robo-AO, the first robotic AO system, to comprehensively demonstrate this laser-only correction. We analyze observations from four years of efficient robotic operation covering 15,000 targets and 42,000 observations, each realizing different seeing conditions. Using an autoguider (or a post-processing software equivalent) and the laser to improve effective seeing independent of the brightness of a target, Robo-AO observations show a 39+/-19% improvement to effective FWHM, without any tip-tilt correction. We also demonstrate that 50% encircled-energy performance without tip-tilt correction remains comparable to diffraction-limited, standard Robo-AO performance. Faint-target science programs primarily limited by 50% encircled-energy (e.g. those employing integral field spectrographs placed behind the AO system) may see significant benefits to sky coverage from employing laser-only AO. △ Less

Submitted 12 November, 2017; originally announced November 2017.

Comments: Accepted for publication in The Astronomical Journal. 7 pages, 6 figures

Simulating the UV Escape Fractions from Molecular Cloud Populations in Star-forming Dwarf and Spiral Galaxies

Authors: Corey S. Howard, Ralph E. Pudritz, Bill E. Harris, Ralf S. Klessen

Abstract: The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) --- Giant Molecular Clouds (GMCs) --- is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, f$_{\text{esc,GMC}}$, from a suite of individual GMC simulations with masses in the range 10$^{4-6}$ M… ▽ More The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) --- Giant Molecular Clouds (GMCs) --- is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, f$_{\text{esc,GMC}}$, from a suite of individual GMC simulations with masses in the range 10$^{4-6}$ M$_{\odot}$ using the adaptive-mesh refinement code FLASH. We find significantly different f$_{\text{esc,GMC}}$ depending on the GMC mass which can reach $>$90% in the evolution of 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ clouds or remain low at $\sim$5% for most of the lifetime of more massive GMCs. All clouds show fluctuations over short, sub-Myr timescales produced by flickering HII regions. We combine our results to calculate the total escape fraction (f$_{\text{esc,tot}}$) from GMC populations in dwarf starburst and spiral galaxies by randomly drawing clouds from a GMC mass distribution (dN/dM$\propto$M$^α$, where $α$ is either -1.5 or -2.5) over fixed time intervals. We find typical f$_{\text{esc,tot}}$ values of 8% for both the dwarf and spiral models. The fluctuations of f$_{\text{esc,tot}}$, however, are much larger for the dwarf models with values as high as 90%. The photons escaping from the 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ GMCs are the dominant contributors to f$_{\text{esc,tot}}$ in all cases. We also show that the accompanying star formation rates (SFRs) of our model ($\sim$2$\times$10$^{-2}$ and 0.73 M$_{\odot}$yr$^{-1}$) are consistent with observations of SFRs in dwarf starburst and spiral galaxies, respectively. △ Less

Submitted 18 December, 2017; v1 submitted 11 October, 2017; originally announced October 2017.

Comments: 15 pages, 7 figures, accepted to MNRAS

Ultra Violet Escape Fractions from Giant Molecular Clouds During Early Cluster Formation

Authors: Corey S. Howard, Ralph E. Pudritz, Ralf S. Klessen

Abstract: The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the Interstellar Medium (ISM), and extragalactic processes, such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction (f$_{esc}$) arising as a consequence of star cluster formation in a turbulent, 10$^6$ M$_{\odot}$ GMC, simulated using the… ▽ More The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the Interstellar Medium (ISM), and extragalactic processes, such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction (f$_{esc}$) arising as a consequence of star cluster formation in a turbulent, 10$^6$ M$_{\odot}$ GMC, simulated using the code FLASH. We make use of sink particles to represent young, star-forming clusters coupled with a radiative transfer scheme to calculate the emergent UV flux. We find that the ionizing photon flux across the cloud boundary is highly variable in time and space due to the turbulent nature of the intervening gas. The escaping photon fraction remains at $\sim$5% for the first 2.5 Myr, followed by two pronounced peaks at 3.25 and 3.8 Myr with a maximum f$_{esc}$ of 30% and 37%, respectively. These peaks are due to the formation of large HII regions, that expand into regions of lower density and some of which reach the cloud surface. However, these phases are short lived and f$_{esc}$ drops sharply as the HII regions are quenched by the central cluster passing through high-density material due to the turbulent nature of the cloud. We find an average f$_{esc}$ of 15% with factor of two variations over 1 Myr timescales. Our results suggest that assuming a single value for f$_{esc}$ from a molecular cloud is in general a poor approximation, and that the dynamical evolution of the system leads to large temporal variation. △ Less

Submitted 10 January, 2017; v1 submitted 8 November, 2016; originally announced November 2016.

Comments: ApJ, Volume 834, Issue 1, article id. 40. 10 pages, 5 figures