Pseudospectral implementation of the Einstein-Maxwell system
Jorge Expósito Patiño, Hannes Robert Rüter, David Hilditch
SciPost Phys. 19, 112 (2025) · published 27 October 2025
- doi: 10.21468/SciPostPhys.19.4.112
- Submissions/Reports
-
Abstract
Electromagnetism plays an important role in a variety of applications in gravity that we wish to investigate. To that end, in this work, we present an implementation of the Maxwell equations within the adaptive-mesh pseudospectral numerical relativity code BAMPS. We perform a thorough analysis of the evolution equations as a first order symmetric hyperbolic system of PDEs. This includes both the construction of the characteristic variables for use in our penalty boundary communication scheme, as well as radiation controlling, constraint preserving outer boundary conditions which, for the first time in a numerical context, are shown to be boundary-stable. After choosing a formulation of the Maxwell constraints that we may solve for initial data, we move on to show a suite of numerical tests. Our simulations, both within the Cowling approximation, and in full non-linear evolution, demonstrate rapid convergence of error with resolution, as well as consistency with known quasinormal decay rates on the Kerr background. Finally we evolve the electrovacuum equations of motion with strong data, a good representation of typical critical collapse runs.
Authors / Affiliations: mappings to Contributors and Organizations
See all Organizations.- 1 2 Jorge Expósito Patiño,
- 1 2 Hannes Robert Rüter,
- 1 2 David Hilditch
- 1 Universidade de Lisboa / University of Lisbon
- 2 Instituto Superior Técnico / Instituto Superior Técnico [IST]
- Fundação para a Ciência e a Tecnologia (through Organization: Fundação para a Ciência e Tecnologia [FCT])
- H2020 Future and Emerging Technologies (FET) (through Organization: European Commission [EC])
- Leibniz-Rechenzentrum / Leibniz Supercomputing Centre [LRZ]
