文献解读 | 柔性纳米间隙超表面：机械弯曲激活太赫兹隧穿非线性（arXiv 2026）

论文简介

这次解读的工作是 Mechanical Activation of Terahertz Tunneling in Metallic Nanogaps，作者为蔚山科学技术院（UNIST）物理系的 Dasom Kim 与 Dai-Sik Kim、仁川国立大学物理系的 Young-Mi Bahk，以及 Dukhyung Lee（现任职于穆罕默德六世理工大学），2026 年 9 月 18 日提交到 arXiv，原文链接为 https://arxiv.org/abs/2609.21506 。

一句话概括它的贡献：把金/聚甲基丙烯酸甲酯/金（Au/PMMA/Au）纳米间隙做在柔性聚对苯二甲酸乙二醇酯（PET）衬底上，用毫米级的宏观弯曲把锥形间隙的上端压到约 1 纳米，原本“造好就定死”的太赫兹隧穿非线性，就成了可以事后用机械形变打开、并用弯曲程度调档的效应。

要看清这件事的分量，先得分清两个尺度。太赫兹波在自由空间里的波长有几百微米，而电子要穿过的势垒只有几纳米。前者属于宏观电磁，后者属于量子输运，通常很难用同一个旋钮同时调节。这篇工作让它们各管一头：太赫兹脉冲负责驱动，机械弯曲负责改变势垒厚度。剩下的任务，就是从远场透射谱里把“隧穿到底有没有发生”读出来。

器件：把太赫兹场关进 30 纳米的缝里

图 1：柔性纳米间隙超表面的制备与几何结构。(a) 由方环 Au 纳米间隙构成、图形化在 250 μm 厚 PET 衬底上的超表面示意；槽尺寸 $l_x=l_y=50$ μm，阵列周期 $d_x=d_y=100$ μm，间隙宽度 $w_{\mathrm{t}}=30$ nm，Au 厚度 $h=160$ nm；入射太赫兹电场 $E_x$ 偏振方向垂直于方环的两条对边。(b) 截面电子显微图显示相邻金属面之间形成的锥形纳米间隙，标尺 100 nm。(c) 制备流程示意：在 PET 上预图形化的 Au 结构上旋涂 PMMA（绿色），随后沉积并压平第二层 Au，形成 PMMA 填充的 Au/PMMA/Au 纳米间隙。

器件读起来并不复杂。一片 250 μm 厚的 PET 上排着周期 100 μm 的方环金结构，边长 50 μm，金层厚 160 nm；方环每条边都断出一道缝，缝里填 PMMA，30 nm 的 PMMA 厚度直接定义了缝的顶部宽度。截面电子显微图里能看出这道缝是锥形的，顶部约 30 nm，越靠近底部越宽，到 80 nm 左右。

缝为什么要做到几十纳米？金属在太赫兹频段很接近理想导体，场几乎进不去金属内部，于是方环正对的两条边之间构成一个电容器。入射电场偏振垂直于这两条边时，电荷在它们表面往复堆积，缝隙里的局域场被抬得很高。算一下尺度会更直观：0.67 THz 对应自由空间波长 $\lambda=c/f\approx447.5$ μm，而单元周期只有 100 μm，约为波长的 0.22 倍。场被塞进远小于波长的体积里，才有机会去驱动纳米尺度的电子输运。

制备流程见于图 1c：先在 PET 上图形化第一层金，旋涂 PMMA 把缝宽定下来，再沉积第二层金并压平，最后用胶带把多余的金粘走。用胶带剥离听起来有点土，好处却实在——它能在宏观面积上得到均匀的金属—绝缘层—金属纳米间隙，同时保住衬底的可弯曲性。可弯曲这件事，是后面所有测量的前提。

弯曲把锥形间隙的上端压窄

图 2：柔性纳米间隙超表面的机械收窄与低场太赫兹响应。(a) 平坦态与弯曲态超表面的截面示意，$L_0=25$ mm 与 $L$ 分别为初始与弯曲后的样品长度。(b) 由 $L$ 计算得到的弯曲半径。(c) 弯曲前后纳米间隙截面的有限元结构力学模拟，色图给出间隙附近的压力分布（模拟中略去 PMMA）。(d) 由模拟得到的上端间隙宽度随弯曲半径的变化。(e) 平坦态（蓝）与弯曲态（$R=12$ mm，红）的低场太赫兹透射谱，弯曲后谐振透射被强烈抑制。

样品原始长度 25 mm。把两端往中间推，它就弯成一段圆弧，弯曲半径 $R$ 由压缩量决定。这里有一个容易忽略的细节：弯曲对样品上下两侧的作用并不对称，PET 衬底受压，贴在它上面的金层被偏心压缩，于是锥形缝的上端比下端收得更快。

结构力学有限元把这一点算了出来（图 2c、图 2d）：上端缝宽随弯曲半径减小而缩小，外推到 $R\approx16$ mm 时在几何上闭合。别忘了这个模型里没有放 PMMA。实验用到的最大弯曲是 $R\approx12$ mm，比几何闭合点还紧，却始终没有出现金属直接接触。缝隙里的 PMMA 承担了最后的间隔，几何闭合并不等于电学短路。这个区别后面还要用到。

弯曲的后果在低场透射谱上一眼可见（图 2e）。平坦样品在 0.67 THz 附近有一个谐振透射峰，来源是缝隙的电容充电；弯到 $R=12$ mm 后，这个峰被明显压下去。原因并不神秘：对向金属面靠近，间隙的等效阻抗随之改变，原来的谐振条件被破坏。到这一步，宏观形变控制纳米结构这件事已经在低场下得到验证。可是低场只看得到电容，隧穿还藏在后面。

电容旁边多出一条隧穿通道

势垒宽的时候，结的行为基本就是电容：电子波函数几乎穿不过 PMMA，电荷堆在金属表面，间隙电压跟着外加场一起涨。势垒一旦窄到纳米量级，情况就变了。波函数能从势垒里透出去一部分，产生有限的隧穿电流，于是间隙电容旁边多出一条并联的隧穿电导（图 3a）。两条通道谁占上风，取决于势垒厚度和电场强度。

势垒本身的样子也不难画（图 3b）。零偏时势垒高度 $\phi_0=4.8$ eV，加上间隙电压后，势垒被倾斜，有效高度下降，电子透射概率继续上升。宽度和电压各贡献一份非线性，隧穿通道因此“越用越通”。

作者用 Simmons 隧穿模型算出了电流密度与间隙电压的关系（图 3c），势垒宽度取了 1、3、5、10、20 nm 五档。结果分化得很清楚：1020 nm 这一档，在考察的电压范围内电流密度都小到可以忽略；压到 13 nm，电流密度提高若干数量级。把 Simmons 公式的指数因子拆开看，这个差距有多大就好理解了。指数因子近似写成 $\exp[-1.025\Delta s(\text{Å})\sqrt{\phi_0(\text{eV})}]$，取 $\phi_0=4.8$ eV：30 nm 势垒对应 $e^{-673.7}$，约 $10^{-293}$；1 nm 势垒对应 $e^{-22.5}\approx1.8\times10^{-10}$。两者差了大约 283 个数量级。所谓“平坦态几乎没有隧穿电流”，在模型里是天然成立的结果。

需要交代一句：真实的残余势垒厚度很难直接测，1 nm 这个值是作者依据对向金表面的粗糙度作的一个代表性假设。图 3c 里的电流密度是模型输出，不是测出来的数。真正被实验测到的是间隙电压和透射谱，模型负责把这两者之间的因果补上。

间隙电压为什么会长到饱和

图 4：强场太赫兹激发下的间隙电压与估算的隧穿响应。(a) 在机械弯曲纳米间隙超表面上进行强场太赫兹透射测量的示意，入射电场偏振垂直于被收窄的间隙侧边。(b) 平坦态（蓝）与弯曲态（$R=12$ mm，红）在 75（实线）、108（点线）、150 kV/cm（虚线）下的间隙电压时域波形。(c) 平坦态与弯曲态的峰值间隙电压随入射场强的变化：平坦态近似线性增长，弯曲态趋于饱和。(d) 把实验测得的峰值间隙电压代入 Simmons 模型算出的隧穿电流密度，弯曲态的局域收窄区近似为间隙宽度 1 nm 的平行 Au/PMMA/Au 结。(e) 对应的间隙有效介电常数虚部（0.67 THz）；弯曲态不断增强的耗散响应与机械收窄势垒中场驱动隧穿的开启一致。

强场实验的做法是：单周期太赫兹脉冲正入射，电场偏振垂直于被收窄的那两条边，峰值场取 75、108、150 kV/cm 三档，平坦态与强弯曲态（$R=12$ mm）各测一遍。

驱动过程可以这样想。磁场在光学厚的金层上感生面电流 $\mathbf{K}0(t)=\hat{\mathbf{n}}\times2\mathbf{H}{\mathrm{inc}}(t)$，它主要由入射脉冲决定，相当于给结外接了一个瞬态的电流源。接下来是一场拉锯：电流给对向金属面充电，间隙电压往上抬；电压一高，隧穿电流又把电荷泄掉，把电压按回去。最终观测到的间隙电压，是充电与放电两者竞争的结果。

从远场波形反推间隙电压用了一个面积平均的处理，写成

$$V_{\mathrm{gap}}(t)=E_{\mathrm{gap}}(t)w_{\mathrm{t}}=\frac{w_{\mathrm{t}}}{\beta}\cdot\frac{4n_{\mathrm{PET}}}{(1+n_{\mathrm{PET}})^2}\cdot\frac{E_{\mathrm{sam}}(t)}{E_{\mathrm{PET}}^{\mathrm{pk}}}\cdot E_0$$

其中 $n_{\mathrm{PET}}=1.78$，纳米间隙的覆盖率 $\beta=2l_yw_{\mathrm{t}}/(p_xp_y)$。代入尺寸可以验一个结论：$\beta=2\times50\ \mu\text{m}\times30\ \text{nm}/(100\ \mu\text{m})^2=3.0\times10^{-4}$，因此 $w_{\mathrm{t}}/\beta=100\ \mu\text{m}$，恰等于 $p_xp_y/(2l_y)$。估算式与间隙宽度无关。正是靠这层解耦，作者才能用名义上的 $w_{\mathrm{t}}$ 从远场波形里推出电压，不必事先知道间隙实际还剩多宽。

平坦态的电压随场强近似线性上涨，150 kV/cm 时峰值约 85 V。拿这个数估一下局域场：$85\ \text{V}/30\ \text{nm}=2.8\times10^{9}$ V/m，约 28 MV/cm，是入射峰值场 150 kV/cm 的 190 倍左右。缝隙确实把场攒了起来。

弯曲态是另一幅样子：同样的场强范围里，峰值电压只从约 10 V 升到约 12 V，几乎是平的。原因还是那场拉锯。势垒被压到 1 nm 量级，隧穿电流开始替电容泄放电荷，电压刚涨一点就被泄掉，自然长不高。把实测电压回代 Simmons 模型，平坦态（按 30 nm 计）的电流密度可以忽略，弯曲态（按 1 nm 计）则随场强明显上升。换算成间隙的有效介电常数虚部 $\mathrm{Im}{\varepsilon_{\mathrm{gap}}}=J_{\mathrm{tun}}/(\varepsilon_0\omega E_{\mathrm{gap}}^{\mathrm{pk}})$，它随入射场增大，意味着耗散通道在增强。

这里还有一个反面论证值得记一笔。假如弯曲直接把两块金压成了欧姆接触，谐振透射会在所有场强下都被压掉，并且与场强无关。实测并非如此：弯曲样品在最低场强下仍保留可测的谐振，场强升高才被进一步压制。这种场依赖关系，只能由有限厚度势垒上的隧穿给出。

用弯曲程度给非线性调档

时域和单点看完了，换到频域更能看出弯曲提供的连续性。图 5a 给出平坦态与强弯曲态的透射谱；把每个状态对各自的 75 kV/cm 谐振透射归一化之后（图 5b），三种弯曲程度排出一条清晰的序列。

图 5：机械激活的非线性太赫兹透射。(a) 由时域波形傅里叶变换得到、并以裸 PET 透射归一的纳米间隙超表面太赫兹透射谱，平坦态与弯曲态分别给出；实线、点线、虚线对应入射场强 75、108、150 kV/cm。(b) 平坦态、中等弯曲（$R=14$ mm）与强弯曲（$R=12$ mm）的归一化谐振透射随入射场强的变化（另含 97 kV/cm 一次测量）；每个弯曲条件下都以各自 75 kV/cm 的谐振透射归一，$T_{\mathrm{norm}}^{\mathrm{res}}(E_0)=T^{\mathrm{res}}(E_0)/T^{\mathrm{res}}(75\ \mathrm{kV/cm})$。

表 1：平坦态与两种弯曲程度下纳米间隙响应的对比。结构模型给出的顶部间隙在强弯曲时外推到约 1 nm，而实验中的残余势垒厚度难以直接测定；归一化谐振透射取 150 kV/cm 入射场强的数据。

状态	顶部间隙	峰值间隙电压	归一化谐振透射（150 kV/cm）	主导机制
平坦（未弯曲）	约 30 nm（PMMA 厚度）	约 85 V，随场强近似线性	约 0.9	电容充电
中等弯曲（$R=14$ mm）	明显收窄	未单独给出	约 0.7	电容为主，隧穿开始参与
强弯曲（$R=12$ mm）	结构模型外推约 1 nm	约 10→12 V，饱和	约 0.2	隧穿放电主导

这张表把三步连成了一条线：弯曲越强，势垒越薄，间隙电压越早饱和，谐振透射被压得越狠。150 kV/cm 下从约 0.9 一路降到约 0.2，降幅足以说明机械形变对非线性强度有相当直接的控制力。换个角度看，平坦态即便把场强从 75 拉到 150 kV/cm，谐振透射几乎不动；强弯曲态在同样的场强变化下掉了八成。同一片器件、同一束脉冲，差别只来自弯曲半径。

意义与展望

这项工作把两个相隔很远的尺度接了起来：一边是毫米级的机械弯曲，一边是发生在皮秒尺度、跨越纳米势垒的电子输运。有意思的是它的读出方式。远场透射谱本身就是探测器，不必给每个纳米间隙单独搭电极，整片超表面的局部结响应能被一次性读出来。对柔性太赫兹器件来说，“无接触读出加事后激活”这个组合很实用：非线性不再由制造一次定死，而是靠弯一弯打开，还能用弯曲程度连续调档。

从原理验证走向实际器件，还有几件事要办。残余势垒厚度目前只能靠界面粗糙度估一个代表值，把界面做得更平整、用更直接的实验手段把厚度测准，隧穿电流密度的绝对量级才能被可靠标定，实验与独立模型的定量对照也才有共同的分母。机械形变的长期稳定性同样需要观察，反复弯折之后间隙几何会不会漂移，决定这类器件能不能长期工作。再往前看，如果电容主导、隧穿参与、欧姆接触这三个区间能在同一片器件上连续而可逆地走通，柔性太赫兹调制器和超快非线性开关就多了一个可以设计的旋钮。

这篇论文最值得带走的一个判断，是量子输运未必需要复杂的电学接触才能被利用。结构尺度做对了，一道几十纳米宽的缝加上一次弯曲，隧穿就能从难以触及的量子效应，变成可测量、可调节的宏观响应。

参考资料

1. Dasom Kim, Dukhyung Lee, Young-Mi Bahk, Dai-Sik Kim. Mechanical Activation of Terahertz Tunneling in Metallic Nanogaps. arXiv:2609.21506（2026-09-18）. https://arxiv.org/abs/2609.21506
2. Hyoung-Taek Lee, Jeonghoon Kim, Joon Sue Lee, Mina Yoon, Hyeong-Ryeol Park. Over 30,000-fold field enhancement of terahertz nanoresonators enabled by rapid inverse design. Nano Letters 23, 11685–11692（2023）. https://arxiv.org/abs/2308.07561
3. Xiaoguang Zhao, Jacob Schalch, Jingdi Zhang, Huseyin R. Seren, Guangwu Duan, Richard D. Averitt, Xin Zhang. Electromechanically Tunable Metasurface Transmission Waveplate at Terahertz Frequencies. Optica 5, 303（2018）. https://arxiv.org/abs/1711.08003
4. John G. Simmons. Generalized Formula for the Electric Tunnel Effect between Similar Electrodes Separated by a Thin Insulating Film. Journal of Applied Physics 34(6), 1793–1803（1963）. https://doi.org/10.1063/1.1702682
5. M. A. Seo, H. R. Park, S. M. Koo, D. J. Park, J. H. Kang, O. K. Suwal 等. Terahertz field enhancement by a metallic nano slit operating beyond the skin-depth limit. Nature Photonics 3, 152–156（2009）. https://doi.org/10.1038/nphoton.2009.22
6. Wikipedia. Quantum tunnelling. https://en.wikipedia.org/wiki/Quantum_tunnelling （2026-10-01 访问）
7. Wikipedia. Terahertz metamaterial. https://en.wikipedia.org/wiki/Terahertz_metamaterial （2026-10-01 访问）