On this page, we calculate the circular dichroism (CD) of a gammadion shaped structure, then optimize the structure dimensions to maximize the CD at 1.1um.

We consider a gammadion shaped periodic structure shown in Fig. 1. In this structure, an aluminum layer is sandwiched by silver layers and the excitation of surface plasmon leads to the enhancement of the circular dichrosim.
Figure 1 Gammadion shaped planar chiral material
The circular dichroism CD is defined by
$$ C D=\leftT_{R}T_{L}\right $$ (1)
where TR and TL is a transmittance when the right and the left circularly polarized plane wave is incident on the device, respectively. To get the transmittance for circularly polarized incidence, we have two alternatives as follows.
1. Use two plane wave sources to generate circularly polarized illumination, as described in the circular polarization page. In this case, two FDTD simulations will be required to get the CD; One for rightcircular polarization and the other for leftcircular polarization. This approach is not used in the associated example simulation file.
2. Use one plane wave source (as in the example simulation file gammadion_dichrosim.fsp). By taking advantage of the fourfold rotational symmetry of the structure, the transmittance can be obtained from a single simulation, as explained below:
The field distributions F (E or H) for circular illumination can be obtained from a single linearly polarized simulation by
$$ \begin{array}{l}{F_{R}(x, y, z)=F_{x}(x, y, z)+F_{y}(x, y, z) e^{j \pi / 2}} \\ {F_{L}(x, y, z)=F_{x}(x, y, z)+F_{y}(x, y, z) e^{j \pi / 2}}\end{array} $$ (2)
where FR (FL) is the field distribution for right (left) circularly polarized incident wave, and Fx (Fy) is the field distribution for a x (y) linearly polarized plane wave. If we assume fourfold rotational symmetry of the structure, the field distribution for ypolarized plane wave is incident on the structure, Fy, is given by that for xpolarized incident plane wave, Fx, as
$$ \begin{array}{l}{F_{x_{y}}(x, y, z)=F_{y_{y} x}\left(y_{z}x_{z}\right)} \\ {F_{y_{} y}(x, y, z)=F_{x_{} x}\left(y_{z}x_{z}\right)} \\ {F_{z_{} y}(x, y, z)=F_{z_{} x}\left(y_{z}x_{z}\right)}\end{array} $$ (3)
where FU_V (u=x, y, z; v=x, y) is the ucomponent of the field distribution for vpolarized incident wave.
Once we get the field distribution FU (U=R or L) for circularly polarized plane wave using the relation Eqs.(2) and (3) from one FDTD simulation (simulation for x or y polarized incident wave), the power traveling down in the substrate over a unit cell is given by
$$ P_{U}=\int_{S_{} \ldots \infty} \frac{1}{2} \operatorname{Re}\left(E_{U} \times H_{U}^{*}\right) \mathrm{d} \mathrm{S} $$ (4)
If we normalize the power P by incident power using script function "sourcepower", we obtain transmittance T as
$$ T_{v}=P_{v} / \text { sourcepower } $$ (5)
The simulation file gammadion_dichrosim.fsp uses a single xpolarized source. The field distribution on a plane under the gammadion structure (in the substrate) is recorded in a power monitor named "T", within the analysis group named "CD". The script in the "Analysis" tab => "Script" tab of this analysis group calculates the transmittance of circular polarizations following the way mentioned above (method 2). From the CD analysis object, you can plot the CD vs wavelength by use of the "visualizer". In the figure below, we can see a peak in the CD around 1.1um.