APP下载

矩阵形式的不变本征算符方法以及几种介观电路的本征频率*

2019-11-28吴泽范洪义

物理学报 2019年22期
关键词:本征拉格朗矩阵

吴泽 范洪义

1) (中国科学技术大学近代物理系,合肥 230026)

2) (中国科学技术大学材料科学与工程系,合肥 230026)

本文把不变本征算符方法(invariant eigen-operator,IEO方法)推广到了基于拉格朗日量的矩阵形式,将以往计算的思路和过程用简约的矩阵形式表示出来,这对大规模复杂多回路的介观电路的计算有着重要的意义.此外用该方法计算了三个L-C介观电路的本征频率,包括存在互感和不存在互感的两种情形.通过计算结果得出了这些电路的相关性质,说明了本征频率只与介观电路本身的元件性质有关.

1 引 言

介观电路的量子化是量子电路领域的重要课题之一.众所周知,经典 L-C 耦合电路的频率可以通过Kirchhoff定律写出久期方程求解得到;尽管量子电路的频率往往跟对应的经典电路频率无异,但是将经典电路的方法直接推广到量子介观电路是要十分谨慎的.其中一个重要的原因是量子化后的介观电路所有的观测量都是算符,而不再是普通的C-数,算符的非对易性质会导致一些与经典电路不一样的结论.所以处理介观量子电路应该使用量子力学的方法,不变本征算符(invariant eigenoperator,IEO)方法就是此类方法之一[1-5].

原始的IEO方法要先写出介观电路的拉式量、确定广义坐标,再计算共轭动量和哈密顿量;接着假设一个参数待定的不变本征算符,代入二阶对易关系求解出参数并确定不变本征算符的具体形式,最后再利用IEO方程求解出介观电路的本征频率[6-8].我们发现每次求解频率都要经历这一系列程式化的流程,而且随着电路规模的增大,该过程将会变得愈加复杂[9].于是本文利用指标表示的方法将这一系列过程收束为简单而优雅的矩阵表示,即提出了基于拉格朗日量的IEO方法的矩阵形式,该形式能够在很大程度上简化IEO方法的计算过程.基于该形式计算了三种 L-C 电路的本征频率.

2 基于拉格朗日量的IEO方法的矩阵表示

假设多回路介观电路的回路电荷为 Qi,电流为 Ii,下标代表不同的回路.则其拉格朗日量可以写为如下的矩阵形式[10,11]

把电荷 Qi当成广义坐标,则对应的共轭动量为

它的逆表示

有了哈密顿量就可以把介观电路量子化,此时的Pi,Qj量子化后应满足正则对易关系(约定 ℏ=1 )

下面利用IEO方法计算频率,设不变本征算符为

t是待定的系数,则对易子:

二阶对易子:

考虑IEO方程

则有

把指标形式换成向量-矩阵的形式:

这意味着 ω2是矩阵 M-1N 的本征值;这样一来只要知道了多回路介观电路的拉格朗日量,就能很方便地利用这个结论求解出介观电路的本征频率.下面举几个例子来验证我们的结论.

3 矩阵IEO方法计算两回路L-C介观电路特征频率

首先考虑图1所示的 L-C介观电路[12-15].

图1 两回路 L- C 介观电路Fig.1.Two-loop L- C mesoscopic circuit.

用IEO方法求解没有耦合的频率,图1电路的拉格朗日量为

对应IEO方法中的矩阵

套用IEO公式

解得

其中根号下的Δ具体形式为

对称情形可以得到简洁的解,令 l1=l2=l,c1=c2=c3=c,则

该结果不依赖回路电流的方向,因为即便上述回路电流反向,也仅仅是改变了矩阵两个反对角元的符号(正号变负号),而求解本征值的时候,这两个负号又恰好抵消.再来看图2所示的电路[16,17].

图2 无互感的 L- C 介观电路Fig.2.L- C mesoscopic circuit without mutual inductance.

图2电路的拉格朗日量为

对应IEO方法中的矩阵

套用IEO公式计算本征值可得本征频率

其中

同时我们发现拉格朗日量中的 I1±I2代表两个回路电路方向一致或是相反,对应M矩阵的非对角项的正负.注意到矩阵的本征值跟矩阵的非对角项的正负无关,所以这并不会影响本征频率,这也是“本征”二字的体现——是电路的固有特征.

此外我们还考虑了存在互感的电路[18],如图3所示,l1和l2的互感为 2m.它的拉格朗日量为

对应IEO方法中的矩阵

同样可得本征频率

其中

同样可以写出对称形式的频率,令 l1=l2=l3=l,c1=c2=c 得

图3 带互感的 L- C 介观电路Fig.3.L- C mesoscopic circuit with mutual inductance 2m.

此时要求 1 -c和l -m同号,ω-才有可能存在,否则只存在一个特征频率.

4 结 论

不变本征算符方法不失为一种计算介观量子电路特征频率的好办法,此外该方法在光-原子体系[19-21]以及固体物理[22]中也有很多应用,而本文对该方法的改进使得其中繁琐的程式化的步骤变成了简单的求解矩阵本征值的过程.考察的三个量子介观电路能够验证这种方法的方便快捷之处,也从这种形式中体现了“本征频率”只与该介观电路本身的属性有关.在该方法的基础上,更多的关于量子介观电路的研究得以顺利展开.

猜你喜欢

本征拉格朗矩阵
Generative Adversarial Network Based Heuristics for Sampling-Based Path Planning
基于本征正交分解的水平轴风力机非定常尾迹特性分析
基于APDL 语言的本征应变法重构激光冲击强化后的残余应力场
这样的完美叫“自私”
拉格朗日的“自私”
拉格朗日的“自私”
探索Euler图的等价命题
这样的完美叫“自私”
多项式理论在矩阵求逆中的应用
矩阵