I am trying to reproduce figures 4f and 4i from the open-access paper: " Launching of hyperbolic phonon-polaritons in h-BN slabs by resonant metal plasmonic antennas" (link below). This paper uses metallic rod antennas to launch polaritons in hBN slab (which is an anisotropic 2D material that supports volume polaritons).
Although the results are close, my simulation shows some extra field normal to the antenna border (indicated in the figure). I have tried to change the mesh resolution, the source type from TFSF to a Gaussian beam, and even the PML boundaries from standard to stabilizes. All these attempts resulted in the same near-field profile.
It is worth saying that I would not expect this extra field since the experimental and simulated data from the paper does not present this artifact.
I would be greatly thanked if anyone could help me with this issue.
Obs: I am using version 8.23 of FDTD solutions.
Antenna_on_hBN_fig_4.fsp (303.2 KB)