Good question!
The example article does not mention anything regarding to this. But it should be from the fact that, such waveguide has limited 3rd dimension. Therefore, its effective or equivalent refractive index becomes smaller, with weak mode confinement in 2D. You can refer to varFDTD physics to get some insight on how to calculate the equivalent refractive index from 3D to 2D: MODE - 2.5D varFDTD solver introduction – Ansys Optics
For SiO2, since it is the environment, it is approximate. Using the above method one may get smaller data or sometimes unphysical value, such as smaller than 1.
Anyway, such 2D example is simply an illustration for such initial design with topology optimization. The final design should be verified with 3D FDTD, and you may use 3D FDTD for further optimization. However it will be much slower than 2D.