-
Quasi-periodic X-ray eruptions years after a nearby tidal disruption event
Authors:
M. Nicholl,
D. R. Pasham,
A. Mummery,
M. Guolo,
K. Gendreau,
G. C. Dewangan,
E. C. Ferrara,
R. Remillard,
C. Bonnerot,
J. Chakraborty,
A. Hajela,
V. S. Dhillon,
A. F. Gillan,
J. Greenwood,
M. E. Huber,
A. Janiuk,
G. Salvesen,
S. van Velzen,
A. Aamer,
K. D. Alexander,
C. R. Angus,
Z. Arzoumanian,
K. Auchettl,
E. Berger,
T. de Boer
, et al. (39 additional authors not shown)
Abstract:
Quasi-periodic Eruptions (QPEs) are luminous bursts of soft X-rays from the nuclei of galaxies, repeating on timescales of hours to weeks. The mechanism behind these rare systems is uncertain, but most theories involve accretion disks around supermassive black holes (SMBHs), undergoing instabilities or interacting with a stellar object in a close orbit. It has been suggested that this disk could b…
▽ More
Quasi-periodic Eruptions (QPEs) are luminous bursts of soft X-rays from the nuclei of galaxies, repeating on timescales of hours to weeks. The mechanism behind these rare systems is uncertain, but most theories involve accretion disks around supermassive black holes (SMBHs), undergoing instabilities or interacting with a stellar object in a close orbit. It has been suggested that this disk could be created when the SMBH disrupts a passing star, implying that many QPEs should be preceded by observable tidal disruption events (TDEs). Two known QPE sources show long-term decays in quiescent luminosity consistent with TDEs, and two observed TDEs have exhibited X-ray flares consistent with individual eruptions. TDEs and QPEs also occur preferentially in similar galaxies. However, no confirmed repeating QPEs have been associated with a spectroscopically confirmed TDE or an optical TDE observed at peak brightness. Here we report the detection of nine X-ray QPEs with a mean recurrence time of approximately 48 hours from AT2019qiz, a nearby and extensively studied optically-selected TDE. We detect and model the X-ray, ultraviolet and optical emission from the accretion disk, and show that an orbiting body colliding with this disk provides a plausible explanation for the QPEs.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
Double "acct": a distinct double-peaked supernova matching pulsational pair-instability models
Authors:
C. R. Angus,
S. E. Woosley,
R. J. Foley,
M. Nicholl,
V. A. Villar,
K. Taggart,
M. Pursiainen,
P. Ramsden,
S. Srivastav,
H. F. Stevance,
T. Moore,
K. Auchettl,
W. B. Hoogendam,
N. Khetan,
S. K. Yadavalli,
G. Dimitriadis,
A. Gagliano,
M. R. Siebert,
A. Aamer,
T. de Boer,
K. C. Chambers,
A. Clocchiatti,
D. A. Coulter,
M. R. Drout,
D. Farias
, et al. (13 additional authors not shown)
Abstract:
We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction wit…
▽ More
We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction with hydrogen-free circumstellar material. The second peak is fainter (M$_{r}$ = -17.29 $\pm$ 0.03 mag), and spectroscopically similar to an evolved stripped-envelope SNe, with strong blended forbidden [Ca II] and [O II] features. No other known double-peak SN exhibits a light curve similar to that of SN 2020acct. We find the likelihood of two individual SNe occurring in the same star-forming region within that time to be highly improbable, while an implausibly fine-tuned configuration would be required to produce two SNe from a single binary system. We find that the peculiar properties of SN2020acct match models of pulsational pair instability (PPI), in which the initial peak is produced by collisions of shells of ejected material, shortly followed by a terminal explosion. Pulsations from a star with a 72 M$_{\odot}$ helium core provide an excellent match to the double-peaked light curve. The local galactic environment has a metallicity of 0.4 Z$_{\odot}$, a level where massive single stars are not expected retain enough mass to encounter the PPI. However, late binary mergers or a low-metallicity pocket may allow the required core mass. We measure the rate of SN 2020acct-like events to be $<3.3\times10^{-8}$ Mpc$^{-3}$ yr$^{-1}$ at z = 0.07, or <0.1% of the total core-collapse SN rate.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
SN 2021foa: The "Flip-Flop" Type IIn / Ibn supernova
Authors:
D. Farias,
C. Gall,
G. Narayan,
S. Rest,
V. A. Villar,
C. R. Angus,
K. Auchettl,
K. W. Davis,
R. Foley,
A. Gagliano,
J. Hjorth,
L. Izzo,
C. D. Kilpatrick,
H . M. L. Perkins,
E. Ramirez-Ruiz,
C. L. Ransome,
Sarangi. A.,
R. Yarza,
D. A. Coulter,
D. O. Jones,
N. Khetan,
A. Rest,
M. R. Siebert,
J. J. Swift,
K. Taggart
, et al. (7 additional authors not shown)
Abstract:
We present a comprehensive analysis of the photometric and spectroscopic evolution of SN~2021foa, unique among the class of transitional supernovae for repeatedly changing its spectroscopic appearance from hydrogen-to-helium-to-hydrogen-dominated (IIn-to-Ibn-to-IIn) within 50 days past peak brightness. The spectra exhibit multiple narrow ($\approx$ 300--600~km~s$^{-1}$) absorption lines of hydroge…
▽ More
We present a comprehensive analysis of the photometric and spectroscopic evolution of SN~2021foa, unique among the class of transitional supernovae for repeatedly changing its spectroscopic appearance from hydrogen-to-helium-to-hydrogen-dominated (IIn-to-Ibn-to-IIn) within 50 days past peak brightness. The spectra exhibit multiple narrow ($\approx$ 300--600~km~s$^{-1}$) absorption lines of hydrogen, helium, calcium and iron together with broad helium emission lines with a full-width-at-half-maximum (FWHM) of $\sim 6000$~km~s$^{-1}$. For a steady, wind-mass loss regime, light curve modeling results in an ejecta mass of $\sim 8$ M$_{\odot}$ and CSM mass below 1 M$_{\odot}$, and an ejecta velocity consistent with the FWHM of the broad helium lines. We obtain a mass-loss rate of $\approx 2$ M$_{\odot} {\rm yr}^{-1}$. This mass-loss rate is three orders of magnitude larger than derived for normal Type II SNe. We estimate that the bulk of the CSM of SN~2021foa must have been expelled within half a year, about 15 years ago. Our analysis suggests that SN~2021foa had a helium rich ejecta which swept up a dense shell of hydrogen rich CSM shortly after explosion. At about 60 days past peak brightness, the photosphere recedes through the dense ejecta-CSM region, occulting much of the red-shifted emission of the hydrogen and helium lines, which results in observed blue-shift ($\sim -3000$~km~s$^{-1}$). Strong mass loss activity prior to explosion, such as those seen in SN~2009ip-like objects and SN~2021foa as precursor emission, are the likely origin of a complex, multiple-shell CSM close to the progenitor star.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Moment-preserving Monte-Carlo Coulomb collision method for particle codes
Authors:
Justin Ray Angus,
Yichen Fu,
Vasily Geyko,
Dave Grote,
David Larson
Abstract:
Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weigh…
▽ More
Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. The efficacy of the model is illustrated with various test problems.
△ Less
Submitted 26 July, 2024;
originally announced July 2024.
-
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.
-
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.
-
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,
D. O. Jones,
C. -C. Lin,
E. A. Magnier,
K. S. Mandel,
T. Matheson
, et al. (7 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 24 July, 2024; v1 submitted 1 April, 2024;
originally announced April 2024.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
$\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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Simulation of FuZE axisymmetric stability using gyrokinetic and extended-MHD models
Authors:
V. I. Geyko,
J. R. Angus,
M. A. Dorf
Abstract:
Axisymmetric ($m=0$) gyrokinetic and extended-MHD simulations of sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The electromagnetic MHD model include…
▽ More
Axisymmetric ($m=0$) gyrokinetic and extended-MHD simulations of sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The electromagnetic MHD model includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. A prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the FuZE experiment and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center ($E{\times} B$) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease of the linear growth rate with an increase in the flow shear, however full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.
△ Less
Submitted 5 November, 2020;
originally announced November 2020.