Simulated normalized changes in the variance in the AFD R-point mode amplitude <u_AFD^2[t]> for the same excitation fluence using the full model presented in equations (1)-(3). The isolated contributions arising from the nonlinear coupling of this mode to the infrared-active phonon Q IR (magenta) and strain Q eta (orange) are also shown.