Related papers: Fast Explicit Solutions for Neutrino-Electron Scat…
In physics, accurately modeling large-scale phenomena such as core-collapse supernovae, (CCSN), and neutron star mergers are computationally challenging and require solving large sets of partial and ordinary differential equations.…
A recently published method for solving the neutrino evolution equation with constant matter density is further refined and used to lay out an exact algorithm for computing oscillation probabilities, which is moderately faster than previous…
In contrast to the prevailing view in the literature, it is shown that even extremely stiff sets of ordinary differential equations may be solved efficiently by explicit methods if limiting algebraic solutions are used to stabilize the…
We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends…
%context {Recent observations of hard X-rays and very high energy gamma-rays from a number of young shell type supernova remnants indicate the importance of detailed quantitative studies of energy spectra of relativistic electrons formed…
In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard…
Understanding nuclear effects is essential for improving the sensitivity of neutrino oscillation measurements. Validating nuclear models solely through neutrino scattering data is challenging due to limited statistics and the broad energy…
This work proposes a general strategy for solving possibly nonlinear problems arising from implicit time discretizations as a sequence of explicit solutions. The resulting sequence may exhibit instabilities similar to those of the base…
We examine the effect of spatial resolution on initial mass ejection in grid-based hydrodynamic simulations of binary neutron star mergers. The subset of the dynamical ejecta with velocities greater than $\sim 0.6$c can generate an…
It is shown that the two-flavor neutrino oscillation equations admit an exact analytic solution for arbitrarily chosen normalized electron neutrino population, provided the electron plasma density is adjusted in a certain way. The…
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the…
Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the…
Theoretical predictions for elastic neutrino-electron scattering have no hadronic or nuclear uncertainties at leading order making this process an important tool for normalizing neutrino flux. However, the process is subject to large…
Based on the Kompaneets approximation, we develop a robust methodology to calculate spectral redistribution via inelastic neutrino-nucleon scattering in the context of core-collapse supernova simulations. The resulting equations conserve…
Integral-equation-based fast direct solvers for electromagnetic scattering can substantially reduce computational costs, especially in the presence of multiple excitations. We recently proposed a new high-frequency fast direct solver…
Precision tests of the standard model are essential for constraining models of new physics. Neutrino-electron elastic scattering offers a clean probe into many electroweak effects that are complimentary to the more canonical measurements…
Explicit stabilized integrators are an efficient alternative to implicit or semi-implicit methods to avoid the severe timestep restriction faced by standard explicit integrators applied to stiff diffusion problems. In this paper, we provide…
The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations.…
Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for…
The use of implicit time-stepping schemes for the numerical approximation of solutions to stiff nonlinear time-evolution equations brings well-known advantages including, typically, better stability behaviour and corresponding support of…