Lin Zschiedrich

Magnetic Metamaterials at Telecommunication and Visible Frequencies

Arrays of gold split rings with a 50-nm minimum feature size and with an LC resonance at 200 THz frequency (1.5 microm wavelength) are fabricated. For normal-incidence conditions, they exhibit a pronounced fundamental magnetic mode, arising from a coupling via the electric component of the incident light. For oblique incidence, a coupling via the magnetic component is demonstrated as well. Moreover, we identify a novel higher-order magnetic resonance at around 370 THz (800 nm wavelength) that evolves out of the Mie resonance for oblique incidence. Comparison with theory delivers good agreement and also shows that the structures allow for a negative magnetic permeability.

Adaptive finite element method for simulation of optical nano structures

We discuss realization, properties and performance of the adaptive finite element approach to the design of nano-photonic components. Central issues are the construction of vectorial finite elements and the embedding of bounded components into the unbounded and possibly heterogeneous exterior. We apply the finite element method to the optimization of the design of a hollow core photonic crystal fiber. Thereby we look at the convergence of the method and discuss automatic and adaptive grid refinement and the performance of higher order elements.

Solving Time-Harmonic Scattering Problems Based on the Pole Condition II: Convergence of the PML Method

In this paper we study the PML method for Helmholtz-type scattering problems with radially symmetric potential. The PML method consists of surrounding the computational domain with a perfectly matched sponge layer. We prove that the approximate solution obtained by the PML method converges exponentially fast to the true solution in the computational domain as the thickness of the sponge layer tends to infinity. This is a generalization of results by Lassas and Somersalo based on boundary integral equation techniques. Here we use techniques based on the pole condition instead. This makes it possible to treat problems without an explicitly known fundamental solution.

Benchmark of FEM, waveguide, and FDTD algorithms for rigorous mask simulation

An extremely fast time-harmonic finite element solver developed for the transmission analysis of photonic crystals was applied to mask simulation problems. The applicability was proven by examining a set of typical problems and by a benchmarking against two established methods (FDTD and a differential method) and an analytical example. The new finite element approach was up to 100 times faster than the competing approaches for moderate target accuracies, and it was the only method which allowed to reach high target accuracies.

Domain decomposition method for Maxwell’s equations: Scattering off periodic structures

Abstract We present a domain decomposition approach for the computation of the electromagnetic field within periodic structures. We use a Schwarz method with transparent boundary conditions at the interfaces of the domains. Transparent boundary conditions are approximated by the perfectly matched layer method (PML). An adaptive strategy to determine optimal PML parameters is developed. Thus we can treat Wood anomalies appearing in periodic structures. We focus on the application to typical EUV lithography line masks. Light propagation within the multilayer stack of the EUV mask is treated analytically. This results in a drastic reduction of the computational costs and allows for the simulation of next generation lithography masks on a standard personal computer.

Numerical solution of nonlocal hydrodynamic Drude model for arbitrary shaped nano-plasmonic structures using Nédélec finite elements

Nonlocal material response distinctively changes the optical properties of nano-plasmonic scatterers and waveguides. It is described by the nonlocal hydrodynamic Drude model, which - in frequency domain - is given by a coupled system of equations for the electric field and an additional polarization current of the electron gas modeled analogous to a hydrodynamic flow. Recent attempt to simulate such nonlocal model using the finite difference time domain method encountered difficulties in dealing with the grad-div operator appearing in the governing equation of the hydrodynamic current. Therefore, in these studies the model has been simplified with the curl-free hydrodynamic current approximation; but this causes spurious resonances. In this paper we present a rigorous weak formulation in the Sobolev spaces H(curl) for the electric field and H(div) for the hydrodynamic current, which directly leads to a consistent discretization based on Nedelecs finite element spaces. Comparisons with the Mie theory results agree well. We also demonstrate the capability of the method to handle any arbitrary shaped scatterer.

JCMmode : An adaptive finite element solver for the computation of leaky modes

We present our simulation tool JCMmode for calculating propagating modes of an optical waveguide. As ansatz functions we use higher order, vectorial elements (Nedelec elements, edge elements). Further we construct transparent boundary conditions to deal with leaky modes even for problems with inhomogeneous exterior domains as for integrated hollow core Arrow waveguides. We have implemented an error estimator which steers the adaptive mesh refinement. This allows the precise computation of singularities near the metals corner of a Plasmon-Polariton waveguide even for irregular shaped metal films on a standard personal computer.

FEM modeling of 3D photonic crystals and photonic crystal waveguides

We present a finite-element simulation tool for calculating light fields in 3D nano-optical devices. This allows to solve challenging problems on a standard personal computer. We present solutions to eigenvalue problems, like Bloch-type eigenvalues in photonic crystals and photonic crystal waveguides, and to scattering problems, like the transmission through finite photonic crystals. The discretization is based on unstructured tetrahedral grids with an adaptive grid refinement controlled and steered by an error-estimator. As ansatz functions we use higher order, vectorial elements (Nedelec, edge elements). For a fast convergence of the solution we make use of advanced multi-grid algorithms adapted for the vectorial Maxwells equations.

Finite Element Simulation of Radiation Losses in Photonic Crystal Fibers

In our work we focus on the accurate computation of light propagation in finite size photonic crystal structures with the finite element method (FEM). We discuss how we utilize numerical concepts like high-order finite elements, transparent boundary conditions and goal-oriented error estimators for adaptive grid refinement in order to compute radiation leakage in photonic crystal fibers and waveguides. Due to the fast convergence of our method we can use it e.g. to optimize the design of photonic crystal structures with respect to geometrical parameters, to minimize radiation losses and to compute attenutation spectra for different geometries. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Advanced finite element method for nano-resonators

Miniaturized optical resonators with spatial dimensions of the order of the wavelength of the trapped light offer prospects for a variety of new applications like quantum processing or construction of meta-materials. Light propagation in these structures is modelled by Maxwells equations. For a deeper numerical analysis one may compute the scattered field when the structure is illuminated or one may compute the resonances of the structure. We therefore address in this paper the electromagnetic scattering problem as well as the computation of resonances in an open system. For the simulation effcient and reliable numerical methods are required which cope with the infinite domain. We use transparent boundary conditions based on the Perfectly Matched Layer Method (PML) combined with a novel adaptive strategy to determine optimal discretization parameters like the thickness of the sponge layer or the mesh width. Further a novel iterative solver for time-harmonic Maxwells equations is presented.