搜索

x

留言板

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

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

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

吴泽 范洪义

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

吴泽, 范洪义
PDF
HTML
导出引用
  • 本文把不变本征算符方法(invariant eigen-operator, IEO方法)推广到了基于拉格朗日量的矩阵形式, 将以往计算的思路和过程用简约的矩阵形式表示出来, 这对大规模复杂多回路的介观电路的计算有着重要的意义. 此外用该方法计算了三个LC介观电路的本征频率, 包括存在互感和不存在互感的两种情形. 通过计算结果得出了这些电路的相关性质, 说明了本征频率只与介观电路本身的元件性质有关.
      通信作者: 范洪义, fhym@ustc.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11775208)资助的课题
    [1]

    Fan H Y, Li C 2004 Phys. Lett. A 321 75

    [2]

    任益充, 范洪义 2013 物理学报 62 156301

    Ren Y C, Fan H Y 2013 Acta Phys. Sin. 62 156301

    [3]

    范洪义, 吴昊, 袁洪春 2011 量子力学的不变本征算符方法 (上海: 上海交通大学出版社) 第175−178页

    Fan H Y, Wu H, Yuan H C 2011 Invariant Eigen-Operator Method in Quantum Mechanics (Shanghai: Shanghai Jiao Tong University Press) pp175−178 (in Chinese)

    [4]

    Fan H Y, Wu H, Xu X F 2005 Int. J. Mod. Phys. B 19 4073

    [5]

    Fan H Y, Tang X B 2007 Commun. Theor. Phys. 47 865

    [6]

    Ren G, Fan H Y 2009 Int. J. Phys. 48 2016

    [7]

    Fan H Y, Wang T T 2007 Int. J. Mod. Phys. B 21 1961

    [8]

    Song T Q, Fan H Y 2006 Mod. Phys. Lett. A 21 451

    [9]

    Liu Y M 2009 Int. J. Theor. Phys. 48 2372

    [10]

    Fan H Y, Tang X B 2006 Commun. Theor. Phys. 46 603

    [11]

    Fan H Y, Tang X B, Hu H P 2008 Commun. Theor. Phys. 50 674

    [12]

    Fan H Y, Wu H 2007 Mod. Phys. Lett. B 21 1751

    [13]

    Fan H Y, Wu H 2008 Commun. Theor. Phys. 49 50

    [14]

    Tang X B, Fan H Y 2010 Int. J. Theor. Phys. 49 877

    [15]

    Fan H Y, Wu H 2009 Int. J. Mod. Phys. B 23 234

    [16]

    Fan H Y, Wang T T, Yang Y L 2006 Int. J. Mod. Phys. B 20 5417

    [17]

    Song T Q, Zhu J Y, Fan H Y 2003 Commun. Theor. Phys. 39 91

    [18]

    Fan H Y, Tang X B 2006 Commun. Theor. Phys. 45 1003

    [19]

    Fan H Y, Gui W J, Gui J Z 2007 Int. J. Mod. Phys. B 21 737

    [20]

    Yu T X, Fan H Y 2010 Commun. Theor. Phys. 53 257

    [21]

    Meng X G, Wang J S, Zhai Y, Fan H Y 2008 Chin. Phys. Lett. 25 1205

    [22]

    张科, 范承玉, 范洪义 2018 物理学报 67 170301

    Zhang K, Fan C Y, Fan H Y, 2018 Acta Phys. Sin. 67 170301

  • 图 1  两回路$ L\text{-}C$介观电路

    Fig. 1.  Two-loop $L\text{-}C$ mesoscopic circuit

    图 2  无互感的$ L\text-C$介观电路

    Fig. 2.  $ L\text{-}C$ mesoscopic circuit without mutual inductance

    图 3  带互感的$ L\text-C$介观电路

    Fig. 3.  $ L\text-C$ mesoscopic circuit with mutual inductance 2m.

  • [1]

    Fan H Y, Li C 2004 Phys. Lett. A 321 75

    [2]

    任益充, 范洪义 2013 物理学报 62 156301

    Ren Y C, Fan H Y 2013 Acta Phys. Sin. 62 156301

    [3]

    范洪义, 吴昊, 袁洪春 2011 量子力学的不变本征算符方法 (上海: 上海交通大学出版社) 第175−178页

    Fan H Y, Wu H, Yuan H C 2011 Invariant Eigen-Operator Method in Quantum Mechanics (Shanghai: Shanghai Jiao Tong University Press) pp175−178 (in Chinese)

    [4]

    Fan H Y, Wu H, Xu X F 2005 Int. J. Mod. Phys. B 19 4073

    [5]

    Fan H Y, Tang X B 2007 Commun. Theor. Phys. 47 865

    [6]

    Ren G, Fan H Y 2009 Int. J. Phys. 48 2016

    [7]

    Fan H Y, Wang T T 2007 Int. J. Mod. Phys. B 21 1961

    [8]

    Song T Q, Fan H Y 2006 Mod. Phys. Lett. A 21 451

    [9]

    Liu Y M 2009 Int. J. Theor. Phys. 48 2372

    [10]

    Fan H Y, Tang X B 2006 Commun. Theor. Phys. 46 603

    [11]

    Fan H Y, Tang X B, Hu H P 2008 Commun. Theor. Phys. 50 674

    [12]

    Fan H Y, Wu H 2007 Mod. Phys. Lett. B 21 1751

    [13]

    Fan H Y, Wu H 2008 Commun. Theor. Phys. 49 50

    [14]

    Tang X B, Fan H Y 2010 Int. J. Theor. Phys. 49 877

    [15]

    Fan H Y, Wu H 2009 Int. J. Mod. Phys. B 23 234

    [16]

    Fan H Y, Wang T T, Yang Y L 2006 Int. J. Mod. Phys. B 20 5417

    [17]

    Song T Q, Zhu J Y, Fan H Y 2003 Commun. Theor. Phys. 39 91

    [18]

    Fan H Y, Tang X B 2006 Commun. Theor. Phys. 45 1003

    [19]

    Fan H Y, Gui W J, Gui J Z 2007 Int. J. Mod. Phys. B 21 737

    [20]

    Yu T X, Fan H Y 2010 Commun. Theor. Phys. 53 257

    [21]

    Meng X G, Wang J S, Zhai Y, Fan H Y 2008 Chin. Phys. Lett. 25 1205

    [22]

    张科, 范承玉, 范洪义 2018 物理学报 67 170301

    Zhang K, Fan C Y, Fan H Y, 2018 Acta Phys. Sin. 67 170301

  • 引用本文:
    Citation:
计量
  • 文章访问数:  1684
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-03
  • 修回日期:  2019-08-30
  • 上网日期:  2019-11-01
  • 刊出日期:  2019-11-20

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

  • 1. 中国科学技术大学近代物理系, 合肥 230026
  • 2. 中国科学技术大学材料科学与工程系, 合肥 230026
  • 通信作者: 范洪义, fhym@ustc.edu.cn
    基金项目: 国家自然科学基金(批准号: 11775208)资助的课题

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

English Abstract

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

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

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

      $ \begin{split} {L} =\, & \frac{1}{2}\left({{I_1}}\quad {{I_2}}\quad {{I_3}}\quad \cdots\right){{M}}\left( {\begin{aligned} {{I_1}}\\ {{I_2}}\\ {{I_3}}\\ \vdots\; \end{aligned}} \right) \\ &- \frac{1}{2}\left({{Q_1}}\quad {{Q_2}}\quad {{Q_3}}\quad \cdots\right){{N}}\left( {\begin{aligned} {{Q_1}}\\ {{Q_2}}\\ {{Q_3}}\\ \vdots \;\; \end{aligned}} \right)\\ \equiv\, & \frac{1}{2}{I_i}{M_{ij}}{I_j} - \frac{1}{2}{Q_i}{N_{ij}}{Q_j}, \end{split} $

      其中重复指标代表求和, $ I_{i} = \dot{Q}_{i} $, M, N都为对称矩阵, 即有:

      $ M_{ij} = M_{ji},\; N_{ij} = N_{ji}, $

      把电荷$ Q_{i} $当成广义坐标, 则对应的共轭动量为

      $ P_{i} = \frac{\partial L}{\partial \dot{Q}_{i}} = M_{ij}I_{j}, $

      它的逆表示

      $ I_{i} = M^{-1}_{ij}P_{j}, $

      这里的$ M_{ij}^{-1} $表示$ { M}^{-1} $的分量. 那么哈密顿量为

      $ \begin{split} {H} \, &= P_{i}\dot{Q}_{i}-{ L} = \frac{1}{2}I_{i}M_{ij}I_{j}+\frac{1}{2}Q_{i}N_{ij}Q_{j}\\ &= \frac{1}{2}P_{i}M_{ij}^{-1}P_{j}+\frac{1}{2}Q_{i}N_{ij}Q_{j} . \end{split} $

      有了哈密顿量就可以把介观电路量子化, 此时的$ P_{i}, Q_{j} $量子化后应满足正则对易关系(约定$ \hbar = 1 $)

      $ \left[Q_{i},P_{j}\right] = {\rm i}\delta_{ij}. $

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

      $ {O_{\rm e}} = \left({{t_1}}\quad {{t_2}}\quad{{t_3}}\quad \cdots\right)\left( {\begin{aligned} {{P_1}}\\ {{P_2}}\\ {{P_3}}\\ \vdots\;\, \end{aligned}} \right) = {t_i}{P_i} ,$

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

      $ \begin{split} \left[ {{O_{\rm e}},2H} \right] \, &= \left[ {{t_i}{P_i},{Q_j}{N_{jk}}{Q_k}} \right]\\ & = {t_i}{N_{jk}}\left( {{Q_j}\left[ {{P_i},{Q_k}} \right] + \left[ {{P_i},{Q_j}} \right]{Q_k}} \right)\\ &= - {\rm i}{t_i}{N_{jk}}\left( {{Q_j}{\delta _{ik}} + {Q_k}{\delta _{ij}}} \right)\\ & = - {\rm i}\left( {{t_i}{N_{ji}}{Q_j} + {t_i}{N_{ik}}{Q_i}} \right)\\ &= - {\rm i}2{t_i}{N_{ij}}{Q_j}, \end{split} $

      二阶对易子:

      $\begin{split} & \frac{1}{{ - 2{\rm i}}}\left[ {\left[ {{O_{\rm e}},2H} \right],2H} \right] \\ = \,& \left[ {{t_i}{N_{ij}}{Q_j},{P_k}M_{kl}^{ - 1}{P_l}} \right]\\ = \,& {t_i}{N_{ij}}M_{kl}^{ - 1}\left( {{P_k}\left[ {{Q_j},{P_l}} \right] + \left[ {{Q_j},{P_k}} \right]{P_l}} \right)\\ = \,& i{t_i}{N_{ij}}M_{kl}^{ - 1}\left( {{P_k}{\delta _{jl}} + {P_l}{\delta _{jk}}} \right)\\ =\,& {\rm i}\left( {{t_i}{N_{ij}}M_{kj}^{ - 1}{P_k} + {t_i}{N_{ij}}M_{jl}^{ - 1}{P_l}} \right)\\ =\,& 2{\rm i}{t_i}{N_{ij}}M_{jk}^{ - 1}{P_k}, \end{split} $

      $ \left[\left[O_{\rm e},H\right],H\right] = t_{i}N_{ij}M^{-1}_{jk}P_{k}, $

      考虑IEO方程

      $ \left[\left[O_{\rm e},H\right],H\right] = \omega^{2}O_{\rm e}, $

      则有

      $ t_{i}N_{ij}M^{-1}_{jk}P_{k} = \omega^{2}t_{k}P_{k}, $

      $ t_{i}N_{ij}M^{-1}_{jk} = \omega^{2}t_{k}\quad\Leftrightarrow\quad M^{-1}_{kj}N_{ji}t_{i} = \omega^{2}t_{i}. $

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

      $ { M}^{-1}{ N}{ {t}} = \omega^{2}{ {t}}. $

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

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

      图  1  两回路$ L\text{-}C$介观电路

      Figure 1.  Two-loop $L\text{-}C$ mesoscopic circuit

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

      $ L = \frac{1}{2}\left( l_{1}I_{1}^{2}+l_{2}I_{2}^{2}\right)-\frac{1}{2}\left[\frac{Q_{1}^{2}}{c_{1}}+\frac{Q_{2}^{2}}{c_{2}}+\frac{\left(Q_{1}+Q_{2}\right)^{2}}{c_{3}}\right], $

      对应IEO方法中的矩阵

      $ {{M}} = \left(\!\!\!{\begin{array}{*{20}{c}} {{l_1}}&0\\ 0&{{l_2}} \end{array}}\!\!\!\right),\;\;\;\;\;{{N}} = \left(\!\!\!{\begin{array}{*{20}{c}} {\dfrac{1}{{{c_1}}} + \dfrac{1}{{{c_3}}}}&{\dfrac{1}{{{c_3}}}}\\ {\dfrac{1}{{{c_3}}}}&{\dfrac{1}{{{c_2}}} + \dfrac{1}{{{c_3}}}} \end{array}}\!\!\!\right). $

      套用IEO公式

      $ { M}^{-1}{ N}{{t}} = \omega^{2}{{t}}, $

      解得

      $ \omega^{2} = \frac{c_{2}c_{3}^{3}l_{1}l_{2}^{2}+c_{1}c_{3}^{3}l_{1}l_{2}\left(c_{3}l_{1}+c_{2}l_{1}+c_{2}l_{2}\right)\pm\sqrt{\varDelta}}{2c_{1}c_{2}c_{3}^{3}l_{1}^{2}l_{2}^{2}}, $

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

      $\begin{split} \varDelta =\, & c_{3}^{4}l_{1}^{2}l_{2}^{2}\left[c_{1}^{2}\left(c_{2}+c_{3}\right)^{2}l_{1}^{2}+c_{2}^{2}\left(c_{1}+c_{3}\right)^{2}l_{2}^{2}\right.\\ & \left.-2c_{1}c_{2}\left(c_{3}^{2}+c_{1}c_{3}+c_{2}c_{3}-c_{1}c_{2}\right)l_{1}l_{2}\right]\end{split}. $

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

      $ \omega_{+} = \sqrt{\frac{3}{cl}},\quad \omega_{-} = \sqrt{\frac{1}{cl}}. $

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

      图  2  无互感的$ L\text-C$介观电路

      Figure 2.  $ L\text{-}C$ mesoscopic circuit without mutual inductance

      图2电路的拉格朗日量为

      $ L = \frac{1}{2}\left[l_{1}I_{1}^{2}+l_{2}I_{2}^{2}+l_{3}\left(I_{1}\pm I_{2}\right)^{2}\right]-\frac{1}{2}\left(\frac{q_{1}^{2}}{c_{1}}+\frac{q_{2}^{2}}{c_{2}}\right) , $

      对应IEO方法中的矩阵

      $ {{M}} = \left( {\begin{array}{*{20}{c}} {{l_1} + {l_3}}&{ \pm {l_3}}\\ { \pm {l_3}}&{{l_2} + {l_3}} \end{array}} \right),\;\;\;\;\;{{N}} = \left( {\begin{array}{*{20}{c}} {\dfrac{1}{{{c_1}}}}&0\\ 0&{\dfrac{1}{{{c_2}}}} \end{array}} \right), $

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

      $ \omega^{2} = \frac{c_{1}l_{1}+c_{2}l_{2}+c_{1}l_{3}+c_{2}l_{3}\pm\sqrt{\varDelta}}{2c_{1}c_{2}\left(l_{2}l_{3}+l_{1}l_{2}+l_{1}l_{3}\right)}, $

      其中

      $\begin{split}\varDelta =\,& \left[c_{1}\left(l_{1}+l_{3}\right)+c_{2}\left(l_{2}+l_{3}\right)\right]^{2}\\ &-4c_{1}c_{2}\left(l_{2}l_{3}+l_{1}l_{2}+l_{1}l_{3}\right)\end{split},$

      对称情形下$ l_{1} = l_{2} = l_{3} = l, c_{1} = c_{2} = c $

      $ \omega_{+} = \sqrt{\frac{3}{cl}},\quad \omega_{-} = \sqrt{\frac{1}{cl}}. $

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

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

      图  3  带互感的$ L\text-C$介观电路

      Figure 3.  $ L\text-C$ mesoscopic circuit with mutual inductance 2m.

      $\begin{split} L =\, & \frac{1}{2}\left[l_{1}I_{1}^{2}+l_{2}I_{2}^{2}+l_{3}\left(I_{1}\pm I_{2}\right)^{2}\pm 2m I_{1}I_{2}\right]\\ &-\frac{1}{2}\left(\frac{q_{1}^{2}}{c_{1}}+\frac{q_{2}^{2}}{c_{2}}\right),\end{split}$

      对应IEO方法中的矩阵

      $ {{M}} = \left(\!\!\!{\begin{array}{*{20}{c}} {{l_1} + {l_3}}&{{l_3} \pm m}\\ {{l_3} \pm m}&{{l_2} + {l_3}} \end{array}}\!\!\!\!\right),\quad {{N}} = \left( {\begin{array}{*{20}{c}} {\dfrac{1}{{{c_1}}}}&0\\ 0&{\dfrac{1}{{{c_2}}}} \end{array}} \right), $

      同样可得本征频率

      $ \omega^{2} = \frac{c_{1}l_{1}\!+\!c_{2}l_{2}\!+\!c_{1}l_{3}\!+\!c_{2}l_{3}\!-\!2c_{1}c_{2}l_{3}\!-\!2c_{1}c_{2}m\pm\sqrt{\varDelta}}{2c_{1}c_{2}\left(l_{2}l_{3}+l_{1}l_{2}+l_{1}l_{3}-2ml_{3}-m^{2}\right)} , $

      其中

      $ \begin{split} \varDelta =\, & \left[c_{1}\left(l_{1}+l_{3}\right)+c_{2}\left(l_{2}+l_{3}\right)-2c_{1}c_{2}l_{3}-2c_{1}c_{2}m\right]^{2}\\ &+4c_{1}c_{2}\left(c_{1}c_{2}\!-\!1\right)\left(l_{2}l_{3}\!+\!l_{1}l_{2}\!+\!l_{1}l_{3}\!-\!2ml_{3}\!-\!m^{2}\right) \end{split}. $

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

      $ \omega_{+} = \sqrt{\frac{c+1}{3c\left(l+m\right)}},\qquad \omega_{-} = \sqrt{\frac{1-c}{c\left(l-m\right)}}, $

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

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

参考文献 (22)

目录

    /

    返回文章
    返回