Abstract

Traditionally, one can calculate the update coefficients of the finite-difference time-domain algorithm at material interfaces by averaging the material properties of both sides, which leads to numerical inaccuracies of the reflection depending on the grid resolution. We propose a novel method to calculate the update coefficients such that the algorithm exactly fulfills the boundary conditions at a frequency of optimization, which allows a significant increase in grid spacing while limiting the numerical error. Using the proposed method, we reduced the computational expenses for the full-wave simulation of vertical-cavity surface-emitting lasers such that large structures can be treated without the need to exploit rotational symmetry. The method is demonstrated with the help of several examples.

© 2003 Optical Society of America

Combined analytical-finite difference time-domain full wave simulation of mode-locked vertical-extended-cavity semiconductor lasers

Philipp Kreuter and Bernd Witzigmann
J. Opt. Soc. Am. B 25(7) 1118-1126 (2008)

Different regimes of low-frequency fluctuations in vertical-cavity surface-emitting lasers

Marc Sciamanna, Cristina Masoller, Neal B. Abraham, Fabien Rogister, Patrice Mégret, and Michel Blondel
J. Opt. Soc. Am. B 20(1) 37-44 (2003)

Simulation of whispering gallery modes in the Mie regime using the nonstandard finite-difference time domain algorithm

Naoki Okada and James B. Cole
J. Opt. Soc. Am. B 27(4) 631-639 (2010)

Optical and thermal finite-difference time-domain model for passively mode-locked surface-emitting lasers

Mayank Bahl, Nicolae C. Panoiu, and Richard M. Osgood, Jr.
J. Opt. Soc. Am. B 26(8) 1558-1568 (2009)

Auxiliary differential equation: efficient implementation in the finite-difference time-domain method

T. O. Körner and W. Fichtner
Opt. Lett. 22(21) 1586-1588 (1997)