Finitedifference timedomain method
Finitedifference timedomain (FDTD) or Yee's method (named after the Chinese American applied mathematician Kane S. Yee, born 1934) is a numerical analysis technique used for modeling computational electrodynamics (finding approximate solutions to the associated system of differential equations). Since it is a timedomain method, FDTD solutions can cover a wide frequency range with a single simulation run, and treat nonlinear material properties in a natural way.
The FDTD method belongs in the general class of gridbased differential numerical modeling methods (finite difference methods). The timedependent Maxwell's equations (in partial differential form) are discretized using centraldifference approximations to the space and time partial derivatives. The resulting finitedifference equations are solved in either software or hardware in a leapfrog manner: the electric field vector components in a volume of space are solved at a given instant in time; then the magnetic field vector components in the same spatial volume are solved at the next instant in time; and the process is repeated over and over again until the desired transient or steadystate electromagnetic field behavior is fully evolved.
Contents
History
Finite difference schemes for timedependent PDEs have been employed for many years in computational fluid dynamics problems,^{[1]} including the idea of using centered finite difference operators on staggered grids in space and time to achieve secondorder accuracy.^{[1]} The novelty of Kane Yee's FDTD scheme, presented in his seminal 1966 paper,^{[2]} was to apply centered finite difference operators on staggered grids in space and time for each electric and magnetic vector field component in Maxwell's curl equations. The descriptor "Finitedifference timedomain" and its corresponding "FDTD" acronym were originated by Allen Taflove in 1980.^{[3]} Since about 1990, FDTD techniques have emerged as primary means to computationally model many scientific and engineering problems dealing with electromagnetic wave interactions with material structures. Current FDTD modeling applications range from nearDC (ultralowfrequency geophysics involving the entire Earthionosphere waveguide) through microwaves (radar signature technology, antennas, wireless communications devices, digital interconnects, biomedical imaging/treatment) to visible light (photonic crystals, nanoplasmonics, solitons, and biophotonics).^{[4]} In 2006, an estimated 2,000 FDTDrelated publications appeared in the science and engineering literature (see Popularity). As of 2013, there are at least 25 commercial/proprietary FDTD software vendors; 13 freesoftware/opensourcesoftware FDTD projects; and 2 freeware/closedsource FDTD projects, some not for commercial use (see External links).
Development of FDTD and Maxwell's equations
An appreciation of the basis, technical development, and possible future of FDTD numerical techniques for Maxwell’s equations can be developed by first considering their history. The following lists some of the key publications in this area.
Partial chronology of FDTD techniques and applications for Maxwell's equations.^{[5]}  

year  event 
1928  Courant, Friedrichs, and Lewy (CFL) publish seminal paper with the discovery of conditional stability of explicit timedependent finite difference schemes, as well as the classic FD scheme for solving secondorder wave equation in 1D and 2D.^{[6]} 
1950  First appearance of von Neumann's method of stability analysis for implicit/explicit timedependent finite difference methods.^{[7]} 
1966  Yee described the FDTD numerical technique for solving Maxwell’s curl equations on grids staggered in space and time.^{[2]} 
1969  Lam reported the correct numerical CFL stability condition for Yee’s algorithm by employing von Neumann stability analysis.^{[8]} 
1975  Taflove and Brodwin reported the first sinusoidal steadystate FDTD solutions of two and threedimensional electromagnetic wave interactions with material structures;^{[9]} and the first bioelectromagnetics models.^{[10]} 
1977  Holland and Kunz & Lee applied Yee’s algorithm to EMP problems.^{[11]}^{[12]} 
1980  Taflove coined the FDTD acronym and published the first validated FDTD models of sinusoidal steadystate electromagnetic wave penetration into a threedimensional metal cavity.^{[3]} 
1981  Mur published the first numerically stable, secondorder accurate, absorbing boundary condition (ABC) for Yee’s grid.^{[13]} 
1982–83  Taflove and Umashankar developed the first FDTD electromagnetic wave scattering models computing sinusoidal steadystate nearfields, farfields, and radar crosssection for two and threedimensional structures.^{[14]}^{[15]} 
1984  Liao et al reported an improved ABC based upon spacetime extrapolation of the field adjacent to the outer grid boundary.^{[16]} 
1985  Gwarek introduced the lumped equivalent circuit formulation of FDTD.^{[17]} 
1986  Choi and Hoefer published the first FDTD simulation of waveguide structures.^{[18]} 
1987–88  Kriegsmann et al and Moore et al published the first articles on ABC theory in IEEE Transactions on Antennas and Propagation.^{[19]}^{[20]} 
1987–88, 1992  Contourpath subcell techniques were introduced by Umashankar et al to permit FDTD modeling of thin wires and wire bundles,^{[21]} by Taflove et al to model penetration through cracks in conducting screens,^{[22]} and by Jurgens et al to conformally model the surface of a smoothly curved scatterer.^{[23]} 
1988  Sullivan et al published the first 3D FDTD model of sinusoidal steadystate electromagnetic wave absorption by a complete human body.^{[24]} 
1988  FDTD modeling of microstrips was introduced by Zhang et al.^{[25]} 
1990–91  FDTD modeling of frequencydependent dielectric permittivity was introduced by Kashiwa and Fukai,^{[26]} Luebbers et al,^{[27]} and Joseph et al.^{[28]} 
1990–91  FDTD modeling of antennas was introduced by Maloney et al,^{[29]} Katz et al,^{[30]} and Tirkas and Balanis.^{[31]} 
1990  FDTD modeling of picosecond optoelectronic switches was introduced by Sano and Shibata,^{[32]} and ElGhazaly et al.^{[33]} 
1992–94  FDTD modeling of the propagation of optical pulses in nonlinear dispersive media was introduced, including the first temporal solitons in one dimension by Goorjian and Taflove;^{[34]} beam selffocusing by Ziolkowski and Judkins;^{[35]} the first temporal solitons in two dimensions by Joseph et al;^{[36]} and the first spatial solitons in two dimensions by Joseph and Taflove.^{[37]} 
1992  FDTD modeling of lumped electronic circuit elements was introduced by Sui et al.^{[38]} 
1993  Toland et al published the first FDTD models of gain devices (tunnel diodes and Gunn diodes) exciting cavities and antennas.^{[39]} 
1993  Aoyagi et al present a hybrid Yee algorithm/scalarwave equation and demonstrate equivalence of Yee scheme to finite difference scheme for electromagnetic wave equation.^{[40]} 
1994  Thomas et al introduced a Norton’s equivalent circuit for the FDTD space lattice, which permits the SPICE circuit analysis tool to implement accurate subgrid models of nonlinear electronic components or complete circuits embedded within the lattice.^{[41]} 
1994  Berenger introduced the highly effective, perfectly matched layer (PML) ABC for twodimensional FDTD grids,^{[42]} which was extended to three dimensions by Katz et al,^{[43]} and to dispersive waveguide terminations by Reuter et al.^{[44]} 
1994  Chew and Weedon introduced the coordinate stretching PML that is easily extended to three dimensions, other coordinate systems and other physical equations.^{[45]} 
1995–96  Sacks et al and Gedney introduced a physically realizable, uniaxial perfectly matched layer (UPML) ABC.^{[46]}^{[47]} 
1997  Liu introduced the pseudospectral timedomain (PSTD) method, which permits extremely coarse spatial sampling of the electromagnetic field at the Nyquist limit.^{[48]} 
1997  Ramahi introduced the complementary operators method (COM) to implement highly effective analytical ABCs.^{[49]} 
1998  Maloney and Kesler introduced several novel means to analyze periodic structures in the FDTD space lattice.^{[50]} 
1998  Nagra and York introduced a hybrid FDTDquantum mechanics model of electromagnetic wave interactions with materials having electrons transitioning between multiple energy levels.^{[51]} 
1998  Hagness et al introduced FDTD modeling of the detection of breast cancer using ultrawideband radar techniques.^{[52]} 
1999  Schneider and Wagner introduced a comprehensive analysis of FDTD grid dispersion based upon complex wavenumbers.^{[53]} 
2000–01  Zheng, Chen, and Zhang introduced the first threedimensional alternatingdirection implicit (ADI) FDTD algorithm with provable unconditional numerical stability.^{[54]}^{[55]} 
2000  Roden and Gedney introduced the advanced convolutional PML (CPML) ABC.^{[56]} 
2000  Rylander and Bondeson introduced a provably stable FDTD  finiteelement timedomain hybrid technique.^{[57]} 
2002  Hayakawa et al and Simpson and Taflove independently introduced FDTD modeling of the global Earthionosphere waveguide for extremely lowfrequency geophysical phenomena.^{[58]}^{[59]} 
2003  DeRaedt introduced the unconditionally stable, “onestep” FDTD technique.^{[60]} 
2008  Ahmed, Chua, Li and Chen introduced the threedimensional locally onedimensional (LOD)FDTD method and proved unconditional numerical stability.^{[61]} 
2008  Taniguchi, Baba, Nagaoka and Ametani introduced a Thin Wire Representation for FDTD Computations for conductive media^{[62]} 
2009  Oliveira and Sobrinho applied the FDTD method for simulating lightning strokes in a power substation^{[63]} 
2010  Chaudhury and Boeuf demonstrated the numerical procedure to couple FDTD and Plasma fluid model for studying Microwave Plasma Interaction.^{[64]} 
2012  Moxley et al developed a generalized finitedifference timedomain quantum method for the Nbody interacting Hamiltonian.^{[65]} 
2013  Moxley et al developed a generalized finitedifference timedomain scheme for solving nonlinear Schrödinger equations.^{[66]} 
2014  Moxley et al developed an implicit generalized finitedifference timedomain scheme for solving nonlinear Schrödinger equations.^{[67]} 
FDTD models and methods
When Maxwell's differential equations are examined, it can be seen that the change in the Efield in time (the time derivative) is dependent on the change in the Hfield across space (the curl). This results in the basic FDTD timestepping relation that, at any point in space, the updated value of the Efield in time is dependent on the stored value of the Efield and the numerical curl of the local distribution of the Hfield in space.^{[2]}
The Hfield is timestepped in a similar manner. At any point in space, the updated value of the Hfield in time is dependent on the stored value of the Hfield and the numerical curl of the local distribution of the Efield in space. Iterating the Efield and Hfield updates results in a marchingintime process wherein sampleddata analogs of the continuous electromagnetic waves under consideration propagate in a numerical grid stored in the computer memory.
This description holds true for 1D, 2D, and 3D FDTD techniques. When multiple dimensions are considered, calculating the numerical curl can become complicated. Kane Yee's seminal 1966 paper proposed spatially staggering the vector components of the Efield and Hfield about rectangular unit cells of a Cartesian computational grid so that each Efield vector component is located midway between a pair of Hfield vector components, and conversely.^{[2]} This scheme, now known as a Yee lattice, has proven to be very robust, and remains at the core of many current FDTD software constructs.
Furthermore, Yee proposed a leapfrog scheme for marching in time wherein the Efield and Hfield updates are staggered so that Efield updates are conducted midway during each timestep between successive Hfield updates, and conversely.^{[2]} On the plus side, this explicit timestepping scheme avoids the need to solve simultaneous equations, and furthermore yields dissipationfree numerical wave propagation. On the minus side, this scheme mandates an upper bound on the timestep to ensure numerical stability.^{[9]} As a result, certain classes of simulations can require many thousands of timesteps for completion.
Using the FDTD method
To implement an FDTD solution of Maxwell's equations, a computational domain must first be established. The computational domain is simply the physical region over which the simulation will be performed. The E and H fields are determined at every point in space within that computational domain. The material of each cell within the computational domain must be specified. Typically, the material is either freespace (air), metal, or dielectric. Any material can be used as long as the permeability, permittivity, and conductivity are specified.
The permittivity of dispersive materials in tabular form cannot be directly substituted into the FDTD scheme. Instead, it can be approximated using multiple Debye, Drude, Lorentz or critical point terms. This approximation can be obtained using open fitting programs^{[68]} and does not necessarily have physical meaning.
Once the computational domain and the grid materials are established, a source is specified. The source can be current on a wire, applied electric field or impinging plane wave. In the last case FDTD can be used to simulate light scattering from arbitrary shaped objects, planar periodic structures at various incident angles,^{[69]}^{[70]} and photonic band structure of infinite periodic structures.^{[71]}^{[72]}
Since the E and H fields are determined directly, the output of the simulation is usually the E or H field at a point or a series of points within the computational domain. The simulation evolves the E and H fields forward in time.
Processing may be done on the E and H fields returned by the simulation. Data processing may also occur while the simulation is ongoing.
While the FDTD technique computes electromagnetic fields within a compact spatial region, scattered and/or radiated far fields can be obtained via neartofarfield transformations.^{[14]}
Strengths of FDTD modeling
Every modeling technique has strengths and weaknesses, and the FDTD method is no different.
 FDTD is a versatile modeling technique used to solve Maxwell's equations. It is intuitive, so users can easily understand how to use it and know what to expect from a given model.
 FDTD is a timedomain technique, and when a broadband pulse (such as a Gaussian pulse) is used as the source, then the response of the system over a wide range of frequencies can be obtained with a single simulation. This is useful in applications where resonant frequencies are not exactly known, or anytime that a broadband result is desired.
 Since FDTD calculates the E and H fields everywhere in the computational domain as they evolve in time, it lends itself to providing animated displays of the electromagnetic field movement through the model. This type of display is useful in understanding what is going on in the model, and to help ensure that the model is working correctly.
 The FDTD technique allows the user to specify the material at all points within the computational domain. A wide variety of linear and nonlinear dielectric and magnetic materials can be naturally and easily modeled.
 FDTD allows the effects of apertures to be determined directly. Shielding effects can be found, and the fields both inside and outside a structure can be found directly or indirectly.
 FDTD uses the E and H fields directly. Since most EMI/EMC modeling applications are interested in the E and H fields, it is convenient that no conversions must be made after the simulation has run to get these values.
Weaknesses of FDTD modeling
 Since FDTD requires that the entire computational domain be gridded, and the grid spatial discretization must be sufficiently fine to resolve both the smallest electromagnetic wavelength and the smallest geometrical feature in the model, very large computational domains can be developed, which results in very long solution times. Models with long, thin features, (like wires) are difficult to model in FDTD because of the excessively large computational domain required. Methods such as Eigenmode Expansion can offer a more efficient alternative as they do not require a fine grid along the zdirection.^{[73]}
 There is no way to determine unique values for permittivity and permeability at a material interface.
 Space and time steps must satisfy the CFL condition, or the leapfrog integration used to solve the partial differential equation is probable to become unstable.
 FDTD finds the E/H fields directly everywhere in the computational domain. If the field values at some distance are desired, it is likely that this distance will force the computational domain to be excessively large. Farfield extensions are available for FDTD, but require some amount of postprocessing.^{[4]}
 Since FDTD simulations calculate the E and H fields at all points within the computational domain, the computational domain must be finite to permit its residence in the computer memory. In many cases this is achieved by inserting artificial boundaries into the simulation space. Care must be taken to minimize errors introduced by such boundaries. There are a number of available highly effective absorbing boundary conditions (ABCs) to simulate an infinite unbounded computational domain.^{[4]} Most modern FDTD implementations instead use a special absorbing "material", called a perfectly matched layer (PML) to implement absorbing boundaries.^{[42]}^{[46]}
 Because FDTD is solved by propagating the fields forward in the time domain, the electromagnetic time response of the medium must be modeled explicitly. For an arbitrary response, this involves a computationally expensive time convolution, although in most cases the time response of the medium (or Dispersion (optics)) can be adequately and simply modeled using either the recursive convolution (RC) technique, the auxiliary differential equation (ADE) technique, or the Ztransform technique. An alternative way of solving Maxwell's equations that can treat arbitrary dispersion easily is the Pseudospectral SpatialDomain method (PSSD), which instead propagates the fields forward in space.
Grid truncation techniques
The most commonly used grid truncation techniques for openregion FDTD modeling problems are the Mur absorbing boundary condition (ABC),^{[13]} the Liao ABC,^{[16]} and various perfectly matched layer (PML) formulations.^{[4]}^{[42]}^{[46]} The Mur and Liao techniques are simpler than PML. However, PML (which is technically an absorbing region rather than a boundary condition per se) can provide ordersofmagnitude lower reflections. The PML concept was introduced by J.P. Berenger in a seminal 1994 paper in the Journal of Computational Physics.^{[42]} Since 1994, Berenger's original splitfield implementation has been modified and extended to the uniaxial PML (UPML), the convolutional PML (CPML), and the higherorder PML. The latter two PML formulations have increased ability to absorb evanescent waves, and therefore can in principle be placed closer to a simulated scattering or radiating structure than Berenger's original formulation.
To reduce undesired numerical reflection from the PML additional back absorbing layers technique can be used.^{[74]}
Popularity
This section possibly contains original research. (August 2013)

Notwithstanding both the general increase in academic publication throughput during the same period and the overall expansion of interest in all Computational electromagnetics (CEM) techniques, there are seven primary reasons for the tremendous expansion of interest in FDTD computational solution approaches for Maxwell’s equations:
 FDTD uses no linear algebra. Being a fully explicit computation, FDTD avoids the difficulties with linear algebra that limit the size of frequencydomain integralequation and finiteelement electromagnetics models to generally fewer than 10^{9} electromagnetic field unknowns.^{[4]} FDTD models with as many as 10^{9} field unknowns have been run; there is no intrinsic upper bound to this number.^{[4]}
 FDTD is accurate and robust. The sources of error in FDTD calculations are well understood, and can be bounded to permit accurate models for a very large variety of electromagnetic wave interaction problems.^{[4]}
 FDTD treats impulsive behavior naturally. Being a timedomain technique, FDTD directly calculates the impulse response of an electromagnetic system. Therefore, a single FDTD simulation can provide either ultrawideband temporal waveforms or the sinusoidal steadystate response at any frequency within the excitation spectrum.^{[4]}
 FDTD treats nonlinear behavior naturally. Being a timedomain technique, FDTD directly calculates the nonlinear response of an electromagnetic system. This allows natural hybriding of FDTD with sets of auxiliary differential equations that describe nonlinearities from either the classical or semiclassical standpoint.^{[4]} One research frontier is the development of hybrid algorithms which join FDTD classical electrodynamics models with phenomena arising from quantum electrodynamics, especially vacuum fluctuations, such as the Casimir effect.^{[4]}^{[75]}
 FDTD is a systematic approach. With FDTD, specifying a new structure to be modeled is reduced to a problem of mesh generation rather than the potentially complex reformulation of an integral equation. For example, FDTD requires no calculation of structuredependent Green functions.^{[4]}
 Parallelprocessing computer architectures have come to dominate supercomputing. FDTD scales with high efficiency on parallelprocessing CPUbased computers, and extremely well on recently developed GPUbased accelerator technology.^{[4]}
 Computer visualization capabilities are increasing rapidly. While this trend positively influences all numerical techniques, it is of particular advantage to FDTD methods, which generate timemarched arrays of field quantities suitable for use in color videos to illustrate the field dynamics.^{[4]}
Taflove has argued that these factors combine to suggest that FDTD will remain one of the dominant computational electrodynamics techniques (as well as potentially other multiphysics problems).^{[4]}
Implementations
There are hundreds of simulation tools that implement FDTD algorithms, many optimized to run on parallelprocessing clusters.
Frederick Moxley suggests further applications with computational quantum mechanics and simulations.^{[76]}
See also
 Computational electromagnetics
 Eigenmode expansion
 Beam propagation method
 Finitedifference frequencydomain
 Scattering Matrix Method
 Discrete dipole approximation
References
 ↑ ^{1.0} ^{1.1} J. von Neumann and RD Richtmyer (March 1950). "A method for the numerical calculation of hydrodynamic shocks". Journal of Applied Physics. 21: 232–237. Bibcode:1950JAP....21..232V. doi:10.1063/1.1699639.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{2.0} ^{2.1} ^{2.2} ^{2.3} ^{2.4} ^{2.5} Kane Yee (1966). "Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media". IEEE Transactions on Antennas and Propagation. 14 (3): 302–307. Bibcode:1966ITAP...14..302Y. doi:10.1109/TAP.1966.1138693.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{3.0} ^{3.1} A. Taflove (1980). "Application of the finitedifference timedomain method to sinusoidal steady state electromagnetic penetration problems" (PDF). IEEE Transactions on Electromagnetic Compatibility. 22 (3): 191–202. doi:10.1109/TEMC.1980.303879.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{4.00} ^{4.01} ^{4.02} ^{4.03} ^{4.04} ^{4.05} ^{4.06} ^{4.07} ^{4.08} ^{4.09} ^{4.10} ^{4.11} ^{4.12} ^{4.13} Allen Taflove and Susan C. Hagness (2005). Computational Electrodynamics: The FiniteDifference TimeDomain Method, 3rd ed. Artech House Publishers. ISBN 1580538320.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Adapted with permission from Taflove and Hagness (2005).
 ↑ Richard Courant, Kurt Otto Friedrichs, and Hans Lewy (1928). "Über die partiellen Differenzengleichungen der mathematischen Physik". Mathematische Annalen (in German). 100 (1): 32–74. Bibcode:1928MatAn.100...32C. doi:10.1007/BF01448839. JFM 54.0486.01. MR 1512478.CS1 maint: multiple names: authors list (link) CS1 maint: unrecognized language (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ G. G. O’Brien, M. A Hyman, and S. Kaplan (1950). "A study of the numerical solution of partial differential equations". Journal of Mathematical Physics. 29 (1): 223–251. MR 0040805.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ DongHoa Lam (1969). "Finite Difference Methods for Electromagnetic Scattering Problems". Mississippi State University, Interaction Notes. 44.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{9.0} ^{9.1} A. Taflove and M. E. Brodwin (1975). "Numerical solution of steadystate electromagnetic scattering problems using the timedependent Maxwell's equations" (PDF). IEEE Transactions on Microwave Theory and Techniques. 23 (8): 623–630. Bibcode:1975ITMTT..23..623T. doi:10.1109/TMTT.1975.1128640.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Taflove and M. E. Brodwin (1975). "Computation of the electromagnetic fields and induced temperatures within a model of the microwaveirradiated human eye" (PDF). IEEE Transactions on Microwave Theory and Techniques. 23 (11): 888–896. Bibcode:1975ITMTT..23..888T. doi:10.1109/TMTT.1975.1128708.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. Holland (1977). "Threde: A freefield EMP coupling and scattering code". IEEE Transactions on Nuclear Science. 24 (6): 2416–2421. Bibcode:1977ITNS...24.2416H. doi:10.1109/TNS.1977.4329229.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ K. S. Kunz and K. M. Lee (1978). "A threedimensional finitedifference solution of the external response of an aircraft to a complex transient EM environment". IEEE Transactions on Electromagnetic Compatibility. 20 (2): 333–341. doi:10.1109/TEMC.1978.303727.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{13.0} ^{13.1} G. Mur (1981). "Absorbing boundary conditions for the finitedifference approximation of the timedomain electromagnetic field equations". IEEE Transactions on Electromagnetic Compatibility. 23 (4): 377–382. doi:10.1109/TEMC.1981.303970.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{14.0} ^{14.1} K. R. Umashankar and A. Taflove (1982). "A novel method to analyze electromagnetic scattering of complex objects" (PDF). IEEE Transactions on Electromagnetic Compatibility. 24 (4): 397–405. doi:10.1109/TEMC.1982.304054.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Taflove and K. R. Umashankar (1983). "Radar cross section of general threedimensional scatterers" (PDF). IEEE Transactions on Electromagnetic Compatibility. 25 (4): 433–440. doi:10.1109/TEMC.1983.304133.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{16.0} ^{16.1} Z. P. Liao, H. L. Wong, B. P. Yang, and Y. F. Yuan (1984). "A transmitting boundary for transient wave analysis". Scientia Sinica a. 27: 1063–1076.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ W. Gwarek (1985). "Analysis of an arbitrarily shaped planar circuit — A timedomain approach". IEEE Transactions on Microwave Theory and Techniques. 33 (10): 1067–1072. Bibcode:1985ITMTT..33.1067G. doi:10.1109/TMTT.1985.1133170.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ D. H. Choi and W. J. Hoefer (1986). "The finitedifference timedomain method and its application to eigenvalue problems". IEEE Transactions on Microwave Theory and Techniques. 34 (12): 1464–1470. Bibcode:1986ITMTT..34.1464C. doi:10.1109/TMTT.1986.1133564.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ G. A. Kriegsmann, A. Taflove, and K. R. Umashankar (1987). "A new formulation of electromagnetic wave scattering using an onsurface radiation boundary condition approach" (PDF). IEEE Transactions on Antennas and Propagation. 35 (2): 153–161. Bibcode:1987ITAP...35..153K. doi:10.1109/TAP.1987.1144062.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ T. G. Moore, J. G. Blaschak, A. Taflove, and G. A. Kriegsmann (1988). "Theory and application of radiation boundary operators" (PDF). IEEE Transactions on Antennas and Propagation. 36 (12): 1797–1812. Bibcode:1988ITAP...36.1797M. doi:10.1109/8.14402.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ K. R. Umashankar, A. Taflove, and B. Beker (1987). "Calculation and experimental validation of induced currents on coupled wires in an arbitrary shaped cavity" (PDF). Antennas and Propagation, IEEE Transactions on. 35 (11): 1248–1257. Bibcode:1987ITAP...35.1248U. doi:10.1109/TAP.1987.1144000.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Taflove, K. R. Umashankar, B. Beker, F. A. Harfoush, and K. S. Yee (1988). "Detailed FDTD analysis of electromagnetic fields penetrating narrow slots and lapped joints in thick conducting screens" (PDF). Antennas and Propagation, IEEE Transactions on. 36 (2): 247–257. Bibcode:1988ITAP...36..247T. doi:10.1109/8.1102.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ T. G. Jurgens, A. Taflove, K. R. Umashankar, and T. G. Moore (1992). "Finitedifference timedomain modeling of curved surfaces" (PDF). Antennas and Propagation, IEEE Transactions on. 40 (4): 357–366. Bibcode:1992ITAP...40..357J. doi:10.1109/8.138836.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ D. M. Sullivan, O. P. Gandhi, and A. Taflove (1988). "Use of the finitedifference timedomain method in calculating EM absorption in man models" (PDF). IEEE Transactions on Biomedical Engineering. 35 (3): 179–186. doi:10.1109/10.1360. PMID 3350546.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ X. Zhang, J. Fang, K. K. Mei, and Y. Liu (1988). "Calculation of the dispersive characteristics of microstrips by the timedomain finitedifference method". IEEE Transactions on Microwave Theory and Techniques. 36 (2): 263–267. Bibcode:1988ITMTT..36..263Z. doi:10.1109/22.3514.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ T. Kashiwa and I. Fukai (1990). "A treatment by FDTD method of dispersive characteristics associated with electronic polarization". Microwave and Optical Technology Letters. 3 (6): 203–205. doi:10.1002/mop.4650030606.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. Luebbers, F. Hunsberger, K. Kunz, R. Standler, and M. Schneider (1990). "A frequencydependent finitedifference timedomain formulation for dispersive materials". Electromagnetic Compatibility, IEEE Transactions on. 32 (3): 222–227. doi:10.1109/15.57116.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. M. Joseph, S. C. Hagness, and A. Taflove (1991). "Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses" (PDF). Optics Letters. 16 (18): 1412–4. Bibcode:1991OptL...16.1412J. doi:10.1364/OL.16.001412. PMID 19776986.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ J. G. Maloney, G. S. Smith, and W. R. Scott, Jr. (1990). "Accurate computation of the radiation from simple antennas using the finitedifference timedomain method". Antennas and Propagation, IEEE Transactions on. 38 (7): 1059–1068. Bibcode:1990ITAP...38.1059M. doi:10.1109/8.55618.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ D. S. Katz, A. Taflove, M. J. PiketMay, and K. R. Umashankar (1991). "FDTD analysis of electromagnetic wave radiation from systems containing horn antennas" (PDF). Antennas and Propagation, IEEE Transactions on. 39 (8): 1203–1212. Bibcode:1991ITAP...39.1203K. doi:10.1109/8.97356.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ P. A. Tirkas and C. A. Balanis (1991). "Finitedifference timedomain technique for radiation by horn antennas". Antennas and Propagation Society International Symposium Digest, IEEE. 3: 1750–1753. doi:10.1109/APS.1991.175196. ISBN 0780301447.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ E. Sano and T. Shibata (1990). "Fullwave analysis of picosecond photoconductive switches". Quantum Electronics, IEEE Journal of. 26 (2): 372–377. Bibcode:1990IJQE...26..372S. doi:10.1109/3.44970.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ S. M. ElGhazaly, R. P. Joshi, and R. O. Grondin (1990). "Electromagnetic and transport considerations in subpicosecond photoconductive switch modeling". Microwave Theory and Techniques, IEEE Transactions on. 38 (5): 629–637. Bibcode:1990ITMTT..38..629E. doi:10.1109/22.54932.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ P. M. Goorjian and A. Taflove (1992). "Direct time integration of Maxwell's equations in nonlinear dispersive media for propagation and scattering of femtosecond electromagnetic solitons" (PDF). Optics Letters. 17 (3): 180–182. Bibcode:1992OptL...17..180G. doi:10.1364/OL.17.000180.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. W. Ziolkowski and J. B. Judkins (1993). "Fullwave vector Maxwell's equations modeling of selffocusing of ultrashort optical pulses in a nonlinear Kerr medium exhibiting a finite response time". Journal of the Optical Society of America B. 10 (2): 186–198. Bibcode:1993JOSAB..10..186Z. doi:10.1364/JOSAB.10.000186.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. M. Joseph, P. M. Goorjian, and A. Taflove (1993). "Direct time integration of Maxwell's equations in 2D dielectric waveguides for propagation and scattering of femtosecond electromagnetic solitons" (PDF). Optics Letters. 18 (7): 491–3. Bibcode:1993OptL...18..491J. doi:10.1364/OL.18.000491. PMID 19802177.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. M. Joseph and A. Taflove (1994). "Spatial soliton deflection mechanism indicated by FDTD Maxwell's equations modeling" (PDF). IEEE Photonics Technology Letters. 2: 1251–1254. Bibcode:1994IPTL....6.1251J. doi:10.1109/68.329654.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ W. Sui, D. A. Christensen, and C. H. Durney (1992). "Extending the twodimensional FDTD method to hybrid electromagnetic systems with active and passive lumped elements". IEEE Transactions on Microwave Theory and Techniques. 40 (4): 724–730. Bibcode:1992ITMTT..40..724S. doi:10.1109/22.127522.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ B. Toland, B. Houshmand, and T. Itoh (1993). "Modeling of nonlinear active regions with the FDTD method". IEEE Microwave and Guided Wave Letters. 3 (9): 333–335. doi:10.1109/75.244870.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Aoyagi, P.H. and Lee, J.F. and Mittra, R. (1993). "A hybrid Yee algorithm/scalarwave equation approach". IEEE Transactions on Microwave Theory and Techniques. IEEE. 41 (9): 1593–1600. Bibcode:1993ITMTT..41.1593A. doi:10.1109/22.245683.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ V. A. Thomas, M. E. Jones, M. J. PiketMay, A. Taflove, and E. Harrigan (1994). "The use of SPICE lumped circuits as subgrid models for FDTD highspeed electronic circuit design" (PDF). IEEE Microwave and Guided Wave Letters. 4: 141–143. doi:10.1109/75.289516.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{42.0} ^{42.1} ^{42.2} ^{42.3} J. Berenger (1994). "A perfectly matched layer for the absorption of electromagnetic waves" (PDF). Journal of Computational Physics. 114 (2): 185–200. Bibcode:1994JCoPh.114..185B. doi:10.1006/jcph.1994.1159.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ D. S. Katz, E. T. Thiele, and A. Taflove (1994). "Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FDTD meshes" (PDF). IEEE Microwave and Guided Wave Letters. 4 (8): 268–270. doi:10.1109/75.311494.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ C. E. Reuter, R. M. Joseph, E. T. Thiele, D. S. Katz, and A. Taflove (1994). "Ultrawideband absorbing boundary condition for termination of waveguiding structures in FDTD simulations" (PDF). IEEE Microwave and Guided Wave Letters. 4 (10): 344–346. doi:10.1109/75.324711.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ W.C. Chew and W.H. Weedon (1994). "A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates". Microwave Optical Technology Letters. 7 (13): 599–604. doi:10.1002/mop.4650071304.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ ^{46.0} ^{46.1} ^{46.2} S. D. Gedney (1996). "An anisotropic perfectly matched layer absorbing media for the truncation of FDTD lattices". IEEE Transactions on Antennas and Propagation. 44 (12): 1630–1639. Bibcode:1996ITAP...44.1630G. doi:10.1109/8.546249.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Z. S. Sacks, D. M. Kingsland, R. Lee, and J. F. Lee (1995). "A perfectly matched anisotropic absorber for use as an absorbing boundary condition". IEEE Transactions on Antennas and Propagation. 43 (12): 1460–1463. Bibcode:1995ITAP...43.1460S. doi:10.1109/8.477075.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Q. H. Liu (1997). "The pseudospectral timedomain (PSTD) method: A new algorithm for solutions of Maxwell's equations". IEEE Antennas and Propagation Society International Symposium Digest. 1: 122–125. doi:10.1109/APS.1997.630102. ISBN 0780341783.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ O. M. Ramahi (1997). "The complementary operators method in FDTD simulations". IEEE Antennas and Propagation Magazine. 39 (6): 33–45. Bibcode:1997IAPM...39...33R. doi:10.1109/74.646801.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ J. G. Maloney and M. P. Kesler (1998). "Analysis of Periodic Structures". Chap. 6 in Advances in Computational Electrodynamics: the FiniteDifference TimeDomain Method, A. Taflove, ed., Artech House, publishers.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. S. Nagra and R. A. York (1998). "FDTD analysis of wave propagation in nonlinear absorbing and gain media". IEEE Transactions on Antennas and Propagation. 46 (3): 334–340. Bibcode:1998ITAP...46..334N. doi:10.1109/8.662652.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ S. C. Hagness, A. Taflove, and J. E. Bridges (1998). "Twodimensional FDTD analysis of a pulsed microwave confocal system for breast cancer detection: Fixedfocus and antennaarray sensors" (PDF). IEEE Transactions on Biomedical Engineering. 45 (12): 1470–1479. doi:10.1109/10.730440. PMID 9835195.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ J. B. Schneider, and C. L. Wagner (1999). "FDTD dispersion revisited: Fasterthanlight propagation". IEEE Microwave and Guided Wave Letters. 9 (2): 54–56. doi:10.1109/75.755044.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ F. Zhen, Z. Chen, and J. Zhang (2000). "Toward the development of a threedimensional unconditionally stable finitedifference timedomain method". Microwave Theory and Techniques, IEEE Transactions on. 48 (9): 1550–1558. Bibcode:2000ITMTT..48.1550Z. doi:10.1109/22.869007.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ F. Zheng and Z. Chen (2001). "Numerical dispersion analysis of the unconditionally stable 3D ADIFDTD method". Microwave Theory and Techniques, IEEE Transactions on. 49 (5): 1006–1009. Bibcode:2001ITMTT..49.1006Z. doi:10.1109/22.920165.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ J. A. Roden and S. D. Gedney (2000). "Convolution PML (CPML): An efficient FDTD implementation of the CFSPML for arbitrary media". Microwave and Optical Technology Letters. 27 (5): 334–339. doi:10.1002/10982760(20001205)27:5<334::AIDMOP14>3.0.CO;2A.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ T. Rylander and A. Bondeson (2000). "Stable FDTDFEM hybrid method for Maxwell's equations". Computer Physics Communications. 125: 75–82. Bibcode:2000CoPhC.125...75R. doi:10.1016/S00104655(99)004634.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ M. Hayakawa and T. Otsuyama (2002). "FDTD analysis of ELF wave propagation in inhomogeneous subionospheric waveguide models". ACES Journal. 17: 239–244.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ J. J. Simpson and A. Taflove (2002). "Twodimensional FDTD model of antipodal ELF propagation and Schumann resonance of the Earth" (PDF). IEEE Antennas and Wireless Propagation Letters. 1 (2): 53–56. Bibcode:2002IAWPL...1...53S. doi:10.1109/LAWP.2002.805123.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ H. De Raedt, K. Michielsen, J. S. Kole, and M. T. Figge (2003). "Solving the Maxwell equations by the Chebyshev method: A onestep finite difference timedomain algorithm". IEEE Transactions on Antennas and Propagation. 51 (11): 3155–3160. arXiv:physics/0208060. Bibcode:2003ITAP...51.3155D. doi:10.1109/TAP.2003.818809.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ I. Ahmed, E. K. Chua, E. P. Li and Z. Chen (2008). "Development of the threedimensional unconditionally stable LODFDTD method". IEEE Transactions on Antennas and Propagation. 56 (11): 3596–3600. Bibcode:2008ITAP...56.3596A. doi:10.1109/TAP.2008.2005544.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Y.Taniguchi, Y. Baba, N. Nagaoka and A. Ametani (2008). "Taniguchi, Y. and Baba, Y. and Nagaoka, N. and Ametani, A.". IEEE Transactions on Antennas and Propagation. IEEE. 56 (10): 3248–3252. Bibcode:2008ITAP...56.3248T. doi:10.1109/TAP.2008.929447.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ R. M. S. de Oliveira and C. L. S. S. Sobrinho (2009). "Computational Environment for Simulating Lightning Strokes in a Power Substation by FiniteDifference TimeDomain Method". IEEE Transactions on Electromagnetic Compatibility. IEEE. 51 (4): 995–1000. doi:10.1109/TEMC.2009.2028879.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ B. Chaudhury and J. P. Boeuf (2010). "Computational Studies of Filamentary Pattern Formation in a High Power Microwave Breakdown Generated Air Plasma". IEEE Transactions on Plasma Science. 38 (9): 2281–2288. Bibcode:2010ITPS...38.2281C. doi:10.1109/TPS.2010.2055893.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ F. I. Moxley III, T. Byrnes, F. Fujiwara, and W. Dai (2012). "A generalized finitedifference timedomain quantum method for the Nbody interacting Hamiltonian". Computer Physics Communications. 183 (11): 2434–2440. Bibcode:2012CoPhC.183.2434M. doi:10.1016/j.cpc.2012.06.012.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ F. I. Moxley III, D. T. Chuss, and W. Dai (2013). "A generalized finitedifference timedomain scheme for solving nonlinear Schrödinger equations". Computer Physics Communications. 184 (8): 1834–1841. Bibcode:2013CoPhC.184.1834M. doi:10.1016/j.cpc.2013.03.006.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Frederick Moxley; et al. (2014). Contemporary Mathematics: Mathematics of Continuous and Discrete Dynamical Systems. American Mathematical Society. ISBN 9780821898628.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ "Fitting of dielectric function".<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ I. Valuev, A. Deinega, and S. Belousov (2008). "Iterative technique for analysis of periodic structures at oblique incidence in the finitedifference timedomain method". Opt. Lett. 33: 1491. Bibcode:2008OptL...33.1491V. doi:10.1364/ol.33.001491.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Aminian and Y. RahmatSamii (2006). "Spectral FDTD: a novel technique for the analysis of oblique incident plane wave on periodic structures". IEEE Trans. Antennas and Propagation. 54: 1818–1825. Bibcode:2006ITAP...54.1818A. doi:10.1109/tap.2006.875484.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Deinega, S. Belousov, and I. Valuev (2009). "Hybrid transfermatrix FDTD method for layered periodic structures". Opt. Lett. 34: 860. Bibcode:2009OptL...34..860D. doi:10.1364/ol.34.000860.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ Y. Hao and R. Mittra (2009). FDTD Modeling of Metamaterials: Theory and Applications. Artech House Publishers.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ D. Gallagher (2008). "Photonics CAD Matures" (PDF). LEOS Newsletter.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ A. Deinega and I. Valuev (2011). "Longtime behavior of PML absorbing boundaries for layered periodic structures". Comp. Phys. Comm. 182: 149–151. Bibcode:2011CoPhC.182..149D. doi:10.1016/j.cpc.2010.06.006.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 ↑ S. G. Johnson, "Numerical methods for computing Casimir interactions," in Casimir Physics (D. Dalvit, P. Milonni, D. Roberts, and F. da Rosa, eds.), vol. 834 of Lecture Notes in Physics, ch. 6, pp. 175–218, Berlin: Springer, June 2011.
 ↑ Hartmut Ruhl, Nils Moscḧuring, and Nina Elkina (2012). "Computational Physics Course 17104 Lecture 9" (PDF). Cite journal requires
journal=
(help)CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
Further reading
The following article in Nature Milestones: Photons illustrates the historical significance of the FDTD method as related to Maxwell's equations:
 David Pile (May 2010). "Milestone 2 (1861) Maxwell's equations". Nature Milestones: Photons. Nature. doi:10.1038/nmat2639. Retrieved 17 June 2010.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
Allen Taflove's interview, "Numerical Solution," in the January 2015 focus issue of Nature Photonics honoring the 150th anniversary of the publication of Maxwell's equations. This interview touches on how the development of FDTD ties into the century and onehalf history of Maxwell's theory of electrodynamics:
The following universitylevel textbooks provide a good general introduction to the FDTD method:
 Karl S. Kunz and Raymond J. Luebbers (1993). The Finite Difference Time Domain Method for Electromagnetics. CRC Press. ISBN 0849386578.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 Allen Taflove and Susan C. Hagness (2005). Computational Electrodynamics: The FiniteDifference TimeDomain Method, 3rd ed. Artech House Publishers. ISBN 1580538320.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 Wenhua Yu, Raj Mittra, Tao Su, Yongjun Liu, and Xiaoling Yang (2006). Parallel FiniteDifference TimeDomain Method. Artech House Publishers. ISBN 1596930853.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
 John B. Schneider (2010). Understanding the FDTD Method. available online.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
External links
Free software/Opensource software FDTD projects:
 FDTD++: advanced, fully featured FDTD software, with included C++ source code, along with sophisticated material models and predefined fits as well as discussion/support forums and email support
 openEMS (Fully 3D Cartesian & Cylindrical graded mesh ECFDTD Solver, written in C++, using a Matlab/OctaveInterface)
 pFDTD (3D C++ FDTD codes developed by SeHeon Kim)
 JFDTD (2D/3D C++ FDTD codes developed for nanophotonics by Jeffrey M. McMahon)
 WOLFSIM (NCSU) (2D)
 Meep (MIT, 2D/3D/cylindrical parallel FDTD)
 (Geo) Radar FDTD
 bigboy (unmaintained, no release files. must get source from cvs)
 toyFDTD
 Parallel (MPI&OpenMP) FDTD codes in C++ (developed by Zs. Szabó)
 FDTD code in Fortran 90
 FDTD code in C for 2D EM Wave simulation
 Angora (3D parallel FDTD software package, maintained by Ilker R. Capoglu)
 GSvit (3D FDTD solver with graphics card computing support, written in C, graphical user interface XSvit available)
 gprMax (Open Source (GPLv3), 3D/2D FDTD modelling code in Python/Cython developed for GPR but can be used for general EM modelling.)
Freeware/Closed source FDTD projects (some not for commercial use):
 EMTL (Electromagnetic Template Library) (Free С++ library for electromagnetic simulations. The current version implements mainly the FDTD).
 CS1 maint: multiple names: authors list
 CS1 maint: unrecognized language
 CS1 errors: missing periodical
 Articles that may contain original research from August 2013
 Numerical software
 Simulation software
 Electromagnetic radiation
 Numerical differential equations
 Computational science
 Computational physics
 Electromagnetism
 Electrodynamics
 Scattering, absorption and radiative transfer (optics)