Nonlinear Dynamics in Quadratic Gravity

We present the first numerically stable nonlinear evolution for quadratic gravity. We first demonstrate our results in spherical symmetry. We explore well-posedness of the respective initial-value problem by simulating randomly perturbed flat-space and black-hole initial data. Our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives. We will also demonstrate current effort to simulate in full 3+1 dimensions using Dendro framework, wavelet adaptive multiresolution code for relativistic astrophysics.