To evaluate the quality of the silica optimized force field with the best fitness score, the amorphous silica structure introduced in the Fogarty et al. was reconstructed. The amorphous silica system included 2000 SiO 2 molecules with an initial density of 2.2 g/cm 3 (Fig. ). The amorphous silica system was energyminimized to eliminate initial bad contacts. The system was then annealed twice between 300K and 4000K. The first annealing loop was performed using NVT ensemble with heating and cooling rate of 37 K/ps. The second annealing loop was performed in NPT ensemble between 300K-4000K using Nose-Hoover thermostat and barostat 1. the heating and cooling rate was 37 K/ps. Finally, the silica system was equilibrated in NPT ensemble using T=300K and P=1.01325 bar for additional 200 ps as the production run. These calculations were performed using our fitted force field and two previous ReaxFF force fields introduced for such systems. The properties of the final configuration of these silica structures are compared in Table . The radial distributions of the final configuration of silica structure equilibrated using our fitted force field for Si-O, O-O and Si-Si are shown in Fig. . These results show good force field quality for silica structure using JAX-ReaxFF.