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

电磁偏置各向异性石墨烯界面的传播矩阵

本站小编 Free考研考试/2021-12-29
摘要:基于电磁场边界条件和相位匹配, 推导出电、磁偏置下呈各向异性的石墨烯导电界面的传播矩阵, 并进一步给出各向异性石墨烯界面的反透射系数解析解; 该传播矩阵耦合了基本的横电波和横磁波极化, 并包括偏置电、磁场的影响. 将跨石墨烯界面传播矩阵嵌入各向同性分层介质传播矩阵, 获得的新传播矩阵可用于解析分析平面电磁波以任意角度入射含各向异性石墨烯界面层状介质时的传播和反透射特性, 并且为分层介质与各向异性导电界面复合结构的相关分析和设计提供了一种非常简单的工具.
关键词: 各向异性石墨烯/
电、磁偏置/
分层介质/
传播矩阵

English Abstract


全文HTML

--> --> -->
1.引　言
石墨烯是由一层碳原子构成的二维蜂窝状晶体结构, 其厚度只有一个原子大小, 约为0.34 nm, 是世界上最薄的半导体材料. 石墨烯出色的力学特性和承载超大电流的能力, 使其成为未来微纳米技术中最有前途的材料之一[1-7], 并且在2004年通过微机械剥离法制备成功[8-11].
相对于其他单原子金属薄膜, 石墨烯在热力学上很稳定, 并且在载流子浓度消失的情况下仍具有最小的导电性[9]; 同时, 石墨烯是一种零带隙半导体, 利用静电和/或静磁偏置的电场和霍尔效应可以调节其导电性[8,10]. 讨论分析石墨烯与电磁波相互作用问题时, 可以将其看成具有导电性的“无限薄”表面[12], 其导电性由表面导电率$\sigma $描述[13]; 微观量子动力学方法推导出了$\sigma $的数学模型, 其中还包括静电和/或静磁偏置情况下的重要影响—在存在偏压和偏磁的情况下, 石墨烯呈各向异性, 其表面导电率为张量${{\sigma}} $[14]. 此外应注意的是, 即使没有偏置场, 在太赫兹(terahertz, THz)频率以上, 石墨烯也具有空间色散效应, 这时亦呈现各向异性特性[15].
本文主要研究在太赫兹频率以下, 存在偏置外静电和静磁场时, 各向异性的“无限薄”石墨烯与平面电磁波的相互作用: 利用电磁场边界条件和相位匹配条件[16], 推导出平面电磁波在各向异性石墨烯平面两侧的跨界面传播矩阵(propagation matrix, PM), 其中耦合了基本的横电(transverse electric, TE)波和横磁(transverse magnetic, TM)波极化, 同时包括偏置电、磁场的影响; 基于该传播矩阵, 进一步给出石墨烯界面的反透射系数解析解; 将跨石墨烯界面传播矩阵嵌入各向同性分层介质传播矩阵[17](TE和TM波合并), 获得的新传播矩阵可用于解析分析平面电磁波以任意角度入射含各向异性石墨烯界面层状介质时的传播和反透射特性. 本文中时谐因子取$\exp \left( { - {\rm{j}}\omega t} \right)$.
2.电、磁偏置各向异性石墨烯界面的传播矩阵
2
2.1.各向异性石墨烯薄层的表面电导率张量
-->

2.1.各向异性石墨烯薄层的表面电导率张量

设一“无限薄”石墨烯层(各向异性)平行于xoy平面, 如图1所示, 该石墨烯层将全空间分开为上、下两部分, 设其分别为均匀各向同性介质${\varepsilon _l}, {\mu _l}$和${\varepsilon _{l + 1}}, {\mu _{l + 1}}$(介电系数$\varepsilon $和磁导系数$\mu $可以是复数, 代表有耗介质).
图 1 “无限薄”石墨烯层
Figure1. “Infinitely thin” graphene sheet.

在不考虑太赫兹情形下石墨烯的空间色散效应[15], 石墨烯层的空间电磁特性可以通过垂直方向的外置偏置电、磁场${E_{{\rm{bias}}}}\hat{{z}}, {B_{{\rm{bias}}}}\hat{{z}}$调节, 一般地, 石墨烯层表面电导率可表示为如下张量形式[14,18-21]:
$\begin{split}& {{\sigma }} = \left[ {\begin{array}{*{20}{c}} {{\sigma _{xx}}}&{{\sigma _{xy}}} \\ {{\sigma _{yx}}}&{{\sigma _{yy}}} \end{array}} \right] , \;\; {\sigma _{xx}} = {\sigma _{yy}} = {\sigma _{\rm{D}}} , \\& {\sigma _{xy}} = - {\sigma _{yx}} = {\sigma _{\rm{O}}} , \;{{\sigma }} = \left[ {\begin{array}{*{20}{c}} {{\sigma _{\rm{D}}}}&{{\sigma _{\rm{O}}}} \\ { - {\sigma _{\rm{O}}}}&{{\sigma _{\rm{D}}}} \end{array}} \right], \end{split}$
其中
$\begin{split} & {\sigma _{\rm{D}}} = {\sigma _{{\rm{DR}}}} + {\rm{j}}{\sigma _{{\rm{DI}}}}= - {\rm{j}}\dfrac{{{e^2}v_{\rm{F}}^2\left| {e{B_{{\rm{bias}}}}} \right|\left( {\omega + {\rm{j}}2\varGamma } \right)\hbar }}{\pi } \\& \cdot \sum\limits_{n = 0}^\infty\! \!\bigg\{ \dfrac{{{n_{\rm{F}}}( {{M_n}} ) \!-\! {n_{\rm{F}}}( {{M_{n + 1}}} ) \!+\! {n_{\rm{F}}}( { - {M_{n + 1}}} ) \!-\! {n_{\rm{F}}}( { - {M_n}} )}}{{{{( {{M_{n + 1}} \!-\! {M_n}} )}^2} \!-\! {{( {\omega \!+\! {\rm{j}}2\varGamma } )}^2}{\hbar ^2}}} \\ & \cdot \Big( {1 - \frac{{{\varDelta ^2}}}{{{M_n}{M_{n + 1}}}}} \Big)\frac{1}{{{M_{n + 1}} - {M_n}}} \\ & + \dfrac{{{n_{\rm{F}}}( { - {M_n}} ) - {n_{\rm{F}}}( {{M_{n + 1}}} ) \!+\! {n_{\rm{F}}}( { - {M_{n + 1}}} ) - {n_{\rm{F}}}( {{M_n}} )}}{{{{( {{M_{n + 1}} + {M_n}})}^2} - {{( {\omega + {\rm{j}}2\varGamma } )}^2}{\hbar ^2}}} \\ & \cdot \left( {1 + \frac{{{\varDelta ^2}}}{{{M_n}{M_{n + 1}}}}} \right)\frac{1}{{{M_{n + 1}} + {M_n}}} \bigg\}, \\[-15pt]\end{split} $
$\begin{split} & {\sigma _{\rm{O}}} = {\sigma _{{\rm{OR}}}} + {\rm{j}}{\sigma _{{\rm{OI}}}} = - \dfrac{{{e^3}v_{\rm{F}}^2{B_{{\rm{bias}}}}}}{\pi } \sum\limits_{n = 0}^\infty \bigg\{ \big[ {n_{\rm{F}}}\left( {{M_n}} \right) \\ & - {n_{\rm{F}}}\left( {{M_{n + 1}}} \right) + {n_{\rm{F}}}\left( { - {M_{n + 1}}} \right) - {n_{\rm{F}}}\left( { - {M_n}} \right) \big] \\ &\cdot\bigg[ \left( \!{1 \!-\! \frac{{{\varDelta ^2}}}{{{M_n}{M_{n + 1}}}}} \right)\!\frac{1}{{{{\left( {{M_{n + 1}} \!-\! {M_n}} \right)}^2} - {{\left( {\omega + {\rm{j}}2\varGamma } \right)}^2}{\hbar ^2}}} \\ & \!+\! \left( {1 + \frac{{{\varDelta ^2}}}{{{M_n}{M_{n + 1}}}}} \right)\frac{1}{{{{\big({{M_{n + 1}} \!+\! {M_n}} \big)}^2} \!-\! {{( {\omega \!+\! {\rm{j}}2\varGamma })}^2}{\hbar ^2}}} \bigg] \bigg\}, \end{split}$
及
${M_n} = \sqrt {{\varDelta ^2} + 2nv_{\rm{F}}^2\left| {e{B_{{\rm{bias}}}}} \right|\hbar } ,$
其中e, $\hbar $是电子电量绝对值、约化普朗克常量, ${v_{\rm{F}}} \cong {10^6}\;{\rm{m}}/{\rm{s}}$为石墨烯中的费米速率, $\varGamma = {1 / {\left( {2\tau } \right)}}$是散射率($\tau $为电子弛豫时间), $\varDelta $表示激子能隙(室温下近似为0), 而费米-狄拉克分布${n_{\rm{F}}}$可表示为
${n_{\rm{F}}}\left( \varepsilon \right) = 1/\big[{{1 + {{\rm{e}}^{{{\left( {\varepsilon - {\mu _{\rm{c}}}} \right)} / {\left( {{k_{\rm{B}}}T} \right)}}}}}}\big],$
其中$\varepsilon $表示能量; ${k_{\rm{B}}}$是玻尔兹曼常数; ${\mu _{\rm{c}}}$是依赖于偏置电场${E_{{\rm{bias}}}}$的化学势, 关系为[15]
$\dfrac{{{\varepsilon _0}\pi {\hbar ^2}v_{\rm{F}}^2}}{e}{E_{{\rm{bias}}}} \!=\! \int_0^{ + \infty } {\varepsilon [ {{n_{\rm{F}}}(\varepsilon) \!-\! {n_{\rm{F}}}( {\varepsilon \!+\! 2{\mu _{\rm{c}}}} )} ]{\rm{d}}\varepsilon } .$
将(5)式代入(6)式, 并通过数值求解可获得化学势与偏置电场关系${\mu _{\rm{c}}}\left( {{E_{{\rm{bias}}}}} \right)$, 如图2所示.
图 2 石墨烯化学势${\mu _{\rm{c}}}\left( {{E_{{\rm{bias}}}}} \right)$与偏置电场${E_{{\rm{bias}}}}$的关系
Figure2. Graphical representation of the relation between the chemical potential ${\mu _{\rm{c}}}\left( {{E_{{\rm{bias}}}}} \right)$ and the electrostatic bias field ${E_{{\rm{bias}}}}$.

应注意, 石墨烯在偏置电场和磁场同时存在时为各向异性, 当偏置电场或/和磁场为0时, 石墨烯为各向同性, 即${\sigma _{\rm{O}}}{{ = }}0$, 而$\sigma = {\sigma _{\rm{D}}}$, 其表达式可参考文献[1].
2
2.2.各向异性石墨烯界面上的电磁场边界条件
-->

2.2.各向异性石墨烯界面上的电磁场边界条件

根据电磁场边界条件[16], 在石墨烯层两侧区域的电、磁场有如下关系[22]:
$ \begin{split}&\hat{{z}}\times \left({{E}}_{l}-{{{E}}}_{l+1}\right)=-{{{J}}}_{\rm{m}}^{\rm{s}}=0\;, \\ &\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}}$为石墨烯表面电流和磁流, 由第二式得
$\begin{split}&\left( {{H_{\left( {l + 1} \right)y}} - {H_{ly}}} \right)\hat{{x}} + \left( {{H_{lx}} - {H_{\left( {l + 1} \right)x}}} \right)\hat{{y}}\\&= {{J}}_{\rm{e}}^{\rm{s}} = J_{{\rm{e}}x}^{\rm{s}}\hat{{x}} + J_{{\rm{e}}y}^{\rm{s}}\hat{{y}}.\end{split}$
又因为
${{J}}_{\rm{e}}^{\rm{s}} = {{\sigma E}},$
即
$J_{{\rm{e}}x}^{\rm{s}} = {\sigma _{\rm{D}}}{E_x} + {\sigma _{\rm{O}}}{E_y},\;\;\;J_{{\rm{e}}y}^{\rm{s}} = - {\sigma _{\rm{O}}}{E_x} + {\sigma _{\rm{D}}}{E_y}.$
代入(8)式得到
$\begin{split}&{H_{\left( {l + 1} \right)x}} - {H_{lx}} = {\sigma _{\rm{O}}}{E_x} - {\sigma _{\rm{D}}}{E_y},\\&{H_{\left( {l + 1} \right)y}} - {H_{ly}} = {\sigma _{\rm{D}}}{E_x} + {\sigma _{\rm{O}}}{E_y}.\end{split}$
而根据边界条件(7)式的第一式可知分界面(石墨烯层)上的电场切向分量是连续的, 即
${E_{\left( {l + 1} \right)x,\left( {l + 1} \right)y}} - {E_{lx,ly}} = 0.$
所以(11)式中等号右边的${E_x}$, ${E_y}$可以取介质l或$l + 1$区域中电场的x, y分量.
2
2.3.各向异性石墨烯界面的前向传播矩阵
-->

2.3.各向异性石墨烯界面的前向传播矩阵

下来我们考虑平面电磁波入射石墨烯界面时, 界面两侧场分量间的关系. 根据电磁波理论, 平面电磁波斜入射界面垂直于z轴的层状各向同性介质时, 在介质的任意l层区域中, 波可按z轴分为TE和TM平面波, 波方程分别为[16]
$ \begin{split}&{E}_{ly}=\left[{A}_{l}^{\rm{TE}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)+{B}_{l}^{\rm{TE}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\right]{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\&{H}_{lx}=\frac{-{k}_{lz}}{\omega {\mu }_{l}}\Big[{A}_{l}^{\rm{TE}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)-{B}_{l}^{\rm{TE}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\Big]\\&\;\;\qquad\times{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\ &{H}_{lz}=\frac{{k}_{x}}{\omega {\mu }_{l}}\Big[{A}_{l}^{\rm{TE}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)+{B}_{l}^{\rm{TE}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\Big]\\&\qquad\;\;\times{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\[-13pt]\end{split}$
和
$ \begin{split}&{H}_{ly}=\left[{A}_{l}^{\rm{TM}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)\!+\!{B}_{l}^{\rm{TM}}{\rm{exp}}\left(\!-{\rm{j}}{k}_{lz}z\right)\right]{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\&{E}_{lx}=\frac{{k}_{lz}}{\omega {\varepsilon }_{l}}\Big[{A}_{l}^{\rm{TM}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)-{B}_{l}^{\rm{TM}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\Big]\\&\qquad\;\;\;\times{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\ &{E}_{lz}=\frac{-{k}_{x}}{\omega {\varepsilon }_{l}}\Big[{A}_{l}^{\rm{TM}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right)+{B}_{l}^{\rm{TM}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\Big]\\&\qquad\;\;\;\times{\rm{exp}}\left({\rm{j}}{k}_{x}x\right), \\[-13pt]\end{split}$
其中$A_l^{{\rm{TE/TM}}}$和$B_l^{{\rm{TE/TM}}}$分别代表区域l中上行和下行TE/TM波的电/磁场幅值.
设石墨烯层位于l和$l + 1$区域分界面, 由边界条件(11)式和(12)式可解得
$\begin{split}& \Big[ A_{\left( {l + 1} \right)}^{{\rm{TE}}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right),B_{\left( {l + 1} \right)}^{{\rm{TE}}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right),A_{\left( {l + 1} \right)}^{{\rm{TM}}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right),B_{\left( {l + 1} \right)}^{{\rm{TM}}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \Big]^{\rm{T}} \\=\; & {{{U}}_G} \cdot \Big[ A_l^{{\rm{TE}}}\exp \left( {{\rm{j}}{k_{lz}}z} \right),B_l^{{\rm{TE}}}\exp \left( { - {\rm{j}}{k_{lz}}z} \right),A_l^{{\rm{TM}}}\exp \left( {{\rm{j}}{k_{lz}}z} \right),B_l^{{\rm{TM}}}\exp \left( { - {\rm{j}}{k_{lz}}z} \right) \Big]^{\rm{T}} ,\\[-15pt] \end{split} $
其中
${{{U}}_G} = \dfrac{1}{2}\left[ {\begin{array}{*{20}{c}} {1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + \kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}{\sigma _{\rm{D}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + \kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}{\sigma _{\rm{D}}}}&{ - \dfrac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}}}{{\kappa _l^{{\rm{TM}}}}}{\sigma _{\rm{O}}}}&{\dfrac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}}}{{\kappa _l^{{\rm{TM}}}}}{\sigma _{\rm{O}}}} \\ {1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - \kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}{\sigma _{\rm{D}}}}&{1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - \kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}{\sigma _{\rm{D}}}}&{\dfrac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}}}{{\kappa _l^{{\rm{TM}}}}}{\sigma _{\rm{O}}}}&{ - \dfrac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}}}{{\kappa _l^{{\rm{TM}}}}}{\sigma _{\rm{O}}}} \\ {{\sigma _{\rm{O}}}}&{{\sigma _{\rm{O}}}}&{1 + p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _l^{{\rm{TM}}}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _l^{{\rm{TM}}}}}} \\ {{\sigma _{\rm{O}}}}&{{\sigma _{\rm{O}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _l^{{\rm{TM}}}}}}&{1 + p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _l^{{\rm{TM}}}}}} \end{array}} \right],$
称为各向异性石墨烯分界面的前向传播矩阵, 式中
$\begin{split}&\kappa _{l}^{\rm{TM}}=\frac{\omega \varepsilon_l }{k_{lz}},\;\; \kappa _{l}^{\rm{TE}}=\frac{\omega \mu_l }{k_{lz}},\;\;p_{\left( {l + 1} \right)l}^{{\rm{TM}}} = \frac{{{\varepsilon _{\left( {l + 1} \right)}}{k_{lz}}}}{{{\varepsilon _l}{k_{\left( {l + 1} \right)z}}}} = \frac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TM}}}}}{{\kappa _l^{{\rm{TM}}}}},\;\;p_{\left( {l + 1} \right)l}^{{\rm{TE}}} = \frac{{{\mu _{\left( {l + 1} \right)}}{k_{lz}}}}{{{\mu _l}{k_{\left( {l + 1} \right)z}}}} = \frac{{\kappa _{\left( {l + 1} \right)}^{{\rm{TE}}}}}{{\kappa _l^{{\rm{TE}}}}}.\end{split}$
由于各向异性, 石墨烯层分界面两侧的TE和TM平面波出现耦合, 单纯TE/TM波入射会产生同时含TE和TM波分量的反射和透射波.
当${\sigma _{\rm{O}}} = 0$时, 石墨烯为各向同性, (15)式可分解为TE, TM波对应的方程
$ \begin{split}&{\left[{A}_{\left(l+1\right)}^{\rm{TE}}{\rm{exp}}\left({\rm{j}}{k}_{\left(l+1\right)z}z\right),{B}_{\left(l+1\right)}^{\rm{TE}}{\rm{exp}}\left(-{\rm{j}}{k}_{\left(l+1\right)z}z\right)\right]}^{\rm{T}}={{V}}^{\rm{TE}}\cdot {\left[{A}_{l}^{\rm{TE}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right),{B}_{l}^{\rm{TE}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\right]}^{\rm{T}}\;, \\ &{{V}}^{\rm{TE}}=\dfrac{1}{2}\left[\begin{array}{cc}1+\dfrac{{k}_{lz}{\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\mu }_{l}}+\dfrac{\omega {\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}}{\sigma }_{\rm{D}}& 1-\dfrac{{k}_{lz}{\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\mu }_{l}}+\dfrac{\omega {\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}}{\sigma }_{\rm{D}}\\ 1-\dfrac{{k}_{lz}{\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\mu }_{l}}-\dfrac{\omega {\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}}{\sigma }_{\rm{D}}& 1+\dfrac{{k}_{lz}{\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\mu }_{l}}-\dfrac{\omega {\mu }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}}{\sigma }_{\rm{D}}\end{array}\right], \end{split}$
$ \begin{split}&{\left[{A}_{\left(l+1\right)}^{\rm{TM}}{\rm{exp}}\left({\rm{j}}{k}_{\left(l+1\right)z}z\right),{B}_{\left(l+1\right)}^{\rm{TM}}{\rm{exp}}\left(-{\rm{j}}{k}_{\left(l+1\right)z}z\right)\right]}^{\rm{T}}={{V}}^{\rm{TM}}\cdot {\left[{A}_{l}^{\rm{TM}}{\rm{exp}}\left({\rm{j}}{k}_{lz}z\right),{B}_{l}^{\rm{TM}}{\rm{exp}}\left(-{\rm{j}}{k}_{lz}z\right)\right]}^{\rm{T}}\;, \\ &{{V}}^{\rm{TM}}=\dfrac{1}{2}\left[\begin{array}{cc}1+\dfrac{{k}_{lz}{\varepsilon }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\varepsilon }_{l}}+\dfrac{{k}_{lz}}{\omega {\varepsilon }_{l}}{\sigma }_{\rm{D}}& 1-\dfrac{{k}_{lz}{\varepsilon }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\varepsilon }_{l}}-\dfrac{{k}_{lz}}{\omega {\varepsilon }_{l}}{\sigma }_{\rm{D}}\\ 1-\dfrac{{k}_{lz}{\varepsilon }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\varepsilon }_{l}}+\dfrac{{k}_{lz}}{\omega {\varepsilon }_{l}}{\sigma }_{\rm{D}}& 1+\dfrac{{k}_{lz}{\varepsilon }_{\left(l+1\right)}}{{k}_{\left(l+1\right)z}{\varepsilon }_{l}}-\dfrac{{k}_{lz}}{\omega {\varepsilon }_{l}}{\sigma }_{\rm{D}}\end{array}\right]\;\;, \end{split}$
其中${{{V}}^{{\rm{TE}}}}$, ${{{V}}^{{\rm{TM}}}}$即为文献[17]中所给出的各向同性石墨烯分界面的前向传播矩阵.
2
2.4.各向异性石墨烯界面的反透射系数
-->

2.4.各向异性石墨烯界面的反透射系数

设石墨烯薄层位于$z = 0$处, 平面波自上往下斜入射, 下来分别讨论TE和TM入射波情形, 此时设入射区为介质1, 即$l = 1$, 则透射区$l + 1$为介质2.
3
2.4.1.TE波入射
-->
2.4.1.TE波入射
此时在1区没有下行的TM波, $B_1^{{\rm{TM}}} = 0$, 将下行即入射的TE波$B_1^{{\rm{TE}}}$看做已知参量, 是入射波电场幅值; 而2区为透射区, 所以没有上行波, $A_2^{{\rm{TE}}} = 0, A_2^{{\rm{TM}}} = 0$. 由(15)式可得
$\left[\! {\begin{array}{*{20}{c}} 0 \\ {B_2^{{\rm{TE}}}} \\ 0 \\ {B_2^{{\rm{TM}}}} \end{array}} \!\right] = \dfrac{1}{2}\!\left[\!\! {\begin{array}{*{20}{c}} {p_{21}^{{\rm{TE}}} + \kappa _2^{{\rm{TE}}}{\sigma _{\rm{D}}} + 1}&{ - p_{21}^{{\rm{TE}}} + \kappa _2^{{\rm{TE}}}{\sigma _{\rm{D}}} + 1}&{ - \dfrac{{\kappa _2^{{\rm{TE}}}}}{{\kappa _1^{{\rm{TM}}}}}{\sigma _{\rm{O}}}}&{\frac{{\kappa _2^{{\rm{TE}}}}}{{\kappa _1^{{\rm{TM}}}}}{\sigma _{\rm{O}}}} \\ { - p_{21}^{{\rm{TE}}} - \kappa _2^{{\rm{TE}}}{\sigma _{\rm{D}}} + 1}&{p_{21}^{{\rm{TE}}} - \kappa _2^{{\rm{TE}}}{\sigma _{\rm{D}}} + 1}&{\dfrac{{\kappa _2^{{\rm{TE}}}}}{{\kappa _1^{{\rm{TM}}}}}{\sigma _{\rm{O}}}}&{ - \dfrac{{\kappa _2^{{\rm{TE}}}}}{{\kappa _1^{{\rm{TM}}}}}{\sigma _{\rm{O}}}} \\ {{\sigma _{\rm{O}}}}&{{\sigma _{\rm{O}}}}&{p_{21}^{{\rm{TM}}} + \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _1^{{\rm{TM}}}}} + 1}&{ - p_{21}^{{\rm{TM}}} - \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _1^{{\rm{TM}}}}} + 1} \\ {{\sigma _{\rm{O}}}}&{{\sigma _{\rm{O}}}}&{ - p_{21}^{{\rm{TM}}} + \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _1^{{\rm{TM}}}}} + 1}&{p_{21}^{{\rm{TM}}} - \dfrac{{{\sigma _{\rm{D}}}}}{{\kappa _1^{{\rm{TM}}}}} + 1} \end{array}}\!\! \right] \cdot \left[\! {\begin{array}{*{20}{c}} {A_1^{{\rm{TE}}}} \\ {B_1^{{\rm{TE}}}} \\ {A_1^{{\rm{TM}}}} \\ 0 \end{array}} \!\right],$
其中
$\begin{split}&\kappa _1^{{\rm{TE}}} = \frac{{\omega {\mu _1}}}{{{k_{1z}}}},\;\kappa _2^{{\rm{TE}}} = \frac{{\omega {\mu _2}}}{{{k_{2z}}}},\;\kappa _1^{{\rm{TM}}} = \frac{{\omega {\varepsilon _1}}}{{{k_{1z}}}},\;\;\kappa _2^{{\rm{TM}}} = \frac{{\omega {\varepsilon _2}}}{{{k_{2z}}}},\;p_{21}^{{\rm{TE/TM}}} = \frac{{\kappa _2^{{\rm{TE/TM}}}}}{{\kappa _1^{{\rm{TE/TM}}}}}.\end{split}$
定义TE波入射石墨烯单层时的反透射系数为
${R^\text{{{TE-TE}}}} \!=\! \dfrac{{A_1^{{\rm{TE}}}}}{{B_1^{{\rm{TE}}}}},\;{R^\text{{{TM}}-{{TE}}}}\! =\! {\eta _1}\dfrac{{A_1^{{\rm{TM}}}}}{{B_1^{{\rm{TE}}}}} \!=\! \sqrt {\dfrac{{{\mu _1}}}{{{\varepsilon _1}}}} \dfrac{{A_1^{{\rm{TM}}}}}{{B_1^{{\rm{TE}}}}},\;\;{T^\text{{{TE-TE}}}} \!=\! \frac{{B_2^{{\rm{TE}}}}}{{B_1^{{\rm{TE}}}}},\;{T^\text{{{TM}}-{{TE}}}} \!=\! {\eta _2}\frac{{B_2^{{\rm{TM}}}}}{{B_1^{{\rm{TE}}}}}\! =\!\! \sqrt {\frac{{{\mu _2}}}{{{\varepsilon _2}}}} \frac{{B_2^{{\rm{TM}}}}}{{B_1^{{\rm{TE}}}}},$
其分别为反射TE波、TM波, 透射TE波、TM波与入射TE波的电场幅值比. 由于$A_1^{{\rm{TM}}}, \;B_2^{{\rm{TM}}}$表示的是TM波的磁场幅值, 所以(22)式中利用各自区域介质的波阻抗将其转换为电场幅值. 求解(20)式可得TE波入射石墨烯单层时的反透射系数
$\begin{split}&{R}^{\rm{\text{TE-TM}}}=\frac{-2{\eta }_{1}{\sigma }_{\rm{O}}{\kappa }_{1}^{\rm{TM}}{\kappa }_{2}^{\rm{TE}}}{\left({\kappa }_{1}^{\rm{TE}}+{\kappa }_{2}^{\rm{TE}}+{\kappa }_{1}^{\rm{TE}}{\kappa }_{2}^{\rm{TE}}{\sigma }_{\rm{D}}\right)\left({\kappa }_{1}^{\rm{TM}}+{\kappa }_{2}^{\rm{TM}}+{\sigma }_{\rm{D}}\right)+{\kappa }_{1}^{\rm{TM}}{\kappa }_{2}^{\rm{TE}}{\eta }_{1}{}^{2}{\sigma }_{\rm{O}}{}^{2}}\;, \\ &{R}^{\rm{\text{TE-TM}}}=\frac{{p}_{21}^{\rm{TE}}-{\kappa }_{2}^{\rm{TE}}{\sigma }_{\rm{D}}-1}{{p}_{21}^{\rm{TE}}+{\kappa }_{2}^{\rm{TE}}{\sigma }_{\rm{D}}+1}+\frac{{\kappa }_{2}^{\rm{TE}}}{{\eta }_{1}{\kappa }_{1}^{\rm{TM}}}\left(\frac{{\sigma }_{\rm{O}}}{{p}_{21}^{\rm{TE}}+{\kappa }_{2}^{\rm{TE}}{\sigma }_{\rm{D}}+1}\right){R}^{\rm{\text{TE-TM}}}\;, \\ &{T}^{\rm{\text{TE-TM}}}=-\frac{{\eta }_{2}}{{\eta }_{1}}{p}_{21}^{\rm{TM}}\cdot {R}^{\rm{\text{TE-TM}}}, \;\;\;\;\;\;\;\;\;{T}^{\rm{\text{TE-TM}}}={R}^{\rm{\text{TE-TM}}}+1\;{.}\end{split}$

3
2.4.2.TM波入射
-->
2.4.2.TM波入射
此时在1区没有下行的TE波, $B_1^{{\rm{TE}}} = 0$, 将下行即入射的TM波$B_1^{{\rm{TM}}}$看做已知参量, 是入射波磁场幅值; 而2区为透射区, 所以没有上行波, $A_2^{{\rm{TE}}} = 0, A_2^{{\rm{TM}}} = 0$. 按照和2.4.1小节类似讨论方法可得TM波入射石墨烯单层时的反透射系数:
$ \begin{split}&{R}^{\rm{\text{TE-TM}}}=\frac{-2{\eta }_{1}{\sigma }_{\rm{O}}{\kappa }_{1}^{\rm{TM}}{\kappa }_{2}^{\rm{TE}}}{\left({\kappa }_{1}^{\rm{TE}}+{\kappa }_{2}^{\rm{TE}}+{\kappa }_{1}^{\rm{TE}}{\kappa }_{2}^{\rm{TE}}{\sigma }_{\rm{D}}\right)\left({\kappa }_{1}^{\rm{TM}}+{\kappa }_{2}^{\rm{TM}}+{\sigma }_{\rm{D}}\right)+{\kappa }_{1}^{\rm{TM}}{\kappa }_{2}^{\rm{TE}}{\eta }_{1}{}^{2}{\sigma }_{\rm{O}}{}^{2}}\;, \\ &{R}^{\rm{\text{TE-TM}}}=\frac{{\kappa }_{2}^{\rm{TM}}-{\kappa }_{1}^{\rm{TM}}+{\sigma }_{\rm{D}}}{{\kappa }_{2}^{\rm{TM}}+{\kappa }_{1}^{\rm{TM}}+{\sigma }_{\rm{D}}}-\frac{{\eta }_{1}{\kappa }_{1}^{\rm{TM}}{\sigma }_{\rm{O}}}{{\kappa }_{1}^{\rm{TM}}+{\kappa }_{2}^{\rm{TM}}+{\sigma }_{\rm{D}}}{R}^{\rm{{\text{TE-TM}}}}\;, \\ &{T}^{\rm{\text{TE-TM}}}=\frac{{\eta }_{1}}{{\eta }_{2}}{R}^{\rm{{\text{TE-TM}}}},\;\;\;\;{T}^{\rm{\text{TE-TM}}}={p}_{21}^{\rm{TM}}\left(1-{R}^{\rm{\text{TE-TM}}}\right)\;, \end{split}$
其中
$\begin{split}&{R^{{\rm{TM \text{-} TM}}}} = \frac{{A_1^{{\rm{TM}}}}}{{B_1^{{\rm{TM}}}}},\\& {R^{{\rm{TE \text{-} TM}}}} = \frac{{A_1^{{\rm{TE}}}}}{{{\eta _1}B_1^{{\rm{TM}}}}} = \sqrt {\frac{{{\varepsilon _1}}}{{{\mu _1}}}} \frac{{A_1^{{\rm{TE}}}}}{{B_1^{{\rm{TM}}}}},~~\\& {T^{{\rm{TM \text{-} TM}}}} = \frac{{B_2^{{\rm{TM}}}}}{{B_1^{{\rm{TM}}}}},\\& {T^{{\rm{TE\text{-} TM}}}} = \frac{{B_2^{{\rm{TE}}}}}{{{\eta _2}B_1^{{\rm{TM}}}}} = \sqrt {\frac{{{\varepsilon _2}}}{{{\mu _2}}}} \frac{{B_2^{{\rm{TE}}}}}{{B_1^{{\rm{TM}}}}},\end{split}$
分别是反射TM波、TE波, 透射TM波、TE波与入射TM波的磁场幅值比, 由于$A_1^{{\rm{TE}}}, \;B_2^{{\rm{TE}}}$表示的是TE波的电场幅值, 所以(25)式中利用各自区域介质的波阻抗将其转换为磁场幅值.
(23)式和(24)式即为平面波入射石墨烯单层时的反透射系数公式(含同极化和交叉极化). 另, 由二式可知, 石墨烯单层的交叉反射系数相同, ${R^{{\rm{TM}} \text{-} {\rm{TE}}}} = {R^{{\rm{TE \text{-} TM}}}}$.
当石墨烯界面两侧为同种介质, 在均匀平面波入射时, 有
$ \begin{split}&{\eta }_{1}={\eta }_{2}=\eta =\sqrt{\mu /\varepsilon },\\&{k}_{1z}={k}_{2z}={k}_{z}=k{\rm{cos}}\theta =\omega \sqrt{\mu \varepsilon }{\rm{cos}}\theta \;, \\&{\kappa }_{1}^{\rm{TE/TM}}={\kappa }_{2}^{\rm{TE/TM}}={\kappa }^{\rm{TE/TM}},\;{p}_{21}^{\rm{TE/TM}}=1,\\&{\kappa }^{\rm{TE}}=\frac{\omega \mu }{{k}_{z}}=\eta {\rm{cos}}\theta,\;{\kappa }^{\rm{TM}}=\frac{\omega \varepsilon }{{k}_{z}}=\frac{{\rm{cos}}\theta }{\eta }\;{.}\end{split}$
代入(23)式和(24)式得到
$ \begin{split}&{R}^{\rm{TM\text{-}TE}}={R}^{\rm{TE\text{-}TM}}=\\&\frac{-2\eta {\sigma }_{\rm{O}}{\rm{cos}}\theta }{4{\rm{cos}}\theta +2\eta \left({{\rm{cos}}}^{2}\theta +1\right){\sigma }_{\rm{D}}+{\eta }^{2}{\rm{cos}}\theta \left({\sigma }_{\rm{D}}{}^{2}+{\sigma }_{\rm{O}}{}^{2}\right)}\;, \\ &{R}^{\rm{TE\text{-}TE}}=\frac{-\eta {\rm{cos}}\theta \left[{\sigma }_{\rm{D}}-\dfrac{{\sigma }_{\rm{O}}}{{\rm{cos}}\theta }{R}^{\rm{TM\text{-}TE}}\right]}{2+\eta {\rm{cos}}\theta {\sigma }_{\rm{D}}},\\&{R}^{\rm{TM\text{-}TM}}=\frac{\eta \left[{\sigma }_{\rm{D}}-{\sigma }_{\rm{O}}{\rm{cos}}\theta {R}^{\rm{TE\text{-}TM}}\right]}{2{\rm{cos}}\theta +\eta {\sigma }_{\rm{D}}}\;, \\ &{T}^{\rm{TE\text{-}TE}}={R}^{\rm{TE\text{-}TE}}+1,\;\;{T}^{\rm{TM\text{-}TM}}=1-{R}^{\rm{TM\text{-}TM}}\;, \\ &{T}^{\rm{TM\text{-}TE}}=-{R}^{\rm{TM\text{-}TE}},\;\;{T}^{\rm{TE\text{-}TM}}={R}^{\rm{TE\text{-}TM}}\;{.}\\[-12pt]\end{split}$
可见, 均匀各向同性介质中石墨烯薄层对入射TE和TM波的交叉极化反射系数以及入射TM波的交叉极化透射系数是相等的, 且与入射TE波的交叉极化透射系数模值相等、相位相差${\text{π}}$. 令上式中$\eta = {\eta _0}$, 即得真空中石墨烯薄层的反透射系数.
3.含各向异性石墨烯界面的层状介质传播矩阵
考虑如图3所示各向同性分层介质, 假设其中某一个或多个界面被嵌入“无限薄”石墨烯, 这时由于各向异性石墨烯的交叉极化效应, 使得该层状结构对入射平面波产生交叉极化方向的反射和透射, 所以我们把透射区t和入射区0的场值关系表示为
图 3 层状介质
Figure3. Stratified medium.

$\begin{split}&~~~~\left[ {\begin{array}{*{20}{c}} 0 \\ {\left( {{T^{{\rm{TE \text{-} TE}}}}B_0^{{\rm{TE}}} + {\eta _t}{T^{{\rm{TE \text{-} TM}}}}B_0^{{\rm{TM}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \\ 0 \\ {\left( {{T^{{\rm{TM \text{-} TM}}}}B_0^{{\rm{TM}}} + \dfrac{{{T^{{\rm{TM}} \text{-} {\rm{TE}}}}B_0^{{\rm{TE}}}}}{{{\eta _t}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \end{array}} \right]\\ &= {{{V}}_{t0}}\!\!\left[\!\!\! {\begin{array}{*{20}{c}} {\left( {{R^{{\rm{TE \text{-} TE}}}}B_0^{{\rm{TE}}} + {\eta _0}{R^{{\rm{TE\text{-} TM}}}}B_0^{{\rm{TM}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {B_0^{{\rm{TE}}}\exp ( - {\rm{j}}{k_{0z}}{z_0})} \\ {\left( {{R^{{\rm{TM \text{-} TM}}}}B_0^{{\rm{TM}}} + \dfrac{{{R^{{\rm{TM}}\text{-} {\rm{TE}}}}B_0^{{\rm{TE}}}}}{{{\eta _0}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {B_0^{{\rm{TM}}}\exp ( - {\rm{j}}{k_{0z}}{z_0})} \end{array}}\!\!\! \right].\end{split}$
上式等号左边和右边的$4 \times 1$矩阵各行依次分别表示透射区t和入射区0中的上行TE波电场、下行TE波电场、上行TM波电场、下行TM波电场. 矩阵${{{V}}_{t0}}$为从入射点${M_0}$到透射点${N_t}$的前向传播矩阵.
当不存在石墨烯界面时, 各向同性分层介质的前向传播矩阵有两组${{V}}_{t0}^{{\rm{TE}}}$和${{V}}_{t0}^{{\rm{TM}}}$, 分别对应TE和TM两种线极化波情形[16,17]. 考虑到石墨烯界面各向异性特性产生的交叉极化效应, 我们首先将这两组传播矩阵整合为一组, 如下
$ \begin{split}&{{{V}}}_{t0}=\left[\begin{array}{cc}{{{V}}}_{t0}^{\rm{TE}}& 0\\ 0& {{{V}}}_{t0}^{\rm{TM}}\end{array}\right],\;\;\;\;\;\;{{{V}}}_{t0}^{\rm{TE/TM}}={{{U}}}_{n}^{\rm{TE/TM}}\cdot {{{V}}}_{n}^{\rm{TE/TM}}\cdots {{{U}}}_{1}^{\rm{TE/TM}}\cdot {{{V}}}_{1}^{\rm{TE/TM}}\cdot {{{U}}}_{0}^{\rm{TE/TM}}\;, \\ &{{{U}}}_{l=0\sim n}^{\rm{TE/TM}}=\frac{1}{2}\left(1+{p}_{\left(l+1\right)l}^{\rm{TE/TM}}\right)\left[\begin{array}{cc}1& {R}_{\left(l+1\right)l}^{\rm{TE/TM}}\\ {R}_{\left(l+1\right)l}^{\rm{TE/TM}}& 1\end{array}\right],\;\;\;\;\;\;{R}_{\left(l+1\right)l}^{\rm{TE/TM}}=\frac{1-{p}_{\left(l+1\right)l}^{\rm{TE/TM}}}{1+{p}_{\left(l+1\right)l}^{\rm{TE/TM}}}\;, \\ &{{{V}}}_{l=1\sim n}^{\rm{TE/TM}}=\left[\begin{array}{cc}{\rm{exp}}\left[{\rm{j}}{k}_{\left(l+1\right)z}\left({z}_{l+1}-{z}_{l}\right)\right]& 0\\ 0& {\rm{exp}}\left[-{\rm{j}}{k}_{\left(l+1\right)z}\left({z}_{l+1}-{z}_{l}\right)\right]\end{array}\right]\;, \end{split}$
即
$ \begin{split}&{{{V}}}_{t0}={{{U}}}_{n}\cdot {{{V}}}_{n}\cdots {{{U}}}_{1}\cdot {{{V}}}_{1}\cdot {{{U}}}_{0}\;, \\ &{{{U}}}_{l=0\sim n}=\left[\begin{array}{cc}{{{U}}}_{l}^{\rm{TE}}& 0\\ 0& {{{U}}}_{l}^{\rm{TM}}\end{array}\right],\\ &{{{V}}}_{l=1\sim n}=\left[\begin{array}{cc}{{{V}}}_{l}^{\rm{TE}}& 0\\ 0& {{{V}}}_{l}^{\rm{TM}}\end{array}\right]\;, \end{split}$
其中$V_l^{{\rm{TE/TM}}}$, ${{U}}_l^{{\rm{TE/TM}}}$分别是l层介质中、跨越区域l和$l + 1$界面的传播矩阵.
假设在区域l和$l + 1$的分界面$z = {z_l}$处嵌入石墨烯, 则将(30)式中的矩阵${{{U}}_l}$替换为(15)式和(16)式中的${{{U}}_G}$.
如果是纯粹的各向同性分层介质, (28)式就退化为TE模和TM模两组传播矩阵, 此时${R^{{\rm{TM}} \text{-} {\rm{TE}}}}, \;{T^{{\rm{TM}}\text{-} {\rm{TE}}}}, \;{R^{{\rm{TE \text{-} TM}}}}, \;{T^{{\rm{TE \text{-} TM}}}}$都是0.
下面我们分别考虑TE和TM线极化波入射时的情形:
2
3.1.TE波入射
-->

3.1.TE波入射

此时入射TE波电场幅值设为1, 而入射TM波磁场幅值为0, 即$B_0^{{\rm{TE}}} = 1$和$B_0^{{\rm{TM}}} = 0$, 则(28)式成为
$\begin{split}&\left[ {\begin{array}{*{20}{c}} 0 \\ {\left( {{T^{{\rm{TE \text{-} TE}}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \\ 0 \\ {\left( {{{{T^{{\rm{TM}} \text{-} {\rm{TE}}}}} / {{\eta _t}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \end{array}} \right] \\&= {{{V}}_{t0}}\left[ {\begin{array}{*{20}{c}} {\left( {{R^{{\rm{TE \text{-} TE}}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {\exp ( - {\rm{j}}{k_{0z}}{z_0})} \\ {\left( {{{{R^{{\rm{TM}} \text{-} {\rm{TE}}}}} / {{\eta _0}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ 0 \end{array}} \right],\end{split}$
可求得
$ \begin{split}&{R}^{\rm{TE\text{-}TE}}=\frac{{v}_{13}{v}_{32}-{v}_{12}{v}_{33}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0})\;, \\ &{T}^{\rm{TE\text{-}TE}} = \left({v}_{22}+\frac{{v}_{13}{v}_{32}-{v}_{12}{v}_{33}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{v}_{21}\right.\\&\qquad\qquad\left.+\frac{{v}_{12}{v}_{31}\!-\!{v}_{11}{v}_{32}}{{v}_{11}{v}_{33}\!-\!{v}_{13}{v}_{31}}{v}_{23}\!\right)\\& \qquad\qquad \times{\rm{exp}}\left({\rm{j}}{k}_{tz}{z}_{n}\!-\!{\rm{j}}{k}_{0z}{z}_{0}\right), \\ &{R}^{\rm{TM\text{-}TE}}={\eta }_{0}\frac{{v}_{12}{v}_{31}-{v}_{11}{v}_{32}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0})\;, \\ &{T}^{\rm{TM\text{-}TE}}={\eta }_{t}\left({v}_{42}+\frac{{v}_{12}{v}_{31}-{v}_{11}{v}_{32}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{v}_{43}\right.\\&\qquad\qquad\left.+\frac{{v}_{13}{v}_{32}\!-\!{v}_{12}{v}_{33}}{{v}_{11}{v}_{33}\!-\!{v}_{13}{v}_{31}}{v}_{41}\!\right) \\ & \qquad \quad \times{\rm{exp}}\left({\rm{j}}{k}_{tz}{z}_{n}\!-\!{\rm{j}}{k}_{0z}{z}_{0}\right), \end{split}$
其中${v_{ij}}$是${{{V}}_{t0}}$的元素. 当没有石墨烯界面时
$\begin{split} {v_{13}} ={}& {v_{14}} = {v_{23}} = {v_{24}} = {v_{31}} = {v_{32}} \\={}& {v_{41}} = {v_{42}} = 0,\end{split}$
代入(32)式得到
$ \begin{split}&{R}^{\rm{TE\text{-}TE}}=-\left(\dfrac{v_{12}}{v_{11}}\right){\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0}),\\&{R}^{\rm{TM\text{-}TE}}=0\;, \\ &{T}^{\rm{TE\text{-}TE}}=\left({v}_{22}-\frac{{v}_{12}}{{v}_{11}}{v}_{21}\right) \\& \qquad \qquad \times {\rm{exp}}\left({\rm{j}}{k}_{tz}{z}_{n}-{\rm{j}}{k}_{0z}{z}_{0}\right),\\&{T}^{\rm{TM\text{-}TE}}=0\;,\\[-10pt] \end{split}$
即为TE波入射一般分层介质反透射计算式[16,17].
2
3.2.TM波入射
-->

3.2.TM波入射

此时入射TM波磁场场幅值设为1, 而入射TE波磁场幅值为0, 即$B_0^{{\rm{TM}}} = 1$和$B_0^{{\rm{TE}}} = 0$, 则(28)式成为
$\begin{split}&\qquad \left[ {\begin{array}{*{20}{c}} 0 \\ {\left( {{\eta _t}{T^{{\rm{TE \text{-} TM}}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \\ 0 \\ {\left( {{T^{{\rm{TM \text{-}TM}}}}} \right)\exp ( - {\rm{j}}{k_{tz}}{z_n})} \end{array}} \right] \\&= {{{V}}_{t0}}\left[ {\begin{array}{*{20}{c}} {\left( {{\eta _0}{R^{{\rm{TE\text{-} TM}}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ 0 \\ {\left( {{R^{{\rm{TM \text{-} TM}}}}} \right)\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {\exp ( - {\rm{j}}{k_{0z}}{z_0})} \end{array}} \right],\end{split}$
可求得
$ \begin{split}&{R}^{\rm{TE\text{-}TM}}=\frac{1}{{\eta }_{0}}\left(\frac{{v}_{13}{v}_{34}-{v}_{14}{v}_{33}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}\right){\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0})\;, \\ &{T}^{\rm{TE\text{-}TM}}=\frac{1}{{\eta }_{t}}\left({v}_{24}+\frac{{v}_{13}{v}_{34}-{v}_{14}{v}_{33}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{v}_{21}\right.\\&\qquad\qquad\left.+\frac{{v}_{14}{v}_{31}\!-\!{v}_{11}{v}_{34}}{{v}_{11}{v}_{33}\!-\!{v}_{13}{v}_{31}}{v}_{23}\!\right)\!{\rm{exp}}\left[{\rm{j}}\left({k}_{tz}{z}_{n}\!-\!{k}_{0z}{z}_{0}\right)\right], \\ &{R}^{\rm{TM\text{-}TM}}=\left(\frac{{v}_{14}{v}_{31}-{v}_{11}{v}_{34}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}\right){\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0})\;, \\ &{T}^{\rm{TM\text{-}TM}}=\left({v}_{44}+\frac{{v}_{14}{v}_{31}-{v}_{11}{v}_{34}}{{v}_{11}{v}_{33}-{v}_{13}{v}_{31}}{v}_{43}\right.\\&\qquad\qquad\left.+\frac{{v}_{13}{v}_{34}\!-\!{v}_{14}{v}_{33}}{{v}_{11}{v}_{33}\!-\!{v}_{13}{v}_{31}}{v}_{41}\!\right)\!{\rm{exp}}\left[{\rm{j}}\left({k}_{tz}{z}_{n}\!-\!{k}_{0z}{z}_{0}\right)\right], \end{split}$
其中${v_{ij}}$是${{{V}}_{t0}}$的元素. 当没有石墨烯界面时, 将(33)式代入(36)式得到
$ \begin{split}&{R}^{\rm{TM\text{-}TM}}=-\frac{{v}_{34}}{{v}_{33}}{\rm{exp}}(-2{\rm{j}}{k}_{0z}{z}_{0}),~~ {R}^{\rm{TE\text{-}TM}}=0\;, \\ &{T}^{\rm{TM\text{-}TM}}=\left({v}_{44}-\frac{{v}_{34}}{{v}_{33}}{v}_{43}\right){\rm{exp}}\left[{\rm{j}}\left({k}_{tz}{z}_{n}-{k}_{0z}{z}_{0}\right)\right],\\&{T}^{\rm{TE\text{-}TM}}=0\;,\\[-10pt] \end{split}$
即为TM波入射一般分层介质反透射计算式[16,17].
此外, 还可以证明(32)式中第三式和(36)式中第二式是相等的(限于篇幅, 证明略), 即含各向异性石墨烯界面的各向同性层状介质对入射TE波和TM波的交叉极化反射系数相等, ${R^{{\rm{TM}} \text{-} {\rm{TE}}}} = {R^{{\rm{TE \text{-}TM}}}}$.
4.算　例
以下算例皆考虑室温条件下$T = 300\;{\rm{K}}$以及石墨烯电子弛豫时间为$\tau = 3\;{\rm{ps}}$.
例 1　电、磁偏置下石墨烯薄层的表面电导率及其屏蔽效率、透射波极化旋转
图4(a)是偏置电场${E_{{\rm{bias}}}} = 0.5$V/nm, 频率$f = 1\;{\rm{GHz}}$时, 石墨烯电导率张量元素${\sigma _{\rm{D}}}, {\sigma _{\rm{O}}}$的实虚部和偏置磁场${B_{{\rm{bias}}}}$($- 0.2\text{—} + 0.2\;{\rm{T}}$)的关系. 图4(a)中清楚体现了${\sigma _{\rm{D}}}, {\sigma _{\rm{O}}}$分别为偏置磁场${B_{{\rm{bias}}}}$的偶函数和奇函数.
图 4 石墨烯电导率张量元素及屏蔽效率随偏置磁场变化　(a) 电导率张量元素; (b) 屏蔽效率
Figure4. Diagonal and off-diagonal components of the graphene conductivity tensor and SE as a function of the applied magnetostatic bias: (a) Components of the tensor; (b) SE.

采用本文给出的传播矩阵, 计算了$f = 1\;{\rm{GHz}}$时石墨烯层的反透射并得到其屏蔽效率[1] (shielding effectiveness, SE)如图4(b)所示. 图4(b)中给出偏置电场${E_{{\rm{bias}}}}$分别为$0, \;0.5, \;2, \;5\;{\rm V/nm}$时, SE与偏置磁场${B_{{\rm{bias}}}}$($- 0.2 \text{—} + 0.2\;{\rm{T}}$)的关系.
当偏置电磁场都不为0时, 石墨烯呈各向异性, 即使入射波为线极化, 在反射和透射波中也会产生交叉极化分量, 从而导致反射、透射波成为椭圆极化波. 图5给出了$f = 1\;{\rm{GHz}}$线极化波垂直入射石墨烯层时, 透射椭圆极化波的极化状态参量$\tau, \;\varepsilon $(在庞加莱球上的参数)[23]与偏置电磁场的关系, 其中图5(a)和图5(b)子图分别是偏置电场${E_{{\rm{bias}}}} = 1\;{\rm{V}}/{\rm{n}}{{\rm{m}}}$、偏置磁场${B_{{\rm{bias}}}} = 0.5\;{\rm{T}}$时极化状态参量与偏置磁场${B_{{\rm{bias}}}}$($- 0.2\text{—} + 0.2\;{\rm{T}}$)、偏置电场${E_{{\rm{bias}}}}$($- 4\text{—} + 4\;{\rm V/nm}$)的关系.
图 5 透射波极化状态参量$\tau, \varepsilon $与偏置磁场、偏置电场关系　(a) 偏置电场${E_{{\rm{bias}}}} = 1\;{\rm{V}} \cdot {\rm{n}}{{\rm{m}}^{{{ - }}1}}$; (b) 偏置磁场${B_{{\rm{bias}}}} = 0.5\;{\rm{T}}$
Figure5. Angles $\tau, \;\varepsilon $ describing the polarization state of the transmitted wave as functions of the magnetostatic bias and the electrostatic bias: (a) ${E_{{\rm{bias}}}} = 1\;{\rm{V}} \cdot {\rm{n}}{{\rm{m}}^{{{ - }}1}}$; (b) ${B_{{\rm{bias}}}} = 0.5\;{\rm{T}}$.

该算例为验证算例, 采用和文献[1]相同参数, 计算结果也和文献一致, 表明本文推导正确, 方法有效.
例 2　电、磁偏置下单石墨烯层的反射、透射特性
应用本文方法在偏置电、磁场为${E_{{\rm{bias}}}} = 0.5\;{\rm V/nm}$, ${B_{{\rm{bias}}}} = 0.1\;{\rm{T}}$条件下, 计算了入射波频率$f = 0-1000\;{\rm{GHz}}$频段内石墨烯的反射、透射系数, 如图6所示. 图6(a)为同极化反、透射系数模值, 随着入射角增大, TE波的同极化反射(黑色线图例)会增强而透射(灰色线图例)减弱, 而TM波则刚好相反(黑色、灰色散点图例); 图6(b)为交叉极化反、透射系数模值, 如前所述((27)式), 石墨烯薄层的交叉极化反射系数和交叉极化透射系数的模值都是相等的, 并由图知, 它们随着入射角增大而减小; 图6(c)给出了TE波入射时交叉极化反、透射系数的相位, 和前面的理论一致, 它们相差$\pi $; 图6(d)是将两种透射系数整合并按照定义计算得到的石墨烯的屏蔽效率SE, 可见对TE波的屏蔽效率要高于对TM波的, 而且入射角越大, 对TE波的屏蔽效率越高, TM波则相反.
图 6 各向异性石墨烯界面的反透射及屏蔽效率随频率变化关系　(a) 同极化反透射系数模值; (b) 交叉极化反透射系数模值; (c) 交叉极化反透射系数相位; (d) 屏蔽效率
Figure6. Reflection/transmission and the SE of an isotropic graphene sheet as functions of frequency: (a) Modulus of co-polarization reflective and transmittance coefficients; (b) modulus of cross-polarization reflective and transmittance coefficients; (c) phase of cross-polarization reflective and transmittance coefficients; (d) SE.

例 3　含各向异性石墨烯界面的分层介质的反射与透射
设有一Si/SiO2周期层状结构如图7所示, 第一层为厚度${d_s}$的SiO2层, 第2—10层交替为厚度${d_{\rm{a}}}$和${d_{\rm{b}}}$的Si和SiO2层, 它们的相对复介电系数${\varepsilon _{{\rm{ra}}}}$和${\varepsilon _{r{\rm{b}}}}$分别为
图 7 Si/SiO2周期层状结构
Figure7. Si/SiO2 periodic layered structure.

$ \begin{split}&{\varepsilon }_{\rm{ra}}=\varepsilon ^{\prime }_{\rm{ra}}+{\rm{j}}{\varepsilon^{\prime\prime}_{\rm{ra}}},\;\;{\varepsilon }^{\prime }_{\rm{ra}}=12,\\& \varepsilon^{\prime\prime}_{\rm{ra}} =\frac{\sigma }{\omega {\varepsilon }_{0}},\;\sigma =2.0\;{\rm{S}}/{\rm{m}}, \\ &{\varepsilon }_{\rm{rb}}={{\varepsilon }^{\prime }}_{\rm{rb}}+{\rm{j}}\varepsilon''_{\rm{rb}},\;\;\varepsilon^{\prime }_{\rm rb} = 4.5,\;\;\varepsilon''_{\rm{rb}}=0.07\;{.}\end{split}$
介质层厚度分别为${d_{\rm{a}}} = 0.25\dfrac{{\lambda _0}} {\sqrt {{{\varepsilon '}_{{\rm{ra}}}}} }, {d_{\rm{b}}} = 0.25\dfrac{{\lambda _0}} {\sqrt {{{\varepsilon '}_{{\rm{rb}}}}} }, {d_{\rm{s}}} = 2{d_{\rm{b}}}$, 其中${\lambda _0} = 0.5\;{\rm{mm}}$.
图8中给出当石墨烯(偏置电、磁场为$E_{\rm bias} = 0.5~\rm V/nm$, ${B_{{\rm{bias}}}} = 0.1\;{\rm{T}}$)位于该层状结构表面、正中间(第5和6层界面处)和底层时对垂直入射平面波的反射和透射. 由于石墨烯的各向异性, 反透射场中出现交叉极化成分, 但相对较弱; 作为对比, 图中用灰线给出在没有石墨烯时该层状结构的反透射, 可以看出, 该结构在大约500—700 GHz 频段有一禁带, 石墨烯界面的加入会使禁带向高频方向偏移, 同时会增强非禁带频段的反射, 削弱非禁带频段的透射; 此外, 在禁带频段, 反透射场中的交叉极化成分几乎为0.
图 8 各向异性石墨烯界面位于不同位置时层状结构对垂直入射波的反透射随频率变化关系　(a) 反射系数模值; (b) 透射系数模值
Figure8. Reflection/transmission of normally incident waves from a periodic layered structure with an isotropic graphene at different interface as functions of frequency: (a) Modulus of reflective coefficient; (b) modulus of transmittance coefficient.

图9是石墨烯位于该层状结构表面时对$\theta = {60^ \circ }$斜入射线极化平面波的反射和透射. 可以明显地看到, TE波的同极化反射要强于TM波的, 而同极化透射比TM波弱; 同时, 两种交叉极化反射系数相同, 和前面理论讨论一致.
图 9 各向异性石墨烯界面位于上表面时层状结构对斜入射波的反透射随频率变化关系　(a) 反射系数模值; (b) 透射系数模值
Figure9. Reflection/transmission of oblique incident waves from a periodic layered structure with an isotropic graphene surface as functions of frequency: (a) Modulus of reflective coefficient; (b) modulus of transmittance coefficient.

5.结　论
基于边界条件和相位匹配推导出在电、磁偏置情形下具有各向异性特征的石墨烯界面的传播矩阵, 并且进一步将其嵌入传统各向同性层状介质传播矩阵序列, 使其能够解析计算含石墨烯界面的分层介质对平面电磁波的反透射. 偏置静电、静磁场影响调节石墨烯表面电导率张量, 从而实现对石墨烯反透射特性及屏蔽性能的控制, 并可能使用偏压的石墨烯作为极化旋转器. 本文方法同样适用于石墨烯为各向同性时情形(偏置静电或/和静磁场为零时), 其可看作一种特例. 此外, 在太赫兹及以上更高频率, 需要考虑空间色散效应引起的各向异性, 此时可采用文中相似的处理方法进行理论上的推导及进一步分析讨论.
相关话题/传播 结构 电磁场 石墨 计算
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • 基于夹层结构的偏振无关1×2定向耦合型解复用器的设计
    摘要:提出一种基于夹层结构的偏振无关1×2定向耦合型解复用器,用于分离1310nm和1550nm两个波长.通过合理选择夹层结构芯区的折射率及波导间隙,可以调节同一波长两个正交偏振模的耦合长度相等,实现偏振无关;通过合理选择夹层结构波导宽度,可以使两个波长分别从不同输出波导端口输出,实现解复用功能.运 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 具有非广延分布电子的碰撞等离子体磁鞘的结构
    摘要:很多关于等离子体鞘层的研究工作都是基于电子满足经典的麦克斯韦速度分布函数,而等离子体中的粒子具有长程电磁相互作用,使用Tsallis提出的非广延分布来描述电子更为恰当.本文建立一个具有非广延分布电子的碰撞等离子体磁鞘模型,理论推导出受非广延参数q影响的玻姆判据,离子马赫数的下限数值会随着参数q ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 多狄拉克锥的二维CrPSe<sub>3</sub>的半金属铁磁性与电子结构<span style="color:#ff0000;">（已撤稿）&
    摘要:采用第一性原理方法和紧束缚模型研究了二维层状结构过渡金属硫族磷化物CrPSe3的磁学特性和电子结构.二维CrPSe3晶格呈现六角蜂窝状结构,声子谱证明其具备较好的动力学稳定性.二维CrPSe3的电子结构呈现出自旋无能隙磁性半金属特性.蒙特卡罗方法模拟出二维CrPSe3的居里温度为224K,分子 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于人工表面等离激元结构的超表面磁镜
    摘要:提出了一种表面粗糙磁镜的概念,该界面由人工表面等离激元结构阵列设计而成.这种人工表面等离激元结构通过周期性地将螺旋金属条插入到介电圆盘构造中以支持强磁偶极共振模式.特别地,对于不同外半径下的螺旋结构,可以通过调节每个结构的螺旋度以支持相同共振频率的磁偶极模式.为此,设计了由多种不同尺寸的人工表 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 双螺线圈射频共振结构增强硅空位自旋传感灵敏度方法
    摘要:针对硅空位自旋磁共振信号射频场非均匀展宽问题,提出并设计了一种双螺线圈射频共振结构,利用双螺线圈平行对称特性,构建射频场均匀区,非均匀性小于0.9%,相比单根直线性结构,均匀性提高了56.889倍.同时,利用射频信号近距离互感耦合共振特性,实现了射频场的增强,相比单线圈结构增强了1.587倍, ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 涂覆石墨烯的非对称椭圆电介质纳米并行线的模式分析
    摘要:设计了一种涂覆石墨烯的非对称椭圆电介质纳米并行线波导.在椭圆柱坐标系中,借助于Mathieu函数和坐标变换,采用多极方法对波导所支持的6个最低阶模式进行了研究,并分析了这些模式的特性与工作波长、石墨烯费米能以及波导结构参数之间的依赖关系.结果表明,调节波导的工作波长、石墨烯的费米能及纳米线之间 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 全固态电池中界面的结构演化和物质输运
    摘要:全固态电池中科学问题的本质在于引入的固态电解质的特性及全新的固-固界面的存在.从构-效关系出发,固-固界面和电解质自身的结构演化与物质输运过程决定了全固态电池的性能.随着固态电解质材料研究的不断丰富,目前全固态电池中的问题主要集中在固-固界面，界面处的组成和结构限制了全固态电池的性能.根据固- ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于分子动力学模拟的铜晶面石墨烯沉积生长机理
    摘要:化学气相沉积法是大面积、高质量石墨烯沉积制备实践中的重要方法.本文采用分子动力学仿真技术,模拟了利用化学气相沉积法在铜(111)晶面制备石墨烯的过程,研究揭示了石墨烯在铜(111)晶面上的微观生长机理.研究结果表明:石墨烯的沉积生长可描述为第一阶段的二元碳、三元碳和碳链形成阶段,以及第二阶段的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 新型铁电拓扑结构的构筑及其亚埃尺度结构特性
    摘要:铁电拓扑结构因其尺寸小而且具有优良的物理特性,有望应用于未来高性能电子器件中.本文从应变、屏蔽和外场等对于铁电材料至关重要的几个外部要素出发,结合薄膜厚度等材料内部参数,针对PbTiO3和BiFeO3这两种典型的铁电材料,简要总结新型铁电拓扑结构的形成及其在外场作用下的演变规律.利用具有亚埃尺 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 铁电纳米结构中奇异极化拓扑畴的研究新进展
    摘要:铁电体中极化拓扑畴(如涡旋畴)有望带来一系列新颖物理现象、新性能和新应用前景(如存储器件应用),从而引起了广泛兴趣.尤其是近年来在铁电纳米结构中发现了一系列有趣的新奇极化拓扑畴态,例如涡旋、中心畴、斯格明子、麦韧(Meron,也有称半子)等,引发了新一轮探索热潮.这些发现为进一步探索其中蕴含的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.