“…For the nonlinear simulations, our approach is as follows: First, the electric fields at the fundamental wavelength are calculated by the Fourier modal method under the assumption that any self-phase modulation can be neglected. Then, we use these fields and Miller’s rule to calculate the third-order nonlinear polarization within the undepleted pump approximation. , The nonlinear polarization acts as a source at the third-harmonic wavelength. ,, However, instead of calculating the direct emission problem, we make use of the reciprocity principle , and calculate the near fields for p- and s-polarized plane waves incident from the far fields at the third-harmonic wavelength. Then, the electric field radiated to the far field at the third harmonic can be calculated as the overlap integral of these near fields and the nonlinear source. , …”