×

Positional dependence of FDTD mode detection in photonic crystal systems. (English) Zbl 1156.82404

Summary: We have developed an algorithm for evaluating the accuracy and reliability of photonic crystal (PhC) simulations, and used it to analyze the influence of excitation and detector placement in finite-difference time-domain algorithm (FDTD) simulations of two canonical PhC systems. In order to perform this computationally expensive analysis, we evaluated the use of filter diagonalization as an alternative to the Fourier Transform for mode detection, and developed a parallelization algorithm to take advantage of the inherent concurrency in simulating periodic systems. A map of locations where mode detection fails was generated, and we show that this is equivalent to a map of the node densities of the system. In addition to the expected high nodal densities at the symmetry areas of each system, we find more difficult to characterize patterns of high nodal density for the higher-order modes. Based on the observed behavior we are able to provide concrete rules to optimize the detection and excitation of modes in FDTD simulations of PhC systems. Although PhCs were studied, the presented strategies and results apply to the much broader class of all computational time-domain problems where Bloch-Floquet boundary conditions are used.

MSC:

82D25 Statistical mechanics of crystals
82C80 Numerical methods of time-dependent statistical mechanics (MSC2010)
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Shih, High-quality-factor photonic crystal heterostructure laser, Applied Physics Letters 89 pp 101104– (2006)
[2] Nakamura, Ultra-fast photonic crystal/quantum dot all-optical switch for future photonic networks, Optics Express 12 (26) pp 6606– (2004)
[3] Johnson, Block-iterative frequency-domain methods for maxwell’s equations in a planewave basis, Optics Express 8 (3) pp 173– (2001)
[4] Moreno, Modelling of discontinuities in photonic crystal waveguides with the multiple multipole method, Physical Review E 663 pp 036618-1– (2002)
[5] Chan, Order-n spectral method for electromagnetic waves, Physical Review B 51 (23) pp 16635– (1995)
[6] Taflove, Computational Electrodynamics: The Finite-Difference Time-Domain Method (2005) · Zbl 0840.65126
[7] Joannopoulos, Photonic Crystals, Molding the Flow of Light (1995) · Zbl 1035.78500
[8] Farjadpour, Improving accuracy by subpixel smoothing in the finite-difference time domain, Optics Letters 31 pp 2872– (2006)
[9] Main J. Use of harmonic inversion techniques in semiclassical quantization and analysis of quantum spectra. arXiv:chao-dyn/9902008v1.
[10] Wall, Extraction, through filter diagonalization, of general quantum eigenvalues or classical normal mode frequencies from a small number of residues of a short-time segment of a signal, Journal of Chemical Physics 102 pp 8011– (1995)
[11] Burr, Balancing accuracy against computation time: 3-d FDTD for nanophotonics device computation, Proceedings of the SPIE 5733 pp 336– (2005)
[12] Johnson SG. harminv, http://ab-initio.mit.edu/wiki/index.php/Harminv.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.