Gravitational-wave memory: waveforms and phenomenology

The non-linear gravitational-wave memory effect is a prediction of general relativity in which test masses are permanently displaced by gravitational radiation. We implement a method for calculating the expected memory waveform from an oscillatory gravitational-wave time series. We use this method to explore the phenomenology of gravitational-wave memory using a numerical relativity surrogate model. Previous methods of calculating the memory have considered only the dominant oscillatory ($\ell=2$, $m=|2|$) mode in the spherical harmonic decomposition or the post-Newtonian expansion. We explore the contribution of higher-order modes and reveal a richer phenomenology than is apparent with $\ell=|m|=2$ modes alone. We also consider the `memory of the memory' in which the memory is, itself, a source of memory, which leads to a small, $O\left(10^{-4}\right)$, correction to the memory waveform. The method is implemented in the python package {\tt\sc GWMemory}, which is made publicly available.