删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

含石墨烯分界面有耗分层介质的传播矩阵

本站小编 Free考研考试/2021-12-29

摘要:给出一种适用于含导电界面的有耗分层介质的传播矩阵方法. 利用相位匹配原理给出斜入射时有耗介质波矢量的实部和虚部, 二者方向不同使得在介质中传播非均匀平面波. 根据边界条件, 推导了跨越石墨烯界面的传播矩阵, 以及“无限薄”石墨烯层的反透射系数解析式. 最终将传播矩阵方法推广应用于含石墨烯界面的有耗分层介质情形, 可用于快速解析分析分层介质与导电界面复合结构的反透射和电波传播特性.
关键词: 石墨烯/
有耗分层介质/
传播矩阵

English Abstract


--> --> -->
2004年, 世界上最薄的半导体材料—石墨烯, 由Novoselov和Geim的团队采用微机械剥离法制备成功[1]. 石墨烯是由单层碳原子构成的二维晶体结构, 具有高电子迁移率、低电阻率、光学透明性、高灵活性及结构稳定等特性[2], 这使得其在光学和光电传感等方面具有广泛的应用前景[3-6], 同时也被认为非常适合于太赫兹光学设备[7-10]. 石墨烯非常薄, 厚度只有约0.34 nm, 这导致其光吸收率很低, 应用受到限制. 为此, 常采用周期掺杂、表面等离子体、光学微腔等方法来增强光学吸收[11-16], 其中利用光子晶体增强石墨烯光吸收是一种常用方法. 2012年, Liu等[17]利用石墨烯和间隔层在一维光子晶体(one-dimensional photonic crystal, 1DPC)表面构造表面缺陷, 导致光的局域化, 使可见光光吸收被增强约4倍. 2013年, Peres和Bludov[18]利用1DPC使石墨烯的THz吸收增强了超过3倍. 2014年, Xie等[19]分析了1DPC表面石墨烯在应力赝磁场作用下的THz吸收. 2018年, Zhang等[20]设计了一种高对比光栅耦合波导结构, 提高了石墨烯在近红外和可见光频段的光吸收.
1DPC本质上是一个周期分层介质, 而单层的石墨烯可以看做是“无限薄”的[21], 因此可以把石墨烯和1DPC构造的结构看成一个含石墨烯分界面的分层介质. 传播矩阵方法[22]是一种分析分层介质中光或电磁波传播特性的解析理论方法, 可以直接计算分层介质对平面电磁波的反透射等. 根据电磁理论, 当分层介质中某层为有耗介质时, 该层中的波矢量为复数, 而在斜入射情形下, 考虑相位匹配, 该层中波矢量的实部和虚部并不平行, 即相位和振幅的传播方向不同, 此时该层中传播的应为非均匀平面波. 本文首先根据相位匹配条件, 推导了有耗分层介质中复数波矢量的实部和虚部, 给出传播矩阵中所需要的纵向波矢量计算式; 然后考虑石墨烯表面电导率, 根据边界条件推导了跨越石墨烯分界面的传播矩阵; 最终, 使得传播矩阵方法可用于解析分析任意入射角情形下含导电界面有耗分层介质的反透射和电波传播特性. 本文中时谐因子取$\exp \left( { - {\rm{j}}\omega t} \right)$.
2
2.1.有耗分层介质中的复数波矢量
-->设平面波从真空入射到$n$层分层介质, 如图1所示, 入射角为$\theta $. 考虑一般情形, 设各层介质的相对介电系数和磁导系数${\varepsilon _{{\rm{r}}l}}, {\mu _{{\rm{r}}l}}\;(l = 1, 2, \cdots, n)$皆为复数,
图 1 分层介质
Figure1. Stratified medium.

${\varepsilon _{{\rm{r}}l}} = {\varepsilon '_{{\rm{r}}l}} + {\rm{j}}{\varepsilon ''_{{\rm{r}}l}},\;{\mu _{{\rm{r}}l}} = {\mu '_{{\rm{r}}l}} + {\rm{j}}{\mu ''_{{\rm{r}}l}}.$
此时, 介质有耗, 介质层中的波矢量${{{k}}_l}$也为复数. 下面来具体分析介质层中的复数波矢量${{{k}}_l}$.
各分层中的复数波矢量${{{k}}_l}$可以写为
${{{k}}_l} = {{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}},$
其中${{{k}}_{{\rm{R}}l}}$${{{k}}_{{\rm{I}}l}}$分别表示${{{k}}_l}$的实部和虚部. 如图1所示, 入射面为$xoz$平面, 故有
$\begin{split}& {{{k}}_{{\rm{R}}l}} = \hat z{k_{{\rm{R}}lz}} + \hat x{k_{{\rm{R}}lx}}, \\& {{{k}}_{{\rm{I}}l}} = \hat z{k_{{\rm{I}}lz}} + \hat x{k_{{\rm{I}}lx}}.\end{split} $
根据界面处的相位匹配(波矢量切向分量连续)条件可得
${k_{{\rm{R}}lx}} = {k_{0x}} = {k_0}\sin \theta,{\rm{ }}{k_{{\rm{I}}lx}} = 0.$
${{{k}}_l}$的虚部${{{k}}_{{\rm{I}}l}}$只有$z$分量, 实部${{{k}}_{{\rm{R}}l}}$$x$分量在每个区域都相同, 都等于区域0中入射波矢量${{{k}}_0}$$x$分量${k_{0 x}}$. 由此可以看到, 平面波斜入射有耗分层介质时, 介质中波矢量的实部和虚部方向不平行, 虚部只有纵向分量. 由于波矢量实部代表相位传播方向, 虚部代表振幅传播方向, 可知此时平面波是非均匀的, 波矢量为
$\begin{split} {{{k}}_l} ={}& \hat z\left( {{k_{{\rm{R}}lz}} + {\rm{j}}{k_{{\rm{I}}l}}} \right) + \hat x{k_{{\rm{R}}lx}} \\ ={}& \hat z\left( {{k_{{\rm{R}}lz}} + {\rm{j}}{k_{{\rm{I}}l}}} \right) + \hat x{k_{0x}}. \end{split} $
由电磁波理论可知在有耗介质中[22]
$\begin{split} {{{k}}_l} \cdot {{{k}}_l} ={}& \omega {\varepsilon _0}{\mu _0}{\varepsilon _{{\rm{r}}l}}{\mu _{{\rm{r}}l}} \\ ={}& k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}} + {\rm{j}}{{\varepsilon ''}_{{\rm{r}}l}}} \right)\left( {{{\mu '}_{{\rm{r}}l}} + {\rm{j}}{{\mu ''}_{{\rm{r}}l}}} \right). \end{split}$
而由(2)和(5)式可得
$\begin{split} {{{k}}_l} \cdot {{{k}}_l}\, &= \left( {{{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}}} \right) \cdot \left( {{{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}}} \right) \\ &= {{{k}}_{{\rm{R}}l}} \cdot {{{k}}_{{\rm{R}}l}} - {{{k}}_{{\rm{I}}l}} \cdot {{{k}}_{{\rm{I}}l}} + 2{\rm{j}}\left( {{{{k}}_{{\rm{R}}l}} \cdot {{{k}}_{{\rm{I}}l}}} \right) \\ & = k_{{\rm{R}}lz}^2 - k_{{\rm{I}}l}^2 + 2{\rm{j}}\left( {{k_{{\rm{R}}lz}} \cdot {k_{{\rm{I}}l}}} \right) + k_{0x}^2.\end{split} $
则根据(6)和(7)式得
$\begin{split} & k_{{\rm{R}}lz}^2 - k_{{\rm{I}}l}^2 = k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}} - {{\varepsilon ''}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}}} \right) - k_0^2{\sin ^2}\theta, \\ & 2{k_{{\rm{R}}lz}} \cdot {k_{{\rm{I}}l}} = k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}}{\mu ''}_{\rm rl} + {{\varepsilon ''}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}}} \right).\end{split} $

$\begin{split} & AA = ({{\varepsilon '}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}} - {{\varepsilon ''}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}}) - {\sin ^2}\theta,\\ & BB = \left( {{{\varepsilon '}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}} + {{\varepsilon ''}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}}} \right), \\ & CC = {{BB} / {AA}}.\end{split} $
求解(8)式得到
$\left\{ \begin{split}& {k_{{\rm{I}} l}} = {k_{{\rm{R}}lz}} = {k_0}\sqrt {\frac{{BB}}{2}},\quad\quad\quad\quad\quad\quad\quad\;\; AA = 0; \\& {k_{{\rm{I}} l}} = {k_0}\sqrt {AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} - 1}}{2}},\\ & {k_{{\rm{R}}lz}} = {k_0}\sqrt {AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} + 1}}{2}},\quad\;\;\;\; AA > 0; \\ & {k_{{\rm{I}} l}} = {k_0}\sqrt { - AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} + 1}}{2}},\\ & {k_{{\rm{R}}lz}} = {k_0}\sqrt { - AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} - 1}}{2}}, \quad\; AA < 0.\end{split} \right.$
综上所述, 考虑到相位匹配, 有耗分层介质中平面波的复数波矢量可通过(10)式计算. 计算得到波矢量虚部${{{k}}_{{\rm{I}}l}}$沿$z$方向, 所以波在介质中沿深度方向衰减. 波矢量实部${{{k}}_{{\rm{R}}l}}$的方向代表波传播方向, 介质中的折射角可以依此计算.
2
2.2.有耗分层介质的传播矩阵
-->对整个分层介质, 在图1中从入射点${M_0}$到透射点${N_{\rm{t}}}$的传播用前向传播矩阵表示为
$\left[ {\begin{array}{*{20}{c}} 0 \\ {T{B_0}\exp ( - {\rm{j}}{k_{{\rm{t}}z}}{z_n})} \end{array}} \right] = {{{V}}_{{\rm{t}}0}}\left[ {\begin{array}{*{20}{c}} {R{B_0}\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {{B_0}\exp ( - {\rm{j}}{k_{0z}}{z_0})} \end{array}} \right],$
其中${B_0}$${M_0}$点处下行波幅值. 整个分层介质总的前向传播矩阵为[22]
${{{V}}_{{\rm{t0}}}} = {{{U}}_n} \cdot {{{V}}_n} \cdots {{{U}}_1} \cdot {{{V}}_1} \cdot {{{U}}_0},$
其中${{{U}}_l}$${{{V}}_l}$分别是跨越$l-l + 1$层界面和第$l$层中的传播矩阵. 对有耗介质, 传播矩阵计算时${\varepsilon _l}, {\mu _l}$取复数, 另外要特别注意${k_{lz}}$, 应该按照前述方法计算.
2
2.3.跨越石墨烯界面的传播矩阵
-->单层的石墨烯可以看做是“无限薄”的[21], 石墨烯薄层两个侧面的表面电导率为${\sigma _{\rm{g}}}\left( {\omega, {\mu _{\rm c}}, \varGamma, T} \right)$, 其中$\omega $, ${\mu _{\rm c}}$, $\varGamma $$T$分别是角频率、化学势、散射率和室温. 由Kobu公式得到石墨烯表面电导率表达式为${\sigma _{\rm{g}}} = {\sigma _{{\rm{intra}}}} + {\sigma _{{\rm{inter}}}}$, ${\sigma _{{\rm{intra}}}}$取决于能带内的能量, 而${\sigma _{{\rm{inter}}}}$取决于能带间的能量. 因为能带内的能量占能量主要部分, 所以只考虑${\sigma _{{\rm{intra}}}}$这一项, 而忽略${\sigma _{{\rm{inter}}}}$项.
${\sigma _{{\rm{intra}}}}$可用Drude表达式表示为
${\sigma _{{\rm{intra}}}}(\omega,{\mu _{\rm{c}}},\varGamma,T) = \frac{{{\sigma _0}}}{{1 - {\rm{j}}\omega \tau }},$
式中
${\sigma _0} = \frac{{{e^2}\tau {k_{\rm{B}}}T}}{{{\text{π}}{\hbar ^2}}}\left\{ {\frac{{{\mu _{\rm{c}}}}}{{{k_{\rm{B}}}T}} + 2\ln \left[ {\exp \left( { - \frac{{{\mu _{\rm{c}}}}}{{{k_{\rm{B}}}T}}} \right) + 1} \right]} \right\},$
为直流电导率, 其中e是电子电量, ${k_{\rm{B}}}$是玻尔兹曼常数, $\hbar $是约化普朗克常量; $\tau = {1 / {\left( {2\varGamma } \right)}}$是电子弛豫时间.
设在分层介质区域$l$$l + 1$的分界面$z = {z_l}$处有一层“无限薄”石墨烯层, 如图1所示, 根据电磁理论, 在介质分界面两侧的电磁场满足边界条件, 所以在分界面$z = {z_l}$处有
$\begin{split} & \hat z \times \left( {{{{E}}_l} - {{{E}}_{l + 1}}} \right) = - {{J}}_{\rm{m}}^{\rm{s}}, \\& \hat z \times \left( {{{{H}}_l} - {{{H}}_{l + 1}}} \right) = {{J}}_{\rm{e}}^{\rm{s}},\end{split} $
其中${{J}}_{\rm{e}}^{\rm{s}}$${{J}}_{\rm{m}}^{\rm{s}}$表示表面电流和磁流. 由于该界面为“无限薄”石墨烯层, 则此处[23]
$\begin{split}& {{J}}_{\rm{m}}^{\rm{s}} = 0, \\& J_{{\rm{e}},x}^{\rm{s}} = {\sigma _{\rm{g}}}{E_x},~~ J_{{\rm{e}},y}^{\rm{s}} = {\sigma _{\rm{g}}}{E_y}.\end{split} $
联合(15)式可得
${E_{l\alpha }}\left( {z = {z_l}^ + } \right) = {E_{\left( {l + 1} \right)\alpha }}\left( {z = {z_l}^ - } \right),{\rm{ }}\alpha = x,y,$
${H_{\left( {l + 1} \right)y}}\left( {z = {z_l}^ - } \right) - {H_{ly}}\left( {z = {z_l}^ + } \right) = {\sigma _{\rm{g}}}{E_x}\left( {z = {z_l}} \right),$
${H_{( {l + 1} )x}}( {z = {z_l}^ - } ) - {H_{lx}}( {z = {z_l}^ + } ) = - {\sigma _{\rm{g}}}{E_y}( {z = {z_l}} ).$
首先考虑TE波情形, 区域l中的总场可写为[22]
$\begin{split}&{E_{ly}} = [ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) + {B_l}\exp ( { - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ), \\& {H_{lx}} = \frac{{ - {k_{lz}}}}{{\omega {\mu _l}}}[ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) \!-\! {B_l}\exp ( {\! - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ), \\& {H_{lz}} = \frac{{{k_x}}}{{\omega {\mu _l}}}[ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) \!+\! {B_l}\exp ( {\! - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ),\end{split} $
式中${A_l}$${B_l}$分别代表区域$l$中的上行波和下行波幅值. 在区域$l$$l + 1$的分界面$z = {z_l}$处, 根据边界条件(17)和(19)式可得
$\begin{split}& {A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) - {B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ =\, & \left( {p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \\ &+ \left( { - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \\[-16pt]\end{split} $
$\begin{split}& \left( {1 - s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ &- \left( {1 + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ = \; &p_{\left( {l + 1} \right)l}^{{\rm{TE}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) - {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right)} \right], \end{split} $
其中
$s_l^{{\rm{TE}}} = \frac{{\omega {\mu _l}{\sigma _{\rm{g}}}}}{{{k_{lz}}}}.$
求解并整理(21)和(22)式可得
$\left[\!\!\!{\begin{array}{*{20}{c}} {{A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}{z_l}} \right)} \\ {{B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}{z_l}} \right)} \end{array}}\!\!\!\right] \!=\! {{U}}_{\rm{g}}^{{\rm{TE}}} \cdot \left[\!\!\!{\begin{array}{*{20}{c}} {{A_l}\exp \left( {{\rm{j}}{k_{lz}}{z_l}} \right)} \\ {{B_l}\exp \left( { - {\rm{j}}{k_{lz}}{z_l}} \right)} \end{array}}\!\!\!\right],$
其中
${{U}}_{\rm{g}}^{{\rm{TE}}} = \frac{1}{2}\left[\!\!\!{\begin{array}{*{20}{c}} {1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{l + 1}^{{\rm{TE}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{l + 1}^{{\rm{TE}}}} \\ {1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - s_{l + 1}^{{\rm{TE}}}}&{1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - s_{l + 1}^{{\rm{TE}}}} \end{array}}\!\!\! \right]$
是穿越$z = {z_l}$处石墨烯分界面从${M_l}$点到${N_{l + 1}}$点的TE波前向传播矩阵.
对TM波情形, 区域$l$中的总场可写为
$\begin{split} &{H_{ly}} = \left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) + {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right), \\ &{E_{lx}} = \frac{{{k_{lz}}}}{{\omega {\varepsilon _l}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \!- \!{B_l}\exp \left( {\! - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right), \\& {E_{lz}} = \frac{{ - {k_x}}}{{\omega {\varepsilon _l}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \!+\! {B_l}\exp \left( {\! - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right).\end{split} $
在区域$l$$l + 1$的分界面$z = {z_l}$处, 根据边界条件(17)和(18)式可得
$\begin{split} &{A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) + {B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ =\, & \left( {1 + s_l^{{\rm{TM}}}} \right){A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \\ & + \left( {1 - s_l^{{\rm{TM}}}} \right){B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \\[-16pt]\end{split} $
$\begin{split} & \left( {1 - s_{\left( {l + 1} \right)}^{{\rm{TM}}}} \right){A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ &+ \left( {1 + s_{\left( {l + 1} \right)}^{{\rm{TM}}}} \right){B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\=\,& {A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) + {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \end{split}$
其中
$s_l^{{\rm{TM}}} = \frac{{{k_{lz}}{\sigma _{\rm{g}}}}}{{\omega {\varepsilon _l}}}.$
求解并整理(27)和(28)式可得
$\left[\!\!\!{\begin{array}{*{20}{c}} {{A_{( {l + 1} )}}\exp ( {{\rm{j}}{k_{( {l + 1} )z}}z} )} \\ {{B_{( {l + 1} )}}\exp ( { - {\rm{j}}{k_{( {l + 1} )z}}z} )}\end{array}} \!\!\!\right] = {{U}}_{\rm{g}}^{{\rm{TM}}} \cdot \left[\!\!\! {\begin{array}{*{20}{c}} {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} )} \\ {{B_l}\exp ( { - {\rm{j}}{k_{lz}}z} )}\end{array}}\!\!\!\right],$
其中
${{U}}_{\rm{g}}^{{\rm{TM}}} = \frac{1}{2}\left[\!\!\!{\begin{array}{*{20}{c}} {1 + p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + s_l^{{\rm{TM}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - s_l^{{\rm{TM}}}} \\ {1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + s_l^{{\rm{TM}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - s_l^{{\rm{TM}}}} \end{array}}\!\!\!\right]$
是穿越$z = {z_l}$处石墨烯分界面从${M_l}$点到${N_{l + 1}}$点的TM波前向传播矩阵.
实际上, 上述方法中所涉及的石墨烯界面可拓展到表面电导率为${\sigma _{\rm{g}}}$的导电界面.
综上, 若分层介质中某界面含“无限薄”导电层, 则可以通过将总传播矩阵${{{V}}_{{\rm{t0}}}}$((12)式)中该界面处的跨界面传播矩阵${{{U}}_{l = 0 - n}}$替换为${{U}}_{\rm{g}}^{{\rm{TE/TM}}}$((25)和(31)式), 从而得到含导电界面的分层介质传播矩阵. 相反, 若令${\sigma _{\rm{g}}} = 0$, 则${{U}}_{\rm{g}}^{{\rm{TE/TM}}}$又可以退化为${{{U}}_{l = 0 - n}}$, 即为一般电介质分界面.
注意到, 由于石墨烯界面的导电特性, 此处TE和TM波的边界条件不对偶, 因而其跨界面传播矩阵在TE和TM波情形下并不对偶.
2
2.4.“无限薄”石墨烯层的反透射系数
-->设“无限薄”石墨烯上、下半空间分别为均匀介质1和2, 在(24)和(30)式中令$l = 1$${A_2} = 0$, 结合(25)式、(31)式和反透射系数定义可以得到“无限薄”石墨烯界面的反透射系数
$\begin{split} {R^{{\rm{TE}}}} = & - \frac{{s_1^{{\rm{TE}}} - s_2^{{\rm{TE}}} + s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{s_1^{{\rm{TE}}} + s_2^{{\rm{TE}}} + s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}} \\ = & - \frac{{{\mu _1}{k_{2z}} - {\mu _2}{k_{1z}} + {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{{\mu _1}{k_{2z}} + {\mu _2}{k_{1z}} + {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}, \\ {\rm{ }}{T^{{\rm{TE}}}} =\; & \frac{{s_1^{{\rm{TE}}} + s_2^{{\rm{TE}}} - s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{2s_1^{{\rm{TE}}}}} \\ &+ {R^{{\rm{TE}}}} \cdot \left( {\frac{{s_1^{{\rm{TE}}} - s_2^{{\rm{TE}}} - s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{2s_1^{{\rm{TE}}}}}} \right)\\= \;& \frac{{{\mu _1}{k_{2z}} + {\mu _2}{k_{1z}} - {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{2{\mu _1}{k_{2z}}}} \\ & + {R^{{\rm{TE}}}} \cdot \left( {\frac{{{\mu _1}{k_{2z}} - {\mu _2}{k_{1z}} - {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{2{\mu _1}{k_{2z}}}}} \right),\end{split} $
$\begin{split} {R^{{\rm{TM}}}} =\, & - \frac{{s_2^{{\rm{TM}}} - s_1^{{\rm{TM}}} - s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{s_2^{{\rm{TM}}} + s_1^{{\rm{TM}}} + s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}\\ =\, & - \frac{{{\varepsilon _1}{k_{2z}} - {\varepsilon _2}{k_{1z}} - {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{{\varepsilon _1}{k_{2z}} + {\varepsilon _2}{k_{1z}} + {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}, \\ {T^{{\rm{TM}}}} =\, & \frac{{s_2^{{\rm{TM}}} + s_1^{{\rm{TM}}} - s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{2s_2^{{\rm{TM}}}}} \\ & + {R^{{\rm{TM}}}} \cdot \left( {\frac{{s_2^{{\rm{TM}}} - s_1^{{\rm{TM}}} + s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{2s_2^{{\rm{TM}}}}}} \right) \\ =\, & \frac{{{\varepsilon _1}{k_{2z}} + {\varepsilon _2}{k_{1z}} - {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{2{\varepsilon _1}{k_{2z}}}} \\ & + {R^{{\rm{TM}}}} \cdot \left( {\frac{{{\varepsilon _1}{k_{2z}}\omega - {\varepsilon _2}{k_{1z}}\omega + {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{2{\varepsilon _1}{k_{2z}}\omega }}} \right).\end{split} $
当垂直入射时, ${k_{1 z}} = {k_1} = {k_0}\sqrt {{\varepsilon _{{\rm{r1}}}}{\mu _{{\rm{r1}}}}}, {k_{2 z}} = {k_2} = $${k_0}\sqrt {{\varepsilon _{{\rm{r2}}}}{\mu _{{\rm{r2}}}}} $, 由(32)和(33)式可得“无限薄”石墨烯界面的垂直反透射系数
$\begin{split}& R = \frac{{{\eta _2} - {\eta _1} - {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}}{{{\eta _2} + {\eta _1} + {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}},\\ & {\eta _n} = \sqrt {\frac{{{\mu _n}}}{{{\varepsilon _n}}}}\;\;\left( {n = 1,2} \right), \\& T = \frac{{2{\eta _2}}}{{{\eta _2} + {\eta _1} + {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}} = 1 + R.\end{split} $
结果和文献[23]一致. 注意到, 在(32), (33)和(34)式中若令${\sigma _{\rm{g}}} = 0$, 则可以自然过渡为一般介质分界面的Fresnel公式.
例1 石墨烯单层的屏蔽效率
设各向同性石墨烯层的参数为T = 300 K, $ {\mu _{\rm{c}}} = 0, \varGamma = 0.11$ MeV, 应用传播矩阵计算其反透射系数, 进而得到屏蔽效率(shielding effectiveness, SE)[21]随频率变化曲线如图2(b)所示(图2(a)为该石墨烯层电导率实虚部与频率的关系曲线). SE的计算结果和文献[21]结果相符, 同时由图2可见, 石墨烯层对TM模的SE低于对TE模的.
图 2 石墨烯电导率和SE (a) 电导率实虚部; (b) SE
Figure2. Graphene conductivity and shielding effectiveness: (a) Real and imaginary parts of conductivity; (b) SE.

以下算例中石墨烯层的参数取
$T = 300\;{\rm{ K}},{\mu _{\rm{c}}} = 0.5\;{\rm{ eV}},\tau = 0.5\;{\rm{ ps}}.$
例2 空气中石墨烯单层的反透射
设平面电磁波入射真空中一单层石墨烯层, 根据(32)和(33)式可计算得到该石墨烯层的反透射系数, 计算结果如图3所示. 图3(a)图3(b)分别是TE和TM波的反透射系数模值, 其中实线、短划线、点线分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时的反射系数, 方形、圆形、星形分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时的透射系数. 可以看到, 随着频率增大, 该石墨烯层的反射系数减小而透射系数增大; 随着入射角增大, TE波反射增强透射减弱, 而TM波相反.
图 3 石墨烯单层的反射和透射系数模值 (a) TE模; (b) TM模
Figure3. Modulus of reflective and transmittance coefficients of a graphene sheet:(a) TE mode; (b) TM mode.

例3 涂覆石墨烯层的半空间有耗介质(碲化镉CdTe)界面
设在碲化镉(CdTe)表面上涂覆着一层石墨烯, CdTe的相对介电系数为[24]
${\varepsilon _{\rm{r}}} = {\varepsilon _\infty } + {\varepsilon _\infty }\left( {\frac{{\omega _{{\rm{LO}}}^2 - \omega _{{\rm{TO}}}^2}}{{\omega _{{\rm{TO}}}^2 - {\omega ^2} - {\rm{j}}\gamma \omega }}} \right),$
其中${\varepsilon _\infty }$, ${\omega _{{\rm{LO}}}}$, ${\omega _{{\rm{TO}}}}$, $\gamma $分别是CdTe的高频极限相对介电常数、纵波光学声子角频率、横波光学声子角频率、碰撞角频率. 本算例取文献[24]给出的测量值
$\begin{split} &{\varepsilon _\infty } = 6.93,\;\;{\omega _{{\rm{TO}}}} = 2{\text{π}} \times 4.25\;{\rm{ THz}},\\ &{\omega _{{\rm{LO}}}} = 2{\text{π}} \times 5.01\;{\rm{ THz}},\;\;\gamma = 2{\text{π}} \times 0.125\;{\rm{ THz}}.\end{split}$
计算得到CdTe的禁带为
${f_{\rm{L}}} < f < {f_{\rm{H}}},~{f_{\rm{L}}} = 4.2548\;{\rm{ THz}},~ {f_{\rm{H}}} = 5.0044\;{\rm{ THz}}.$
根据(32)和(33)式可计算该CdTe界面的反透射系数, 计算结果如图4图5所示. 图4图5分别是TE波和TM波情形, 图(a)和(b)分别是反射、透射系数模值, 各图中实线、短划线、点线分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时不含石墨烯涂层时的情形, 方形、圆形、星形分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时含石墨烯涂层时的情形. 由图4图5可见, 石墨烯的高导电特性使得其涂覆CdTe表面后, 反射增强而透射减弱, 并且在禁带以下的低频段尤为明显.
图 4 含石墨烯涂层CdTe半空间的反透射系数模值(TE模) (a) 反射系数; (b) 透射系数
Figure4. Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TE mode): (a) Reflective coefficient; (b) transmittance coefficient.

图 5 含石墨烯涂层CdTe半空间的反透射系数模值(TM模) (a) 反射系数; (b) 透射系数
Figure5. Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TM mode): (a) Reflective coefficient; (b) transmittance coefficient.

图6图7是该CdTe半空间界面TE和TM入射波的反透射光场分布, 图(a)和(b)分别是界面不含和含有石墨烯涂层时的情形, 从图6图7可见, 石墨烯涂层不改变反透射角, 但增强了表面反射, 起到了一定屏蔽作用.
图 6 CdTe半空间的反透射光场(TE模) (a) 无石墨烯涂层; (b) 含石墨烯涂层
Figure6. Optical field of reflection and transmission coefficients of CdTe half-space (TE mode): (a) Without graphene coating; (b) with graphene coating.

图 7 CdTe半空间的反透射光场(TM模) (a) 无石墨烯涂层; (b) 含石墨烯涂层
Figure7. Optical field of reflection and transmission coefficients of CdTe half-space (TM mode): (a) Without graphene coating; (b) with graphene coating.

例4 涂覆石墨烯层的Si/SiO2周期结构型一维光子晶体
有一含SiO2缺陷层的Si/SiO2周期结构型1DPC, 如图8所示. 该周期结构周期数$N = 8$, 介质a (Si)和介质b (SiO2)参数为
图 8 Si/SiO2周期结构型1DPC
Figure8. Si/SiO2 1DPC with periodic structure.

$\begin{split}& {\varepsilon _{{\rm{ra}}}} = {{\varepsilon }'_{{\rm{ra}}}} + {\rm{j}}{\varepsilon ''_{{\rm{ra}}}},\;\; {{\varepsilon }'_{{\rm{ra}}}} = 12,\;\;{{\varepsilon }''_{{\rm{ra}}}} = 0, \\& {\varepsilon _{{\rm{rb}}}} = {{\varepsilon }'_{{\rm{rb}}}} + {\rm{j}}{{\varepsilon}''_{{\rm{rb}}}},\;\; {{\varepsilon }'_{{\rm{rb}}}} = 4.5,\;\; {{\varepsilon }''_{{\rm{rb}}}} = 0.07.\end{split} $
Si, SiO2层和缺陷层的厚度分别为${d_{\rm{a}}} = {{\lambda _0}} / {\sqrt {{{\varepsilon}'_{{\rm{ra}}}}} },$$ {d_{\rm{b}}} = {{{\lambda _0}} / {\sqrt {{{\varepsilon}'_{{\rm{rb}}}}} }}, \;{d_{\rm{s}}} = {{{d_{\rm{b}}}}/ {0.85}}$, 其中${\lambda _0} = $ 20 μm.
图9是平面波垂直入射时该光子晶体的反透射系数, 其中圆圈表示无石墨烯涂层时的反透射系数, 作为对比用虚线给出了无耗情形(${\varepsilon ''_{{\rm{rb}}}} = 0$)时的反透射系数, 可以看到, 有耗时反透射都有所减弱. 图9中实心和空心星形分别代表石墨烯涂层位于光子晶体表面(第1个界面位置)和底层(最后一个界面位置)时的情形, 可以看到, 在禁带下方低频段, 石墨烯涂层的加入使反射增强, 透射减弱, 在禁带内位于底层的石墨烯涂层对反射系数几乎没有影响而位于表面的石墨烯涂层对反射系数影响较大, 使之减弱.
图 9 含石墨烯涂层Si/SiO2周期结构1DPC的反透射系数 (a) 反射系数; (b) 透射系数
Figure9. Modulus of reflective and transmittance coefficients of Si/SiO2 1DPC containing graphene sheet: (a) Reflective coefficient; (b) transmittance coefficient.

图10是平面波垂直入射时该光子晶体的吸收率曲线, 石墨烯涂层增强了光子晶体的吸收, 而且当石墨烯涂层位于表面时, 石墨烯和间隔层构成的表面缺陷导致光的局域化, 使得THz光吸收被增强[18], 在禁带及以上频段, 吸收率显著提高; 当石墨烯涂层位于底层时, 吸收率的提高不明显, 而且主要集中在禁带下方低频段.
图 10 含石墨烯界面Si/SiO2周期结构1DPC的吸收率
Figure10. Absorbance of Si/SiO2 1DPC containing graphene sheet.

图11图12分别是TE和TM波情形下, 该光子晶体的吸收率随频率和入射角变化的伪色彩图, 其中子图(a), (b)和(c)分别是无石墨烯涂层、涂层位于表面和底层时的情形, 整体上看, 位于表面的涂层对光吸收的增强效应更为明显, 而且对TE波的吸收要强于TM波.
图 11 含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TE模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层
Figure11. Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

图 12 含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TM模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层
Figure12. Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

根据相位匹配, 讨论给出了平面波斜入射分层介质时, 有耗层中波矢量实部和虚部的计算方法. 有耗层中复数波矢量的实部和虚部不平行, 传播的平面波非均匀. 代表波振幅传播方向的波矢量虚部沿纵向方向, 即波沿深度方向衰减. 波矢量实部代表相位传播方向, 折射角依此计算. 依据边界条件, 推导了跨越石墨烯界面的传播矩阵, 以及“无限薄”石墨烯层的反透射系数解析式. 推广了传播矩阵方法, 使之可应用于包括含导电界面的有耗分层介质反透射和电波传播特性的计算和分析, 为分层介质与导电界面复合结构的设计和验证提供了支撑.
相关话题/传播 结构 计算 光子 光学

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于超强耦合量子点-纳米机械振子系统的全光学质量传感
    摘要:提出一种复合量子点-纳米机械振子系统,该系统以半导体芯片为基底,量子点嵌入倒置半导体圆锥纳米线的底端,通过光学抽运-探测技术来驱动量子点-纳米机械振子系统,研究该系统中的相干光学特性.通过探测吸收谱给出确定机械振子频率和量子点-纳米机械振子耦合强度的全光学方法.此外,基于该系统理论上提出一种在 ...
    本站小编 Free考研考试 2021-12-29
  • 空位及氮掺杂二维ZnO单层材料性质:第一性原理计算与分子轨道分析
    摘要:采用基于密度泛函理论的第一性原理计算方法,系统地研究了带缺陷的二维类石墨烯结构的ZnO(graphenelike-ZnO,g-ZnO)的几何结构、电子结构、磁性性质和吸收光谱性质.研究的缺陷类型包括锌原子空位(VZn_g-ZnO)、氧原子空位(VO_g-ZnO)、氮原子取代氧原子(NO_g-Z ...
    本站小编 Free考研考试 2021-12-29
  • 结构参数对N极性面GaN/InAlN高电子迁移率晶体管性能的影响
    摘要:基于漂移-扩散传输模型、费米狄拉克统计模型以及Shockley-Read-Hall复合模型等,通过自洽求解薛定谔方程、泊松方程以及载流子连续性方程,模拟研究了材料结构参数对N极性面GaN/InAlN高电子迁移率晶体管性能的影响及其物理机制.结果表明,增加GaN沟道层的厚度(5—15nm)与In ...
    本站小编 Free考研考试 2021-12-29
  • 氧化石墨烯的结构稳定性及硝酸催化作用的第一性原理研究
    摘要:利用第一性原理方法,采用超软赝势库系统研究了硝酸熏蒸石墨烯得到的氧化石墨烯结构的稳定性及电子结构.基于石墨烯正交元胞的2×2超胞模型建立相应的正交晶系硝酸熏蒸氧化石墨烯模型,包含15个碳原子和2个氧原子.结果表明熏蒸后包含碳氧双键的氧化石墨烯结构为能量较低的稳定结构,与实验报道一致.力学稳定性 ...
    本站小编 Free考研考试 2021-12-29
  • 基于凹槽结构抑制AlGaN/GaN高电子迁移率晶体管电流崩塌效应
    摘要:基于双脉冲技术,研究了GaN缓冲层陷阱对AlGaN/GaN高电子迁移率晶体管电流崩塌效应的影响.结果表明,栅边缘漏侧的电场峰值使得沟道电子跃迁至缓冲层,并被缓冲层中的陷阱俘获是造成电流崩塌的主要原因之一.提出了势垒层局部凹槽结构,降低了栅边缘漏侧的电场峰值,使电场分布更加均匀,改善了器件的电流 ...
    本站小编 Free考研考试 2021-12-29
  • 宣布式单光子源宣布效率的宣布测量基相关性
    摘要:对非线性晶体参量下转换过程所产生的纠缠光子对的一路光子进行测量,“宣布”式地在另一路获得单光子,是制备单光子源的一种重要途径.其中,核心的问题是如何得到较高的宣布效率.本文以I类相位匹配偏硼酸钡(β-BaB2O4,BBO)非线性晶体参量下转换过程中所产生的偏振纠缠光子对为例,通过研究其纠缠特性 ...
    本站小编 Free考研考试 2021-12-29
  • 基于全保偏光纤结构的主振荡脉冲非线性放大系统
    摘要:提出了基于全保偏光纤结构的主振荡脉冲非线性放大系统,该系统由基于半导体可饱和吸收镜锁模的直线型光纤振荡器、二级放大结构脉冲非线性光纤放大器和具有负色散的单模传导光纤的脉冲压缩器构成.通过此系统获得了中心波长为1560nm,重复频率为200MHz的超短激光脉冲,脉冲半高全宽为44fs,单脉冲能量 ...
    本站小编 Free考研考试 2021-12-29
  • 碱金属和碱土金属掺杂二维GaN材料电磁特性的第一性原理计算
    摘要:基于密度泛函理论和投影缀加波赝势方法,采用广义梯度近似算法研究了碱金属(Li,Na,K和Rb)和碱土金属(Be,Mg和Sr)掺杂二维GaN单层的电子结构和磁学性质.研究表明,除Be原子位于GaN单层平面内之外,其余掺杂原子均略微隆起于平面.通过比较七种掺杂体系在不同环境下的形成能,发现在富N环 ...
    本站小编 Free考研考试 2021-12-29
  • 石墨烯莫尔超晶格体系的拓扑性质及光学研究进展
    摘要:当两个晶格常数或晶格转角不同的二维材料叠加在一起时会出现周期性的莫尔条纹结构,这种莫尔超晶格形成了一个新的二维周期势,可以大大改变原有体系的物理性质.最近石墨烯与石墨烯、石墨烯与六方氮化硼形成的莫尔超晶格提供了一个非常有趣的体系,在该体系中石墨烯的电子能带结构发生了根本性的改变,在原本的能谷处 ...
    本站小编 Free考研考试 2021-12-29
  • 拓扑光子学研究进展
    摘要:受凝聚态中拓扑相和拓扑相变概念的启发,一种基于拓扑能带论的新的研究领域——拓扑光子学正在兴起,它突破了传统基于实空间光场叠加原理和倒空间固体能带色散理论的光场调控思想,提供了一种新颖的光场调控机制和丰富的输运和光操控性质.例如,背散射抑制且缺陷免疫的边界输运特性、自旋轨道依赖的选择传输特性、高 ...
    本站小编 Free考研考试 2021-12-29