Showing 1–50 of 83 results for author: Angus, R

Discovery and Extensive Follow-Up of SN 2024ggi, a nearby type IIP supernova in NGC 3621

Authors: Ting-Wan Chen, Sheng Yang, Shubham Srivastav, Takashi J. Moriya, Stephen J. Smartt, Sofia Rest, Armin Rest, Hsing Wen Lin, Hao-Yu Miao, Yu-Chi Cheng, Amar Aryan, Chia-Yu Cheng, Morgan Fraser, Li-Ching Huang, Meng-Han Lee, Cheng-Han Lai, Yu Hsuan Liu, Aiswarya Sankar. K, Ken W. Smith, Heloise F. Stevance, Ze-Ning Wang, Joseph P. Anderson, Charlotte R. Angus, Thomas de Boer, Kenneth Chambers , et al. (23 additional authors not shown)

Abstract: We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o… ▽ More We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o- and r-band light curves show a rapid rise of 3.3 magnitudes in 13.7 hours, much faster than SN 2023ixf (another recent, nearby, and well-observed SN II). Between 13.8 and 18.8 hours after explosion SN 2024ggi became bluer, with u-g colour dropping from 0.53 to 0.15 mag. The rapid blueward evolution indicates a wind shock breakout (SBO) scenario. No hour-long brightening expected for the SBO from a bare stellar surface was detected during our observations. The classification spectrum, taken 17 hours after the SN explosion, shows flash features of high-ionization species such as Balmer lines, He I, C III, and N III. Detailed light curve modeling reveals critical insights into the properties of the circumstellar material (CSM). Our favoured model has an explosion energy of 2 x 10^51 erg, a mass-loss rate of 10^-3 solar_mass/yr (with an assumed 10 km/s wind), and a confined CSM radius of 6 x 10^14 cm. The corresponding CSM mass is 0.4 solar_mass. Comparisons with SN 2023ixf highlight that SN 2024ggi has a smaller CSM density, resulting in a faster rise and fainter UV flux. The extensive dataset and the involvement of citizen astronomers underscore that a collaborative network is essential for SBO searches, leading to more precise and comprehensive SN characterizations. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures in manuscript, 6 pages in appendix, submitted to ApJL

The Extremely Metal-Poor SN 2023ufx: A Local Analog to High-Redshift Type II Supernovae

Authors: Michael A. Tucker, Jason Hinkle, Charlotte R. Angus, Katie Auchettl, Willem B. Hoogendam, Benjamin Shappee, Christopher S. Kochanek, Chris Ashall, Thomas de Boer, Kenneth C. Chambers, Dhvanil D. Desai, Aaron Do, Michael D. Fulton, Hua Gao, Joanna Herman, Mark Huber, Chris Lidman, Chien-Cheng Lin, Thomas B. Lowe, Eugene A. Magnier, Bailey Martin, Paloma Minguez, Matt Nicholl, Miika Pursiainen, S. J. Smartt , et al. (4 additional authors not shown)

Abstract: We present extensive observations of the Type II supernova (SN II) 2023ufx which is likely the most metal-poor SN II observed to-date. It exploded in the outskirts of a low-metallicity ($Z_{\rm host} \sim 0.1~Z_\odot$) dwarf ($M_g = -13.23\pm0.15$~mag; $r_e\sim 1$~kpc) galaxy. The explosion is luminous, peaking at $M_g\approx -18.5~$mag, and shows rapid evolution. The $r$-band (pseudo-bolometric)… ▽ More We present extensive observations of the Type II supernova (SN II) 2023ufx which is likely the most metal-poor SN II observed to-date. It exploded in the outskirts of a low-metallicity ($Z_{\rm host} \sim 0.1~Z_\odot$) dwarf ($M_g = -13.23\pm0.15$~mag; $r_e\sim 1$~kpc) galaxy. The explosion is luminous, peaking at $M_g\approx -18.5~$mag, and shows rapid evolution. The $r$-band (pseudo-bolometric) light curve has a shock-cooling phase lasting 20 (17) days followed by a 19 (23)-day plateau. The entire optically-thick phase lasts only $\approx 55~$days following explosion, indicating that the red supergiant progenitor had a thinned H envelope prior to explosion. The early spectra obtained during the shock-cooling phase show no evidence for narrow emission features and limit the pre-explosion mass-loss rate to $\dot{M} \lesssim 10^{-3}~\rm M_\odot$/yr. The photospheric-phase spectra are devoid of prominent metal absorption features, indicating a progenitor metallicity of $\lesssim 0.1~Z_\odot$. The semi-nebular ($\sim 60-130~$d) spectra reveal weak Fe II, but other metal species typically observed at these phases (Ti II, Sc II, Ba II) are conspicuously absent. The late-phase optical and near-infrared spectra also reveal broad ($\approx 10^4~\rm{km}~\rm s^{-1}$) double-peaked H$α$, P$β$, and P$γ$ emission profiles suggestive of a fast outflow launched during the explosion. Outflows are typically attributed to rapidly-rotating progenitors which also prefer metal-poor environments. This is only the second SN II with $\lesssim 0.1~Z_\odot$ and both exhibit peculiar evolution, suggesting a sizable fraction of metal-poor SNe II have distinct properties compared to nearby metal-enriched SNe II. These observations lay the groundwork for modeling the metal-poor SNe II expected in the early Universe. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 17 pages, 15 figures and 3 tables in main text, an additional 5 pages, 4 figures, and 2 tables in the appendix. Submitted to ApJ, comments welcome. All data will be made publicly available upon publication

Anomaly Detection and Approximate Similarity Searches of Transients in Real-time Data Streams

Authors: P. D. Aleo, A. W. Engel, G. Narayan, C. R. Angus, K. Malanchev, K. Auchettl, V. F. Baldassare, A. Berres, T. J. L. de Boer, B. M. Boyd, K. C. Chambers, K. W. Davis, N. Esquivel, D. Farias, R. J. Foley, A. Gagliano, C. Gall, H. Gao, S. Gomez, M. Grayling, C. -C. Lin, E. A. Magnier, K. S. Mandel, T. Matheson, S. I. Raimundo , et al. (5 additional authors not shown)

Abstract: We present LAISS (Lightcurve Anomaly Identification and Similarity Search), an automated pipeline to detect anomalous astrophysical transients in real-time data streams. We deploy our anomaly detection model on the nightly ZTF Alert Stream via the ANTARES broker, identifying a manageable $\sim$1-5 candidates per night for expert vetting and coordinating follow-up observations. Our method leverages… ▽ More We present LAISS (Lightcurve Anomaly Identification and Similarity Search), an automated pipeline to detect anomalous astrophysical transients in real-time data streams. We deploy our anomaly detection model on the nightly ZTF Alert Stream via the ANTARES broker, identifying a manageable $\sim$1-5 candidates per night for expert vetting and coordinating follow-up observations. Our method leverages statistical light-curve and contextual host-galaxy features within a random forest classifier, tagging transients of rare classes (spectroscopic anomalies), of uncommon host-galaxy environments (contextual anomalies), and of peculiar or interaction-powered phenomena (behavioral anomalies). Moreover, we demonstrate the power of a low-latency ($\sim$ms) approximate similarity search method to find transient analogs with similar light-curve evolution and host-galaxy environments. We use analogs for data-driven discovery, characterization, (re-)classification, and imputation in retrospective and real-time searches. To date we have identified $\sim$50 previously known and previously missed rare transients from real-time and retrospective searches, including but not limited to: SLSNe, TDEs, SNe IIn, SNe IIb, SNe Ia-CSM, SNe Ia-91bg-like, SNe Ib, SNe Ic, SNe Ic-BL, and M31 novae. Lastly, we report the discovery of 325 total transients, all observed between 2018-2021 and absent from public catalogs ($\sim$1% of all ZTF Astronomical Transient reports to the Transient Name Server through 2021). These methods enable a systematic approach to finding the "needle in the haystack" in large-volume data streams. Because of its integration with the ANTARES broker, LAISS is built to detect exciting transients in Rubin data. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 43 pages (67 pages with Appendix), 15 figures, submitted to ApJ

zoomies: A tool to infer stellar age from vertical action in Gaia data

Authors: Sheila Sagear, Adrian M. Price-Whelan, Sarah Ballard, Yuxi, Lu, Ruth Angus, David W. Hogg

Abstract: Stellar age measurements are fundamental to understanding a wide range of astronomical processes, including galactic dynamics, stellar evolution, and planetary system formation. However, extracting age information from Main Sequence stars is complicated, with techniques often relying on age proxies in the absence of direct measurements. The Gaia data releases have enabled detailed studies of the d… ▽ More Stellar age measurements are fundamental to understanding a wide range of astronomical processes, including galactic dynamics, stellar evolution, and planetary system formation. However, extracting age information from Main Sequence stars is complicated, with techniques often relying on age proxies in the absence of direct measurements. The Gaia data releases have enabled detailed studies of the dynamical properties of stars within the Milky Way, offering new opportunities to understand the relationship between stellar age and dynamics. In this study, we leverage high-precision astrometric data from Gaia DR3 to construct a stellar age prediction model based only on stellar dynamical properties; namely, the vertical action. We calibrate two distinct, hierarchical stellar age--vertical action relations, first employing asteroseismic ages for red giant branch stars, then isochrone ages for main-sequence turn-off stars. We describe a framework called "zoomies" based on this calibration, by which we can infer ages for any star given its vertical action. This tool is open-source and intended for community use. We compare dynamical age estimates from "zoomies" with ages derived from other techniques for a sample of open clusters and main-sequence stars with asteroseismic age measurements. We also compare dynamical age estimates for stellar samples from the Kepler, K2, and TESS exoplanet transit surveys. While dynamical age relations are associated with large uncertainty, they are generally mass-independent and depend on homogeneously measured astrometric data. These age predictions are uniquely useful for large-scale demographic investigations, especially in disentangling the relationship between planet occurrence, metallicity, and age for low-mass stars. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 23 pages, 12 figures, 2 tables. Submitted to AAS Journals

Final Moments II: Observational Properties and Physical Modeling of CSM-Interacting Type II Supernovae

Authors: W. V. Jacobson-Galán, L. Dessart, K. W. Davis, C. D. Kilpatrick, R. Margutti, R. J. Foley, R. Chornock, G. Terreran, D. Hiramatsu, M. Newsome, E. Padilla Gonzalez, C. Pellegrino, D. A. Howell, A. V. Filippenko, J. P. Anderson, C. R. Angus, K. Auchettl, K. A. Bostroem, T. G. Brink, R. Cartier, D. A. Coulter, T. de Boer, M. R. Drout, N. Earl, K. Ertini , et al. (30 additional authors not shown)

Abstract: We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δt < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist… ▽ More We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δt < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist on a characteristic timescale ($t_{\rm IIn}$) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early-time IIn-like features in addition to 35 "comparison" SNe with no evidence of early-time IIn-like features, all with ultraviolet observations. The total sample consists of 50 unpublished objects with 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both $t_{\rm IIn}$ and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through matching of peak multiband absolute magnitudes, rise times, $t_{\rm IIn}$ and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: $t_{\rm IIn} \approx 3.8[\dot{M}/(0.01 \textrm{M}_{\odot} \textrm{yr}^{-1})]$ days. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 58 pages, 24 figures, submitted to ApJ. Supplementary figures available on Github (https://github.com/wynnjacobson-galan/Flash_Spectra_Sample). Data release following publication

Methods for the detection of stellar rotation periods in individual TESS sectors and results from the Prime mission

Authors: Isabel L. Colman, Ruth Angus, Trevor David, Jason Curtis, Soichiro Hattori, Yuxi Lucy Lu

Abstract: For ongoing studies of the role of rotation in stellar evolution, we require large catalogs of rotation periods for testing and refining gyrochronology. While there is a wealth of data from the Kepler and K2 missions, TESS presents both an opportunity and a challenge: despite its all-sky coverage, rotation periods remain hard to detect. We analyzed individual TESS sectors to detect short-period st… ▽ More For ongoing studies of the role of rotation in stellar evolution, we require large catalogs of rotation periods for testing and refining gyrochronology. While there is a wealth of data from the Kepler and K2 missions, TESS presents both an opportunity and a challenge: despite its all-sky coverage, rotation periods remain hard to detect. We analyzed individual TESS sectors to detect short-period stellar rotation, using only parameters measured from light curves for a robust and unbiased method of evaluating detections. We used random forest classifiers for vetting, trained on a large corpus of period measurements in KELT data from the Oelkers et al. (2018) catalog and using TESS full-frame image light curves generated by eleanor (Feinstein et al. 2019). Finally, using data from the first 26 sectors of TESS, we analyzed 432,704 2-minute cadence single-sector light curves for FGKM dwarfs. We detected 16,800 periods in individual sector light curves, covering 10,909 distinct targets, and we present a catalog of the median period for each target as measured by a Lomb-Scargle periodogram. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 23 pages, 15 figures, accepted for publication in The Astronomical Journal

In this Day and Age: An Empirical Gyrochronology Relation for Partially and Fully Convective Single Field Stars

Authors: Yuxi Lu, Ruth Angus, Daniel Foreman-Mackey, Soichiro Hattori

Abstract: Gyrochronology, the field of age-dating stars using mainly their rotation periods and masses, is ideal for inferring the ages of individual main-sequence stars. However, due to the lack of physical understanding of the complex magnetic fields in stars, gyrochronology relies heavily on empirical calibrations that require consistent and reliable stellar age measurements across a wide range of period… ▽ More Gyrochronology, the field of age-dating stars using mainly their rotation periods and masses, is ideal for inferring the ages of individual main-sequence stars. However, due to the lack of physical understanding of the complex magnetic fields in stars, gyrochronology relies heavily on empirical calibrations that require consistent and reliable stellar age measurements across a wide range of periods and masses. In this paper, we obtain a sample of consistent ages using the gyro-kinematic age-dating method, a technique to calculate the kinematics ages of stars. Using a Gaussian Process model conditioned on ages from this sample (~ 1 - 14 Gyr) and known clusters (0.67 - 3.8 Gyr), we calibrate the first empirical gyrochronology relation that is capable of inferring ages for single, main-sequence stars between 0.67 Gyr to 14 Gyr. Cross-validating and testing results suggest our model can infer cluster and asteroseismic ages with an average uncertainty of just over 1 Gyr. With this model, we obtain gyrochronology ages for ~ 100,000 stars within 1.5 kpc of the Sun with period measurements from Kepler and ZTF, and 384 unique planet host stars. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Submitted to AJ. Missing citations welcome

SN 2022oqm: A Bright and Multi-peaked Calcium-rich Transient

Authors: S. Karthik Yadavalli, V. Ashley Villar, Luca Izzo, Yossef Zenati, Ryan J. Foley, J. Craig Wheeler, Charlotte R. Angus, Dominik Bánhidi, Katie Auchettl, Barna Imre Bíró, Attila Bódi, Zsófia Bodola, Thomas de Boer, Kenneth C. Chambers, Ryan Chornock, David A. Coulter, István Csányi, Borbála Cseh, Srujan Dandu, Kyle W. Davis, Connor Braden Dickinson, Diego Farias, Joseph Farah, Christa Gall, Hua Gao , et al. (38 additional authors not shown)

Abstract: We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheri… ▽ More We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheric-phase spectra, and strong forbidden calcium emission starting 38 days after discovery. SN 2022oqm has a relatively high peak luminosity (MB = -17 mag) for (CaRTs), making it an outlier in the population. We determine that three power sources are necessary to explain the light curve (LC), with each corresponding to a distinct peak. The first peak is powered by an expanding blackbody with a power law luminosity, suggesting shock cooling by circumstellar material (CSM). Subsequent LC evolution is powered by a double radioactive decay model, consistent with two sources of photons diffusing through optically thick ejecta. From the LC, we derive an ejecta mass and 56Ni mass of ~0.6 solar masses and ~0.09 solar masses. Spectroscopic modeling suggests 0.6 solar masses of ejecta, and with well-mixed Fe-peak elements throughout. We discuss several physical origins for SN 2022oqm and find either a surprisingly massive white dwarf progenitor or a peculiar stripped envelope model could explain SN 2022oqm. A stripped envelope explosion inside a dense, hydrogen- and helium-poor CSM, akin to SNe Icn, but with a large 56Ni mass and small CSM mass could explain SN 2022oqm. Alternatively, helium detonation on an unexpectedly massive white dwarf could also explain SN 2022oqm. △ Less

Submitted 4 April, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 35 pages, 17 figures, 7 tables, Accepted for Publication in ApJ

SN 2023emq: a flash-ionised Ibn supernova with possible CIII emissio

Authors: M. Pursiainen, G. Leloudas, S. Schulze, P. Charalampopoulos, C. R. Angus, J. P. Anderson, F. Bauer, T. -W. Chen, L. Galbany, M. Gromadzki, C. P. Gutiérrez, C. Inserra, J. Lyman, T. E. Müller-Bravo, M. Nicholl, S. J. Smartt, L. Tartaglia, P. Wiseman, D. R. Young

Abstract: SN 2023emq is a fast-evolving transient initially classified as a rare Type Icn supernova (SN), interacting with a H- and He-free circumstellar medium (CSM) around maximum light. Subsequent spectroscopy revealed the unambiguous emergence of narrow He lines, confidently placing SN 2023emq in the more common Type Ibn class. Photometrically SN 2023emq has several uncommon properties regardless of its… ▽ More SN 2023emq is a fast-evolving transient initially classified as a rare Type Icn supernova (SN), interacting with a H- and He-free circumstellar medium (CSM) around maximum light. Subsequent spectroscopy revealed the unambiguous emergence of narrow He lines, confidently placing SN 2023emq in the more common Type Ibn class. Photometrically SN 2023emq has several uncommon properties regardless of its class, including its extreme initial decay (faster than > 90% of Ibn/Icn SNe) and sharp transition in the decline rate from 0.20 mag/d to 0.07 mag/d at +20 d. The bolometric light curve can be modelled as CSM interaction with 0.32M_Sun of ejecta and 0.12M_Sun of CSM, with 0.006M_Sun of nickel, as expected of fast interacting SNe. Furthermore, broad-band polarimetry at +8.7 days (P = 0.55 +/- 0.30%) is consistent with spherical symmetry. A discovery of a transitional Icn/Ibn SN would be unprecedented and would give valuable insights into the nature of mass loss suffered by the progenitor just before death, but we favour an interpretation that SN 2023emq is a type Ibn SN that exhibited flash-ionised features in the earliest spectrum, as the features are not an exact match with other SNe Icn to date. However, the feature at 5700Å, in the region of C III and N II emission, is significantly stronger in SN 2023emq than in the few other flash-ionised Type Ibn SNe, and if it is related to C III, it possibly implies a continuum of properties between the two classes. △ Less

Submitted 27 November, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted to ApJL on 22/11/2023

An abrupt change in the stellar spin-down law at the fully convective boundary

Authors: Yuxi Lu, Victor See, Louis Amard, Ruth Angus, Sean P. Matt

Abstract: The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective sta… ▽ More The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective stars exhibit very similar period-activity relations, hinting that dynamos generated by stars with and without radiative cores hold similar properties. Here, using kinematic ages, we discover an abrupt change in the stellar spin-down law across the fully convective boundary. We found that fully convective stars exhibit a higher angular momentum loss rate, corresponding to a torque that is $\sim$ 2.25 times higher for a given angular velocity than partially convective stars around the fully convective boundary. This requires a dipole field strength that is larger by a factor of $\sim$2.5, a mass loss rate that is $\sim$4.2 times larger, or some combination of both of those factors. Since stellar-wind torques depend primarily on large-scale magnetic fields and mass loss rates, both of which derive from magnetic activity, the observed abrupt change in spin-down law suggests that the dynamos of partially and fully convective stars may be fundamentally different △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: Under review at Nature Astronomy

SN 2023ixf in Messier 101: Photo-ionization of Dense, Close-in Circumstellar Material in a Nearby Type II Supernova

Authors: W. V. Jacobson-Galan, L. Dessart, R. Margutti, R. Chornock, R. J. Foley, C. D. Kilpatrick, D. O. Jones, K. Taggart, C. R. Angus, S. Bhattacharjee, L. A. Braff, D. Brethauer, A. J. Burgasser, F. Cao, C. M. Carlile, K. C. Chambers, D. A. Coulter, E. Dominguez-Ruiz, C. B. Dickinson, T. de Boer, A. Gagliano, C. Gall, H. Gao, E. L. Gates, S. Gomez , et al. (43 additional authors not shown)

Abstract: We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) l… ▽ More We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for $\sim$8 days with respect to first light, at which time Doppler broadened features from the fastest SN ejecta form, suggesting a reduction in CSM density at $r \gtrsim 10^{15}$ cm. The early-time light curve of SN2023ixf shows peak absolute magnitudes (e.g., $M_{u} = -18.6$ mag, $M_{g} = -18.4$ mag) that are $\gtrsim 2$ mag brighter than typical type II supernovae, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light curve and multi-epoch spectral models from the non-LTE radiative transfer code CMFGEN and the radiation-hydrodynamics code HERACLES suggests dense, solar-metallicity, CSM confined to $r = (0.5-1) \times 10^{15}$ cm and a progenitor mass-loss rate of $\dot{M} = 10^{-2}$ M$_{\odot}$yr$^{-1}$. For the assumed progenitor wind velocity of $v_w = 50$ km s$^{-1}$, this corresponds to enhanced mass-loss (i.e., ``super-wind'' phase) during the last $\sim$3-6 years before explosion. △ Less

Submitted 21 August, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJL

Flaring Latitudes in Ensembles of Low Mass Stars

Authors: Ekaterina Ilin, Ruth Angus, Rodrigo Luger, Brett M. Morris, Florian U. Jehn

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

Submitted 3 June, 2023; originally announced June 2023.

Comments: 15 pages, 10 figures. Accepted to MNRAS. Repository with source code: https://github.com/ekaterinailin/flare-locations-ensembles-science Repository with data: https://zenodo.org/record/7996929

Flight of the Bumblebee: the Early Excess Flux of Type Ia Supernova 2023bee revealed by $TESS$, $Swift$ and Young Supernova Experiment Observations

Authors: Qinan Wang, Armin Rest, Georgios Dimitriadis, Ryan Ridden-harper, Matthew R. Siebert, Mark Magee, Charlotte R. Angus, Katie Auchettl, Kyle W. Davis, Ryan J. Foley, Ori D. Fox, Sebastian Gomez, Jacob E. Jencson, David O. Jones, Charles D. Kilpatrick, Justin D. R. Pierel, Anthony L. Piro, Abigail Polin, Collin A. Politsch, César Rojas-bravo, Melissa Shahbandeh, V. Ashley Villar, Yossef Zenati, C. Ashall, Kenneth C. Chambers , et al. (19 additional authors not shown)

Abstract: We present high-cadence ultraviolet through near-infrared observations of the Type Ia supernova (SN Ia) 2023bee in NGC~2708 ($D = 32 \pm 3$ Mpc), finding excess flux in the first days after explosion relative to the expected power-law rise from an expanding fireball. This deviation from typical behavior for SNe Ia is particularly obvious in our 10-minute cadence $TESS$ light curve and $Swift$ UV d… ▽ More We present high-cadence ultraviolet through near-infrared observations of the Type Ia supernova (SN Ia) 2023bee in NGC~2708 ($D = 32 \pm 3$ Mpc), finding excess flux in the first days after explosion relative to the expected power-law rise from an expanding fireball. This deviation from typical behavior for SNe Ia is particularly obvious in our 10-minute cadence $TESS$ light curve and $Swift$ UV data. Compared to a few other normal SNe Ia with detected early excess flux, the excess flux in SN 2023bee is redder in the UV and less luminous. We present optical spectra of SN 2023bee, including two spectra during the period where the flux excess is dominant. At this time, the spectra are similar to those of other SNe Ia but with weaker Si II, C II and Ca II absorption lines, perhaps because the excess flux creates a stronger continuum. We compare the data to several theoretical models that have been proposed to explain the early flux excess in SNe Ia. Interaction with either a nearby companion star or close-in circumstellar material is expected to produce a faster evolution than seen in the data. Radioactive material in the outer layers of the ejecta, either from a double detonation explosion or simply an explosion with a $^{56}$Ni clump near the surface, can not fully reproduce the evolution either, likely due to the sensitivity of early UV observable to the treatment of the outer part of ejecta in simulation. We conclude that no current model can adequately explain the full set of observations. We find that a relatively large fraction of nearby, bright SNe Ia with high-cadence observations have some amount of excess flux within a few days of explosion. Considering potential asymmetric emission, the physical cause of this excess flux may be ubiquitous in normal SNe Ia. △ Less

Submitted 19 November, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 21 pages, 12 figures. Accepted by the astrophysical journal

The broad-lined Type-Ic supernova SN 2022xxf with extraordinary two-humped light curves

Authors: H. Kuncarayakti, J. Sollerman, L. Izzo, K. Maeda, S. Yang, S. Schulze, C. R. Angus, M. Aubert, K. Auchettl, M. Della Valle, L. Dessart, K. Hinds, E. Kankare, M. Kawabata, P. Lundqvist, T. Nakaoka, D. Perley, S. I. Raimundo, N. L. Strotjohann, K. Taguchi, Y. -Z. Cai, P. Charalampopoulos, Q. Fang, M. Fraser, C. P. Gutierrez , et al. (38 additional authors not shown)

Abstract: We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical… ▽ More We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical and near-infrared photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high signal-to-noise-ratio spectroscopy were obtained within 130 days, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN 2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km s$^{-1}$) are present in the spectra from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, in comparison to the typical strong narrow lines of CSM-interacting SNe, for example, Type IIn, Ibn, and Icn, but some are readily noticeable at late times such as in Mg I $λ$5170 and [O I] $λ$5577. Unusually, the near-infrared spectra show narrow line peaks in a number of features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. (Abstract abridged) △ Less

Submitted 14 August, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted version

Journal ref: A&A 678, A209 (2023)

An unlikely survivor: a low-density hot Neptune orbiting a red giant star

Authors: Samuel Grunblatt, Nicholas Saunders, Daniel Huber, Daniel Thorngren, Shreyas Vissapragada, Stephanie Yoshida, Kevin Schlaufman, Steven Giacalone, Mason MacDougall, Ashley Chontos, Emma Turtelboom, Corey Beard, Joseph M. Akana Murphy, Malena Rice, Howard Isaacson, Ruth Angus, Andrew W. Howard

Abstract: Hot Neptunes, gaseous planets smaller than Saturn ($\sim$ 3-8 R$_\oplus$) with orbital periods less than 10 days, are rare. Models predict this is due to high-energy stellar irradiation stripping planetary atmospheres over time, often leaving behind only rocky planetary cores. We present the discovery of a 6.2 R$_\oplus$(0.55 R$_\mathrm{J}$), 19.2 M$_\oplus$(0.060 M$_\mathrm{J}$) planet transiting… ▽ More Hot Neptunes, gaseous planets smaller than Saturn ($\sim$ 3-8 R$_\oplus$) with orbital periods less than 10 days, are rare. Models predict this is due to high-energy stellar irradiation stripping planetary atmospheres over time, often leaving behind only rocky planetary cores. We present the discovery of a 6.2 R$_\oplus$(0.55 R$_\mathrm{J}$), 19.2 M$_\oplus$(0.060 M$_\mathrm{J}$) planet transiting a red giant star every 4.21285 days. The old age and high equilibrium temperature yet remarkably low density of this planet suggests that its gaseous envelope should have been stripped by high-energy stellar irradiation billions of years ago. The present day planet mass and radius suggest atmospheric stripping was slower than predicted. Unexpectedly low stellar activity and/or late-stage planet inflation could be responsible for the observed properties of this system. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Comments: 54 pages, 16 figures, submitted to Science, comments welcome!

Detecting Solar System Analogs through Joint Radial Velocity/Astrometric Surveys

Authors: Daniel A. Yahalomi, Ruth Angus, David N. Spergel, Daniel Foreman-Mackey

Abstract: Earth-mass exoplanets on year-long orbits and cool gas giants (CGG) on decade-long orbits lie at the edge of current detection limits. The Terra Hunting Experiment (THE) will take nightly radial velocity (RV) observations on HARPS3 of at least 40 bright nearby G and K dwarfs for 10 years, with a target 1$σ$ measurement error of $\sim$0.3 m/s, in search of exoplanets that are Earth-like in mass and… ▽ More Earth-mass exoplanets on year-long orbits and cool gas giants (CGG) on decade-long orbits lie at the edge of current detection limits. The Terra Hunting Experiment (THE) will take nightly radial velocity (RV) observations on HARPS3 of at least 40 bright nearby G and K dwarfs for 10 years, with a target 1$σ$ measurement error of $\sim$0.3 m/s, in search of exoplanets that are Earth-like in mass and temperature. However, RV observations can only provide minimum mass estimates, due to the mass-inclination degeneracy. Astrometric observations of these same stars, with sufficient precision, could break this degeneracy. Gaia will soon release $\sim$100-200 astrometric observations of the THE stars with a 10 year baseline and $\sim$34.2 $μ$as 1$σ$ along-scan measurement error. The Nancy Grace Roman Space Telescope will be capable of precision astrometry using its wide field imager (target $\sim$5-20 $μ$as 1$σ$ measurement error for bright stars) and could extend the astrometric observational baseline to $\sim$25 years. We simulate and model an observing program that combines data from these three telescopes. We find that (1) THE RVs and Gaia astrometry can detect Earth-like and CGG-like exoplanets around bright Sun-like stars at 10 parsecs and that (2) adding Roman astrometry improves the detection precision for CGG masses and periods by a factor up to $\sim$10 and $\sim$4, respectively. Such a survey could provide insight into the prevalence of Solar System analogs, exoplanet architectures reminiscent of the mass and orbital separation hierarchy of our Solar System, for the nearest Sun-like stars. △ Less

Submitted 13 October, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 21 pages, 10 figures. Revised based on comments from anonymous reviewer at AAS Journals. Code available at https://github.com/dyahalomi/rv_and_astrometry

Evidence for the volatile-rich composition of a 1.5-$R_\oplus$ planet

Authors: Caroline Piaulet, Björn Benneke, Jose M. Almenara, Diana Dragomir, Heather A. Knutson, Daniel Thorngren, Merrin S. Peterson, Ian J. M. Crossfield, Eliza M. -R. Kempton, Daria Kubyshkina, Andrew W. Howard, Ruth Angus, Howard Isaacson, Lauren M. Weiss, Charles A. Beichman, Jonathan J. Fortney, Luca Fossati, Helmut Lammer, P. R. McCullough, Caroline V. Morley, Ian Wong

Abstract: The population of planets smaller than approximately $1.7~R_\oplus$ is widely interpreted as consisting of rocky worlds, generally referred to as super-Earths. This picture is largely corroborated by radial-velocity (RV) mass measurements for close-in super-Earths but lacks constraints at lower insolations. Here we present the results of a detailed study of the Kepler-138 system using 13 Hubble an… ▽ More The population of planets smaller than approximately $1.7~R_\oplus$ is widely interpreted as consisting of rocky worlds, generally referred to as super-Earths. This picture is largely corroborated by radial-velocity (RV) mass measurements for close-in super-Earths but lacks constraints at lower insolations. Here we present the results of a detailed study of the Kepler-138 system using 13 Hubble and Spitzer transit observations of the warm-temperate $1.51\pm0.04~R_\oplus$ planet Kepler-138 d ($T_{\mathrm{eq, A_B=0.3}}$~350 K) combined with new Keck/HIRES RV measurements of its host star. We find evidence for a volatile-rich "water world" nature of Kepler-138 d, with a large fraction of its mass contained in a thick volatile layer. This finding is independently supported by transit timing variations, RV observations ($M_d=2.1_{-0.7}^{+0.6}~M_\oplus$), as well as the flat optical/IR transmission spectrum. Quantitatively, we infer a composition of $11_{-4}^{+3}$\% volatiles by mass or ~51% by volume, with a 2000 km deep water mantle and atmosphere on top of a core with an Earth-like silicates/iron ratio. Any hypothetical hydrogen layer consistent with the observations ($<0.003~M_\oplus$) would have swiftly been lost on a ~10 Myr timescale. The bulk composition of Kepler-138 d therefore resembles those of the icy moons rather than the terrestrial planets in the solar system. We conclude that not all super-Earth-sized planets are rocky worlds, but that volatile-rich water worlds exist in an overlapping size regime, especially at lower insolations. Finally, our photodynamical analysis also reveals that Kepler-138 c ($R_c=1.51 \pm 0.04~R_\oplus$, $M_c=2.3_{-0.5}^{+0.6}~M_\oplus$) is a slightly warmer twin of Kepler-138 d, i.e., another water world in the same system, and we infer the presence of Kepler-138 e, a likely non-transiting planet at the inner edge of the habitable zone. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: Published in Nature Astronomy. 4 main figures, 10 extended data figures, 13 supplementary figures. 4 tables

The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae

Authors: P. D. Aleo, K. Malanchev, S. Sharief, D. O. Jones, G. Narayan, R. J. Foley, V. A. Villar, C. R. Angus, V. F. Baldassare, M. J. Bustamante-Rosell, D. Chatterjee, C. Cold, D. A. Coulter, K. W. Davis, S. Dhawan, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, J. Hjorth, M. E. Huber, W. V. Jacobson-Galán, C. D. Kilpatrick, D. Langeroodi , et al. (58 additional authors not shown)

Abstract: We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from… ▽ More We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 21 February, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJS; 64 pages; 35 figures; 10 tables

SN 2022ann: A type Icn supernova from a dwarf galaxy that reveals helium in its circumstellar environment

Authors: K. W. Davis, K. Taggart, S. Tinyanont, R. J. Foley, V. A. Villar, L. Izzo, C. R. Angus, M. J. Bustamante-Rosell, D. A. Coulter, N. Earl, D. Farias, J. Hjorth, M. E. Huber, D. O. Jones, P. L. Kelly, C. D. Kilpatrick, D. Langeroodi, H. -Y. Miao, C. M. Pellegrino, E. Ramirez-Ruiz, C. L. Ransome, S. Rest, S. N. Sharief, M. R. Siebert, G. Terreran , et al. (43 additional authors not shown)

Abstract: We present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of 800 km/s; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outfl… ▽ More We present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of 800 km/s; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outflowing slower than a typical Wolf-Rayet wind velocity of $>$1000 km/s. We identify helium in NIR spectra obtained two weeks after maximum and in optical spectra at three weeks, demonstrating that the CSM is not fully devoid of helium. We never detect broad spectral features from SN ejecta, including in spectra extending to the nebular phase, a unique characteristic among SNe~Icn. Compared to other SNe Icn, SN 2022ann has a low luminosity, with a peak o-band absolute magnitude of -17.7, and evolves slowly. We model the bolometric light curve and find it is well-described by 1.7 M_Sun of SN ejecta interacting with 0.2 M_sun of CSM. We place an upper limit of 0.04 M_Sun of Ni56 synthesized in the explosion. The host galaxy is a dwarf galaxy with a stellar mass of 10^7.34 M_Sun (implied metallicity of log(Z/Z_Sun) $\approx$ 0.10) and integrated star-formation rate of log(SFR) = -2.20 M_sun/yr; both lower than 97\% of the galaxies observed to produce core-collapse supernovae, although consistent with star-forming galaxies on the galaxy Main Sequence. The low CSM velocity, nickel and ejecta masses, and likely low-metallicity environment disfavour a single Wolf-Rayet progenitor star. Instead, a binary companion star is likely required to adequately strip the progenitor before explosion and produce a low-velocity outflow. The low CSM velocity may be indicative of the outer Lagrangian points in the stellar binary progenitor, rather than from the escape velocity of a single Wolf-Rayet-like massive star. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 21 pages, 13 figures, to be submitted to MNRAS

TESS Giants Transiting Giants III: An eccentric warm Jupiter supports a period-eccentricity relation for giant planets transiting evolved stars

Authors: Samuel K. Grunblatt, Nicholas Saunders, Ashley Chontos, Soichiro Hattori, Dimitri Veras, Daniel Huber, Ruth Angus, Malena Rice, Katelyn Breivik, Sarah Blunt, Steven Giacalone, Jack Lubin, Howard Isaacson, Andrew W. Howard, David R. Ciardi, Boris S. Safonov, Ivan A. Strakhov, David W. Latham, Allyson Bieryla, George R. Ricker, Jon M. Jenkins, Peter Tenenbaum, Avi Shporer, Edward H. Morgan, Veselin Kostov , et al. (5 additional authors not shown)

Abstract: The fate of planets around rapidly evolving stars is not well understood. Previous studies have suggested that relative to the main sequence population, planets transiting evolved stars ($P$ $<$ 100 d) tend to have more eccentric orbits. Here we present the discovery of TOI-4582 b, a 0.94 $\pm$ 0.12 R$_\mathrm{J}$, 0.53 $\pm$ 0.05 M$_\mathrm{J}$ planet orbiting an intermediate-mass subgiant star e… ▽ More The fate of planets around rapidly evolving stars is not well understood. Previous studies have suggested that relative to the main sequence population, planets transiting evolved stars ($P$ $<$ 100 d) tend to have more eccentric orbits. Here we present the discovery of TOI-4582 b, a 0.94 $\pm$ 0.12 R$_\mathrm{J}$, 0.53 $\pm$ 0.05 M$_\mathrm{J}$ planet orbiting an intermediate-mass subgiant star every 31.034 days. We find that this planet is also on a significantly eccentric orbit ($e$ = 0.51 $\pm$ 0.05). We then compare the population of planets found transiting evolved (log$g$ $<$ 3.8) stars to the population of planets transiting main sequence stars. We find that the rate at which median orbital eccentricity grows with period is significantly higher for evolved star systems than for otherwise similar main sequence systems, particularly for systems with only one planet detected. In general, we observe that mean planet eccentricity $<e>$ = $a$ + $b$log$_{10}$($P$) for the evolved population with a single transiting planet where $a$ = (-0.18 $\pm$ 0.08) and $b$ = (0.38 $\pm$ 0.06), significantly distinct from the main sequence planetary system population. This trend is seen even after controlling for stellar mass and metallicity. These systems do not appear to represent a steady evolution pathway from eccentric, long-period planetary orbits to circular, short period orbits, as orbital model comparisons suggest inspiral timescales are uncorrelated with orbital separation or eccentricity. Characterization of additional evolved planetary systems will distinguish effects of stellar evolution from those of stellar mass and composition. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: 21 pages, 11 figures, favorably reviewed by AAS Journals

Bridging the gap -- the disappearance of the intermediate period gap for fully convective stars, uncovered by new ZTF rotation periods

Authors: Yuxi Lu, Jason L. Curtis, Ruth Angus, Trevor J. David, Soichiro Hattori

Abstract: The intermediate period gap, discovered by Kepler, is an observed dearth of stellar rotation periods in the temperature-period diagram at $\sim$ 20 days for G dwarfs and up to $\sim$ 30 days for early-M dwarfs. However, because Kepler mainly targeted solar-like stars, there is a lack of measured periods for M dwarfs, especially those at the fully convective limit. Therefore it is unclear if the in… ▽ More The intermediate period gap, discovered by Kepler, is an observed dearth of stellar rotation periods in the temperature-period diagram at $\sim$ 20 days for G dwarfs and up to $\sim$ 30 days for early-M dwarfs. However, because Kepler mainly targeted solar-like stars, there is a lack of measured periods for M dwarfs, especially those at the fully convective limit. Therefore it is unclear if the intermediate period gap exists for mid- to late-M dwarfs. Here, we present a period catalog containing 40,553 rotation periods (9,535 periods $>$ 10 days), measured using the Zwicky Transient Facility (ZTF). To measure these periods, we developed a simple pipeline that improves directly on the ZTF archival light curves and reduces the photometric scatter by 26%, on average. This new catalog spans a range of stellar temperatures that connect samples from Kepler with MEarth, a ground-based time domain survey of bright M-dwarfs, and reveals that the intermediate period gap closes at the theoretically predicted location of the fully convective boundary ($G_{\rm BP} - G_{\rm RP} \sim 2.45$ mag). This result supports the hypothesis that the gap is caused by core-envelope interactions. Using gyro-kinematic ages, we also find a potential rapid spin-down of stars across this period gap. △ Less

Submitted 12 October, 2022; originally announced October 2022.

A fast rising tidal disruption event from a candidate intermediate mass black hole

Authors: C. R. Angus, V. F. Baldassare, B. Mockler, R. J. Foley, E. Ramirez-Ruiz, S. I. Raimundo, K. D. French, K. Auchettl, H. Pfister, C. Gall, J. Hjorth, M. R. Drout, K. D. Alexander, G. Dimitriadis, T. Hung, D. O. Jones, A. Rest, M. R. Siebert, K. Taggart, G. Terreran, S. Tinyanont, C. M. Carroll, L. DeMarchi, N. Earl, A. Gagliano , et al. (14 additional authors not shown)

Abstract: Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emit… ▽ More Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emitted from accreting gaseous discs close to the BHs. Without this light, they are difficult to detect. Tidal disruption events (TDEs), the luminous flares produced when a star strays close to a BH and is shredded, are a direct way to probe massive BHs. The rise times of these flares theoretically correlate with the BH mass. Here we present AT2020neh, a fast rising TDE candidate, hosted by a dwarf galaxy. AT2020neh can be described by the tidal disruption of a main sequence star by a 10$^{4.7} - 10^{5.9} M_{\odot}$ BH. We find the observable rate of fast rising nuclear transients like AT2020neh to be rare, at $\lesssim 2 \times 10^{-8}$ events Mpc$^{-3}$ yr$^{-1}$. Finding non-accreting BHs in dwarf galaxies is important to determine how prevalent BHs are within these galaxies, and constrain models of BH formation. AT2020neh-like events may provide a galaxy-independent method of measuring IMBH masses. △ Less

Submitted 5 September, 2022; v1 submitted 31 August, 2022; originally announced September 2022.

Comments: Accepted for publication in Nature Astronomy

wdwarfdate: A Python Package to Derive Bayesian Ages of White Dwarfs

Authors: Rocio Kiman, Siyi Xu, Jacqueline K. Faherty, Jonathan Gagne, Ruth Angus, Timothy D. Brandt, Sarah L. Casewell, Kelle L. Cruz

Abstract: White dwarfs have been successfully used as cosmochronometers in the literature, however their reach has been limited in comparison to their potential. We present wdwarfdate, a publicly available Python package to derive the Bayesian age of a white dwarf, based on its effective temperature (Teff) and surface gravity (logg). We make this software easy to use with the goal of transforming the usage… ▽ More White dwarfs have been successfully used as cosmochronometers in the literature, however their reach has been limited in comparison to their potential. We present wdwarfdate, a publicly available Python package to derive the Bayesian age of a white dwarf, based on its effective temperature (Teff) and surface gravity (logg). We make this software easy to use with the goal of transforming the usage of white dwarfs as cosmochronometers into an accessible tool. The code estimates the mass and cooling age of the white dwarf, as well as the mass and main-sequence age of the progenitor star, allowing for a determination of the total age of the object. We test the reliability of the method by estimating the parameters of white dwarfs from previous studies, and find agreement with the literature within measurement errors. By analyzing the limitation of the code we find a typical uncertainty of 10% on the total age when both input parameters have uncertainties of 1%, and an uncertainty of 25% on the total age when Teff has an uncertainty of 10% and logg of 1%. Furthermore, wdwarfdate assumes single star evolution, and can be applied to calculate the total age of a white dwarf with parameters in the range 1,500<Teff<90,000 K and 7.9<logg<9.3. Finally, the code assumes a uniform mixture of C/O in the core and single star evolution, which is reliable in the range of white dwarf masses 0.45-1.1 Msun (7.73<logg<8.8). △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 19 pages, 12 figures, 2 tables, Accepted for publication in AJ

The 3D Galactocentric velocities of Kepler stars: marginalizing over missing RVs

Authors: Ruth Angus, Adrian M. Price-Whelan, Joel C. Zinn, Megan Bedell, Yuxi, Lu, Daniel Foreman-Mackey

Abstract: Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to dire… ▽ More Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities, however one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three-quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: Accepted for publication in AAS Journals

Hubble constant and nuclear equation of state from kilonova spectro-photometric light curves

Authors: M. A. Pérez-García, L. Izzo, D. Barba, M. Bulla, A. Sagués-Carracedo, E. Pérez, C. Albertus, S. Dhawan, F. Prada, A. Agnello, C. R. Angus, S. H. Bruun, C. del Burgo, C. Dominguez-Tagle, C. Gall, A. Goobar, J. Hjorth, D. Jones, A. R. López-Sánchez, J. Sollerman

Abstract: The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected dur… ▽ More The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy $r-$process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the $10.4$ m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600-10000 Angstron spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and b) measure the Hubble constant, $H_0$, with precision improved by up to 40$\%$, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of $H_0$. △ Less

Submitted 14 December, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 21 pages, 23 figures, matches published version in A&A 666, A67 (2022)

Journal ref: A&A 666, A67 (2022)

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature-Period Distribution

Authors: Trevor J. David, Ruth Angus, Jason L. Curtis, Jennifer L. van Saders, Isabel L. Colman, Gabriella Contardo, Yuxi Lu, Joel C. Zinn

Abstract: We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature-period distribution for main-sequence stars with temperatures exceeding $\sim 5500$K. The long-period pileup is well-described by a curve of constant Rossby number, with a critical valu… ▽ More We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature-period distribution for main-sequence stars with temperatures exceeding $\sim 5500$K. The long-period pileup is well-described by a curve of constant Rossby number, with a critical value of $\mathrm{Ro_{crit}} \lesssim 2$. The long-period pileup was predicted by van Saders et al. (2019) as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages ($\sim 2-6$Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well-described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. (2020) to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well-described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. (2013a) in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a non-uniform star-formation history in the Kepler field. △ Less

Submitted 10 May, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted to ApJ. 29 pages, 21 figures. The data and code required to reproduce this work is available at http://github.com/trevordavid/rossby-ridge

The Circumstellar Environments of Double-Peaked, Calcium-strong Supernovae 2021gno and 2021inl

Authors: Wynn Jacobson-Galán, Padma Venkatraman, Raffaella Margutti, David Khatami, Giacomo Terreran, Ryan J. Foley, Rodrigo Angulo, Charlotte R. Angus, Katie Auchettl, Peter K. Blanchard, Alexey Bobrick, Joe S. Bright, Cirilla D. Couch, David A. Coulter, Karoli Clever, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Sierra A. Dodd, David O. Jones, Jessica Johnson, Charles D. Kilpatrick, Nandita Khetan, Zhisen Lai, Danial Langeroodi , et al. (20 additional authors not shown)

Abstract: We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock co… ▽ More We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock cooling emission (SCE) and/or shock interaction with circumstellar material (CSM). The primary peak in SN 2021gno is coincident with luminous, rapidly decaying X-ray emission ($L_x = 5 \times 10^{41}$ erg s$^{-1}$) detected by Swift-XRT at $δt = 1$ day after explosion, this observation being the second ever detection of X-rays from a calcium-strong transient. We interpret the X-ray emission from SN 2021gno in the context of shock interaction with dense CSM that extends to $r < 3 \times 10^{14}$ cm. Based on modeling of the SN 2021gno X-ray spectrum, we calculate a CSM mass range of $M_{\rm CSM} = (0.3 - 1.6) \times 10^{-3}$ M$_{\odot}$ and particle densities of $n = (1-4) \times 10^{10}$ cm$^{-3}$. Radio non-detections of SN 2021gno indicate a low-density environment at larger radii ($r > 10^{16}$ cm) and a progenitor mass loss rate of $\dot{M} < 10^{-4}$ M$_{\odot}$ yr$^{-1}$, for $v_w = 500$ km s$^{-1}$. For radiation derived from SCE, modeling of the primary light curve peak in both SNe indicates an extended progenitor envelope mass and radius of $M_e = 0.02 - 0.05$ M$_{\odot}$ and $R_e = 30 - 230$ R$_{\odot}$. The explosion properties of SNe 2021gno and 2021inl suggest progenitor systems containing either a low-mass massive star or a white dwarf (WD), the former being unlikely for either object given the lack of star formation at both explosion sites. Furthermore, the progenitor environments of both SNe are consistent with explosion models for low-mass hybrid He/C/O WD + C/O WD binaries. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 33 pages, 19 figures. Submitted to ApJ. Comments welcome!

Systematic errors on optical-SED stellar mass estimates for galaxies across cosmic time and their impact on cosmology

Authors: A. Paulino-Afonso, S. González-Gaitán, L. Galbany, A. M. Mourão, C. R. Angus, M. Smith, J. P. Anderson, J. D. Lyman, H. Kuncarayakti, M. A. Rodrigues

Abstract: Studying galaxies at different cosmic epochs entails several observational effects that need to be taken into account to compare populations across a large time span in a consistent manner. We use a sample of 166 nearby galaxies that hosted type Ia supernovae (SNe Ia) and have been observed with the integral field spectrograph MUSE through the AMUSING survey. Here, we present a study of the system… ▽ More Studying galaxies at different cosmic epochs entails several observational effects that need to be taken into account to compare populations across a large time span in a consistent manner. We use a sample of 166 nearby galaxies that hosted type Ia supernovae (SNe Ia) and have been observed with the integral field spectrograph MUSE through the AMUSING survey. Here, we present a study of the systematic errors and bias in the host stellar mass with increasing redshifts that are generally overlooked in SNe Ia cosmological analyses. We simulate observations at different redshifts (0.1<z<2.0) using four photometric bands (griz, similar to the Dark Energy Survey-SN program) to then estimate the host galaxy properties across cosmic time. We find that stellar masses are systematically underestimated as we move towards higher redshifts, due mostly to different rest-frame wavelength coverage, with differences reaching 0.3 dex at z~1. We have used the newly derived corrections as a function of redshift to correct the stellar masses of a known sample of SN Ia hosts and derive cosmological parameters. We show that these corrections have a small impact on the derived cosmological parameters. The most affected is the value of the mass step $Δ_M$, which is reduced by $\sim$0.004 (6% lower). The dark energy equation of state parameter $w$ changes by $Δw \sim $0.006 (0.6% higher) and the value of $Ω_m$ increases at most by 0.001 ($\sim$0.3%), all within the derived uncertainties of the model. While the systematic error found in the estimate of the host stellar mass does not significantly affect the derived cosmological parameters, it is an important source of a systematic error that one should correct for as we enter a new era of precision cosmology. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 12 pages, 9 figures, accepted for publication on A&A

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

SN 2018bsz: a Type I superluminous supernova with aspherical circumstellar material

Authors: M. Pursiainen, G. Leloudas, E. Paraskeva, A. Cikota, J. P. Anderson, C. R. Angus, S. Brennan, M. Bulla, E. Camacho-Iñiguez, P. Charalampopoulos, T. -W. Chen, M. Delgado Mancheño, M. Fraser, C. Frohmaier, L. Galbany, C. P. Gutiérrez, M. Gromadzki, C. Inserra, J. Maund, T. E. Müller-Bravo, S. Muñoz Torres, M. Nicholl, F. Onori, F. Patat, P. J. Pessi , et al. (4 additional authors not shown)

Abstract: We present a spectroscopic analysis of Type I superluminous supernova (SLSN-I), SN 2018bsz. While it closely resembles SLSNe-I, the multi-component H$α$ line appearing at $\sim30$ d post-maximum is the most atypical. The H$α$ is characterised by two emission components, one at $+3000$ km/s and a second at $-7500$ km/s, with a third, near-zero velocity component appearing after a delay. The blue an… ▽ More We present a spectroscopic analysis of Type I superluminous supernova (SLSN-I), SN 2018bsz. While it closely resembles SLSNe-I, the multi-component H$α$ line appearing at $\sim30$ d post-maximum is the most atypical. The H$α$ is characterised by two emission components, one at $+3000$ km/s and a second at $-7500$ km/s, with a third, near-zero velocity component appearing after a delay. The blue and central components can be described by Gaussian profiles of intermediate width, but the red component is significantly broader and Lorentzian. The blue component evolves towards lower velocity before fading at $100$ d post-peak, concurrently with a light curve break. Multi-component profiles are observed in other hydrogen lines including Pa$β$, and in lines of Ca II and He I. Spectropolarimetry obtained before (10.2 d) and after (38.4 d) the appearance of the H lines show a large shift on the Stokes $Q$ -- $U$ plane consistent with SN 2018bsz undergoing radical changes in its geometry. Assuming the SN is almost unpolarised at 10.2 d, the continuum polarisation at 38.4 d reaches $P \sim1.8\%$ implying a highly asymmetric configuration. We propose that the observed evolution of SN 2018bsz can be explained by highly aspherical CSM. After the SN explosion, the CSM is quickly overtaken by the ejecta, but as the photosphere starts to recede, the different CSM regions re-emerge producing the peculiar line profiles. Based on the first appearance of H$α$, we can constrain the distance of the CSM to be less than $430$ AU, or even lower ($<87$ AU) if the pre-peak plateau is related to an eruption that created the CSM. The presence of CSM has been inferred for other SLSNe-I. However, it is not clear whether the rare properties of SN 2018bsz can be generalised for SLSNe-I or whether they are the result of an uncommon evolutionary path, possibly involving a binary companion. △ Less

Submitted 29 June, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: Accepted to A&A on 22/06/2022

Journal ref: A&A 666, A30 (2022)

TESS Giants Transiting Giants II: The hottest Jupiters orbiting evolved stars

Authors: Samuel K. Grunblatt, Nicholas Saunders, Meng Sun, Ashley Chontos, Melinda Soares-Furtado, Nora Eisner, Filipe Pereira, Thaddeus Komacek, Daniel Huber, Karen Collins, Gavin Wang, Chris Stockdale, Samuel N. Quinn, Rene Tronsgaard, George Zhou, Grzegorz Nowak, Hans J. Deeg, David R. Ciardi, Andrew Boyle, Malena Rice, Fei Dai, Sarah Blunt, Judah Van Zandt, Corey Beard, Joseph M. Akana Murphy , et al. (20 additional authors not shown)

Abstract: Giant planets on short-period orbits are predicted to be inflated and eventually engulfed by their host stars. However, the detailed timescales and stages of these processes are not well known. Here we present the discovery of three hot Jupiters (P $<$ 10 d) orbiting evolved, intermediate-mass stars ($M_\star$ $\approx$ 1.5 M$_\odot$, 2 R$_\odot$ $<$ $R_\star < $ 5 R$_\odot$). By combining \tess p… ▽ More Giant planets on short-period orbits are predicted to be inflated and eventually engulfed by their host stars. However, the detailed timescales and stages of these processes are not well known. Here we present the discovery of three hot Jupiters (P $<$ 10 d) orbiting evolved, intermediate-mass stars ($M_\star$ $\approx$ 1.5 M$_\odot$, 2 R$_\odot$ $<$ $R_\star < $ 5 R$_\odot$). By combining \tess photometry with ground-based photometry and radial velocity measurements, we report masses and radii for these three planets between 0.4 and 1.8 M$_\mathrm{J}$ and 0.8 and 1.8 R$_\mathrm{J}$. \planet has the shortest period (P=\period) of any planet discovered around a red giant star to date. Both \planettwo and \planetthree appear to be inflated, but \planet does not show any sign of inflation. The large radii and relatively low masses of \planettwo and \planetthree place them among the lowest density hot Jupiters currently known, while \planet is conversely one of the highest. All three planets have orbital eccentricities below 0.2. The large spread in radii for these systems implies that planet inflation has a complex dependence on planet mass, radius, incident flux, and orbital properties. We predict that \planet has the shortest orbital decay timescale of any planet currently known, but do not detect any orbital decay in this system. Transmission spectroscopy of \planettwo would provide a favorable opportunity for the detection of water, carbon dioxide and carbon monoxide features in the atmosphere of a planet orbiting an evolved star, and could yield new information about planet formation and atmospheric evolution. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 22 pages, 15 figures, accepted for publication in the Astronomical Journal

SN 2020kyg and the rates of faint Iax Supernovae from ATLAS

Authors: Shubham Srivastav, S. J. Smartt, M. E. Huber, K. C. Chambers, C. R. Angus, T. -W. Chen, F. P. Callan, J. H. Gillanders, O. R. McBrien, S. A. Sim, M. Fulton, J. Hjorth, K. W. Smith, D. R. Young, K. Auchettl, J. P. Anderson, G. Pignata, T. J. L. de Boer, C. -C. Lin, E. A. Magnier

Abstract: We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M… ▽ More We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M$_{\odot}$ and a low kinetic energy of $E \approx 0.05 \pm 0.02 \times 10^{51}$ erg. We construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) events within 60 Mpc at $12^{+14}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+17}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for only $0.9^{+1.1}_{-0.5}\%$ of the Ia rate (a 2$σ$ upper limit of approximately 3\%). We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates. △ Less

Submitted 18 January, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted to MNRAS after minor revision

Is the high-energy neutrino event IceCube-200530A associated with a hydrogen-rich superluminous supernova?

Authors: Tetyana Pitik, Irene Tamborra, Charlotte R. Angus, Katie Auchettl

Abstract: The Zwicky Transient Facility (ZTF) follow-up campaign of alerts released by the IceCube Neutrino Observatory has led to the likely identification of the transient AT2019fdr as the source of the neutrino event IC200530A. AT2019fdr was initially suggested to be a tidal disruption event in a Narrow-Line Seyfert 1 galaxy. However, the combination of its spectral properties, color evolution, and featu… ▽ More The Zwicky Transient Facility (ZTF) follow-up campaign of alerts released by the IceCube Neutrino Observatory has led to the likely identification of the transient AT2019fdr as the source of the neutrino event IC200530A. AT2019fdr was initially suggested to be a tidal disruption event in a Narrow-Line Seyfert 1 galaxy. However, the combination of its spectral properties, color evolution, and feature-rich light curve suggests that AT2019fdr may be a Type IIn superluminous supernova. In the latter scenario, IC200530A may have been produced via inelastic proton-proton collisions between the relativistic protons accelerated at the forward shock and the cold protons of the circumstellar medium. Here, we investigate this possibility and find that at most $4.6\times 10^{-2}$ muon neutrino and antineutrino events are expected to be detected by the IceCube Neutrino Observatory within $394$ days of discovery in the case of excellent discrimination of the atmospheric background. After correcting for the Eddington bias, which occurs when a single cosmic neutrino event is adopted to infer the neutrino emission at the source, we conclude that IC200530A may originate from the hydrogen-rich superluminous supernova AT2019fdr. △ Less

Submitted 3 March, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 16 pages, including 10 figures. Improved modeling for neutrino production, conclusions unchanged, matches version accepted for publication in ApJ

Journal ref: Astrophys. J. 929 (2022) 163

$\textit{TESS}$ Giants Transiting Giants I: A Non-inflated Hot Jupiter Orbiting a Massive Subgiant

Authors: Nicholas Saunders, Samuel K. Grunblatt, Daniel Huber, Karen A. Collins, Eric L. N. Jensen, Andrew Vanderburg, Rafael Brahm, Andrés Jordán, Néstor Espinoza, Thomas Henning, Melissa J. Hobson, Samuel N. Quinn, George Zhou, R. Paul Butler, Lisa Crause, Rudi B. Kuhn, K. Moses Mogotsi, Coel Hellier, Ruth Angus, Soichiro Hattori, Ashley Chontos, George R. Ricker, Jon M. Jenkins, Peter Tenenbaum, David W. Latham , et al. (5 additional authors not shown)

Abstract: While the population of confirmed exoplanets continues to grow, the sample of confirmed transiting planets around evolved stars is still limited. We present the discovery and confirmation of a hot Jupiter orbiting TOI-2184 (TIC 176956893), a massive evolved subgiant ($M_\star= 1.53 \pm 0.12 M_\odot$, $R_\star= 2.90 \pm 0.14 R_\odot$) in the $\textit{TESS}$ Southern Continuous Viewing Zone. The pla… ▽ More While the population of confirmed exoplanets continues to grow, the sample of confirmed transiting planets around evolved stars is still limited. We present the discovery and confirmation of a hot Jupiter orbiting TOI-2184 (TIC 176956893), a massive evolved subgiant ($M_\star= 1.53 \pm 0.12 M_\odot$, $R_\star= 2.90 \pm 0.14 R_\odot$) in the $\textit{TESS}$ Southern Continuous Viewing Zone. The planet was flagged as a false positive by the $\textit{TESS}$ Quick-Look Pipeline due to periodic systematics introducing a spurious depth difference between even and odd transits. Using a new pipeline to remove background scattered light in $\textit{TESS}$ Full Frame Image (FFI) data, we combine space-based $\textit{TESS}$ photometry, ground-based photometry, and ground-based radial velocity measurements to report a planet radius of $R_p= 1.017 \pm 0.051 R_J$ and mass of $M_p= 0.65 \pm 0.16 M_J$. For a planet so close to its star, the mass and radius of TOI-2184b are unusually well matched to those of Jupiter. We find that the radius of TOI-2184b is smaller than theoretically predicted based on its mass and incident flux, providing a valuable new constraint on the timescale of post-main-sequence planet inflation. The discovery of TOI-2184b demonstrates the feasibility of detecting planets around faint ($\textit{TESS}$ magnitude $>12$) post-main sequence stars and suggests that many more similar systems are waiting to be detected in the $\textit{TESS}$ FFIs, whose confirmation may elucidate the final stages of planetary system evolution. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 20 pages, 13 figures

The unpopular Package: a Data-driven Approach to De-trend TESS Full Frame Image Light Curves

Authors: Soichiro Hattori, Daniel Foreman-Mackey, David W. Hogg, Benjamin T. Montet, Ruth Angus, T. A. Pritchard, Jason L. Curtis, Bernhard Schölkopf

Abstract: The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the c… ▽ More The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the causal pixel model method. Under the assumption that shared flux variations across multiple distant pixels are likely to be systematics, unpopular removes these common (i.e., popular) trends by modeling the systematics in a given pixel's light curve as a linear combination of light curves from many other distant pixels. To prevent overfitting we employ ridge regression and a train-and-test framework where the data points being de-trended are separated from those used to obtain the model coefficients. We also allow for simultaneous fitting with a polynomial model to capture any long-term astrophysical trends. We validate our method by de-trending different sources (e.g., supernova, tidal disruption event, exoplanet-hosting star, fast rotating star) and comparing our light curves to those obtained by other pipelines when appropriate. We also show that unpopular is able to preserve sector-length astrophysical signals, allowing for the extraction of multi-sector light curves from the FFI data. The unpopular source code and tutorials are freely available online. △ Less

Submitted 4 April, 2022; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: 35 pages, 13 figures, accepted for publication in The Astronomical Journal

The Early Phases of Supernova 2020pni: Shock-Ionization of the Nitrogen-Enriched Circumstellar Material

Authors: G. Terreran, W. V. Jacobson-Galan, J. H. Groh, R. Margutti, D. L. Coppejans, G. Dimitriadis, C. D. Kilpatrick, D. J. Matthews, M. R. Siebert, C. R. Angus, T. G. Brink, A. V. Filippenko, R. J. Foley, D. O. Jones, S. Tinyanont, C. Gall, H. Pfister, Y. Zenati, Z. Ansari, K. Auchettl, K. El-Badry, E. A. Magnier, W. Zheng

Abstract: We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a pro… ▽ More We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a progenitor mass-loss rate of $\dot{M}=(3.5-5.3)\times10^{-3}$ M$_{\odot}$ yr$^{-1}$ (assuming a wind velocity of $v_w=200\,\textrm{km}\,\textrm{s}^{-1}$), estimated at a radius of $R_{\rm in}=2.5\times10^{14}\,\rm{cm}$. In addition, we find that the progenitor of SN 2020pni was enriched in helium and nitrogen (relative abundances in mass fractions of 0.30$-$0.40, and $8.2\times10^{-3}$, respectively). Radio upper limits are also consistent with a dense CSM, and a mass-loss rate of $\dot M>5 \times 10^{-4}\,\rm{M_{\odot}\,yr^{-1}}$. During the first 4 days after first light, we also observe an increase in velocity of the hydrogen lines (from $\sim 250\,\textrm{km}\,\textrm{s}^{-1}$ to $\sim 1000\,\textrm{km}\,\textrm{s}^{-1}$), which suggests a complex CSM. The presence of dense and confined CSM, as well as its inhomogeneous structure, suggest a phase of enhanced mass loss of the progenitor of SN 2020pni during the last year before explosion. Finally, we compare SN 2020pni to a sample of other shock-photoionization events. We find no evidence of correlations among the physical parameters of the explosions and the characteristics of the CSM surrounding the progenitors of these events. This favors the idea that the mass-loss experienced by massive stars during their final years could be governed by stochastic phenomena, and that, at the same time, the physical mechanisms responsible for this mass-loss must be common to a variety of different progenitors. △ Less

Submitted 1 March, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 19 pages, 14 figures

Age-Dating Red Giant Stars Associated with Galactic Disk and Halo Substructures

Authors: Samuel K. Grunblatt, Joel C. Zinn, Adrian M. Price-Whelan, Ruth Angus, Nicholas Saunders, Marc Hon, Amalie Stokholm, Earl P. Bellinger, Sarah L. Martell, Benoit Mosser, Emily Cunningham, Jamie Tayar, Daniel Huber, Jakob Lysgaard Rørsted, Victor Silva Aguirre

Abstract: The vast majority of Milky Way stellar halo stars were likely accreted from a small number ($\lesssim$3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained, relying predominantly on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for star… ▽ More The vast majority of Milky Way stellar halo stars were likely accreted from a small number ($\lesssim$3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained, relying predominantly on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for stars associated with galactic substructures to more directly constrain the merger history of the Galaxy. By combining kinematic, asteroseismic, and spectroscopic data where available, we infer stellar ages for a sample of 10 red giant stars that were kinematically selected to be associated with the stellar halo, a subset of which are associated with the Gaia-Enceladus-Sausage halo substructure, and compare their ages to 3 red giant stars in the Galactic disk. Despite systematic differences in both absolute and relative ages determined by this work, age rankings of stars in this sample are robust. Passing the same observable inputs to multiple stellar age determination packages, we measure a weighted average age for the Gaia-Enceladus-Sausage stars in our sample of 8 $\pm$ 3 (stat.) $\pm$ 1 (sys.) Gyr. We also determine hierarchical ages for the populations of Gaia-Enceladus-Sausage, in situ halo and disk stars, finding a Gaia-Enceladus-Sausage population age of 8.0$^{+3.2}_{-2.3}$ Gyr. Although we cannot distinguish hierarchical population ages of halo or disk structures with our limited data and sample of stars, this framework should allow distinct characterization of Galactic substructures using larger stellar samples and additional data available in the near future. △ Less

Submitted 10 September, 2021; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: 28 pages, 21 figures. Published by AAS Journals

Journal ref: ApJ 916 (2021) Issue 2, id.88, 19 pp

Evaluating Rotation Periods of M Dwarfs Across the Ages

Authors: Mark Popinchalk, Jacqueline K. Faherty, Rocio Kiman, Jonathan Gagné, Jason L. Curtis, Ruth Angus, Kelle L. Cruz, Emily L. Rice

Abstract: In this work we examine M dwarf rotation rates at a range of ages to establish benchmarks for Mdwarf gyrochronology. This work includes a sample of 713 spectroscopically-classified M0-M8 dwarfs with new rotation rates measured from K2 light curves. We analyzed data and recover rotation rates for 179 of these objects. We add these to rotation rates for members of clusters with known ages (5-700 Myr… ▽ More In this work we examine M dwarf rotation rates at a range of ages to establish benchmarks for Mdwarf gyrochronology. This work includes a sample of 713 spectroscopically-classified M0-M8 dwarfs with new rotation rates measured from K2 light curves. We analyzed data and recover rotation rates for 179 of these objects. We add these to rotation rates for members of clusters with known ages (5-700 Myr), as well as objects assumed to have field ages ($>$1 Gyr). We use Gaia DR2 parallax and (G-GRP) photometry to create color-magnitude diagrams to compare objects across samples. We use color-period plots to analyze the period distributions across age, as well as incorporate Halpha equivalent width and tangential velocity where possible to further comment on age dependence. We find that the age of transition from rapid to slow rotation in clusters, which we define as an elbow in the period-color plots, depends on spectral type. Later spectral types transition at older ages: M4 for Praesepe at 700 Myr, one of the oldest clusters for which M dwarf rotation rates have been measured. The transition from active to inactive Halpha equivalent width also occurs at this elbow, as objects transition from rapid rotation to the slowly rotating sequence. Redder or smaller stars remain active at older ages. Finally, using Gaia kinematics we find evidence for rotation stalling for late Ms in the field sample, suggesting the transition happens much later for mid to late-type M dwarfs. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: 30 pages, 20 figures

Calibration of the Halpha Age-Activity relation for M dwarfs

Authors: Rocio Kiman, Jacqueline K. Faherty, Kelle L. Cruz, Jonathan Gagné, Ruth Angus, Sarah J. Schmidt, Andrew W. Mann, Daniella C. Bardalez Gagliuffi, Emily Rice

Abstract: In this work, we calibrate the relationship between Halpha emission and M dwarf ages. We compile a sample of 892 M dwarfs with Halpha equivalent width (HaEW) measurements from the literature that are either co-moving with a white dwarf of known age (21 stars) or in a known young association (871 stars). In this sample we identify 7 M dwarfs that are new candidate members of known associations. By… ▽ More In this work, we calibrate the relationship between Halpha emission and M dwarf ages. We compile a sample of 892 M dwarfs with Halpha equivalent width (HaEW) measurements from the literature that are either co-moving with a white dwarf of known age (21 stars) or in a known young association (871 stars). In this sample we identify 7 M dwarfs that are new candidate members of known associations. By dividing the stars into active and inactive categories according to their HaEW and spectral type (SpT), we find that the fraction of active dwarfs decreases with increasing age, and the form of the decline depends on SpT. Using the compiled sample of age-calibrators we find that HaEW and fractional Halpha luminosity (LHaLbol) decrease with increasing age. HaEW for SpT<M7 decreases gradually up until ~1Gyr. For older ages, we found only two early M dwarfs which are both inactive and seem to continue the gradual decrease. We also found 14 mid-type out of which 11 are inactive and present a significant decrease of HaEW, suggesting that the magnetic activity decreases rapidly after ~1Gyr. We fit LHaLbol versus age with a broken power-law and find an index of -0.11+0.02-0.01 for ages <~776Myr. The index becomes much steeper at older ages however a lack of field age-calibrators leaves this part of the relation far less constrained. Finally, from repeated independent measurements for the same stars we find that 94% of these has a level of HaEW variability <=5A at young ages (<1Gyr). △ Less

Submitted 2 April, 2021; originally announced April 2021.

Comments: 30 pages, 11 figures, 6 Tables. Accepted for publication in AJ

Gyro-Kinematic Ages for around 30,000 Kepler Stars

Authors: Yuxi, Lu, Ruth Angus, Jason L. Curtis, Trevor J. David, Rocio Kiman

Abstract: Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, a… ▽ More Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, and crucially K and M dwarfs, using their rotation and kinematic properties. We apply the simple assumption that the velocity dispersion of stars increases over time and adopt an age--velocity--dispersion relation (AVR) to estimate average stellar ages for groupings of coeval stars. We calculate the vertical velocity dispersion of stars in bins of absolute magnitude, temperature, rotation period, and Rossby number and then convert velocity dispersion to kinematic age via an AVR. Using this method, we estimate gyro-kinematic ages for 29,949 Kepler stars with measured rotation periods. We are able to estimate ages for clusters and asteroseismic stars with an RMS of 1.22 Gyr and 0.26 Gyr respectively. With our Astraea machine learning algorithm, which predicts rotation periods, we suggest a new selection criterion (a weight of 0.15) to increase the size of the McQuillian et al. (2014) catalog of Kepler rotation periods by up to 25%. Using predicted rotation periods, we estimated gyro-kinematic ages for stars without measured rotation periods and found promising results by comparing 12 detailed age--element abundance trends with literature values. △ Less

Submitted 2 February, 2021; originally announced February 2021.

A decade of radial-velocity monitoring of Vega and new limits on the presence of planets

Authors: Spencer A. Hurt, Samuel N. Quinn, David W. Latham, Andrew Vanderburg, Gilbert A. Esquerdo, Michael L. Calkins, Perry Berlind, Ruth Angus, Christian A. Latham, George Zhou

Abstract: We present an analysis of 1524 spectra of Vega spanning 10 years, in which we search for periodic radial velocity variations. A signal with a periodicity of 0.676 days and a semi-amplitude of ~10 m/s is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The… ▽ More We present an analysis of 1524 spectra of Vega spanning 10 years, in which we search for periodic radial velocity variations. A signal with a periodicity of 0.676 days and a semi-amplitude of ~10 m/s is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The timescale of evolution of these features can provide insight into the mechanism that sustains the weak magnetic fields in normal A type stars. Modeling the radial velocities with a Gaussian process using a quasi-periodic kernel suggests that the characteristic spot evolution timescale is ~180 days, though we cannot exclude the possibility that it is much longer. Such long timescales may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega's photometric rotational modulation for the first time, with a total amplitude of only 10 ppm, and a comparison of the spectroscopic and photometric amplitudes suggest the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial velocity injection recovery tests exclude the presence of transiting planets larger than 2 Earth radii and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial velocity signal with a period of 2.43 days and a semi-amplitude of 6 m/s. If caused by an orbiting companion, its minimum mass would be ~20 Earth masses; because of Vega's pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 17 pages, 8 figures, 5 tables; accepted for publication in The Astronomical Journal

Stellar Rotation in the K2 Sample: Evidence for Modified Spindown

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

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

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

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

TOI-1259Ab -- a gas giant planet with 2.7% deep transits and a bound white dwarf companion

Authors: David V. Martin, Kareem El-Badry, Vedad Kunovac Hodžić, Amaury H. M. J. Triaud, Ruth Angus, Jessica Birky, Daniel Foreman-Mackey, Christina Hedges, Benjamin T. Montet, Simon J. Murphy, Alexandre Santerne, Keivan G. Stassun, Alexander P. Stephan, Ji Wang, Paul Benni, Vadim Krushinsky, Nikita Chazov, Nikolay Mishevskiy, Carl Ziegler, Abderahmane Soubkiou, Zouhair Benkhaldoun, Douglas A. Caldwell, Karen Collins, Christopher E. Henze, Natalia M. Guerrero , et al. (11 additional authors not shown)

Abstract: We present TOI-1259Ab, a 1.0 Rjup gas giant planet transiting a 0.71 Rsun K-dwarf on a 3.48 day orbit. The system also contains a bound white dwarf companion TOI-1259B with a projected distance of approximately 1600 AU from the planet host. Transits are observed in nine TESS sector and are 2.7 per cent deep - among the deepest known - making TOI-1259Ab a promising target for atmospheric characteri… ▽ More We present TOI-1259Ab, a 1.0 Rjup gas giant planet transiting a 0.71 Rsun K-dwarf on a 3.48 day orbit. The system also contains a bound white dwarf companion TOI-1259B with a projected distance of approximately 1600 AU from the planet host. Transits are observed in nine TESS sector and are 2.7 per cent deep - among the deepest known - making TOI-1259Ab a promising target for atmospheric characterization. Our follow-up radial velocity measurements indicate a variability of semiamplitude K = 71 m/s, implying a planet mass of 0.44 Mjup. By fitting the spectral energy distribution of the white dwarf we derive a total age of 4.08 (+1.21 -0.53) Gyr for the system. The K-dwarf's light curve reveals a rotational variability with a period of 28 days, which implies a gyrochronology age broadly consistent with the white dwarf's total age. △ Less

Submitted 1 July, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

Comments: Accepted to MNRAS. Some structural changes from first arXiv version but no significant changes to results. One figure got a bit prettier

Systematics-insensitive Periodogram for finding periods in TESS observations of long-period rotators

Authors: Christina Hedges, Ruth Angus, Geert Barentsen, Nicholas Saunders, Benjamin T. Montet, Michael Gully-Santiago

Abstract: NASA's TESS mission \citep{tess} has produced high precision photometry of millions of stars to the community. The majority of TESS observations have a duration of $\approx$27 days, corresponding to a single observation during a TESS sector. A small subset of TESS targets are observed for multiple sectors, with approximately 1-2\% of targets falling in the Continuous Viewing Zone (CVZ) during the… ▽ More NASA's TESS mission \citep{tess} has produced high precision photometry of millions of stars to the community. The majority of TESS observations have a duration of $\approx$27 days, corresponding to a single observation during a TESS sector. A small subset of TESS targets are observed for multiple sectors, with approximately 1-2\% of targets falling in the Continuous Viewing Zone (CVZ) during the prime mission \citep{yield}, where targets are observed continuously for a year. These targets are highly valuable for extracting long period rotation rates, which can be linked to stellar ages. We present a pip installable Python tool for extracting long period rotation rates in the TESS CVZ, while simultaneously mitigating instrument systematics. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 3 pages, 1 figure

Journal ref: Res. Notes AAS 4 220 (2020)

WASP-107b's density is even lower: a case study for the physics of planetary gas envelope accretion and orbital migration

Authors: Caroline Piaulet, Björn Benneke, Ryan A. Rubenzahl, Andrew W. Howard, Eve J. Lee, Daniel Thorngren, Ruth Angus, Merrin Peterson, Joshua E. Schlieder, Michael Werner, Laura Kreidberg, Tareq Jaouni, Ian J. M. Crossfield, David R. Ciardi, Erik A. Petigura, John Livingston, Courtney D. Dressing, Benjamin J. Fulton, Charles Beichman, Jessie L. Christiansen, Varoujan Gorjian, Kevin K. Hardegree-Ullman, Jessica Krick, Evan Sinukoff

Abstract: With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4-year Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and… ▽ More With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4-year Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of its gas envelope. We reveal that WASP-107b's mass is only 1.8 Neptune masses ($M_b = 30.5 \pm 1.7$ $M_\oplus$). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of $> 85$% unless it is substantially inflated. The corresponding core mass of $<4.6$ $M_\oplus$ at 3$σ$ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at $\gtrsim 1$ AU onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day $a_b = 0.0566 \pm 0.0017$ AU. Beyond WASP-107b, we also detect a second more massive planet ($M_c \sin i = 0.36 \pm 0.04$ $M_{J}$) on a wide eccentric orbit ($e_c = 0.28 \pm 0.07$) which may have influenced the orbital migration and spin-orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system's formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication in AJ, 16 pages, 7 figures

Evolution of the Exoplanet Size Distribution: Forming Large Super-Earths Over Billions of Years

Authors: Trevor J. David, Gabriella Contardo, Angeli Sandoval, Ruth Angus, Yuxi, Lu, Megan Bedell, Jason L. Curtis, Daniel Foreman-Mackey, Benjamin J. Fulton, Samuel K. Grunblatt, Erik A. Petigura

Abstract: The radius valley, a bifurcation in the size distribution of small, close-in exoplanets, is hypothesized to be a signature of planetary atmospheric loss. Such an evolutionary phenomenon should depend on the age of the star-planet system. In this work, we study the temporal evolution of the radius valley using two independent determinations of host star ages among the California-Kepler Survey (CKS)… ▽ More The radius valley, a bifurcation in the size distribution of small, close-in exoplanets, is hypothesized to be a signature of planetary atmospheric loss. Such an evolutionary phenomenon should depend on the age of the star-planet system. In this work, we study the temporal evolution of the radius valley using two independent determinations of host star ages among the California-Kepler Survey (CKS) sample. We find evidence for a wide and nearly empty void of planets in the period-radius diagram at the youngest system ages ($\lesssim$2-3 Gyr) represented in the CKS sample. We show that the orbital period dependence of the radius valley among the younger CKS planets is consistent with that found among those planets with asteroseismically determined host star radii. Relative to previous studies of preferentially older planets, the radius valley determined among the younger planetary sample is shifted to smaller radii. This result is compatible with an atmospheric loss timescale on the order of gigayears for progenitors of the largest observed super-Earths. In support of this interpretation, we show that the planet sizes which appear to be unrepresented at ages $\lesssim$2-3 Gyr are likely to correspond to planets with rocky compositions. Our results suggest the size distribution of close-in exoplanets, and the precise location of the radius valley, evolves over gigayears. △ Less

Submitted 23 March, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: Accepted to The Astronomical Journal. 46 pages, 27 figures

From core collapse to superluminous: The rates of massive stellar explosions from the Palomar Transient Factory

Authors: C. Frohmaier, C. R. Angus, M. Vincenzi, M. Sullivan, M. Smith, P. E. Nugent, S. B. Cenko, A. Gal-Yam, S. R. Kulkarni, N. M. Law, R. M. Quimby

Abstract: We present measurements of the local core collapse supernova (SN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light curve realizations coupled with the detailed PTF survey detection efficiencies to forward-model the SN rates in PTF. Using a sample of 86 core collapse SNe, including 26 stripped-envelope SNe (SESNe… ▽ More We present measurements of the local core collapse supernova (SN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light curve realizations coupled with the detailed PTF survey detection efficiencies to forward-model the SN rates in PTF. Using a sample of 86 core collapse SNe, including 26 stripped-envelope SNe (SESNe), we show that the overall core collapse SN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$ at $ \langle z \rangle = 0.028$, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times10^{-5}\, \text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using 8 events at $z{\le}0.2$ of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star-formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of $\sim1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of $\sim 1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass ($\mathrm{log} M_{*} < 9.5 \mathrm{M}_\odot$) galaxies, which are the only environments that host SLSN-I in our sample, we measure a SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass and show that the specific rates of all core collapse SNe decrease with increasing stellar mass. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 17 pages, 10 figures, submitted to MNRAS

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Authors: D. O. Jones, R. J. Foley, G. Narayan, J. Hjorth, M. E. Huber, P. D. Aleo, K. D. Alexander, C. R. Angus, K. Auchettl, V. F. Baldassare, S. H. Bruun, K. C. Chambers, D. Chatterjee, D. L. Coppejans, D. A. Coulter, L. DeMarchi, G. Dimitriadis, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, T. Hung, L. Izzo, W. V. Jacobson-Galán , et al. (46 additional authors not shown)

Abstract: Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS tele… ▽ More Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date. △ Less

Submitted 5 January, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: ApJ, in press; more information at https://yse.ucsc.edu/

When Do Stalled Stars Resume Spinning Down? Advancing Gyrochronology with Ruprecht 147

Authors: Jason Lee Curtis, Marcel A. Agüeros, Sean P. Matt, Kevin R. Covey, Stephanie T. Douglas, Ruth Angus, Steven H. Saar, Ann Marie Cody, Andrew Vanderburg, Nicholas M. Law, Adam L. Kraus, David W. Latham, Christoph Baranec, Reed Riddle, Carl Ziegler, Mikkel N. Lund, Guillermo Torres, Søren Meibom, Victor Silva Aguirre, Jason T. Wright

Abstract: Recent measurements of rotation periods ($P_\text{rot}$) in the benchmark open clusters Praesepe (670 Myr), NGC 6811 (1 Gyr), and NGC 752 (1.4 Gyr) demonstrate that, after converging onto a tight sequence of slowly rotating stars in mass$-$period space, stars temporarily stop spinning down. These data also show that the duration of this epoch of stalled spin-down increases toward lower masses. To… ▽ More Recent measurements of rotation periods ($P_\text{rot}$) in the benchmark open clusters Praesepe (670 Myr), NGC 6811 (1 Gyr), and NGC 752 (1.4 Gyr) demonstrate that, after converging onto a tight sequence of slowly rotating stars in mass$-$period space, stars temporarily stop spinning down. These data also show that the duration of this epoch of stalled spin-down increases toward lower masses. To determine when stalled stars resume spinning down, we use data from the $K2$ mission and the Palomar Transient Factory to measure $P_\text{rot}$ for 58 dwarf members of the 2.7-Gyr-old cluster Ruprecht 147, 39 of which satisfy our criteria designed to remove short-period or near-equal-mass binaries. Combined with the $Kepler$ $P_\text{rot}$ data for the approximately coeval cluster NGC 6819 (30 stars with $M_\star > 0.85$ M$_\odot$), our new measurements more than double the number of $\approx$2.5 Gyr benchmark rotators and extend this sample down to $\approx$0.55 M$_\odot$. The slowly rotating sequence for this joint sample appears relatively flat (22 $\pm$ 2 days) compared to sequences for younger clusters. This sequence also intersects the $Kepler$ intermediate period gap, demonstrating that this gap was not created by a lull in star formation. We calculate the time at which stars resume spinning down, and find that 0.55 M$_\odot$ stars remain stalled for at least 1.3 Gyr. To accurately age-date low-mass stars in the field, gyrochronology formulae must be modified to account for this stalling timescale. Empirically tuning a core$-$envelope coupling model with open cluster data can account for most of the apparent stalling effect. However, alternative explanations, e.g., a temporary reduction in the magnetic braking torque, cannot yet be ruled out. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 51 pages and 21 figures. Machine-readable tables for Ruprecht 147 and the Benchmark Clusters catalogs are included in arxiv source

Astraea: Predicting Long Rotation Periods with 27-Day Light Curves

Authors: Yuxi, Lu, Ruth Angus, Marcel A. ~Agüeros, Kirsten Blancato, Melissa Ness, Jason L. ~Curtis, Sam Grunblatt

Abstract: The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for… ▽ More The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27-day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4-year light curves with 13% uncertainty overall (and a 9% uncertainty for periods > 30 days). By training on 27-day Kepler light curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods > 30 days). After training this tool on these 27-day Kepler light curve segments, we applied \texttt{Astraea} to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics. △ Less

Submitted 7 August, 2020; originally announced August 2020.

Data-driven derivation of stellar properties from photometric time series data using convolutional neural networks

Authors: Kirsten Blancato, Melissa Ness, Daniel Huber, Yuxi Lu, Ruth Angus

Abstract: Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars… ▽ More Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of \textit{Kepler} data, we recover stellar properties, including surface gravity for red giant stars (with an uncertainty of $\lesssim$ 0.06 dex), and rotation period for main sequence stars (with an uncertainty of $\lesssim$ 5.2 days, and unbiased from $\approx$5 to 40 days). Shortening the \textit{Kepler} data to a 27-day TESS-like baseline, we recover stellar properties with a small decrease in precision, $\sim$0.07 dex for log $g$ and $\sim$5.5 days for $P_{\rm rot}$, unbiased from $\approx$5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal feature engineering, and can be generalized to other surveys and datasets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 33 pages, 10 figures, 2 tables. Submitted to ApJ