Path integral computation of the lowest order modes in arbitrary shaped inhomogeneous waveguides

A general numerical algorithm for the computa-tion of the fundamental modes and related cutoff wavenumbers in arbitrary shaped inhomogeneous waveguides is presented. The method exploits a generalized Donsker-Ka˘ c formula to ex-press the lowest order modes in terms of asymptotic generalized Wiener-Itˆ o integrals, whose computation is carried out by means of Monte Carlo methods. Comparison with known solutions and computational budget indicate that the proposed method is indeed accurate, versatile, as well as computationally efficient