Fixing the dynamical evolution of self-interacting vector fields

Abstract

Numerical simulations of the Cauchy problem for self-interacting massive vector fields often face instabilities and apparent pathologies. We explicitly demonstrate that these issues, previously reported in the literature, are actually due to the breakdown of the well-posedness of the initial-value problem. This is akin to shortcomings observed in scalar-tensor theories when derivative self-interactions are included. Building on previous work done for k-essence, we characterize the well-posedness breakdowns, differentiating between Tricomi and Keldysh-like behaviors. We show that these issues can be avoided by ``fixing the equations'', enabling stable numerical evolutions in spherical symmetry. Additionally, we show that for a class of vector self-interactions, no Tricomi-type breakdown takes place. Finally, we investigate initial configurations for the massive vector field which lead to gravitational collapse and the formation of black holes.

13 pages, 6 figures, refs added, matches PRD version

References in corpus (13)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• GW190521: A Binary Black Hole Merger with a Total Mass of $150M_{\odot }$
• Echoes of ECOs: gravitational-wave signatures of exotic compact objects and of quantum corrections at the horizon scale
• Light ring stability in ultra-compact objects
• Fixing extensions to General Relativity in the non-linear regime
• Spherical Boson Stars as Black Hole mimickers
• Supermassive black holes or boson stars? Hair counting with gravitational wave detectors
• Numerical relativity for Horndeski gravity
• Gravitational collapse of charged scalar fields
• Shock Formation in Lovelock Theories
• Proca-Higgs balls and stars in a UV completion for Proca self-interactions
• Wave propagation and shock formation in the most general scalar-tensor theories
• Efficient implementation of finite volume methods in Numerical Relativity