Abstract

The one-step leapfrog alternating-direction-implicit finite-difference time-domain (ADI-FDTD) method is reformulated for simulating general electrically dispersive media. It models material dispersive properties with equivalent polarization currents. These currents are then solved with the auxiliary differential equation (ADE) and then incorporated into the one-step leapfrog ADI-FDTD method. The final equations are presented in the form similar to that of the conventional FDTD method but with second-order perturbation. The adapted method is then applied to characterize (a) electromagnetic wave propagation in a rectangular waveguide loaded with a magnetized plasma slab, (b) transmission coefficient of a plane wave normally incident on a monolayer graphene sheet biased by a magnetostatic field, and (c) surface plasmon polaritons (SPPs) propagation along a monolayer graphene sheet biased by an electrostatic field. The numerical results verify the stability, accuracy and computational efficiency of the proposed one-step leapfrog ADI-FDTD algorithm in comparison with analytical results and the results obtained with the other methods.

© 2013 Optical Society of America

1. Introduction

Due to its simplicity in directly discretizing Maxwell’s curl equations in time domain and its ability for obtaining a wideband response with a single run of simulation, the Finite-Difference Time-Domain method (FDTD) has been widely used for modeling various isotropic and anisotropic dispersive media and structures for RF to Optical applications [110]. In general, the FDTD schemes for modeling dispersive media can be divided into several categories: auxiliary differential equation (ADE) method [11], Z-transform technique [12], and discrete convolution of the dispersion relation (DC-DR) method [13]. Among them, the ADE method is simple and applicable for handling linear and nonlinear dispersive media, and its capability can be enhanced greatly by using the complex- conjugate pole-residue pair (CC-PR), critical points (CP), and rational-fraction dispersion (RFM) techniques developed recently [1416]. Different from the ADE method, the DC-DR method proposed in [17,18] can model all dispersive media with a single general formulation. It is based on the discretization of the dispersion relation with the convolutions that include current density convolution (JEC), recursive convolution (RC), and piecewise linear recursive convolution (PLRC) [19,20]. However, their numerical implementations are more complex than that of the ADE with little accuracy improvement.

Many FDTD methods were developed in the past few years, such as EJ collocated FDTD [21], Runge–Kutta exponential time differencing formulation (RKETD) FDTD [22], matrix exponential (ME) FDTD [23], and exponential time differencing (ETD) FDTD [24], etc. Unfortunately, they cannot be easily used to model the general dispersive media.

On the other hand, in order to overcome the Courant–Friedrich–Levy (CFL) of the conventional FDTD method, one-step leapfrog ADI-FDTD method has been developed [25]. It is originated from the conventional ADI-FDTD but with no mid-time computation needed [26]. Therefore, it is similar to the conventional FDTD and has higher computational efficiency than the other unconditionally stable FDTD methods [27]. Recently, both convolutional perfectly matched layer (CPML) and Mur absorbing boundary condition for the leapfrog ADI-FDTD have been developed to simulate open structure problems [28,29]. While in [30], the leapfrog ADI-FDTD method for periodic structures was presented. More recently, the divergence properties of the leapfrog ADI-FDTD was studied in [31], and it has been further extended to simulating the lossy media [32]. In addition, the leapfrog ADI-FDTD has also been adapted to modeling magnetized plasma [33].

In this paper, we further extend the leapfrog ADI-FDTD method to simulating the general dispersive media based on the newly developed ADE method presented in [17,18]. Its successful applications are demonstrated by simulating the magnetized plasma, monolayer graphene sheet biased by a magnetostatic field, and surface plasmon polaritons (SPPs) propagating along the graphene sheet biased by an electrostatic field. Evidently, it can also be employed to simulate various graphene-based nonreciprocal RF, THz nanoelectonic and nanophotonic devices and components.

2. The proposed leapfrog ADI-FDTD method for general dispersive media

For a generalized electrically anisotropic media, the Ampere’s law in frequency domain can be written as

jωε0εreE(ω)=×H(ω)
Where E=[Ex,Ey,Ez]T, H=[Hx,Hy,Hz]T, ε0is the permittivity in free space, and εre is the relative permittivity tensor, i.e.
εre(ω)=[εxx(ω)εxy(ω)εxz(ω)εyx(ω)εyy(ω)εyz(ω)εzx(ω)εzy(ω)εzz(ω)]
Therefore, the x-component of ε0εre(ω)E(ω)can be written as
[ε0εre(ω)E(ω)]x=ε0Ex(ω)+ε0(εxx1)Ex(ω)+ε0εxyEy(ω)+ε0εxzEz(ω)=ε0Ex(ω)+Pxx(ω)+Pxy(ω)+Pxz(ω)
wherePxx(ω), Pxy(ω) and Pxz(ω)are the polarization variables [18]. Then, by substituting Eq. (2) into Eq. (1a), we have
jωε0Ex(ω)+jω[Pxx(ω)+Pxy(ω)+Pxz(ω)]=jωε0Ex(ω)+Jx(ω)=[×H(ω)]x
whereJx(ω)is the x-component of the polarization current density, and
Jx(ω)=jω[Pxx(ω)+Pxy(ω)+Pxz(ω)]
Similarly, we also have
Jy(ω)=jω[Pyx(ω)+Pyy(ω)+Pyz(ω)]
Jz(ω)=jω[Pzx(ω)+Pzy(ω)+Pzz(ω)].
Therefore, Eq. (1a) can be expressed as
jωε0E(ω)=×H(ω)-J(ω).
Converting Eq. (5) into the time domain by usingjω/t, we have
dE/dt=[(AB)HJ]/ε0
Where J=[JxJyJz]T,

A=[00/y/z000/x0],B=[0/z000/x/y00]. (6b,c)

Similarly, the Faraday’s induction law can also be rewritten as

dH/dt=(BA)E/μ0
where μ0 is the permeability in free space.

By applying the ADI process to Eq. (6a) and Eq. (7), we can obtain the ADI-FDTD formulations for the dispersive media as:

For the sub-step of n to n + 1/2,

En+1/2=En+(0.5Δt/ε0)(AHn+1/2BHnJn+p1)
Hn+1/2=Hn+(0.5Δt/μ0)(BEn+1/2AEn)

For the sub-step of n + 1/2 to n + 1,

En+1=En+1/2+(0.5Δt/ε0)(AHn+1/2BHn+1Jn+p2)
Hn+1=Hn+1/2+(0.5Δt/μ0)(BEn+1/2AEn+1)
where p1 and p2 are the time index within single time step, and they are within the ranges of [0, 1/2] and [1/2, 1], respectively. With the algebraic manipulations used in [25], the one-step leapfrog ADI-FDTD method with the polarization current density can be derived as follows.

First, by substituting Eq. (8b) into Eq. (8a), we have

En+1/2=En+Δt2ε0(AHn+Δt2μ0ABEn+1/2Δt2μ0AAEnBHnJn+p1).
By substituting Eq. (8d) into Eq. (8c), and replacing n with n-1, we obtain
En=En1/2+Δt2ε0(AHnΔt2μ0ABEn+1/2+Δt2μ0AAEnBHnJn+p21).
Adding Eq. (9) and Eq. (10) on the both corresponding sides, we obtain the iterative equation for the electric field in the leapfrog ADI-FDTD method as

(IΔt24μ0ε0AB)En+1/2=(IΔt24μ0ε0AB)En1/2+Δtε0(AHnBHn)Δt2ε0(Jn+p1+Jn+p21).

Similarly, an iterative equation for the magnetic field can be derived as follows. From Eq. (8a), it has

En=En+1/2(0.5Δt/ε0)(AHn+1/2BHnJn+p1).
By substituting Eq. (12) into Eq. (8b), we obtain
Hn+1/2=Hn+Δt2μ0(BEn+1/2+Δt2ε0AAHn+1/2Δt2ε0ABHnAEn+1/2Δt2ε0AJn+p1).
By substituting Eq. (8c) into Eq. (8d), we also obtain
Hn+1=Hn+1/2+Δt2μ0(BEn+1/2Δt2ε0AAHn+1/2+Δt2ε0ABHn+1AEn+1/2+Δt2ε0AJn+p2).
Then, add Eq. (13) and Eq. (14) on their both sides, the update equation for the magnetic field in the leapfrog ADI-FDTD method can be obtained as

(IΔt24μ0ε0AB)Hn+1=(IΔt24μ0ε0AB)Hn+Δtμ0(BEn+1/2AEn+1/2)+Δt24μ0ε0A(Jn+p2Jn+p1).

In practice, it is not easy to compute the cross term Δt2/4μ0ε0A(Jn+p2Jn+p1) in Eq. (15). To simplify it, we set p1 = p2 = 1/2 as

(IΔt24μ0ε0AB)En+1/2=(IΔt24μ0ε0AB)En1/2+Δtε0(AHnBHn)Δtε0Jn
(IΔt24μ0ε0AB)Hn+1=(IΔt24μ0ε0AB)Hn+Δtμ0(BEn+1/2AEn+1/2).
Here the Taylor series approximation Jn1/2+Jn+1/2=2Jn+O(Δt2) is used. Equations (16a) and (16b) are the equations of the leapfrog ADI-FDTD method applicable for modeling general dispersive media. The method is efficient since no mid-time step computation is needed with less memory consumed.

The above derivations of Eqs. (16a) and (16b) appear quite tedious. In fact, both of them can be reformulated in the form similar to that of the explicit FDTD method by two auxiliary variables e and h such that

En+1/2=En1/2+Δtε0(AHnBHn)Δtε0Jn+en
Hn+1=Hn+Δtμ0(BEn+1/2AEn+1/2)+hn+1/2
where
en=Δt24μ0ε0AB(En+1/2En1/2),hn+1/2=Δt24μ0ε0AB(Hn+1Hn). (17c,d)
Equations (17a) and (17b) show that the leapfrog ADI-FDTD can be regarded as the conventional FDTD method, with the second-order perturbation terms en and hn+1/2 added.

The above leapfrog ADI-FDTD formulas are applicable for various dispersive media as demonstrated below.

2.1 Lossy media

A lossy medium can be considered as a simple dispersive medium as studied in [32]. For simplicity, we choose its conductivity and equivalent permittivity as

σ=[σx000σy000σz],εre(ω)=[1+σx/jωε00001+σy/jωε00001+σz/jωε0]. (18a,b)
By substituting Eq. (18) into Eq. (2) and Eq. (4), we have
Jn=σEn=σ(En+1/2+En1/2)/2.
Furthermore, by substituting Eq. (19) into Eq. (16a), we obtain

(IΔt24μ0ε0AB+Δt2ε0σ)En+1/2=(IΔt24μ0ε0ABΔt2ε0σ)En1/2+Δtε0(AHnBHn).

It should be noted that Eq. (20) and Eq. (16b) have the same form as those in [32], and its unconditionally stable property has been studied in [34].

2.2 Magnetized plasma

Assume that the applied biasing magnetic field is in the z-axis direction, we then have [18]

εre(ω)=[εd(ω)jεg(ω)0jεg(ω)εd(ω)000εz(ω)]
εd(ω)=1+ωp2+νcωp2/(jω)(νc2+ωb2)+2jνcωω2,εg(ω)=ωbωp2/(jω)(νc2+ωb2)+2jνcωω2 (21b,c)
εz(ω)=1ωp2jνcω+ω2
where ωp is the plasma frequency, ωb is the cyclotron frequency, and νc is the electron collision rate. By substituting Eqs. (21a)-(21d) into Eq. (4), we obtain
Jx(ω)=ε0ωp2jω+ε0νcωp2(νc2+ωb2)+2jνcωω2Ex+ε0ωbωp2(νc2+ωb2)+2jνcωω2Ey.
As discussed in [18], the first and second terms on the right-hand side of Eq. (22) can be replaced by two auxiliary variables Jxx and Jxy, i.e. Jxn=Jxxn+Jxyn. Both Jxx(ω) and Jxy(ω) have the general form as
J(ω)=a0+a1jωb0+b1jω+b2(jω)2E(ω)
wherea0, a1, b0, b1and b2 are the real coefficients. Converting Eq. (23) into the time domain, we have
b0J+b1dJdt+b2d2Jdt2=a0E+a1dEdt.
By discretizing and approximating the derivatives around the time step n-1, we obtain
b0Jn1+b1JnJn22Δt+b2Jn2Jn1+Jn2Δt2=a0En1/2+En3/22+a1En1/2En3/2Δt.
and then, one iterative equation for J is derived as

Jn=4b2Δt2b0Δt2b1Δt+2b2Jn1+b1Δt2b2b1Δt+2b2Jn2+a0Δt2+2a1Δtb1Δt+2b2En1/2+a0Δt22a1Δtb1Δt+2b2En3/2.

2.3 Graphene biased by a magnetostatic field

The monolayer graphene is biased by a magnetostatic field Bo in the z-axis direction, and it is placed in the x-y plane. Then, its volume conductivityσcan be calculated from its surface conductivity σsby usingσ=σs/Δz. In the frequency up to 10 THz, its σscan be described by a 2-D Drude model [35], withεre(ω)=I+σ/(jωε0), and

εre(ω)=[εd(ω)jεg(ω)0jεg(ω)εd(ω)0000]
εd(ω)=1+σ0+νσ0/(jω)(ν2+ωc2)+2jνωω2
εg(ω)=ωbσ0/(jω)(ν2+ωc2)+2jνωω2
σ0=2e2τkBTνπ2ε0Δzln(2coshμc2kBT)
where ν is the phenomenological scattering rate,ωcis the cyclotron frequency and ωc=eB0vF2/μc, e is the electron charge, vF is the Fermi velocity, kB is the Boltzmann constant, μc is the chemical potential, T is the operating temperature, and is the reduced Planck’s constant.

According to Eqs. (27a)-(27d), it is found that theεre(ω)of graphene sheet can be obtained by replacing the parameters of magnetized plasma in Eq. (21) withωp2σ0, νcν, ωbωc, andεz(ω)=0. Therefore, the magnetized graphene can also be simulated using Eq. (16a), Eq. (16b), and Eq. (26).

2.4 Graphene biased by an electrostatic field

If the monolayer graphene sheet is only biased by an electrostatic fieldE0, its relative permittivity can be obtained by setting ωc=0 in Eq. (27), i.e.

εd(ω)=1+[σ0/(jω)]/(ν+jω)
and εg(ω)=0. With the iterative equation Eq. (26) used, both Jx and Jy can be obtained easily.

In our simulation below, the chemical potential μcof graphene is calculated by

E0=[e/(ε0π2vF2)]0ε[fd(ε)fd(ε+2μc)]dε
wherefd(ε)is the Femi-Dirac distribution function.

3. Numerical results and discussion

To demonstrate the capability of our developed leapfrog ADI-FDTD algorithm for treating general dispersive media, three numerical experiments were performed: electromagnetic wave propagation in a PEC waveguide loaded with a magnetized plasma slab [33], electromagnetic wave transmission through a magnetized two-dimensional monolayer graphene sheet, and the surface plasmon polaritons along a two-dimensional graphene sheet.

3.1 The PEC Waveguide Loaded with Magnetized Plasma

Figure 1 shows the rectangular waveguide of 40mm × 40mm, inserted with a magnetized plasma slab of 5 mm in width. The waveguide is terminated by the 10-layer CPMLs at both ends [33]. In such a structure, both propagating and evanescent modes exist.

 

Fig. 1 A rectangular waveguide of 40mm × 40mm loaded with a magnetized plasma slab of 5 mm in width.

Download Full Size | PPT Slide | PDF

In our simulation, we set ∆x = ∆y = ∆z = ∆ = 1mm. A modulated Gaussian pulse ofcos(2πf0t)exp[4π(tt0)2/τ2]was used to excite electromagnetic field, where f0 = 10 GHz,Δt=Δ/3c0CFLN,τ=160Δt/CFLN, and t0=τ. CFLN is the CFL number which measures the ratio of the time step to the CFL limit. The TE10-mode was excited and its cutoff frequency is 3.75GHz. The plasma frequency ωp = 2π × 10 Grad/s, ν = 10 GHz, and ωb = 10 Grad/s. An observation point is set at the waveguide center, and it is ten cells away from the magnetized plasma slab. Figure 2 shows the recorded Ez -component at the observation point verses time with different CFLNs. It is seen that for CFLN = 1, the recorded Ez –component agrees well with that of the conventional FDTD. As the CFLN increases from 1 to 3, the errors are increased slightly, but they are in an acceptable range. To check the stability of our developed algorithm, we run the simulation to one million time steps with the time step CFLN = 2 and 3, respectively; the recorded field was found to be always stable.

 

Fig. 2 The Recorded Ez –component versus time with CFLN = 1, 2, and 3; the result obtained with the conventional FDTD method is also plotted for comparison.

Download Full Size | PPT Slide | PDF

3.2 Magnetized graphene sheet

When a graphene sheet is biased by a magnetostatic field B0, it produces gyrotropy [35]. By applying the proposed method, we can predict the transmission properties of an infinite graphene sheet at both THz and optical bands.

In our simulation, a two-dimensional graphene sheet was placed at the center of z-direction in the x-y plane as in Fig. 3(a). We chose the maximum frequency fmax = c0/λmin = 2 THz, the FDTD cells ∆x = ∆y = ∆z = ∆ = λmin/20, ν = 20THz, T = 300K, μc = 0.1 eV, and Δt=Δ/(3c0). The overall computational domain consists of 40 × 40 × 180 cells, including 10-layer CPML [28]. The normally incident plane wave is polarized in the x-direction and described by

Exinc(t)=exp[4π(tt0)2/τ2]
where τ=80Δt and t0=2τ. The observation point was set at the center of x-y plane and 4 cells away from the graphene sheet. The analytical transmission coefficient of the infinite graphene sheet can be calculated as [35]:
Tanalytic(f)=2|2+jω[εdω1]η0|2+|ε0εgω2η0|2|(2+jω[εdω1]η0)2+(ε0εgω2η0)2|
where η0is the wave impedance in the free space. The transmission coefficient in the FDTD simulation was computed with
Tleapfrog(f)=|FFT(Ex,ob(t))|2+|FFT(Ey,ob(t))|2|FFT(Ein(t))|
where FFT represents the fast Fourier transformation.

 

Fig. 3 (a) The computational domain where a plane wave normally incident on an infinite graphene sheet. (b) The computed transmission coefficient |T| of the plane wave normally incident on an infinite graphene sheet biased by a magnetostatic field B0 = 1 T. The analytical result is also plotted for comparison.

Download Full Size | PPT Slide | PDF

Figure 3(b) shows the computed transmission coefficient using the leapfrog ADI-FDTD with B0 = 1T and CFLN = 1, 2, and 3, respectively. The analytical result is also plotted for comparison. It is seen that, when CFLN = 1, the transmission coefficient |T| and the analytical solution agrees very well up to 2 THz. When CFLN = 2, the value of |T| still agrees well with the analytical solution at low frequencies. When CFLN = 3, the errors become much larger at high frequencies. This means that the leapfrog ADI-FDTD is quite suitable for computing the transmission coefficient at low frequencies.

Figures 4(a) and 4(b) show the snapshots of Ey – component at time t = 2.5t0 for CFLN = 1 and 3, respectively. The field is recorded at 10 cells away from the CPML in the y-z plane. It is seen that the field distribution computed with CFLN = 3 is slightly different from that of CFLN = 1. On the other hand, it is noted that, as the incident EM wave is polarized in the x-axis direction, the Ey – component remains zero before it interacts with the graphene sheet. Therefore, Fig. 4 shows the gyrotropic property of the magnetized graphene sheet.

 

Fig. 4 Ey-component snapshot at t = 2.5t0 for (a) CFLN = 1 and (b) CFLN = 3, respectively. The field is recorded at 10 cells away from the CPML in the y-z plane.

Download Full Size | PPT Slide | PDF

3.3 SPPs propagating along the graphene sheet

It is known that whenRe[εd(ω)]<0, the graphene sheet acts as a thin metal film, and a TM-mode SPPs can be excited. The properties of SPPs have been studied using the finite-difference frequency-domain (FDFD) and the conventional FDTD [36,37], respectively. However, as it is an atomically thin structure, both methods may take very long simulation time because of the small cell sizes and small time steps used. Here, we fast solve it using the proposed method with a large CFLN.

In our simulation, we chose ∆y = ∆z = ∆ = 50 nm, ν = 1 THz, T = 300 K, μc = 0.15 eV,and Δt=Δ/(3c0). The overall computational domain, as shown in Fig. 5, consists of 820 × 180 cells, including 10-layer CPML [31]. A single point source, placed one cell over the graphene sheet and one cell away from the CPML, is used to excite the SPPs, and the exciting magnetic current source isJm,x(t)=sin(2πf0t)with f0 = 6 THz. Therefore, the grid resolution is λ0/∆ = c0/(f0∆) = 1000, and for such a fine structure the large CFLN is needed so as to reduce the simulation time.

 

Fig. 5 Computational domain for the simulation of SPPs on the graphene sheet.

Download Full Size | PPT Slide | PDF

Figures 6(a)-6(d) show snapshots of the Ez – and Ey –components at time t = 2.887 ps. In Figs. 6(a) and 6(b), the numerical results were obtained with the leapfrog ADI-FDTD with CFLN = 10, and in Figs. 6(c) and 6(d), they were computed using the conventional ADE-FDTD [17] with CFLN = 1. It is seen that the results obtained with the leapfrog ADI-FDTD agree well with that of the conventional FDTD. Its computational time is about 46 seconds for 3000 time steps with CFLN = 10, while the conventional FDTD uses 142 seconds for 30000 time steps. Our simulation platform is the Lenovo PC with Intel Dual Core P8700 of 2.53 GHz and RAM of 8GB.

 

Fig. 6 Snapshots of the normalized electric field at time t = 2.887 ps. (a) Ez: leapfrog ADI-FDTD with CFLN = 10; (b) Ey: leapfrog ADI-FDTD with CFLN = 10; (c) Ez: conventional ADE-FDTD with CFLN = 1 [17]; and (d) Ey: conventional ADE-FDTD with CFLN = 1 [17].

Download Full Size | PPT Slide | PDF

4. Conclusion

In this paper, the one-step leapfrog ADI-FDTD method is adapted for simulating general dispersive media. In the method, dispersive properties are modelled with equivalent polarization currents which are then solved by the newly developed ADE formulations. Numerical experiments have shown that the proposed method can simulate the conventional magnetized plasma as well as magnetized graphene sheet accurately. In particular, it can be used to simulate the SPPs on a monolayer graphene sheet and the simulation time was reduced greatly without sacrificing the accuracy.

Acknowledgments

This work was supported in part by the NSFC under Grant 61171037, the Specialized Research Fund for the Doctoral Program of Higher Education under Grant 20100101110065, the State Key Lab of MOI of Zhejiang University, the National Basic Research Program under Grant 2009CB320204 of China, and the Natural Science and Engineering Research Council of Canada through its Discovery Grant program.

References and links

1. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, 2005).

2. C. Lundgren, R. Lopez, J. Redwing, and K. Melde, “FDTD modeling of solar energy absorption in silicon branched nanowires,” Opt. Express 21(9), 329–400 (2013). [PubMed]  

3. E. H. Khoo, I. Ahmed, R. S. M. Goh, K. H. Lee, T. G. G. Hung, and E. P. Li, “Efficient analysis of mode profiles in elliptical microcavity using dynamic-thermal electron-quantum medium FDTD method,” Opt. Express 21(5), 5910–5923 (2013). [CrossRef]   [PubMed]  

4. D. T. Nguyen and R. A. Norwood, “Label-free, single-object sensing with a microring resonator: FDTD simulation,” Opt. Express 21(1), 49–59 (2013). [CrossRef]   [PubMed]  

5. I. R. Çapoğlu, A. Taflove, and V. Backman, “Computation of tightly-focused laser beams in the FDTD method,” Opt. Express 21(1), 87–101 (2013). [CrossRef]   [PubMed]  

6. S. Buil, J. Laverdant, B. Berini, P. Maso, J. P. Hermier, and X. Quélin, “FDTD simulations of localization and enhancements on fractal plasmonics nanostructures,” Opt. Express 20(11), 11968–11975 (2012). [CrossRef]   [PubMed]  

7. C. M. Dissanayake, M. Premaratne, I. D. Rukhlenko, and G. P. Agrawal, “FDTD modeling of anisotropic nonlinear optical phenomena in silicon waveguides,” Opt. Express 18(20), 21427–21448 (2010). [CrossRef]   [PubMed]  

8. G. Singh, K. Ravi, Q. Wang, and S. T. Ho, “Complex-envelope alternating-direction-implicit FDTD method for simulating active photonic devices with semiconductor/solid-state media,” Opt. Lett. 37(12), 2361–2363 (2012). [CrossRef]   [PubMed]  

9. G. Singh, E. L. Tan, and Z. N. Chen, “Split-step finite-difference time-domain method with perfectly matched layers for efficient analysis of two-dimensional photonic crystals with anisotropic media,” Opt. Lett. 37(3), 326–328 (2012). [CrossRef]   [PubMed]  

10. S. Zhao, “High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media,” Opt. Lett. 36(16), 3245–3247 (2011). [CrossRef]   [PubMed]  

11. W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997). [CrossRef]  

12. D. Sullivan, “Nonlinear FDTD formulations using Z transforms,” IEEE Trans. Microw. Theory Tech. 43(3), 676–682 (1995). [CrossRef]  

13. R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992). [CrossRef]  

14. M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006). [CrossRef]  

15. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006). [CrossRef]   [PubMed]  

16. L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012). [CrossRef]  

17. M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009). [CrossRef]  

18. A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013). [CrossRef]  

19. L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006). [CrossRef]  

20. D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996). [CrossRef]  

21. Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010). [CrossRef]  

22. S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008). [CrossRef]  

23. D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009). [CrossRef]  

24. S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005). [CrossRef]  

25. S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009). [CrossRef]  

26. F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000). [CrossRef]  

27. E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antenn. Propag. 56(1), 170–177 (2008). [CrossRef]  

28. X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012). [CrossRef]  

29. T. H. Gan and E. T. Tan, “Mur absorbing boundary condition for 2-D leapfrog ADI-FDTD method,” in proceedings of IEEE Conference on Asia-Pacific Antennas and Propagation, 3–4 (2012). [CrossRef]  

30. Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013). [CrossRef]  

31. T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012). [CrossRef]  

32. T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).

33. X. H. Wang, W. Y. Yin, and Z. Chen, “One-step leapfrog ADI-FDTD method for anisotropic magnetized plasma,” in proceedings of IEEE International Microwave Symposium, 1–4 (2013).

34. J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).

35. D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012). [CrossRef]  

36. B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012). [CrossRef]  

37. H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, 2005).
  2. C. Lundgren, R. Lopez, J. Redwing, and K. Melde, “FDTD modeling of solar energy absorption in silicon branched nanowires,” Opt. Express 21(9), 329–400 (2013).
    [PubMed]
  3. E. H. Khoo, I. Ahmed, R. S. M. Goh, K. H. Lee, T. G. G. Hung, and E. P. Li, “Efficient analysis of mode profiles in elliptical microcavity using dynamic-thermal electron-quantum medium FDTD method,” Opt. Express 21(5), 5910–5923 (2013).
    [Crossref] [PubMed]
  4. D. T. Nguyen and R. A. Norwood, “Label-free, single-object sensing with a microring resonator: FDTD simulation,” Opt. Express 21(1), 49–59 (2013).
    [Crossref] [PubMed]
  5. I. R. Çapoğlu, A. Taflove, and V. Backman, “Computation of tightly-focused laser beams in the FDTD method,” Opt. Express 21(1), 87–101 (2013).
    [Crossref] [PubMed]
  6. S. Buil, J. Laverdant, B. Berini, P. Maso, J. P. Hermier, and X. Quélin, “FDTD simulations of localization and enhancements on fractal plasmonics nanostructures,” Opt. Express 20(11), 11968–11975 (2012).
    [Crossref] [PubMed]
  7. C. M. Dissanayake, M. Premaratne, I. D. Rukhlenko, and G. P. Agrawal, “FDTD modeling of anisotropic nonlinear optical phenomena in silicon waveguides,” Opt. Express 18(20), 21427–21448 (2010).
    [Crossref] [PubMed]
  8. G. Singh, K. Ravi, Q. Wang, and S. T. Ho, “Complex-envelope alternating-direction-implicit FDTD method for simulating active photonic devices with semiconductor/solid-state media,” Opt. Lett. 37(12), 2361–2363 (2012).
    [Crossref] [PubMed]
  9. G. Singh, E. L. Tan, and Z. N. Chen, “Split-step finite-difference time-domain method with perfectly matched layers for efficient analysis of two-dimensional photonic crystals with anisotropic media,” Opt. Lett. 37(3), 326–328 (2012).
    [Crossref] [PubMed]
  10. S. Zhao, “High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media,” Opt. Lett. 36(16), 3245–3247 (2011).
    [Crossref] [PubMed]
  11. W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997).
    [Crossref]
  12. D. Sullivan, “Nonlinear FDTD formulations using Z transforms,” IEEE Trans. Microw. Theory Tech. 43(3), 676–682 (1995).
    [Crossref]
  13. R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992).
    [Crossref]
  14. M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
    [Crossref]
  15. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
    [Crossref] [PubMed]
  16. L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
    [Crossref]
  17. M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009).
    [Crossref]
  18. A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
    [Crossref]
  19. L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006).
    [Crossref]
  20. D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996).
    [Crossref]
  21. Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010).
    [Crossref]
  22. S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008).
    [Crossref]
  23. D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009).
    [Crossref]
  24. S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005).
    [Crossref]
  25. S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
    [Crossref]
  26. F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
    [Crossref]
  27. E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antenn. Propag. 56(1), 170–177 (2008).
    [Crossref]
  28. X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
    [Crossref]
  29. T. H. Gan and E. T. Tan, “Mur absorbing boundary condition for 2-D leapfrog ADI-FDTD method,” in proceedings of IEEE Conference on Asia-Pacific Antennas and Propagation, 3–4 (2012).
    [Crossref]
  30. Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
    [Crossref]
  31. T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012).
    [Crossref]
  32. T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).
  33. X. H. Wang, W. Y. Yin, and Z. Chen, “One-step leapfrog ADI-FDTD method for anisotropic magnetized plasma,” in proceedings of IEEE International Microwave Symposium, 1–4 (2013).
  34. J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).
  35. D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012).
    [Crossref]
  36. B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
    [Crossref]
  37. H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
    [Crossref]

2013 (7)

C. Lundgren, R. Lopez, J. Redwing, and K. Melde, “FDTD modeling of solar energy absorption in silicon branched nanowires,” Opt. Express 21(9), 329–400 (2013).
[PubMed]

E. H. Khoo, I. Ahmed, R. S. M. Goh, K. H. Lee, T. G. G. Hung, and E. P. Li, “Efficient analysis of mode profiles in elliptical microcavity using dynamic-thermal electron-quantum medium FDTD method,” Opt. Express 21(5), 5910–5923 (2013).
[Crossref] [PubMed]

D. T. Nguyen and R. A. Norwood, “Label-free, single-object sensing with a microring resonator: FDTD simulation,” Opt. Express 21(1), 49–59 (2013).
[Crossref] [PubMed]

I. R. Çapoğlu, A. Taflove, and V. Backman, “Computation of tightly-focused laser beams in the FDTD method,” Opt. Express 21(1), 87–101 (2013).
[Crossref] [PubMed]

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).

2012 (10)

D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012).
[Crossref]

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012).
[Crossref]

T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

S. Buil, J. Laverdant, B. Berini, P. Maso, J. P. Hermier, and X. Quélin, “FDTD simulations of localization and enhancements on fractal plasmonics nanostructures,” Opt. Express 20(11), 11968–11975 (2012).
[Crossref] [PubMed]

G. Singh, K. Ravi, Q. Wang, and S. T. Ho, “Complex-envelope alternating-direction-implicit FDTD method for simulating active photonic devices with semiconductor/solid-state media,” Opt. Lett. 37(12), 2361–2363 (2012).
[Crossref] [PubMed]

G. Singh, E. L. Tan, and Z. N. Chen, “Split-step finite-difference time-domain method with perfectly matched layers for efficient analysis of two-dimensional photonic crystals with anisotropic media,” Opt. Lett. 37(3), 326–328 (2012).
[Crossref] [PubMed]

2011 (1)

2010 (2)

C. M. Dissanayake, M. Premaratne, I. D. Rukhlenko, and G. P. Agrawal, “FDTD modeling of anisotropic nonlinear optical phenomena in silicon waveguides,” Opt. Express 18(20), 21427–21448 (2010).
[Crossref] [PubMed]

Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010).
[Crossref]

2009 (3)

M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009).
[Crossref]

D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009).
[Crossref]

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

2008 (2)

E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antenn. Propag. 56(1), 170–177 (2008).
[Crossref]

S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008).
[Crossref]

2006 (3)

L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006).
[Crossref]

M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
[Crossref]

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
[Crossref] [PubMed]

2005 (1)

S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005).
[Crossref]

2000 (1)

F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
[Crossref]

1997 (1)

W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997).
[Crossref]

1996 (1)

D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996).
[Crossref]

1995 (1)

D. Sullivan, “Nonlinear FDTD formulations using Z transforms,” IEEE Trans. Microw. Theory Tech. 43(3), 676–682 (1995).
[Crossref]

1992 (1)

R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992).
[Crossref]

Agrawal, G. P.

Ahmed, I.

Al-Jabr, A.

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009).
[Crossref]

Alsunaidi, M.

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009).
[Crossref]

Alvarez, J.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Angulo, L. D.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Antonsen, T. M.

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

Backman, V.

Berini, B.

Botton, M.

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

Buil, S.

Caloz, C.

D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012).
[Crossref]

Çapoglu, I. R.

Chen, B.

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

Chen, J.

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

Chen, Z.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
[Crossref]

Chen, Z. N.

Cooke, S. J.

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

Dissanayake, C. M.

Dutton, R.

M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
[Crossref]

Etchegoin, P. G.

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
[Crossref] [PubMed]

Fan, S.

M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
[Crossref]

Gan, T. H.

T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012).
[Crossref]

T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).

Gao, J. Y.

J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).

Garcia, S. G.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Goh, R. S. M.

Guo, Y.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

Han, L.

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

Han, M.

M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
[Crossref]

Heh, D. Y.

D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009).
[Crossref]

Hermier, J. P.

Ho, S. T.

Huang, S.

S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005).
[Crossref]

Huang, W. P.

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

Hung, T. G. G.

Hunsberger, F.

R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992).
[Crossref]

Kelley, D. F.

D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996).
[Crossref]

Khee, T.

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

Khoo, E. H.

Laverdant, J.

Le Ru, E. C.

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
[Crossref] [PubMed]

Lee, K. H.

Levush, B.

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

Li, E. P.

Li, F.

S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005).
[Crossref]

Li, K.

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

Li, X.

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

Lin, H.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Liu, S.

S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008).
[Crossref]

S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008).
[Crossref]

Lopez, R.

Luebbers, R. J.

D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996).
[Crossref]

R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992).
[Crossref]

Lundgren, C.

Mao, Y. F.

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

Martin, R. G.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Maso, P.

Melde, K.

Meyer, M.

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
[Crossref] [PubMed]

Nguyen, D. T.

Norwood, R. A.

Ooi, B. S.

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

Pantoja, M. F.

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

Premaratne, M.

Quélin, X.

Rappaport, C. M.

W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997).
[Crossref]

Ravi, K.

Redwing, J.

Rukhlenko, I. D.

Simpson, J.

Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010).
[Crossref]

Singh, G.

Sounas, D. L.

D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012).
[Crossref]

Sullivan, D.

D. Sullivan, “Nonlinear FDTD formulations using Z transforms,” IEEE Trans. Microw. Theory Tech. 43(3), 676–682 (1995).
[Crossref]

Taflove, A.

Tan, E. L.

G. Singh, E. L. Tan, and Z. N. Chen, “Split-step finite-difference time-domain method with perfectly matched layers for efficient analysis of two-dimensional photonic crystals with anisotropic media,” Opt. Lett. 37(3), 326–328 (2012).
[Crossref] [PubMed]

T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).

T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012).
[Crossref]

D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009).
[Crossref]

E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antenn. Propag. 56(1), 170–177 (2008).
[Crossref]

Tang, J. Z.

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

Teng, J.

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

Wang, B.

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

Wang, J.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

Wang, Q.

Wang, X. H.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

Weedon, W. H.

W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997).
[Crossref]

Xia, J. L.

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

Xu, L. J.

L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006).
[Crossref]

Yin, W. Y.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

Yu, Y.

Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010).
[Crossref]

Yu, Y. Q.

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

Yuan, N. C.

L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006).
[Crossref]

Yuan, X.

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

Zhang, J.

F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
[Crossref]

Zhang, X.

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

Zhao, S.

Zheng, F.

F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
[Crossref]

Zheng, H. X.

J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).

Zhou, D.

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

Appl. Phys. Lett. (1)

B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012).
[Crossref]

IEEE Antennas Wirel. Propag. Lett. (3)

Y. F. Mao, B. Chen, J. L. Xia, J. Chen, and J. Z. Tang, “Application of the leapfrog ADI-FDTD method to periodic structures,” IEEE Antennas Wirel. Propag. Lett. 12, 599–602 (2013).
[Crossref]

S. Liu and S. Liu, “Runge-Kutta exponential time differencing FDTD method for anisotropic magnetized plasma,” IEEE Antennas Wirel. Propag. Lett. 7, 306–309 (2008).
[Crossref]

L. J. Xu and N. C. Yuan, “FDTD formulations for scattering from 3-D anisotropic magnetized plasma objects,” IEEE Antennas Wirel. Propag. Lett. 5(1), 335–338 (2006).
[Crossref]

IEEE J. Light. Tech. (1)

L. Han, D. Zhou, K. Li, X. Li, and W. P. Huang, “A rational-fraction dispersion model for efficient simulation of dispersive material in FDTD method,” IEEE J. Light. Tech. 30(13), 2216–2225 (2012).
[Crossref]

IEEE Microw. Wirel. Compon. Lett. (4)

M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” IEEE Microw. Wirel. Compon. Lett. 16(3), 119–121 (2006).
[Crossref]

D. Y. Heh and E. L. Tan, “FDTD modeling for dispersive media using matrix exponential method,” IEEE Microw. Wirel. Compon. Lett. 19(2), 53–55 (2009).
[Crossref]

S. Huang and F. Li, “FDTD implementation for magnetoplasma medium using exponential time differencing,” IEEE Microw. Wirel. Compon. Lett. 15(3), 183–185 (2005).
[Crossref]

H. Lin, M. F. Pantoja, L. D. Angulo, J. Alvarez, R. G. Martin, and S. G. Garcia, “FDTD modeling of graphene devices using complex conjugate dispersion material model,” IEEE Microw. Wirel. Compon. Lett. 22(12), 612–614 (2012).
[Crossref]

IEEE Photon. Technol. Lett. (1)

M. Alsunaidi and A. Al-Jabr, “A general ADE-FDTD algorithm for the simulation of dispersive structures,” IEEE Photon. Technol. Lett. 21(12), 817–819 (2009).
[Crossref]

IEEE Trans. Antenn. Propag. (7)

A. Al-Jabr, M. Alsunaidi, T. Khee, and B. S. Ooi, “A simple FDTD algorithm for simulating EM-wave propagation in general dispersive anisotropic material,” IEEE Trans. Antenn. Propag. 61(3), 1321–1326 (2013).
[Crossref]

W. H. Weedon and C. M. Rappaport, “A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media,” IEEE Trans. Antenn. Propag. 45(3), 401–410 (1997).
[Crossref]

R. J. Luebbers and F. Hunsberger, “FDTD Nth-order dispersive media,” IEEE Trans. Antenn. Propag. 40(11), 1297–1301 (1992).
[Crossref]

T. H. Gan and E. L. Tan, “Analysis of the divergence properties for the three-dimensional leapfrog ADI-FDTD method,” IEEE Trans. Antenn. Propag. 60(12), 5801–5808 (2012).
[Crossref]

D. F. Kelley and R. J. Luebbers, “Piecewise linear recursive convolution for dispersive media using FDTD,” IEEE Trans. Antenn. Propag. 44(6), 792–797 (1996).
[Crossref]

Y. Yu and J. Simpson, “An EJ collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma,” IEEE Trans. Antenn. Propag. 58(2), 469–478 (2010).
[Crossref]

E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antenn. Propag. 56(1), 170–177 (2008).
[Crossref]

IEEE Trans. Electromagn. Compat. (1)

X. H. Wang, W. Y. Yin, Y. Q. Yu, Z. Chen, J. Wang, and Y. Guo, “A convolutional perfect matched layer (CPML) for one-step leapfrog ADI-FDTD method and its applications to EMC problems,” IEEE Trans. Electromagn. Compat. 54(5), 1066–1076 (2012).
[Crossref]

IEEE Trans. Microw. Theory Tech. (3)

F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three dimensional unconditionally stable finite-different time-domain method,” IEEE Trans. Microw. Theory Tech. 48(9), 1550–1558 (2000).
[Crossref]

D. Sullivan, “Nonlinear FDTD formulations using Z transforms,” IEEE Trans. Microw. Theory Tech. 43(3), 676–682 (1995).
[Crossref]

D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microw. Theory Tech. 60(4), 901–914 (2012).
[Crossref]

Int. J. Numer. Model. (1)

S. J. Cooke, M. Botton, T. M. Antonsen, and B. Levush, “A leapfrog formulation of the 3D ADI-FDTD algorithm,” Int. J. Numer. Model. 22(2), 187–200 (2009).
[Crossref]

J. Chem. Phys. (1)

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125(16), 164705 (2006).
[Crossref] [PubMed]

Opt. Express (6)

Opt. Lett. (3)

Progress In Electromagnetics Research Letters (1)

J. Y. Gao and H. X. Zheng, “One-Step leapfrog ADI-FDTD method for lossy media and its stability analysis,” Progress In Electromagnetics Research Letters 40, 49–60 (2013).

Progress In Electromagnetics Research M (1)

T. H. Gan and E. L. Tan, “Unconditionally stable leapfrog ADI-FDTD method for lossy media,” Progress In Electromagnetics Research M 26, 173–186 (2012).

Other (3)

X. H. Wang, W. Y. Yin, and Z. Chen, “One-step leapfrog ADI-FDTD method for anisotropic magnetized plasma,” in proceedings of IEEE International Microwave Symposium, 1–4 (2013).

T. H. Gan and E. T. Tan, “Mur absorbing boundary condition for 2-D leapfrog ADI-FDTD method,” in proceedings of IEEE Conference on Asia-Pacific Antennas and Propagation, 3–4 (2012).
[Crossref]

A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, 2005).

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1

A rectangular waveguide of 40mm × 40mm loaded with a magnetized plasma slab of 5 mm in width.

Fig. 2
Fig. 2

The Recorded Ez –component versus time with CFLN = 1, 2, and 3; the result obtained with the conventional FDTD method is also plotted for comparison.

Fig. 3
Fig. 3

(a) The computational domain where a plane wave normally incident on an infinite graphene sheet. (b) The computed transmission coefficient |T| of the plane wave normally incident on an infinite graphene sheet biased by a magnetostatic field B0 = 1 T. The analytical result is also plotted for comparison.

Fig. 4
Fig. 4

Ey-component snapshot at t = 2.5t0 for (a) CFLN = 1 and (b) CFLN = 3, respectively. The field is recorded at 10 cells away from the CPML in the y-z plane.

Fig. 5
Fig. 5

Computational domain for the simulation of SPPs on the graphene sheet.

Fig. 6
Fig. 6

Snapshots of the normalized electric field at time t = 2.887 ps. (a) Ez: leapfrog ADI-FDTD with CFLN = 10; (b) Ey: leapfrog ADI-FDTD with CFLN = 10; (c) Ez: conventional ADE-FDTD with CFLN = 1 [17]; and (d) Ey: conventional ADE-FDTD with CFLN = 1 [17].

Equations (47)

Equations on this page are rendered with MathJax. Learn more.

jω ε 0 ε re E(ω)=×H(ω)
ε re (ω)=[ ε xx (ω) ε xy (ω) ε xz (ω) ε yx (ω) ε yy (ω) ε yz (ω) ε zx (ω) ε zy (ω) ε zz (ω) ]
[ ε 0 ε re (ω)E(ω) ] x = ε 0 E x (ω)+ ε 0 ( ε xx 1) E x (ω)+ ε 0 ε xy E y (ω)+ ε 0 ε xz E z (ω) = ε 0 E x (ω)+ P xx (ω)+ P xy (ω)+ P xz (ω)
jω ε 0 E x (ω)+jω[ P xx (ω)+ P xy (ω)+ P xz (ω)]=jω ε 0 E x (ω)+ J x (ω)= [ ×H(ω) ] x
J x (ω)=jω[ P xx (ω)+ P xy (ω)+ P xz (ω)]
J y (ω)=jω[ P yx (ω)+ P yy (ω)+ P yz (ω)]
J z (ω)=jω[ P zx (ω)+ P zy (ω)+ P zz (ω)].
jω ε 0 E(ω)=×H(ω) -J(ω).
dE/dt=[ (AB)HJ ]/ ε 0
A=[ 0 0 /y /z 0 0 0 /x 0 ], B=[ 0 /z 0 0 0 /x /y 0 0 ].
dH/dt=(BA)E/ μ 0
E n+1/2 = E n +(0.5Δt/ ε 0 )(A H n+1/2 B H n J n+ p 1 )
H n+1/2 = H n +(0.5Δt/ μ 0 )(B E n+1/2 A E n )
E n+1 = E n+1/2 +(0.5Δt/ ε 0 )(A H n+1/2 B H n+1 J n+ p 2 )
H n+1 = H n+1/2 +(0.5Δt/ μ 0 )(B E n+1/2 A E n+1 )
E n+1/2 = E n + Δt 2 ε 0 (A H n + Δt 2 μ 0 AB E n+1/2 Δt 2 μ 0 AA E n B H n J n+ p 1 ).
E n = E n1/2 + Δt 2 ε 0 (A H n Δt 2 μ 0 AB E n+1/2 + Δt 2 μ 0 AA E n B H n J n+ p 2 1 ).
(I Δ t 2 4 μ 0 ε 0 AB) E n+1/2 =(I Δ t 2 4 μ 0 ε 0 AB) E n1/2 + Δt ε 0 (A H n B H n ) Δt 2 ε 0 ( J n+ p 1 + J n+ p 2 1 ).
E n = E n+1/2 (0.5Δt/ ε 0 )(A H n+1/2 B H n J n+ p 1 ).
H n+1/2 = H n + Δt 2 μ 0 (B E n+1/2 + Δt 2 ε 0 AA H n+1/2 Δt 2 ε 0 AB H n A E n+1/2 Δt 2 ε 0 A J n+ p 1 ).
H n+1 = H n+1/2 + Δt 2 μ 0 (B E n+1/2 Δt 2 ε 0 AA H n+1/2 + Δt 2 ε 0 AB H n+1 A E n+1/2 + Δt 2 ε 0 A J n+ p 2 ).
(I Δ t 2 4 μ 0 ε 0 AB) H n+1 =(I Δ t 2 4 μ 0 ε 0 AB) H n + Δt μ 0 (B E n+1/2 A E n+1/2 ) + Δ t 2 4 μ 0 ε 0 A( J n+ p 2 J n+ p 1 ).
(I Δ t 2 4 μ 0 ε 0 AB) E n+1/2 =(I Δ t 2 4 μ 0 ε 0 AB) E n1/2 + Δt ε 0 (A H n B H n ) Δt ε 0 J n
(I Δ t 2 4 μ 0 ε 0 AB) H n+1 =(I Δ t 2 4 μ 0 ε 0 AB) H n + Δt μ 0 (B E n+1/2 A E n+1/2 ).
E n+1/2 = E n1/2 + Δt ε 0 (A H n B H n ) Δt ε 0 J n + e n
H n+1 = H n + Δt μ 0 (B E n+1/2 A E n+1/2 )+ h n+1/2
e n = Δ t 2 4 μ 0 ε 0 AB( E n+1/2 E n1/2 ), h n+1/2 = Δ t 2 4 μ 0 ε 0 AB( H n+1 H n ).
σ=[ σ x 0 0 0 σ y 0 0 0 σ z ], ε re (ω)=[ 1+ σ x /jω ε 0 0 0 0 1+ σ y /jω ε 0 0 0 0 1+ σ z /jω ε 0 ].
J n =σ E n =σ( E n+1/2 + E n1/2 )/2.
(I Δ t 2 4 μ 0 ε 0 AB+ Δt 2 ε 0 σ) E n+1/2 =(I Δ t 2 4 μ 0 ε 0 AB Δt 2 ε 0 σ) E n1/2 + Δt ε 0 (A H n B H n ).
ε re (ω)=[ ε d (ω) j ε g (ω) 0 j ε g (ω) ε d (ω) 0 0 0 ε z (ω) ]
ε d (ω)=1+ ω p 2 + ν c ω p 2 /(jω) ( ν c 2 + ω b 2 )+2j ν c ω ω 2 , ε g (ω)= ω b ω p 2 /(jω) ( ν c 2 + ω b 2 )+2j ν c ω ω 2
ε z (ω)=1 ω p 2 j ν c ω+ ω 2
J x (ω)= ε 0 ω p 2 jω+ ε 0 ν c ω p 2 ( ν c 2 + ω b 2 )+2j ν c ω ω 2 E x + ε 0 ω b ω p 2 ( ν c 2 + ω b 2 )+2j ν c ω ω 2 E y .
J(ω)= a 0 + a 1 jω b 0 + b 1 jω+ b 2 (jω) 2 E(ω)
b 0 J+ b 1 dJ dt + b 2 d 2 J d t 2 = a 0 E+ a 1 dE dt .
b 0 J n1 + b 1 J n J n2 2Δt + b 2 J n 2 J n1 + J n2 Δ t 2 = a 0 E n1/2 + E n3/2 2 + a 1 E n1/2 E n3/2 Δt .
J n = 4 b 2 Δt2 b 0 Δ t 2 b 1 Δt+2 b 2 J n1 + b 1 Δt2 b 2 b 1 Δt+2 b 2 J n2 + a 0 Δ t 2 +2 a 1 Δt b 1 Δt+2 b 2 E n1/2 + a 0 Δ t 2 2 a 1 Δt b 1 Δt+2 b 2 E n3/2 .
ε re (ω)=[ ε d (ω) j ε g (ω) 0 j ε g (ω) ε d (ω) 0 0 0 0 ]
ε d (ω)=1+ σ 0 +ν σ 0 /(jω) ( ν 2 + ω c 2 )+2jνω ω 2
ε g (ω)= ω b σ 0 /(jω) ( ν 2 + ω c 2 )+2jνω ω 2
σ 0 = 2 e 2 τ k B Tν π 2 ε 0 Δz ln( 2cosh μ c 2 k B T )
ε d (ω)=1+[ σ 0 /(jω)]/(ν+jω)
E 0 =[e/( ε 0 π 2 v F 2 )] 0 ε[ f d ( ε ) f d ( ε+2 μ c ) ]dε
E x inc (t)=exp[4π (t t 0 ) 2 / τ 2 ]
T analytic (f)= 2 | 2+jω[ ε d ω1] η 0 | 2 + | ε 0 ε g ω 2 η 0 | 2 | ( 2+jω[ ε d ω1] η 0 ) 2 + ( ε 0 ε g ω 2 η 0 ) 2 |
T leapfrog (f)= | FFT( E x,ob (t)) | 2 + | FFT( E y,ob (t)) | 2 | FFT( E in (t)) |

Metrics