Processing math: 10%

Search

Article

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Study of magnetoplasmons in graphene rings with two-dimensional finite element method

Wang Wei-Hua

Wang Wei-Hua. Study of magnetoplasmons in graphene rings with two-dimensional finite element method. Acta Phys. Sin., 2023, 72(8): 087301. doi: 10.7498/aps.72.20222467
Citation: Wang Wei-Hua. Study of magnetoplasmons in graphene rings with two-dimensional finite element method. Acta Phys. Sin., 2023, 72(8): 087301. doi: 10.7498/aps.72.20222467

Study of magnetoplasmons in graphene rings with two-dimensional finite element method

Wang Wei-Hua
Article Text (iFLYTEK Translation)
PDF
HTML
Get Citation
  • Graphene plasmons are important collective excitations in graphene, which play a key role in determining the optical properties of graphene. They have quite lots of unique features in comparison with classical plasmons in noble metals. Of them, the active tunability is the most attractive, which is realized by external gating (equivalently electric field). As is well known, graphene also has strong magnetic response (e.g. room temperature quantum Hall effect), so magnetic field can act as another degree of freedom for actively tuning graphene plasmons, with the new quasi particles being so-called graphene magneto-plasmons. Because of the two-dimensional nature of graphene, the numerical studies (or full wave simulations) of graphene magneto-plasmons are usually carried out through a three-dimensional approximation, e.g. treating two-dimensional graphene as a very thin three-dimensional film. Actually, this treatment takes quite some time and requires high memory consumption. Herein, starting from Coulomb law and charge conservation law, we propose an alternative numerical method, namely, two-dimensional finite element method, to solve this problem. All the calculations are now performed in two-dimensional graphene plane, and the usual three-dimensional approximation is not required. To characterize the excitations of graphene magneto-plasmons, the eigenvalue loss spectrum is introduced. Based on this method, graphene magneto-plasmons in graphene rings of four kinds are investigated. The strongest magneto-optic effect is observed in circular ring, which is consistent with its highest rotational symmetry. In all the rings, the lowest dipolar graphene magneto-plasmon always supports symmetric mode splitting, which can be further modified by the interaction between inner edge and outer edge of ring. As the hole size is very small, the edge current confined to the outer edge dominates, and that confined to the inner edge can be ignored; while increasing the hole size, the interaction between these two edges increases, which results in the reduction of the symmetric mode splitting; when the hole size is larger than a critical value, the symmetric mode splitting will disappear.
      PACS:
      73.20.Mf(Collective excitations (including excitons, polarons, plasmons and other charge-density excitations))
      78.20.Ls(Magneto-optical effects)
      78.67.Wj(Optical properties of graphene)
      71.45.Gm(Exchange, correlation, dielectric and magnetic response functions, plasmons)
      Corresponding author: Wang Wei-Hua, wh.wang@outlook.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 12174440).

    等离激元源于自由电子的集体振荡, 是固体中的一种元激发, 具有很多独特的光学性质, 如亚波长约束、局域场增强与突破衍射极限[1-3], 可极大地提升光与物质的相互作用; 作为兼具电子与光子优势的新型载体, 其在物理、化学、材料和生物等领域得到了高度的关注, 广泛应用于光电探测[4-6]、生物传感[7-9]、高分辨成像[10-12]、催化反应[13-15]等. 与其他准粒子(如声子)类似, 振动频率与弛豫时间是等离激元的两个重要指标, 当下的研究多集中于贵金属体系(如Au和Ag), 前者可通过几何结构、尺寸来调节, 以满足不同工作频率的需求, 而后者则严重受制于固有的欧姆损耗[16,17], 相干时间与传输距离较短, 成为其在器件应用方面的主要障碍. 要从根本上解决这一问题, 必须挖掘更好的等离激元材料[18,19], 以替代传统贵金属. 人们尝试了Al[20], Na[21], Cu[22], 以及掺杂的半导体[23]等, 情况有所改善, 但更多有价值的探索仍有待进一步开展.

    2004年, 二维(two-dimensional, 2D)材料石墨烯首次由Novoselov等[24]通过机械剥离方法成功制备, 展现出了优异的力学[25-27]、热学[28-30]、光学[31-34]及电学[35-38]特性, 迅即成为理论与应用研究的焦点, 更是掀起了2D材料研究的热潮[39-41]. 作为一种零带隙的半金属, 中性石墨烯不支持等离激元共振, 但当其费米面偏离狄拉克(Dirac)点时[36,37], 便提供了该准粒子激发所需的内在条件[42-45], 且较贵金属有一系列独特的优势: 1) 电磁场高度局域于单原子层的厚度处[46-48]; 2) 共振频率可通过外置电压动态调控[49,50]; 3) 欧姆损耗较小(缘于极高的载流子迁移率) [51-53]; 4) 工作在稀缺的太赫兹及红外波段[54-57]. 相关的工作众多, 汇聚成石墨烯等离激元学(graphene plasmonics)这一新兴的学科方向, 也极大地推动了2D材料光学的发展[58-60]. 毫无疑问, 这里外置电压(即电场)赋予的动态可调性打破了传统等离激元仅可依靠结构调控的桎梏, 是石墨烯等离激元最引人注目的亮点, 为多频段、多功能、可切换光电器件的设计提供了潜在的解决方案. 然而, 该体系中是否还蕴藏着其他调控自由度?

    值得注意的是, 与传统的二维电子气(two-dimensional electron gases, 2DEGs)相比, 线性能带结构的石墨烯具有更强的磁场响应, 很多极限条件下的物理效应在常规条件下便可以观测, 如量子霍尔效应[61-63]、分数量子霍尔效应[64-66]和量子自旋霍尔效应[67,68]等; 同时, 体系的光与物质相互作用也受磁场强烈地调制, 即磁光效应(magneto-optic effect, MOE), 比如显著增强的法拉第旋转[69-72]. 鉴于此, 电子的集体振荡以单粒子回旋运动为基础, 其准粒子俗称磁等离激元[73-75], 是MOE与等离激元优势的结合体[76], 而贵金属不支持MOE, 因此和电场一样, 磁场也是石墨烯等离激元独特的调控自由度. 有趣的是, 这两个自由度的调控机理差异很大: 电场改变的是电子的浓度, 而其单粒子的性质不变; 磁场改变的是单粒子的运动状态, 而电子的浓度不变; 两种方式并行不悖, 可灵活搭配以实现最优的等离激元调控[77].

    当前, 石墨烯磁等离激元(graphene magneto-plasmon, GMP)的研究多开展于较为简单的石墨烯体系, 如无限大[78-82]、半无限大[83-85]和圆盘结构[86,87]等. 此外, 磁场的存在打破了时间反演对称性, GMP很自然地与非互易[88,89]、拓扑[90-92]等物理现象密切相关. 考虑到不可避免且尤为重要的界面效应, 我们曾探讨了石墨烯半无限薄片[93]、圆盘[94]、正多边形盘[95]、圆环[96,97]中GMP的模式劈裂与杂化等基本物性, 但对界面效应仍缺乏较为全面的认识, 更多深入的探究诸如内外界面的耦合、GMP高阶态的激发等仍有待开展.

    鉴于GMP研究大多采用三维(three-dimensional, 3D)近似, 即将2D石墨烯等效成厚度很薄的3D块材, 仿真计算需消耗大量资源. 本文提出一种可行的替代方案, 不再需要3D近似, 所有的计算只限于2D石墨烯面内, 极大地节省了资源、提升了计算效率. 利用该2D计算方法, 讨论了4类石墨烯环中GMP的激发, 分析了内外边界的相互作用对GMP模式劈裂的影响.

    媒质的光学响应由其微观的单粒子跃迁(single-particle transition, SPT)决定, 但可通过宏观的电磁参量如介电函数或电导率表征. MOE体系中SPT发生在不同的朗道能级(Landau level)之间[98-102], 传统抛物线形的2DEG中第n个朗道能级的能量为En=(n+1/2)ωc, 这里的回旋共振频率(cyclotron resonance frequency, CRF) ωc=eB/m, m为电子的有效质量; 而石墨烯中第n个朗道能级的能量可表示为En=±nωc, 这里ωc=2vF/lB, 费米速度vF=106m/s, 磁长度lB=/(eB), 即ωcB, 而不是像传统2DEG中的线性关系; 但在n很大时(对应高掺杂情况), 石墨烯和2DEG中由于朗道能级不同而导致的SPT差异被逐渐抹平, 石墨烯的磁光面电导率也可用经典的德鲁德(Drude)模型表述[94]:

    \begin{split} \bar {\bar{{\boldsymbol{\sigma}} }}=\;&\left( \begin{array}{cc}{\sigma }_{xx}& {\sigma }_{xy}\\ -{\sigma }_{xy}& {\sigma }_{yy}\end{array} \right)\\=\;&\dfrac{{e}^{2}{E}_{{\rm{F}}}}{{\rm{\pi }}{\hslash }^{2}\left[{\left(\omega +{\rm{i}}{\tau }^{-1}\right)}^{2}-{\omega }_{{\rm{c}}}^{2}\right]}\\ &\times\left( \begin{array}{cc}{\rm{i}}\left(\omega +{\rm{i}}{\tau }^{-1}\right)& {\omega }_{{\rm{c}}}\\ -{\omega }_{c}& {\rm{i}}\left(\omega +{\rm{i}}{\tau }^{-1}\right)\end{array} \right)\text{. } \end{split} (1)

    这里 {\omega }_{{\rm{c}}}=eB/{m}^{*}=eB{v}_{{\rm{F}}}^{2}/{E}_{{\rm{F}}} 为等效的CRF, 而不是石墨烯体系严格的CRF, \tau =\mu {E}_{{\rm{F}}}/\left(e{v}_{{\rm{F}}}^{2}\right) 为弛豫时间, 是表征损耗的唯象参数, 系数 \mu 是载流子迁移率. 本文中将这些参数设定为常数(除非单独说明), 如磁场 B=5 T, 费米能级 {E}_{{\rm{F}}}=0.6 eV, 以及载流子迁移率μ = 10000 cm2·V–1·s–1.

    基于上述磁光面电导率, 体系的GMP可通过麦克斯韦方程组求解, 一般处理的方法首先是将2D石墨烯近似成非常薄的3D块材[94-97], 比如其厚度是面内尺寸的1%甚至更小, 然后在任意电磁仿真软件中进行建模. 但是这样处理往往有一定的代价, 就是人为引入的这个极小厚度会使剖分的3D网格数量急剧增加, 耗费很多的计算资源与时间. 本文在准静态近似下, 以库仑定律和电荷守恒定律为基础, 构建自洽的2D面积分方程, 并在面内剖分2D网格进行求解. 此前, 我们曾用该方法计算无磁场情况下石墨烯及双层石墨烯体系的等离激元[103-105], 而这里GMP的计算不能简单照搬先前的方案, 需要考虑面电导率中非对角项的贡献. 自洽的2D面积分方程为

    \phi \left(\boldsymbol{r}\right)={\phi }^{{\rm{e}}{\rm{x}}{\rm{t}}}\left(\boldsymbol{r}\right)+\frac{1}{4{\rm{\pi }}{\varepsilon }_{{\rm{s}}}}\int \frac{\rho \left({\boldsymbol{r}}'\right)}{\left|\boldsymbol{r}-{\boldsymbol{r}}'\right|}{\rm{d}}{\boldsymbol{r}}', \tag{2a}
    \rho \left(\boldsymbol{r}\right)=\frac{{\rm{i}}{\sigma }_{xx}}{\omega }{\nabla }^{2}\phi \left(\boldsymbol{r}\right). \tag{2b}

    这里的积分和微分运算都只在石墨烯面内, 例如 \nabla ={\partial }_{x}\widehat{\boldsymbol{x}}+{\partial }_{y}\widehat{\boldsymbol{y}} , {\varepsilon }_{{\rm{s}}}是周围环境的平均介质常数, 如石墨烯两侧介质的介电常数分别为 {\varepsilon }_{1} {\varepsilon }_{2} , 则{\varepsilon }_{{\rm{s}}}= \left({\varepsilon }_{1}+{\varepsilon }_{2}\right)/2. (2a)式是库仑定律的积分形式, 这一点很直观, 需要说明的是, 由这个方程可以求得空间任意位置的电势 \phi \left(\boldsymbol{r}\right) , 而不限于石墨烯面内, 但电荷 \rho \left(\boldsymbol{r}\right) 只分布于石墨烯面内, 因而(2b)式中的电势 \phi \left(\boldsymbol{r}\right) 只限于石墨烯面内, 即石墨烯面内的 \phi \left(\boldsymbol{r}\right) \rho \left(\boldsymbol{r}\right) 可通过这个方程组自洽求解. (2b)式来自电荷守恒定律, 但为什么其只与面电导率的对角项有关. 解开这个疑团, 需先从石墨烯的面电流入手:

    \begin{split} \boldsymbol{J}=\;&{J}_{x}\widehat{\boldsymbol{x}}+{J}_{y}\widehat{\boldsymbol{y}}=\left({\sigma }_{xx}{E}_{x}+{\sigma }_{xy}{E}_{y}\right)\widehat{\boldsymbol{x}} \\ & +\left(-{\sigma }_{xy}{E}_{x}+{\sigma }_{xx}{E}_{y}\right)\widehat{\boldsymbol{y}}\\ =\;&-{\sigma }_{xx}\nabla \phi -{\sigma }_{xy}\nabla \times \left[\phi \left(\boldsymbol{r}\right)\widehat{\boldsymbol{z}}\right]{\rm{. }} \end{split} (3)

    这里已经利用了 {\sigma }_{xx}={\sigma }_{yy} ; 由此, 不难看出 -{\sigma }_{xy}\nabla \times \left[\phi \left(\boldsymbol{r}\right)\widehat{\boldsymbol{z}}\right] 是一个矢量函数的旋度, 其散度恰好为零, 即如(2b)式所示, \nabla \cdot \boldsymbol{J} 只由 -{\sigma }_{xx}\nabla \phi 贡献.

    对于任意形状的石墨烯结构, (2a)式和(2b)式需要离散化求解, 这里沿用之前发展的2D有限元方法(2D finite element method, 2DFEM) [104], 以每个格点上的 \phi \left(\boldsymbol{r}\right) \rho \left(\boldsymbol{r}\right) 作为基函数. (2a)式可直接给出关于 \phi \left(\boldsymbol{r}\right) \rho \left(\boldsymbol{r}\right) 之间关系的第一个矩阵方程, 比如

    \left(\begin{array}{c}{\phi }_{1}\\ \vdots\\ {\phi }_{N}\end{array}\right)=\left(\begin{array}{c}{\phi }_{1}^{{\rm{e}}{\rm{x}}{\rm{t}}}\\ \vdots\\ {\phi }_{N}^{{\rm{e}}{\rm{x}}{\rm{t}}}\end{array}\right)+{{{\boldsymbol{M}}}}\left(\begin{array}{c}{\rho }_{1}\\ \vdots\\ {\rho }_{N}\end{array}\right)\text{, } (4)

    N 为格点的数目. 问题的关键在于(2b)式的处理, 需将其转化为弱形式, 即在其两边乘以试探函数 \phi \left(\boldsymbol{r}\right) 再积分,

    \displaystyle\int \phi \left(\boldsymbol{r}\right)\rho \left(\boldsymbol{r}\right){\rm{d}}\boldsymbol{r}=\dfrac{{\rm{i}}{\sigma }_{xx}}{\omega }\displaystyle\int \phi \left(\boldsymbol{r}\right){\nabla }^{2}\phi \left(\boldsymbol{r}\right){\rm{d}}\boldsymbol{r} .

    利用恒等式\nabla \cdot \left(f\boldsymbol{g}\right)=\nabla f\cdot \boldsymbol{g}+f\nabla \cdot \boldsymbol{g}, 等式右边可进行如下约化:

    \begin{split} &\displaystyle\int \nabla \cdot \left[\phi \left(\boldsymbol{r}\right){\sigma }_{xx}\nabla \phi \left(\boldsymbol{r}\right)\right]{\rm{d}}\boldsymbol{r} - {\sigma }_{xx}\displaystyle\int \nabla \phi \left(\boldsymbol{r}\right) \cdot \nabla \phi \left(\boldsymbol{r}\right)\\ =\;& \displaystyle\int {{\boldsymbol{n}}} \cdot \left[\phi \left(\boldsymbol{r}\right){\sigma }_{xx}\nabla \phi \left(\boldsymbol{r}\right)\right]{\rm{d}}\boldsymbol{l} - {\sigma }_{xx}\displaystyle\int \nabla \phi \left(\boldsymbol{r}\right) \cdot \nabla \phi \left(\boldsymbol{r}\right){\rm{d}}\boldsymbol{r}. \end{split} (5)

    等式右边第一项是沿着边界的线积分, n为边界上每处的外法线方向. 要完全求解这个方程必须考虑合理的边界条件, 按照惯例可使法向电流为0,{\boldsymbol{n}}\cdot{\boldsymbol{J}} = 0, 即对于电势 \phi \left(\boldsymbol{r}\right) 的诺伊曼边界条件(Neumann boundary condition). 在无磁场的情况下, 由于 \boldsymbol{J}= -{\sigma }_{xx}\nabla \phi \left(\boldsymbol{r}\right) 第一项线积分正好为0; 但有磁场时, 结合(3)式容易得到

    \begin{split} & \displaystyle \int \phi \left(\boldsymbol{r}\right)\rho \left(\boldsymbol{r}\right){\rm{d}}\boldsymbol{r} \\ =\;&-\dfrac{{\rm{i}}{\sigma }_{xx}}{\omega }\displaystyle \int \nabla \phi \left(\boldsymbol{r}\right)\cdot \nabla \phi \left(\boldsymbol{r}\right){\rm{d}}\boldsymbol{r}\\ &-\dfrac{{\rm{i}}{\sigma }_{xy}}{\omega }\displaystyle \int \phi \left(\boldsymbol{r}\right)\boldsymbol{n}\cdot \nabla \times \left[\phi \left(\boldsymbol{r}\right)\widehat{\boldsymbol{z}}\right]{\rm{d}}\boldsymbol{l}. \end{split} (6)

    该方程离散化的最终形式为

    \begin{split} &\left(\begin{array}{ccc}{\phi }_{1}& \cdots & {\phi }_{N}\end{array}\right)\left(\begin{array}{c}{\rho }_{1}\\ \vdots\\ {\rho }_{N}\end{array}\right)=\left(\begin{array}{ccc}{\phi }_{1}& \cdots & {\phi }_{N}\end{array}\right)\\ &\times\left[-\frac{{\rm{i}}{\sigma }_{xx}}{\omega }\boldsymbol{P}-\frac{{\rm{i}}{\sigma }_{xy}}{\omega }\boldsymbol{Q}\right]\left(\begin{array}{c}{\phi }_{1}\\ \vdots\\ {\phi }_{N}\end{array}\right) \text{, } \end{split}

    \left(\begin{array}{c}{\rho }_{1}\\ \vdots\\ {\rho }_{N}\end{array}\right)=\left[-\dfrac{{\rm{i}}{\sigma }_{xx}}{\omega }\boldsymbol{P}-\dfrac{{\rm{i}}{\sigma }_{xy}}{\omega }\boldsymbol{Q}\right]\left(\begin{array}{c}{\phi }_{1}\\ \vdots\\ {\phi }_{N}\end{array}\right). (7)

    (4)式和(7)式构成格点上电势 \phi \left(\boldsymbol{r}\right) 与电荷 \rho \left(\boldsymbol{r}\right) 的一组自洽方程, 将(7)式代入(4)式可得

    \left[1+\dfrac{{\rm{i}}{\sigma }_{xx}}{\omega }\boldsymbol{M}\boldsymbol{P}+\dfrac{{\rm{i}}{\sigma }_{xy}}{\omega }\boldsymbol{M}\boldsymbol{Q}\right]\left(\begin{array}{c}{\phi }_{1}\\ \vdots\\ {\phi }_{N}\end{array}\right)=\left(\begin{array}{c}{\phi }_{1}^{{\rm{e}}{\rm{x}}{\rm{t}}}\\ \vdots\\ {\phi }_{N}^{{\rm{e}}{\rm{x}}{\rm{t}}}\end{array}\right). (8)

    很显然, 方程左边括号里的部分对应于体系等效的介电函数(dielectric function), 找出其为零的解, 即给出体系自持(self-sustaining)的本征GMP模式, 但考虑到耗散及其矩阵形式, 该方法并不具备可操作性; 这里将使用其本征值损失谱方法, 具体如下: 对任一频率 \omega , 对角化介电函数矩阵, 得到 N 个本征值 \left({\varepsilon }_{1}, \cdots {, \varepsilon }_{N}\right) 及其本征态, 挑出其中实部趋于零且虚部最小的那个本征值(通过损失函数 -{\rm{I}}{\rm{m}}\left({\varepsilon }_{n}^{-1}\right) 的极值), 记为该频率处的谱值, 而在频谱上谱值的峰即对应自持的GMP模式, 如果峰值处对应多个相同的本征值, 则表明该模式为多重简并. 在实施过程中, 矩阵 \boldsymbol{M} , \boldsymbol{P} 以及 \boldsymbol{Q} 只由石墨烯结构决定, 因此只需求解一次, 而频率 \omega 的依赖则完全包含在系数 {\rm{i}}{\sigma }_{xx}/\omega {\rm{i}}{\sigma }_{xy}/\omega 中. 除了本征值损失谱, 该2DFEM还方便计算电磁波散射谱、吸收谱等, 对于入射场 \left({\phi }_{1}^{{\rm{e}}{\rm{x}}{\rm{t}}}, \cdots {, \phi }_{N}^{{\rm{e}}{\rm{x}}{\rm{t}}}\right) , 求解(8)式可得电势分布 \left({\phi }_{1}, \cdots {, \phi }_{N}\right) . 结合(7)式, 求得激发电荷分布 \left({\rho }_{1}, \cdots {, \rho }_{N}\right) , 接着计算偶极矩其虚部则对应电磁波吸收; 这里的入射场即每个格点上电势的初始值可随意设置, 如平面波、点源、电子束等, 以满足个性化的计算需求.

    通过上述计算方法, 探讨石墨烯圆环、方环、方孔圆环、圆孔方环中GMP的激发. 不失一般性, 首先考虑环的直径或边长 {D}_{{\rm{R}}}=100 nm, 孔的直径或边长 {D}_{{\rm{H}}}=50 nm, 其结构如图1所示. 借助开源软件Gmsh[106], 采用三角形微元对结构进行离散化, 在此微元基础上开展微分、积分运算的细节, 请参阅此前的工作[104]. 在实际的计算中, 为了保证结果的收敛性和精度, 往往需要较为精细的网格. 值得注意的是, 最佳的方案并非是直接设置很小的网格(即很小的三角形微元), 而是设置一般尺寸的三角形微元, 然后对每一个三角形微元进行拆分细化(refine by splitting), 即将一个三角形微元拆分成多个三角形微元. 这样一来, 每个三角形的顶点(对应(8)式中的格点)会被更多的微元共用, 即通过引入更少的格点产生更多的微元, 前者决定(8)式中介电函数矩阵的维度, 后者决定计算的精度, 如图1(a)所示的圆环, 5312个格点构成10240个微元, 此时计算结果已完全收敛. 该方案可极大地提高计算效率, 广泛应用于有限元计算.

    图 1 石墨烯圆环(a)、方环(b)、方孔圆环(c)、圆孔方环(d)网格划分示意图, 环的直径或边长均为100 nm, 孔的直径或边长均为50 nm\r\nFig. 1. Schematic diagrams of mesh generation for graphene circular ring (a), square ring (b), circular ring with square hole (c), and square ring with circular hole (d). The diameteror side length of ring is 100 nm, and that of hole is 50 nm.
    图 1  石墨烯圆环(a)、方环(b)、方孔圆环(c)、圆孔方环(d)网格划分示意图, 环的直径或边长均为100 nm, 孔的直径或边长均为50 nm
    Fig. 1.  Schematic diagrams of mesh generation for graphene circular ring (a), square ring (b), circular ring with square hole (c), and square ring with circular hole (d). The diameteror side length of ring is 100 nm, and that of hole is 50 nm.

    为了检验算法的可靠性, 首先对比2DFEM与常见3D有限元(3DFEM)的计算结果. 简单起见, 考查石墨烯圆盘中GMP激发, 即图1(a)中孔直径为零的情形. 至于3DFEM, 借助全波仿真软件COMSOL来开展. 构建球形计算区域, 通过完美匹配层(perfectly matched layer, PML)截断, 石墨烯圆盘的厚度设置为0.5 nm, 置于区域中心, 线极化电磁波正入射, 消光截面 {\sigma }_{{\rm{e}}{\rm{x}}{\rm{t}}} 用于表征体系光学响应(详情参阅此前工作[95-97]). 这里, {\sigma }_{{\rm{e}}{\rm{x}}{\rm{t}}} 可借助光学定理求解, 将其除以圆盘的面积 S , 即 {\sigma }_{{\rm{e}}{\rm{x}}{\rm{t}}}/S , 就得到了消光系数. 该3DFEM在GMP研究中最为常用, 但只适用于探究含有偶极矩的低阶模式, 比如偶极GMP, 其他模式的 {\sigma }_{{\rm{e}}{\rm{x}}{\rm{t}}} 均接近于0, 俗称暗模式. 鉴于此, 不妨比较偶极GMP的计算结果(图2). 相比于2DFEM的结果, 3DFEM的谱线出现了细微的红移, 但在磁场诱发的对称模式劈裂(symmetric mode splitting, SMS)方面表现一致(如劈裂的大小 {\omega }_{+}-{\omega }_{-} 相同). 图2(b)中劈裂之后的峰高(红色)是蓝色峰高的一半, 源于线极化本征态可分解为等量左/右旋圆极化本征态的线性叠加. 图2(b)中虚线与蓝色曲线相交部分, 给出了磁场 B=0 时模式的半峰宽, 对应模式的损耗, 其由辐射损耗和材料吸收损耗构成. 虽然2DFEM由于采用准静态近似而忽略了辐射损耗, 但两种方法给出的半峰宽几乎一样(包括蓝色曲线与红色曲线的峰), 表明材料吸收损耗主导, 这是等离激元共振体系的特性, 在深亚波长GMP体系则尤为明显.

    图 2 石墨烯圆盘的偶极GMP谱线, 数据分别来自2DFEM (a)和3DFEM (b), 其中虚线标记了蓝色曲线半峰的高度\r\nFig. 2. Dipolar GMP spectra in graphene disk, obtained from 2DFEM (a) and 3DFEM (b), respectively. The dashed lines mark the position of half maximum of blue curves.
    图 2  石墨烯圆盘的偶极GMP谱线, 数据分别来自2DFEM (a)和3DFEM (b), 其中虚线标记了蓝色曲线半峰的高度
    Fig. 2.  Dipolar GMP spectra in graphene disk, obtained from 2DFEM (a) and 3DFEM (b), respectively. The dashed lines mark the position of half maximum of blue curves.

    与3DFEM对比, 2DFEM框架下的本征值损失谱方法不依赖于模式的偶极矩, 可方便探究各阶GMP的行为. 如图3(a)所示, 角动量 l\ne 0 的模式都会产生SMS, 即每组GMP( {\omega }_{-} {\omega }_{+} , 红色曲线的峰值)对称分布于磁场为0时等离激元({\omega }_{{\rm{p}}}, 蓝色曲线的峰值)的两侧, 如 l=1 的偶极模式、 l=2 的四极模式等; l=0 的模式, 俗称轴对称呼吸模式[103], 没有发生SMS, 其共振频率与磁场为0时等离激元的频率基本一致. 这个结果可借助一个简单的图像来理解. 磁场作用于载流子的洛伦兹力宏观表现为沿边界传输的环形电流, 其环向与磁场构成右手定则, 不妨假定为顺时针方向. 对于 l\ne 0 的模式, 其共振电流也沿环向, 且两个简并态 \pm l 对应的环向恰好相反, 因此顺时针方向模式因与边界环流同向而缩短运动周期, 导致频率增大, 但逆时针方向模式因与边界环流反向而延长运动周期, 导致频率减小, 这里的增减量均等于边界环流即SMS. 对于 l=0 的模式, 其共振电流沿径向, 边界环流不会对其产生加速或减速效果, 且该情形下 {\omega }_{{\rm{p}}}\gg {\omega }_{{\rm{c}}} , 边界环流对于快速的径向电流可做准静态近似, 即忽略其对后者的影响. 该图像不局限于圆形, 适用于任意形状体系中GMP激发的描述, 详情参阅此前的工作[95].

    图 3 石墨烯圆盘(a)和方块(b)的本征值损失谱, 圆盘和方块的直径与边长均为100 nm\r\nFig. 3. Eigenvalue loss spectra of graphene disk (a) and graphene square (b). The diameter of disk and the side length of square are both 100 nm.
    图 3  石墨烯圆盘(a)和方块(b)的本征值损失谱, 圆盘和方块的直径与边长均为100 nm
    Fig. 3.  Eigenvalue loss spectra of graphene disk (a) and graphene square (b). The diameter of disk and the side length of square are both 100 nm.

    图3(b)所示, 对于石墨烯方块, 也有5组GMP模式出现在相似的频率范围, 但相比圆盘, 整体发生红移, 且仅有2组伴随着SMS, 被第1条和第5条虚线标注. 从上述分析中可知, 这两个模式类似于圆盘中偶极和四极模式, 而其他3个模式则由径向振动主导. GMP模式上的差异将会影响它们之间的耦合, 圆形或方形作为内边界或外边界时作用也将不同, 因而4类石墨烯环中GMP的值需要进一步深入探究.

    在给定的频率范围内, 数值计算的结果如图4所示, 不难发现4类石墨烯环支持的GMP激发有很大不同, 归纳为4点. 1) 圆环拥有最多的4组GMP经历SMS, 体现最强的磁场效应, 这与其体系最高的对称性一致; 2) 方孔圆环仅有2组GMP发生SMS, 因其面积最小, 更多的模式被推至高频; 3) 最低阶偶极GMP均支持SMS, 这由体系的 {C}_{4} 旋转对称性决定; 4) SMS的大小呈不均匀分布, 而非越低阶GMP的SMS越大, 这由环的内外边界对SMS不同的贡献所致, 在均匀的磁场中, 内外边界处的环流恰好反向, 即内边界的出现会抑制SMS, 而抑制的程度取决于GMP在内外边界处的局域场分布.

    图 4 4类石墨烯环的本征值损失谱, 椭圆标记了发生SMS的GMP模式\r\nFig. 4. Eigenvalue loss spectra of graphene rings of 4 kinds. The ellipses mark the symmetric mode splitting of those GMPs.
    图 4  4类石墨烯环的本征值损失谱, 椭圆标记了发生SMS的GMP模式
    Fig. 4.  Eigenvalue loss spectra of graphene rings of 4 kinds. The ellipses mark the symmetric mode splitting of those GMPs.

    上述讨论中提及, 4类石墨烯环中最低阶的偶极GMP都会发生SMS, 但是鉴于内外边界的反向环流, 其耦合作用将极大地调制GMP激发, 这个过程可通过变化孔的尺寸展示. 如图5所示, 在孔很小时(如10 nm), 4类环都具有最大的SMS, 且圆环与方孔圆环、方环与圆孔方环的 {\omega }_{+} {\omega }_{-} 都分别相等. 表明该情形下, 外边界的作用主导, 内边界的作用可忽略; 随着孔的尺寸增加, {\omega }_{+} {\omega }_{-} 集体红移, 这由其对称的模式分布(即内外边界的激发电荷同号)决定, 与此同时它们的差值, 即SMS明显变小. 值得一提的是, 内外边界的电荷异号, 即典型的反对称模式, 归为径向振动, 不会产生SMS; 在孔很大时(如大于60 nm), 内外边界的反向环流相互抵消, {\omega }_{+} {\omega }_{-} 重合, SMS消失. 图6展示了孔的尺寸为60 nm时, 4类石墨烯环中偶极GMP的电势分布图, 印证了关于SMS消失的讨论. 此外, 鉴于该场分布, 亦可通过(6)式判断SMS效应, 如此时内外边界的电势已非常接近, 鉴于其相反的法线方向即对Q矩阵的贡献恰好相反, 导致最终的Q矩阵趋于0.

    图 5 4类石墨烯环中偶极GMP模式频率随孔尺寸的演化, 圆圈和菱形分别代表$ {\omega }_{+} $与$ {\omega }_{-} $模式\r\nFig. 5. Resonance frequencies of dipolar GMP as a function of the size of inner hole. The circles and rhombuses are the results of $ {\omega }_{+} $ and $ {\omega }_{-} $, respectively.
    图 5  4类石墨烯环中偶极GMP模式频率随孔尺寸的演化, 圆圈和菱形分别代表 {\omega }_{+} {\omega }_{-} 模式
    Fig. 5.  Resonance frequencies of dipolar GMP as a function of the size of inner hole. The circles and rhombuses are the results of {\omega }_{+} and {\omega }_{-} , respectively.
    图 6 4类石墨烯环中偶极GMP的电势分布图, 孔的尺寸为60 nm\r\nFig. 6. Potential distribution of the dipolar GMP in graphene rings of four kinds. The diameter or side length of each ring is 60 nm.
    图 6  4类石墨烯环中偶极GMP的电势分布图, 孔的尺寸为60 nm
    Fig. 6.  Potential distribution of the dipolar GMP in graphene rings of four kinds. The diameter or side length of each ring is 60 nm.

    石墨烯微纳结构中的GMP是决定体系磁光性质的重要电磁元激发, 对其开展研究大多仍沿用传统的3D建模方案, 即将2D石墨烯近似成厚度很小的3D块材, 代价是要消耗极大的计算资源. 本文基于库仑定律和电荷连续方程提出用于计算石墨烯体系中GMP的2DFEM, 通过本征值损失谱表征GMP的模式分布. 借助2DFEM, 探讨了4类石墨烯环中GMP 物性, 圆环表现出最强的磁光响应, 与其最高的对称性一致; 最低阶的偶极GMP都支持SMS, 当孔的尺寸很小时, 其对SMS的影响可以忽略, 随着孔的尺寸变大, 内外边界的相互作用增强并最终导致SMS消失. 该2DFEM具有一定的普适性, 可用于任意形状石墨烯结构中GMP的计算, 亦可拓展到其他2D材料, 为GMP的深入研究提供高效便捷的方案.

    [1]

    Barnes W L, Dereux A, Ebbesen T W 2003 Nature 424 824Google Scholar

    [2]

    Schuller J A, Barnard E S, Cai W, Jun Y C, White J S, Brongersma M L 2010 Nat. Mater. 9 193Google Scholar

    [3]

    Gramotnev D K, Bozhevolnyi S I 2010 Nat. Photon. 4 83Google Scholar

    [4]

    Knight M W, Sobhani H, Nordlander P, Halas N J 2011 Science 332 702Google Scholar

    [5]

    Goykhman I, Desiatov B, Khurgin J, Shappir J, Levy U 2011 Nano Lett. 11 2219Google Scholar

    [6]

    Senanayake P, Hung C, Shapiro J, Lin A, Liang B, Williams B S, Huffaker D L 2011 Nano Lett. 11 5279Google Scholar

    [7]

    Anker J N, Hall W P, Lyandres O, Shah N C, Zhao J, Van Duyne R P 2008 Nat. Mater. 7 442Google Scholar

    [8]

    Kabashin A V, Evans P, Pastkovsky S, Hendren W, Wurtz G A, Atkinson R, Pollard R, Podolskiy V A, Zayats A V 2009 Nat. Mater. 8 867Google Scholar

    [9]

    Zhang S, Bao K, Halas N J, Xu H, Nordlander P 2011 Nano Lett. 11 1657Google Scholar

    [10]

    Fang N, Lee H, Sun C, Zhang X 2005 Science 308 534Google Scholar

    [11]

    Liu Z, Steele J M, Srituravanich W, Pikus Y, Sun C, Zhang X 2005 Nano Lett. 5 1726Google Scholar

    [12]

    Kawata S, Inouye Y, Verma P 2009 Nat. Phton. 3 388Google Scholar

    [13]

    Christopher P, Xin H, Marimuthu A, Linic S 2012 Nat. Mater. 11 1044Google Scholar

    [14]

    Zhou L, Swearer D F, Zhang C, Robatjazi H, Zhao H, Henderson L, Dong L, Christopher P, Carter E A, Nordlander P, Halas N J 2018 Science 362 69Google Scholar

    [15]

    周利, 王取泉 2019 物理学报 68 147301Google Scholar

    Zhou L, Wang Q Q 2019 Acta Phys. Sin. 68 147301Google Scholar

    [16]

    Boltasseva A, Atwater H A 2011 Science 331 290Google Scholar

    [17]

    Khurgin J B 2015 Nat. Nanotech. 10 2Google Scholar

    [18]

    West P R, Ishii S, Naik G V, Emani N K, Shalaev V M, Boltasseva A 2010 Laser Photonics Rev. 4 795Google Scholar

    [19]

    Naik G V, Shalaev V M, Boltasseva A 2013 Adv. Mater. 25 3264Google Scholar

    [20]

    Knight M W, King N S, Liu L, Everitt H O, Nordlander P, Halas N J 2014 ACS Nano 8 834Google Scholar

    [21]

    Loa I, Syassen K, Monaco G, Vanko G, Krisch M, Hanfland M 2011 Phys. Rev. Lett. 107 086402Google Scholar

    [22]

    Fedyanin D Y, Yakubovsky D I, Kirtaev R V, Volkov V S 2016 Nano Lett. 16 326Google Scholar

    [23]

    Taliercio T, Biagioni P 2019 Nanophotonics 8 949Google Scholar

    [24]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsov A A 2004 Science 306 666Google Scholar

    [25]

    Frank I W, Tanenbaum D M, Van der Zande A M, McEuen P L 2007 J. Vac. Sci. Technol. B 25 2558Google Scholar

    [26]

    Rafiee M A, Rafiee J, Wang Z, Song H, Yu Z, Koratkar N 2009 ACS Nano 3 3884Google Scholar

    [27]

    Young R J, Kinloch I A, Gong L, Novoselov K S 2012 Comps. Sci. Technol. 72 1459Google Scholar

    [28]

    Balandin A A, Ghosh S, Bao W, Calizo I, Teweldebrhan D, Miao F, Lau C N 2008 Nano Lett. 8 902Google Scholar

    [29]

    Cai W, Moore A L, Zhu Y, Li X, Chen S, Shi L, Ruoff R S 2010 Nano Lett. 10 1645Google Scholar

    [30]

    Pop E, Varshney V, Roy A K 2012 MRS Bull. 37 1273Google Scholar

    [31]

    Nair R R, Blake P, Grigorenko A N, Novoselov K S, Booth T J, Stauber T, Peres N M R, Geim A K 2008 Science 320 1308Google Scholar

    [32]

    Wang F, Zhang Y, Tian C, Girit C, Zettl A, Crommie M, Shen Y R 2008 Science 320 206Google Scholar

    [33]

    Mueller T, Xia F, Avouris P 2010 Nat. Photon. 4 297Google Scholar

    [34]

    Liu M, Yin X, Ulin-Avila E, Geng B, Zentgraf T, Ju L, Wang F, Zhang X 2011 Nature 474 64Google Scholar

    [35]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Katsnelson M I, Grigorieva I V, Dubonos S V, Firsov A A 2005 Nature 438 197Google Scholar

    [36]

    Geim A K, Novoselov K S 2007 Nat. Mater. 6 183Google Scholar

    [37]

    Castro Neto A H, Guinea F, Peres N M R, Novoselov K S, Geim A K 2009 Rev. Mod. Phys. 81 109Google Scholar

    [38]

    Das Sarma S, Adam S, Hwang E H, Rossi E 2011 Rev. Mod. Phys. 83 407Google Scholar

    [39]

    Mas-Balleste R, Gomez-Navarro C, Gomez-Herrero J, Zamora F 2011 Nanoscale 3 20Google Scholar

    [40]

    Fiori G, Bonaccorso F, Iannaccone G, Palacios T, Neumaier D, Seabaugh A, Bnerjee S K, Colombo L 2014 Nat. Nanotechonl. 9 768Google Scholar

    [41]

    Novoselov K S, Mishchenko A, Carvalho A, Castro Neto A H 2016 Science 353 aac9439Google Scholar

    [42]

    Wunsch B, Stauber T, Sols F, Guinea F 2006 New J. Phys. 8 318Google Scholar

    [43]

    Hwang E H, Das Sarma S 2007 Phys. Rev. B 75 205418Google Scholar

    [44]

    Polini M, Asgari R, Borghi G, Barlas Y, Pereg-Barnea T, MacDonald A H 2008 Phys. Rev. B 77 081411Google Scholar

    [45]

    Jablan M, Buljan H, Soljacic M 2009 Phys. Rev. B 80 245435Google Scholar

    [46]

    Koppens F H L, Chang D E, de Abajo F J G 2011 Nano Lett. 11 3370Google Scholar

    [47]

    Brar V W, Jang M S, Sherrott M, Lopez J J, Atwater H A 2013 Nano Lett. 13 2541Google Scholar

    [48]

    de Abajo F J G 2014 ACS Photon. 1 135Google Scholar

    [49]

    Chen J, Badioli M, Alonso-Gonzalez P, Thongrattanasiri S, Huth F, Osmond J, Spasenovic M, Centeno A, Pesquera A, Godignon P, Elorza A Z, Camara N, de Abajo F J G, Hillenbrand R, Koppens F H L 2012 Nature 487 77Google Scholar

    [50]

    Fei Z, Rodin A S, Andreev G O, Bao W, Mcleod A S, Wagner M, Zhang L, Zhao Z, Thiemens M, Dominguez G, Fogler M M, Castro Neto A H, Lau C N, Keilmann F, Basov D N 2012 Nature 487 82Google Scholar

    [51]

    Morozov S V, Novoselov K S, Katsnelson M I, Schedin F, Elias D C, Jaszczak J A, Geim A K 2008 Phys. Rev. Lett. 100 016602Google Scholar

    [52]

    Du X, Skachko I, Barker A, Andrei E Y 2008 Nat. Nanotechnol. 3 491Google Scholar

    [53]

    Bolotin K I, Sikes K J, Hone J, Stormer H L, Kim P 2008 Phys. Rev. Lett. 101 096802Google Scholar

    [54]

    Ju L, Geng B, Horng J, Girit C, Martin M, Hao Z, Bechtel H A, Liang X, Zettl A, Shen Y R, Wang F 2011 Nature Nanotechnol. 6 630Google Scholar

    [55]

    Yan H, Li X, Chandra B, Tulevski G, Wu Y, Freitag M, Zhu W, Avouris P, Xia F 2012 Nature Nanotechnol. 7 330Google Scholar

    [56]

    Low T, Avouris P 2014 ACS Nano 8 1086Google Scholar

    [57]

    吴晨晨, 郭相东, 胡海, 杨晓霞, 戴庆 2019 物理学报 68 148103Google Scholar

    Wu C C, Guo X D, Hu H, Yang X X, Dai Q 2019 Acta Phys. Sin. 68 148103Google Scholar

    [58]

    Eda G, Maier S A 2013 ACS Nano 7 5660Google Scholar

    [59]

    Xia F, Wang H, Xiao D, Dubey M, Ramasubramaniam A 2014 Nat. Photon. 8 899Google Scholar

    [60]

    Li Y, Li Z, Cheng C, Shan H, Zheng L, Fang Z 2017 Adv. Sci. 4 1600430Google Scholar

    [61]

    Zhang Y, Tan Y, Stormer H L, Kim P 2005 Nature 438 201Google Scholar

    [62]

    Novoselov K S, Jiang Z, Zhang Y, Morozov S V, Stormer H L, Zeitler U, Maan J C, Boebinger G S, Kim P, Geim A K 2007 Science 315 1379Google Scholar

    [63]

    Guinea F, Katsnelson M I, Geim A K 2010 Nat. Phys. 6 30Google Scholar

    [64]

    Du X, Skachko I, Duerr F, Luican A, Andrei E Y 2009 Nature 462 192Google Scholar

    [65]

    Bolotin K L, Ghahari F, Shulman M D, Stormer H L, Kim P 2009 Nature 462 196Google Scholar

    [66]

    Dean C R, Young A F, Cadden-Zimansky P, Wang L, Ren H, Watanabe K, Taniguchi T, Kim P, Hone J, Shepard K L 2011 Nat. Phys. 7 693Google Scholar

    [67]

    Kane C L, Mele E J 2005 Phys. Rev. Lett. 95 226801Google Scholar

    [68]

    Ferreira A, Rappoport T G, Cazalilla M A, Castro Neto A H 2014 Phys. Rev. Lett. 112 066601Google Scholar

    [69]

    Crassee I, Levallois J, Walter A L, Ostler M, Bostwick A, Rotenberg E, Seyller T, ven der Marel D, Kuzmenko A B 2011 Nat. Phys. 7 48Google Scholar

    [70]

    Fallahi A, Perruisseau-Carrier J 2012 Appl. Phys. Lett. 101 231605Google Scholar

    [71]

    Sounas D L, Skulason H S, Nguyen H V, Guermoune A, Slaj M, Szkopek T, Caloz C 2013 Appl. Phys. Lett. 102 191901Google Scholar

    [72]

    Shimano R, Yumoto G, Yoo J Y, Matsunaga R, Tanabe S, Hibino H, Morimoto T, Aoki H 2013 Nat. Commun. 4 1841Google Scholar

    [73]

    Fetter A L 1985 Phys. Rev. B 32 7676Google Scholar

    [74]

    Mast D B, Dahm A J, Fetter A L 1985 Phys. Rev. Lett. 54 1706Google Scholar

    [75]

    Wu J, Hawrylak P, Eliasson G, Quinn J J 1986 Phys. Rev. B 33 7091Google Scholar

    [76]

    Armelles G, Cebollada A, Garcia-Martin A, Gonzalez M U 2013 Adv. Opt. Mater. 1 10Google Scholar

    [77]

    Poumirol J M, Liu P Q, Slipchenko T M, Nikitin A Y, Martin-Moreno L, Faist J, Kuzmenko A B 2017 Nat. Commun. 8 14626Google Scholar

    [78]

    Bychkov Y A, Martinez G 2008 Phys. Rev. B 77 125417Google Scholar

    [79]

    Berman O L, Gumbs G, Lozovik Y E 2008 Phys. Rev. B 78 085401Google Scholar

    [80]

    Wu J, Chen S, Roslyak O, Gumbs G, Lin M 2011 ACS Nano 5 1026Google Scholar

    [81]

    Crassee I, Orlita M, Potemski M, Walter A L, Ostler M, Seyller T H, Gaponenko I, Chen J, Kuzmenko A B 2012 Nano Lett. 12 2470Google Scholar

    [82]

    Ferreira A, Peres N M R, Castro Neto A H 2012 Phys. Rev. B 85 205426Google Scholar

    [83]

    Mishchenko E G, Shyton A V, Silvestrov P G 2010 Phys. Rev. Lett. 104 156806Google Scholar

    [84]

    Petkovic I, Williams F I B, Bennaceur K, Portier F, Roche P, Glattli D C 2013 Phys. Rev. Lett. 110 016801Google Scholar

    [85]

    Sokolik A A, Lozovik 2019 Phys. Rev. B 100 125409Google Scholar

    [86]

    Yan H, Li Z, Li X, Zhu W, Avouris P, Xia F 2012 Nano Lett. 12 3766Google Scholar

    [87]

    Kumada N, Roulleau P, Roche B, Hashisaka M, Hibino H, Petkovic I, Glattli D C 2014 Phys. Rev. Lett. 113 266601Google Scholar

    [88]

    Lin X, Xu Y, Zhang B, Hao R, Chen H, Li E 2013 New J. Phys. 15 113003Google Scholar

    [89]

    Chamanara N, Sounas D, Caloz C 2013 Opt. Express 21 11248Google Scholar

    [90]

    Jin D, Lu L, Wang Z, Fang C, Joannopoulos J D, Soljacic M, Fu L, Fang N 2016 Nat. Commun. 7 13486Google Scholar

    [91]

    Pan D, Yu R, Xu H, de Abajo F J G 2017 Nat. Commun. 8 1243Google Scholar

    [92]

    Jin D, Christensen T, Soljacic M, Fang N, Lu L, Zhang X 2017 Phys. Rev. Lett. 118 245301Google Scholar

    [93]

    Wang W, Kinaret J M, Apell S P 2012 Phys. Rev. B 85 235444Google Scholar

    [94]

    Wang W, Apell S P, Kinaret J M 2012 Phys. Rev. B 86 125450Google Scholar

    [95]

    Jiao N, Kang S, Han K, Shen X, Wang W 2019 Phys. Rev. B 99 195447Google Scholar

    [96]

    Jiao N, Kang S, Han K, Shen X, Wang W 2021 Phys. Rev. B 103 085405Google Scholar

    [97]

    Wang N, Ding L, Wang W 2022 Phys. Rev. B 105 235435Google Scholar

    [98]

    Gusynin V P, Sharapov S G, Carbotte J P 2007 J. Phys. Condens. Matter 19 026222Google Scholar

    [99]

    Gusynin V P, Sharapov S G, Carbotte J P 2007 Phys. Rev. B 75 165407Google Scholar

    [100]

    Roldan R, Fuchs J, Goerbig M O 2009 Phys. Rev. B 80 085408Google Scholar

    [101]

    Morimoto T, Hatsugai Y, Aoki H 2009 Phys. Rev. Lett. 103 116803Google Scholar

    [102]

    Yang C, Peeters F M, Xu W 2010 Phys. Rev. B 82 205428Google Scholar

    [103]

    Wang W, Apell S P, Kinaret J M 2011 Phys. Rev. B 84 085423Google Scholar

    [104]

    Wang W, Christensen T, Jauho A P, Thygesen K S, Wubs M, Mortensen N A 2015 Sci. Rep. 5 9535Google Scholar

    [105]

    Wang W, Xiao S, Mortensen N A 2016 Phys. Rev. B 93 165407Google Scholar

    [106]

    Geuzaine C, Remacle J F 2009 Int. J. Numer. Meth. Eng. 79 1309Google Scholar

  • 图 1  石墨烯圆环(a)、方环(b)、方孔圆环(c)、圆孔方环(d)网格划分示意图, 环的直径或边长均为100 nm, 孔的直径或边长均为50 nm

    Figure 1.  Schematic diagrams of mesh generation for graphene circular ring (a), square ring (b), circular ring with square hole (c), and square ring with circular hole (d). The diameteror side length of ring is 100 nm, and that of hole is 50 nm.

    图 2  石墨烯圆盘的偶极GMP谱线, 数据分别来自2DFEM (a)和3DFEM (b), 其中虚线标记了蓝色曲线半峰的高度

    Figure 2.  Dipolar GMP spectra in graphene disk, obtained from 2DFEM (a) and 3DFEM (b), respectively. The dashed lines mark the position of half maximum of blue curves.

    图 3  石墨烯圆盘(a)和方块(b)的本征值损失谱, 圆盘和方块的直径与边长均为100 nm

    Figure 3.  Eigenvalue loss spectra of graphene disk (a) and graphene square (b). The diameter of disk and the side length of square are both 100 nm.

    图 4  4类石墨烯环的本征值损失谱, 椭圆标记了发生SMS的GMP模式

    Figure 4.  Eigenvalue loss spectra of graphene rings of 4 kinds. The ellipses mark the symmetric mode splitting of those GMPs.

    图 5  4类石墨烯环中偶极GMP模式频率随孔尺寸的演化, 圆圈和菱形分别代表 {\omega }_{+} {\omega }_{-} 模式

    Figure 5.  Resonance frequencies of dipolar GMP as a function of the size of inner hole. The circles and rhombuses are the results of {\omega }_{+} and {\omega }_{-} , respectively.

    图 6  4类石墨烯环中偶极GMP的电势分布图, 孔的尺寸为60 nm

    Figure 6.  Potential distribution of the dipolar GMP in graphene rings of four kinds. The diameter or side length of each ring is 60 nm.

  • [1]

    Barnes W L, Dereux A, Ebbesen T W 2003 Nature 424 824Google Scholar

    [2]

    Schuller J A, Barnard E S, Cai W, Jun Y C, White J S, Brongersma M L 2010 Nat. Mater. 9 193Google Scholar

    [3]

    Gramotnev D K, Bozhevolnyi S I 2010 Nat. Photon. 4 83Google Scholar

    [4]

    Knight M W, Sobhani H, Nordlander P, Halas N J 2011 Science 332 702Google Scholar

    [5]

    Goykhman I, Desiatov B, Khurgin J, Shappir J, Levy U 2011 Nano Lett. 11 2219Google Scholar

    [6]

    Senanayake P, Hung C, Shapiro J, Lin A, Liang B, Williams B S, Huffaker D L 2011 Nano Lett. 11 5279Google Scholar

    [7]

    Anker J N, Hall W P, Lyandres O, Shah N C, Zhao J, Van Duyne R P 2008 Nat. Mater. 7 442Google Scholar

    [8]

    Kabashin A V, Evans P, Pastkovsky S, Hendren W, Wurtz G A, Atkinson R, Pollard R, Podolskiy V A, Zayats A V 2009 Nat. Mater. 8 867Google Scholar

    [9]

    Zhang S, Bao K, Halas N J, Xu H, Nordlander P 2011 Nano Lett. 11 1657Google Scholar

    [10]

    Fang N, Lee H, Sun C, Zhang X 2005 Science 308 534Google Scholar

    [11]

    Liu Z, Steele J M, Srituravanich W, Pikus Y, Sun C, Zhang X 2005 Nano Lett. 5 1726Google Scholar

    [12]

    Kawata S, Inouye Y, Verma P 2009 Nat. Phton. 3 388Google Scholar

    [13]

    Christopher P, Xin H, Marimuthu A, Linic S 2012 Nat. Mater. 11 1044Google Scholar

    [14]

    Zhou L, Swearer D F, Zhang C, Robatjazi H, Zhao H, Henderson L, Dong L, Christopher P, Carter E A, Nordlander P, Halas N J 2018 Science 362 69Google Scholar

    [15]

    周利, 王取泉 2019 物理学报 68 147301Google Scholar

    Zhou L, Wang Q Q 2019 Acta Phys. Sin. 68 147301Google Scholar

    [16]

    Boltasseva A, Atwater H A 2011 Science 331 290Google Scholar

    [17]

    Khurgin J B 2015 Nat. Nanotech. 10 2Google Scholar

    [18]

    West P R, Ishii S, Naik G V, Emani N K, Shalaev V M, Boltasseva A 2010 Laser Photonics Rev. 4 795Google Scholar

    [19]

    Naik G V, Shalaev V M, Boltasseva A 2013 Adv. Mater. 25 3264Google Scholar

    [20]

    Knight M W, King N S, Liu L, Everitt H O, Nordlander P, Halas N J 2014 ACS Nano 8 834Google Scholar

    [21]

    Loa I, Syassen K, Monaco G, Vanko G, Krisch M, Hanfland M 2011 Phys. Rev. Lett. 107 086402Google Scholar

    [22]

    Fedyanin D Y, Yakubovsky D I, Kirtaev R V, Volkov V S 2016 Nano Lett. 16 326Google Scholar

    [23]

    Taliercio T, Biagioni P 2019 Nanophotonics 8 949Google Scholar

    [24]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsov A A 2004 Science 306 666Google Scholar

    [25]

    Frank I W, Tanenbaum D M, Van der Zande A M, McEuen P L 2007 J. Vac. Sci. Technol. B 25 2558Google Scholar

    [26]

    Rafiee M A, Rafiee J, Wang Z, Song H, Yu Z, Koratkar N 2009 ACS Nano 3 3884Google Scholar

    [27]

    Young R J, Kinloch I A, Gong L, Novoselov K S 2012 Comps. Sci. Technol. 72 1459Google Scholar

    [28]

    Balandin A A, Ghosh S, Bao W, Calizo I, Teweldebrhan D, Miao F, Lau C N 2008 Nano Lett. 8 902Google Scholar

    [29]

    Cai W, Moore A L, Zhu Y, Li X, Chen S, Shi L, Ruoff R S 2010 Nano Lett. 10 1645Google Scholar

    [30]

    Pop E, Varshney V, Roy A K 2012 MRS Bull. 37 1273Google Scholar

    [31]

    Nair R R, Blake P, Grigorenko A N, Novoselov K S, Booth T J, Stauber T, Peres N M R, Geim A K 2008 Science 320 1308Google Scholar

    [32]

    Wang F, Zhang Y, Tian C, Girit C, Zettl A, Crommie M, Shen Y R 2008 Science 320 206Google Scholar

    [33]

    Mueller T, Xia F, Avouris P 2010 Nat. Photon. 4 297Google Scholar

    [34]

    Liu M, Yin X, Ulin-Avila E, Geng B, Zentgraf T, Ju L, Wang F, Zhang X 2011 Nature 474 64Google Scholar

    [35]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Katsnelson M I, Grigorieva I V, Dubonos S V, Firsov A A 2005 Nature 438 197Google Scholar

    [36]

    Geim A K, Novoselov K S 2007 Nat. Mater. 6 183Google Scholar

    [37]

    Castro Neto A H, Guinea F, Peres N M R, Novoselov K S, Geim A K 2009 Rev. Mod. Phys. 81 109Google Scholar

    [38]

    Das Sarma S, Adam S, Hwang E H, Rossi E 2011 Rev. Mod. Phys. 83 407Google Scholar

    [39]

    Mas-Balleste R, Gomez-Navarro C, Gomez-Herrero J, Zamora F 2011 Nanoscale 3 20Google Scholar

    [40]

    Fiori G, Bonaccorso F, Iannaccone G, Palacios T, Neumaier D, Seabaugh A, Bnerjee S K, Colombo L 2014 Nat. Nanotechonl. 9 768Google Scholar

    [41]

    Novoselov K S, Mishchenko A, Carvalho A, Castro Neto A H 2016 Science 353 aac9439Google Scholar

    [42]

    Wunsch B, Stauber T, Sols F, Guinea F 2006 New J. Phys. 8 318Google Scholar

    [43]

    Hwang E H, Das Sarma S 2007 Phys. Rev. B 75 205418Google Scholar

    [44]

    Polini M, Asgari R, Borghi G, Barlas Y, Pereg-Barnea T, MacDonald A H 2008 Phys. Rev. B 77 081411Google Scholar

    [45]

    Jablan M, Buljan H, Soljacic M 2009 Phys. Rev. B 80 245435Google Scholar

    [46]

    Koppens F H L, Chang D E, de Abajo F J G 2011 Nano Lett. 11 3370Google Scholar

    [47]

    Brar V W, Jang M S, Sherrott M, Lopez J J, Atwater H A 2013 Nano Lett. 13 2541Google Scholar

    [48]

    de Abajo F J G 2014 ACS Photon. 1 135Google Scholar

    [49]

    Chen J, Badioli M, Alonso-Gonzalez P, Thongrattanasiri S, Huth F, Osmond J, Spasenovic M, Centeno A, Pesquera A, Godignon P, Elorza A Z, Camara N, de Abajo F J G, Hillenbrand R, Koppens F H L 2012 Nature 487 77Google Scholar

    [50]

    Fei Z, Rodin A S, Andreev G O, Bao W, Mcleod A S, Wagner M, Zhang L, Zhao Z, Thiemens M, Dominguez G, Fogler M M, Castro Neto A H, Lau C N, Keilmann F, Basov D N 2012 Nature 487 82Google Scholar

    [51]

    Morozov S V, Novoselov K S, Katsnelson M I, Schedin F, Elias D C, Jaszczak J A, Geim A K 2008 Phys. Rev. Lett. 100 016602Google Scholar

    [52]

    Du X, Skachko I, Barker A, Andrei E Y 2008 Nat. Nanotechnol. 3 491Google Scholar

    [53]

    Bolotin K I, Sikes K J, Hone J, Stormer H L, Kim P 2008 Phys. Rev. Lett. 101 096802Google Scholar

    [54]

    Ju L, Geng B, Horng J, Girit C, Martin M, Hao Z, Bechtel H A, Liang X, Zettl A, Shen Y R, Wang F 2011 Nature Nanotechnol. 6 630Google Scholar

    [55]

    Yan H, Li X, Chandra B, Tulevski G, Wu Y, Freitag M, Zhu W, Avouris P, Xia F 2012 Nature Nanotechnol. 7 330Google Scholar

    [56]

    Low T, Avouris P 2014 ACS Nano 8 1086Google Scholar

    [57]

    吴晨晨, 郭相东, 胡海, 杨晓霞, 戴庆 2019 物理学报 68 148103Google Scholar

    Wu C C, Guo X D, Hu H, Yang X X, Dai Q 2019 Acta Phys. Sin. 68 148103Google Scholar

    [58]

    Eda G, Maier S A 2013 ACS Nano 7 5660Google Scholar

    [59]

    Xia F, Wang H, Xiao D, Dubey M, Ramasubramaniam A 2014 Nat. Photon. 8 899Google Scholar

    [60]

    Li Y, Li Z, Cheng C, Shan H, Zheng L, Fang Z 2017 Adv. Sci. 4 1600430Google Scholar

    [61]

    Zhang Y, Tan Y, Stormer H L, Kim P 2005 Nature 438 201Google Scholar

    [62]

    Novoselov K S, Jiang Z, Zhang Y, Morozov S V, Stormer H L, Zeitler U, Maan J C, Boebinger G S, Kim P, Geim A K 2007 Science 315 1379Google Scholar

    [63]

    Guinea F, Katsnelson M I, Geim A K 2010 Nat. Phys. 6 30Google Scholar

    [64]

    Du X, Skachko I, Duerr F, Luican A, Andrei E Y 2009 Nature 462 192Google Scholar

    [65]

    Bolotin K L, Ghahari F, Shulman M D, Stormer H L, Kim P 2009 Nature 462 196Google Scholar

    [66]

    Dean C R, Young A F, Cadden-Zimansky P, Wang L, Ren H, Watanabe K, Taniguchi T, Kim P, Hone J, Shepard K L 2011 Nat. Phys. 7 693Google Scholar

    [67]

    Kane C L, Mele E J 2005 Phys. Rev. Lett. 95 226801Google Scholar

    [68]

    Ferreira A, Rappoport T G, Cazalilla M A, Castro Neto A H 2014 Phys. Rev. Lett. 112 066601Google Scholar

    [69]

    Crassee I, Levallois J, Walter A L, Ostler M, Bostwick A, Rotenberg E, Seyller T, ven der Marel D, Kuzmenko A B 2011 Nat. Phys. 7 48Google Scholar

    [70]

    Fallahi A, Perruisseau-Carrier J 2012 Appl. Phys. Lett. 101 231605Google Scholar

    [71]

    Sounas D L, Skulason H S, Nguyen H V, Guermoune A, Slaj M, Szkopek T, Caloz C 2013 Appl. Phys. Lett. 102 191901Google Scholar

    [72]

    Shimano R, Yumoto G, Yoo J Y, Matsunaga R, Tanabe S, Hibino H, Morimoto T, Aoki H 2013 Nat. Commun. 4 1841Google Scholar

    [73]

    Fetter A L 1985 Phys. Rev. B 32 7676Google Scholar

    [74]

    Mast D B, Dahm A J, Fetter A L 1985 Phys. Rev. Lett. 54 1706Google Scholar

    [75]

    Wu J, Hawrylak P, Eliasson G, Quinn J J 1986 Phys. Rev. B 33 7091Google Scholar

    [76]

    Armelles G, Cebollada A, Garcia-Martin A, Gonzalez M U 2013 Adv. Opt. Mater. 1 10Google Scholar

    [77]

    Poumirol J M, Liu P Q, Slipchenko T M, Nikitin A Y, Martin-Moreno L, Faist J, Kuzmenko A B 2017 Nat. Commun. 8 14626Google Scholar

    [78]

    Bychkov Y A, Martinez G 2008 Phys. Rev. B 77 125417Google Scholar

    [79]

    Berman O L, Gumbs G, Lozovik Y E 2008 Phys. Rev. B 78 085401Google Scholar

    [80]

    Wu J, Chen S, Roslyak O, Gumbs G, Lin M 2011 ACS Nano 5 1026Google Scholar

    [81]

    Crassee I, Orlita M, Potemski M, Walter A L, Ostler M, Seyller T H, Gaponenko I, Chen J, Kuzmenko A B 2012 Nano Lett. 12 2470Google Scholar

    [82]

    Ferreira A, Peres N M R, Castro Neto A H 2012 Phys. Rev. B 85 205426Google Scholar

    [83]

    Mishchenko E G, Shyton A V, Silvestrov P G 2010 Phys. Rev. Lett. 104 156806Google Scholar

    [84]

    Petkovic I, Williams F I B, Bennaceur K, Portier F, Roche P, Glattli D C 2013 Phys. Rev. Lett. 110 016801Google Scholar

    [85]

    Sokolik A A, Lozovik 2019 Phys. Rev. B 100 125409Google Scholar

    [86]

    Yan H, Li Z, Li X, Zhu W, Avouris P, Xia F 2012 Nano Lett. 12 3766Google Scholar

    [87]

    Kumada N, Roulleau P, Roche B, Hashisaka M, Hibino H, Petkovic I, Glattli D C 2014 Phys. Rev. Lett. 113 266601Google Scholar

    [88]

    Lin X, Xu Y, Zhang B, Hao R, Chen H, Li E 2013 New J. Phys. 15 113003Google Scholar

    [89]

    Chamanara N, Sounas D, Caloz C 2013 Opt. Express 21 11248Google Scholar

    [90]

    Jin D, Lu L, Wang Z, Fang C, Joannopoulos J D, Soljacic M, Fu L, Fang N 2016 Nat. Commun. 7 13486Google Scholar

    [91]

    Pan D, Yu R, Xu H, de Abajo F J G 2017 Nat. Commun. 8 1243Google Scholar

    [92]

    Jin D, Christensen T, Soljacic M, Fang N, Lu L, Zhang X 2017 Phys. Rev. Lett. 118 245301Google Scholar

    [93]

    Wang W, Kinaret J M, Apell S P 2012 Phys. Rev. B 85 235444Google Scholar

    [94]

    Wang W, Apell S P, Kinaret J M 2012 Phys. Rev. B 86 125450Google Scholar

    [95]

    Jiao N, Kang S, Han K, Shen X, Wang W 2019 Phys. Rev. B 99 195447Google Scholar

    [96]

    Jiao N, Kang S, Han K, Shen X, Wang W 2021 Phys. Rev. B 103 085405Google Scholar

    [97]

    Wang N, Ding L, Wang W 2022 Phys. Rev. B 105 235435Google Scholar

    [98]

    Gusynin V P, Sharapov S G, Carbotte J P 2007 J. Phys. Condens. Matter 19 026222Google Scholar

    [99]

    Gusynin V P, Sharapov S G, Carbotte J P 2007 Phys. Rev. B 75 165407Google Scholar

    [100]

    Roldan R, Fuchs J, Goerbig M O 2009 Phys. Rev. B 80 085408Google Scholar

    [101]

    Morimoto T, Hatsugai Y, Aoki H 2009 Phys. Rev. Lett. 103 116803Google Scholar

    [102]

    Yang C, Peeters F M, Xu W 2010 Phys. Rev. B 82 205428Google Scholar

    [103]

    Wang W, Apell S P, Kinaret J M 2011 Phys. Rev. B 84 085423Google Scholar

    [104]

    Wang W, Christensen T, Jauho A P, Thygesen K S, Wubs M, Mortensen N A 2015 Sci. Rep. 5 9535Google Scholar

    [105]

    Wang W, Xiao S, Mortensen N A 2016 Phys. Rev. B 93 165407Google Scholar

    [106]

    Geuzaine C, Remacle J F 2009 Int. J. Numer. Meth. Eng. 79 1309Google Scholar

  • [1] Duan Yu, Dai Xiao-Kang, Wu Chen-Chen, Yang Xiao-Xia. Tunable acoustic graphene plasmon enhanced nano-infrared spectroscopy. Acta Physica Sinica, 2024, 73(13): 138101. doi: 10.7498/aps.73.20240489
    [2] Zhou Kun, Ma Hao-Yue, Sun Xi-Xian, Wu Xiao-Hu. Active tuning hBN phonon polaritons and spontaneous emission rates based on VO2 and graphene. Acta Physica Sinica, 2023, 72(7): 074201. doi: 10.7498/aps.72.20222167
    [3] Wang Jian, Zhang Chao-Yue, Yao Zhao-Yu, Zhang Chi, Xu Feng, Yang Yuan. A method of rapidly designing graphene-based terahertz diffusion surface. Acta Physica Sinica, 2021, 70(3): 034102. doi: 10.7498/aps.70.20201034
    [4] Li Ze-Yu, Jiang Qu-Han, Ma Teng-Zhou, Yuan Ying-Hao, Chen Lin. Multi-parameter tunable phase transition based terahertz graphene plasmons and its application. Acta Physica Sinica, 2021, 70(22): 224202. doi: 10.7498/aps.70.20210445
    [5] Liu Liang, Han De-Zhuan, Shi Lei. Plasmonic band structures and its applications. Acta Physica Sinica, 2020, 69(15): 157301. doi: 10.7498/aps.69.20200193
    [6] Hu Bao-Jing, Huang Ming, Li Peng, Yang Jing-Jing. Multiband plasmon-induced transparency based on nanometals-graphene hybrid model. Acta Physica Sinica, 2020, 69(17): 174201. doi: 10.7498/aps.69.20200200
    [7] Zhao Cheng-Xiang, Qie Yuan, Yu Yao, Ma Rong-Rong, Qin Jun-Fei, Liu Yan. Enhanced optical absorption of graphene by plasmon. Acta Physica Sinica, 2020, 69(6): 067801. doi: 10.7498/aps.69.20191645
    [8] Cui Shu-Wen, Li Lu, Wei Lian-Jia, Qian Ping. Theoretical study of density functional of confined CO oxidation reaction between bilayer graphene. Acta Physica Sinica, 2019, 68(21): 218101. doi: 10.7498/aps.68.20190447
    [9] Wu Chen-Chen, Guo Xiang-Dong, Hu Hai, Yang Xiao-Xia, Dai Qing. Graphene plasmon enhanced infrared spectroscopy. Acta Physica Sinica, 2019, 68(14): 148103. doi: 10.7498/aps.68.20190903
    [10] Chen Hao, Zhang Xiao-Xia, Wang Hong, Ji Yue-Hua. Near-infrared absorption of graphene-metal nanostructure based on magnetic polaritons. Acta Physica Sinica, 2018, 67(11): 118101. doi: 10.7498/aps.67.20180196
    [11] Deng Hong-Mei, Huang Lei, Li Jing, Lu Ye, Li Chuan-Qi. Tunable unidirectional surface plasmon polariton coupler utilizing graphene-based asymmetric nanoantenna pairs. Acta Physica Sinica, 2017, 66(14): 145201. doi: 10.7498/aps.66.145201
    [12] Lu Ya-Xin, Ma Ning. The coupled electromagnetic field effects on quantum magnetic oscillations of graphene. Acta Physica Sinica, 2016, 65(2): 027502. doi: 10.7498/aps.65.027502
    [13] Zhang Chao-Jie, Zhou Ting, Du Xin-Peng, Wang Tong-Biao, Liu Nian-Hua. Enhancement of quantum friction via coupling of surface phonon polariton and graphene plasmons. Acta Physica Sinica, 2016, 65(23): 236801. doi: 10.7498/aps.65.236801
    [14] Yang Jing-Jing, Li Jun-Jie, Deng Wei, Cheng Cheng, Huang Ming. Transmission mode of a single layer graphene and its performance in the detection of the vibration spectrum of gas molecular. Acta Physica Sinica, 2015, 64(19): 198102. doi: 10.7498/aps.64.198102
    [15] Sheng Shi-Wei, Li Kang, Kong Fan-Min, Yue Qing-Yang, Zhuang Hua-Wei, Zhao Jia. Tooth-shaped plasmonic filter based on graphene nanoribbon. Acta Physica Sinica, 2015, 64(10): 108402. doi: 10.7498/aps.64.108402
    [16] Zhang Bao-Lei, Wang Jia-Xu, Xiao Ke, Li Jun-Yang. Quasi-static finite element calculation of interaction between graphene and nanoprobe. Acta Physica Sinica, 2014, 63(15): 154601. doi: 10.7498/aps.63.154601
    [17] Deng Wei-Yin, Zhu Rui, Deng Wen-Ji. Electronic state of the limited graphene. Acta Physica Sinica, 2013, 62(8): 087301. doi: 10.7498/aps.62.087301
    [18] Zou Wei-Bo, Zhou Jun, Jin Li, Zhang Hao-Peng. Properties of localized surface plasmon resonance of gold nanoshell pairs. Acta Physica Sinica, 2012, 61(9): 097805. doi: 10.7498/aps.61.097805
    [19] Sun Hong-Xiang, Xu Bai-Qiang, Wang Ji-Jun, Xu Gui-Dong, Xu Chen-Guang, Wang Feng. Numerical simulation of laser-generated Rayleigh wave by finite element method on viscoelastic materials. Acta Physica Sinica, 2009, 58(9): 6344-6350. doi: 10.7498/aps.58.6344
    [20] Feng Yong-Ping, Cui Jun-Zhi, Deng Ming-Xiang. The two-scale finite element computation for thermoelastic problem in periodic perforated domain. Acta Physica Sinica, 2009, 58(13): 327-S337. doi: 10.7498/aps.58.327
Metrics
  • Abstract views:  5001
  • PDF Downloads:  127
Publishing process
  • Received Date:  31 December 2022
  • Accepted Date:  06 February 2023
  • Available Online:  17 February 2023
  • Published Online:  20 April 2023

/

返回文章
返回