Spectral computation of highly oscillatory integral equations in laser theory

Abstract We are concerned in this paper with the numerical computation of the spectra of highly oscillatory integrals that arise in laser simulations. Discretised using the modified Fourier basis, the spectral problem for the integral equation is converted into two independent infinite systems of linear equations whose unknowns are the coefficients of the modified Fourier functions, namely the cosine and shifted sine functions, respectively. Each ( m , n ) entry of the resulting coefficient matrices can be represented exactly by expressions involving the error function with an argument that involves the oscillatory parameter ω and the numbers m and n . Moreover, considering the behaviour of the error function for a large argument, the asymptotics for each entry are analysed for large ω or for large m and n and this enables efficient truncation of the infinite systems. Numerical experiments are provided to illustrate the effectiveness of this method.

[1]  Arieh Iserles,et al.  The computation of the spectra of highly oscillatory Fredholm integral operators , 2011 .

[2]  A. Siegman,et al.  Unstable optical resonators. , 1974, Applied Optics.

[3]  Arieh Iserles,et al.  From high oscillation to rapid approximation I: Modified Fourier expansions , 2008 .

[4]  Byeong-Cheol Lee,et al.  Numerical simulation of laser resonators , 2004 .

[5]  K. Atkinson.  The Numerical Solution of Integral Equations of the Second Kind , 1997 .

[6]  Michael V Berry,et al.  Fractal modes of unstable lasers with polygonal and circular mirrors , 2001 .

[7]  Arieh Iserles,et al.  The Fox-Li operator as a test and a spur for Wiener-Hopf theory , 2012 .

[8]  Arieh Iserles,et al.  On the singular values and eigenvalues of the Fox–Li and related operators , 2010 .

[9]  D. Yakubovich,et al.  Completeness and spectral synthesis of nonselfadjoint one-dimensional perturbations of selfadjoint operators , 2012, 1212.5965.

[10]  V. A. Il'in,et al.  Properties of Spectral Expansions Corresponding to Non-Self-Adjoint Differential Operators , 2003 .

[11]  A. G. Fox,et al.  Resonant modes in a maser interferometer , 1961 .

[12]  Oscar F. Bandtlow,et al.  Explicit upper bounds for the spectral distance of two trace class operators , 2015, 1510.05638.

[13]  Paul Horwitz,et al.  Asymptotic theory of unstable resonator modes , 1973 .

[14]  Michael V Berry,et al.  Theory of unstable laser modes : edge waves and fractality , 2001 .

[15]  Daan Huybrechs,et al.  Computing highly oscillatory integrals , 2017 .

[16]  L. Trefethen,et al.  Spectra and pseudospectra : the behavior of nonnormal matrices and operators , 2005 .

[17]  M. Sugiura.  Unitary representations and harmonic analysis : an introduction , 1990 .

[18]  V. A. Svich,et al.  Eigenmodes diffraction losses of marginally unstable semispherical resonator , 2000, Proceedings of LFNM'2000. 2nd International Workshop on Laser and Fiber-Optical Networks Modeling (Cat. No.00EX419).