石墨烯的介电常数一般是多少？

推荐一个小工具，可以很方便地计算石墨烯的电导率、介电常数、折射率。

下载链接：https:// opt-simul.com/archives/ 705

以下是正文：

石墨烯 的电导率、介电常数、折射率等电磁参数主要与频率 \omega 、化学势 \mu_{\rm c} 、散射率 \Gamma 、温度 T 、层数 N_{\rm g} 、单层石墨烯厚度 t_{\rm g} 等因素有关。其电导率可以用 Kubo公式 来计算，然后可以根据电导率进一步计算 介电常数 和 折射率 。

1、Kubo公式

描述石墨烯电导率的Kubo公式由 带内电导率 \sigma_{\rm intra} 和 带间电导率 \sigma_{\rm inter} 两部分组成，总电导率是两者的和。当石墨烯层数 N_{\rm g} 不大时，总电导率可以写为 ^{[1] [2] [3] [4] [5] [6]} ：

\sigma(\omega,\mu_{\rm c},\Gamma,T)=N_{\rm g}(\sigma_{\rm intra}+\sigma_{\rm inter}) (1a)

其中

\sigma_{\rm intra}=\frac{-{\rm i}e^2}{\pi\hbar^2(\omega+2{\rm i} \Gamma)}\int_{0}^{\infty}\xi\left[ \frac{\partial f_d(\xi)}{\partial \xi}-\frac{\partial f_d(-\xi)}{\partial \xi} \right]{\rm d}\xi (1b)

\sigma_{\rm inter}=\frac{{\rm i}e^2(\omega+2{\rm i}\Gamma)}{\pi\hbar^2}\int_{0}^{\infty}\frac{f_d(-\xi)-f_d(\xi)}{(\omega+2{\rm i}\Gamma)^2-4(\xi/\hbar)^2}{\rm d}\xi (1c)

f_{\rm d}(\xi)=\left[ {\rm exp}\left( \frac{\xi-\mu_{\rm c}}{k_{\rm B}T} \right)+1 \right]^{-1} 是 费米狄拉克分布 (1d)

在计算需要注意单位统一：如果角频率 \omega 的单位取[rad/s]，那么 \Gamma 的单位也必须是[rad/s]；如果光子能量 \hbar\omega 的单位取[eV]，那么 \hbar\Gamma 的单位也必须是[eV]；在 费米狄拉克分布 中，如果 \mu_{\rm c} 的单位取[eV]，那么 \xi 和 k_{\rm B}T 的单位也必须是[eV]。为了计算方便，建议所有能量单位都取为[eV]。另外还需注意，石墨烯是二维材料，上面Eq.(1)算出来的电导率是 表面电导率 ，单位是[S]，而不是[S/m]。

在得到电导率后，利用公式 ^[7]

\epsilon_{//}=1+\frac{{\rm i}\sigma}{\omega\epsilon_0 N_{\rm g} t_{\rm g}} 和 n_{//}=\sqrt{\epsilon_{//}} (2a)

分别得到石墨烯的 切向相对介电常数 和 切向折射率 ，而石墨烯的 法向相对介电常数 和 法向折射率 都是1，即：

\epsilon_{\bot}=1 和 n_{\bot}=1 (2b)

Kubo公式 对紫外到太赫兹波段都适用，但是我们可以看到 Eq.(1) 还是比较复杂的，要正确计算这些式子并不容易。因此后来人们推导出一些近似公式来简化Kubo公式，例如：

2、L. A. Falkovsky近似的第一种公式 ^{[4] [5]}

\sigma(\omega,\mu_{\rm c},\Gamma,T)=N_{\rm g}(\sigma_{\rm intra}+\sigma_{\rm inter}) (3a)

其中

\sigma_{\rm intra}=\frac{2{\rm i}e^2 k_{\rm B} T}{\pi \hbar^2 (\omega+2{\rm i}\Gamma)} {\rm ln}\left[ 2 {\rm cosh} \left( \frac{mu_{\rm c}}{2k_{\rm B} T} \right) \right] (3b)

\sigma_{\rm inter}=\frac{e^2}{4\hbar}\left[ G\left( \frac{\hbar\omega}{2} \right)-\frac{\rm i}{2\pi}{\rm ln}\frac{(\hbar\omega+2\mu_{\rm c})^2}{(\hbar\omega-2\mu_{\rm c})^2+(2k_{\rm B}T)^2} \right] (3c)

G(\xi)=\frac{{\rm sinh}[\xi/(k_{\rm B}T)]}{{\rm cosh}[\mu_{\rm c}/(k_{\rm B}T)]+{\rm cosh}[\xi/(k_{\rm B}T)]} (3d)

3、L. A. Falkovsky近似的第二种公式 ^{[5] [8]}

也就是将L. A. Falkovsky第一种公式中的 G\left( \hbar\omega/2 \right) 进一步简化成

G\left( \frac{\hbar\omega}{2} \right)=\frac{1}{2}+\frac{1}{\pi}{\rm arctan}\left( \frac{\hbar\omega-2\mu_{\rm c}}{2k_{\rm B} T} \right) (4)

4、G. W. Hanson近似的公式 ^[9]

\sigma(\omega,\mu_{\rm c},\Gamma,T)=N_{\rm g}(\sigma_{\rm intra}+\sigma_{\rm inter}) (5a)

其中

\sigma_{\rm intra}=\frac{{\rm i}e^2 k_{\rm B} T}{\pi \hbar^2 (\omega+2{\rm i}\Gamma)} \left[ \frac{mu_{\rm c}}{k_{\rm B}T}+2{\rm ln}\left( 1+{\rm e}^{-\frac{mu_{\rm c}}{k_{\rm B}T}} \right) \right] (5b)

\sigma_{\rm inter}=\frac{{\rm i}e^2}{4\pi\hbar}{\rm ln}\left[ \frac{2|\mu_{\rm c}|-(\omega+2{\rm i}\Gamma)\hbar}{2|\mu_{\rm c}|+(\omega+2{\rm i}\Gamma)\hbar} \right] (5c)

5、Drude-like近似 ^{[10] [11]}

\sigma(\omega,\mu_{\rm c},\Gamma,T)=N_{\rm g}\frac{e^2\mu_{\rm c}}{\pi\hbar^2}\frac{\rm i}{\omega+2{\rm i}\Gamma} (6)

这是我算出来的结果与文献 ^[12] 中的图进行对比，两者符合得非常好：

上面5种公式的Matlab计算代码 ^[13] ：https:// opt-simul.com/archives/ 691

有需要的朋友还可以从我的主页下载光学仿真案例（包括COMSOL案例和Lumerical FDTD案例），这些案例都是我原创的，大部分是SCI论文复现：

参考

1. ^ V. P. Gusynin, S. G. Sharapov and J. P. Carbotte, J. Phys.: Condens. Matter 19, 026222 (2007). https://doi.org/10.1088/0953-8984/19/2/026222
2. ^ V. P. Gusynin, S. G. Sharapov and J. P. Carbotte, Phys. Rev. B 75, 165407 (2007). https://doi.org/10.1103/PhysRevB.75.165407
3. ^ L. A. Falkovsky and A. A. Varlamov, Eur. Phys. J. B 56, 281 (2007). https://doi.org/10.1140/epjb/e2007-00142-3
4. ^ ^{a b} L. A. Falkovsky and S. S. Pershoguba, Phys. Rev. B 76, 153410 (2007). https://doi.org/10.1103/PhysRevB.76.153410
5. ^ ^{a b c} L. A. Falkovsky, Journal of Physics: Conference Series 129, 012004 (2008). https://doi.org/10.1088/1742-6596/129/1/012004
6. ^ G. W. Hanson, J. Appl. Phys. 103, 064302 (2008). https://doi.org/10.1063/1.2891452
7. ^ Supporting online material for A. Vakil, N. Engheta, Science 332, 1291 (2011). https://doi.org/10.1126/science.1202691
8. ^ Y. Yao, M. A. Kats, P. Genevet, et. al., Nano Lett. 13, 1257 (2013). https://doi.org/10.1021/nl3047943
9. ^ G. W. Hanson, J. Appl. Phys. 103, 064302 (2008). https://doi.org/10.1063/1.2891452
10. ^ N. M. R. Peres, J. M. B. Lopes dos Santos, and T. Stauber, Phys. Rev. B 76, 073412 (2007). https://doi.org/10.1103/PhysRevB.76.073412
11. ^ P. A. D Gonçalves, N. M. R Peres, An introduction to graphene plasmonics (World Scientific Publishing 2016). https://doi.org/10.1142/9948
12. ^ R. Alaee, M. Farhat, C. Rockstuhl, et. al., Opt. Express 20, 28017 (2012). https://doi.org/10.1364/OE.20.028017
13. ^ 石墨烯光学参数Matlab代码 https://www.90pan.com/b2103414