BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Philipp Moesta

Philipp Moesta contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-ph.IMphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2026arXiv

thornado+Flash-X: A Hybrid DG-IMEX and Finite-Volume Framework for Neutrino-Radiation Hydrodynamics in Core-Collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement (AMR)-based multiphysics simulation framework Flash-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $O(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino-matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino-electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved with a hybrid DG-FV representation of the fluid variables and operator-split evolution in Flash-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within Flash-X. thornado's neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How loud are neutron star mergers?

We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $κ^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRHydro: A new open source general-relativistic magnetohydrodynamics code for the Einstein Toolkit

We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein Toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the Toolkit builds upon previous releases and implements the evolution of relativistic magnetised fluids in the ideal MHD limit in fully dynamical spacetimes using the same shock-capturing techniques previously applied to hydrodynamical evolution. In order to maintain the divergence-free character of the magnetic field, the code implements both hyperbolic divergence cleaning and constrained transport schemes. We present test results for a number of MHD tests in Minkowski and curved spacetimes. Minkowski tests include aligned and oblique planar shocks, cylindrical explosions, magnetic rotors, Alfvén waves and advected loops, as well as a set of tests designed to study the response of the divergence cleaning scheme to numerically generated monopoles. We study the code's performance in curved spacetimes with spherical accretion onto a black hole on a fixed background spacetime and in fully dynamical spacetimes by evolutions of a magnetised polytropic neutron star and of the collapse of a magnetised stellar core. Our results agree well with exact solutions where these are available and we demonstrate convergence. All code and input files used to generate the results are available on http://einsteintoolkit.org. This makes our work fully reproducible and provides new users with an introduction to applications of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accurate Simulations of Binary Black-Hole Mergers in Force-Free Electrodynamics

We provide additional information on our recent study of the electromagnetic emission produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. As anticipated in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~ 100 times smaller than the total one, which is predominantly quadrupolar. We here discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L ~ Omega^{10/3-8/3}, while the collimated one scales as L ~ Omega^{5/3-6/3}, thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black-hole horizons as probes of black-hole dynamics I: post-merger recoil in head-on collisions

The understanding of strong-field dynamics near black-hole horizons is a long-standing and challenging prob- lem in general relativity. Recent advances in numerical relativity and in the geometric characterization of black- hole horizons open new avenues into the problem. In this first paper in a series of two, we focus on the analysis of the recoil occurring in the merger of binary black holes, extending the analysis initiated in [1] with Robinson- Trautman spacetimes. More specifically, we probe spacetime dynamics through the correlation of quantities defined at the black-hole horizon and at null infinity. The geometry of these hypersurfaces responds to bulk gravitational fields acting as test screens in a scattering perspective of spacetime dynamics. Within a 3 + 1 approach we build an effective-curvature vector from the intrinsic geometry of dynamical-horizon sections and correlate its evolution with the flux of Bondi linear momentum at large distances. We employ this setup to study numerically the head-on collision of nonspinning black holes and demonstrate its validity to track the qualita- tive aspects of recoil dynamics at infinity. We also make contact with the suggestion that the antikick can be described in terms of a "slowness parameter" and how this can be computed from the local properties of the horizon. In a companion paper [2] we will further elaborate on the geometric aspects of this approach and on its relation with other approaches to characterize dynamical properties of black-hole horizons.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black-hole horizons as probes of black-hole dynamics II: geometrical insights

In a companion paper [1], we have presented a cross-correlation approach to near-horizon physics in which bulk dynamics is probed through the correlation of quantities defined at inner and outer spacetime hypersurfaces acting as test screens. More specifically, dynamical horizons provide appropriate inner screens in a 3+1 setting and, in this context, we have shown that an effective-curvature vector measured at the common horizon produced in a head-on collision merger can be correlated with the flux of linear Bondi-momentum at null infinity. In this paper we provide a more sound geometric basis to this picture. First, we show that a rigidity property of dynamical horizons, namely foliation uniqueness, leads to a preferred class of null tetrads and Weyl scalars on these hypersurfaces. Second, we identify a heuristic horizon news-like function, depending only on the geometry of spatial sections of the horizon. Fluxes constructed from this function offer refined geometric quantities to be correlated with Bondi fluxes at infinity, as well as a contact with the discussion of quasi-local 4-momentum on dynamical horizons. Third, we highlight the importance of tracking the internal horizon dual to the apparent horizon in spatial 3-slices when integrating fluxes along the horizon. Finally, we discuss the link between the dissipation of the non-stationary part of the horizon's geometry with the viscous-fluid analogy for black holes, introducing a geometric prescription for a "slowness parameter" in black-hole recoil dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the detectability of dual jets from binary black holes

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.

0 citations0 reviewsNo code linkedOpen access