Computational Photonics E-Book

Salah Obayya

• Herausgeber: John Wiley & Sons
• Kategorie: Wissenschaft und neue Technologien
• Sprache: Englisch

Seitenzahl: 419

Veröffentlichungsjahr: 2011

Beschreibung

This book explores the state-of-the art in computational modelling techniques for photonic devices In this book, the author provides a comprehensive coverage of modern numerical modelling techniques for designing photonic devices for use in modern optical telecommunications systems. In addition the book presents the state-of-the-art in computational photonics techniques, covering methods such as full-vectorial finite-element beam propagation, bidirectional beam propagation, complex-envelope alternative direction implicit finite difference time domain, multiresolution time domain, and finite volume time domain. The book guides the reader through the concepts of modelling, analysing, designing and optimising the performance of a wide range of photonic devices by building their own numerical code using these methods. Key Features: * Provides a thorough presentation of the state-of-the art in computational modelling techniques for photonics * Contains broad coverage of both frequency- and time-domain techniques to suit a wide range of photonic devices * Reviews existing commercial software packages for photonics * Presents the advantages and disadvantages of the different modelling techniques as well as their suitability for various photonic devices * Shows the reader how to model, analyse, design and optimise the performance of a wide range of photonic devices by building their own numerical code using these methods * Accompanying website contains the numerical examples representing the numerical techniques in this book, as well as several design examples (href="http://www.wiley.com/go/obayya_computational">http://www.wiley.com/go/obayya_computational) This book will serve as an invaluable reference for researchers, optical telecommunications engineers, engineers in the photonics industry. PhD and MSc students undertaking courses in the areas of photonics and optical telecommunications will also find this book of interest.

Leseprobe

Contents

Preface

Acknowledgements

1 Introduction

1.1 Photonics: The Countless Possibilities of Light Propagation

1.2 Modelling Photonics

2 Full-Vectorial Beam Propagation Methods

2.1 Introduction

2.2 Overview of the Beam Propagation Methods

2.3 Maxwell’s Equations

2.4 Magnetic-Field Formulation of the Wave Equation

2.5 Electric-Field Formulation of the Wave Equation

2.6 Perfectly Matched Layer

2.7 Finite-Element Analysis

2.8 Derivation of BPM Equations

2.9 Imaginary-Distance BPM: Mode Solver

3 Assessment of the Full-Vectorial Beam Propagation Method

3.1 Introduction

3.2 Analysis of Rectangular Waveguide

3.3 Photonic Crystal Fibre

3.4 Liquid-Crystal-Based Photonic Crystal Fibre

3.5 Electro-Optical Modulators

3.6 Switches 61 References

4 Bidirectional Beam Propagation Method

4.1 Introduction

4.2 Optical Waveguide Discontinuity Problem

4.3 Finite-Element Analysis of Discontinuity Problems

4.4 Derivation of Finite-Element Matrices

4.5 Application of Taylor’s Series Expansion

4.6 Computation of Reflected, Transmitted and Radiation Waves

4.7 Optical Fibre Facet Problem

4.8 Finite-Element Analysis of Optical Fibre Facets

4.9 Iterative Analysis of Multiple Discontinuities

4.10 Numerical Assessment

5 Complex-Envelope Alternating-Direction-Implicit Finite-Difference Time-Domain Method

5.1 Introduction

5.2 Maxwell’s Equations

5.3 Brief History of the Finite-Difference Time-Domain (FDTD) Method

5.4 Finite-Difference Time-Domain (FDTD) Method

5.5 Alternating-Direction-Implicit FDTD (ADI-FDTD): Beyond the Courant Limit

5.6 Complex-Envelope ADI-FDTD (CE-ADI-FDTD)

5.7 Perfectly Matched Layer (PML) Boundary Conditions

5.8 Uniaxial Perfectly Matched Layer (UPML) Absorbing Boundary Condition

5.9 PML Parameters

6 Finite-Volume Time-Domain Method

6.1 Introduction

6.2 Numerical Analysis

6.3 The UPWIND Scheme for Flux Calculation

6.4 Nondiffusive Scheme for the Flux Calculation

6.5 2D Formulation of the FVTD Method

6.6 Boundary Conditions

6.7 Nonlinear Optics

6.8 Nonlinear Optical Interactions

6.9 Extension of the FDTD Method to Nonlinear Problems

7 Numerical Analysis of Linear and Nonlinear PhC-Based Devices

7.1 Introduction

7.2 FVTD Method Assessment: PhC Cavity

7.3 FVTD Method Assessment: PhC Waveguide

7.4 FVTD Method Assessment: PBG T-Branch

7.5 PhC Multimode Resonant Cavity

7.6 FDTD Analysis of Nonlinear Devices

7.7 FVTD Analysis of Nonlinear Photonic Crystal Wires

8 Multiresolution Time Domain

8.1 Introduction

8.2 MRTD Basics

8.3 MRTD Update Scheme

8.4 Scaling MRTD

9 MRTD Analysis of PhC Devices

9.1 Introduction

9.2 UPML-MRTD: Test and Code Validation

9.3 MRTD versus FDTD for the Analysis of Linear Photonic Devices

10 MRTD Analysis of SHG PhC Devices

10.1 Introduction

10.2 Second Harmonic Generation in Optics

10.3 Extended S-MRTD for SHG Analysis

10.4 SHG in PhC-Waveguides

10.5 Selective SHG in Compound PhC-Based Structures

10.6 New Design for Selective SHG: PhC Microcavity Coupling

11 Dispersive Nonlinear MRTD for SHG Applications

11.1 Introduction

11.2 Dispersion Analysis

11.3 SHG-MRTD Scheme for Dispersive Materials

11.4 Simulation Results

Index

This edition first published 2011© 2011 John Wiley & Sons, Ltd

Registered officeJohn Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex, PO19 8SQ, United Kingdom

For details of our global editorial offices, for customer services and for information about how to apply for permission to reuse the copyright material in this book please see our website at www.wiley.com.

The right of the author to be identified as the author of this work has been asserted in accordance with the Copyright, Designs and Patents Act 1988.

All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording or otherwise, except as permitted by the UK Copyright, Designs and Patents Act 1988, without the prior permission of the publisher.

Wiley also publishes its books in a variety of electronic formats. Some content that appears in print may not be available in electronic books.

Designations used by companies to distinguish their products are often claimed as trademarks. All brand names and product names used in this book are trade names, service marks, trademarks or registered trademarks of their respective owners. The publisher is not associated with any product or vendor mentioned in this book. This publication is designed to provide accurate and authoritative information in regard to the subject matter covered. It is sold on the understanding that the publisher is not engaged in rendering professional services. If professional advice or other expert assistance is required, the services of a competent professional should be sought.

Library of Congress Cataloging-in-Publication Data

Obayya, Salah.Computational photonics/Salah Obayya.p. cm.Includes bibliographical references and index.ISBN 9780470688939 (cloth)1. Optoelectronic devices–Mathematical models. 2. Photonics–Mathematics. I. Title.TK8304.O23 2010621.36–dc222010019376

A catalogue record for this book is available from the British Library.

Print ISBN: 978-0-470-68893-9 (H/B)ePDF ISBN: 9780470667071oBook ISBN: 9780470667064

All Praise is due to Allah, and peace and blessings be upon. ProphetMuhammad and upon his family and his Companions.This work is dedicated to my mother and my late father who raised me up tothe good values of hardworking with full devotion.I would like to dedicate this work also to my wife, Mona, whose love, support,patience, understanding are beyond any scope. In a word, to Monawho means everything to me.

Preface

Photonics as a field today covers a huge range from communications to science and technology applications, including laser manufacturing, biological and chemical sensing, display technology and optical computing.

For a number of reasons ranging from a higher bandwidth to the inexpensiveness of optical-fibre materials, components based on electric currents are more and more complemented or replaced by technology based on light. Relatively cheap, compact and highly sensitive photonic devices have already been commercialised for a variety of areas, and are expected to be basic blocks for all-optical integrated circuits. As the fabrication of photonic devices is still costly and essentially time consuming, the development of science and technology has grown parallel to the development of accurate numerical techniques able to perform “modelling” of the devices. Besides analysis of the structures, the numerical simulation of the performance of the device makes it possible to perform operations of optimisation and design at an early stage with significant cost savings. The field of numerical modelling has become so important that developing the innovative potential of optics and photonics relies today on sophisticated simulation techniques.

Up to now, there have been many books in the market that deal with numerical modelling techniques. However, the majority only concentrate on one or two numerical methods. This book, on the other hand, covers a comprehensive state-of-the art of many modelling possibilities, ranging from the most innovative techniques in the frequency domain (e.g. the bidirectional beam propagation method), to ones in time domain (e.g. the multiresolution time domain and finite volume time domain methods), for future photonic devices.

Through an extensive number of simulation models and solutions for generic and more specific problems, this book focuses on equipping the reader with the sophisticated concepts of computational modelling in a new easy way to build his/her own codes. Furthermore, it comes with a CD that shows samples of codes, and some examples to enable the reader to understand various computational modelling cases.

This book consists of 11 chapters.

The introduction to the book is given in Chapter 1, where a general overview of modern photonics and recent advances in computational modelling for new photonic devices are given.

From Chapter 2 to Chapter 4, beam propagation method (BPM)-based numerical techniques in the frequency domain are considered. Among them, Chapter 2 presents the governing equations for the full-vectorial BPM, including the finite-element analysis, the perfectly matched layer (PML) scheme for the treatment of boundary conditions and the imaginary-distance BPM that provides the mode solver. Full-vectorial BPM is then assessed in Chapter 3, where modal analysis of rectangular waveguides is given, together with analysis of photonic crystal fibres and liquid-crystal-based photonic crystal fibres.

The bidirectional BPM is presented in Chapter 4. The chapter focuses on the optical waveguide discontinuity problem, giving the formulation of a numerical method able to perform its simulation. After deriving the governing equations, the bidirectional BPM is assessed in several optical examples.

From Chapter 5 onwards, numerical techniques belonging to the area of the time-domain schemes are presented. Chapter 5 starts with a recent innovative modification of the conventional finite-difference time-domain (FDTD) method: the complex-envelope alternating-direction-implicit FDTD. After derivation of the fundamental equations to solve Maxwell’s equations, assessment of the technique is given in the specific problem of photonic crystal cavities.

In Chapter 6, the finite-volume time-domain method is proposed as a novel alternative technique to FDTD for the study of electromagnetic problems. The scheme is presented in detail, providing mathematical formulations, analysis of numerical stability and dispersion, and an efficient scheme for the treatment of the boundary conditions. The FVTD is then assessed in Chapter 7 where both linear and nonlinear devices are investigated.

Further modelling possibilities are given in Chapter 8 with the multiresolution time domain (MRTD) method. Basic concepts of the multiresolution analyis are given and its application to the solution of the Maxwell’s equations is explained. An accurate and innovative extension of the method to the analysis of second-order nonlinear effects is given. Assessment of the MRTD scheme is presented in Chapter 9, providing code validation in linear photonic devices. Assessment of the technique for the study of second-order nonlinear photonic devices follows in Chapter 10. Chapter 11 shows improvement of the nonlinear-MRTD scheme with the inclusion of auxiliary differential equations that enable the accurate analysis of the linear dispersion of the media. Chapter 11 also presents validation of the code with the generation of second harmonics in a planar waveguide and in a one-dimensional photonic crystal.

Salah Obayya

Acknowledgements

The author would like to acknowledge the people whose help and efforts played an important role in making this book a reality. In particular, I express sincere gratitude to some of my group members; namely Dr. Rosa Letizia, Dr. Domenico Pinto and Dr. Mohamed Farahat.

I wish to mention also the students I taught through the years. Their questions inspired me to undertake this project in the first place.

The author would like to thank the John Wiley team, Sarah Tilley, Jasmine Chang and Jo Tyszka for their patience and support throughout this project.

Last but not least, the author would like to express his gratitude to Prof H. A. Elmikati, Dean of Faculty of Engineering, Mansoura University, whom I owe a lot.

1

Introduction

1.1 Photonics: The Countless Possibilities of Light Propagation

Following the advances in semiconductor physics that have allowed us to fully exploit the conducting properties of certain materials, thereby initiating the transistor revolution in electronics, in the last few decades a new frontier has opened up. The goal in this case is to control the optical properties of materials. An enormous range of technological developments would become possible if we could engineer materials that respond to light waves over a desired range of frequencies by perfectly reflecting them, or allowing them to propagate only in certain directions, or confining them within a specified volume.

Optical solutions to engineering problems are being increasingly found in many fields of application, including medicine, communication, entertainment, sensing and homeland security. Already, fibre optic cables, which simply guide light, have revolutionised the telecommunications industry. Laser engineering, high-speed computing and spectroscopy are just a few of the fields next in line to reap the benefits of the advances in optical materials. Photonics as a field covers today a huge range, from communications to science and technology applications, including laser manufacturing, biological and chemical sensing, display technology and optical computing.

For a number of reasons, ranging from a higher bandwidth to the inexpensiveness of optical fibre materials, components based on electric currents are more and more complemented or replaced by technology based on light. Relatively cheap, compact and highly sensitive photonic devices have already been commercialised for a variety of areas, and are expected to be basic blocks for all-optical integrated circuits, truly revolutionising the way we live.

All-optical systems are believed to be the best solution for the realisation of devices capable of meeting the high-performance characteristics required today and by next-generation telecommunications. In recent decades, with the advent of photonic crystals (PhCs) technology, the goal of a system capable of all-optical signal processing has made a giant step towards becoming a reality. PhC technology has brought the possibility of manipulating the light in a way that is basically not possible with conventional optical technology. Moreover, this capability of manipulating light is also obtained with an efficiency that is far greater than that obtained with conventional optical devices. However, the complex geometry of devices realised with PhC technology is still a challenge for fabrication technology, although significant progress has been made in this field, making it possible to fabricate new devices exploiting such technology. On the other hand, the mass production of these devices is still far from becoming a reality, and this can pose a serious threat to the development of PhC technology and, most of all, to the development of all-optical systems for ultra-fast telecommunication applications. In this aspect, numerical methods can be of vital help. Generally speaking, the possibility to predict the performance of a device before its practical fabrication can be a key factor towards the development of new and innovative technologies. With the exponential growth of the computational capabilities of modern personal computers (PCs), numerical methods have experienced a sudden expansion in all fields of engineering. This has made possible the development of new applications at a rate that was unthinkable just few decades ago.

Up to now, there have been many books on the market that treat numerical modelling techniques. However, the majority only concentrate on one or two numerical methods. The main goal of this book is to present a comprehensive state-of-the art of the latest modelling possibilities for the analysis and design of innovative devices for next-generation optical communications. Through a thorough analysis of the latest advances in numerical methods, weaknesses and drawbacks will be addressed in order to develop new techniques capable of dealing with simulations of optical devices in an innovative fashion, so as to increase either the accuracy or the efficiency. This aim will be reached through a deep theoretical analysis which will involve the physics behind the functionality of the optical device more directly in the core of the numerical techniques. In this way, the obtained numerical tools will be considered to be more state-of-the-art computational techniques than mere ‘number-crunching’ algorithms for the solution of Maxwell’s equations. The gain in terms of efficiency can be usefully employed for the extension of the methods to dealing with three-dimensional (3D) problems without requiring a prohibitive amount of computational resources, making it possible to obtain simulation results closer to reality.

The information in the book will be essential for the reader to perform and excel in the following:

Understanding the physics of light propagation in various photonic devices.Understanding Maxwell’s equations, governing the propagation phenomena.Grasping basic concepts of general numerical modelling techniques.Deriving the equations underlying the numerical modelling techniques.Applying the gained knowledge of computational modelling to modelling photonic devices.

These functions range from introductory (understanding the basic concepts) to advanced levels (application of the concepts to the numerical modelling of photonic devices).

1.2 Modelling Photonics

The last decade has witnessed dramatic progress and interest in micro-and nanofabrication techniques of complex photonic devices. In almost all cases, an accurate quantitative theoretical modelling of these devices has to be based on advanced computational techniques that solve the corresponding, numerically very large linear, nonlinear or coupled partial differential equations.

Photonics is especially suitable for computation because Maxwell’s equations are practically exact, the relevant material properties are well known and the length scales are not too small. Therefore, an exciting aspect of this field is that quantitative theoretical predictions can be made from first principles, without any questionable assumptions or simplifications. The results of such computations have consistently agreed with experiments. This makes it possible and preferable to optimise the design of photonic devices on a computer at an early stage, prior the actual fabrication. The computer becomes the pre-laboratory.

1.2.1 History of Computational Modelling

Many standard numerical techniques for the solution of partial differential equations have been applied to electromagnetics (EM), and each has its own particular strengths and weaknesses. High-quality ‘black-box’ software is widely available, including free, open-source programs. Indeed, computational photonics has matured so much that many are familiar only with the general principles and capabilities of the different tools.

In the last decades, a wide variety of numerical techniques have been developed and utilised for the design and optimisation of optical devices. Most of them have been successfully employed for the simulation of novel devices, mostly in two-dimensional (2D) space domains. Despite their powerful numerical capabilities, all these techniques possess drawbacks which pose limitations in the range of applicability of the methods to specific classes of problems. These drawbacks become more severe once the extension of these techniques to full-vectorial 3D space domains is considered. These limitations are mainly due to the huge growth in the computational burden required for 3D problems. Furthermore, if nonlinear phenomena are also taken into account, the limitations appear to be even more stringent.

Typically, the existing different numerical schemes can be divided into two main categories: time-domain and frequency-domain methods.

1.2.2 Frequency Domain

Before 1960, the principal approaches in the area of frequency-domain (FD) based numerical techniques involved closed-form and infinite-series analytical solutions, with numerical results from these analyses obtained using mechanical calculators. After 1960, the increasing availability of programmable electronic digital computers permitted such FD approaches to increase significantly in sophistication. Researchers were provided with a whole new range of capabilities offered by powerful high-level programming languages, such as Fortran, rapid random-access storage of large arrays of numbers, and computational speeds orders of magnitude faster than possible with mechanical calculators. In this period, the principal computational techniques for Maxwell’s equations included high-frequency asymptotic methods and integral equations. However, these FD methods have some difficulties and drawbacks. For instance, while asymptotic analyses are well suited for modelling the scattering properties of large electrical shapes, such approaches are find it difficult to deal with nonmetallic material compositions and the volumetric complexity of a structure. On the other hand, integral equation methods can deal with material and structural complexity, however their need to construct and solve systems of linear equations limits the electrical size of possible models, especially those requiring detailed treatment of geometric details within a volume.

Although significant progress has been made in solving the ultra-large systems of equations generated by these FD integral equations, the capabilities of even the latest of such technologies cannot keep up with many volumetrically complex structures of recent engineering interest. This also holds for FD finite-element techniques, which generate sparse rather than dense matrices. Moreover, properties of material such as nonlinearities cannot be easily incorporated into the FD solutions of Maxwell’s equations, which is a severe constraint, as research today is very active in the fields of active electromagnetic/electronic and electromagnetic/quantum-optical systems, such as high-speed digital circuits, and microwave and millimetre-wave amplifiers and lasers.

1.2.3 Time Domain

Since the arrival of the digital computer, which has profoundly changed the possibilities, time-domain (TD) modelling has offered efficient and flexible techniques to study computational electromagnetic propagation in linear and nonlinear optics.

There are two basic reasons for the success of TD over FD modelling: computational efficiency and problem requirements. Generally, when broadband information is analysed, a TD approach is intrinsically a more immediate choice because it provides a transient response whose bandwidth is limited only by the frequency content of the source, and the time and space sampling adopted in the numerical approach. Moreover, the computational efficiency of TD is also derived from its natural ability to adapt to parallel computer architectures. In addition, the TD approach can usually model problems involving time-varying media and components in a more straightforward way.

A representative example of the rapid growth in TD research is the popularity of the finite-difference time-domain (FDTD) method developed by Taflove.

There are general steps that need to be carried out in order to perform numerical modelling in the time domain:

Develop time-dependent integral or Maxwell’s curl equations.Discretise the equations in space and in time by means of an appropriate grid in space, and suitable basis and testing functions.Derive a set of equations that relate unknown with known quantities (starting from an initial value that usually is given by the source field).Generate a numerical solution of this initial-value problem in space and time.

1.2.4 Chapter Overview

Following this Introduction, from Chapter 2 to Chapter 4, beam propagation method (BPM) based numerical techniques in the frequency domain are considered. Amongst them, Chapter 2 presents the governing equations for the full-vectorial BPM, including the finite-element analysis, the perfectly matched layer (PML) scheme for the treatment of boundary conditions and the imaginary-distance BPM that provides the mode solver. Full-vectorial BPM is then assessed in Chapter 3, where modal analysis of rectangular waveguides is given together with analysis of photonic crystal fibres and liquid-crystal-based photonic crystal fibres.

The bidirectional BPM is presented in Chapter 4. The chapter focuses on the optical waveguide discontinuity problem, giving the formulation of a numerical method able to perform its simulation. After deriving the governing equations, the bidirectional BPM is assessed in several optical examples.

From Chapter 5, numerical techniques belonging to the area of time-domain schemes are presented. Chapter 5 starts with a recent innovative modification of the conventional finite-difference time-domain (FDTD) method: complex-envelope alternating-direction-implicit FDTD. After derivation of the fundamental equations to solve Maxwell’s equations, assessment of the technique is driven by the specific problem of photonic crystal cavities.

In Chapter 6, the finite-volume time-domain (FVTD) method is proposed as a novel alternative technique to FDTD for the study of electromagnetic problems. The scheme is presented in detail, providing mathematical formulations, analysis of numerical stability and dispersion, and an efficient scheme for the treatment of the boundary conditions. The FVTD is then assessed in Chapter 7, where both linear and nonlinear devices are investigated.

Further modelling possibilities are given in Chapter 8, with the multiresolution time domain (MRTD) method. Basic concepts of the multiresolution analysis are given and its application to the solution of Maxwell’s equations is explained. An accurate and innovative extension of the method to the analysis of second-order nonlinear effects is given. Assessment of the MRTD scheme is presented in Chapter 9, providing code validation in linear photonic devices. Assessment of the technique for the study of second-order nonlinear photonic devices follows in Chapter 10. Chapter 11 shows improvement of the nonlinear-MRTD scheme, with the inclusion of auxiliary differential equations (ADE) that enable the accurate analysis of linear dispersion of the media. Chapter 11 also presents validation of the code with generation of second harmonics in a planar waveguide and in a one-dimensional (1D) photonic crystal.

1.2.5 Overview of Commercial Software for Photonics

The continuous advancement of microwave circuits and electromagnetic devices towards increased functionality and performance requires simultaneous development of modelling tools that are able to keep up with the growing level of sophistication. In order to get a deeper insight into the field of computational electromagnetics (CEM) and grasp what is available in the market, the most widely used commercial packages are listed as follows:

COMSOL Multiphysics (formerly FEMLAB) is a finite element analysis, solver and simulation software package for various physics and engineering applications, especially coupled phenomena, or multiphysics. COMSOL Multiphysics also offers an extensive interface to MATLAB and its toolboxes for a large variety of programming, pre-processing and post-processing possibilities.FIMMWAVE is a generic full-vectorial mode-finder for waveguide structures. FIMMWAVE combines both methods based on semi-analytical techniques with other more numerical methods such as finite difference or finite element.CST MICROWAVE STUDIO® (CST MWS) offers five solver modules; the Transient, Eigenmode, Frequency Domain, ‘Resonant: Fast S-Parameter’, ‘Resonant: S-Parameter, Fields’ (formerly known as Modal Analysis), and the Integral Equation Solver, each offering distinct advantages in their own domains. There are numerical advantages offered by the method used in most of the solvers, the finite integration technique (FIT).CrystalWave is a design environment for the layout and design of integrated optics components optimised for the design of photonic crystal structures. It is based on both FDTD and finite-element frequency-domain (FEFD) simulators and includes a masque file generator carefully optimised for planar photonic crystal structures.Optiwave is a suite of engineering design tools. Amongst these, there are OptiFDTD, which is based on the FDTD algorithm with second-order numerical accuracy and the most advanced boundary condition – the uniaxial perfectly matched layer (UPML) boundary condition, and OptiBPM, which is based on the BPM.RSoft’s Photonic Component Design Suite allows the design and the simulation of both passive and active photonic devices for optical communications, optoelectronics and semiconductor manufacturing. FullWAVE is a simulation tool for studying the propagation of light in a wide variety of photonic structures, including integrated and fibre-optic waveguide devices, as well as circuits and nanophotonic devices, such as photonic crystals. The software employs the FDTD method for the full-vector simulation of photonic structures. BandSOLVE is a design tool for the calculation of photonic band structures for all photonic crystal (PC) devices which employs the plane wave expansion (PWE) algorithm.

Indeed, commercial EM simulation programs available today provide powerful design tools, however no single numerical method provides a universal solution, and commercial codes often fail to accurately simulate high-end problems found in cutting-edge research. Therefore, research in computational EM is still essential to keep up with the increasing complexity of devices throughout the EM spectrum and has to develop parallel to fabrication technologies to allow the full establishment of the new photonic solutions in the market.

2

Full-Vectorial Beam Propagation Methods

2.1 Introduction

Numerical simulations play an important role for the design and modelling of guided-wave optoelectronic devices. There are various modelling methods in which not only a full-vector model, but also an approximate scalar model, are used. In this chapter, an overview of beam propagation methods (BPMs) [1–4] is introduced. In addition, the formulation of the wave equations in terms of the electric and magnetic fields is included. Moreover, this chapter includes an introduction to finite-element analysis, followed by the derivation of the finite-element BPM. Also, the formulation of the imaginary-distance full-vectorial finite element BPM scheme proposed in [5] is extended further to fully treat the vectorial complex modes.

2.2 Overview of the Beam Propagation Methods

Analysis and simulation of electromagnetic wave propagation are essential in the modelling and design of optical waveguide devices. The BPM [1–4] has been one of the most popular techniques for modelling and simulation of such optical devices. The major concept of the BPM is the development of a formula that permits the propagation of an initial field distribution along the axial direction by steps of sufficiently small length, as shown in [4].