-
Gammapy: A Python package for gamma-ray astronomy
Authors:
Axel Donath,
Régis Terrier,
Quentin Remy,
Atreyee Sinha,
Cosimo Nigro,
Fabio Pintore,
Bruno Khélifi,
Laura Olivera-Nieto,
Jose Enrique Ruiz,
Kai Brügge,
Maximilian Linhoff,
Jose Luis Contreras,
Fabio Acero,
Arnau Aguasca-Cabot,
David Berge,
Pooja Bhattacharjee,
Johannes Buchner,
Catherine Boisson,
David Carreto Fidalgo,
Andrew Chen,
Mathieu de Bony de Lavergne,
José Vinícius de Miranda Cardoso,
Christoph Deil,
Matthias Füßling,
Stefan Funk
, et al. (20 additional authors not shown)
Abstract:
In this article, we present Gammapy, an open-source Python package for the analysis of astronomical $γ$-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from different $γ$-ray instruments for many analysis scenarios. Gammapy complies w…
▽ More
In this article, we present Gammapy, an open-source Python package for the analysis of astronomical $γ$-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from different $γ$-ray instruments for many analysis scenarios. Gammapy complies with several well-established data conventions in high-energy astrophysics, providing serialized data products that are interoperable with other software packages. Starting from event lists and instrument response functions, Gammapy provides functionalities to reduce these data by binning them in energy and sky coordinates. Several techniques for background estimation are implemented in the package to handle the residual hadronic background affecting $γ$-ray instruments. After the data are binned, the flux and morphology of one or more $γ$-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal, and spatial models. Estimation of flux points, likelihood profiles, and light curves is also supported. After describing the structure of the package, we show, using publicly available $γ$-ray data, the capabilities of Gammapy in multiple traditional and novel $γ$-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and power are displayed in a final multi-instrument example, where datasets from different instruments, at different stages of data reduction, are simultaneously fitted with an astrophysical flux model.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
Authors:
The Astropy Collaboration,
Adrian M. Price-Whelan,
Pey Lian Lim,
Nicholas Earl,
Nathaniel Starkman,
Larry Bradley,
David L. Shupe,
Aarya A. Patil,
Lia Corrales,
C. E. Brasseur,
Maximilian Nöthe,
Axel Donath,
Erik Tollerud,
Brett M. Morris,
Adam Ginsburg,
Eero Vaher,
Benjamin A. Weaver,
James Tocknell,
William Jamieson,
Marten H. van Kerkwijk,
Thomas P. Robitaille,
Bruce Merry,
Matteo Bachetti,
H. Moritz Günther,
Thomas L. Aldcroft
, et al. (111 additional authors not shown)
Abstract:
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as…
▽ More
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
Four Small Planets Buried in K2 Systems: What Can We Learn for TESS?
Authors:
Christina Hedges,
Nicholas Saunders,
Geert Barentsen,
Jeffrey L. Coughlin,
Josè Vinícius de Miranda Cardoso,
Veselin B. Kostov,
Jessie Dotson,
Ann Marie Cody
Abstract:
The Kepler, K2, and Transiting Exoplanet Survey Satellite (TESS) missions have provided a wealth of confirmed exoplanets, benefiting from a huge effort from the planet-hunting and follow-up community. With careful systematics mitigation, these missions provide precise photometric time series, which enable detection of transiting exoplanet signals. However, exoplanet hunting can be confounded by se…
▽ More
The Kepler, K2, and Transiting Exoplanet Survey Satellite (TESS) missions have provided a wealth of confirmed exoplanets, benefiting from a huge effort from the planet-hunting and follow-up community. With careful systematics mitigation, these missions provide precise photometric time series, which enable detection of transiting exoplanet signals. However, exoplanet hunting can be confounded by several factors, including instrumental noise, search biases, and host star variability. In this Letter, we discuss strategies to overcome these challenges using newly emerging techniques and tools. We demonstrate the power of new, fast open-source community tools (e.g., lightkurve, starry, celerite, exoplanet), and discuss four high signal-to-noise ratio (S/N) exoplanets that showcase specific challenges present in planet detection: K2-43c, K2-168c, K2-198c, and K2-198d. These planets have been undetected in several large K2 planet searches, despite having transit signals with S/N > 10. Two of the planets discussed here are new discoveries. In this work we confirm all four as true planets. Alongside these planet systems, we discuss three key challenges in finding small transiting exoplanets. The aim of this Letter is to help new researchers understand where planet detection efficiency gains can be made, and to encourage the continued use of K2 archive data. The considerations presented in this Letter are equally applicable to Kepler, K2, and TESS, and the tools discussed here are available for the community to apply to improve exoplanet discovery and fitting.
△ Less
Submitted 18 July, 2019;
originally announced July 2019.
-
eleanor: An open-source tool for extracting light curves from the TESS Full-Frame Images
Authors:
Adina D. Feinstein,
Benjamin T. Montet,
Daniel Foreman-Mackey,
Megan E. Bedell,
Nicholas Saunders,
Jacob L. Bean,
Jessie L. Christiansen,
Christina Hedges,
Rodrigo Luger,
Daniel Scolnic,
Jose Vinicius de Miranda Cardoso
Abstract:
During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture ph…
▽ More
During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture photometry is recorded every two minutes. However, TESS is also recording and delivering Full-Frame Images (FFIs) of each detector at a 30 minute cadence. We have created an open-source tool, eleanor, to produce light curves for objects in the TESS FFIs. Here, we describe the methods used in eleanor to produce light curves that are optimized for planet searches. The tool performs background subtraction, aperture and PSF photometry, decorrelation of instrument systematics, and cotrending using principal component analysis. We recover known transiting exoplanets in the FFIs to validate the pipeline and perform a limited search for new planet candidates in Sector 1. Our tests indicate that eleanor produces light curves with significantly less scatter than other tools that have been used in the literature. Cadence-stacked images, and raw and detrended eleanor light curves for each analyzed star will be hosted on MAST, with planet candidates on ExoFOP-TESS as Community TESS Objects of Interest (CTOIs). This work confirms the promise that the TESS FFIs will enable the detection of thousands of new exoplanets and a broad range of time domain astrophysics.
△ Less
Submitted 13 June, 2019; v1 submitted 21 March, 2019;
originally announced March 2019.
-
K2 Observations of SN 2018oh Reveal a Two-Component Rising Light Curve for a Type Ia Supernova
Authors:
G. Dimitriadis,
R. J. Foley,
A. Rest,
D. Kasen,
A. L. Piro,
A. Polin,
D. O. Jones,
A. Villar,
G. Narayan,
D. A. Coulter,
C. D. Kilpatrick,
Y. -C. Pan,
C. Rojas-Bravo,
O. D. Fox,
S. W. Jha,
P. E. Nugent,
A. G. Riess,
D. Scolnic,
M. R. Drout,
G. Barentsen,
J. Dotson,
M. Gully-Santiago,
C. Hedges,
A. M. Cody,
T. Barclay
, et al. (125 additional authors not shown)
Abstract:
We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unus…
▽ More
We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.
△ Less
Submitted 25 November, 2018;
originally announced November 2018.
-
Photometric and Spectroscopic Properties of Type Ia Supernova 2018oh with Early Excess Emission from the $Kepler$ 2 Observations
Authors:
W. Li,
X. Wang,
J. Vinkó,
J. Mo,
G. Hosseinzadeh,
D. J. Sand,
J. Zhang,
H. Lin,
T. Zhang,
L. Wang,
J. Zhang,
Z. Chen,
D. Xiang,
L. Rui,
F. Huang,
X. Li,
X. Zhang,
L. Li,
E. Baron,
J. M. Derkacy,
X. Zhao,
H. Sai,
K. Zhang,
L. Wang,
D. A. Howell
, et al. (140 additional authors not shown)
Abstract:
Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry…
▽ More
Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $Δ$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.
△ Less
Submitted 25 November, 2018;
originally announced November 2018.
-
A Discrete Set of Possible Transit Ephemerides for Two Long Period Gas Giants Orbiting HIP 41378
Authors:
Juliette C. Becker,
Andrew Vanderburg,
Joseph E. Rodriguez,
Mark Omohundro,
Fred C. Adams,
Keivan G. Stassun,
Xinyu Yao,
Joel Hartman,
Joshua Pepper,
Gaspar Bakos,
Geert Barentsen,
Thomas G. Beatty,
Waqas Bhatti,
Ashley Chontos,
Andrew Collier Cameron,
Coel Hellier,
Daniel Huber,
David James,
Rudolf B. Kuhn,
Michael B. Lund,
Don Pollacco,
Robert J. Siverd,
Daniel J. Stevens,
Jose Vinicius de Miranda Cardoso,
Richard West
Abstract:
In 2015, K2 observations of the bright (V = 8.9, K = 7.7) star HIP 41378 revealed a rich system of at least five transiting exoplanets, ranging in size from super-Earths to gas giants. The 2015 K2 observations only spanned 74.8 days, and the outer three long-period planets in the system were only detected with a single transit, so their orbital periods and transit ephemerides could not be determin…
▽ More
In 2015, K2 observations of the bright (V = 8.9, K = 7.7) star HIP 41378 revealed a rich system of at least five transiting exoplanets, ranging in size from super-Earths to gas giants. The 2015 K2 observations only spanned 74.8 days, and the outer three long-period planets in the system were only detected with a single transit, so their orbital periods and transit ephemerides could not be determined at that time. Here, we report on 50.8 days of new K2 observations of HIP 41378 from summer 2018. These data reveal additional transits of the long-period planets HIP 41378 d and HIP 41378 f, yielding a set of discrete possible orbital periods for these two planets. We identify the most probable orbital periods for these two planets using our knowledge of the planets' transit durations, the host star's properties, the system's dynamics, and data from the ground-based HATNet, KELT, and WASP transit surveys. Targeted photometric follow-up during the most probable future transit times will be able to determine the planets' orbital periods, and will enable future observations with facilities like the James Webb Space Telescope. The methods developed herein to determine the most probable orbital periods will be important for long-period planets detected by the Transiting Exoplanet Survey Satellite, where similar period ambiguities will frequently arise due to the telescope's survey strategy.
△ Less
Submitted 26 June, 2019; v1 submitted 27 September, 2018;
originally announced September 2018.
-
Seeing Double: ASASSN-18bt Exhibits a Two-Component Rise in the Early-Time K2 Light Curve
Authors:
B. J. Shappee,
T. W. -s. Holoien,
M. R. Drout,
K. Auchettl,
M. D. Stritzinger,
C. S. Kochanek,
K. Z. Stanek,
E. Shaya,
G. Narayan,
J. S. Brown,
S. Bose,
D. Bersier,
J. Brimacombe,
Ping Chen,
Subo Dong,
S. Holmbo,
B. Katz,
J. A. Munnoz,
R. L. Mutel,
R. S. Post,
J. L. Prieto,
J. Shields,
D. Tallon,
T. A. Thompson,
P. J. Vallely
, et al. (88 additional authors not shown)
Abstract:
On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the A…
▽ More
On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.
△ Less
Submitted 23 November, 2018; v1 submitted 30 July, 2018;
originally announced July 2018.
-
K2SUPERSTAMP: The release of calibrated mosaics for the {\em Kepler/K2} Mission
Authors:
Ann Marie Cody,
Geert Barentsen,
Christina Hedges,
Michael Gully-Santiago,
José Vinícius de Miranda Cardoso
Abstract:
We describe the release of a new High Level Science Product (HLSP) available at the MAST archive. The HLSP, called K2Superstamp, consists of a series of FITS images for four open star clusters observed by the K2 Mission using so-called "superstamp" pixel masks: M35, the $\sim$150 Myr old open cluster observed during K2 Campaign 0, M67, the solar-age, solar-metallicity benchmark cluster observed du…
▽ More
We describe the release of a new High Level Science Product (HLSP) available at the MAST archive. The HLSP, called K2Superstamp, consists of a series of FITS images for four open star clusters observed by the K2 Mission using so-called "superstamp" pixel masks: M35, the $\sim$150 Myr old open cluster observed during K2 Campaign 0, M67, the solar-age, solar-metallicity benchmark cluster observed during Campaign 5, Ruprecht 147, the $\sim$3 Gyr-old open cluster observed during Campaign 7, and the Lagoon Nebula (M8/NGC 6530), the high-mass star-forming region observed during Campaign 9. While the data for these regions have long been served on MAST, until now they were only available as a disconnected set of smaller Target Pixel Files (TPFs) because the spacecraft stored these observations in small chunks. As a result, these regions have hitherto been ignored by many lightcurve and planet search pipelines. With this new release, we have stitched these TPFs together into spatially contiguous FITS images (one per cadence) to make their scientific analysis easier. In addition, each image has been fit with an accurate WCS solution so that you may locate any object of interest via its right ascension and declination. We describe here the process of stitching and astrometric calibration.
△ Less
Submitted 18 February, 2018;
originally announced February 2018.