Cut-off frequency and dominant eigenfunction computation in complex dielectric geometries via donsker-kač formula and Monte Carlo method