FAST BPM METHOD TO ANALYZE NONLINEAR OPTICAL DEVICES
|
|
- Abraham Stafford
- 5 years ago
- Views:
Transcription
1 FAST BPM METHOD TO ANALYZE NONLINEAR OPTICAL DEVICES J. de Oliva Rubio 1 and I. Molina Fernandez 1 1 Departamento de Ingenierıa de Comunicaciones Escuela Tecnica Superior de Ingenieros de Telecomunicacion Universidad de Malaga 9071 Malaga, Spain Recei ed February 1998 ABSTRACT: A new finite-difference BPM method based on the application of the Runge Kutta ( RK) technique altogether with transparent boundary conditions ( TBCs) is proposed. Numerical experiments are carried out to erify the compatibility of TBCs with this method, and to compare its performance with the standard Crank Nicholson scheme when simulating nonlinear optical de ices John Wiley & Sons, Inc. Microwave Opt Technol Lett 18: , Key words: nonlinear beam-propagation method; nonlinear optical wa eguide; transparent boundary conditions; Runge Kutta method I. INTRODUCTION Simulation of integrated optical devices is commonly performed by means of scalar beam-propagation methods Ž BPM. based upon the resolution of the Fresnel wave equation 1. Although, historically, BPM techniques based on the FFT were first introduced and exploited, today, finite-difference schemes Ž FDBPM. are the most widely used due to their superior ability to analyze strongly guiding and nonlinear structures. To date, the most popular FDBPM method is based upon the discretization of the scalar Fresnel equation by means of the implicit Crank Nicholson Ž CN. finite-difference scheme 1 together with the transparent boundary conditions Ž TBCs.. This method, which will be referred to as CN FDBPM in this paper, leads to a sparse nonlinear system of equations which must be solved in each propagation step. Explicit finite-difference BPM methods have also been developed 3, 4 which avoid the time-consuming algorithms for the solution of the nonlinear system of equations. These methods are not unconditionally stable and, to our knowledge, they have not been used in conjunction with TBCs, so absorbent boundary conditions Ž ABCs. were used, which considerably increase the numerical magnitude of the problem and reduce the method efficiency. In this paper, a fourth-order Runge Kutta Ž RK. explicit method is developed for the solution of the Fresnel equation which allows TBCs to be successfully applied, thus drastically reducing the computational effort. This method relies on the discretization of the transversal operator with a second-order centered finite-difference scheme Žthe same as in the CN algorithm., leading to a system of nonlinear ordinary differential equations in the longitudinal variable, which is then solved by means of the RK technique. As this method employs a finite-difference discretization for the transversal operator and an RK strategy for the longitudinal one, it will be called in this paper the RK FDBPM technique. Although, for the sake of simplicity, only the scalar problem has been analyzed in this paper, since the proposed method is transversally identical to the standard CN FDBPM technique, its extension to already existing semivectorial and vectorial formulations can easily be performed. The RK method has been previously used for the resolution of the nonlinear Fresnel equation 5, but, to our knowledge, this is the first time that this method has been applied in conjunction with TBCs, and that a detailed comparison between RK and CN FDBPM has been performed. Assessment of the proposed technique is carried out in two steps. First, to study the behavior of TBCs within the RK FDBPM algorithm, two previously proposed problems for testing the boundary conditions are analyzed: diffraction and oblique incidence on a boundary of a Gaussian beam propagating through a linear homogeneous medium. Second, the ability of the RK FDBPM technique to analyze nonlinear structures is established and compared with the standard CN FDBPM. This is carried out by simulating the propagation of a three-dimensional Gaussian beam in an optical slab with a nonlinear Kerr-type saturable substrate 6, and studying the functioning modes of a nonlinear directional coupler 7. A rigorous comparison of the CN FDBPM and RK FDBPM methods is obtained by choosing mesh sizes which guarantee the same numerical dispersion in both cases 8. II. THE NONLINEAR CN FDBPM METHOD The application of the finite differences implicit in CN scheme plus the TBCs to the resolution of the scalar Fresnel equation produces a nonlinear system of equations that must be solved at every propagation step. The fastest way to solve it is to linearize the system using the known value of the field at every point of the previous propagation step to estimate the nonlinear refraction index in the actual propagation step. If better accuracy is desired, it is possible to iterate, using the new field obtained, to make a better estimation of the nonlinear refraction index and recalculate the field. In the described technique, a system of linear equations must be solved at every iteration of the nonlinear method. The system solution may be achieved by LU decomposition, but this is a very expensive method in terms of memory and CPU time, so it is preferable to exploit the sparse characteristic of the system matrix using, for example, the Gauss Seidel Ž GS. method 9. This method is the one that has been used in this paper to perform the CN FDBPM algorithm. The use of the Crank Nicholson scheme together with the Gauss Seidel method, for the usual mesh sizes in nonlinear optical problems, reduces the memory requirements at least two orders of magnitude and the consumed CPU time one order of magnitude, depending upon the accuracy imposed on the GS method. In most cases of interest, one single iteration of the non-linear method is enough to get sufficiently accurate results. The reason is that the nonlinear coefficient of the usually employed materials is very small, as happens with the propagation step, so that the error made in the estimation of the nonlinear coefficient by means of the field in the previous propagation step is very small. Regardless of the improvement introduced by the GS method, the need to solve at least one linear system at every propagation step considerably increases the computation time, especially when very fine meshes are used to reduce numerical dispersion. III. THE NONLINEAR RK FDBPM METHOD WITH TBCS The computational effort of the previously described technique can be reduced if an explicit scheme is employed; however, the stable nonlinear explicit methods proposed to 418 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
2 date have the disadvantage of being incompatible with the TBCs 4, 8. Absorbent boundary conditions considerably increase the transversal magnitude of the numerical problem, and thus the CPU time and memory needed. To have an explicit stable nonlinear method compatible with the TBCs, the application of the fourth-order Runge Kutta method to the solution of the scalar Fresnel equation has been developed. Under the slowly varying approximation, the propagation of optical beams in a nonlinear Kerr saturable-type waveguide is governed by the Fresnel equation A A A Ž. jk N k n x, y, z, A N 0 0 A z x y Ž 1. Ž. where n x, y, z, E is the nonlinear refractive index given by A L NL n Žx, y, z, E. n Ž x, y, z. C Ž x, y, z. Ž. 1 A and A ' NL E is the envelope of the transversal electric field normalized to the medium nonlinear coefficient. In Ž., nl is the linear part of the refractive index, CNL accounts for the spatial distribution of the nonlinear Kerr saturable effect, and is 1 for saturable media and 0 for nonsaturable ones. The first step in the developed technique is to discretize the transversal partial derivatives using a second-order centered finite-difference scheme. In doing so, the Fresnel differential equation is transformed into the following system of nonlinear ordinary differential equations: da 1 A A A dz jk0 N x p, q p 1, q p, q p 1, q A A A p, q 1 p, q p, q 1 y Ž. k n N A Ž 3. 0 p, q p, q which is well suited to the application of the fourth-order Runge Kutta method. To do it, the transversal mesh points are consecutively numbered from left to right and from down to up; the resulting transversal operator is labeled FŽz, AŽ z.., A A 1, 1, A 1,,..., A 1, Nx,..., A, 1,..., A, Nx, A Ny,1, A Ny,,..., A Ny, Nx being the vector whose components A k are the values Ap, q arranged as described, i.e., with k p qn x, Nx being the width and Ny the height of the transverse mesh. In doing so, Eq. Ž. 3 becomes Ž. da z dz where the transversal operator is Ž Ž.. Ž. F z,a z 4 1 FŽ z,a. MŽz, A. A, Ž 5. jk N 0 M being the pentadiagonal matrix obtained from the transversal discretization. For the mesh points interior to the contour, its elements are 1 1 Mk, k 1, Mk, k Nx, x y Ž. Mk, k k0 np, qž z, Ap, q z. N, x y Ž 6. while for the mesh points at the contour, the field is calculated with the TBCs that relate the value of the field at any point in the contour with the value of the field at the adjacent point interior to the contour. For example, at any point on the right boundary, the TBC condition is Ž. Ž. A A exp jk x 7 Nx 1, q Nx q x where ANx 1, q is the value of the field at a point on the boundary, ANx q is the field at the adjacent left point, and the value of k x is calculated from the previous step field at the points Ž N 1, q. and Ž N, q. from the expression x x Ž. Ž. A A exp jk x. 8 Nx, q Nx 1, q x To get the field AŽ z z. from AŽ z. 0 0, the RK algorithm is applied; notice that the evaluation of the right-hand side of the nonlinear system of ordinary differential equations is easily performed by multiplying a sparse matrix Ž M. by a vector, Ž A. which is a fast operation, involving fewer than 5NxNy multiplications. IV. NUMERICAL RESULTS Numerical simulations have been carried out to study the behavior of the TBCs within the RK FDBPM algorithm and the performance of the method when analyzing nonlinear structures. The TBCs behavior has been tested by comparing simulation results with analytical solutions for the diffraction, and oblique incidence on the computational boundary, of a Gaussian beam propagating through a linear homogeneous medium. The mesh size in both transverse coordinates and the propagation step for these experiments are shown in the figure captions. In Figure 1, analytical and numerical results for the diffraction of a Gaussian beam propagating in the z-direction are plotted together, and it can be seen that they are almost indistinguishable. Figure shows a comparison of the normalized power in the computational window obtained by the numerical and analytical methods. These results, which are comparable to those obtained in, show that the behavior of the RK FDBPM dealing with the outflow of the beam power by the four boundaries at the same time is excellent; furthermore, the behavior is fairly good even when the fraction of power remaining in the computational window is very small. In Figure 3, the incidence on a boundary of a Gaussian beam propagating with a tilt angle of.5 with respect to the z-axis is shown and compared with analytical results for different propagation distances. It can be seen that the beam, which propagates almost diffractionless, passes through the boundary almost completely without distortion. Figure 4 shows a comparison between the numerically and analytically calculated power of the beam in the computational window MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
3 Figure 1 Gaussian beam diffraction in a linear homogeneous medium, 1550 nm, x y z 1 m, N 1, spot size 0.5 m. numerical RK FDBPM with TBCs, analytical Figure Variation with distance of the normalized power in the computational window for Gaussian beam diffraction with the same data as Figure 1. analytical, numerical as a function of the propagation distance. It can be seen that very good agreement between both results is obtained up to the distance where 75% of the initial power has been radiated. After that distance, the numerically calculated power in the computational window is greater than it should be due to reflection from the TBCs. However, it must be noted that the numerically calculated power tends asymptotically to zero with the propagation distance, resulting in a reflection coefficient of , calculated as the remaining normalized power in the computational window when the beam has completely passed through the boundary Ž z 4000 m.. Obviously, the reflection coefficient obtained with this technique Ž TBCs and RK FDBPM. is not as good as the one obtained with the TBCs and the CN FDBPM for which they were originally developed, but, as will be shown in the next paragraphs, it can be enough for many practical situations. From previous results, the correct behavior of the TBCs when applied to the proposed RK FDBPM technique has been stated; however, where the proposed method shows its superior performance is when dealing with nonlinear structures where an accurate description of dispersion and nonlinearity must be achieved. To compare the behavior of the CN FDBPM and the RK FDBPM techniques when analyzing nonlinear structures, two previously reported problems have been simulated: the propagation of a three-dimensional Gaussian beam in an optic slab with nonlinear saturable Kerr-type substrate 6, and the study of the nonlinear directional coupler working modes 7. The excitation and slab s refractive indexes are the same as those of the references, and are reproduced in the figure captions for convenience. Figure 5 shows the result of propagation of a Gaussian beam in a slab with a nonlinear saturable substrate 6 for z 70. Figure 5Ž. a has been obtained by means of the CN FDBPM technique, while Figure 5Ž b. has been obtained with the RK FDBPM technique. In both cases, the mesh sizes were x 0.075, y 0.3, and z 0.01, which have been chosen to keep the numerical dispersion under % for at least 95% of the bandwidth of the Gaussian excitation 8. In this experiment, the computational effort of the CN FDBPM technique is twice the effort of the RK FDBPM Žthe CPU times on an HP workstation were 4h, 38 min for the CN FDBPM and, 44 min for the RK FDBPM.. Comparing both figures, it is clear that the CN solution presents a higher penetration of the beam in the substrate, while in the RK solution, a greater amount of energy remains trapped by the film. To find which of the previous results is better, the same simulation has been repeated with both techniques with reduced mesh sizes Žnumerical dispersion below 0.06%, x 0.05, y 0.05, and z 0.01., both methods yielding similar solutions. The RK FDBPM technique presented in Figure 6 Ž CPU time 1 h, 50 min., showing good agreement with previously published data 6. As compared with 40 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
4 Figure 3 Gaussian beam propagation through a TBC, 1550 nm, x y z 1 m, N 1, spot size 10 m. numerical Ž RK FDBPM with TBCs., analytical Figure 4 Variation with distance of the normalized power in the computational window for Gaussian beam propagation through a TBC with the same data as in Figure 3. analytical, numerical the plots of Figure 5Ž. a and Ž. b, it can be stated that, when dealing with nonlinear structures, the RK FDBPM technique gives more accurate results and is more computationally efficient than the standard CN FDBPM technique with equal mesh sizes. This result is due to the fact that the nonlinearity is solved with fourth-order precision in the RK algorithm, while it is only second order in the CN method. To further validate the performance of the proposed technique when dealing with nonlinear structures, a directional coupler with a nonlinear substrate 7 has been analyzed. In Figures 7Ž. a and Ž. b, propagation in this device, when excited with a Gaussian field profile closely resembling the mode trapped in a single arm of the coupler, is presented. Figure 7Ž. a corresponds to simulation results obtained by means of the CN FDBPM technique Ž CPU time 17 h, 1 min., while Figure 7Ž b. has been obtained with the RK FDBPM technique Ž CPU time 3 h, 6 min.. Transversal mesh sizes have been fixed to k0 x 0.1, k0 y 0.1, k0 z in both cases, which assures that the numerical dispersion is Figure 5 Results of the propagation of a Gaussian beam in a nonlinear slab at z 70. Ž. a CN FDBPM with TBCs. Ž. b RK FDBPM with TBCs. x 0.075, y 0.3, z 0.01, A , x m, w0 197 m, a 1, N 1.644, d Excitation: Ax,y,z 0 A Ž. exp Ž x x. Ž y a. w Refractive index: n c.647 x 0 Ž. n x d n.647 A Ž1 A. s x d n x, y, A f MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
5 Figure 6 Gaussian beam propagation in a slab with nonlinear substrate. RK FDBPM with reduced mesh size: x 0.05, y 0.05, z Excitation and refraction index are the same as those in Figure 5. under % for 95% of the excitation bandwidth. When compared with previously published results Ž 7 reprinted in Fig. 8 with permission., it can easily be seen that the RK FDBPM exhibits much better agreement than CN FDBPM. As the numerical dispersion has been kept very low and the same for both simulations, the better performance of the RK technique can only be explained as a consequence of the nonlinear fourth-order precision of the RK technique. V. CONCLUSIONS An RK FDBPM method to analyze nonlinear 3-D optical devices under the scalar slowly varying envelope approximation has been presented, and its performance has been compared with the standard Crank Nicholson technique. The method is based on the finite-difference discretization in the transversal coordinates and on the application of the fourthorder Runge Kutta technique in the propagation coordinate, and this allows transparent boundary conditions Ž TBCs. to be included. Two different sets of numerical experiments have been carried out: 1. to determine the compatibility of the RK FDBPM method with the TBCs,. to establish the ability of the proposed method to deal with nonlinear devices. From the first set of experiments, it is shown that, although the reflection coefficient of the TBCs obtained with the proposed technique is not as good as the one attainable within the CN FDBPM, the TBCs in conjunction with the RK FDBPM lead to a continuous outflow of energy from the computational window, giving a reflection coefficient that can be sufficiently small for many practical situations. On the other hand, the reduction of the size of the computational window allowed by the use of TBCs greatly increases the numerical efficiency of the proposed RK FDBPM algorithm. The second set of experiments clearly evidences the greater accuracy and computational efficiency of the RK FDBPM Figure 7 Gaussian beam propagation through the left arm of a nonlinear directional coupler. Ž. a CN FDBPM with TBCs. Ž. b RK FDBPM with TBCs. k0 x 0.1, k0 y 0.1, k0 z 0.003, N 0.5, k0 1, A , k0x0 1.5, k0 0.7, k , k d 1.5, k a 1. Excitation: Ax,y,z 0A Ž. exp Ž x x. 4 y 4. Refractive index: n Žx, y, A. ½.5 A d a x d a 0.5 A d a x d a 4 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
6 Propagation Method, IEEE J. Quantum Electron., Vol. 8, Jan. 199, pp Y. Chung and N. Dagli, An Explicit Finite Difference Beam Propagation Method: Application to Semiconductor Rib Waveguide Y-Junction Analysis, Electron. Lett., Vol. 6, No. 11, 1990, pp Y. Chung and N. Dagli, Analysis of z-invariant and z-variant Semiconductor Rib Waveguides by Explicit Finite Difference Beam Propagation Method with Nonuniform Mesh Configuration, IEEE J. Quantum Electron., Vol. 7, Oct. 1991, pp H. Liu and W. Wang, Beam Propagation Analysis of the Non- Linear Tapered Optical Waveguide, IEEE Microwa e Guided Wa e Lett., Vol. 5, Feb. 1995, pp D. Mihalache, D. M. Baboiu, D. Mazilu, L. Torner, and J. P. Torres, Gaussian-Beam Excitation and Stability of Three-Dimensional Non-Linear Guided Waves, J. Opt. Soc. Amer. B, Vol. 11, No. 7, 1994, pp A. B. Aceves, A. D. Capobianco, B. Constantini, C. De Angelis, and G. F. Nalesso, Beam Dynamics in Non-Linear Coupled Slab Waveguides: Three-Dimensional Variational Analysis, J. Opt. Soc. Amer. B, Vol. 11, July 1994, pp J. de Oliva Rubio and I. Molina Fernandez, Tecnicas de Diferencias Finitas para el Analisis de Estructuras Opticas Tridimensionales No-Lineales, Actas del XI Simp. Nacional de la Union Cientıfica Internacional del Radio Ž URSI., Madrid, Spain, Sept. 1996, Vol., pp W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes: The Art of Scientific Computing, Cambridge University Press, London John Wiley & Sons, Inc. CCC NUMERICAL DISPERSION IN THE FINITE-ELEMENT METHOD USING THREE-DIMENSIONAL EDGE ELEMENTS Figure 8 Results obtained in 7 for the Gaussian beam propagation through the left arm of a nonlinear directional coupler. Reprinted with permission of the OSA and the authors. when analyzing nonlinear devices, which can lead to an 80% reduction of the computational time with respect to the standard CN FDBPM technique. From the presented results, it can be stated that the proposed technique can be a very fast and inexpensive simulation tool for the analysis and design of slightly powerradiating 3-D nonlinear optical devices with small requirements on the reflection coefficient from the numerical boundaries. ACKNOWLEDGMENT The authors wish to thank the Optical Society of America and the authors of 7 for their permission to reprint Figure 8, and especially Prof. Capobianco for providing them with high-quality PostScript copies of the figure. This work was supported by the Spanish CICYT, Project TIC C REFERENCES 1. H. P. Nolting and R. Marz, Results of Benchmark Tests for Different Numerical BMP Algorithms, J. Lightwa e Technol., Vol. 13, Feb. 199, pp G. R. Hadley, Transparent Boundary Condition for the Beam Gregory S. Warren 1 and Waymond R. Scott, Jr. 1 1 School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta, Georgia Recei ed 8 January 1998 ABSTRACT: The discretization inherent in the ector finite-element method results in the numerical dispersion of a propagating wa e. The numerical dispersion of a time-harmonic plane wa e propagating through an infinite, three-dimensional, finite-element mesh composed of hexahedral and tetrahedral edge elements is in estigated in this work. The effects on the numerical dispersion of the propagation direction of the wa e and the electrical size of the elements are in estigated. The numerical dispersion of the tetrahedral edge elements is found to be dependent upon the polarization of the plane wa e propagating through the mesh. In addition, the dispersion of the tetrahedral elements is significantly smaller than the dispersion of the hexahedral edge elements. Both elements are found to ha e a phase error that con erges at the rate of O[( h ) ] John Wiley & Sons, Inc. Microwave Opt Technol Lett 18: 43 49, Key words: finite-element method; numerical dispersion; wa e propagation I. INTRODUCTION The finite-element method is a popular technique in computational electromagnetics for solving three-dimensional vector field problems. For this type of application, edge ele- MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 18, No. 6, August
THE beam propagation method (BPM) is at present the
JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 17, NO. 11, NOVEMBER 1999 2389 Three-Dimensional Noniterative Full-Vectorial Beam Propagation Method Based on the Alternating Direction Implicit Method Yu-li Hsueh,
More informationEVALUATION of the modal characteristics of an optical
JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 23, NO. 5, MAY 2005 1947 Implicit Yee-Mesh-Based Finite-Difference Full-Vectorial Beam-Propagation Method Junji Yamauchi, Member, IEEE, Member, OSA, Takanori Mugita,
More informationModeling liquid-crystal devices with the three-dimensional full-vector beam propagation method
214 J. Opt. Soc. Am. A/ Vol. 23, No. 8/ August 26 Wang et al. Modeling liquid-crystal devices with the three-dimensional full-vector beam propagation method Qian Wang, Gerald Farrell, and Yuliya Semenova
More informationNUMERICAL methods are indispensable tools for the
IEEE JOURNAL OF QUANTUM ELECTRONICS, VOL. 45, NO. 2, FEBRUARY 2009 117 A Numerical Approach for Full-Vectorial Analysis of 3-D Guided Wave Structures With Multiple and Strong Longitudinal Discontinuities
More informationBEAM propagation methods (BPM s) are very powerful
IEEE JOURNAL OF QUANTUM ELECTRONICS, VOL. 33, NO. 3, MARCH 1997 367 Efficient Semivectorial Mode Solvers Frank Wijnands, Thomas Rasmussen, Hugo J. W. M. Hoekstra, Jørn H. Povlsen, and René M. de Ridder,
More informationNumerical Assessment of Finite Difference Time Domain and Complex-Envelope Alternating-Direction-Implicit Finite-Difference-Time-Domain
Proceedings of the Federated Conference on Computer Science and Information Systems pp. 255 260 ISBN 978-83-60810-22-4 Numerical Assessment of Finite Difference Time Domain and Complex-Envelope Alternating-Direction-Implicit
More informationGaussian beam diffraction in inhomogeneous media of cylindrical symmetry
Optica Applicata, Vol. XL, No. 3, 00 Gaussian beam diffraction in inhomogeneous media of cylindrical symmetry PAWEŁ BERCZYŃSKI, YURI A. KRAVTSOV, 3, GRZEGORZ ŻEGLIŃSKI 4 Institute of Physics, West Pomeranian
More informationModeling microlenses by use of vectorial field rays and diffraction integrals
Modeling microlenses by use of vectorial field rays and diffraction integrals Miguel A. Alvarez-Cabanillas, Fang Xu, and Yeshaiahu Fainman A nonparaxial vector-field method is used to describe the behavior
More informationBeam propagation method for waveguide device simulation
1/29 Beam propagation method for waveguide device simulation Chrisada Sookdhis Photonics Research Centre, Nanyang Technological University This is for III-V Group Internal Tutorial Overview EM theory,
More informationSolitons. Nonlinear pulses and beams
Solitons Nonlinear pulses and beams Nail N. Akhmediev and Adrian Ankiewicz Optical Sciences Centre The Australian National University Canberra Australia m CHAPMAN & HALL London Weinheim New York Tokyo
More informationStable Implicit Scheme for TM Transient Scattering from 2D Conducting Objects Using TD-EFIE
Progress In Electromagnetics Research Letters, Vol. 73, 99 104, 2018 Stable Implicit Scheme for TM Transient Scattering from 2D Conducting Objects Using TD-EFIE Qiang Wang 1, 2, Lixin Guo 1, *,PengjuYang
More informationBeam Propagation Method Solution to the Seminar Tasks
Beam Propagation Method Solution to the Seminar Tasks Matthias Zilk The task was to implement a 1D beam propagation method (BPM) that solves the equation z v(xz) = i 2 [ 2k x 2 + (x) k 2 ik2 v(x, z) =
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Lecture No. # 09 Electromagnetic Wave Propagation Inhomogeneous Plasma (Refer Slide Time: 00:33) Today, I
More informationFinite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations
Finite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations Wolfgang Freude and Jan Brosi Institute of High-Frequency and Quantum Electronics (IHQ), University of Karlsruhe, Germany
More informationA new method for sensitivity analysis of photonic crystal devices
A new method for sensitivity analysis of photonic crystal devices Georgios Veronis, Robert W. Dutton, and Shanhui Fan Department of Electrical Engineering, Stanford University, Stanford, California 9435
More informationVector diffraction theory of refraction of light by a spherical surface
S. Guha and G. D. Gillen Vol. 4, No. 1/January 007/J. Opt. Soc. Am. B 1 Vector diffraction theory of refraction of light by a spherical surface Shekhar Guha and Glen D. Gillen* Materials and Manufacturing
More informationSoliton trains in photonic lattices
Soliton trains in photonic lattices Yaroslav V. Kartashov, Victor A. Vysloukh, Lluis Torner ICFO-Institut de Ciencies Fotoniques, and Department of Signal Theory and Communications, Universitat Politecnica
More informationPart 1: Fano resonances Part 2: Airy beams Part 3: Parity-time symmetric systems
Lecture 3 Part 1: Fano resonances Part 2: Airy beams Part 3: Parity-time symmetric systems Yuri S. Kivshar Nonlinear Physics Centre, Australian National University, Canberra, Australia http://wwwrsphysse.anu.edu.au/nonlinear/
More informationTM-Radiation From an Obliquely Flanged Parallel-Plate Waveguide
1534 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 50, NO. 11, NOVEMBER 2002 TM-Radiation From an Obliquely Flanged Parallel-Plate Waveguide Jae Yong Kwon, Member, IEEE, Jae Wook Lee, Associate Member,
More informationOptimized split-step method for modeling nonlinear pulse propagation in fiber Bragg gratings
448 J. Opt. Soc. Am. B/ Vol. 25, No. 3/ March 2008 Z. Toroker and M. Horowitz Optimized split-step method for modeling nonlinear pulse propagation in fiber Bragg gratings Zeev Toroker* and Moshe Horowitz
More informationOptical Nonlinearities at the Interface between Glass and Liquid Crystal 1
Optical Nonlinearities at the Interface between Glass and Liquid Crystal 1 PEDRO MENENDEZ-VALDES and JOSE A. MARTIN-PEREDA Universidad Politecnica de Madrid, Escuela Tecnica Superior de Ingenieros de Teiecomunicacion,
More informationAmbiguity of optical coherence tomography measurements due to rough surface scattering
Ambiguity of optical coherence tomography measurements due to rough surface scattering Y. Ashtamker, 1 V Freilikher, 1,* and J C Dainty 2 1 Department of Physics, Bar-Ilan University, Ramat Gan 52900,
More informationFAST METHODS FOR EVALUATING THE ELECTRIC FIELD LEVEL IN 2D-INDOOR ENVIRONMENTS
Progress In Electromagnetics Research, PIER 69, 247 255, 2007 FAST METHODS FOR EVALUATING THE ELECTRIC FIELD LEVEL IN 2D-INDOOR ENVIRONMENTS D. Martinez, F. Las-Heras, and R. G. Ayestaran Department of
More information1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationSource Free Surface x
Finite-dierence time-domain model for elastic waves in the ground Christoph T. Schroeder and Waymond R. Scott, Jr. School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta,
More informationGuided and Leaky Modes of Planar Waveguides: Computation via High Order Finite Elements and Iterative Methods
PERS ONLNE, VOL. 6, NO. 7, 669 Guided and Leaky Modes of Planar Waveguides: Computation via High Order Finite Elements and terative Methods D. Stowell and J. Tausch Southern Methodist University, USA Abstract
More informationTHE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION
Progress In Electromagnetics Research Letters, Vol. 5, 137 149, 2008 THE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION Y.-L. Li, M.-J.
More informationResearch Article Solution of the Porous Media Equation by a Compact Finite Difference Method
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2009, Article ID 9254, 3 pages doi:0.55/2009/9254 Research Article Solution of the Porous Media Equation by a Compact Finite Difference
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationEnhanced backscattering due to total internal reflection at a dielectric-air interface
JOURNAL OF MODERN OPTICS, 1995, VOL. 42, NO. 2, 257-269 Enhanced backscattering due to total internal reflection at a dielectric-air interface R. E. LUNAt, E. R. MENDEZ Centro de Investigacion Cientifica
More informationComputational Electromagnetics Definitions, applications and research
Computational Electromagnetics Definitions, applications and research Luis E. Tobón Pontificia Universidad Javeriana Seminario de investigación Departamento de Electrónica y Ciencias de la Computación
More informationComputer simulation of cylindrical laser beam self-focusing in a plasma
Computer simulation of cylindrical laser beam self-focusing in a plasma D. Subbarao a *, Karuna Batra a, b, Manik Bali c, Sugata Mitra b a Fusion Studies Program, Plasma Sc. and Tech. Group, Centre for
More informationVectorial Analysis of the Tapered Dielectric Waveguides
Proceedings of the 5th WSEAS Int. Conf. on Electronics, Hardware, Wireless and Optical Communications, Madrid, Spain, February 5-7, 006 (pp03-08) Vectorial Analysis of the Tapered Dielectric Waveguides
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More informationNear-field diffraction of irregular phase gratings with multiple phase-shifts
References Near-field diffraction of irregular phase gratings with multiple phase-shifts Yunlong Sheng and Li Sun Center for optics, photonics and laser (COPL), University Laval, Quebec City, Canada, G1K
More information4.4 Microstrip dipole
4.4 Microstrip dipole Basic theory Microstrip antennas are frequently used in today's wireless communication systems. Thanks to their low profile, they can be mounted to the walls of buildings, to the
More informationLight Localization in Left-Handed Media
Vol. 112 (2007) ACTA PHYSICA POLONICA A No. 4 Proceedings of the 3rd Workshop on Quantum Chaos and Localisation Phenomena Warsaw, Poland, May 25 27, 2007 Light Localization in Left-Handed Media M. Rusek,
More informationPAPER Fast Algorithm for Solving Matrix Equation in MoM Analysis of Large-Scale Array Antennas
2482 PAPER Fast Algorithm for Solving Matrix Equation in MoM Analysis of Large-Scale Array Antennas Qiang CHEN, Regular Member, Qiaowei YUAN, Nonmember, and Kunio SAWAYA, Regular Member SUMMARY A new iterative
More informationExact analytical Helmholtz bright and dark solitons
Exact analytical Helmholtz bright and dark solitons P. CHAMORRO POSADA Dpto. Teoría de la Señal y Comunicaciones e Ingeniería Telemática ETSI Telecomunicación Universidad de Valladolid, Spain G. S. McDONALD
More informationEXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC WAVELENGTH SCALE COMPONENTS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux Eds c TU Delft, The Netherlands, 2006 EXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC
More informationDifferential transfer-matrix method for solution of one-dimensional linear nonhomogeneous optical structures
S. Khorasani and K. Mehrany Vol. 20, No. 1/January 2003/J. Opt. Soc. Am. B 91 Differential transfer-matrix method for solution of one-dimensional linear nonhomogeneous optical structures Sina Khorasani
More informationFourier Decomposition Analysis af Anisotropic Inhomogeneous Dielectric Waveguide Structures
University of Pennsylvania ScholarlyCommons Departmental Papers (ESE) Department of Electrical & Systems Engineering August 2007 Fourier Decomposition Analysis af Anisotropic Inhomogeneous Dielectric Waveguide
More informationELECTROMAGNETIC SCATTERING FROM A CHIRAL- COATED NIHILITY CYLINDER
Progress In Electromagnetics Research Letters, Vol. 18, 41 5, 21 ELECTROMAGNETIC SCATTERING FROM A CHIRAL- COATED NIHILITY CYLINDER S. Ahmed and Q. A. Naqvi Department of Electronics Quaid-i-Azam University
More informationPropagation dynamics of abruptly autofocusing Airy beams with optical vortices
Propagation dynamics of abruptly autofocusing Airy beams with optical vortices Yunfeng Jiang, 1 Kaikai Huang, 1,2 and Xuanhui Lu 1, * 1 Institute of Optics, Department of Physics, Zhejiang University,
More informationLight bullets and dynamic pattern formation in nonlinear dissipative systems
Light bullets and dynamic pattern formation in nonlinear dissipative systems Philippe Grelu Laboratoire de Physique de l Université de Bourgogne, Unité Mixte de Recherche 5027 du Centre National de Recherche
More informationOptical component modelling and circuit simulation using SERENADE suite
Optical component modelling and circuit simulation using SERENADE suite Laurent Guilloton, Smail Tedjini, Tan-Phu Vuong To cite this version: Laurent Guilloton, Smail Tedjini, Tan-Phu Vuong. Optical component
More informationA RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES
Progress In Electromagnetics Research, PIER 22, 197 212, 1999 A RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES M. Akbari, M. Shahabadi, and K. Schünemann Arbeitsbereich Hochfrequenztechnik Technische
More informationPRINCIPLES OF PHYSICAL OPTICS
PRINCIPLES OF PHYSICAL OPTICS C. A. Bennett University of North Carolina At Asheville WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION CONTENTS Preface 1 The Physics of Waves 1 1.1 Introduction
More informationCanalization of Sub-wavelength Images by Electromagnetic Crystals
Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August 22-26 37 Canalization of Sub-wavelength Images by Electromagnetic Crystals P. A. Belov 1 and C. R. Simovski 2 1 Queen Mary
More informationNumerical Analysis of Low-order Modes in Thermally Diffused Expanded Core (TEC) Fibers
Proceedings of the 4th WSEAS Int. Conference on Electromagnetics, Wireless and Optical Communications, Venice, Italy, November 2-22, 26 26 Numerical Analysis of Low-order Modes in Thermally Diffused Expanded
More informationPROCEEDINGS OF SPIE. Scaling rules for the design of a narrow-band grating filter at the focus of a free-space beam
PROCEEDINGS OF SPIE SPIEDigitalLibrary.org/conference-proceedings-of-spie Scaling rules for the design of a narrow-band grating filter at the focus of a free-space beam Eanuel Bonnet, Alain Cachard, Alexandre
More informationResonant mode flopping in modulated waveguiding structures
Resonant mode flopping in modulated waveguiding structures Yaroslav V. Kartashov, Victor A. Vysloukh, and Lluis Torner ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, and Universitat
More informationarxiv: v1 [physics.optics] 30 Mar 2010
Analytical vectorial structure of non-paraxial four-petal Gaussian beams in the far field Xuewen Long a,b, Keqing Lu a, Yuhong Zhang a,b, Jianbang Guo a,b, and Kehao Li a,b a State Key Laboratory of Transient
More informationProgress In Electromagnetics Research, PIER 52, , 2005 FDTD ANALYSIS OF MICROSTRIP PATCH ANTENNA COVERED BY PLASMA SHEATH
Progress In Electromagnetics Research, PIER 52, 173 183, 25 FDTD ANALYSIS OF MICROSTRIP PATCH ANTENNA COVERED BY PLASMA SHEATH Z. H. Qian and R. S. Chen Department of Communication Engineering Nanjing
More informationSupplementary Information Supplementary Figures
Supplementary Information Supplementary Figures Supplementary Figure S1. Schematic defining the relevant physical parameters in our grating coupler model. The schematic shows a normally-incident light
More informationAngular responses of the first diffracted order in over-modulated volume diffraction gratings
Angular responses of the first diffracted order in over-modulated volume diffraction gratings C. Neipp 1, M. L. Alvarez 1, S. Gallego 2, M. Ortuño 2, J. Sheridan 3, I. Pascual 2 and A. Beléndez 1 1 Departamento
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationTime Domain Modeling of All-Optical Switch based on PT-Symmetric Bragg Grating
Time Domain Modeling of All-Optical Switch based on PT-Symmetric Bragg Grating Sendy Phang 1, Ana Vukovic 1, Hadi Susanto 2, Trevor M. Benson 1, and Phillip Sewell 1 1 School of Electrical and Electronic
More informationPhysics 3312 Lecture 7 February 6, 2019
Physics 3312 Lecture 7 February 6, 2019 LAST TIME: Reviewed thick lenses and lens systems, examples, chromatic aberration and its reduction, aberration function, spherical aberration How do we reduce spherical
More informationWAVEGUIDES FILLED WITH BILAYERS OF DOUBLE- NEGATIVE (DNG) AND DOUBLE-POSITIVE (DPS) METAMATERIALS
Progress In Electromagnetics Research B, Vol., 75 9, WAVEGUIDES FILLED WITH BILAYERS OF DOUBLE- NEGATIVE (DNG) AND DOUBLE-POSITIVE (DPS) METAMATERIALS E. Cojocaru * Department of Theoretical Physics, Horia
More informationSelf-assembled nanostructures for antireflection optical coatings
Self-assembled nanostructures for antireflection optical coatings Yang Zhao 1, Guangzhao Mao 2, and Jinsong Wang 1 1. Deaprtment of Electrical and Computer Engineering 2. Departmentof Chemical Engineering
More informationEfficiency analysis of diffractive lenses
86 J. Opt. Soc. Am. A/ Vol. 8, No. / January 00 Levy et al. Efficiency analysis of diffractive lenses Uriel Levy, David Mendlovic, and Emanuel Marom Faculty of Engineering, Tel-Aviv University, 69978 Tel-Aviv,
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationFIBER Bragg gratings are important elements in optical
IEEE JOURNAL OF QUANTUM ELECTRONICS, VOL. 40, NO. 8, AUGUST 2004 1099 New Technique to Accurately Interpolate the Complex Reflection Spectrum of Fiber Bragg Gratings Amir Rosenthal and Moshe Horowitz Abstract
More informationANALYTICAL POLE RESIDUE CALCULATION IN SPECTRAL METHOD OF MOMENTS FORMULATIONS FOR PERIODIC STRUCTURES
Progress In Electromagnetics Research M, Vol. 13, 83 94, 2010 ANALYTICAL POLE RESIDUE CALCULATION IN SPECTRAL METHOD OF MOMENTS FORMULATIONS FOR PERIODIC STRUCTURES Y. Kaganovsky Department of Physical
More informationAnalysis of Photonic Band Structure in 1-D Photonic Crystal using PWE and FDTD Method
P P IJISET - International Journal of Innovative Science, Engineering & Technology, Vol. Issue 8, August 05. Analysis of Photonic Band Structure in -D Photonic Crystal using PWE and FDTD Method Pooja ChhokerP
More informationStability and instability of solitons in inhomogeneous media
Stability and instability of solitons in inhomogeneous media Yonatan Sivan, Tel Aviv University, Israel now at Purdue University, USA G. Fibich, Tel Aviv University, Israel M. Weinstein, Columbia University,
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationOne-Dimensional Numerical Solution of the Maxwell-Minkowski Equations
Tamkang Journal of Science and Engineering, Vol. 12, No. 2, pp. 161168 (2009) 161 One-Dimensional Numerical Solution of the Maxwell-Minkowski Equations Mingtsu Ho 1 and Yao-Han Chen 2 1 Department of Electronic
More informationFAST AND ACCURATE RADAR CROSS SECTION COM- PUTATION USING CHEBYSHEV APPROXIMATION IN BOTH BROAD FREQUENCY BAND AND ANGULAR DOMAINS SIMULTANEOUSLY
Progress In Electromagnetics Research Letters, Vol. 13, 121 129, 2010 FAST AND ACCURATE RADAR CROSS SECTION COM- PUTATION USING CHEBYSHEV APPROXIMATION IN BOTH BROAD FREQUENCY BAND AND ANGULAR DOMAINS
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
More informationAPPLICATION OF BILAYER ANISOTROPIC STRUC- TURES FOR DESIGNING LOW-PASS FILTERS AND PO- LARIZERS
Progress In Electromagnetics Research M, Vol. 29, 95 108, 2013 APPLICATION OF BILAYER ANISOTROPIC STRUC- TURES FOR DESIGNING LOW-PASS FILTERS AND PO- LARIZERS Amir Raeesi *, Ali Abdolali, and Hossein Mirzaei
More informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
More informationStudy of Propagating Modes and Reflectivity in Bragg Filters with AlxGa1-xN/GaN Material Composition
Study of Propagating Modes and Reflectivity in Bragg Filters with AlxGa1-xN/GaN Material Composition Sourangsu Banerji Department of Electronics & Communication Engineering, RCC Institute of Information
More informationUnited Nations Educational, Scientific and Cultural Organization and International Atomic Energy Agency
Available at: http://publications.ictp.it IC /2010/046 United Nations Educational, Scientific and Cultural Organization and International Atomic Energy Agency THE ABDUS SALAM INTERNATIONAL CENTRE FOR THEORETICAL
More informationComparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 1014 1020 Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses Jian-Bo Chen, Yan Shen,
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationOptimum Access Waveguide Width for 1xN Multimode. Interference Couplers on Silicon Nanomembrane
Optimum Access Waveguide Width for 1xN Multimode Interference Couplers on Silicon Nanomembrane Amir Hosseini 1,*, Harish Subbaraman 2, David Kwong 1, Yang Zhang 1, and Ray T. Chen 1,* 1 Microelectronic
More informationFocal shift in vector beams
Focal shift in vector beams Pamela L. Greene The Institute of Optics, University of Rochester, Rochester, New York 1467-186 pgreene@optics.rochester.edu Dennis G. Hall The Institute of Optics and The Rochester
More informationUltra-narrow-band tunable laserline notch filter
Appl Phys B (2009) 95: 597 601 DOI 10.1007/s00340-009-3447-6 Ultra-narrow-band tunable laserline notch filter C. Moser F. Havermeyer Received: 5 December 2008 / Revised version: 2 February 2009 / Published
More informationGray spatial solitons in nonlocal nonlinear media
Gray spatial solitons in nonlocal nonlinear media Yaroslav V. Kartashov and Lluis Torner ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, and Universitat Politecnica de Catalunya, 08860
More informationA SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new
More informationA Novel Design of Photonic Crystal Lens Based on Negative Refractive Index
PIERS ONLINE, VOL. 4, NO. 2, 2008 296 A Novel Design of Photonic Crystal Lens Based on Negative Refractive Index S. Haxha 1 and F. AbdelMalek 2 1 Photonics Group, Department of Electronics, University
More informationLight as a Transverse Wave.
Waves and Superposition (Keating Chapter 21) The ray model for light (i.e. light travels in straight lines) can be used to explain a lot of phenomena (like basic object and image formation and even aberrations)
More informationEFFICIENTLY pumping a fiber over a distance of a few
IEEE JOURNAL OF QUANTUM ELECTRONICS, VOL. 40, NO. 9, SEPTEMBER 2004 1301 Slab Delivery of Incoherent Pump Light to Double-Clad Fiber Amplifiers: Numerical Simulations Patrick Kano, Dmitrii Kouznetsov,
More informationPublication II Wiley Periodicals. Reprinted by permission of John Wiley & Sons.
Publication II Ilkka Laakso and Tero Uusitupa. 2008. Alternative approach for modeling material interfaces in FDTD. Microwave and Optical Technology Letters, volume 50, number 5, pages 1211-1214. 2008
More informationAn Efficient Method to Simulate the Pulse Propagation and Switching Effects of a Fiber Bragg Grating
An Efficient Method to Simulate the Pulse Propagation and Switching Effects of a Fiber ragg Grating F. Emami, Member IAENG, A. H. Jafari, M. Hatami, and A. R. Keshavarz Abstract In this paper we investigated
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationarxiv:physics/ v1 [physics.optics] 25 Jun 1998
arxiv:physics/9806043v [physics.optics] 5 Jun 998 Nonlinear phase shift without cascaded second-order processes and third order nonlinearity V.P. Drachev, S.V. Perminov Institute of Semiconductor Physics,
More informationEvaluation of Second Order Nonlinearity in Periodically Poled
ISSN 2186-6570 Evaluation of Second Order Nonlinearity in Periodically Poled KTiOPO 4 Crystal Using Boyd and Kleinman Theory Genta Masada Quantum ICT Research Institute, Tamagawa University 6-1-1 Tamagawa-gakuen,
More informationSpatial Intensity Nonuniformities of an OMEGA Beam due to Nonlinear Beam Propagation
Spatial Intensity Nonuniformities of an OMEGA Beam due to Nonlinear Beam Propagation Several applications require that a laser beam maintain a high degree of wavefront quality after propagating a long
More informationA COMPACT POLARIZATION BEAM SPLITTER BASED ON A MULTIMODE PHOTONIC CRYSTAL WAVEGUIDE WITH AN INTERNAL PHOTONIC CRYSTAL SECTION
Progress In Electromagnetics Research, PIER 103, 393 401, 2010 A COMPACT POLARIZATION BEAM SPLITTER BASED ON A MULTIMODE PHOTONIC CRYSTAL WAVEGUIDE WITH AN INTERNAL PHOTONIC CRYSTAL SECTION Y. C. Shi Centre
More informationB. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea
Progress In Electromagnetics Research, PIER 77, 111 120, 2007 ANALYSIS OF TRANSIENT ELECTROMAGNETIC SCATTERING WITH PLANE WAVE INCIDENCE USING MOD-FDM B. H. Jung Department of Information and Communication
More informationTwo-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO2 gas
Two-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO gas M. H. Mahdieh 1, and B. Lotfi Department of Physics, Iran University of Science and Technology,
More informationWhite light generation and amplification using a soliton pulse within a nano-waveguide
Available online at www.sciencedirect.com Physics Procedia 00 (009) 000 000 53 57 www.elsevier.com/locate/procedia Frontier Research in Nanoscale Science and Technology White light generation and amplification
More informationThe concept of perfect electromagnetic conductor (PEMC) has been defined by medium conditions of the form [1, 2]
Progress In Electromagnetics Research B, Vol. 5, 169 183, 2008 REFLECTION AND TRANSMISSION OF WAVES AT THE INTERFACE OF PERFECT ELECTROMAGNETIC CONDUCTOR PEMC I. V. Lindell and A. H. Sihvola Electromagnetics
More informationStep index planar waveguide
N. Dubreuil S. Lebrun Exam without document Pocket calculator permitted Duration of the exam: 2 hours The exam takes the form of a multiple choice test. Annexes are given at the end of the text. **********************************************************************************
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationSimulations of the nonlinear Helmholtz equation: arrest of beam collapse, nonparaxial solitons and counter-propagating beams
Simulations of the nonlinear Helmholtz equation: arrest of beam collapse, nonparaxial solitons and counter-propagating beams G. Baruch, G. Fibich and Semyon Tsynkov 2 School of Mathematical Sciences, Tel
More informationTHE PROPAGATION AND CUTOFF FREQUENCIES OF THE RECTANGULAR METALLIC WAVEGUIDE PAR- TIALLY FILLED WITH METAMATERIAL MULTILAYER SLABS
Progress In Electromagnetics Research M, Vol. 9, 35 40, 2009 THE PROPAGATION AND CUTOFF FREQUENCIES OF THE RECTANGULAR METALLIC WAVEGUIDE PAR- TIALLY FILLED WITH METAMATERIAL MULTILAYER SLABS D. Zhang
More information