Showing 1–50 of 69 results for author: Huppenkothen, D

A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437$\unicode{x2013}$4715

Authors: Devarshi Choudhury, Tuomo Salmi, Serena Vinciguerra, Thomas E. Riley, Yves Kini, Anna L. Watts, Bas Dorsman, Slavko Bogdanov, Sebastien Guillot, Paul S. Ray, Daniel J. Reardon, Ronald A. Remillard, Anna V. Bilous, Daniela Huppenkothen, James M. Lattimer, Nathan Rutherford, Zaven Arzoumanian, Keith C. Gendreau, Sharon M. Morsink, Wynn C. G. Ho

Abstract: We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Pa… ▽ More We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Parkes Pulsar Timing Array (PPTA) (Reardon et al. 2024), and we use NICER background models to constrain the non-source background, cross-checking with data from XMM-Newton. We assume two distinct hot emitting regions, and various parameterized hot region geometries that are defined in terms of overlapping circles; while simplified, these capture many of the possibilities suggested by detailed modeling of return current heating. For the preferred model identified by our analysis we infer a mass of $M = 1.418 \pm 0.037$ M$_\odot$ (largely informed by the PPTA mass prior) and an equatorial radius of $R = 11.36^{+0.95}_{-0.63}$ km, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. This radius favors softer dense matter equations of state and is highly consistent with constraints derived from gravitational wave measurements of neutron star binary mergers. The hot regions are inferred to be non-antipodal, and hence inconsistent with a pure centered dipole magnetic field. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 27 pages, 16 figures, + appendices. Accepted in ApJL. Data files and reproduction package available in Zenodo

The Radius of the High Mass Pulsar PSR J0740+6620 With 3.6 Years of NICER Data

Authors: Tuomo Salmi, Devarshi Choudhury, Yves Kini, Thomas E. Riley, Serena Vinciguerra, Anna L. Watts, Michael T. Wolff, Zaven Arzoumanian, Slavko Bogdanov, Deepto Chakrabarty, Keith Gendreau, Sebastien Guillot, Wynn C. G. Ho, Daniela Huppenkothen, Renee M. Ludlam, Sharon M. Morsink, Paul S. Ray

Abstract: We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gra… ▽ More We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gravitational mass are $12.49_{-0.88}^{+1.28}$ km and $2.073_{-0.069}^{+0.069}$ $M_\odot$ respectively, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles, with an estimated systematic error $\lesssim 0.1$ km. This result was obtained using the best computationally feasible sampler settings providing a strong radius lower limit but a slightly more uncertain radius upper limit. The inferred radius interval is also close to the $R=12.76_{-1.02}^{+1.49}$ km obtained by Dittmann et al. 2024, when they require the radius to be less than $16$ km as we do. The results continue to disfavor very soft equations of state for dense matter, with $R<11.15$ km for this high mass pulsar excluded at the $95\,\%$ probability. The results do not depend significantly on the assumed cross-calibration uncertainty between NICER and XMM-Newton. Using simulated data that resemble the actual observations, we also show that our pipeline is capable of recovering parameters for the inferred models reported in this paper. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures (2 of which are figure sets), 2 tables, accepted for publication in ApJ

FlopPITy: Enabling self-consistent exoplanet atmospheric retrievals with machine learning

Authors: Francisco Ardévol Martínez, Michiel Min, Daniela Huppenkothen, Inga Kamp, Paul I. Palmer

Abstract: Interpreting the observations of exoplanet atmospheres to constrain physical and chemical properties is typically done using Bayesian retrieval techniques. Because these methods require many model computations, a compromise is made between model complexity and run time. Reaching this compromise leads to the simplification of many physical and chemical processes (e.g. parameterised temperature stru… ▽ More Interpreting the observations of exoplanet atmospheres to constrain physical and chemical properties is typically done using Bayesian retrieval techniques. Because these methods require many model computations, a compromise is made between model complexity and run time. Reaching this compromise leads to the simplification of many physical and chemical processes (e.g. parameterised temperature structure). Here we implement and test sequential neural posterior estimation (SNPE), a machine learning inference algorithm, for exoplanet atmospheric retrievals. The goal is to speed up retrievals so they can be run with more computationally expensive atmospheric models, such as those computing the temperature structure using radiative transfer. We generate 100 synthetic observations using ARCiS (ARtful Modeling Code for exoplanet Science, an atmospheric modelling code with the flexibility to compute models in varying degrees of complexity) and perform retrievals on them to test the faithfulness of the SNPE posteriors. The faithfulness quantifies whether the posteriors contain the ground truth as often as we expect. We also generate a synthetic observation of a cool brown dwarf using the self-consistent capabilities of ARCiS and run a retrieval with self-consistent models to showcase the possibilities that SNPE opens. We find that SNPE provides faithful posteriors and is therefore a reliable tool for exoplanet atmospheric retrievals. We are able to run a self-consistent retrieval of a synthetic brown dwarf spectrum using only 50,000 forward model evaluations. We find that SNPE can speed up retrievals between $\sim2\times$ and $\geq10\times$ depending on the computational load of the forward model, the dimensionality of the observation, and the signal-to-noise ratio of the observation. We make the code publicly available for the community on Github. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted for publication at A&A

The Sonified Hertzsprung-Russell Diagram

Authors: Daniela Huppenkothen, Juan Pampin, James R. A. Davenport, James Wenlock

Abstract: Understanding the physical properties of stars, and putting these properties into the context of stellar evolution, is a core challenge in astronomical research. A key visualization in studying stellar evolution is the Hertzsprung-Russell diagram (HRD), organizing data about stellar luminosity and colour into a form that is informative about stellar structure and evolution. However, connecting the… ▽ More Understanding the physical properties of stars, and putting these properties into the context of stellar evolution, is a core challenge in astronomical research. A key visualization in studying stellar evolution is the Hertzsprung-Russell diagram (HRD), organizing data about stellar luminosity and colour into a form that is informative about stellar structure and evolution. However, connecting the HRD with other sources of information, including stellar time series, is an outstanding challenge. Here we present a new method to turn stellar time series into sound. This method encodes physically meaningful features such that auditory comparisons between sonifications of different stars preserve astrophysical differences between them. We present an interactive multimedia version of the HRD that combines both visual and auditory components and that allows exploration of different types of stars both on and off the main sequence through both visual and auditory media. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Comments: 8 pages, 5 figures; accepted for publication in the proceedings of "The 28th International Conference on Auditory Display (ICAD 2023) - Special Session on Astronomical Data Sonification"

Discovery of spin-phase dependent QPOs in the super-critical accretion regime from the X-ray pulsar RX J0440.9+4431

Authors: Christian Malacaria, Daniela Huppenkothen, Oliver J. Roberts, Lorenzo Ducci, Enrico Bozzo, Peter Jenke, Colleen A. Wilson-Hodge, Maurizio Falanga

Abstract: RX J0440.9+4431 is an accreting X-ray pulsar (XRP) that remained relatively unexplored until recently, when major X-ray outburst activity enabled more in-depth studies. Here, we report on the discovery of ${\sim}0.2$ Hz quasi-periodic oscillations (QPOs) from this source observed with $Fermi$-GBM. The appearance of QPOs in RX J0440.9+4431 is thricely transient, that is, QPOs appear only above a ce… ▽ More RX J0440.9+4431 is an accreting X-ray pulsar (XRP) that remained relatively unexplored until recently, when major X-ray outburst activity enabled more in-depth studies. Here, we report on the discovery of ${\sim}0.2$ Hz quasi-periodic oscillations (QPOs) from this source observed with $Fermi$-GBM. The appearance of QPOs in RX J0440.9+4431 is thricely transient, that is, QPOs appear only above a certain luminosity, only at certain pulse phases (namely corresponding to the peak of its sine-like pulse profile), and only for a few oscillations at time. We argue that this newly discovered phenomenon (appearance of thricely transient QPOs -- or ATTO) occurs if QPOs are fed through an accretion disk whose inner region viscosity is unstable to mass accretion rate and temperature variations. Such variations are triggered when the source switches to the super-critical accretion regime and the emission pattern changes. We also argue that the emission region configuration is likely responsible for the observed QPOs spin-phase dependence. △ Less

Submitted 14 November, 2023; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: 5 + 2 appendix pages. Accepted on A&A. Proofs version

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Authors: D. Huppenkothen, M. Ntampaka, M. Ho, M. Fouesneau, B. Nord, J. E. G. Peek, M. Walmsley, J. F. Wu, C. Avestruz, T. Buck, M. Brescia, D. P. Finkbeiner, A. D. Goulding, T. Kacprzak, P. Melchior, M. Pasquato, N. Ramachandra, Y. -S. Ting, G. van de Ven, S. Villar, V. A. Villar, E. Zinger

Abstract: Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr… ▽ More Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Society

Dense Forests of Microshots in Bursts from FRB 20220912A

Authors: Danté M. Hewitt, Jason W. T. Hessels, Omar S. Ould-Boukattine, Pragya Chawla, Ismaël Cognard, Akshatha Gopinath, Lucas Guillemot, Daniela Huppenkothen, Kenzie Nimmo6, Mark P. Snelders

Abstract: We report on exceptionally bright bursts (>400 Jy ms) detected from the repeating fast radio burst source FRB 20220912A using the Nançay Radio Telescope (NRT), as part of the ECLAT (Extragalactic Coherent Light from Astrophysical Transients) monitoring campaign. These bursts exhibit extremely luminous, broadband, short-duration structures (~ 16 microseconds), which we term 'microshots' and which c… ▽ More We report on exceptionally bright bursts (>400 Jy ms) detected from the repeating fast radio burst source FRB 20220912A using the Nançay Radio Telescope (NRT), as part of the ECLAT (Extragalactic Coherent Light from Astrophysical Transients) monitoring campaign. These bursts exhibit extremely luminous, broadband, short-duration structures (~ 16 microseconds), which we term 'microshots' and which can be especially well studied in the NRT data given the excellent signal-to-noise and dynamic range (32-bit samples). The estimated peak flux density of the brightest microshot is 450 Jy. We show that the microshots are clustered into dense 'forests', by modelling them as Weibull distributions and obtaining Weibull shape parameters of approximately 0.5. Our polarimetric analysis reveals that the bursts are nearly 100% linearly polarised; have < 10% circular polarisation fractions; a near-zero average rotation measure of 0.10(6) rad/m^2; and varying polarisation position angles over the burst duration. For one of the bursts, we analyse raw voltage data from simultaneous observations with the Westerbork RT-1 single 25-m dish. These data allow us to measure the scintillation bandwidth, 0.30(3) MHz, and to probe the bursts on (sub-)microsecond timescales. Some important nuances related to dedispersion are also discussed. We propose that the emission mechanism for the broadband microshots is potentially different from the emission mechanism of the broader burst components which still show a residual drift of a few hundred MHz/ms after correcting for dispersion using the microshots. We discuss how the observed emission is phenomenologically analogous to different types of radio bursts from the Sun. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: submitted to MNRAS, 19 pages, 13 figures

An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451

Authors: Serena Vinciguerra, Tuomo Salmi, Anna L. Watts, Devarshi Choudhury, Thomas E. Riley, Paul S. Ray, Slavko Bogdanov, Yves Kini, Sebastien Guillot, Deepto Chakrabarty, Wynn C. G. Ho, Daniela Huppenkothen, Sharon M. Morsink, Zorawar Wadiasingh

Abstract: In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposi… ▽ More In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyse, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and jointly analyse also XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multi-modal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of $\sim [1.4 \mathrm{M_\odot}, 11.5$ km] and $\sim [1.7 \mathrm{M_\odot}, 14.5$ km], depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 35 pages, 11 figures, 7 tables

Journal ref: ApJ, 961, 62 (2024)

Swift Deep Galactic Plane Survey Classification of Swift J170800$-$402551.8 as a Candidate Intermediate Polar Cataclysmic Variable

Authors: B. O'Connor, E. Gogus, J. Hare, K. Mukai, D. Huppenkothen, J. Brink, D. A. H. Buckley, A. Levan, M. G. Baring, R. Stewart, C. Kouveliotou, P. Woudt, E. Bellm, S. B. Cenko, P. A. Evans, J. Granot, C. Hailey, F. Harrison, D. Hartmann, A. J. van der Horst, L. Kaper, J. A. Kennea, S. B. Potter, P. O. Slane, D. Stern , et al. (2 additional authors not shown)

Abstract: Here, we present the results of our multi-wavelength campaign aimed at classifying \textit{Swift} J170800$-$402551.8 as part of the \textit{Swift} Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with \textit{Swift}, \textit{NICER}, \textit{XMM-Newton}, \textit{NuSTAR}, and the Southern African Large Telescope (SALT), as well as multi-wavelength archival obse… ▽ More Here, we present the results of our multi-wavelength campaign aimed at classifying \textit{Swift} J170800$-$402551.8 as part of the \textit{Swift} Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with \textit{Swift}, \textit{NICER}, \textit{XMM-Newton}, \textit{NuSTAR}, and the Southern African Large Telescope (SALT), as well as multi-wavelength archival observations from \textit{Gaia}, VPHAS, and VVV. The source displays a periodicity of 784 s in our \textit{XMM-Newton} observation. The X-ray spectrum (\textit{XMM-Newton} and \textit{NuSTAR}) can be described by thermal bremsstrahlung radiation with a temperature of $kT$\,$\approx$\,$30$ keV. The phase-folded X-ray lightcurve displays a double-peaked, energy-dependent pulse-profile. We used \textit{Chandra} to precisely localize the source, allowing us to identify and study the multi-wavelength counterpart. Spectroscopy with SALT identified a Balmer H$α$ line, and potential HeI lines, from the optical counterpart. The faintness of the counterpart ($r$\,$\approx$\,$21$ AB mag) favors a low-mass donor star. Based on these criteria, we classify \textit{Swift} J170800$-$402551.8 as a candidate intermediate polar cataclysmic variable, where the spin period of the white dwarf is 784 s. △ Less

Submitted 28 August, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: Accepted in MNRAS

The Swift Deep Galactic Plane Survey (DGPS) Phase-I Catalog

Authors: B. O'Connor, C. Kouveliotou, P. A. Evans, N. Gorgone, A. J. van Kooten, S. Gagnon, H. Yang, M. G. Baring, E. Bellm, P. Beniamini, J. Brink, D. A. H. Buckley, S. B. Cenko, O. D. Egbo, E. Gogus, J. Granot, C. Hailey, J. Hare, F. Harrison, D. Hartmann, A. J. van der Horst, D. Huppenkothen, L. Kaper, O. Kargaltsev, J. A. Kennea , et al. (8 additional authors not shown)

Abstract: The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March… ▽ More The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March 2017 and May 2021. The Survey is complete to depth $L_X$\,$>$\,$10^{34}$ erg s$^{-1}$ to the edge of the Galaxy. The main Survey goal is to produce a rich sample of new X-ray sources and transients, while also covering a broad discovery space. Here, we introduce the Survey strategy and present a catalog of sources detected during Phase-I observations. In total, we identify 928 X-ray sources, of which 348 are unique to our X-ray catalog. We report on the characteristics of sources in our catalog and highlight sources newly classified and published by the DGPS team. △ Less

Submitted 21 September, 2023; v1 submitted 25 June, 2023; originally announced June 2023.

Comments: Accepted in ApJS. This is the final version

Quasi-Periodic Peak Energy Oscillations in X-ray Bursts from SGR J1935+2154

Authors: Oliver J. Roberts, Matthew G. Baring, Daniela Huppenkothen, Ersin Gogus, Yuki Kaneko, Chryssa Kouveliotou, Lin Lin, Alexander J. van der Horst, George Younes

Abstract: Magnetars are young neutron stars powered by the strongest magnetic fields in the Universe (10$^{13-15}$ G). Their transient X-ray emission usually manifests as short (a few hundred milliseconds), bright, energetic ($\sim$ 10$^{40-41}$ erg) X-ray bursts. Since its discovery in 2014, magnetar J1935+2154 has become one of the most prolific magnetars, exhibiting very active bursting episodes, and oth… ▽ More Magnetars are young neutron stars powered by the strongest magnetic fields in the Universe (10$^{13-15}$ G). Their transient X-ray emission usually manifests as short (a few hundred milliseconds), bright, energetic ($\sim$ 10$^{40-41}$ erg) X-ray bursts. Since its discovery in 2014, magnetar J1935+2154 has become one of the most prolific magnetars, exhibiting very active bursting episodes, and other fascinating events such as pulse timing anti-glitches and Fast Radio Bursts. Here, we present evidence for possible 42 Hz (24 ms) quasi-periodic oscillations in the $νF_ν$ spectrum peak energy (Ep) identified in a unique burst detected with the Fermi Gamma-ray Burst Monitor in January 2022. While quasi-periodic oscillations have been previously reported in the intensity of magnetar burst lightcurves, quasi-periodic oscillations in the Ep have not. We also find an additional event from the same outburst that appears to exhibit similar character in Ep, albeit of lower statistical quality. For these two exceptional transients, such Ep oscillations can be explained by magnetospheric density and pressure perturbations. For burst-emitting plasma consisting purely of $e^+e^-$ pairs, these acoustic modes propagate along a highly magnetized flux tube of length up to around $L\sim 130$ neutron star radii, with $L$ being lower if ions are present in the emission zone. Detailed time-resolved analyses of other magnetar bursts are encouraged to evaluate the rarity of these events and their underlying mechanisms. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Journal ref: ApJL 956 L27 (2023)

Identification of 1RXS J165424.6-433758 as a polar cataclysmic variable

Authors: B. O'Connor, J. Brink, D. A. H. Buckley, K. Mukai, C. Kouveliotou, E. Gogus, S. B. Potter, P. Woudt, A. Lien, A. Levan, O. Kargaltsev, M. G. Baring, E. Bellm, S. B. Cenko, P. A. Evans, J. Granot, C. Hailey, F. Harrison, D. Hartmann, A. J. van der Horst, D. Huppenkothen, L. Kaper, J. A. Kennea, P. O. Slane, D. Stern , et al. (3 additional authors not shown)

Abstract: We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity… ▽ More We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity $L_X=(6.5\pm0.8)\times10^{31}$ erg s$^{-1}$ at a \textit{Gaia} distance of 460 pc. Spectroscopy with the Southern African Large Telescope (SALT) revealed a flat continuum dominated by emission features, demonstrating an inverse Balmer decrement, the $\lambda4640$ Bowen blend, almost a dozen HeI lines, and HeII $\lambda4541$, $\lambda4686$ and $λ5411$. Our high-speed photometry demonstrates a preponderance of flickering and flaring episodes, and revealed the orbital period of the system, $P_\textrm{orb}=2.87$ hr, which fell well within the cataclysmic variable (CV) period gap between $2-3$ hr. These features classify 1RXS J165424.6-433758 as a nearby polar magnetic CV. △ Less

Submitted 10 September, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

Fermi-GBM Discovery of GRB 221009A: An Extraordinarily Bright GRB from Onset to Afterglow

Authors: S. Lesage, P. Veres, M. S. Briggs, A. Goldstein, D. Kocevski, E. Burns, C. A. Wilson-Hodge, P. N. Bhat, D. Huppenkothen, C. L. Fryer, R. Hamburg, J. Racusin, E. Bissaldi, W. H. Cleveland, S. Dalessi, C. Fletcher, M. M. Giles, B. A. Hristov, C. M. Hui, B. Mailyan, C. Malacaria, S. Poolakkil, O. J. Roberts, A. von Kienlin, J. Wood , et al. (115 additional authors not shown)

Abstract: We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing ana… ▽ More We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{γ,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{γ,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow. △ Less

Submitted 12 July, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 26 pages 7 figures - accepted for publication in ApJL

Practical Guidance for Bayesian Inference in Astronomy

Authors: Gwendolyn M. Eadie, Joshua S. Speagle, Jessi Cisewski-Kehe, Daniel Foreman-Mackey, Daniela Huppenkothen, David E. Jones, Aaron Springford, Hyungsuk Tak

Abstract: In the last two decades, Bayesian inference has become commonplace in astronomy. At the same time, the choice of algorithms, terminology, notation, and interpretation of Bayesian inference varies from one sub-field of astronomy to the next, which can lead to confusion to both those learning and those familiar with Bayesian statistics. Moreover, the choice varies between the astronomy and statistic… ▽ More In the last two decades, Bayesian inference has become commonplace in astronomy. At the same time, the choice of algorithms, terminology, notation, and interpretation of Bayesian inference varies from one sub-field of astronomy to the next, which can lead to confusion to both those learning and those familiar with Bayesian statistics. Moreover, the choice varies between the astronomy and statistics literature, too. In this paper, our goal is two-fold: (1) provide a reference that consolidates and clarifies terminology and notation across disciplines, and (2) outline practical guidance for Bayesian inference in astronomy. Highlighting both the astronomy and statistics literature, we cover topics such as notation, specification of the likelihood and prior distributions, inference using the posterior distribution, and posterior predictive checking. It is not our intention to introduce the entire field of Bayesian data analysis -- rather, we present a series of useful practices for astronomers who already have an understanding of the Bayesian "nuts and bolts" and wish to increase their expertise and extend their knowledge. Moreover, as the field of astrostatistics and astroinformatics continues to grow, we hope this paper will serve as both a helpful reference and as a jumping off point for deeper dives into the statistics and astrostatistics literature. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 10 pages, 4 figures

Mapping the X-ray variability of GRS1915+105 with machine learning

Authors: Benjamin J. Ricketts, James F. Steiner, Cecilia Garraffo, Ronald A. Remillard, Daniela Huppenkothen

Abstract: Black hole X-ray binary systems (BHBs) contain a close companion star accreting onto a stellar-mass black hole. A typical BHB undergoes transient outbursts during which it exhibits a sequence of long-lived spectral states, each of which is relatively stable. GRS 1915+105 is a unique BHB that exhibits an unequaled number and variety of distinct variability patterns in X-rays. Many of these patterns… ▽ More Black hole X-ray binary systems (BHBs) contain a close companion star accreting onto a stellar-mass black hole. A typical BHB undergoes transient outbursts during which it exhibits a sequence of long-lived spectral states, each of which is relatively stable. GRS 1915+105 is a unique BHB that exhibits an unequaled number and variety of distinct variability patterns in X-rays. Many of these patterns contain unusual behaviour not seen in other sources. These variability patterns have been sorted into different classes based on count rate and color characteristics by Belloni et al (2000). In order to remove human decision-making from the pattern-recognition process, we employ an unsupervised machine learning algorithm called an auto-encoder to learn what classifications are naturally distinct by allowing the algorithm to cluster observations. We focus on observations taken by the Rossi X-ray Timing Explorer's Proportional Counter Array. We find that the auto-encoder closely groups observations together that are classified as similar under the Belloni et al (2000) system, but that there is reasonable grounds for defining each class as made up of components from 3 groups of distinct behaviour. △ Less

Submitted 16 May, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 20 pages, 27 figures. For associated projection code to view the interactive 3D plot, see https://github.com/bjricketts/grs1915-auto-encoder

The Future of Hackathon Research and Practice

Authors: Jeanette Falk, Alexander Nolte, Daniela Huppenkothen, Marion Weinzierl, Kiev Gama, Daniel Spikol, Erik Tollerud, Neil Chue Hong, Ines Knäpper, Linda Bailey Hayden

Abstract: Hackathons are time-bounded collaborative events which have become a global phenomenon adopted by both researchers and practitioners in a plethora of contexts. Hackathon events are generally used to accelerate the development of, for example, scientific results and collaborations, communities, and innovative prototypes addressing urgent challenges. As hackathons have been adopted into many differe… ▽ More Hackathons are time-bounded collaborative events which have become a global phenomenon adopted by both researchers and practitioners in a plethora of contexts. Hackathon events are generally used to accelerate the development of, for example, scientific results and collaborations, communities, and innovative prototypes addressing urgent challenges. As hackathons have been adopted into many different contexts, the events have also been adapted in numerous ways corresponding to the unique needs and situations of organizers, participants and other stakeholders. While these interdisciplinary adaptions, in general affords many advantages - such as tailoring the format to specific needs - they also entail certain challenges, specifically: 1) limited exchange of best practices, 2) limited exchange of research findings, and 3) larger overarching questions that require interdisciplinary collaboration are not discovered and remain unaddressed. We call for interdisciplinary collaborations to address these challenges. As a first initiative towards this, we performed an interdisciplinary collaborative analysis in the context of a workshop at the Lorentz Center, Leiden in December 2021. In this paper, we present the results of this analysis in terms of six important areas which we envision to contribute to maturing hackathon research and practice: 1) hackathons for different purposes, 2) socio-technical event design, 3) scaling up, 4) making hackathons equitable, 5) studying hackathons, and 6) hackathon goals and how to reach them. We present these areas in terms of the state of the art and research proposals and conclude the paper by suggesting next steps needed for advancing hackathon research and practice. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 20 pages, 3 figures, 1 table

Fourier Domain

Authors: Matteo Bachetti, Daniela Huppenkothen

Abstract: The changes in brightness of an astronomical source as a function of time are key probes into that source's physics. Periodic and quasi-periodic signals are indicators of fundamental time (and length) scales in the system, while stochastic processes help uncover the nature of turbulent accretion processes. A key method of studying time variability is through Fourier methods, the decomposition of t… ▽ More The changes in brightness of an astronomical source as a function of time are key probes into that source's physics. Periodic and quasi-periodic signals are indicators of fundamental time (and length) scales in the system, while stochastic processes help uncover the nature of turbulent accretion processes. A key method of studying time variability is through Fourier methods, the decomposition of the signal into sine waves, which yields a representation of the data in frequency space. With the extension into \textit{spectral timing} the methods built on the Fourier transform can not only help us characterize (quasi-)periodicities and stochastic processes, but also uncover the complex relationships between time, photon energy and flux in order to help build better models of accretion processes and other high-energy dynamical physics. In this Chapter, we provide a broad, but practical overview of the most important relevant methods. △ Less

Submitted 6 August, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 51 pages, 13 figures. Published In: Bambi, C., Santangelo, A. (eds) Handbook of X-ray and Gamma-ray Astrophysics. Springer, Singapore. This revised version adds a Conclusions section, corrects a few typos (including equations) and rewrites a few sentences for clarity. Acknowledgements include a few more people who contributed comments and discussions on the previous version

Journal ref: In: Bambi, C., Santangelo, A. (eds) Handbook of X-ray and Gamma-ray Astrophysics. Springer, Singapore (2023)

Searching for quasi-periodic oscillations in astrophysical transients using Gaussian processes

Authors: M. Hübner, D. Huppenkothen, P. D. Lasky, A. R. Inglis, C. Ick, D. W. Hogg

Abstract: Analyses of quasi-periodic oscillations (QPOs) are important to understanding the dynamic behaviour in many astrophysical objects during transient events like gamma-ray bursts, solar flares, magnetar flares and fast radio bursts. Astrophysicists often search for QPOs with frequency-domain methods such as (Lomb-Scargle) periodograms, which generally assume power-law models plus some excess around t… ▽ More Analyses of quasi-periodic oscillations (QPOs) are important to understanding the dynamic behaviour in many astrophysical objects during transient events like gamma-ray bursts, solar flares, magnetar flares and fast radio bursts. Astrophysicists often search for QPOs with frequency-domain methods such as (Lomb-Scargle) periodograms, which generally assume power-law models plus some excess around the QPO frequency. Time-series data can alternatively be investigated directly in the time domain using Gaussian Process (GP) regression. While GP regression is computationally expensive in the general case, the properties of astrophysical data and models allow fast likelihood strategies. Heteroscedasticity and non-stationarity in data have been shown to cause bias in periodogram-based analyses. Gaussian processes can take account of these properties. Using GPs, we model QPOs as a stochastic process on top of a deterministic flare shape. Using Bayesian inference, we demonstrate how to infer GP hyperparameters and assign them physical meaning, such as the QPO frequency. We also perform model selection between QPOs and alternative models such as red noise and show that this can be used to reliably find QPOs. This method is easily applicable to a variety of different astrophysical data sets. We demonstrate the use of this method on a range of short transients: a gamma-ray burst, a magnetar flare, a magnetar giant flare, and simulated solar flare data. △ Less

Submitted 25 May, 2022; originally announced May 2022.

Comments: 24 pages, 19 figures, submitted to The Astrophysical Journal

Phase-resolved spectroscopy of a quasi-periodic oscillation in the black hole X-ray binary GRS 1915+105 with NICER and NuSTAR

Authors: Edward Nathan, Adam Ingram, Jeroen Homan, Daniela Huppenkothen, Phil Uttley, Michiel van der Klis, Sara Motta, Diego Altamirano, Matthew Middleton

Abstract: Quasi-periodic oscillations (QPOs) are often present in the X-ray flux from accreting stellar-mass black holes (BHs). If they are due to relativistic (Lense-Thirring) precession of an inner accretion flow which is misaligned with the disc, the iron emission line caused by irradiation of the disc by the inner flow will rock systematically between red and blue shifted during each QPO cycle. Here we… ▽ More Quasi-periodic oscillations (QPOs) are often present in the X-ray flux from accreting stellar-mass black holes (BHs). If they are due to relativistic (Lense-Thirring) precession of an inner accretion flow which is misaligned with the disc, the iron emission line caused by irradiation of the disc by the inner flow will rock systematically between red and blue shifted during each QPO cycle. Here we conduct phase-resolved spectroscopy of a $\sim2.2$ Hz type-C QPO from the BH X-ray binary GRS 1915+105, observed simultaneously with NICER and NuSTAR. We apply a tomographic model in order to constrain the QPO phase-dependent illumination profile of the disc. We detect the predicted QPO phase-dependent shifts of the iron line centroid energy, with our best fit featuring an asymmetric illumination profile ($>2σ$ confidence). The observed line energy shifts can alternatively be explained by the spiral density waves of the accretion-ejection instability model. However we additionally measure a significant ($>3σ$) modulation in reflection fraction, strongly favouring a geometric QPO origin. We infer that the disc is misaligned with previously observed jet ejections, which is consistent with the model of a truncated disc with an inner precessing hot flow. However our inferred disc inner radius is small ($r_\text{in}{\sim} 1.4 GM/c^2$). For this disc inner radius, Lense-Thirring precession cannot reproduce the observed QPO frequency. In fact, this disc inner radius is incompatible with the predictions of all well-studied QPO models in the literature. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 26 pages, 16 figures, accepted to MNRAS

Orbital decay in M82 X-2

Authors: Matteo Bachetti, Marianne Heida, Thomas Maccarone, Daniela Huppenkothen, Gian Luca Israel, Didier Barret, Murray Brightman, McKinley Brumback, Hannah P. Earnshaw, Karl Forster, Felix Fürst, Brian W. Grefenstette, Fiona A. Harrison, Amruta D. Jaodand, Kristin K. Madsen, Matthew Middleton, Sean N. Pike, Maura Pilia, Juri Poutanen, Daniel Stern, John A. Tomsick, Dominic J. Walton, Natalie Webb, Jörn Wilms

Abstract: M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenar… ▽ More M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenarios. In this Paper, we follow the orbit of the neutron star for seven years, measure the decay of the orbit ($\dot{P}_{orb}/{P}_{orb}\approx-8\cdot10^{-6}\mathrm{yr}^{-1}$), and argue that this orbital decay is driven by extreme mass transfer of more than 150 times the mass transfer limit set by the Eddington luminosity. If this is true, the mass available to the accretor is more than enough to justify its luminosity, with no need for beaming. This also strongly favors models where the accretor is a highly-magnetized neutron star. △ Less

Submitted 29 August, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 18 pages, 6 figures, accepted for publication in ApJ

Journal ref: ApJ 937, 125 (2022)

Characterizing Sparse Asteroid Light Curves with Gaussian Processes

Authors: Christina Willecke Lindberg, Daniela Huppenkothen, R. Lynne Jones, Bryce T. Bolin, Mario Juric, V. Zach Golkhou, Eric C. Bellm, Andrew J. Drake, Matthew J. Graham, Russ R. Laher, Ashish A. Mahabal, Frank J. Masci, Reed Riddle, Kyung Min Shin

Abstract: In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to ov… ▽ More In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to overcome these sampling limitations, we introduce a flexible model based on Gaussian Processes to enable Bayesian parameter inference of asteroid time series data. This model is designed to be flexible and extensible, and can model multiple asteroid properties such as the rotation period, light curve amplitude, changing pulse profile, and magnitude changes due to the phase angle evolution at the same time. Here, we focus on the inference of rotation periods. Based on both simulated light curves and real observations from the Zwicky Transient Facility, we show that the new model reliably infers rotational periods from sparsely sampled light curves, and generally provides well-constrained posterior probability densities for the model parameters. We propose this framework as an intermediate method between fast, but very limited period detection algorithms and much more comprehensive, but computationally expensive shape modeling based on ray-tracing codes. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 27 pages, 18 figures, accepted for publication in AJ, associated software available at https://github.com/dirac-institute/asterogap/tree/v0.1

Identification of an X-ray Pulsar in the BeXRB system IGR J18219$-$1347

Authors: B. O'Connor, E. Gogus, D. Huppenkothen, C. Kouveliotou, N. Gorgone, L. J. Townsend, A. Calamida, A. Fruchter, D. A. H. Buckley, M. G. Baring, J. A. Kennea, G. Younes, Z. Arzoumanian, E. Bellm, S. B. Cenko, K. Gendreau, J. Granot, C. Hailey, F. Harrison, D. Hartmann, L. Kaper, A. Kutyrev, P. O. Slane, D. Stern, E. Troja , et al. (3 additional authors not shown)

Abstract: We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore… ▽ More We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore, by comparing the infrared counterpart's spectral energy distribution to known BeXRBs, we confirm this classification and set a distance of approximately $10-15$ kpc for the source. The source's broadband X-ray spectrum ($1.5-50$ keV) is described by an absorbed power-law with photon index $Γ$\,$\sim$\,$0.5$ and cutoff energy at $\sim$\,$13$ keV. △ Less

Submitted 1 February, 2022; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Revised manuscript. Accepted to ApJ

Light curve fingerprints: an automated approach to the extraction of X-ray variability patterns with feature aggregation -- an example application to GRS 1915+105

Authors: Jakub K. Orwat-Kapola, Antony J. Bird, Adam B. Hill, Diego Altamirano, Daniela Huppenkothen

Abstract: Time series data mining is an important field of research in the era of "Big Data". Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light curve characterisation that employs a pipeline consisting of a neural network with a Long-Short Term Memory Variational Autoencoder architecture and… ▽ More Time series data mining is an important field of research in the era of "Big Data". Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light curve characterisation that employs a pipeline consisting of a neural network with a Long-Short Term Memory Variational Autoencoder architecture and a Gaussian mixture model. The pipeline performs extraction and aggregation of features from light curve segments into feature vectors of fixed length which we refer to as light curve "fingerprints". This representation can be readily used as input of down-stream machine learning algorithms. We demonstrate the proposed method on a data set of Rossi X-ray Timing Explorer observations of the galactic black hole X-ray binary GRS 1915+105, which was chosen because of its observed complex X-ray variability. We find that the proposed method can generate a representation that characterises the observations and reflects the presence of distinct classes of GRS 1915+105 X-ray flux variability. We find that this representation can be used to perform efficient classification of light curves. We also present how the representation can be used to quantify the similarity of different light curves, highlighting the problem of the popular classification system of GRS 1915+105 observations, which does not account for intermediate class behaviour. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Pitfalls of periodograms: The non-stationarity bias in the analysis of quasi-periodic oscillations

Authors: Moritz Huebner, Daniela Huppenkothen, Paul D. Lasky, Andrew R. Inglis

Abstract: Quasi-periodic oscillations (QPOs) are an important key to understand the dynamic behavior of astrophysical objects during transient events like gamma-ray bursts, solar flares, and magnetar flares. Searches for QPOs often use the periodogram of the time series and perform spectral density estimation using a Whittle likelihood function. However, the Whittle likelihood is only valid if the time seri… ▽ More Quasi-periodic oscillations (QPOs) are an important key to understand the dynamic behavior of astrophysical objects during transient events like gamma-ray bursts, solar flares, and magnetar flares. Searches for QPOs often use the periodogram of the time series and perform spectral density estimation using a Whittle likelihood function. However, the Whittle likelihood is only valid if the time series is stationary since the frequency bins are otherwise not statistically independent. We show that if time series are non-stationary, the significance of QPOs can be highly overestimated and estimates of the central frequencies and QPO widths can be overconstrained. The effect occurs if the QPO is only present for a fraction of the time series and the noise level is varying throughout the time series. This can occur for example if background noise from before or after the transient is included in the time series or if the low-frequency noise profile varies strongly over the time series. We confirm the presence of this bias in previously reported results from solar flare data and show that significance can be highly overstated. Finally, we provide some suggestions that help identify if an analysis is affected by this bias. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: 19 pages, 11 Figures

Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference

Authors: D. Huppenkothen, M. Bachetti

Abstract: Because many of our X-ray telescopes are optimized towards observing faint sources, observations of bright sources like X-ray binaries in outburst are often affected by instrumental biases. These effects include dead time and photon pile-up, which can dramatically change the statistical inference of physical parameters from these observations. While dead time is difficult to take into account in a… ▽ More Because many of our X-ray telescopes are optimized towards observing faint sources, observations of bright sources like X-ray binaries in outburst are often affected by instrumental biases. These effects include dead time and photon pile-up, which can dramatically change the statistical inference of physical parameters from these observations. While dead time is difficult to take into account in a statistically consistent manner, simulating dead time-affected data is often straightforward. This structure makes the issue of inferring physical properties from dead time-affected observations fall into a class of problems common across many scientific disciplines. There is a growing number of methods to address them under the name of Simulation-Based Inference (SBI), aided by new developments in density estimation and statistical machine learning. In this paper, we introduce SBI as a principled way to infer variability properties from dead time-affected light curves. We use Sequential Neural Posterior Estimation to estimate the posterior probability for variability properties. We show that this method can recover variability parameters on simulated data even when dead time is variable, and present results of an application of this approach to NuSTAR observations of the galactic black hole X-ray binary GRS 1915+105. △ Less

Submitted 23 November, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 21 pages, 24 figures; accepted for publication in MNRAS; code available at https://github.com/dhuppenkothen/DeadTimeSBI and data products at http://doi.org/10.5281/zenodo.4670161

Photometric Classifications of Evolved Massive Stars: Preparing for the Era of Webb and Roman with Machine Learning

Authors: Trevor Z. Dorn-Wallenstein, James R. A. Davenport, Daniela Huppenkothen, Emily M. Levesque

Abstract: In the coming years, next-generation space-based infrared observatories will significantly increase our samples of rare massive stars, representing a tremendous opportunity to leverage modern statistical tools and methods to test massive stellar evolution in entirely new environments. Such work is only possible if the observed objects can be reliably classified. Spectroscopic observations are infe… ▽ More In the coming years, next-generation space-based infrared observatories will significantly increase our samples of rare massive stars, representing a tremendous opportunity to leverage modern statistical tools and methods to test massive stellar evolution in entirely new environments. Such work is only possible if the observed objects can be reliably classified. Spectroscopic observations are infeasible with more distant targets, and so we wish to determine whether machine learning methods can classify massive stars using broadband infrared photometry. We find that a Support Vector Machine classifier is capable of coarsely classifying massive stars with labels corresponding to hot, cool, and emission line stars with high accuracy, while rejecting contaminating low mass giants. Remarkably, 76\% of emission line stars can be recovered without the need for narrowband or spectroscopic observations. We classify a sample of ${\sim}2500$ objects with no existing labels, and identify fourteen candidate emission line objects. Unfortunately, despite the high precision of the photometry in our sample, the heterogeneous origins of the labels for the stars in our sample severely inhibits our classifier from distinguishing classes of stars with more granularity. Ultimately, no large and homogeneously labeled sample of massive stars currently exists. Without significant efforts to robustly classify evolved massive stars -- which is feasible given existing data from large all-sky spectroscopic surveys -- shortcomings in the labeling of existing data sets will hinder efforts to leverage the next-generation of space observatories. △ Less

Submitted 22 March, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 23 pages, 16 figures, 6 tables. Accepted for publication in ApJ. Comments welcome

Rapid Spectral Variability of a Giant Flare from a Magnetar in NGC 253

Authors: O. J. Roberts, P. Veres, M. G. Baring, M. S. Briggs, C. Kouveliotou, E. Bissaldi, G. Younes, S. I. Chastain, J. J. DeLaunay, D. Huppenkothen, A. Tohuvavohu, P. N. Bhat, E. Gogus, A. J. van der Horst, J. A. Kennea, D. Kocevski, J. D. Linford, S. Guiriec, R. Hamburg, C. A. Wilson-Hodge, E. Burns

Abstract: Magnetars are slowly-rotating neutron stars with extremely strong magnetic fields ($10^{13-15}$ G), episodically emitting $\sim100$ ms long X-ray bursts with energies of $\sim10^{40-41}$ erg. Rarely, they produce extremely bright, energetic giant flares that begin with a short ($\sim0.2$ s), intense flash, followed by fainter, longer lasting emission modulated by the magnetar spin period (typicall… ▽ More Magnetars are slowly-rotating neutron stars with extremely strong magnetic fields ($10^{13-15}$ G), episodically emitting $\sim100$ ms long X-ray bursts with energies of $\sim10^{40-41}$ erg. Rarely, they produce extremely bright, energetic giant flares that begin with a short ($\sim0.2$ s), intense flash, followed by fainter, longer lasting emission modulated by the magnetar spin period (typically 2-12 s), thus confirming their origin. Over the last 40 years, only three such flares have been observed in our local group; they all suffered from instrumental saturation due to their extreme intensity. It has been proposed that extra-galactic giant flares likely constitute a subset of short gamma-ray bursts, noting that the sensitivity of current instrumentation prevents us from detecting the pulsating tail, while the initial bright flash is readily observable out to distances $\sim 10-20$ Mpc. Here, we report X- and gamma-ray observations of GRB 200415A, which exhibits a rapid onset, very fast time variability, flat spectra and significant sub-millisecond spectral evolution. These attributes match well with those expected for a giant flare from an extra-galactic magnetar, noting that GRB 200415A is directionally associated with the galaxy NGC 253 ($\sim$3.5 Mpc away). The detection of $\sim3$ MeV photons provides definitive evidence for relativistic motion of the emitting plasma. The observed rapid spectral evolution can naturally be generated by radiation emanating from such rapidly-moving gas in a rotating magnetar. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Extending the $Z^2_n$ and $H$ statistics to generic pulsed profiles

Authors: Matteo Bachetti, Maura Pilia, Daniela Huppenkothen, Scott M. Ransom, Stefano Curatti, Alessandro Ridolfi

Abstract: The search for astronomical pulsed signals within noisy data, in the radio band, is usually performed through an initial Fourier analysis to find "candidate" frequencies and then refined through the folding of the time series using trial frequencies close to the candidate. In order to establish the significance of the pulsed profiles found at these trial frequencies, pulsed profiles are evaluated… ▽ More The search for astronomical pulsed signals within noisy data, in the radio band, is usually performed through an initial Fourier analysis to find "candidate" frequencies and then refined through the folding of the time series using trial frequencies close to the candidate. In order to establish the significance of the pulsed profiles found at these trial frequencies, pulsed profiles are evaluated with a chi-squared test, to establish how much they depart from a null hypothesis where the signal is consistent with a flat distribution of noisy measurements. In high-energy astronomy, the chi-squared statistic has widely been replaced by the $Z^2_n$ statistic and the H-test as they are more sensitive to extra information such as the harmonic content of the pulsed profile. The $Z^2_n$ statistic and H-test were originally developed for the use with "event data", composed of arrival times of single photons, leaving it unclear how these methods could be used in radio astronomy. In this paper, we present a version of the $Z^2_n$ statistic and H-test for pulse profiles with Gaussian uncertainties, appropriate for radio or even optical pulse profiles. We show how these statistical indicators provide better sensitivity to low-significance pulsar candidates with respect to the usual chi-squared method, and a straightforward way to discriminate between pulse profile shapes. Moreover, they provide an additional tool for Radio Frequency Interference (RFI) rejection. △ Less

Submitted 4 March, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 15 pages, 5 figures

Journal ref: ApJ 909, 33 (2021)

Simultaneous magnetic polar cap heating during a flaring episode from the magnetar 1RXS J170849.0-400910

Authors: George Younes, Matthew G. Baring, Chryssa Kouveliotou, Zorawar Wadiasingh, Daniela Huppenkothen, Alice K. Harding

Abstract: During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of… ▽ More During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: Accepted for publication in the Astrophysical Journal Letters

Discovery and Identification of MAXI J1621-501 as a Type I X-ray Burster with a Super-Orbital Period

Authors: Nicholas M. Gorgone, Chryssa Kouveliotou, Hitoshi Negoro, Ralph A. M. J. Wijers, Enrico Bozzo, Sylvain Guiriec, Peter Bult, Daniela Huppenkothen, Ersin Gogus, Arash Bahramian, Jamie Kennea, Justin D. Linford, James Miller-Jones, Matthew G. Baring, Paz Beniamini, Deepto Chakrabarty, Jonathan Granot, Charles Hailey, Fiona A. Harrison, Dieter H. Hartmann, Wataru Iwakiri, Lex Kaper, Erin Kara, Simona Mazzola, Katsuhiro Murata , et al. (4 additional authors not shown)

Abstract: MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identi… ▽ More MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $Γ$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods. △ Less

Submitted 22 August, 2019; v1 submitted 9 August, 2019; originally announced August 2019.

Comments: 24 pages, Accepted for publication in The Astrophysical Journal (8/21/2019)

Astro2020: Training the Future Generation of Computational Researchers

Authors: Gurtina Besla, Daniela Huppenkothen, Nicole Lloyd-Ronning, Evan Schneider, Peter Behroozi, Blakesley Burkhart, C. K. Chan, Seth A. Jacobson, Sarah Morrison, Hai Ah Nam, Smadar Naoz, Annika Peter, Enrico Ramirez-Ruiz

Abstract: The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics. The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper: State of the Profession Consideration

Entrofy Your Cohort: A Data Science Approach to Candidate Selection

Authors: D. Huppenkothen, B. McFee, L. Norén

Abstract: Selecting a cohort from a set of candidates is a common task within and beyond academia. Admitting students, awarding grants, choosing speakers for a conference are situations where human biases may affect the make-up of the final cohort. We propose a new algorithm, Entrofy, designed to be part of a larger decision making strategy aimed at making cohort selection as just, quantitative, transparent… ▽ More Selecting a cohort from a set of candidates is a common task within and beyond academia. Admitting students, awarding grants, choosing speakers for a conference are situations where human biases may affect the make-up of the final cohort. We propose a new algorithm, Entrofy, designed to be part of a larger decision making strategy aimed at making cohort selection as just, quantitative, transparent, and accountable as possible. We suggest this algorithm be embedded in a two-step selection procedure. First, all application materials are stripped of markers of identity that could induce conscious or sub-conscious bias. During blind review, the committee selects all applicants, submissions, or other entities that meet their merit-based criteria. This often yields a cohort larger than the admissible number. In the second stage, the target cohort can be chosen from this meritorious pool via a new algorithm and software tool. Entrofy optimizes differences across an assignable set of categories selected by the human committee. Criteria could include gender, academic discipline, experience with certain technologies, or other quantifiable characteristics. The Entrofy algorithm yields the computational maximization of diversity by solving the tie-breaking problem with provable performance guarantees. We show how Entrofy selects cohorts according to pre-determined characteristics in simulated sets of applications and demonstrate its use in a case study. This cohort selection process allows human judgment to prevail when assessing merit, but assigns the assessment of diversity to a computational process less likely to be beset by human bias. Importantly, the stage at which diversity assessments occur is fully transparent and auditable with Entrofy. Splitting merit and diversity considerations into their own assessment stages makes it easier to explain why a given candidate was selected or rejected. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 22 pages, 4 figures, submitted to PLOS One. The accompanying software is available at https://github.com/dhuppenkothen/entrofy

Introducing Bayesian Analysis with $\text{m&m's}^\circledR$: an active-learning exercise for undergraduates

Authors: Gwendolyn Eadie, Daniela Huppenkothen, Aaron Springford, Tyler McCormick

Abstract: We present an active-learning strategy for undergraduates that applies Bayesian analysis to candy-covered chocolate $\text{m&m's}^\circledR$. The exercise is best suited for small class sizes and tutorial settings, after students have been introduced to the concepts of Bayesian statistics. The exercise takes advantage of the non-uniform distribution of $\text{m&m's}^\circledR~$ colours, and the di… ▽ More We present an active-learning strategy for undergraduates that applies Bayesian analysis to candy-covered chocolate $\text{m&m's}^\circledR$. The exercise is best suited for small class sizes and tutorial settings, after students have been introduced to the concepts of Bayesian statistics. The exercise takes advantage of the non-uniform distribution of $\text{m&m's}^\circledR~$ colours, and the difference in distributions made at two different factories. In this paper, we provide the intended learning outcomes, lesson plan and step-by-step guide for instruction, and open-source teaching materials. We also suggest an extension to the exercise for the graduate-level, which incorporates hierarchical Bayesian analysis. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: Accepted to the Journal of Statistics Education (in press); 15 pages, 7 figures

Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

Authors: A. Siemiginowska, G. Eadie, I. Czekala, E. Feigelson, E. B. Ford, V. Kashyap, M. Kuhn, T. Loredo, M. Ntampaka, A. Stevens, A. Avelino, K. Borne, T. Budavari, B. Burkhart, J. Cisewski-Kehe, F. Civano, I. Chilingarian, D. A. van Dyk, G. Fabbiano, D. P. Finkbeiner, D. Foreman-Mackey, P. Freeman, A. Fruscione, A. A. Goodman, M. Graham , et al. (27 additional authors not shown)

Abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New met… ▽ More Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years

Authors: Paul S. Ray, Zaven Arzoumanian, David Ballantyne, Enrico Bozzo, Soren Brandt, Laura Brenneman, Deepto Chakrabarty, Marc Christophersen, Alessandra DeRosa, Marco Feroci, Keith Gendreau, Adam Goldstein, Dieter Hartmann, Margarita Hernanz, Peter Jenke, Erin Kara, Tom Maccarone, Michael McDonald, Michael Nowak, Bernard Phlips, Ron Remillard, Abigail Stevens, John Tomsick, Anna Watts, Colleen Wilson-Hodge , et al. (134 additional authors not shown)

Abstract: We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniqu… ▽ More We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science. △ Less

Submitted 8 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

The Zwicky Transient Facility: Science Objectives

Authors: Matthew J. Graham, S. R. Kulkarni, Eric C. Bellm, Scott M. Adams, Cristina Barbarino, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Patrick R. Brady, S. Bradley Cenko, Chan-Kao Chang, Michael W. Coughlin, Kishalay De, Gwendolyn Eadie, Tony L. Farnham, Ulrich Feindt, Anna Franckowiak, Christoffer Fremling, Avishay Gal-yam, Suvi Gezari, Shaon Ghosh, Daniel A. Goldstein, V. Zach Golkhou, Ariel Goobar, Anna Y. Q. Ho , et al. (92 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single… ▽ More The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 26 pages, 7 figures, Published in PASP Focus Issue on the Zwicky Transient Facility

The Zwicky Transient Facility: System Overview, Performance, and First Results

Authors: Eric C. Bellm, Shrinivas R. Kulkarni, Matthew J. Graham, Richard Dekany, Roger M. Smith, Reed Riddle, Frank J. Masci, George Helou, Thomas A. Prince, Scott M. Adams, C. Barbarino, Tom Barlow, James Bauer, Ron Beck, Justin Belicki, Rahul Biswas, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Valery Brinnel, Tim Brooke, Brian Bue, Mattia Bulla, Rick Burruss, S. Bradley Cenko , et al. (91 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation… ▽ More The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: Published in PASP Focus Issue on the Zwicky Transient Facility (https://dx.doi.org/10.1088/1538-3873/aaecbe). 21 Pages, 12 Figures

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 131, Issue 995, pp. 018002 (2019)

Stingray: A Modern Python Library For Spectral Timing

Authors: D. Huppenkothen, M. Bachetti, A. L. Stevens, S. Migliari, P. Balm, O. Hammad, U. M. Khan, H. Mishra, H. Rashid, S. Sharma, R. V. Blanco, E. M. Ribeiro

Abstract: This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build… ▽ More This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. Here, we describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well-tested, with a test coverage of currently 95%, and is accompanied by extensive API documentation and a set of step-by-step tutorials. △ Less

Submitted 9 August, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 16 pages, 6 figures; published in ApJ Vol 881, No 1

Physics and astrophysics of strong magnetic field systems with eXTP

Authors: Andrea Santangelo, Silvia Zane, Hua Feng, RenXin Xu, Victor Doroshenko, Enrico Bozzo, Ilaria Caiazzo, Francesco Coti Zelati, Paolo Esposito, Denis González-Caniulef, Jeremy Heyl, Daniela Huppenkothen, Gianluca Israel, ZhaoSheng Li, Lin Lin, Roberto Mignani, Nanda Rea, Mauro Orlandini, Roberto Taverna, Hao Tong, Roberto Turolla, Cristina Baglio, Federico Bernardini, Niccoló Bucciantini, Marco Feroci , et al. (16 additional authors not shown)

Abstract: In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Conso… ▽ More In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

NuSTAR Hard X-ray View of Low-luminosity Active Galactic Nuclei: High-energy Cutoff and Truncated Thin Disk

Authors: George Younes, Andrew Ptak, Luis C. Ho, Fu-Guo Xie, Yuichi Terasima, Feng Yuan, Daniela Huppenkothen, Mihoko Yukita

Abstract: We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combinatio… ▽ More We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combination of a hot thermal plasma model, a power-law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with a $3σ$ reflection fraction upper-limits $R<0.3$ and $R<0.18$ for NGC 3998 and NGC 4579, respectively; (2) the 6-7 keV line complex in NGC 4579 could either be fit with a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line, or 3 relatively narrow lines, which includes contribution from Fe XXVI; (3) NGC 4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy $E_{\rm cut}=107_{-18}^{+27}$ keV in NGC 3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) NGC 3998 spectrum is consistent with a Comptonization model with either a sphere ($τ\approx3\pm1$) or slab ($τ\approx1.2\pm0.6$) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGN, other LLAGN, and hot accretion flow models. △ Less

Submitted 2 December, 2018; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures, 4 tables, accepted for publication in ApJ

Constraining the limiting brightness temperature and Doppler factors for the largest sample of radio bright blazars

Authors: I. Liodakis, T. Hovatta, D. Huppenkothen, S. Kiehlmann, W. Max-Moerbeck, A. C. S. Readhead

Abstract: Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called {\it Magnetron} to model the light curves of 1029 sources observed by the Owens V… ▽ More Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called {\it Magnetron} to model the light curves of 1029 sources observed by the Owens Valley Radio Observatory's 40-m telescope as a series of flares with an exponential rise and decay, and estimate their variability brightness temperature. Our analysis allows us to place the most stringent constraints on the equipartition brightness temperature i.e., the maximum achieved intrinsic brightness temperature in beamed sources which we found to be $\rm \langle T_{eq}\rangle=2.78\times10^{11}K\pm26\%$. Using our findings we estimated the variability Doppler factor for the largest sample of blazars increasing the number of available estimates in the literature by almost an order of magnitude. Our results clearly show that $γ$-ray loud sources have faster and higher amplitude flares than $γ$-ray quiet sources. As a consequence they show higher variability brightness temperatures and thus are more relativistically beamed, with all of the above suggesting a strong connection between the radio flaring properties of the jet and $γ$-ray emission. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 14 pages, 8 figures, accepted for publication in APJ

The first tidal disruption flare in ZTF: from photometric selection to multi-wavelength characterization

Authors: Sjoert van Velzen, Suvi Gezari, S. Bradley Cenko, Erin Kara, James C. Miller-Jones, Tiara Hung, Joe Bright, Nathaniel Roth, Nadejda Blagorodnova, Daniela Huppenkothen, Lin Yan, Eran Ofek, Jesper Sollerman, Sara Frederick, Charlotte Ward, Matthew J. Graham, Rob Fender, Mansi M. Kasliwal, Chris Canella, Robert Stein, Matteo Giomi, Valery Brinnel, Jakob Santen, Jakob Nordin, Eric C. Bellm , et al. (15 additional authors not shown)

Abstract: We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observati… ▽ More We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST. △ Less

Submitted 7 February, 2019; v1 submitted 7 September, 2018; originally announced September 2018.

Comments: accepted to ApJ, updated with recent Swift data

Detection of non-thermal X-ray emission in the lobes and jets of Cygnus A

Authors: M. N. de Vries, M. W. Wise, D. Huppenkothen, P. E. J. Nulsen, B. Snios, M. J. Hardcastle, M. Birkinshaw, D. M. Worrall, R. T. Duffy, B. R. McNamara

Abstract: We present a spectral analysis of the lobes and X-ray jets of Cygnus A, using more than 2 Ms of $\textit{Chandra}$ observations. The X-ray jets are misaligned with the radio jets and significantly wider. We detect non-thermal emission components in both lobes and jets. For the eastern lobe and jet, we find 1 keV flux densities of $71_{-10}^{+10}$ nJy and $24_{-4}^{+4}$ nJy, and photon indices of… ▽ More We present a spectral analysis of the lobes and X-ray jets of Cygnus A, using more than 2 Ms of $\textit{Chandra}$ observations. The X-ray jets are misaligned with the radio jets and significantly wider. We detect non-thermal emission components in both lobes and jets. For the eastern lobe and jet, we find 1 keV flux densities of $71_{-10}^{+10}$ nJy and $24_{-4}^{+4}$ nJy, and photon indices of $1.72_{-0.03}^{+0.03}$ and $1.64_{-0.04}^{+0.04}$ respectively. For the western lobe and jet, we find flux densities of $50_{-13}^{+12}$ nJy and $13_{-5}^{+5}$ nJy, and photon indices of $1.97_{-0.10}^{+0.23}$ and $1.86_{-0.12}^{+0.18}$ respectively. Using these results, we modeled the electron energy distributions of the lobes as broken power laws with age breaks. We find that a significant population of non-radiating particles is required to account for the total pressure of the eastern lobe. In the western lobe, no such population is required and the low energy cutoff to the electron distribution there needs to be raised to obtain pressures consistent with observations. This discrepancy is a consequence of the differing X-ray photon indices, which may indicate that the turnover in the inverse-Compton spectrum of the western lobe is at lower energies than in the eastern lobe. We modeled the emission from both jets as inverse-Compton emission. There is a narrow region of parameter space for which the X-ray jet can be a relic of an earlier active phase, although lack of knowledge about the jet's electron distribution and particle content makes the modelling uncertain. △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 21 pages, 8 figures

Rotational phase dependence of magnetar bursts

Authors: C. Elenbaas, A. L. Watts, D. Huppenkothen

Abstract: The trigger for the short bursts observed in $γ$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have… ▽ More The trigger for the short bursts observed in $γ$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have been mixed. At the same time, it is not obvious a priori that bursts from a localized active region would actually give rise to a detectable phase-dependence, taking into account issues such as geometry, relativistic effects, and intrinsic burst properties such brightness and duration. In this paper, we build a simple theoretical model to investigate the circumstances under which the latter effects could affect detectability of a dependence of burst emission on rotational phase. We find that even for strongly phase-dependent emission, inferred burst properties may not show a rotational phase dependence depending on the geometry of the system and the observer. Furthermore, the observed properties of bursts with durations short as 10-20% of the spin period can vary strongly depending on the rotational phase at which the burst was emitted. We also show that detectability of a rotational phase dependence depends strongly on the minimum number of bursts observed, and find that existing burst samples may simply be too small to rule out a phase dependence. △ Less

Submitted 6 February, 2018; originally announced February 2018.

Comments: 15 pages, 17 figures

APO Time Resolved Color Photometry of Highly-Elongated Interstellar Object 1I/'Oumuamua

Authors: Bryce T. Bolin, Harold A. Weaver, Yanga R. Fernandez, Carey M. Lisse, Daniela Huppenkothen, R. Lynne Jones, Mario Juric, Joachim Moeyens, Charles A. Schambeau, Colin T. Slater, Zeljko Ivezic, Andrew J. Connolly

Abstract: We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparabl… ▽ More We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017). △ Less

Submitted 29 January, 2018; v1 submitted 13 November, 2017; originally announced November 2017.

Comments: 13 pages, 6 figures, 1 table, accepted for publication in the Astrophysical Journal Letters, data are available at: https://github.com/dhuppenkothen/CometPeriodSearch and https://zenodo.org/record/1068392

Hack Weeks as a model for Data Science Education and Collaboration

Authors: Daniela Huppenkothen, Anthony Arendt, David W. Hogg, Karthik Ram, Jake VanderPlas, Ariel Rokem

Abstract: Across almost all scientific disciplines, the instruments that record our experimental data and the methods required for storage and data analysis are rapidly increasing in complexity. This gives rise to the need for scientific communities to adapt on shorter time scales than traditional university curricula allow for, and therefore requires new modes of knowledge transfer. The universal applicabi… ▽ More Across almost all scientific disciplines, the instruments that record our experimental data and the methods required for storage and data analysis are rapidly increasing in complexity. This gives rise to the need for scientific communities to adapt on shorter time scales than traditional university curricula allow for, and therefore requires new modes of knowledge transfer. The universal applicability of data science tools to a broad range of problems has generated new opportunities to foster exchange of ideas and computational workflows across disciplines. In recent years, hack weeks have emerged as an effective tool for fostering these exchanges by providing training in modern data analysis workflows. While there are variations in hack week implementation, all events consist of a common core of three components: tutorials in state-of-the-art methodology, peer-learning and project work in a collaborative environment. In this paper, we present the concept of a hack week in the larger context of scientific meetings and point out similarities and differences to traditional conferences. We motivate the need for such an event and present in detail its strengths and challenges. We find that hack weeks are successful at cultivating collaboration and the exchange of knowledge. Participants self-report that these events help them both in their day-to-day research as well as their careers. Based on our results, we conclude that hack weeks present an effective, easy-to-implement, fairly low-cost tool to positively impact data analysis literacy in academic disciplines, foster collaboration and cultivate best practices. △ Less

Submitted 31 October, 2017; originally announced November 2017.

Comments: 15 pages, 2 figures, submitted to PNAS, all relevant code available at https://github.com/uwescience/HackWeek-Writeup

No Time for Dead Time: Use the Fourier Amplitude Differences to Normalize Dead-time-affected Periodograms

Authors: Matteo Bachetti, Daniela Huppenkothen

Abstract: Dead time affects many of the instruments used in X-ray astronomy, by producing a strong distortion in power density spectra. This can make it difficult to model the aperiodic variability of the source or look for quasi-periodic oscillations. Whereas in some instruments a simple a priori correction for dead-time-affected power spectra is possible, this is not the case for others such as NuSTAR, wh… ▽ More Dead time affects many of the instruments used in X-ray astronomy, by producing a strong distortion in power density spectra. This can make it difficult to model the aperiodic variability of the source or look for quasi-periodic oscillations. Whereas in some instruments a simple a priori correction for dead-time-affected power spectra is possible, this is not the case for others such as NuSTAR, where the dead time is non-constant and long (~2.5 ms). Bachetti et al. 2015 suggested the cospectrum obtained from light curves of independent detectors within the same instrument as a possible way out, but this solution has always only been a partial one: the measured rms was still affected by dead time, because the width of the power distribution of the cospectrum was modulated by dead time in a frequency-dependent way. In this Letter we suggest a new, powerful method to normalize cospectra and, with some caveats, even power density spectra. Our approach uses the difference of the Fourier amplitudes from two independent detectors to characterize and filter out the effect of dead time. This method is crucially important for the accurate modelling of periodograms derived from instruments affected by dead time on board current missions like NuSTAR and ASTROSAT, but also future missions such as IXPE △ Less

Submitted 31 January, 2018; v1 submitted 27 September, 2017; originally announced September 2017.

Comments: 8 pages, 5 figures, Published on ApJL on 2018 January 30

Journal ref: Matteo Bachetti and Daniela Huppenkothen 2018 ApJL 853, L21

On the Statistical Properties of Cospectra

Authors: D. Huppenkothen, M. Bachetti

Abstract: In recent years, the cross spectrum has received considerable attention as a means of characterising the variability of astronomical sources as a function of wavelength. While much has been written about the statistics of time and phase lags, the cospectrum has only recently been understood as means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as we… ▽ More In recent years, the cross spectrum has received considerable attention as a means of characterising the variability of astronomical sources as a function of wavelength. While much has been written about the statistics of time and phase lags, the cospectrum has only recently been understood as means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as well as a method of characterizing the coherent variability in two wavelength ranges on different time scales. In this paper, we lay out the statistical foundations of the cospectrum, starting with the simplest case of detecting a periodic signal in the presence of white noise. This case is especially relevant for detecting faint X-ray pulsars in detectors heavily affected by instrumental effects, including NuSTAR, Astrosat and IXPE. We show that the statistical distributions of both single and averaged cospectra differ considerably from those for standard periodograms. While a single cospectrum follows a Laplace distribution exactly, averaged cospectra are approximated by a Gaussian distribution only for more than ~30 averaged segments, dependent on the number of trials. We provide an instructive example of a quasi-periodic oscillation in NuSTAR and show that applying standard periodogram statistics leads to underestimated tail probabilities for period detection. We also demonstrate the application of these distributions to a NuSTAR observation of the X-ray pulsar Hercules X-1. △ Less

Submitted 20 December, 2017; v1 submitted 27 September, 2017; originally announced September 2017.

Comments: 12 pages, 8 figures, submitted to ApJ; reproducible code available at https://github.com/dhuppenkothen/cospectra-paper ; Comments welcome

Magnetar giant flare high-energy emission

Authors: C. Elenbaas, D. Huppenkothen, C. Omand, A. L. Watts, E. Bissaldi, I. Caiazzo, J. Heyl

Abstract: High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investig… ▽ More High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the giant flare. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The $RHESSI$ observations of the December 2004 flare modestly favor this latter picture. We assess the prospects for the Fermi Gamma-Ray Space Telescope to detect and characterize a similar high energy component in a future giant flare. Such a detection should help to resolve some of the outstanding issues. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 20 pages, 14 figures

X-ray and radio observations of the magnetar SGR J1935+2154 during its 2014, 2015, and 2016 outbursts

Authors: George Younes, Chryssa Kouveliotou, Amruta Jaodand, Matthew G. Baring, Alexander J. van der Horst, Alice K. Harding, Jason W. T. Hessels, Neil Gehrels, Ramandeep Gill, Daniela Huppenkothen, Jonathan Granot, Ersin Göğüş, Lin Lin

Abstract: We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $Γ=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Cha… ▽ More We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $Γ=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $μ$Jy for the 4.6 GHz observations and 7 $μ$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 17 pages, 7 figures, 5 tables, submitted to ApJ, comments welcome