Incorporating the Havriliak-Negami dielectric model in the FD-TD method

We derive and analyze an efficient algorithm to incorporate the anomalously dispersive Havriliak-Negami dielectric model of induced polarization in the Finite-difference time-domain (FD-TD) method. Our algorithm implements this dielectric model, which in the time-domain involves fractional derivatives and fractional differential operators, with a preset error over the desired computational time interval [0,T_comp] and correctly takes into account the singularity at t=0⁺ of the corresponding time-domain dielectric susceptibility. The overall algorithm is shown to be second-order accurate in space and time, and to obey the standard FD-TD stability condition. Numerical experiments confirm our analysis.