MATLAB Implementation of 2D FDTD with PML Boundary Conditions

We are developing a MATLAB implementation of the 2D Finite-Difference Time-Domain method, designed to simulate light propagation through various materials. The program implements Perfectly Matched Layer boundaries, a critical component that effectively simulates electromagnetic wave behavior in open space environments by minimizing unwanted reflections. Our implementation includes material modeling capabilities that describe optical properties of diverse substances such as metals, glass, and plastics through appropriate dielectric constant and conductivity parameters. The core algorithm employs Yee's grid discretization scheme with staggered electric and magnetic field components, updating fields through iterative time-stepping using Maxwell's curl equations. We continuously enhance the program's accuracy in predicting light-matter interactions, incorporating improved dispersion models and optimization techniques to advance research and applications in optical science and engineering.