The simulations were performed in Lumerical using finite-difference time-domain. The material models were provided by the software's material libraries.