An official website of the United States government
ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.
more »
« less
KilonovaNet : Surrogate models of kilonova spectra with conditional variational autoencoders
ABSTRACT Detailed radiative transfer simulations of kilonova spectra play an essential role in multimessenger astrophysics. Using the simulation results in parameter inference studies requires building a surrogate model from the simulation outputs to use in algorithms requiring sampling. In this work, we present kilonovanet, an implementation of conditional variational autoencoders (cVAEs) for the construction of surrogate models of kilonova spectra. This method can be trained on spectra directly, removing overhead time of pre-processing spectra, and greatly speeds up parameter inference time. We build surrogate models of three state-of-the-art kilonova simulation data sets and present in-depth surrogate error evaluation methods, which can in general be applied to any surrogate construction method. By creating synthetic photometric observations from the spectral surrogate, we perform parameter inference for the observed light-curve data of GW170817 and compare the results with previous analyses. Given the speed with which kilonovanet performs during parameter inference, it will serve as a useful tool in future gravitational wave observing runs to quickly analyse potential kilonova candidates.
more »
« less
- Award ID(s):
- 2122312
- PAR ID:
- 10353486
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 516
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 1137 to 1148
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We present a coherent, re-usable python framework building on the CosmoPower emulator code for high-accuracy calculations of cosmological observables with Einstein–Boltzmann codes. For detailed statistical analyses, such codes require high computing power, making parameter space exploration costly, especially for beyond-$$\Lambda$$CDM analyses. Machine learning-enabled emulators of Einstein–Boltzmann codes are becoming an increasingly popular solution to this problem. To enable generation, sharing, and use of emulators for inference, we define standards for robustly describing, packaging, and distributing them. We present software for easily performing these tasks in an automated and replicable manner and provide examples and guidelines for generating emulators and wrappers for using them in popular cosmological inference codes. We demonstrate our framework with a suite of high-accuracy emulators for the CAMB code’s calculations of CMB $$C_\ell$$, $P(k)$, background evolution, and derived parameter quantities. We show these emulators are accurate enough for analysing both $$\Lambda$$CDM and a set of extension models ($$N_{\rm eff}$$, $$\sum m_\nu$$, $$w_0 w_a$$) with stage-IV observatories, recovering the original high-accuracy spectra to tolerances well within the cosmic variance uncertainties. We show our emulators also recover cosmological parameters in a simulated cosmic-variance limited experiment, finding results well within $$0.1 \sigma$$ of the input cosmology, while requiring $$\lesssim 1/50$$ of the evaluation time.more » « less
-
While Bayesian inference is the gold standard for uncertainty quantification and propagation, its use within physical chemistry encounters formidable computational barriers. These bottlenecks are magnified for modeling data with many independent variables, such as X-ray/neutron scattering patterns and electromagnetic spectra. To address this challenge, we employ local Gaussian process (LGP) surrogate models to accelerate Bayesian optimization over these complex thermophysical properties. The time-complexity of the LGPs scales linearly in the number of independent variables, in stark contrast to the computationally expensive cubic scaling of conventional Gaussian processes. To illustrate the method, we trained a LGP surrogate model on the radial distribution function of liquid neon and observed a 1,760,000-fold speed-up compared to molecular dynamics simulation, beating a conventional GP by three orders-of-magnitude. We conclude that LGPs are robust and efficient surrogate models poised to expand the application of Bayesian inference in molecular simulations to a broad spectrum of experimental data.more » « less
-
Abstract Real-time ranking of optical transient candidates during gravitational-wave and multimessenger follow-up is challenging when only sparse early-time, multi-band photometry is available. We presentKilonovaSCORER, an open-source framework for scoring and ranking in this regime. It quantifies the consistency of each candidate with a physically motivated kilonova model grid in absolute magnitude space using two complementary per-observation metrics,Ptail,KNeandPnear,KNe. These are aggregated into a cumulative ranking score via inverse-variance weighting in logit space, naturally accounting for heterogeneous observational uncertainties across bands and epochs. A sequential Approximate Bayesian Computation diagnostic tracks photometric consistency across epochs, penalizing candidates whose temporal evolution is incompatible with kilonova expectations. We validate the framework on AT 2017gfo and SN 2025ulz, and test it against supernova simulations under a realistic Rubin/LSST Target-of-Opportunity strategy. The framework recovers kilonova candidates with high confidence while ruling out supernova contaminants within five days of the gravitational-wave trigger. In our LSST ToO simulations, median cumulative scores for thermonuclear and core-collapse supernova contaminants fall to zero by 3–4 days post-trigger, whereas kilonova medians remain ≳0.4.KilonovaSCORERsupports real-time workflows for ToO teams and LSST alert brokers, integrates with follow-up coordination platforms such as the Tool for Rapid Object Vetting and Examination, and is publicly available athttps://github.com/phelipedarc/KilonovaSCORER/tree/main.more » « less
-
Patient-specific cardiovascular simulations are an integral part of cardiovascular research and are increasingly finding use in clinical practice. A primary factor hindering large cohort studies and clinical utilization of cardiovascular simulations is the time-consuming process of constructing simulation-suitable computer models from medical image data. Methods have been proposed to automate the model construction process, but are tailored to either vascular or cardiac models. We herein propose a framework, MeshGrow, that is capable of combining automated modeling of cardiac and extended vascular territories. Namely, we demonstrate that this method can be used to reconstruct the main cardiac chambers along with the aorta and its main sub-branches, and returns a simulation-suitable mesh with a defined aortic valve surface and other surfaces needed for boundary condition specification. MeshGrow combines two different machine learning based techniques to address the specific challenges involved with cardiac and vascular model construction. This results in a two-stage approach; (a) meshing the cardiac structures and (b) growing the vasculature out from it. We present results of our method on five CT test datasets and compare predictions against state-of-the-art benchmark methods, as well as ground truth manually constructed models. Results show that MeshGrow achieves higher metric scores than benchmark methods on the test data. Additionally, we showcase the applicability of the method by running three-dimensional computational fluid dynamics simulations with two of the test cases. With this work, we demonstrate a method that can automate the construction of combined cardiac and vascular models for patient-specific cardiovascular hemodynamic simulation.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1093/mnras/stac2342
