content_image/438387.jpg content_image/438388.jpg content_image/438382.jpg content_image/438381.jpg
<img src='content_image/438387.jpg'>

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. 

<img src='content_image/438388.jpg'>

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] 

<img src='content_image/438382.jpg'>

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. 

## B. Generalized Incoherent Intermediate Scattering Function 

In this section we focus on the time-dependent gener- alized intermediate incoherent scattering function [11] 

<img src='content_image/438381.jpg'>

## A. Partial Structure Factor 

## III. RESULTS 

