skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: redback : a Bayesian inference software package for electromagnetic transients
ABSTRACT Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce redback; a Bayesian inference software package for electromagnetic transients. redback provides an object-orientated python interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. redback also provides a simple interface for downloading and processing data from various catalogues such as Swift and FINK. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky coverage or a hypothetical user-constructed survey or a generic transient for target-of-opportunity observations with different telescopes. As a demonstration of its capabilities, we show how redback can be used to jointly fit the spectrum and photometry of a kilonova, enabling a more powerful, holistic probe into the properties of a transient. We also showcase general examples of how redback can be used as a tool to simulate transients for realistic surveys, fit models to real, simulated, or private data, multimessenger inference with gravitational waves, and serve as an end-to-end software toolkit for parameter estimation and interpreting the nature of electromagnetic transients.  more » « less
Award ID(s):
2011759
PAR ID:
10548237
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
531
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
1203 to 1227
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Astrophysical transients with rapid developments on subhour timescales are intrinsically rare. Due to their short durations, events like stellar superflares, optical flashes from gamma-ray bursts, and shock breakouts from young supernovae are difficult to identify on timescales that enable spectroscopic follow-up. This paper presents the Evryscope Fast Transient Engine (EFTE), a new data reduction pipeline that is designed to provide low-latency transient alerts from the Evryscopes—a north–south pair of ultra-wide-field telescopes with an instantaneous footprint covering 38% of the entire sky—and tools for building long-term light curves from Evryscope data.EFTEleverages the optical stability of the Evryscopes by using a simple direct image subtraction routine that is suited to continuously monitoring the transient sky at a cadence of a minute. Candidates are produced within the base Evryscope 2 minute cadence for 98.5% of images, and internally filtered usingvetnet, a convolutional neural network real–bogus classifier.EFTEprovides an extensible and robust architecture for transient surveys probing similar timescales, and serves as the software test bed for the real-time analysis pipelines and public data distribution systems for the Argus Array, a next-generation all-sky observatory with a data rate 62 times higher than that of Evryscope. 
    more » « less
  2. Kumar, Sudhir (Ed.)
    Phylogenetic models have become increasingly complex, and phylogenetic data sets have expanded in both size and richness. However, current inference tools lack a model specification language that can concisely describe a complete phylogenetic analysis while remaining independent of implementation details. We introduce a new lightweight and concise model specification language, ‘LPhy’, which is designed to be both human and machine-readable. A graphical user interface accompanies ‘LPhy’, allowing users to build models, simulate data, and create natural language narratives describing the models. These narratives can serve as the foundation for manuscript method sections. Additionally, we present a command-line interface for converting LPhy-specified models into analysis specification files (in XML format) compatible with the BEAST2 software platform. Collectively, these tools aim to enhance the clarity of descriptions and reporting of probabilistic models in phylogenetic studies, ultimately promoting reproducibility of results. 
    more » « less
  3. null (Ed.)
    ABSTRACT octo-tiger is an astrophysics code to simulate the evolution of self-gravitating and rotating systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. octo-tiger is currently optimized to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs (WDs) or main-sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a ‘correction’ algorithm. This code uses hpx parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, WD binary simulations. Results are compared to analytical solutions, when known, and to other grid-based codes such as flash. We also compute the interaction between two WDs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure octo-tiger’s scaling properties up to a core count of ∼80 000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients. 
    more » « less
  4. Hambleton, J. P. (Ed.)
    Discrete element method (DEM) has been widely applied to simulate granular soil behavior. However, traditional DEM uses sphere clusters to approximate realistic particles, which is computationally demanding when simulating many particles. This study explores the use of physics engine, a platform developed for simulating physical processes in video games, to simulate realistic particles. This paper compares realistic particle simulation methodologies using physics engine and discrete element method, including contact models, parameter settings, computational speeds, and simulation results. The results show that the physics engine and DEM achieve similar simulation outputs, while the physics engine runs significantly faster than DEM, because PhysX uses both CPUs (central processing units) and GPUs (graphics processing units) of computers, triangular face tessellations to represent realistic particles, and a simplified contact model to accelerate simulations. This study provides geo-mechanicians and DEM modelers with one more option for them to consider when they simulate realistic particles. 
    more » « less
  5. null (Ed.)
    Abstract Standard procedures for capture–mark–recapture modelling (CMR) for the study of animal demography include running goodness-of-fit tests on a general starting model. A frequent reason for poor model fit is heterogeneity in local survival among individuals captured for the first time and those already captured or seen on previous occasions. This deviation is technically termed a transience effect. In specific cases, simple, uni-state CMR modeling showing transients may allow researchers to assess the role of these transients on population dynamics. Transient individuals nearly always have a lower local survival probability, which may appear for a number of reasons. In most cases, transients arise due to permanent dispersal, higher mortality, or a combination of both. In the case of higher mortality, transients may be symptomatic of a cost of first reproduction. A few studies working at large spatial scales actually show that transients more often correspond to survival costs of first reproduction rather than to permanent dispersal, bolstering the interpretation of transience as a measure of costs of reproduction, since initial detections are often associated with first breeding attempts. Regardless of their cause, the loss of transients from a local population should lower population growth rate. We review almost 1000 papers using CMR modeling and find that almost 40% of studies fitting the searching criteria (N = 115) detected transients. Nevertheless, few researchers have considered the ecological or evolutionary meaning of the transient phenomenon. Only three studies from the reviewed papers considered transients to be a cost of first reproduction. We also analyze a long-term individual monitoring dataset (1988–2012) on a long-lived bird to quantify transients, and we use a life table response experiment (LTRE) to measure the consequences of transients at a population level. As expected, population growth rate decreased when the environment became harsher while the proportion of transients increased. LTRE analysis showed that population growth can be substantially affected by changes in traits that are variable under environmental stochasticity and deterministic perturbations, such as recruitment, fecundity of experienced individuals, and transient probabilities. This occurred even though sensitivities and elasticities of these parameters were much lower than those for adult survival. The proportion of transients also increased with the strength of density-dependence. These results have implications for ecological and evolutionary studies and may stimulate other researchers to explore the ecological processes behind the occurrence of transients in capture–recapture studies. In population models, the inclusion of a specific state for transients may help to make more reliable predictions for endangered and harvested species. 
    more » « less