An improved self-gravity treatment in general relativistic smoothed particle hydrodynamics
Abstract
General relativistic smoothed particle hydrodynamics (GRSPH) on a fixed black hole (BH) spacetime is suitable for fluid bodies that are much less massive than the BH, but the self-gravity of the fluid is added as an extra term in the equation of motion. Standard GRSPH implementations add this as a Newtonian source term to the momentum equation, which is no longer valid when the fluid moves at a sizeable fraction of the speed of light, close to the BH. We develop an improved prescription by deriving the weak gravitational field of each fluid element in its local rest frame and transforming it to the global frame where the momentum equation is solved. The resulting force uses rest-frame distances and softening lengths, accounts for the time dilation of the receiving particle, and its motion relative to the source particle. The derivation applies to any weakly self-gravitating flow on a fixed BH background, provided the body is small compared with its distance from the BH. Implemented in a tree-based GRSPH code, the prescription reproduces the self-gravity of a boosted polytrope with a median error below 1\% up to a Lorentz factor of 2, and keeps a binary of two polytropes stable over ten orbits. As an application, we study tidal disruption events (TDEs) of solar-type stars by supermassive BHs. The treatment of self-gravity shifts the disruption threshold and changes both the remnant core mass and the debris energy distribution. Newtonian self-gravity leaves substantial surviving cores in encounters that our prescription predicts to be nearly complete disruptions.
Related
- An improvement to the spatial discretisation in the Astrophysical Weighted Particle Hydrodynamics method
- General-relativistic neutrino-radiation hydrodynamics of crust ejection in magnetar giant flares
- How particle physics stops ultra-relativistic bubbles and shapes their gravitational waves
- Guide field effects on particle acceleration in three-dimensional relativistic magnetic reconnection
- Are E>B regions essential for particle injection in relativistic magnetic reconnection?
BibTeX
@misc{banerjee2026improved,
title = {{An improved self-gravity treatment in general relativistic smoothed particle hydrodynamics}},
author = {Pritam Banerjee and Kowsona Chakraborty and Niles Mondal and Tapobrata Sarkar},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.27775},
eprint = {2610.27775},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.27775},
}
RIS
TY - UNPB
TI - An improved self-gravity treatment in general relativistic smoothed particle hydrodynamics
AU - Pritam Banerjee
AU - Kowsona Chakraborty
AU - Niles Mondal
AU - Tapobrata Sarkar
PY - 2026
DA - 2026/10/20
PB - arXiv
DO - 10.48550/arxiv.2610.27775
UR - https://arxiv.org/abs/2610.27775
AB - General relativistic smoothed particle hydrodynamics (GRSPH) on a fixed black hole (BH) spacetime is suitable for fluid bodies that are much less massive than the BH, but the self-gravity of the fluid is added as an extra term in the equation of motion. Standard GRSPH implementations add this as a Newtonian source term to the momentum equation, which is no longer valid when the fluid moves at a sizeable fraction of the speed of light, close to the BH. We develop an improved prescription by deriving the weak gravitational field of each fluid element in its local rest frame and transforming it to the global frame where the momentum equation is solved. The resulting force uses rest-frame distances and softening lengths, accounts for the time dilation of the receiving particle, and its motion relative to the source particle. The derivation applies to any weakly self-gravitating flow on a fixed BH background, provided the body is small compared with its distance from the BH. Implemented in a tree-based GRSPH code, the prescription reproduces the self-gravity of a boosted polytrope with a median error below 1\% up to a Lorentz factor of 2, and keeps a binary of two polytropes stable over ten orbits. As an application, we study tidal disruption events (TDEs) of solar-type stars by supermassive BHs. The treatment of self-gravity shifts the disruption threshold and changes both the remnant core mass and the debris energy distribution. Newtonian self-gravity leaves substantial surviving cores in encounters that our prescription predicts to be nearly complete disruptions.
ER -