Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user$varphi$enics: Vainshtein screening with the finite element method
149
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Within the landscape of modified theories of gravity, progress in understanding the behaviour of, and developing tests for, screening mechanisms has been hindered by the complexity of the field equations involved, which are nonlinear in nature and characterised by a large hierarchy of scales. This is especially true of Vainshtein screening, where the fifth force is suppressed by high-order derivative terms which dominate within a radius much larger than the size of the source, known as the Vainshtein radius. In this work, we present the numerical code $varphi$enics, building on the FEniCS library, to solve the full equations of motion from two theories of interest for screening: a model containing high-order derivative operators in the equation of motion and one characterised by nonlinear self-interactions in two coupled scalar fields. We also include functionalities that allow the computation of higher-order operators of the scalar fields in post-processing, enabling us to check that the profiles we find are consistent solutions within the effective field theory. These two examples illustrate the different challenges experienced when trying to simulate such theories numerically, and we show how these are addressed within this code. The examples in this paper assume spherical symmetry, but the techniques may be straightforwardly generalised to asymmetric configurations. This article therefore also provides a worked example of how the finite element method can be employed to solve the screened equations of motion. $varphi$enics is publicly available and can be adapted to solve other theories of screening.
rate research
Read More
We study the Vainshtein mechanism in the context of slowly rotating stars in scalar-tensor theories. While the Vainshtein screening is well established for spherically symmetric spacetimes, we examine its validity in the axisymmetric case for slowly rotating sources. We show that the deviations from the general relativity solution are small in the weak-field approximation outside the star: the solution for the frame-dragging function is the same as in general relativity at leading order. Moreover, in most cases the corrections are suppressed by powers of the Vainshtein radius provided that the screening operates in spherical symmetry. Outside the Vainshtein radius, the frame dragging function receives corrections that are not suppressed by the Vainshtein radius, but which are still subleading. This suggests that the Vainshtein mechanism in general can be extended to slowly rotating stars and that it works analogously to the static case inside the Vainshtein radius. We also study relativistic stars and show that for some theories the frame-dragging function in vacuum does not receive corrections at all, meaning that the screening is perfect outside the star.
We investigate the wave effects of gravitational waves (GWs) using numerical simulations with the finite element method (FEM) based on the publicly available code {it deal.ii}. We robustly test our code using a point source monochromatic spherical wave. We examine not only the waveform observed by a local observer but also the global energy conservation of the waves. We find that our numerical results agree very well with the analytical predictions. Based on our code, we study the scattering of GWs by compact objects. Using monochromatic waves as the input source, we find that if the wavelength of GWs is much larger than the Schwarzschild radius of the compact object, the amplitude of the total scattered GWs does not change appreciably due to the strong diffraction effect, for an observer far away from the scatterer. This finding is consistent with the results reported in the literature. However, we also find that, near the scatterer, not only the amplitude of the scattered waves is very large, comparable to that of the incident waves, but also the phase of the GWs changes significantly due to the interference between the scattered and incident waves. As the evolution of the phase of GWs plays a crucial role in the matched filtering technique in extracting GW signals from the noisy background, our findings suggest that wave effects should be taken into account in the data analysis in the future low-frequency GW experiments, if GWs are scattered by nearby compact objects in our local environment.
The Vainshtein screening mechanism relies on nonlinear interaction terms becoming dominant close to a compact source. However, theories displaying this mechanism are generally understood to be low-energy theories: it is unclear that operators emerging from UV completion do not interfere with terms inducing Vainshtein screening. In this work, we find a set of interacting massive Galileon theories that exhibit Vainshtein screening; examining potential UV completions of these theories, we determine that the screening does not survive the extension. We find that neglecting operators when integrating out a heavy field is non-trivial, and either care must be taken to ensure that omitted terms are small for the whole domain, or one is forced to work solely with the UV theory. We also comment on massive deformations of the familiar Wess-Zumino Galileons.
We study the screening mechanism in the most general scalar-tensor theories that leave gravitational waves unaffected and are thus compatible with recent LIGO/Virgo observations. Using the effective field theory of dark energy approach, we consider the general action for perturbations beyond linear order, focussing on the quasi-static limit. When restricting to the subclass of theories that satisfy the gravitational wave constraints, the fully nonlinear effective Lagrangian contains only three independent parameters. One of these, $beta_1$, is uniquely present in degenerate higher-order theories. We compute the two gravitational potentials for a spherically symmetric matter source and we find that for $beta_1 ge 0$ they decrease as the inverse of the distance, as in standard gravity, while the case $beta_1 < 0$ is ruled out. For $beta_1 > 0$, the two potentials differ and their gravitational constants are not the same on the inside and outside of the body. Generically, the bound on anomalous light bending in the Solar System constrains $beta_1 lesssim 10^{-5}$. Standard gravity can be recovered outside the body by tuning the parameters of the model, in which case $beta_1 lesssim 10^{-2}$ from the Hulse-Taylor pulsar.
Gravitational theories differing from General Relativity may explain the accelerated expansion of the Universe without a cosmological constant. However, to pass local gravitational tests, a screening mechanism is needed to suppress, on small scales, the fifth force driving the cosmological acceleration. We consider the simplest of these theories, i.e. a scalar-tensor theory with first-order derivative self-interactions, and study isolated (static and spherically symmetric) non-relativistic and relativistic stars. We produce screened solutions and use them as initial data for non-linear numerical evolutions in spherical symmetry. We find that these solutions are stable under large initial perturbations, as long as they do not cause gravitational collapse. When gravitational collapse is triggered, the characteristic speeds of the scalar evolution equation diverge, even before apparent black-hole or sound horizons form. This casts doubts on whether the dynamical evolution of screened stars may be predicted in these effective field theories.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
