-
Quantum simulation serves as a significant tool for studying and understanding novel phenomena in the quantum world. In recent years, it has be realized that apart from quantum platforms, classical systems like photonic crystals, phononic crystals, and mechanical oscillators can also be used to simulate quantum models by analogizing the Schrödinger equation. Among these systems, electrical circuits have emerged as a promising simulation platform owing to their low cost, technological maturity, and ease of scalability, successfully simulating numerous important quantum phenomena. Meanwhile, non-Hermitian physics breaks the Hermiticity of systems’ Hamiltonians in traditional quantum mechanics, providing a fresh perspective for understanding the physics of quantum systems, particularly open quantum systems. Non-Hermitian systems, owing to their manifestation of unique phenomena absent in Hermitian systems, have become emerging research subjects in various fields of physics. However, many non-Hermitian phenomena require specialized configurations that pose relatively high technical thresholds on quantum platforms. For instance, the non-Hermitian skin effect typically requires systems to possess non-reciprocal hopping between lattice sites. Therefore, utilizing flexible electrical circuits to simulate non-Hermitian physics becomes a natural choice. This paper provides a short review of the current experimental progress in simulating non-Hermitian lattice models by using electrical circuits. It offers a brief introduction to the relevant knowledge of non-Hermitian physics, including mathematical concepts and novel phenomena, as well as the simulation theory of electrical circuits, including the mapping theory of the lattice models, the introduction of non-Hermiticity, and the measurement of physical quantities. The aim is to provide readers with a reference for better understanding or engagement in related researches, thus promoting further development in this field. 1. 引 言
量子模拟的概念最早由著名物理学家费曼提出[1], 以解决复杂量子系统无法用经典计算机模拟的问题, 从而更好地理解奇异的量子世界. 随着低温、超导等极端技术的发展, 人造量子平台(如冷原子[2-4]、离子阱[5-7]、超导量子比特[8-10]等)表现出系统纯净、可控性强等优势, 成功模拟了许多重要的量子现象. 然而, 量子模拟平台对技术条件要求苛刻且容易受环境影响而发生退相干[11], 导致实验成本很高. 近些年, 研究者们发现主导经典系统的物态方程在一定条件下可以与量子系统所遵循的薛定谔方程相对应[12], 因此, 经典系统(如光子晶体[13-18]、声子晶体[19-24]、机械振子[25-30]等)同样可以用来模拟量子现象, 并且具有成本低廉、技术成熟和扩展性强等特点. 尤其是近期兴起的经典电路系统[31,32], 原则上可以模拟任意维度和任意边界条件下具有任意格点间跃迁的量子紧束缚模型(即格点模型). 利用经典电路, 人们已经成功模拟了许多量子现象[33-49], 比如拓扑边缘态[35,39,46,48]以及高阶拓扑角态[36,40,43,47]等.
另一方面, 非厄米系统作为量子开放系统的一种有效描述[12,50,51], 本身带来许多不同于传统厄米系统的独特现象, 比如复能谱的出现、宇称-时间反演对称(parity-time-reversal symmetry, PT对称)破缺[52-56]、传统体边对应关系(bulk-boundary correspondence)的失效[57-70]、非厄米动力学[71-73]等, 已经成为当下凝聚态领域中一个新兴的研究热点. 鉴于经典电路对量子厄米拓扑系统的成功模拟, 人们自然也希望用它模拟非厄米系统, 以期更好地研究和理解新奇的非厄米物理现象. 实验上, 研究者们在利用经典电路模拟非厄米物理方面已经取得了很大进展[32], 比如成功模拟具有非互易跃迁的Su-Schriefer-Heeger (SSH)模型[41], 观测到由增益/损耗(gain/loss)诱导的非厄米拓扑边缘态[74]等. 本文将聚焦于经典电路对非厄米格点模型的模拟, 对当前的实验进展进行综述, 为相关研究提供参考, 以推动该领域进一步发展.
本文的剩余部分大致安排如下: 第2节简要介绍非厄米物理中的一些数学知识和新奇现象; 第3节介绍经典电路模拟的理论基础; 第4节概述当下经典电路模拟非厄米格点模型的实验进展; 第5节进行总结.
2. 非厄米物理简介
非厄米物理的研究对象既可以是量子系统也可以是经典系统, 其特征是系统的性质可以通过有效的非厄米矩阵进行描述[12]. 系统的非厄米性通常来源于系统与环境之间的耦合, 比如系统与环境之间的能量交换, 对系统的测量等[50,51].
早期的研究主要关注于PT对称的非厄米系统[52,53], 因为这类系统在特定参数下具有类似厄米系统的纯实数能谱, 以保证态的演化不发散或不消逝. 同时, 此类系统能谱的实复转变对应于本征态的PT对称破缺, 其转变点即为异常点(exceptional point, EP)[75]. 随后, 人们建立了非厄米系统的一般性理论—非厄米量子力学[76], 给出了描述非厄米系统的基本数学范式.
近些年, 传统体边对应关系在非厄米拓扑系统的失效引起新一轮对非厄米物理研究的浪潮. 在重建非厄米体边对应关系的过程中, 人们逐渐发现一些非厄米系统所特有的现象, 比如非厄米趋肤效应[57-59,77-79], 也建立起一些新的非厄米理论, 比如非布洛赫理论[59,80-88]、非厄米拓扑分类[89-92]等. 这些研究在理论上揭示出非厄米系统不同于厄米系统的独特性质[12,86,92]. 近期随着非厄米领域的不断发展, 人们开始将非厄米理论应用于对开放系统的研究中[63,93-101].
本节根据理解相关实验的需要, 简要介绍一些实验中所涉及的非厄米理论的基本数学知识以及非厄米系统所特有的新奇现象.
2.1 非厄米理论的数学基础
2.1.1 非厄米矩阵及双正交基
非厄米系统通常可以用非厄米矩阵H (相当于传统量子力学中系统的哈密顿量)来描述, 其非厄米性表现为
$ H^\dagger \neq H $ . 对于可对角化的非厄米矩阵而言(不可对角化的情况随后介绍), 其本征值分解如下[102]:$$ {S}^{-1}{HS}={\varLambda}, $$ (1) 其中, Λ为对角矩阵, 其对角项
$ \{E_n\} $ 为本征值(相当于传统量子力学中系统的本征能量), 可以为任意复数; S为相似矩阵, 其中的列向量被称为H的右本征矢, 记作$ \{|\psi_n^{(\mathrm{r})}\rangle\} $ , 而${S}^{-1} $ 中的行向量被称为H的左本征矢, 记作$ \{\langle \psi_n^{(\mathrm{l})}|\} $ . 将(1)式写成本征方程的形式:$$ H|{\psi} _{n}^{(\mathrm{r})}\rangle = E_n |{\psi} _{n}^{(\mathrm{r})}\rangle,\; \; H^\dagger|{\psi} _{n}^{(\mathrm{l})}\rangle = E^*_n |{\psi} _{n}^{(\mathrm{l})}\rangle, $$ (2) 表示非厄米系统的定态薛定谔方程. 由相似矩阵的性质可知, 左右本征矢之间满足双正交归一关系:
$$ \langle {\psi} _{m}^{(\mathrm{l})}|{\psi} _{n}^{(\mathrm{r})}\rangle = \delta _{mn}, $$ (3) 从而具有如下完备性:
$$ \sum\limits_n | {\psi _n^{({\rm{r}})}} \rangle \langle {\psi _n^{({\rm{l}})}} | = \mathrm{1}. $$ (4) 因此, 左右矢可构成非厄米矩阵的双正交基(biorthogonal basis)[76].
特别地, 对于厄米矩阵(
$ {H}^{\dagger}={H}$ )而言, 其本征值分解中的相似矩阵变为酉矩阵($ {S}^{-1}={S}^{\dagger} $ ), 即左右矢之间满足$ | {\psi _n^{({\rm{l}})}} \rangle =| {\psi _n^{({\rm{r}})}} \rangle $ 的关系, 因此双正交基退化为大家所熟知的正交基, 相应的本征值$ E_n $ 也变为实数. 不同于厄米矩阵的本征矢, 非厄米矩阵对应不同本征值的右本征矢通常不正交(左本征矢类似), 即$ \langle {\psi} _{m}^{(\mathrm{r})}|{\psi} _{n}^{(\mathrm{r})}\rangle \neq 0 \; (m\ne n) $ .2.1.2 缺陷矩阵及EP
当非厄米矩阵不可对角化(被称为缺陷矩阵)时, 无法对其进行本征值分解, 取而代之的是更一般的约当分解(Jordan decomposition)[12,103]:
$${S}^{-1}{HS}={J}. $$ (5) 其中, S仍为相似矩阵,
$$ {J}=\begin{pmatrix}\ddots & & 0 \\ & J_s(E) & \\ 0 & & \ddots\end{pmatrix} $$ (6) 是在相似变换下最接近完全对角化的块对角矩阵, 其对角块
$ J_{s}(E) $ 具有如下形式: 对角元均为E, 上次对角元均为1, 其他为0. 具有这样形式的J被称为H的约当标准型(Jordan canonical form), 其中$ J_{s}(E) $ 为第s个约当块(Jordan block). 每个约当块$ J_s(E) $ 都有且仅有一个本征矢$ (1, 0, \cdots)^\mathrm{T} $ , E为相应的本征值. 因此, 当存在约当块的维度$ p_{s}(E) = \dim{J_{s}(E)} > 1 $ 时, H即为缺陷矩阵. 具有相同本征值E的约当块的个数$ \alpha(E) $ 即为缺陷矩阵H在E处的简并度. 特别地, 当所有约当块的维度均为1时, 约当分解则退化为本征值分解, 即H不再是缺陷矩阵.如果调节系统参数, 使描述非厄米系统的非厄米矩阵恰好为缺陷矩阵, 此参数即为非厄米系统的EP[75]. 不同于厄米矩阵的简并点, 在EP处不仅本征值重合, 本征矢也会部分或全部合并(coalesce), 即H的左右本征矢无法提供完备的双正交基, 此时需要用广义本征矢(generalized eigenvector)来补足. 通常大家定义在本征值E处的EP阶数为
$ p(E) = \displaystyle\sum\nolimits_{s = 1}^{\alpha(E)}p_s(E) $ , 但此定义只适合对EP的粗糙描述, 并不能反映在此系统参数下是否发生了本征矢合并以及合并的细节. 只有给出本征值E对应的每一个约当块的维度信息$ p_s(E) $ , 才能更好地表征EP的性质.为了度量系统参数离EP的远近程度, 可以定义平均的相刚度(phase rigidity)[104]:
$$ r(z) = \frac{1}{N} \sum\limits_{n = 1}^{N} \frac{\langle \psi_n^{(\mathrm{l})}(z) | \psi_n^{(\mathrm{r})}(z) \rangle }{\sqrt{\langle \psi_n^{(\mathrm{l})}(z) | \psi_n^{(\mathrm{l})}(z) \rangle\langle \psi_n^{(\mathrm{r})}(z) | \psi_n^{(\mathrm{r})}(z) \rangle}}, $$ (7) 其中, N为可对角化矩阵
$ H(z) $ 的维度,$ | \psi_n^{(\mathrm{l/r})}(z)\rangle $ 为其在参数z处的第n个左/右本征矢. 当系统参数接近EP值$ z_0 $ 时, 相刚度$ r(z\rightarrow z_0)\rightarrow 0 $ ; 对于厄米矩阵, 左右本征矢互为复共轭, 所以相刚度$ r(z) = 1 $ , 即厄米矩阵不存在EP.2.2 非厄米系统中的新奇现象
2.2.1 PT对称破缺及赝厄米
在传统的量子力学中, 系统的哈密顿量为厄米算符(对应于厄米矩阵), 其本征谱全为实数, 反映系统能量为实数的物理事实. 1998年, Bender和Boettcher[52,53]发现, PT对称的非厄米哈密顿量(对应于非厄米矩阵)同样可以具有全实能谱. 这里的P和T分别表示空间反演和时间反演. 当非厄米强度(比如增益/损耗强度) γ比较弱时, 系统的所有本征态均具有PT对称性, 从而具有全实的能谱, 此时系统处于PT对称相, 任何量子态在此系统下均具有稳定的动力学. 当非厄米强度很强时, 系统的部分或全部本征态不再具有PT对称性, 其能谱也出现复数, 此时系统处于PT对称破缺相, 量子态在其中的演化通常会发散或消逝. 因此, 在PT对称的非厄米系统里, 存在能谱由全实到复数的转变, 被称为PT转变(PT transition), 如图1所示. 在转变点处部分或全部本征态会发生自发性PT对称破缺. 此转变点即为前面所提到的EP. PT转变是厄米系统所没有的.
其实, 不只有PT对称的非厄米系统可以存在全实能谱. Mostafazadeh[105]在2002年的研究发现, 每一个具有全实能谱的哈密顿量都是赝厄米的(pseudo-Hermitian), 而PT对称只是赝厄米的一种特殊形式. 如果存在一个厄米的可逆算符η, 使得系统的哈密顿量H满足
$ H^{\dagger} = \eta H \eta^{-1} $ , 则称H是赝厄米的[12,105]. 如果η可以取单位算符, 则H退化为厄米的. 赝厄米哈密顿量的能谱一定具备以下性质之一: 1) 全实能谱; 2) 能谱以复共轭的形式成对出现, 且互为复共轭的能谱的简并度相同.如果存在η, 它可以进一步写成
$ \eta = OO^\dagger $ 的形式, 其中O为线性可逆算符, 则赝厄米哈密顿量H一定具有全实能谱[12,106]. 这是非厄米哈密顿量具有全实能谱的充分必要条件. 而PT对称仅能保证非厄米哈密顿量具有产生全实能谱的可能性, 既不是其具有全实能谱的充分条件, 也不是必要条件.2.2.2 传统体边对应关系的失效及非厄米趋肤效应
体边对应关系是被体能隙保护的厄米拓扑系统遵循的一个基本原则, 它描述了系统的体态拓扑不变量与拓扑边缘态之间的关联. 然而, 此原则在某些非厄米系统中并不成立, 表现为开边界条件下的能谱和体态与周期边界下有很大的不同[57,58], 如图2(a)所示. 这是因为此类非厄米系统对边界的选择表现出很强的敏感性, 开边界与周期边界的同一系统在热力学极限下并不等价[57,58,107].
非厄米趋肤效应恰是这种边界敏感性的体现. 它具体表现为非厄米系统的体态在开边界条件下呈指数型地聚集在某一边界, 如图2(c)所示. 非厄米趋肤效应使具有体周期性的非厄米系统的体态丢失了布洛赫态的特性, 与厄米系统中体态弥散在全域的情形完全不同.
基于对非厄米趋肤效应的观察, Yao等[60]在一维非厄米SSH模型中建立起非布洛赫理论(non-Bloch theory), 成功重塑了此非厄米系统中的体边对应关系, 引起了后续的广泛研究[80-88]. 在非布洛赫理论中, 对于热力学极限下的非厄米系统, 周期边界下的哈密顿量
$ H(k) $ 被开边界下的$ H(\beta) $ 所取代. 这里, 原先定义在布里渊区上的晶格动量k被扩展为一个复变量$ \beta = r {\rm e}^{{\rm i}k} $ , 它在复平面的集合被称为广义布里渊区(generalized Brillouin zone). 广义布里渊区通常为一个闭合路径, 如图2(b)所示, 其与原点的距离r反映系统在开边界条件下体态的趋肤性质:$ r = 1 $ 代表布洛赫态, 即没有非厄米趋肤效应;$ r < 1 $ 和$ r > 1 $ 分别对应趋向于不同边界的体趋肤态[60,81]. 类比厄米情形[108], 可以用非布洛赫态在广义布里渊区中定义非布洛赫的拓扑不变量, 从而重构非厄米拓扑系统的体边对应关系, 即非布洛赫拓扑不变量与非厄米拓扑边缘态之间的关联: 拓扑不变量为0表示不存在拓扑边缘态的拓扑平庸相, 非0表示存在拓扑边缘态的非厄米拓扑相.例如, 对于具有手征对称性的一维非厄米系统, 其非布洛赫缠绕数可定义为[59]
$$ w = \frac{1}{2 \pi \mathrm{i}} \oint_{\rm{GBZ}} \mathrm{tr}[{q}^{-1}(\beta) {\rm d} { q}(\beta)], $$ (8) 其中, 矩阵
${q}(\beta) $ 由哈密顿量$ H(\beta) $ 的Q-矩阵来定义:$$ \begin{split} {Q}(\beta) =\;& \sum\limits_{n\in N_-} \left[\varGamma|u^{\rm{(r)}}_n(\beta)\rangle\langle u^{\rm{(l)}}_n(\beta)|\varGamma\right.\\ &\left.-|u^{\rm{(r)}}_n(\beta)\rangle\langle u^{\rm{(l)}}_n(\beta)|\right] \\ =\;& \left(\begin{array}{ccc} 0 & {q}(\beta) \\ {q}^{-1}(\beta) &0 \end{array}\right). \end{split} $$ (9) 由于系统具有手征对称性
$ \{H(\beta), \varGamma\} = 0 $ (Γ为相应的手征算符), 导致能谱关于零点对称, 因此可以将$ H(\beta) $ 的本征态$ \{|u^{\rm{(r, l)}}_n(\beta)\rangle\} $ (n为能带指标)按本征能量$ \{E_n(\beta), -E_n(\beta)\} $ 划分成两个子空间$ N_+ $ 和$ N_-$ . 以手征算符Γ的本征态为基, 手征对称的Q-矩阵便可写成(9)式中第2行的反对角形式.又如, 二维非厄米系统的非布洛赫陈数可定义为[60,109]
$$ C = \frac{1}{2 \pi} \oint_{\rm{GBZ}} \mathrm{tr}[F_{xy}({\boldsymbol{\beta}})] {\rm d}\beta_x {\rm d}\beta_y, $$ (10) 其中(
$ \partial_{\mu = x, y}\equiv\partial/\partial \beta_\mu $ ),$$ \begin{split} F_{xy}({\boldsymbol{\beta}}) =\;& \partial_x A_y({\boldsymbol{\beta}})-\partial_y A_x({\boldsymbol{\beta}}) \\ & -{\rm i}[A_x({\boldsymbol{\beta}}),A_y({\boldsymbol{\beta}})], \end{split} $$ (11) $$ A_\mu^{nm}({\boldsymbol{\beta}}) = \mathrm{i}\langle u_n^{({\rm{l}})}({\boldsymbol{\beta}}) | \partial_\mu u_m^{({\rm{r}})}({\boldsymbol{\beta}})\rangle $$ (12) 分别为非布洛赫贝里曲率(Berry curvature)和贝里联络(Berry connection).
以上用非布洛赫本征态定义的拓扑不变量, 虽然能很好地反映非厄米系统在开边界条件下的拓扑相变, 包括能隙的关闭以及拓扑边缘态的产生, 但并没有体现出非厄米趋肤态本身的拓扑性质. 借助于非厄米系统的能谱一般为复数的特性, 可以定义能量缠绕数[12,67]:
$$ w_{C,E_{{\rm{b}}}} = \frac{1}{2 \pi \mathrm{i}} \oint_{C} \frac{\rm d}{{\rm d} z} \ln\det \left[H(z)-E_ \mathrm{b}\right] {\rm d} z, $$ (13) 其中, C为积分回路, 可以是布里渊区(BZ)或广义布里渊区(GBZ), 也可以是其他周期参数空间;
$ E_\mathrm{b} $ 是能量复平面内的基准能量. 对于一维单带非厄米系统, 周期边界下的能量缠绕数$ w_{{\rm{BZ}}, E_{{\rm{b}}}} $ 可以反映相应开边界下趋肤态的性质:$ w_{{\rm{BZ}}, E_{{\rm{b}}}}\ne 0 $ 表明存在非厄米趋肤效应, 其正负反映趋肤态的趋肤方向[67]. 这种周期边界下能谱的拓扑性质与开边界下非厄米趋肤态的对应关系是非厄米拓扑系统所独有的.2.2.3 EP诱导的高灵敏度及分数级数展开
与传统的简并点(此处仅是能量重合而态并不合并)不同, EP能使非厄米系统的能量对微扰产生更灵敏的响应. 对于一个p阶EP, 能谱对微扰z的响应ε在一定条件下最大可以达到
$ \varepsilon\propto z^{1/p} $ 的量级, 而传统的简并点仅为$ \varepsilon\propto z $ [110], 如图3所示. 基于这一特性, EP可以用于制造高灵敏度的传感器, 其最早方案由Wiersig[111]于2014年提出.EP产生高灵敏度的原因可以通过数学上的皮瑟级数(Puiseux series)来理解. 已知非厄米系统
$ H(z) $ (参数$ z\in {\mathbb{C}} $ )在$ z = 0 $ 处有一个能量为$ E_0 $ 的p阶EP. 当系统偏离EP时, 假设能量完全劈裂为p支不同的能量函数$ E_h(z)\; (h = 0, \cdots, p-1) $ 且它们恰好构成以$ z = 0 $ 为$ p-1 $ 阶支点的p叶黎曼面(即这p支能量函数之间在$ z = 0 $ 附近是连续解析的, 且函数值绕p圈才能回到初始值), 则能量函数在EP附近展开有如下皮瑟级数[92,112]:$$ E_h(z) = E_0+\sum\limits_{n = 1}^{\infty} \alpha_n {\rm e}^{{\rm i}2\pi n h/p}z^{n/p}, $$ (14) 其中
$ \alpha_n $ 是展开系数. 因为$ |z|\ll 1 $ , 所以能量在EP附近的劈裂由领头项$ z^{1/p} $ 主导. 皮瑟级数是分数级数, 当EP阶数$ p > 1 $ 时, 能量劈裂对参数偏离的响应比传统的简并点(对应于泰勒级数)要大. 这就是EP能导致高灵敏度的来源.值得注意的是, 展开式(14)假设了p阶EP恰好是其附近能量函数的
$ p-1 $ 阶支点, 但这些能量函数也可以形成α个各自连续解析的函数族, 每族函数以p阶EP作为其$ p_s-1 $ 阶支点(满足$ \displaystyle \sum\nolimits_{s = 1}^\alpha p_s = p $ ). 此时, 每一族能量函数的展开式都具有(14)式的形式, 只是将p换成$ p_s $ , 这会导致能量对参数偏离的响应相对较弱[103]. 因此, 能谱对微扰的响应并非只由EP的阶数决定, 还跟具体的微扰形式有关.3. 经典电路模拟的理论基础
3.1 经典电路对格点模型的映射
利用经典电路模拟物理现象的基本逻辑为, 基于电路的基尔霍夫定律, 通过合理近似, 建立起电路中描述某物理量的运动方程与模拟对象所遵循的运动方程之间的映射关系, 以达到通过观测此物理量在电路中的行为从而获知模拟对象相关性质的目的. 由于基尔霍夫方程的离散属性, 经典电路主要用于模拟量子力学中的紧束缚模型(即格点模型). 所以, 本节主要介绍用于映射格点模型薛定谔方程的3种电路理论: 拉普拉斯形式(Laplacian formalism)[35]、刘维尔形式(Liouvillian formalism)[113]以及耦合模理论(coupled mode theory)[114-116].
3.1.1 拉普拉斯形式
对于任意线性电路构成的网络, 元件从节点i到节点j (有方向性)的特性可以通过线性的等效电阻
$ R_{ij} $ 、等效电感$ L_{ij} $ 和等效电容$ C_{ij} $ 来描述. 如果元件的等效值依赖于节点i与j之间的方向, 则称此元件为非互易的(nonreciprocal). 定义节点的外部输入电流$ I_i(t) $ 和对地电压$ V_i(t) $ , 根据基尔霍夫定律可得矩阵形式的线性微分方程[117]:$$ \frac{{\rm d}}{{\rm d}t}{\cal{I}}(t) = {\cal{C}} _\mathrm{C} \frac{{\rm d}^2}{{\rm d}t^2}{\cal{V}}(t) +{\cal{C}} _\mathrm{R} \frac{{\rm d}}{{\rm d}t}{\cal{V}}(t) +{\cal{C}} _\mathrm{L} {\cal{V}}(t), $$ (15) 其中,
$ {\cal{I}}(t) $ 和$ {\cal{V}}(t) $ 分别为$ I_i(t) $ 和$ V_i(t) $ 的矢量形式, 系数矩阵:$$ \begin{split} & ({\cal{C}}_\mathrm{C})_{ii} = C_{i{\rm{g}}}+\sum\limits_{k\ne i}C_{ik},\; \; \, ({\cal{C}}_\mathrm{C})_{ij} = -C_{ij}\; (i\ne j), \\ &({\cal{C}}_\mathrm{R})_{ii} = R_{i{\rm{g}}}^{-1}+\sum\limits_{l\ne i}R^{-1}_{ik},\; ({\cal{C}}_\mathrm{R})_{ij} = -R^{-1}_{ij}\; (i\ne j), \\ &({\cal{C}}_\mathrm{L})_{ii} = L_{i{\rm{g}}}^{-1}+\sum\limits_{k\ne i}L^{-1}_{ik},\; \; ({\cal{C}}_\mathrm{L})_{ij} = -L^{-1}_{ij}\; (i\ne j), \end{split} $$ (16) 分别具有 [电容]、[电阻]–1和 [电感]–1的量纲, g表示对地.
对电路注入驱动频率为
$ \omega\in{\mathbb{R}} $ 的交变电流$ {\cal{I}}(t) = {\cal{I}}{\rm e}^{\mathrm{{\rm i}}\omega t} $ , 电路节点将具有稳定的电压响应$ {\cal{V}}(t) = {\cal{V}} {\rm e}^{\mathrm{{\rm i}} \omega t} $ , 其电流幅$ {\cal{I}} $ 和电压幅$ {\cal{V}} $ 之间满足如下关系:$$ {\cal{I}} = \left(\mathrm{i} \omega {\cal{C}} _\mathrm{C}+{\cal{C}}_\mathrm{R} -\mathrm{i} {\cal{C}} _\mathrm{L}/\omega \right){\cal{V}} \equiv{\cal{J}}(\omega){\cal{V}}. $$ (17) 式中
$ {\cal{J}}(\omega) $ 具有 [导纳] 的量纲, 被称为导纳矩阵(admittance matrix)或电路拉普拉斯量(circuit Laplacian)[35]. 它是驱动频率ω的函数, 其矩阵元$ {\cal{J}}_{ij} $ 表示节点i到j的导纳. 其逆矩阵$ {\cal{{\boldsymbol G}}} = {\cal{{\boldsymbol J}}}^{-1} $ 被称为电路格林函数(circuit Green function), 具有[阻抗] 的量纲, 反映了节点电压$ {\cal{V}} $ 对节点输入电流$ {\cal{I}} $ 的响应. 有时为了方便起见, 也可将导纳矩阵分解为三部分:$$ {\cal{{\boldsymbol J}}} = {\cal{{\boldsymbol A}}}+{\cal{{\boldsymbol D}}}+{\cal{{\boldsymbol W}}}, $$ (18) 其中,
$ {\cal{A}} $ 代表导纳矩阵的非对角部分,$ {\cal{D}} $ 和$ {\cal{W}} $ 分别代表对角部分中的节点间耦合部分和对地部分, 其矩阵元分别为$$ \begin{split} &{\cal{A}}_{ij} = (1-\delta_{ij})\left[-\mathrm{i} \omega C_{ij}-R_{ij}^{-1}-(\mathrm{i} \omega L_{ij})^{-1}\right], \\ &{\cal{D}}_{ij} = \delta_{ij}\sum\limits_{k\ne i}\left[\mathrm{i} \omega C_{ik}+R_{ik}^{-1}+(\mathrm{i} \omega L_{ik})^{-1}\right], \\ &{\cal{W}}_{ij} = \delta_{ij}\left[\mathrm{i} \omega C_{i{\rm{g}}}+R_{i {\rm{g}}}^{-1}+(\mathrm{i} \omega L_{i{\rm{g}}})^{-1}\right], \end{split} $$ (19) 其中
$ \delta_{ij} $ 为Kronecker函数.如果将导纳矩阵
$ {\cal{J}} $ 看作格点模型的哈密顿量, 则电路节点对应格点的位置,$ {\cal{J}} $ 的对角元和非对角元分别对应格点上的在位势能和格点间的跃迁振幅, 而导纳矩阵的本征方程,$$ {\cal{J}}v^{({\rm{r}})}_n = j_nv_n^{({\rm{r}})}, \; \; {\cal{J}}^\dagger v^{({\rm{l}})}_n = j^*_n v_n^{({\rm{l}})}, $$ (20) 便可用于模拟格点模型的定态薛定谔方程(2)式. 这里的
$ j_n $ 表示导纳矩阵$ {\cal{J}} $ 的第n个本征值,$ v_n^{({\rm{l, r}})} $ 为相应的左右本征矢. 因此,$ j_n $ 对应格点模型的本征能量$ E_n $ , 而$ v_n^{({\rm{l, r}})} $ 对应格点模型的本征态$ |\psi_n^{({\rm{l, r}})}\rangle $ . 需要注意的是, 在此对应关系中, 驱动频率ω是作为模拟参数而存在的, 并非直接对应于格点模型的本征能量.特别地, 当无外界输入电流(即
$ {\cal{I}} = 0 $ )时, (17)式变为本征方程$ {\cal{J}}(\omega){\cal{V}} = 0 $ , 其存在非平庸解的条件为$$ \det\left[{\cal{J}}(\widetilde{\omega})\right] = \prod\nolimits_n j_n(\widetilde{\omega}) = 0,$$ (21) 其中使条件成立的
$ \widetilde{\omega} $ 即为电路的本征频率, 非平庸解$ {\cal{V}}(\widetilde{\omega}) $ 为本征频率$ \widetilde{\omega} $ 所对应的电路的本征电压. 条件(21)式说明本征频率$ \widetilde{\omega} $ 至少使一个导纳矩阵的本征值为零, 即存在n使$ j_{n}(\widetilde{\omega}) = 0 $ . 换而言之, 导纳矩阵$ {\cal{J}}(\omega) $ 的本征值$ j_n(\omega) $ 作为ω的连续函数族, 与$ j = 0 $ 的交点所对应的频率即为电路的本征频率$ \widetilde{\omega} $ , 而相应的右本征矢$ v_n^{({\rm{r}})}(\widetilde{\omega}) $ 即为电路的本征电压. 注意以上结论并不限于本征频率是实数, 对于复本征频率同样适用.由(17)式可知, 当电路中仅含有被动(passive)电阻时(即纯耗散电路),
$ {\cal{J}} $ 是厄米的, 而仅含有被动电容和电感时(即纯振荡电路),$ {\cal{J}} $ 是反厄米的, 即$ {\cal{J}}^\dagger = -{\cal{J}} $ . 因此, 通常用$ \mathrm{i} {\cal{J}} $ 建立与格点模型哈密顿量H的映射, 用仅含有被动电容和电感的电路模拟厄米格点模型[33,35,45], 而通过引入被动电阻以及主动(active)元件来模拟非厄米格点模型[41,74,118-129]. 只要对电路元件及排布进行合理的设计并采取合适的驱动频率, 原则上可以用导纳矩阵模拟任意格点模型的稳态性质, 包括本征能和本征态以及由它们定义的各种物理量.3.1.2 刘维尔形式
如果要模拟格点模型的动力学, 需要将描述电路动力学的(15)式与描述格点模型动力学的含时薛定谔方程
$ \mathrm{i}\partial_t |\psi(t)\rangle = H|\psi(t)\rangle $ ($ \hbar = 1 $ )相联系. 由于含时薛定谔方程是齐次方程, 所以这里考虑无外界电流输入的情况, 即$ {\cal{I}}(t) = 0 $ . 此时, (15)式变为关于时间的二阶齐次微分方程而含时薛定谔方程是一阶的, 因此, 需要通过变量替换将(15)式降为一阶:$$ \begin{split} &\frac{{\rm d}}{{\rm d}t}\; {\cal{V}}(t) \equiv {\cal{W}}(t), \\ &\frac{{\rm d}}{{\rm d}t}{\cal{W}}(t) = -{\cal{C}} _\mathrm{C}^{-1} {\cal{C}} _\mathrm{R} {\cal{W}}(t) -{\cal{C}} _\mathrm{C}^{-1} {\cal{C}} _\mathrm{L} {\cal{V}}(t). \end{split} $$ (22) 定义新变量
$ \varPhi(t) = \left[{\cal{V}}(t), {\cal{W}}(t) \right]^{\rm{T}} $ , (22)式可改写为如下刘维尔方程的形式[113]:$$ \frac{{\rm d}}{{\rm d}t}\varPhi(t) = {\cal{L}}\varPhi(t), \; \; {\cal{L}} = \begin{pmatrix} 0&{\cal{E}} \\ -{\cal{C}} _\mathrm{C}^{-1} {\cal{C}} _\mathrm{L}&-{\cal{C}} _\mathrm{C}^{-1} {\cal{C}} _\mathrm{R} \end{pmatrix}. $$ (23) 这里的
$ {\cal{L}} $ 被称为电路刘维尔量(circuit Liouvillian),$ {\cal{E}} $ 为单位矩阵. 刘维尔方程将微分方程(15)从二阶降为一阶, 但变量的个数扩大为原来的2倍.显然, 电路的刘维尔方程类似含时薛定谔方程,
$ \mathrm{i}{\cal{L}} $ 具有哈密顿量的地位. 此方程具有形为$ \varPhi(t) = \varPhi {\rm e}^{-\mathrm{{\rm i}} \widetilde{\omega} t} $ 的稳态解, 其中$ \mathrm{i}{\cal{L}} $ 的本征值$ \widetilde{\omega} $ 即为电路的本征频率, 模拟格点模型的本征能$ E_n $ [122,130]; 矢量$ \varPhi(t) $ 反映电路中电压随时间的演化情况, 模拟格点模型中量子态的动力学[131].3.1.3 耦合模理论
因为格点模型直接来源于凝聚态物理中的紧束缚近似(即将不同的原子轨道弱耦合在一起), 所以同样可以用电路重构类似的近似过程: 将多个具有独自振荡频率的RLC电路, 通过某种形式的弱耦合连接起来, 从而模拟格点模型. 这种近似方法被称为耦合模理论[114-116].
为了清楚地展示耦合模理论的近似过程, 这里选取通过电感元件间的互感M耦合的两个RLC回路作为例子[115], 如图4所示. 根据基尔霍夫电压定律, 可得
$$ \mathrm{i}\omega (L I_{\rm{l,r}}+MI_{\rm{r,l}})+ R_{\rm{l,r}}I_{\rm{l,r}}+\frac{I_{\rm{l,r}}}{ {\rm{i}}\omega C} = 0, $$ (24) 其中, 电阻
$ R_{\rm{l, r}} = \mp R $ (关于负电阻的实现可参见3.2节), 且设参数$ R, L, C, M $ 均为正数. 定义特征量$ \omega_0 = 1/\sqrt{LC} $ ,$ \gamma = R\sqrt{C/L} $ 和$ \mu = M/L $ , 则(24)式可改写为$$ \begin{pmatrix} 1-\dfrac{\omega_0^2}{\omega^2}+\mathrm{i} \dfrac{\gamma\omega_0}{\omega}& \mu\\ \mu & 1-\dfrac{\omega_0^2}{\omega^2}-\mathrm{i} \dfrac{\gamma\omega_0}{\omega} \end{pmatrix} \begin{pmatrix} I_{\rm{l}}\\ I_{\rm{r}} \end{pmatrix} = 0. $$ (25) 当考虑弱耦合的情况时, 即
$$ \mu\ll 1, \; \; \gamma\ll 1,\; \; \frac{\omega-\omega_0}{\omega_0}\ll1, $$ (26) (25)式可化为
$$ \begin{pmatrix} \omega-\omega_0+\mathrm{i}\dfrac{\gamma \omega_0}{2} & \dfrac{\mu \omega_0}{2}\\ \dfrac{\mu\omega_0}{2} & \omega-\omega_0-\mathrm{i}\dfrac{\gamma \omega_0}{2} \end{pmatrix} \begin{pmatrix} I_{\rm{l}}\\ I_{\rm{r}} \end{pmatrix} = 0. $$ (27) 此方程等价于振荡电流
$ I_{\rm{l, r}}(t) = I_{\rm{l, r}} {\rm e}^{\mathrm{{\rm i}} \omega t} $ 所满足的耦合模方程:$$ \mathrm{{\rm i}}\frac{{\rm d}}{{\rm d}t} \begin{pmatrix} I_{\rm{l}}(t)\\ I_{\rm{r}}(t) \end{pmatrix} = {\cal{H}} \begin{pmatrix} I_{\rm{l}}(t)\\ I_{\rm{r}}(t) \end{pmatrix}, $$ (28) 其中,
$$ {\cal{H}} = \omega_0 \begin{pmatrix} -1+\mathrm{i}\gamma /2 & \mu /2\\ \mu/2 & -1-\mathrm{i}\gamma/2 \end{pmatrix}. $$ (29) 显然, (28)式与含时薛定谔方程形式一致, 可以用来模拟相应的动力学以及稳态性质.
值得注意的是, 电路的耦合模方程(28)与刘维尔方程(23)的形式相同, 但却没有增加变量的个数, 这是因为耦合模理论在弱耦合近似下将(25)式中ω的平方项降成了(27)式中的线性项, 等价于将基尔霍夫方程(15)中关于时间的二阶导数降到了一阶. 这便是耦合模理论能直接模拟含时薛定谔方程的原因. 由于弱耦合近似(26)式的限制, 耦合模理论适用范围较窄. 目前在电路系统中主要用于高灵敏EP传感器的设计[114-116,130].
3.2 非厄米在经典电路中的引入
为了用经典电路模拟非厄米格点模型, 需要在以上电路理论中引入等效的非厄米项, 主要包括格点上的复在位势和格点间的非互易跃迁, 分别对应于哈密顿量矩阵中的复对角元和不满足复共轭关系的非对角元, 它们均可以使描述系统的矩阵失去厄米性.
3.2.1 损耗和增益
由于电阻天然的耗散属性, 在电路中引入非厄米最自然且最简单的方式就是加入电阻. 根据电路的拉普拉斯形式, 在电路中添加电阻或等效的负电阻可以模拟非厄米格点模型中的损耗或增益, 使得
$ \mathrm{i} {\cal{J}}\neq \left(\mathrm{i} {\cal{J}}\right)^{\dagger} $ . 具体地, 由(16)式和(17)式可知, 节点对地的电阻$ R_{n{\rm{g}}} $ 使$ \mathrm{i} {\cal{J}} $ 的对角元产生虚部, 可以模拟格点的复在位势; 节点间的电阻$ R_{nm} $ 则会在对角元和非对角元上同时引入虚部, 额外模拟格点间的复跃迁[118].由此可知, 电阻的引入产生损耗, 但如要产生增益, 需引入等效的负电阻. 文献[133]给出了利用运算放大器(简称运放)实现一般负阻抗(负电阻是其特殊情况)的方案, 如图5(a)所示, 其中包括对地端口(上)和两端口间(下)的负阻抗方案, 其阻抗值可由基尔霍夫电流定律分别求得
$$ Z_{i{\rm{g}}}\equiv\frac{V_i}{I_i} = -Z,\; Z_{ij}\equiv\frac{V_{i}-V_{j}}{I_{i}} = -Z = Z_{ji}, $$ (30) 其中, Z为需要变负值的目标阻抗. 图中运放的等效增益函数可用简单的加法器和乘法器实现. 应用负电阻的典型实例是对PT对称破缺的电路模拟[113,134,135].
3.2.2 非互易跃迁
相比于增益/损耗, 非互易跃迁的实现更具有挑战性. 在交流源的驱动下, 被动元件(如电阻、电容和电感)的阻抗/导纳值都与测量方向无关, 即被动元件是互易的(reciprocal). 为了使电路元件具有非互易性, 从而实现对格点间非互易跃迁的模拟, 必须引入主动元件打破其互易性.
实验上常见的方案是使用电流型负阻抗变换器(negative impedance converter with current inversion, INIC)[136], 它由运放和若干线性元件组成, 其原理如图5(b)所示. 由基尔霍夫定律易得
$$ \begin{pmatrix} I_{i}\\ I_{j} \end{pmatrix} = Z^{-1}\begin{pmatrix} -\beta& \beta \\ -1&1 \end{pmatrix} \begin{pmatrix} V_{i}\\ V_{j} \end{pmatrix}, $$ (31) 其中
$ \beta = Z_-/Z_+ $ . 当$ \beta = 1 $ 时, 不同端口测量的导纳反号, 即$ Y _{ij}\equiv Z^{-1}_{ij} = I_{i}/(V_{i}-V_{j}) = -Z^{-1} = -Y_{ji} $ . 将此INIC与被动元件X (设导纳为$ Y_{\rm{X}} $ )并联后再置于节点间, 便实现非互易跃迁的效果, 即$ Y'_{ij} = Y_{\rm{X}}-Z^{-1} $ 和$ Y'_{ji} = Y_{\rm{X}}+Z^{-1} $ . 在用经典电路模拟格点模型方面, INIC最早用于打破时间反演对称性[137], 之后被应用到模拟具有非互易跃迁的非厄米格点模型的实验中, 观测到了与非厄米趋肤效应相关的现象[41,119,121,123,125,128,138]. 当INIC对地连接时, 其效果与负电阻类似, 也常被用作模拟具有复在位势能的非厄米模型[74,114,122].此外, 电压跟随器(voltage follower)[136]也可以用来实现非互易跃迁, 如图5(c)所示. 根据理想运放正负两输入端口的电压相同的特点, 当负输入端口与输出端口相连时, 输出端口的电压
$ V_0 $ 则与正输入端口的电压$ V_{j} $ 相同, 从而实现电压跟随. 由此可得,$ Y_{ij} = I_{i}/(V_{i}-V_{j}) = Z^{-1} $ . 又由于理想运放的输入阻抗趋于无穷大, 导致$ I_{j}\sim0 $ , 所以$ Y_{ji} = I_{j}/(V_{j}-V_{i})\sim 0 $ . 显然,$ Y_{ij}\ne Y_{ji} $ . 将电压跟随器与被动元件并联同样能实现非互易跃迁的效果[124,139], 尤其是对单向跃迁(unidirectional hopping)的实现[129], 相比INIC更为简单.3.3 经典电路中物理量的测量与表征
建立好经典电路对非厄米格点模型的映射后, 则需要考虑如何通过电路测量反映模型中的物理性质. 在非厄米系统中, 大家主要关心的是复能谱和本征态以及由它们所构造的各种物理量等静态性质, 或者是量子态演化等动力学性质. 因此, 本节将概述经典电路实验中常用的测量方法与物理量之间的对应.
3.3.1 导纳矩阵
根据电路的拉普拉斯形式, 测量出导纳谱
$ \{j_n\} $ 便可得到相应格点模型能谱$ \{E_n\} $ 的信息. 但导纳谱通常无法被直接测量, 因为根据导纳矩阵的定义(17)式, 节点电压$ {\cal{V}} $ 对节点输入电流$ {\cal{I}} $ 的响应直接反映的是电路格林函数$ {\cal{G}} $ . 所以, 可以先测出电路格林函数, 再根据$ {\cal{J}} = {\cal{G}}^{-1} $ 得到导纳矩阵, 进而将其对角化得到导纳谱.下面介绍如何测量电路格林函数. 假设电路有N个节点, 如果仅在第j节点注入频率为ω的交变电流
$ I_j{\rm e}^{\mathrm{{\rm i}}\omega t} $ , 由$ {\cal{V}} = {\cal{G}}{\cal{I}} $ 可得稳态响应电压的幅矢量为$$ {\cal{V}} = {\cal{G}} \begin{pmatrix} 0 \\ \vdots \\ I_j\\ \vdots \\ 0 \end{pmatrix} = \begin{pmatrix} {\cal{G}}_{1j} \\ \vdots \\ {\cal{G}} _{Nj} \end{pmatrix} I_j. $$ (32) 于是通过测量每个节点的稳态电压幅
$ V_i $ , 便可得到电路格林函数第j列的矩阵元$ {\cal{G}} _{ij} = V_i/I_j $ . 当以相同的频率对电路中所有节点遍历上述电流驱动, 利用频谱分析仪测量每次驱动下所有节点的稳态电压幅, 即可完全重构电路格林函数, 进而得到导纳矩阵[38]. 这个过程所需要的总的测量次数为$ N^2 $ . 而对于具有平移不变性的周期性电路而言, 仅需遍历一个电路原胞即可, 测量次数可以减少到$ N^2/N_{\rm{c}} $ , 其中$ N_{\rm{c}} $ 表示原胞数[41]. 将得到的导纳矩阵$ {\cal{J}} $ 对角化, 便可得到本征值$ j_n $ 和左右本征矢$ v_n^{({\rm{l, r}})} $ 以及由它们所构成的其他静态物理量; 本征频率$ \tilde{\omega} $ 也可以通过$ {\rm{det}}\left[{\cal{J}}(\widetilde{\omega})\right] = 0 $ 得到.在实验中更高效的方法是, 利用矢量网络分析仪(vector network analyzer, VNA)对电路各节点的散射信号进行扫描, 再将包含散射信息的S参数矩阵(S-parameter matrix)变换到导纳矩阵, 便可计算出电路的导纳谱和本征频率谱[114,115,130,132,134,139-142].
另外, 如果仅是定性地表征导纳矩阵本征矢的特性, 也可以直接测量各个节点对交变电流的电压响应[41]. 根据下式
$$ {\cal{V}} = {\cal{G}}{\cal{I}} = \sum\limits_n j_n^{-1} v_n^{(\mathrm{r})}\left[ v_{n}^{(\mathrm{l})\ast} {\cal{I}} \right]$$ (33) 可知, 电压响应
$ {\cal{V}} $ 是所有导纳矩阵右本征矢的某种线性叠加. 据此可以从电压在节点间的分布粗略判断系统是否发生类似非厄米趋肤效应或安德森局域化的现象[117].3.3.2 两点间阻抗
电路实验中也经常通过测量两点间阻抗(two-point impedance)来表征所关心的模拟量的性质. 节点i到j的阻抗
$ Z_{ij} $ 可以通过输入交变电流$ I_i{\rm e}^{{\rm i}\omega t} $ 并测量两节点的响应电压$ V_{i}{\rm e}^{{\rm i}\omega t} $ 和$ V_{j}{\rm e}^{{\rm i}\omega t} $ 得到. 利用(20)式中导纳矩阵$ {\cal{J}} $ 的本征值$ j_n $ 和本征矢$ v_n^{({\rm{l, r}})} $ , 可以将两点间阻抗表示为[143]$$ \begin{split} Z_{ij} \equiv\;& \frac{V_i-V_j}{I_i} = \sum\limits_k \frac{{\cal{G}}_{ik}I_k-{\cal{G}}_{jk}I_k}{I_i} \\ =\;& \sum\limits_n\frac{\big[v_{ni}^{(\mathrm{l})\ast}+ v_{nj}^{(\mathrm{l})\ast}I_j/I_i\big]\big[ v_{ni}^{(\mathrm{r})} -v_{nj}^{(\mathrm{r})}\big]}{ j_n}. \end{split} $$ (34) 其中,
$ v_{ni}^{\rm{(l, r)}} $ 代表$ v_n^{({\rm{l, r}})} $ 的第i个矩阵元. (34)式说明, 当驱动频率ω接近电路的某个本征频率$ \widetilde{\omega} $ 时, 阻抗$ Z_{ij} $ 将由值很小的$ j_n(\omega\sim \widetilde{\omega}) $ 项主导从而趋于发散, 即$ Z_{ij}(\omega) $ 在$ \widetilde{\omega} $ 处出现共振峰. 更简单的方案是测量对地阻抗, 即令(34)式中的$ v_{n{\rm{g}}}^{({\rm{l, r)}}} = 0 $ , 可得$$ Z_{i{\rm{g}}} \equiv \dfrac{V_i}{I_i} = \sum\limits_n j_n^{-1} v_{ni}^{(\mathrm{l})\ast} v_{ni}^{(\mathrm{r})}. $$ (35) 这样可以通过对每个节点扫频测得两点间阻抗的共振峰的位置, 从而确定电路的本征频率谱. 固定本征频率
$ \widetilde{\omega} $ , 测量对地阻抗$ Z_{i{\rm{g}}}(\widetilde{\omega}) $ 随节点i的变化也反映了导纳本征矢$ v_{ni}^{({\rm{l, r}})}(\widetilde{\omega}) $ (即电路的本征电压)在节点间的分布特点, 常被应用于拓扑边缘态和非厄米趋肤效应[35,74,119,128,144,145]以及高阶拓扑角态[126,127,146]的测量. 不过需要注意的是, 本方法只适用于本征频率谱为近实谱(即本征频率的虚部很小)的情况, 因为扫描频率ω只能是实数.3.3.3 动力学
前两种测量主要模拟的是格点模型的静态性质, 对于动力学性质的模拟, 可直接利用电路的刘维尔方程(23)[131,137,147]或耦合模方程(28)[114,115,130], 类比格点模型的含时薛定谔方程. 此类测量非常简单, 只要按模拟需求给电路各节点注入相应的初始电压或电流, 随后测量其随时间的变化即可.
利用动力学也能验证格点模型的部分静态性质. 比如利用导纳矩阵根据
$ {\rm{det}}\left[{\cal{J}}(\widetilde{\omega})\right] = 0 $ 已得到了系统的本征频率$ \widetilde{\omega} $ (一般为复数), 通过动力学的办法可以观测到此本征频率且判断其正确性. 方法是, 将与本征频率$ \widetilde{\omega} $ 相应的本征电压$ v_0 $ 设为电路的初态, 随后撤掉外部源, 则初态将会按照刘维尔方程(23)进行演化, 其任意时刻的电压应满足:$$ {\cal{V}}(t) = v_0 {\rm e}^{\mathrm{{\rm i}} \widetilde{\omega} t} = v_0 {\rm e}^{-{\rm{Im}}(\widetilde{\omega})t}{\rm e}^{\mathrm{{\rm i}}{\rm{Re}}(\widetilde{\omega})t}. $$ (36) 因此, 可根据电压随时间的振荡频率和幅值的放大/衰减程度分别拟合出本征频率
$ \widetilde{\omega} $ 的实部和虚部, 与理论相比较[122,135,148].虽然电压的动力学行为并不是遵从于导纳矩阵
$ {\cal{J}}(\widetilde{\omega}) $ 的本征值$ j_n(\widetilde{\omega}) $ , 而是依赖于本征频率$ \widetilde{\omega} $ , 但是由于本征频率谱与导纳谱结构上的相似性, 可以用本征电压的动力学定性的反映导纳矩阵本征矢的动力学, 进而模拟相关的行为, 例如拓扑边缘态和非厄米趋肤效应的动力学行为[137,147]以及场相干的演化[127]等.4. 经典电路模拟非厄米格点模型的实验进展
目前, 人们利用经典电路在实验上已经成功模拟和观测了非厄米格点模型的许多性质. 本节将从PT对称破缺、多种非厄米趋肤效应、其他非厄米拓扑态以及EP在传感器上的应用这4个方面介绍相关的实验进展.
4.1 PT对称破缺
PT对称破缺是非厄米系统所独有的现象. 在经典电路实验中, 人们将增益电路与损耗电路耦合形成二聚体(dimer)电路[113,133,135,148]或者在二聚体电路中进一步添加一个中性电路构成的三聚体(trimer)电路[130,141], 观测电路的本征频率谱随增益/损耗强度的变化, 从而对PT对称破缺进行研究. 本节将从不同的电路模型出发对相关实验进展进行介绍.
最早在电路平台上研究PT对称破缺的是Schindler等[113]. 他们在2011年设计了一个由一对RLC电路耦合形成的二聚体电路 [图6(a)][133], 其中增益和损耗部分分别由接地运放和电阻实现. 他们根据基尔霍夫定律写出描述电量和电流变化的刘维尔方程(23)式, 建立起与含时薛定谔方程的对应关系, 在实验上观测本征频率随增益/损耗强度的变化. 如图6(b)所示, 在
$ \gamma/\gamma_{\mathrm{PT}} > 1 $ 时, 观测到复本征频率出现, 对应于PT对称破缺. 2017年, 该研究组[148]基于Floquet理论对一系列由EP分隔开的PT对称相与PT对称破缺相进行了研究. 他们在增益电路中用光电池(photocell)连接金属-氧化物场效应晶体管(MOSFET)的源极与电感, 并将MOSFET的漏极接入振荡电路. 当光电池两端的电压降低时, 其电流会流入电路, 从而实现增益. 在损耗电路中使用接地的光电池连接电感, 其电流由于接地会耗散, 从而实现损耗. 利用变容二极管(varactor diode)耦合增益电路与损耗电路, 实现对电路耦合强度的控制.2020年, Wang等[135]利用接地的INIC实现等效负阻抗, 搭建了由两个RLC回路组成的对称二聚体电路, 在频谱中观测到一阶和二阶EP以及PT对称相和PT对称破缺相. 次年, Zhou等[115]同样搭建了由两个RLC回路耦合的二聚体电路, 基于耦合模理论(25)式对格点模型进行模拟, 研究了微扰对PT对称系统的影响. 实验将VNA的输入端连接到电路实现等效的负阻抗, 同时将另一端连接到到读取器(reader)中测量频率. VNA自身带有阻抗, 连接输入源端口到闭合回路中可以输入电能, 在电路中可以起到增益的作用.
随着人们对于PT对称认识的加深, 对其模拟的电路平台不仅局限于二聚体电路. 三聚体电路由于中性电路的加入具有更强的可扩展性, 同样是重要的模拟平台. 2019年, Sakhdari等[141]在原来的二聚体电路之间添加了一个LC振荡电路, 并利用VNA实现等效负阻抗并同时进行频率测量. 他们比较了二聚体与三聚体电路的实验结果, 观测到三聚体频谱中接近发散EP (divergent exceptional points, DEP)的频率偏移程度较二聚体更大. 2023年, Yin等[130]采用了与Zhou等[115]类似的方法, 利用VNA同时实现负阻抗与测量, 构建了一个由3个RLC回路耦合的三聚体电路, 并观测到环绕周期为3的三阶EP.
除了二聚体或三聚体电路之外, 2021年, Stegmaier等[120]构建了一个类SSH模型的电路来研究PT对称性. 他们根据电路的拉普拉斯形式(17)式, 利用导纳矩阵建立起电路与格点模型之间的对应关系. 他们在电路的每个节点上添加接地电阻以实现格点上的增益和损耗, 并测量了不同增益和损耗强度下的导纳谱. 通过这些测量, 他们观察到了PT对称、PT对称破缺和反PT对称这3种情况. 该实验为研究PT对称破缺提供了一个新的平台, 并为进一步将EP应用在传感器上提供了新的思路.
4.2 多种非厄米趋肤效应
非互易跃迁往往可以诱导出非厄米趋肤效应, 从而导致传统体边对应关系的失效. 在经典电路的实验中, 研究者们通常在电路节点间加入INIC[41,74,121,123,127,138,145,146]或电压跟随器[119,124,129]等主动元件实现对非互易跃迁的模拟, 通过电路中电压响应的分布观测非厄米趋肤效应. 目前, 利用经典电路已实现对非厄米系统中多种非厄米趋肤效应及相关奇异现象的研究和模拟, 包括一维非厄米趋肤效应和传统体边对应关系的失效[41,119]、非布洛赫波的演化和广义布里渊区[127]、高维非厄米趋肤效应[118]、高阶非厄米趋肤效应[121,124,128]、多体非厄米趋肤效应[123]、临界非厄米趋肤效应[125]以及尺寸依赖的非厄米趋肤效应[129]等.
2020年, Helbig等[41]首次在电路实验中观测到传统体边对应关系的失效与非厄米趋肤效应的出现. 他们研究了一维非厄米SSH模型, 在电路节点间加入INIC实现对等效的非互易跃迁的模拟, 并通过电路的拉普拉斯形式(17)式建立起导纳矩阵与SSH模型的对应关系. 他们利用测量节点电压响应的方法(33)式, 对电路中所有节点施加了电流驱动, 测量各节点上的电压响应, 观测到节点上的电压响应局域在电路的一侧边界, 即在该电路中实现了非厄米趋肤效应. 次年, Liu等[119]利用电压跟随器构建了具有非互易跃迁的SSH电路. 他们观察到由传统体边对应关系失效引起的电路本征频率谱的变化. 与此同时, 他们基于广义布里渊区的理论对能量缠绕数(13)式进行了计算, 验证了这种新的体边对应关系.
非厄米趋肤效应本质是一种非布洛赫效应. 在经典电路平台上研究者们已经实现了非布洛赫波的演化及对广义布里渊区的验证. 2023年, Wu等[127]构建了如图7(a)所示的二维电路来观测非布洛赫动力学. 实验中, 他们利用INIC实现了格点间等效的非互易跃迁, 利用傅里叶变换场扫描(Fourier-transformed field scan)的方法对节点的电压分布进行测量, 将其进行拉普拉斯变换(Laplace transform)后得到在广义布里渊区的分布 [图7(b)], 结果与理论一致. 由此验证了广义布里渊区在描述非布洛赫动力学上的有效性.
人们在经典电路实验里还研究了更多影响非厄米趋肤效应的因素. 2022年, Zhang等[125]研究了一个由两条Hatano-Nelson链耦合成的类梯子模型. 他们在实验中测量了在开边界和周期边界条件下各节点的阻抗, 发现当节点数增多到一定程度时, 周期边界的阻抗与开边界的阻抗比值明显大于1, 进一步表明边界连接对远距离的阻抗响应具有非局域的影响, 即具有临界非厄米趋肤效应.
2023年, Su等[129]在电路平台上对单向跃迁的格点模型进行了研究. 他们利用电压跟随器在电路中模拟单向跃迁, 并在不同的节点数量下测量了导纳谱研究发现随着节点数目的增多, 集中在一侧的局域化程度会降低, 实现了对尺寸依赖的非厄米趋肤效应的观测. 同年, 他们也对非厄米趋肤效应与安德森局域化之间的竞争进行了研究[149]. 实验仅采用了线性元件, 通过赋予元件幂指数变化的参数值以模拟格点间的非互易跃迁, 成功实现了非厄米趋肤效应. 在引入准无序后成功观测到非厄米趋肤效应与安德森局域化的竞争, 发现当无序强度足够大时, 本征态从非厄米趋肤相转变为局域相.类似地, Tang等[146]基于二维非互易蜂窝(honeycomb)电路对非厄米趋肤效应与角态之间的竞争关系进行了研究. 实验中, 他们通过调整INIC里电容的大小以改变非互易强度, 并对电路节点的电压响应进行测量, 发现随着非互易强度的增强, 节点电压会从电路两端逐渐转移到同一端, 呈现出角态的特征.
随着对非厄米趋肤效应研究的不断深入以及经典电路模拟方法的不断发展, 研究者们开始将对非厄米趋肤效应的经典电路模拟推向高维、高阶以及多体等更为复杂的系统. 2021年, Zou等[121]基于电路的拉普拉斯形式(17)式, 分别建立了电路的导纳矩阵与二维格点模型以及三维格点模型之间的对应关系. 对杂化高阶拓扑趋肤效应(hybrid higher-order topological skin effect)进行了研究. 杂化高阶拓扑趋肤效应是在高维系统中, 由非厄米趋肤效应以及拓扑局域态相互影响所形成的一种拓扑现象[150]. 实验中, 他们在电路耦合处引入INIC实现等效的非互易跃迁, 继而测量各节点上的电压响应分布, 通过调整INIC的连接方向实现不同的杂化方式, 观测到二阶杂化趋肤-拓扑效应(即电压集中分布在二维平面中的一组对角节点)、二阶杂化趋肤-趋肤效应(即电压分布在二维平面中的其中一个角节点上)、三阶杂化趋肤-趋肤-趋肤效应(即电压集中分布在三维立体中的角节点上)以及三阶杂化趋肤-拓扑-拓扑效应(即电压分布在一对平面上的对角节点上). 2022年, Deng等[145]研究了一个一维N次方根非厄米SSH模型 [图7(c)], 其中N的大小可以通过调整原胞内电容电感等元器件的数目实现. 实验中, 他们选择了不同的N对节点电压进行测量, 观测到电压分布集中在电路节点的一侧, 即出现了非厄米趋肤效应, 并发现非厄米趋肤效应随N的增大变得更加明显. 同年, 他们还对一个一维三粒子Bose-Hubbard模型进行了研究[123]. 他们将此模型的本征态在希尔伯特空间中重新表示为Fock态的形式后, 写出系统的薛定谔方程, 并根据基尔霍夫定律直接构造电路模拟该薛定谔方程. 该团队测量了节点上的阻抗分布, 发现阻抗值集中在一个节点上, 观察到了多体系统中由阻抗表征的非厄米趋肤效应, 这被称为非厄米聚集效应(non-Hermitian aggregation effect). 2023年, 他们在三维电路系统中基于INIC设计了一个拓扑开关, 也可以用于研究非厄米趋肤效应[138].
2022年, Shang等[124]研究了一个二维的非互易格点模型, 他们在实验中计算了二维的非布洛赫能量缠绕数, 发现由其表征的二阶非厄米趋肤效应. 2023年, Zhu等[126]构建了一个具有二阶手性的二维电路. 他们利用标准电阻引入非厄米项, 重构电路格林函数以测量每个节点的电压响应, 观测到阻抗的响应集中在边界局域态上, 展示了二维二阶非厄米趋肤效应.
基于经典电路平台同样可以实现互易非厄米趋肤效应. 2020年, Hofmann等[118]通过在二维电路中的对角节点连接电阻引入非厄米, 研究了一个具有π-通量的格点模型. 实验中, 他们测量并计算了导纳矩阵本征矢 [由(20)式定义] 的倒参与率(inverse participation ratio), 对动量空间中
$ k_y $ 方向上的局域态进行了研究, 观察到在$ k_y = \pi/2, \; k_y = 3\pi/2 $ 处的局域化程度最高, 与理论预测的非厄米π-通量模型呈现出互易趋肤效应时的特征相同.4.3 其他非厄米拓扑态
在用经典电路模拟非厄米趋肤效应的同时, 人们也对其他非厄米拓扑物态进行了模拟, 包括单纯由增益/耗损诱导的非厄米拓扑态[74]、非厄米连续谱中的束缚态(bound state in the continuum)[144]、杂化高阶趋肤拓扑态[121]、高阶拓扑角态[122,127]、N次方根拓扑相[145]、二阶手性系统中的拓扑态及动力学[126]、非厄米拓扑开关[138]以及非厄米Hopf束[128]等.
2020年, Liu等[74]首次实现了单纯由增益/损耗诱导的一维非厄米拓扑模型 [图8(a)]. 他们在电路 [图8(b)] 中引入INIC实现等效负电阻, 通过调控接地电阻
$ R_0 $ 以及代表增益(损耗)的正(负)电阻$ R_1 (R_2) $ , 可观测到不同的拓扑相. 同年, Li 等[144]利用电路研究了一个由两条SSH链耦合构成的梯子模型, 在两条链之间引入耦合电阻使其变成非厄米模型. 实验中, 他们测量了特定节点上的阻抗随输入频率的变化, 发现阻抗的峰值处于特定频率, 对应于非厄米连续谱中的束缚态.2021年, Zou等[121]设计了一个二维格点模型, 对杂化二阶趋肤-拓扑态进行了研究. 实验中, 他们设计了3种不同的INIC, 并通过改变节点间INIC的连接方向实现对趋肤-趋肤和趋肤-拓扑两种态的模拟. 通过测量不同输入频率下电路节点的电压响应, 观测到在杂化二阶趋肤-拓扑或趋肤-趋肤情况下的角态、边缘态以及体态. 杂化趋肤-拓扑态可被用于拓扑开关的设计上[151], 可以通过控制电路内的能量传递实现对趋肤效应的“开关”[121]. 2023年, Zhang等[138]在电路实验中实现拓扑开关. 他们在节点间连接INIC实现对非互易跃迁的模拟, 观测到在三维系统中趋肤-拓扑-趋肤效应以及趋肤-拓扑-拓扑效应之间的转变, 表明利用INIC设计拓扑开关的可行性.
2022年, Wu等[43]在研究非厄米二维电路的基础上, 设计了一个模拟二阶拓扑绝缘体的二维拓扑电路 [图9(a)][122]. 实验中, 他们通过电阻引入非厄米项, 并调整两种电阻的大小关系以诱导不同的拓扑态. 通过对电路本征频率以及节点电压响应的测量, 他们观察到带隙中的一维边缘态 [图9(b)] 以及零维角态 [图9(c)]. 此外, 他们还测量了在电路的角态、边缘态和体态的阻抗响应随驱动频率的变化 [图9(d)].
2023年, Zhu等[126]构建了一个具有二阶手性的二维电路. 他们利用标准电阻引入非厄米项, 观测到阻抗响应局域在电路节点平面的角、边和体三种情况, 分别表征二阶手性二维格点模型中的角态、边缘态和体态.
4.4 EP在传感器上的应用
EP在传感器上的应用与PT对称破缺的实验研究密切相关. 实现EP传感器的电路一般也采取PT对称二聚体电路[114,115,132,140,142]或三聚体电路[116,152]的结构. 实验中通常采用测得的反射谱(reflection spectra)来对传感器的灵敏度进行表征[116,132,140,142,152], 反射谱的响应越大, 即反射谱曲线中波谷所处位置越深 [图10(b)], 传感器越灵敏.
2018年, Chen等[140]成功将EP应用于无线压力传感器(wireless pressure sensor)中, 该传感器通过电容的阻值对压力进行表征. 在增益电路中, 他们改用了可以调参的电感与损耗电路进行耦合, 由此实现了更一般情况下的PTX对称电路, 这里的X表示一种互易的标度操作(reciprocal-scaling operation). 实验中, 他们选择EP附近和非EP处的多个不同耦合强度来对反射谱进行测量与对比, 观测到在相同频段内, 选在EP附近的耦合强度, 反射会更强, 由此表明EP可以增强传感器的灵敏度. 在实际压力传感测试中, 实验人员利用VNA分别对PTX电路、PT电路以及未应用EP的普通传感器的反射谱进行了测量, 在相同的频段内, PTX电路和PT电路都可以增强反射. 由于PTX电路相比于PT电路具有不同的电路特征模式, 可以放大特定频率的反射系数, 具有更大的操作空间.
2019年, Dong等[114]在动物表皮层植入用于测量生物体征(压力, 温度等)的传感器, 并利用读取器对感应频率的偏移进行了测量, 通过频率偏移变化曲线对测量的物理量进行表征. 实验观测到, 随耦合强度变化的频率偏移在EP附近明显增高. 该团队将基于EP设计的传感器的灵敏度提升至当时极限的3.2倍. 2020年, Zhou等[132]比较了传统的被动无线传感器与基于PT对称性的被动无线传感器在灵敏度上的差别 [图10(a)], 根据实验中VNA测得的反射谱 [图10(b)], 在增大传感器距离的同时, 响应仍然能够保持在较强的范围内, 即EP可提升传感器的控制距离. 2021年, Zeng等[116]将三阶EP应用到传感器中, 设计了一个三线圈耦合的电路. 与二阶EP传感器相比, 实验测得的反射谱反映出高阶EP对灵敏度的增强更为显著. 2022年, Yin等[142]利用运放在增益侧实现了等效的负阻抗, 进而设计了一个不需要引入VNA的无线传感器. 得益于运放的非线性饱和效应, 该传感器可以通过观测读取器对电容的阻值进行实时读出. 同年, 他们在二聚体电路的基础上加入了一个由电感与可变电容构成的振荡电路, 设计了一个三聚体电路 [图10(c)][152]. 相比于二聚体电路, 三聚体电路测得的反射谱具有更强的响应. 选取多组阻值的电容进行实验, 三聚体电路测得的反射谱稳定在较强的响应区间, 比二聚体电路具有更高的品质因子(Q factor).
SSH电路同样可应用在EP传感器上. 2023年, Yuan等[139]在实验上采用具有非互易跃迁的SSH模型作为平台设计传感器, 原胞内通过电压跟随器实现格点间的非互易跃迁. 该传感器根据电容的阻值以及测得的频率偏移可以对位置、角度以及水位变化高度等物理量进行表征.
非厄米系统中的增益和损耗不可避免地会引入噪声[153]. 最近, 研究者们注意到在利用EP传感器进行信号放大的过程中往往伴随着噪声的增强, 因此亟待解决的问题是如何提高EP传感器的信噪比. 噪声增强的来源可分为两种: 一是物理上的, 来源于产生EP过程中本征矢的重合; 另一种是技术上的, 来源于EP传感器实现中的放大机制[154]. 物理学家尝试从多个角度解决噪声对EP传感器的影响. 理论上, Langbein[155]在2018年讨论了EP传感器的量子有限信噪比与扰动, Zhang等[153]在2019年提出计算EP传感器信噪比的量子噪声理论; 实验上, Kononchuk等[154]在2022年设计了PT对称式的机电EP传感器等.
5. 总 结
本文通过简单介绍非厄米物理的基本理论和新奇现象, 结合经典电路模拟的理论基础, 概述了当下经典电路模拟非厄米格点模型的实验进展. 目前的实验进展主要涵盖对PT对称破缺、非厄米趋肤效应以及各类新奇非厄米拓扑物态的经典电路模拟. 实验的模拟对象也从简单的一维单体格点模型逐渐推广至高维多体等更复杂的格点模型. 此外, 基于经典电路对非厄米新奇效应的应用也是重要的研究方向, 比如高灵敏度的EP传感器等. 相比于量子模拟平台的高技术门槛, 经典电路的低成本、易扩展等特点使其成为模拟和研究丰富的非厄米新奇现象的首选平台, 表现出广阔的研究和应用前景.
[1] Feynman R P 1982 Int. J. Theor. Phys. 21 467Google Scholar
[2] Greiner M, Mandel O, Esslinger T, Hänsch T W, Bloch I 2002 Nature 415 39Google Scholar
[3] Billy J, Josse V, Zuo Z, Bernard A, Hambrecht B, Lugan P, Clement D, Sanchez-Palencia L, Bouyer P, Aspect A 2008 Nature 453 891Google Scholar
[4] Zhang D W, Zhu Y Q, Zhao Y X, Yan H, Zhu S L 2018 Adv. Phys. 67 253Google Scholar
[5] Friedenauer A, Schmitz H, Glueckert J T, Porras D, Schaetz T 2008 Nat. Phys. 4 757Google Scholar
[6] Barreiro J T, Müller M, Schindler P, Nigg D, Monz T, Chwalla M, Hennrich M, Roos C F, Zoller P, Blatt R 2011 Nature 470 486Google Scholar
[7] Blatt R, Roos C F 2012 Nat. Phys. 8 277Google Scholar
[8] You J Q, Shi X F, Hu X, Nori F 2010 Phys. Rev. B 81 014505Google Scholar
[9] Georgescu I, Ashhab S, Nori F 2014 Rev. Mod. Phys. 86 153Google Scholar
[10] Barends R, Kelly J, Megrant A, Veitia A, Sank D, Jeffrey E, White T C, Mutus J, Fowler A G, Campbell B, Chen Y, Chen Z, Chiaro B, Dunsworth A, Neill C, O’Malley P, Roushan P, Vainsencher A, Wenner J, Korotkov A N, Cleland A N, Martinis J M 2014 Nature 508 500Google Scholar
[11] Buluta I, Nori F 2009 Science 326 108Google Scholar
[12] Ashida Y, Gong Z P, Ueda M 2020 Adv. Phys. 69 249Google Scholar
[13] Raghu S, Haldane F D M 2008 Phys. Rev. A 78 033834Google Scholar
[14] Wang Z, Chong Y, Joannopoulos J D, Soljacic M 2009 Nature 461 772Google Scholar
[15] Rechtsman M C, Zeuner J M, Plotnik Y, Lumer Y, Podolsky D, Dreisow F, Nolte S, Segev M, Szameit A 2013 Nature 496 196Google Scholar
[16] Hafezi M, Mittal S, Fan J, Migdall A, Taylor J M 2013 Nat. Photonics 7 1001Google Scholar
[17] Yan Q, Liu R, Yan Z, Liu B, Chen H, Wang Z, Lu L 2018 Nat. Phys. 14 461Google Scholar
[18] Ozawa T, Price H M, Amo A, Goldman N, Hafezi M, Lu L, Rechtsman M C, Schuster D, Simon J, Zilberberg O, Carusotto I 2019 Rev. Mod. Phys. 91 015006Google Scholar
[19] Hussein M I, Leamy M J, Ruzzene M 2014 Appl. Mech. Rev. 66 040802Google Scholar
[20] Wang P, Lu L, Bertoldi K 2015 Phys. Rev. Lett. 115 104302Google Scholar
[21] Xiao Y X, Ma G, Zhang Z Q, Chan C T 2017 Phys. Rev. Lett. 118 166803Google Scholar
[22] Li F, Huang X, Lu J, Ma J, Liu Z 2018 Nat. Phys. 14 30Google Scholar
[23] Li X X, Li B, Li P B 2020 Phys. Rev. Res. 2 013121Google Scholar
[24] Li X X, Li P B, Li H R, Gao H, Li F L 2021 Phys. Rev. Res. 3 013025Google Scholar
[25] Kane C L, Lubensky T C 2014 Nat. Phys. 10 39Google Scholar
[26] Chen B G g, Upadhyaya N, Vitelli V 2014 Proc. Natl. Acad. Sci. 111 13004Google Scholar
[27] Nash L M, Kleckner D, Read A, Vitelli V, Turner A M, Irvine W T M 2015 Proc. Natl. Acad. Sci. 112 14495Google Scholar
[28] Chaunsali R, Kim E, Thakkar A, Kevrekidis P G, Yang J 2017 Phys. Rev. Lett. 119 024301Google Scholar
[29] Brandenbourger M, Locsin X, Lerner E, Coulais C 2019 Nat. Commun. 10 4608Google Scholar
[30] Ghatak A, Brandenbourger M, van Wezel J, Coulais C 2020 Proc. Natl. Acad. Sci. 117 29561Google Scholar
[31] 罗开发, 余睿 2019 物理学报 68 220305Google Scholar
Luo K F, Yu R 2019 Acta Phys. Sin. 68 220305Google Scholar
[32] 刘硕, 张霜, 崔铁军 2021 中国光学 14 736Google Scholar
Liu S, Zhang S, Cui T J 2021 Chinese Optics 14 736Google Scholar
[33] Albert V V, Glazman L I, Jiang L 2015 Phys. Rev. Lett. 114 173902Google Scholar
[34] Ningyuan J, Owens C, Sommer A, Schuster D, Simon J 2015 Phys. Rev. X 5 021031Google Scholar
[35] Lee C H, Imhof S, Berger C, Bayer F, Brehm J, Molenkamp L W, Kiessling T, Thomale R 2018 Commun. Phys. 1 39Google Scholar
[36] Imhof S, Berger C, Bayer F, Brehm J, Molenkamp L W, Kiessling T, Schindler F, Lee C H, Greiter M, Neupert T, Thomale R 2018 Nat. Phys. 14 925Google Scholar
[37] Ezawa M 2018 Phys. Rev. B 98 201402Google Scholar
[38] Helbig T, Hofmann T, Lee C H, Thomale R, Imhof S, Molenkamp L W, Kiessling T 2019 Phys. Rev. B 99 161114Google Scholar
[39] Liu S, Gao W, Zhang Q, Ma S, Zhang L, Liu C, Xiang Y J, Cui T J, Zhang S 2019 Research 2019 8609875Google Scholar
[40] Bao J, Zou D, Zhang W, He W, Sun H, Zhang X 2019 Phys. Rev. B 100 201406Google Scholar
[41] Helbig T, Hofmann T, Imhof S, Abdelghany M, Kiessling T, Molenkamp L W, Lee C H, Szameit A, Greiter M, Thomale R 2020 Nat. Phys. 16 747Google Scholar
[42] Lee C H, Sutrisno A, Hofmann T, Helbig T, Liu Y, Ang Y S, Ang L K, Zhang X, Greiter M, Thomale R 2020 Nat. Commun. 11 4385Google Scholar
[43] Wu J, Huang X, Lu J, Wu Y, Deng W, Li F, Liu Z 2020 Phys. Rev. B 102 104109Google Scholar
[44] Kotwal T, Moseley F, Stegmaier A, Imhof S, Brand H, Kießling T, Thomale R, Ronellenfitsch H, Dunkel J 2021 Proc. Natl. Acad. Sci. 118 e2106411118Google Scholar
[45] Dong J, Juričić V, Roy B 2021 Phys. Rev. Res. 3 023056Google Scholar
[46] Yang Y, Zhu D, Hang Z, Chong Y 2021 Sci. China: Phys., Mech. Astron. 64 257011Google Scholar
[47] Zhang W, Zou D, Pei Q, He W, Bao J, Sun H, Zhang X 2021 Phys. Rev. Lett. 126 146802Google Scholar
[48] Wu J, Wang Z, Biao Y, Fei F, Zhang S, Yin Z, Hu Y, Song Z, Wu T, Song F, Yu R 2022 Nat. Electron. 5 635Google Scholar
[49] Wang Z, Zeng X T, Biao Y, Yan Z, Yu R 2023 Phys. Rev. Lett. 130 057201Google Scholar
[50] Breuer H P, Petruccione F 2002 The Theory of Open Quantum Systems (Oxford: Oxford University Press
[51] Daley A J 2014 Adv. Phys. 63 77Google Scholar
[52] Bender C M, Boettcher S 1998 Phys. Rev. Lett. 80 5243Google Scholar
[53] Bender C M 2007 Rep. Prog. Phys. 70 947Google Scholar
[54] Longhi S 2019 Opt. Lett. 44 5804Google Scholar
[55] Pan J S, Yi W, Gong J 2021 Commun. Phys. 4 261Google Scholar
[56] Moca C P, Sticlet D, Dóra B, Zaránd G 2023 Phys. Rev. B 107 115111Google Scholar
[57] Lee T E 2016 Phys. Rev. Lett. 116 133903Google Scholar
[58] Kunst F K, Edvardsson E, Budich J C, Bergholtz E J 2018 Phys. Rev. Lett. 121 026808Google Scholar
[59] Yao S, Wang Z 2018 Phys. Rev. Lett. 121 086803Google Scholar
[60] Yao S, Song F, Wang Z 2018 Phys. Rev. Lett. 121 136802Google Scholar
[61] Edvardsson E, Kunst F K, Bergholtz E J 2019 Phys. Rev. B 99 081302Google Scholar
[62] Ghatak A, Das T 2019 J. Phys.: Condens. Matter 31 263001Google Scholar
[63] Song F, Yao S, Wang Z 2019 Phys. Rev. Lett. 123 170401Google Scholar
[64] Song F, Yao S, Wang Z 2019 Phys. Rev. Lett. 123 246801Google Scholar
[65] Borgnia D S, Kruchkov A J, Slager R J 2020 Phys. Rev. Lett. 124 056802Google Scholar
[66] Xiao L, Deng T, Wang K, Zhu G, Wang Z, Yi W, Xue P 2020 Nat. Phys. 16 761Google Scholar
[67] Zhang K, Yang Z, Fang C 2020 Phys. Rev. Lett. 125 126402Google Scholar
[68] Yang Z, Zhang K, Fang C, Hu J 2020 Phys. Rev. Lett. 125 226402Google Scholar
[69] Okuma N, Sato M 2021 Phys. Rev. Lett. 126 176601Google Scholar
[70] Liang Q, Xie D, Dong Z, Li H, Li H, Gadway B, Yi W, Yan B 2022 Phys. Rev. Lett. 129 070401Google Scholar
[71] Milburn T J, Doppler J, Holmes C A, Portolan S, Rotter S, Rabl P 2015 Phys. Rev. A 92 052124Google Scholar
[72] Zhong Q, Christodoulides D N, Khajavikhan M, Makris K G, El-Ganainy R 2018 Phys. Rev. A 97 020105Google Scholar
[73] Zhang K L, Wang P, Song Z 2019 Phys. Rev. A 99 042111Google Scholar
[74] Liu S, Ma S, Yang C, Zhang L, Gao W, Xiang Y J, Cui T J, Zhang S 2020 Phys. Rev. Appl. 13 014047Google Scholar
[75] Heiss W D 2004 J. Phys. A: Math. Gen. 37 2455Google Scholar
[76] Moiseyev N 2011 Non-Hermitian Quantum Mechanics (Cambridge: Cambridge University Press
[77] Okugawa R, Takahashi R, Yokomizo K 2020 Phys. Rev. B 102 241202Google Scholar
[78] Zhang X, Tian Y, Jiang J H, Lu M H, Chen Y F 2021 Nat. Commun. 12 5377Google Scholar
[79] Zhang X, Zhang T, Lu M H, Chen Y F 2022 Adv. Phys.: X 7 2109431Google Scholar
[80] Deng T S, Yi W 2019 Phys. Rev. B 100 035102Google Scholar
[81] Yokomizo K, Murakami S 2019 Phys. Rev. Lett. 123 066404Google Scholar
[82] Zhu X, Wang H, Gupta S K, Zhang H, Xie B, Lu M, Chen Y 2020 Phys. Rev. Res. 2 013280Google Scholar
[83] Yokomizo K, Murakami S 2020 Prog. Theor. Exp. Phys. 2020 12A102Google Scholar
[84] Xiao L, Deng T, Wang K, Wang Z, Yi W, Xue P 2021 Phys. Rev. Lett. 126 230402Google Scholar
[85] Wang K, Li T, Xiao L, Han Y, Yi W, Xue P 2021 Phys. Rev. Lett. 127 270602Google Scholar
[86] Yokomizo K 2022 Non-Bloch Band Theory of Non Hermitian Systems (Singapore: Springer Nature
[87] Yokomizo K, Murakami S 2023 Phys. Rev. B 107 195112Google Scholar
[88] Shi K, Tian M, Sun F X, Zhang W 2023 Phys. Rev. B 107 205154Google Scholar
[89] Zhou H, Lee J Y 2019 Phys. Rev. B 99 235112Google Scholar
[90] Kawabata K, Shiozaki K, Ueda M, Sato M 2019 Phys. Rev. X 9 041015Google Scholar
[91] Kawabata K, Okuma N, Sato M 2020 Phys. Rev. B 101 195147Google Scholar
[92] Bergholtz E J, Budich J C, Kunst F K 2021 Rev. Mod. Phys. 93 015005Google Scholar
[93] Liu C H, Zhang K, Yang Z, Chen S 2020 Phys. Rev. Res. 2 043167Google Scholar
[94] He P, Liu Y G, Wang J T, Zhu S L 2022 Phys. Rev. A 105 023311Google Scholar
[95] Niu X, Li J, Wu S L, Yi X X 2023 Phys. Rev. A 108 032214Google Scholar
[96] Nakanishi Y, Sasamoto T 2022 Phys. Rev. A 105 022219Google Scholar
[97] Li T, Zhang Y S, Yi W 2022 Phys. Rev. B 105 125111Google Scholar
[98] Purkayastha A 2022 Phys. Rev. A 105 062204Google Scholar
[99] Roccati F, Palma G M, Ciccarello F, Bagarello F 2022 Open Syst. Inf. Dyn. 29 2250004Google Scholar
[100] Niu X, Wang J 2023 arXiv: 2211.04233 [quant-ph
[101] Liu Y G, Chen S 2023 Phys. Rev. B 107 134307Google Scholar
[102] Strang G 2023 Introduction to Linear Algebra (6th Ed.) (Wellesley-Cambridge Press
[103] Demange G, Graefe E M 2011 J. Phys. A: Math. Theor. 45 025303Google Scholar
[104] Rotter I 2009 J. Phys. A: Math. Theor 42 153001Google Scholar
[105] Mostafazadeh A 2002 J. Math. Phys. 43 205Google Scholar
[106] Mostafazadeh A 2002 J. Math. Phys. 43 3944Google Scholar
[107] Xiong Y 2018 J. Phys. Commun. 2 035043Google Scholar
[108] Ryu S, Schnyder A P, Furusaki A, Ludwig A W W 2010 New J. Phys 12 065010Google Scholar
[109] Xu Z C, Zhou Z, Cheng E, Lang L J, Zhu S L 2022 Sci. China: Phys., Mech. Astron. 65 283011Google Scholar
[110] Chen W, Kaya Özdemir Ş, Zhao G, Wiersig J, Yang L 2017 Nature 548 192Google Scholar
[111] Wiersig J 2014 Phys. Rev. Lett. 112 203901Google Scholar
[112] Kato T 1995 Perturbation Theory for Linear Operators (Berlin, Heidelberg: Springer
[113] Schindler J, Li A, Zheng M C, Ellis F M, Kottos T 2011 Phys. Rev. A 84 040101Google Scholar
[114] Dong Z, Li Z, Yang F, Qiu C W, Ho J S 2019 Nat. Electron. 2 335Google Scholar
[115] Zhou B B, Wang L F, Dong L, Huang Q A 2021 J. Phys. Commun. 5 045010Google Scholar
[116] Zeng C, Zhu K, Sun Y, Li G, Guo Z, Jiang J, Li Y, Jiang H, Yang Y, Chen H 2021 New J. Phys 23 063008Google Scholar
[117] 成恩宏, 郎利君 2022 物理学报 71 160301Google Scholar
Cheng E H, Lang L J 2022 Acta Phys. Sin. 71 160301Google Scholar
[118] Hofmann T, Helbig T, Schindler F, Salgo N, Brzezińska M, Greiter M, Kiessling T, Wolf D, Vollhardt A, Kabaši A, Lee C H, Bilušić A, Thomale R, Neupert T 2020 Phys. Rev. Res. 2 023265Google Scholar
[119] Liu S, Shao R, Ma S, Zhang L, You O, Wu H, Xiang Y J, Cui T J, Zhang S 2021 Research 2 5608038Google Scholar
[120] Stegmaier A, Imhof S, Helbig T, Hofmann T, Lee C H, Kremer M, Fritzsche A, Feichtner T, Klembt S, Höfling S, Boettcher I, Fulga I C, Ma L, Schmidt O G, Greiter M, Kiessling T, Szameit A, Thomale R 2021 Phys. Rev. Lett. 126 215302Google Scholar
[121] Zou D, Chen T, He W, Bao J, Lee C H, Sun H, Zhang X 2021 Nat. Commun. 12 7201Google Scholar
[122] Wu J, Huang X, Yang Y, Deng W, Lu J, Deng W, Liu Z 2022 Phys. Rev. B 105 195127Google Scholar
[123] Zhang W, Di F, Yuan H, Wang H, Zheng X, He L, Sun H, Zhang X 2022 Phys. Rev. B 105 195131Google Scholar
[124] Shang C, Liu S, Shao R, Han P, Zang X, Zhang X, Salama K N, Gao W, Lee C H, Thomale R, Manchon A, Zhang S, Cui T J, Schwingenschlögl U 2022 Adv. Sci. 9 2202922Google Scholar
[125] Zhang X, Zhang B, Zhao W, Lee C H 2022 arXiv: 2211.09152 [cond-mat.mes-hall
[126] Zhu P, Sun X Q, Hughes T L, Bahl G 2023 Nat. Commun. 14 720Google Scholar
[127] Wu M, Zhao Q, Kang L, Weng M, Chi Z, Peng R, Liu J, Werner D H, Meng Y, Zhou J 2023 Phys. Rev. B 107 064307Google Scholar
[128] Kim Y, Park H C, Kyung M, Lee K, Ryu J W, You O, Zhang S, Min B, Park M J 2023 arXiv: 2303.13721 [cond-mat.mes-hall
[129] Su L, Guo C X, Wang Y, Li L, Ruan X, Du Y, Chen S, Zheng D 2023 Chin. Phys. B 32 038401Google Scholar
[130] Yin K, Hao X, Huang Y, Zou J, Ma X, Dong T 2023 Phys. Rev. Applied 20 L021003Google Scholar
[131] Ezawa M 2019 Phys. Rev. B 99 121411Google Scholar
[132] Zhou B B, Deng W J, Wang L F, Dong L, Huang Q A 2020 Phys. Rev. Appl. 13 064022Google Scholar
[133] Schindler J, Lin Z, Lee J M, Ramezani H, Ellis F M, Kottos T 2012 J. Phys. A: Math. Theor. 45 444029Google Scholar
[134] Ramezani H, Schindler J, Ellis F M, Günther U, Kottos T 2012 Phys. Rev. A 85 062122Google Scholar
[135] Wang T, Fang J, Xie Z, Dong N, Joglekar Y N, Wang Z, Li J, Luo L 2020 Eur. Phys. J. D 74 167Google Scholar
[136] Chen W K 2002 The Circuits and Filters Handbook (CRC Press
[137] Hofmann T, Helbig T, Lee C H, Greiter M, Thomale R 2019 Phys. Rev. Lett. 122 247702Google Scholar
[138] Zhang H, Chen T, Li L, Lee C H, Zhang X 2023 Phys. Rev. B 107 085426Google Scholar
[139] Yuan H, Zhang W, Zhou Z, Wang W, Pan N, Feng Y, Sun H, Zhang X 2023 Adv. Sci. 2301128Google Scholar
[140] Chen P Y, Sakhdari M, Hajizadegan M, Cui Q, Cheng M M C, El-Ganainy R, Alù A 2018 Nat. Electron. 1 297Google Scholar
[141] Sakhdari M, Hajizadegan M, Zhong Q, Christodoulides D N, El-Ganainy R, Chen P Y 2019 Phys. Rev. Lett. 123 193901Google Scholar
[142] Yin K, Huang Y, Ma C, Hao X, Gao X, Ma X, Dong T 2022 Appl. Phys. Lett. 120 194101Google Scholar
[143] Wu F Y 2004 J. Phys. A: Math. Gen. 37 6653Google Scholar
[144] Li Z, Wu J, Huang X, Lu J, Li F, Deng W, Liu Z 2020 Appl. Phys. Lett. 116 263501Google Scholar
[145] Deng W, Chen T, Zhang X 2022 Phys. Rev. Res. 4 033109Google Scholar
[146] Tang C, Yang H, Song L, Yao X, Yan P, Cao Y 2023 Phys. Rev. B 108 035410Google Scholar
[147] Lang L J, Weng Y, Zhang Y, Cheng E, Liang Q 2021 Phys. Rev. B 103 014302Google Scholar
[148] Chitsazi M, Li H, Ellis F M, Kottos T 2017 Phys. Rev. Lett. 119 093901Google Scholar
[149] Su L, Jiang H, Wang Z, Chen S, Zheng D 2023 Phys. Rev. B 107 184108Google Scholar
[150] Lee C H, Li L, Gong J 2019 Phys. Rev. Lett. 123 016805Google Scholar
[151] Li L, Lee C H, Gong J 2020 Phys. Rev. Lett. 124 250402Google Scholar
[152] Yin K, Huang Y, Yin W, Hao X, Ma X, Dong T 2022 Phys. Rev. Appl. 18 064020Google Scholar
[153] Zhang M, Sweeney W, Hsu C W, Yang L, Stone A D, Jiang L 2019 Phys. Rev. Lett. 123 180501Google Scholar
[154] Kononchuk R, Cai J, Ellis F, Thevamaran R, Kottos T 2022 Nature 607 697Google Scholar
[155] Langbein W 2018 Phys. Rev. A 98 023805Google Scholar
-
图 1 能量E的实部(实线)和虚部(虚线)随非厄米强度γ的变化. 点线处为PT转变点, 其左侧为PT对称相(白色区域), 右侧为PT对称破缺相(灰色区域)
Figure 1. The real (solid lines) and imaginary (dashed lines) parts of the energy E versus the strength γ of the non-Hermiticity. The dotted line indicates the PT transition point, to the left side of which is the PT symmetric phase (white region) and to the right side of which is the PT-broken phase (gray region).
图 2 (a)非厄米SSH模型[59]分别在开边界(粉色)和周期边界(灰色)条件下的能谱E在复平面的示意图; (b)与(a)中能谱相对应的布里渊区(灰色)和广义布里渊区(粉色)的示意图, 其中β的定义见正文; (c)开边界条件下拓扑边缘态(红色)和趋肤态(灰色)的在实空间的几率分布
$|\psi_i|^2$ 示意图, i为格点标记Figure 2. (a) The sketch of the energy spectra in complex plane for the non-Hermitian SSH model in Ref. [59] respectively under open (pink) and periodic (gray) boundary conditions; (b) the sketch of the Brillouin zone (black) and the generalized Brillouin zone (pink) corresponding to the spectra with the same colors in (a), where the definition of β can be referred to in the main text; (c) the sketch of probability distribution
$|\psi_i|^2$ of the topological end state (red) and the skin bulk states (gray) in real space under open boundary conditions, where i is the site index.图 3 (a)能谱ε随微扰z在二阶EP附近劈裂的示意图, 具有
$\varepsilon \propto z^{1/2}$ 的形式[110]; (b)能谱ε随微扰z在传统的二重简并点附近劈裂的示意图, 具有$\varepsilon \propto z$ 的形式[110]Figure 3. (a) The sketch of energy spectra ε versus the perturbation z around a two-order EP, satisfying
$\varepsilon \propto z^{1/2}$ [110]; (b) the sketch of energy spectra ε versus the perturbation z around a traditional two-fold degenerate point, satisfying$\varepsilon \propto z$ [110].图 5 (a)负阻抗的电路原理图[133], 上下图分别表示对地端口和自由双端口的情况; (b) INIC的电路原理图[137]; (c)电压跟随器的电路原理图[136]. 以上所用运放的等效增益函数可用加法器和乘法器实现
Figure 5. (a) The schematic circuit for negative impedance[133], where the upper and lower panels represent the one-port and two-port cases, respectively; (b) the schematic circuit for INIC[136]; (c) the schematic circuit for a voltage follower[136]. The equivalent gain functions of the operational amplifiers used above can be implemented using adders and multipliers.
图 6 (a) PT对称的二聚体电路, 电路中的增益和损耗部分通过电容或互感耦合; (b) 实验测量的本征频率随着增益/损耗参数
$\gamma/\gamma_{\rm{PT}}$ 的变换, 在$\gamma/\gamma_{\rm{PT}}=1$ 时发生PT对称破缺. 图来源于文献[113], 版权属于美国物理学会Figure 6. (a) Circuit diagram of a PT-symmetric dimer, where the gain and loss parts are capacitively or inductively coupled; (b) experimentally measured eigenfrequencies versus the gain/loss parameter
$\gamma/\gamma_{\rm{PT}}$ , where the PT symmetry breaking occurs at$\gamma/\gamma_{\rm{PT}}=1$ . All figures are adapted from Ref.[113] with the copyright © 2011 by the American Physical Society.图 7 (a)非布洛赫演化的拓扑电路; (b)实验测得的电压分布经拉普拉斯变换后在广义布里渊区中的等频分布; (c)四次方根非厄米SSH模型的电路图. 图(a)和(b)来源于文献[127], 图(c)来源于文献[145]. 版权属于美国物理学会
Figure 7. (a) Topolectrical circuit for the non-Bloch dynamics; (b) the isofrequency contour of the measured voltage distribution through the Laplace transform in the GBZ; (c) circuit diagram of a 4th-root non-Hermitian SSH model. Subfigures (a) and (b) are adapted from Ref. [127], and Figure (c) from Ref. [145]. Copyright © 2023 and 2022 by the American Physical Society.
图 8 (a)单纯由增益/损耗诱导的非厄米拓扑模型; (b)实现(a)中模型的电路图. 图均来源于文献[74], 版权属于美国物理学会
Figure 8. (a) A non-Hermitian topological model whose topology is purely induced by gain/loss; (b) circuit diagram for the realization of the model in (a). All figures are adapted from Ref.[74] with the copyright © 2020 by the American Physical Society.
图 9 (a)非厄米二阶拓扑电路的原胞示意图; (b)能隙间的一维边缘态, 彩色和圆圈分别表示实验测量和理论计算的结果; (c)实验测量的零维角态的电压分布; (d)实验测量的角态、边缘态和体态的阻抗响应随驱动频率的变化. 图均来源于文献[122], 版权属于美国物理学会
Figure 9. (a) Circuit diagram of a unit cell of the non-Hermitian second-order topological electric circuit; (b) one-dimensional gapped edge states, where the color map and the blue circles represent the data from the experiment and the theoretical calculation, respectively; (c) experimental voltage distributions of the zeroth dimensional corner states; (d) experimental impedance responses of corner states, edge states, and bulk states to the driving frequency. All figures are adapted from Ref. [122] with the copyright © 2022 by the American Physical Society.
图 10 (a)传统的被动无线传感器与基于PT对称性的被动无线传感器的电路图; (b)基于PT对称性的被动无线传感器的反射谱的测量结果; (c)具有PT对称性的三阶无线传感器的电路图. 图(a)和(b)来源于文献[132], 图(c)来源于文献[152], 版权属于美国物理学会
Figure 10. (a) Circuit diagrams of a conventional passive wireless sensor and a PT-symmetry-based passive wireless sensor; (b) measured reflection spectra for the PT-symmetry-based passive wireless sensor; (c) circuit diagram of a PT symmetric third-order wireless sensing system. Subfigures (a) and (b) are adapted from Ref. [132] and subfigure (c) from Ref. [152], with copyright © 2020 and 2022 by the American Physical Society.
-
[1] Feynman R P 1982 Int. J. Theor. Phys. 21 467Google Scholar
[2] Greiner M, Mandel O, Esslinger T, Hänsch T W, Bloch I 2002 Nature 415 39Google Scholar
[3] Billy J, Josse V, Zuo Z, Bernard A, Hambrecht B, Lugan P, Clement D, Sanchez-Palencia L, Bouyer P, Aspect A 2008 Nature 453 891Google Scholar
[4] Zhang D W, Zhu Y Q, Zhao Y X, Yan H, Zhu S L 2018 Adv. Phys. 67 253Google Scholar
[5] Friedenauer A, Schmitz H, Glueckert J T, Porras D, Schaetz T 2008 Nat. Phys. 4 757Google Scholar
[6] Barreiro J T, Müller M, Schindler P, Nigg D, Monz T, Chwalla M, Hennrich M, Roos C F, Zoller P, Blatt R 2011 Nature 470 486Google Scholar
[7] Blatt R, Roos C F 2012 Nat. Phys. 8 277Google Scholar
[8] You J Q, Shi X F, Hu X, Nori F 2010 Phys. Rev. B 81 014505Google Scholar
[9] Georgescu I, Ashhab S, Nori F 2014 Rev. Mod. Phys. 86 153Google Scholar
[10] Barends R, Kelly J, Megrant A, Veitia A, Sank D, Jeffrey E, White T C, Mutus J, Fowler A G, Campbell B, Chen Y, Chen Z, Chiaro B, Dunsworth A, Neill C, O’Malley P, Roushan P, Vainsencher A, Wenner J, Korotkov A N, Cleland A N, Martinis J M 2014 Nature 508 500Google Scholar
[11] Buluta I, Nori F 2009 Science 326 108Google Scholar
[12] Ashida Y, Gong Z P, Ueda M 2020 Adv. Phys. 69 249Google Scholar
[13] Raghu S, Haldane F D M 2008 Phys. Rev. A 78 033834Google Scholar
[14] Wang Z, Chong Y, Joannopoulos J D, Soljacic M 2009 Nature 461 772Google Scholar
[15] Rechtsman M C, Zeuner J M, Plotnik Y, Lumer Y, Podolsky D, Dreisow F, Nolte S, Segev M, Szameit A 2013 Nature 496 196Google Scholar
[16] Hafezi M, Mittal S, Fan J, Migdall A, Taylor J M 2013 Nat. Photonics 7 1001Google Scholar
[17] Yan Q, Liu R, Yan Z, Liu B, Chen H, Wang Z, Lu L 2018 Nat. Phys. 14 461Google Scholar
[18] Ozawa T, Price H M, Amo A, Goldman N, Hafezi M, Lu L, Rechtsman M C, Schuster D, Simon J, Zilberberg O, Carusotto I 2019 Rev. Mod. Phys. 91 015006Google Scholar
[19] Hussein M I, Leamy M J, Ruzzene M 2014 Appl. Mech. Rev. 66 040802Google Scholar
[20] Wang P, Lu L, Bertoldi K 2015 Phys. Rev. Lett. 115 104302Google Scholar
[21] Xiao Y X, Ma G, Zhang Z Q, Chan C T 2017 Phys. Rev. Lett. 118 166803Google Scholar
[22] Li F, Huang X, Lu J, Ma J, Liu Z 2018 Nat. Phys. 14 30Google Scholar
[23] Li X X, Li B, Li P B 2020 Phys. Rev. Res. 2 013121Google Scholar
[24] Li X X, Li P B, Li H R, Gao H, Li F L 2021 Phys. Rev. Res. 3 013025Google Scholar
[25] Kane C L, Lubensky T C 2014 Nat. Phys. 10 39Google Scholar
[26] Chen B G g, Upadhyaya N, Vitelli V 2014 Proc. Natl. Acad. Sci. 111 13004Google Scholar
[27] Nash L M, Kleckner D, Read A, Vitelli V, Turner A M, Irvine W T M 2015 Proc. Natl. Acad. Sci. 112 14495Google Scholar
[28] Chaunsali R, Kim E, Thakkar A, Kevrekidis P G, Yang J 2017 Phys. Rev. Lett. 119 024301Google Scholar
[29] Brandenbourger M, Locsin X, Lerner E, Coulais C 2019 Nat. Commun. 10 4608Google Scholar
[30] Ghatak A, Brandenbourger M, van Wezel J, Coulais C 2020 Proc. Natl. Acad. Sci. 117 29561Google Scholar
[31] 罗开发, 余睿 2019 物理学报 68 220305Google Scholar
Luo K F, Yu R 2019 Acta Phys. Sin. 68 220305Google Scholar
[32] 刘硕, 张霜, 崔铁军 2021 中国光学 14 736Google Scholar
Liu S, Zhang S, Cui T J 2021 Chinese Optics 14 736Google Scholar
[33] Albert V V, Glazman L I, Jiang L 2015 Phys. Rev. Lett. 114 173902Google Scholar
[34] Ningyuan J, Owens C, Sommer A, Schuster D, Simon J 2015 Phys. Rev. X 5 021031Google Scholar
[35] Lee C H, Imhof S, Berger C, Bayer F, Brehm J, Molenkamp L W, Kiessling T, Thomale R 2018 Commun. Phys. 1 39Google Scholar
[36] Imhof S, Berger C, Bayer F, Brehm J, Molenkamp L W, Kiessling T, Schindler F, Lee C H, Greiter M, Neupert T, Thomale R 2018 Nat. Phys. 14 925Google Scholar
[37] Ezawa M 2018 Phys. Rev. B 98 201402Google Scholar
[38] Helbig T, Hofmann T, Lee C H, Thomale R, Imhof S, Molenkamp L W, Kiessling T 2019 Phys. Rev. B 99 161114Google Scholar
[39] Liu S, Gao W, Zhang Q, Ma S, Zhang L, Liu C, Xiang Y J, Cui T J, Zhang S 2019 Research 2019 8609875Google Scholar
[40] Bao J, Zou D, Zhang W, He W, Sun H, Zhang X 2019 Phys. Rev. B 100 201406Google Scholar
[41] Helbig T, Hofmann T, Imhof S, Abdelghany M, Kiessling T, Molenkamp L W, Lee C H, Szameit A, Greiter M, Thomale R 2020 Nat. Phys. 16 747Google Scholar
[42] Lee C H, Sutrisno A, Hofmann T, Helbig T, Liu Y, Ang Y S, Ang L K, Zhang X, Greiter M, Thomale R 2020 Nat. Commun. 11 4385Google Scholar
[43] Wu J, Huang X, Lu J, Wu Y, Deng W, Li F, Liu Z 2020 Phys. Rev. B 102 104109Google Scholar
[44] Kotwal T, Moseley F, Stegmaier A, Imhof S, Brand H, Kießling T, Thomale R, Ronellenfitsch H, Dunkel J 2021 Proc. Natl. Acad. Sci. 118 e2106411118Google Scholar
[45] Dong J, Juričić V, Roy B 2021 Phys. Rev. Res. 3 023056Google Scholar
[46] Yang Y, Zhu D, Hang Z, Chong Y 2021 Sci. China: Phys., Mech. Astron. 64 257011Google Scholar
[47] Zhang W, Zou D, Pei Q, He W, Bao J, Sun H, Zhang X 2021 Phys. Rev. Lett. 126 146802Google Scholar
[48] Wu J, Wang Z, Biao Y, Fei F, Zhang S, Yin Z, Hu Y, Song Z, Wu T, Song F, Yu R 2022 Nat. Electron. 5 635Google Scholar
[49] Wang Z, Zeng X T, Biao Y, Yan Z, Yu R 2023 Phys. Rev. Lett. 130 057201Google Scholar
[50] Breuer H P, Petruccione F 2002 The Theory of Open Quantum Systems (Oxford: Oxford University Press
[51] Daley A J 2014 Adv. Phys. 63 77Google Scholar
[52] Bender C M, Boettcher S 1998 Phys. Rev. Lett. 80 5243Google Scholar
[53] Bender C M 2007 Rep. Prog. Phys. 70 947Google Scholar
[54] Longhi S 2019 Opt. Lett. 44 5804Google Scholar
[55] Pan J S, Yi W, Gong J 2021 Commun. Phys. 4 261Google Scholar
[56] Moca C P, Sticlet D, Dóra B, Zaránd G 2023 Phys. Rev. B 107 115111Google Scholar
[57] Lee T E 2016 Phys. Rev. Lett. 116 133903Google Scholar
[58] Kunst F K, Edvardsson E, Budich J C, Bergholtz E J 2018 Phys. Rev. Lett. 121 026808Google Scholar
[59] Yao S, Wang Z 2018 Phys. Rev. Lett. 121 086803Google Scholar
[60] Yao S, Song F, Wang Z 2018 Phys. Rev. Lett. 121 136802Google Scholar
[61] Edvardsson E, Kunst F K, Bergholtz E J 2019 Phys. Rev. B 99 081302Google Scholar
[62] Ghatak A, Das T 2019 J. Phys.: Condens. Matter 31 263001Google Scholar
[63] Song F, Yao S, Wang Z 2019 Phys. Rev. Lett. 123 170401Google Scholar
[64] Song F, Yao S, Wang Z 2019 Phys. Rev. Lett. 123 246801Google Scholar
[65] Borgnia D S, Kruchkov A J, Slager R J 2020 Phys. Rev. Lett. 124 056802Google Scholar
[66] Xiao L, Deng T, Wang K, Zhu G, Wang Z, Yi W, Xue P 2020 Nat. Phys. 16 761Google Scholar
[67] Zhang K, Yang Z, Fang C 2020 Phys. Rev. Lett. 125 126402Google Scholar
[68] Yang Z, Zhang K, Fang C, Hu J 2020 Phys. Rev. Lett. 125 226402Google Scholar
[69] Okuma N, Sato M 2021 Phys. Rev. Lett. 126 176601Google Scholar
[70] Liang Q, Xie D, Dong Z, Li H, Li H, Gadway B, Yi W, Yan B 2022 Phys. Rev. Lett. 129 070401Google Scholar
[71] Milburn T J, Doppler J, Holmes C A, Portolan S, Rotter S, Rabl P 2015 Phys. Rev. A 92 052124Google Scholar
[72] Zhong Q, Christodoulides D N, Khajavikhan M, Makris K G, El-Ganainy R 2018 Phys. Rev. A 97 020105Google Scholar
[73] Zhang K L, Wang P, Song Z 2019 Phys. Rev. A 99 042111Google Scholar
[74] Liu S, Ma S, Yang C, Zhang L, Gao W, Xiang Y J, Cui T J, Zhang S 2020 Phys. Rev. Appl. 13 014047Google Scholar
[75] Heiss W D 2004 J. Phys. A: Math. Gen. 37 2455Google Scholar
[76] Moiseyev N 2011 Non-Hermitian Quantum Mechanics (Cambridge: Cambridge University Press
[77] Okugawa R, Takahashi R, Yokomizo K 2020 Phys. Rev. B 102 241202Google Scholar
[78] Zhang X, Tian Y, Jiang J H, Lu M H, Chen Y F 2021 Nat. Commun. 12 5377Google Scholar
[79] Zhang X, Zhang T, Lu M H, Chen Y F 2022 Adv. Phys.: X 7 2109431Google Scholar
[80] Deng T S, Yi W 2019 Phys. Rev. B 100 035102Google Scholar
[81] Yokomizo K, Murakami S 2019 Phys. Rev. Lett. 123 066404Google Scholar
[82] Zhu X, Wang H, Gupta S K, Zhang H, Xie B, Lu M, Chen Y 2020 Phys. Rev. Res. 2 013280Google Scholar
[83] Yokomizo K, Murakami S 2020 Prog. Theor. Exp. Phys. 2020 12A102Google Scholar
[84] Xiao L, Deng T, Wang K, Wang Z, Yi W, Xue P 2021 Phys. Rev. Lett. 126 230402Google Scholar
[85] Wang K, Li T, Xiao L, Han Y, Yi W, Xue P 2021 Phys. Rev. Lett. 127 270602Google Scholar
[86] Yokomizo K 2022 Non-Bloch Band Theory of Non Hermitian Systems (Singapore: Springer Nature
[87] Yokomizo K, Murakami S 2023 Phys. Rev. B 107 195112Google Scholar
[88] Shi K, Tian M, Sun F X, Zhang W 2023 Phys. Rev. B 107 205154Google Scholar
[89] Zhou H, Lee J Y 2019 Phys. Rev. B 99 235112Google Scholar
[90] Kawabata K, Shiozaki K, Ueda M, Sato M 2019 Phys. Rev. X 9 041015Google Scholar
[91] Kawabata K, Okuma N, Sato M 2020 Phys. Rev. B 101 195147Google Scholar
[92] Bergholtz E J, Budich J C, Kunst F K 2021 Rev. Mod. Phys. 93 015005Google Scholar
[93] Liu C H, Zhang K, Yang Z, Chen S 2020 Phys. Rev. Res. 2 043167Google Scholar
[94] He P, Liu Y G, Wang J T, Zhu S L 2022 Phys. Rev. A 105 023311Google Scholar
[95] Niu X, Li J, Wu S L, Yi X X 2023 Phys. Rev. A 108 032214Google Scholar
[96] Nakanishi Y, Sasamoto T 2022 Phys. Rev. A 105 022219Google Scholar
[97] Li T, Zhang Y S, Yi W 2022 Phys. Rev. B 105 125111Google Scholar
[98] Purkayastha A 2022 Phys. Rev. A 105 062204Google Scholar
[99] Roccati F, Palma G M, Ciccarello F, Bagarello F 2022 Open Syst. Inf. Dyn. 29 2250004Google Scholar
[100] Niu X, Wang J 2023 arXiv: 2211.04233 [quant-ph
[101] Liu Y G, Chen S 2023 Phys. Rev. B 107 134307Google Scholar
[102] Strang G 2023 Introduction to Linear Algebra (6th Ed.) (Wellesley-Cambridge Press
[103] Demange G, Graefe E M 2011 J. Phys. A: Math. Theor. 45 025303Google Scholar
[104] Rotter I 2009 J. Phys. A: Math. Theor 42 153001Google Scholar
[105] Mostafazadeh A 2002 J. Math. Phys. 43 205Google Scholar
[106] Mostafazadeh A 2002 J. Math. Phys. 43 3944Google Scholar
[107] Xiong Y 2018 J. Phys. Commun. 2 035043Google Scholar
[108] Ryu S, Schnyder A P, Furusaki A, Ludwig A W W 2010 New J. Phys 12 065010Google Scholar
[109] Xu Z C, Zhou Z, Cheng E, Lang L J, Zhu S L 2022 Sci. China: Phys., Mech. Astron. 65 283011Google Scholar
[110] Chen W, Kaya Özdemir Ş, Zhao G, Wiersig J, Yang L 2017 Nature 548 192Google Scholar
[111] Wiersig J 2014 Phys. Rev. Lett. 112 203901Google Scholar
[112] Kato T 1995 Perturbation Theory for Linear Operators (Berlin, Heidelberg: Springer
[113] Schindler J, Li A, Zheng M C, Ellis F M, Kottos T 2011 Phys. Rev. A 84 040101Google Scholar
[114] Dong Z, Li Z, Yang F, Qiu C W, Ho J S 2019 Nat. Electron. 2 335Google Scholar
[115] Zhou B B, Wang L F, Dong L, Huang Q A 2021 J. Phys. Commun. 5 045010Google Scholar
[116] Zeng C, Zhu K, Sun Y, Li G, Guo Z, Jiang J, Li Y, Jiang H, Yang Y, Chen H 2021 New J. Phys 23 063008Google Scholar
[117] 成恩宏, 郎利君 2022 物理学报 71 160301Google Scholar
Cheng E H, Lang L J 2022 Acta Phys. Sin. 71 160301Google Scholar
[118] Hofmann T, Helbig T, Schindler F, Salgo N, Brzezińska M, Greiter M, Kiessling T, Wolf D, Vollhardt A, Kabaši A, Lee C H, Bilušić A, Thomale R, Neupert T 2020 Phys. Rev. Res. 2 023265Google Scholar
[119] Liu S, Shao R, Ma S, Zhang L, You O, Wu H, Xiang Y J, Cui T J, Zhang S 2021 Research 2 5608038Google Scholar
[120] Stegmaier A, Imhof S, Helbig T, Hofmann T, Lee C H, Kremer M, Fritzsche A, Feichtner T, Klembt S, Höfling S, Boettcher I, Fulga I C, Ma L, Schmidt O G, Greiter M, Kiessling T, Szameit A, Thomale R 2021 Phys. Rev. Lett. 126 215302Google Scholar
[121] Zou D, Chen T, He W, Bao J, Lee C H, Sun H, Zhang X 2021 Nat. Commun. 12 7201Google Scholar
[122] Wu J, Huang X, Yang Y, Deng W, Lu J, Deng W, Liu Z 2022 Phys. Rev. B 105 195127Google Scholar
[123] Zhang W, Di F, Yuan H, Wang H, Zheng X, He L, Sun H, Zhang X 2022 Phys. Rev. B 105 195131Google Scholar
[124] Shang C, Liu S, Shao R, Han P, Zang X, Zhang X, Salama K N, Gao W, Lee C H, Thomale R, Manchon A, Zhang S, Cui T J, Schwingenschlögl U 2022 Adv. Sci. 9 2202922Google Scholar
[125] Zhang X, Zhang B, Zhao W, Lee C H 2022 arXiv: 2211.09152 [cond-mat.mes-hall
[126] Zhu P, Sun X Q, Hughes T L, Bahl G 2023 Nat. Commun. 14 720Google Scholar
[127] Wu M, Zhao Q, Kang L, Weng M, Chi Z, Peng R, Liu J, Werner D H, Meng Y, Zhou J 2023 Phys. Rev. B 107 064307Google Scholar
[128] Kim Y, Park H C, Kyung M, Lee K, Ryu J W, You O, Zhang S, Min B, Park M J 2023 arXiv: 2303.13721 [cond-mat.mes-hall
[129] Su L, Guo C X, Wang Y, Li L, Ruan X, Du Y, Chen S, Zheng D 2023 Chin. Phys. B 32 038401Google Scholar
[130] Yin K, Hao X, Huang Y, Zou J, Ma X, Dong T 2023 Phys. Rev. Applied 20 L021003Google Scholar
[131] Ezawa M 2019 Phys. Rev. B 99 121411Google Scholar
[132] Zhou B B, Deng W J, Wang L F, Dong L, Huang Q A 2020 Phys. Rev. Appl. 13 064022Google Scholar
[133] Schindler J, Lin Z, Lee J M, Ramezani H, Ellis F M, Kottos T 2012 J. Phys. A: Math. Theor. 45 444029Google Scholar
[134] Ramezani H, Schindler J, Ellis F M, Günther U, Kottos T 2012 Phys. Rev. A 85 062122Google Scholar
[135] Wang T, Fang J, Xie Z, Dong N, Joglekar Y N, Wang Z, Li J, Luo L 2020 Eur. Phys. J. D 74 167Google Scholar
[136] Chen W K 2002 The Circuits and Filters Handbook (CRC Press
[137] Hofmann T, Helbig T, Lee C H, Greiter M, Thomale R 2019 Phys. Rev. Lett. 122 247702Google Scholar
[138] Zhang H, Chen T, Li L, Lee C H, Zhang X 2023 Phys. Rev. B 107 085426Google Scholar
[139] Yuan H, Zhang W, Zhou Z, Wang W, Pan N, Feng Y, Sun H, Zhang X 2023 Adv. Sci. 2301128Google Scholar
[140] Chen P Y, Sakhdari M, Hajizadegan M, Cui Q, Cheng M M C, El-Ganainy R, Alù A 2018 Nat. Electron. 1 297Google Scholar
[141] Sakhdari M, Hajizadegan M, Zhong Q, Christodoulides D N, El-Ganainy R, Chen P Y 2019 Phys. Rev. Lett. 123 193901Google Scholar
[142] Yin K, Huang Y, Ma C, Hao X, Gao X, Ma X, Dong T 2022 Appl. Phys. Lett. 120 194101Google Scholar
[143] Wu F Y 2004 J. Phys. A: Math. Gen. 37 6653Google Scholar
[144] Li Z, Wu J, Huang X, Lu J, Li F, Deng W, Liu Z 2020 Appl. Phys. Lett. 116 263501Google Scholar
[145] Deng W, Chen T, Zhang X 2022 Phys. Rev. Res. 4 033109Google Scholar
[146] Tang C, Yang H, Song L, Yao X, Yan P, Cao Y 2023 Phys. Rev. B 108 035410Google Scholar
[147] Lang L J, Weng Y, Zhang Y, Cheng E, Liang Q 2021 Phys. Rev. B 103 014302Google Scholar
[148] Chitsazi M, Li H, Ellis F M, Kottos T 2017 Phys. Rev. Lett. 119 093901Google Scholar
[149] Su L, Jiang H, Wang Z, Chen S, Zheng D 2023 Phys. Rev. B 107 184108Google Scholar
[150] Lee C H, Li L, Gong J 2019 Phys. Rev. Lett. 123 016805Google Scholar
[151] Li L, Lee C H, Gong J 2020 Phys. Rev. Lett. 124 250402Google Scholar
[152] Yin K, Huang Y, Yin W, Hao X, Ma X, Dong T 2022 Phys. Rev. Appl. 18 064020Google Scholar
[153] Zhang M, Sweeney W, Hsu C W, Yang L, Stone A D, Jiang L 2019 Phys. Rev. Lett. 123 180501Google Scholar
[154] Kononchuk R, Cai J, Ellis F, Thevamaran R, Kottos T 2022 Nature 607 697Google Scholar
[155] Langbein W 2018 Phys. Rev. A 98 023805Google Scholar
Catalog
Metrics
- Abstract views: 4639
- PDF Downloads: 236