光子晶体介质柱的讨论

矩形介质柱的$k$空间相对介电常数计算：

$$\alpha (G)=\{\begin{array}{l}{\alpha }_{\mathrm{A}}f+{\alpha }_{\mathrm{B}}(1-f),G=0\\ \left({\alpha }_{\mathrm{A}}-{\alpha }_{\mathrm{B}}\right)P(G),G\ne 0\end{array}$$
其中，$f$是填充率，等于介质柱横截面积/二维正格子原胞面积，${\alpha }_{\mathrm{A}}$是背景的相对介电常数，${\alpha }_{\mathrm{B}}$是介质的相对介电常数，$G$在这里是指倒格矢，但在编程中不是完全套倒格矢进去，$P(G)$在这里叫结构因子，对于矩形介质柱，它指代如下：
$$P(\mathbf{G})=f\frac{\sin \left({\tilde{G}}_x{l}_x\right)}{{\tilde{G}}_x{l}_x}\frac{\sin \left({\tilde{G}}_y{l}_y\right)}{{\tilde{G}}_y{l}_y}$$
其中，$2l_x$和$2l_y$分别是矩形柱子的长度和宽度，$f$是填充率，其中还有：
$$\{\begin{array}{l}{\tilde{G}}_x=G_x\cos \theta +G_y\sin \theta \\ {\tilde{G}}_y=-G_x\sin \theta +G_y\cos \theta \end{array}$$
其中$\theta$是我的另一篇博文里说的输入的参数之一，代表的意义如下：

这个是正空间里的图像。

正六边形介质柱的$k$空间相对介电常数计算：

同上面一样，有：
$$\alpha (G)=\{\begin{array}{l}{\alpha }_{\mathrm{A}}f+{\alpha }_{\mathrm{B}}(1-f),G=0\\ \left({\alpha }_{\mathrm{A}}-{\alpha }_{\mathrm{B}}\right)P(G),G\ne 0\end{array}$$
其中各参数和上面意义相同，结构因子为：

$$P(G)=\{\begin{array}{ll}\frac{2f}{3G_{x}b}\left[\sin \left({\tilde{G}}_{x}b\right)+\frac{1-\cos \left({\tilde{G}}_{x}b\right)}{{\tilde{G}}_{x}b}\right],& {\tilde{G}}_x\ne 0,{\tilde{G}}_y=0\\ \frac{f}{3G_{x}b}\left[\frac{1-\cos \left(2{\tilde{G}}_{x}b\right)}{2{\tilde{G}}_{x}b}+\sin \left(2{\tilde{G}}_{x}b\right)\right],& \left|{\tilde{G}}_x\right|=\left|{\tilde{G}}_y\right|/\sqrt{3},{\tilde{G}}_y\ne 0,\\ \frac{f}{G_{y}b}\left[\frac{\cos \left(\frac{{\tilde{G}}_{y}b}{\sqrt{3}}-{\tilde{G}}_{x}b\right)-\cos \left(\frac{2{\tilde{G}}_{y}b}{\sqrt{3}}\right)}{\left({\tilde{G}}_y+\sqrt{3}{\tilde{G}}_x\right)b}+\frac{\cos \left(\frac{{\tilde{G}}_{y}b}{\sqrt{3}}+{\tilde{G}}_{x}b\right)-\cos \left(\frac{2{\tilde{G}}_{y}b}{\sqrt{3}}\right)}{\left({\tilde{G}}_y-\sqrt{3}{\tilde{G}}_x\right)b}\right],\left|{\tilde{G}}_x\right|\ne \left|{\tilde{G}}_y\right|/\sqrt{3},{\tilde{G}}_y\ne 0& \end{array}$$
其中，b是六边形的六个边宽度。而且也有：
$$\{\begin{array}{l}{\tilde{G}}_x=G_x\cos \theta +G_y\sin \theta \\ {\tilde{G}}_y=-G_x\sin \theta +G_y\cos \theta \end{array}$$
其中$\theta$是我的另一篇博文里说的输入的参数之一，代表的意义如下：

对于填充率：
$$f=4l_x{l}_y/a^2\text{ and }f=3\sqrt{3}{b}^2/2a^2$$
参考文献：
Hui Y , Guo-Hong Y , Yong-Jun C . Effects of rotating noncircular scatterers on spin-wave band gaps of two-dimensional magnonic crystals[J]. Chinese Physics B, 2014(9):423-427.

正三角形介质柱的$k$空间相对介电常数计算：

${\epsilon }^{-1}\left(\mathbf{G}-{\mathbf{G}}^{\prime }\right)$ represents the inverse of the matrix $\epsilon \left(\mathbf{G}-{\mathbf{G}}^{\prime }\right).$
同时也有：
$$\epsilon (\mathbf{G})=\{\begin{array}{ll}f\cdot {\epsilon }_{\mathrm{A}}+(1-f){\epsilon }_{\mathrm{B}}& \mathbf{G}=0\\ \left({\epsilon }_{\mathrm{A}}-{\epsilon }_{\mathrm{B}}\right)I(\mathbf{G})& \mathbf{G}\ne 0\end{array}$$
其中，其中，$f$是填充率，等于介质柱横截面积/二维正格子原胞面积，${\epsilon }_{\mathrm{A}}$是背景的相对介电常数，${\epsilon }_{\mathrm{B}}$是介质的相对介电常数，其中$I(\mathbf{G})$等于：
$$I(\mathbf{G})=\frac{1}{s}{\int }_{s_{\mathrm{r}}}{\mathrm{e}}^{-\mathrm{i}G\cdot \mathrm{r}}\mathrm{d}\mathbf{r}$$
这是一个普遍适用的积分。这个积分是很重要的，它的积分区域是在一个单元中棒的横截面积上积分，$s$是代表正格子原胞面积。它是面积积分。
对于正三角形：

其中，$l$是三角形每条边的长度。$G_x$和$G_y$是倒格矢$G$的$x$和$y$分量，而且有：

$$\left(\begin{array}{l}{\hat{G}}_x\\ {\hat{G}}_y\end{array}\right)=\left(\begin{array}{cc}\cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{array}\right)\left(\begin{array}{l}{G}_x\\ G_y\end{array}\right)$$
$\theta$代表的意义如下：

参考文献：
Kuang W , Hou Z , Liu Y , et al. The bandgap of a photonic crystal with triangular dielectric rods in a honeycomb lattice[J]. Journal of Optics A Pure & Applied Optics, 2005, 7(10):525.

椭圆形介质柱的$k$空间相对介电常数计算：

同样有：
$$\epsilon (\mathbf{G})=\{\begin{array}{ll}f\cdot {\epsilon }_{\mathrm{A}}+(1-f){\epsilon }_{\mathrm{B}}& \mathbf{G}=0\\ \left({\epsilon }_{\mathrm{A}}-{\epsilon }_{\mathrm{B}}\right)P(\mathbf{G})& \mathbf{G}\ne 0\end{array}$$
各参数代表的含义和上述相同。
$$\begin{array}{rl}& P(\vec{G})=\frac{2F_d{J}_1\left(\sqrt{{\left(L_x{G}_x^{\prime }\right)}^2+{\left(L_y{G}_y^{\prime }\right)}^2}\right)}{\sqrt{{\left(L_x{G}_x^{\prime }\right)}^2+{\left(L_y{G}_y^{\prime }\right)}^2}}\\ & \text{ where }\\ & \left(\begin{array}{l}{G}_x^{\prime }\\ G_y^{\prime }\end{array}\right)=\left(\begin{array}{cc}\cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{array}\right)\left(\begin{array}{l}{G}_x\\ G_y\end{array}\right)\end{array}$$
其中，其中Lx是椭圆沿x方向的半轴，Ly是椭圆沿y方向的半轴，$\theta$是椭圆的夹角。$J_1$是一阶贝塞尔函数，$F_d=$ $\pi L_x{L}_y/S_{\text{cell }}$，$S_{\text{cell }}$是正格子原胞面积。

参考文献：
Feng Y , Wu F , Zhong H , et al. Defect modes created by an elliptic defect in two-dimensional triangular photonic crystals[J]. Solid State Communications, 2007, 142(4):223-227.