FIG. 1: Schematic sketch of the protocol for the simulation runs. 20 independent initial configurations are obtained via one long simulation run at 6000 K. For each independent config- uration the system is quenched instantaneously and then fully equilibrated at temperatures T i = 3760 K and 5000 K, followed by instantaneous quenches to the temperatures T f = 3250 K, 3000 K, 2750 K, and 2500 K. At each T $_{f}$, time-dependent cor- relation functions are determined for different waiting times t w . Temperature is kept constant at T f by coupling the system to a Nos´ e-Hoover thermostat. The thermostat is switched off after 0.33 ns, followed by the continuation of the simulations in the microcanonical ensemble for 33 ns. the initial temperatures T$_{i}$ = 5000 K (for 16 . 35 ns) and T$_{i}$ = 3760 K (for 32 . 7 ns), the system was quenched in- stantaneously to T$_{f}$ ∈ { 2500 K , 2750 K , 3000 K , 3250 K } . To disturb the dynamics minimally, we used a Nos´e- Hoover thermostat [45, 46] instead of a stochastic heat bath to keep the temperature at T$_{f}$ constant. A velocity Verlet algorithm was used to integrate the Nos´e-Hoover equations of motion (see Appendix A) with a time step of 1.02 fs. After 0 . 33 ns the Nos´e-Hoover thermostat was switched off and the simulation was continued in the NVE ensemble for 33 ns using a time step of 1 . 6 fs. Whereas previous simulations used instead the NVT ensemble for the whole simulation run, we chose to switch to the NVE ensemble to minimize any influence on the dynamics due to the chosen heat bath algorithm. For the comparison with previous simulations and to check for the lack of a temperature drift, we show in Fig. 2 for exemplatory simulation runs the temperature T = 2 E kin 3 N as a function of time where E$_{kin}$ is the time averaged kinetic energy with fluctuations as indicated with error bars. We find that even after switching off the heat bath [53] there is no temperature drift for T = 3250 K and T = 3000 K and for T = 2500 K and T = 2750 K there is only a slight temper- ature drift which is of the same order as the temperature fluctuations and the drift occurs only for t ⪅ 0 . 6 ns. For all times t ⪆ 0 . 6 ns and for all investigated temperatures there is no temperature drift and thus the comparison with previous simulations valid. In all following, we investigate how the structure and dynamics of the system depend on the waiting time t$_{w}$ elapsed after the quench from T$_{i}$ to T$_{f}$ . We varied the waiting time in the range 0 ns ≤ t$_{w}$ ≤ 23 . 98 ns. Figure 3 shows for the temperature quench T$_{i}$ = 5000 K to T$_{f}$ = 2500 K the partial structure factors [11] where r$_{i}$ and r$_{j}$ are the positions of particles i and j of species α, β = O , Si. The partial structure factors for all other ( T$_{i}$, T$_{f}$ ) combinations are very similar. Although Fig. 3 shows S ( q, t$_{w}$ ) for the largest investigated temper- ature quench, there is only a slight t$_{w}$ -dependence for very short waiting times t$_{w}$ ≤ 0 . 33 ns and almost no t$_{w}$ - dependence for t$_{w}$ ≥ 0 . 33 ns. In this section we focus on the time-dependent gener- alized intermediate incoherent scattering function [11] B. Generalized Incoherent Intermediate Scattering Function C$_{q}$ ( t$_{w}$, t$_{w}$ + t ) = 1 N$_{α}$ 〈 N$_{α}$ ∑ j =1 e i q$_{·}$ ( r$_{j}$ ( t w + t ) $_{−}$r$_{j}$ ( t w ) ) 〉 (3) S$_{αβ}$ ( q, t$_{w}$ ) = 1 N 〈 N$_{α}$ ∑ i =1 N$_{β}$ ∑ j =1 e i q$_{·}$ ( r$_{i}$ ( t w ) $_{−}$r$_{j}$ ( t w ) ) 〉 (2) FIG. 2: (Color online) Temperature T = 2 E$_{kin}$ 3 N as a function of time t for T i = 5000 K , T f = 2500 K , 2750 K , 3000 K , 3250 K shown in each case for the 11th independent simulation run. E kin includes a time average and error bars indicate the corre- sponding fluctuations. III. RESULTS A. Partial Structure Factor 3 w t T NVT NVE (Nose Hoover) (Velocity Verlet) 6000 K  5000 K  3760 K  3250 K  3000 K  2750 K  2500 K  T$_{i}$ T$_{f}$ measurement waiting time t NVT (stochastic) NVT (stochastic) 20 initial configurations  10 -2 10 -1 10 0 10 1 t  [ns] 2000 2500 3000 3500 T  [K] T$_{f}$=3250  K T$_{f}$=2500  K