The five significant right singular vectors (RSVs, circles with vertical error bars) from the singular value decomposition (SVD) analysis of the experimental data. The error bars were estimated using the method described in Section 2-3 of the SI, where the LSVs in Supplementary Fig. 5a were used to obtain the time-dependent coefficients (S RSV) as parameterized variables with standard errors. We conducted a global fit analysis with the sum of four exponentials, each with time constants of 130 fs, 1.2 ps, 26 ps (yellow shades), and 292 ps (grey shades to imply that the last dynamics happens slowly after the latest time delay), convoluted with a Gaussian IRF function with a width of 172 fs.