カーボンナノチューブナノワイヤ共振器における幾何学的位相効果の探究

概要

幾何学的位相は、振動の波動関数における追加的な位相発展であり、広範な科学技術分野への応用が期待されています。本研究では、カーボンナノチューブ（CNT）ベースのナノワイヤ共振器における圧縮状態での幾何学的位相の特性を、不変演算子法を用いて理論的に解析しました。複雑な時間依存ハミルトン系の処理に有用な線形不変演算子を導入し、幾何学的位相の解析的公式を導出しました。この公式を用いて、関連するパラメータ設定に基づく幾何学的位相の時間発展を詳細に分析しました。圧縮パラメータが幾何学的位相の発展に及ぼす影響を検討した結果、幾何学的位相は概して振動し、その包絡線は時間とともに増大することがわかりました。古典的振動振幅、減衰因子、駆動力振幅などのパラメータが大きいほど、幾何学的位相の増大率が大きくなります。特に、系の角周波数が共鳴角周波数に近づくと、幾何学的位相が時間に対して非常に鋭く増大することを確認しました。ナノワイヤ共振器における幾何学的位相の特性解明は、ナノワイヤ振動のトポロジカルな特徴を理解する上で不可欠です。

はじめに

CNTベースナノワイヤ、半導体ナノワイヤ、グラフェン、浮遊粒子など、最小スケール共振器の力学振動は10年以上にわたりナノサイエンスの主要研究テーマです。外部周期駆動力によるナノワイヤ共振器の電気力学的振動について、理論・実験両面で活発な研究が進められています。特にCNTベースナノワイヤ共振器は、周囲からの微小摂動に対する高感度と高品質因子を備え、ナノスケール力学デバイスとして注目されています。浮遊CNTナノワイヤ共振器は、電磁波、微小力、質量、温度、ノイズなど幅広い物理量測定装置として有望視されています。

ナノワイヤ振動の量子位相発展の解析は、系の基礎的特徴を理論的に解明する上で必要です。CNTナノワイヤ共振器の量子振動状態において、通常の動力学的位相に加えて幾何学的位相が補助的な位相発展として現れます。幾何学的位相は量子状態のアホロノミーであり、多様な物理学分野への応用が可能です。幾何学的位相の解析は、ナノワイヤのナノ特性（共鳴プロファイル、強い量子振動、歪み緩和機構、ディラック磁気プラズモンの出現、アハロノフ・ボーム振動のトポロジーなど）の特性評価に潜在的に採用できます。

非断熱ダイナミクスに関連する幾何学的位相の研究は、正確なシミュレーション技術の発展に必要なナノメカニカルシステムへの洞察を提供します。量子状態の制備・操作・検出は量子技術の重要な要素です。本研究の目的は、ナノワイヤ振動の量子状態における幾何学的位相の時間発展を明らかにすることです。CNTナノワイヤ振動のメカニズムを理解するため、古典的類似状態である圧縮状態における幾何学的位相の時間発展を調べます。圧縮状態の利点は、一方の四元数成分の不確定性を犠牲に他方を大幅に低減できる点にあり、これはコヒーレント状態では不可能です。特に、共鳴が幾何学的位相に及ぼす影響を分析します。共鳴エネルギーは非共鳴状態のエネルギーと大きく異なるため、波動関数のトポロジカルな振る舞いは非自明で、通常状態から大きく逸脱する可能性があります。物理パラメータと圧縮パラメータの変化が幾何学的位相の発展に及ぼす影響も厳密に解析します。幾何学的位相は力学系に普遍的に存在し、量子計算、強度干渉計、フォトニックマルチタスク、量子センシングプロトコル、波動安定性測定など現代技術への応用が期待されます。

系のハミルトニアンには、系の減衰と外部駆動力に関連する時間関数が含まれるため、時間依存ハミルトン系（TDHS）の一種です。TDHSの量子力学的問題は近年まで広く研究されてきました。TDHSのハミルトニアンに含まれる時間関数は、多くの場合、正準変数の関数から分離できないため、シュレディンガー方程式を解くための従来の変数分離法は適用できません。この困難を克服するために開発された強力な代替手法が、LewisとRiesenfeldによって導入された不変演算子法です。この手法は、TDHSの量子解を導出する際に非常に有用な数学的ツールです。カオス的粒子散乱、時間変化媒質での光伝播、駆動トラップ電子の制御、量子ナノ電子回路の非古典性など、TDHSで記述される多くの量子力学的問題がこの手法で調べられています。TDHSの量子力学的処理には、ユニタリ変換法、リー代数法、ハミルトニアン推定法など多様な手法も存在します。

本系がTDHSであるため、不変演算子法を用いて量子解を求めます。消滅演算子で表される線形不変演算子を導入します。系の時間依存性により消滅・生成演算子が時間で表されますが、これらのラダー演算子を用いてコヒーレント状態と圧縮状態の両方を得られます。圧縮状態の波動関数を用いて、系の幾何学的位相を解析的に評価します。多様なパラメータ選択に基づく描像により、幾何学的位相の時間発展を詳細に分析します。

手法

幾何学的位相を調べるには、まずナノワイヤ先端の古典的運動方程式を設定する必要があります。幾何学的位相がTDHSの量子波動発展に現れるため、扱う特定の量子状態での波動関数を導出する必要があります。前述の通り圧縮状態を考慮します。TDHSの多様な量子状態（圧縮状態を含む）の波動関数は、不変演算子法から得られます。

有効質量mを持つ浮遊カーボンナノチューブの曲げモードにおける時間依存振幅xの運動方程式は以下で与えられます[1]：

$$ \ddot{x}+\left(\frac{\omega_{0}}{Q} +\eta x^{2}\right) \dot{x}+\left(\omega_{0}^{2}+\beta x^{2}\right) x = f_{\mathrm{d}}\cos (\omega t), \quad (1) $$

ここで、ω₀は共鳴角周波数、Qは品質因子、f_dは静電駆動力をmで割った値、ηは非線形減衰係数、βはダフィングパラメータです。ナノワイヤ長に対する先端変位が十分小さいと仮定すると、式(1)の非線形項を無視でき[2]：

$$ \ddot{x}+\frac{\omega_{0}}{Q} \dot{x}+\omega_{0}^{2} x = f_{\mathrm{d}}\cos (\omega t). \quad (2) $$

式(2)を導く系のハミルトニアンは：

$$ \hat{H}= e^{-\gamma t} \frac{\hat{p}^{2}}{2m} +\frac{1}{2}me^{\gamma t} \left[\omega_{0}^{2} \hat{x}^{2} - 2f_{\mathrm{d}}\cos (\omega t)\hat{x}\right], \quad (3) $$

ここでγ = ω₀/Qです。式(2)の古典解は、余関数X_c(t)と特殊解X_p(t)で構成され：

$$ X_{c}(t)=X_{c,0}e^{-\gamma t/2}\cos(\Omega t+\varphi), \quad (4) $$

$$ X_{p}(t) =X_{p,0}\cos (\omega t - \delta), \quad (5) $$

ここでX_{c,0}は定数、Ω = √(ω₀² - γ²/4)、φは任意の位相、

$$ X_{p,0}=\frac{f_{\mathrm{d}}}{\sqrt{\left(\omega_{0}^{2} -\omega^{2}\right)^{2} + \gamma^{2} \omega^{2}}}, \quad (6) $$

$$ \delta = \tan^{-1} \frac{\gamma \omega}{ \omega_{0}^{2} -\omega^{2}}. \quad (7) $$

運動量空間での古典解も同様で、余関数はP_c(t) = m e^{γt} \dot{X}_c(t)、特殊解はP_p(t) = m e^{γt} \dot{X}_p(t)です。系の幾何学的位相を調べるには、まず量子解を導出する必要があります。式(3)のハミルトニアンが時間に明示的に依存することに注意してください。このような時間変化系の量子解を導出するため、不変演算子法[30, 31]を用います。系の不変演算子ÎはLiouville-von Neumann方程式から導出されます：

dÎ/dt = ∂Î/∂t + [Î,Ĥ]/(iℏ) = 0。

式(3)を代入して厳密に評価すると、線形不変演算子[34]が得られます：

$$ \hat{I} = \hat{A} e^{i\Omega t}, \quad (8) $$

ここでÂは消滅演算子で：

$$ \hat{A} = \left(2\hbar m\Omega\right)^{-1/2} \left[ m \left(\Omega+ i\frac{\gamma}{2} \right) e^{\gamma t/2}\left[\hat{x}-X_{p}(t)\right] +ie^{-\gamma t/2} \left[\hat{p}-P_{p}(t)\right] \right]. \quad (9) $$

式(9)のエルミート随伴†が生成演算子です。

Âの固有値方程式は：

$$ \hat{A} |A \rangle = A |A \rangle. \quad (10) $$

これを評価すると固有値は：

$$ A(t) = A(0) e^{-i\Omega t}, \quad (11) $$

ここでA(0)=A₀e^{-iφ}、A₀ = [mΩ/(2ℏ)]^{1/2}X_{c,0}です。(12)

コヒーレント状態|A⟩はÂの固有状態ですが、圧縮状態は演算子B̂の固有状態です：

$$ \hat{B} = \mu \hat{A} + \nu \hat{A}^{\dagger}, \quad (13) $$

ここでμとνは以下を満たす複素変数です：

$$ |\mu|^{2} - |\nu|^{2} =1. \quad (14) $$

B̂の固有値方程式を：

$$ \hat{B} |B \rangle = B |B \rangle, \quad (15) $$

と書くと、|B⟩は圧縮状態です。この方程式を構成空間で解くと：

$$ \langle {x}|B\rangle = \sqrt[4]{\frac{m \Omega e^{\gamma t}}{\hbar\pi(\mu-\nu)(\mu^{*}-\nu^{*})}} \exp \left\{- \frac{1}{\hbar (\mu-\nu)} \left[\frac{1}{2} m e^{\gamma t}\left((\mu+\nu)\Omega +\frac{i\gamma}{2}(\mu-\nu)\right)\left[x-X_{p}(t)\right]^{2} - \left[iP_{p}(t)(\mu-\nu)+ \left(2\hbar m \Omega e^{\gamma t}\right)^{1/2} (\mu A+\nu A^{*}) \right]\left[x-X_{p}(t)\right] -\frac{|A|^{2}+A^{2}}{2(\mu-\nu)(\mu^{*}-\nu^{*})} \right\}. \quad (16) $$

以上で圧縮状態での波動関数が式(16)として導出されました。この系の量子的特徴は、この波動関数の解析的記述に基づいて明確化できます。μ=1、ν=0とすると式(16)はコヒーレント状態の波動関数（式(10)の構成空間での固有状態）に帰着します。式(16)の波動関数を次節で用いて、圧縮状態での幾何学的位相を導出します。

結果と考察

量子波動発展の位相には、動力学的位相に加えて幾何学的位相が含まれることはよく知られています。幾何学的位相は、Berryによって1984年に[12]、パラメータ空間での循環的・断熱的進化をする系で初めて発見されました。量子力学の断熱定理によれば、パラメータ空間での循環的進化における量子状態の瞬間固有状態は後にも同じ状態に留まりますが、追加の量子位相蓄積、すなわちBerry位相が生じます。Berry位相を、非断熱・非循環・非ユニタリな量子系の進化も含む形で一般化したものが幾何学的位相です。

圧縮状態での幾何学的位相は：

$$ \gamma_{G}(t) = \int_{0}^{t} \langle B(t') |i\frac{\partial}{\partial t'}| B(t') \rangle dt' +\gamma_{G}(0). \quad (17) $$

構成空間での波動関数の時間微分は：

$$ \frac{\partial \langle {x}|B\rangle}{\partial t} = \left\{ f_{1}(t) \left[x-X_{p}(t)\right]^{2} + f_{2}(t) \left[x-X_{p}(t)\right] + f_{3}(t) \right\} \langle {x}|B\rangle, \quad (18) $$

ここで、

$$ f_{1}(t) = - \frac{m\gamma e^{\gamma t}}{2\hbar (\mu-\nu)} \left((\mu+\nu)\Omega + \frac{i\gamma}{2}(\mu-\nu) \right), \quad (19) $$

$$ f_{2}(t) = \frac{1}{\hbar (\mu-\nu)}\left[ \left((\mu+\nu)\Omega + \frac{i\gamma}{2}(\mu-\nu) \right) P_{p}(t) -i m e^{\gamma t} \left[\omega_{0}^{2} X_{p}(t) - f_{\mathrm{d}} \cos(\omega t)\right](\mu-\nu) +\left(2\hbar m \Omega e^{\gamma t}\right)^{1/2} \left(\frac{\gamma}{2}\left(\mu A + \nu A^{*}\right)-i\Omega \left(\mu A - \nu A^{*}\right) \right) \right], \quad (20) $$

$$ f_{3}(t) = \frac{\gamma}{4}-\frac{1}{\hbar m e^{\gamma t}(\mu-\nu)} \left[iP_{p}(t)(\mu-\nu) + \left(2\hbar m\Omega e^{\gamma t}\right)^{1/2} \left(\mu A+\nu A^{*}\right)\right] P_{p}(t)+ \frac{i\Omega A^{2}}{(\mu-\nu)\left(\mu^{*}-\nu^{*}\right)}. \quad (21) $$

式(18)を式(17)に代入してさらに評価すると：

$$ \gamma_{G}(t) = \int_{0}^{t} dt' \left[ A_{0}^{2}\left(\frac{\gamma^{2}}{4\Omega}+\Omega + g_{1} \sin\left[2\left(\Omega t'+\varphi\right)\right] +g_{2} \cos\left[2\left(\Omega t'+\varphi\right)\right] \right) -A_{0}\left[ g_{3}(t') \sin\left(\Omega t'+\varphi\right) +g_{4}(t') \cos\left(\Omega t'+\varphi\right) \right]+ g_{5}(t') \right] +\gamma_{G}(0), \quad (22) $$

ここで、

$$ g_{1} = \frac{\gamma}{2} + \frac{i\Omega \left(\mu\nu^{*}-\mu^{*}\nu\right)}{(\mu-\nu)\left(\mu^{*}-\nu^{*}\right)}, \quad (23) $$

$$ g_{2} = \frac{\gamma^{2}}{4\Omega}+\Omega\frac{2|\nu|^{2}- \left(\mu\nu^{*}+\mu^{*}\nu\right)}{(\mu-\nu) \left(\mu^{*}-\nu^{*}\right)}, \quad (24) $$

$$ g_{3}(t) = \left(\frac{2\Omega}{m\hbar e^{\gamma t}} \right)^{1/2}P_{p}(t), \quad (25) $$

$$ g_{4}(t) = \frac{1}{\sqrt{2\hbar\Omega}}\left(\frac{\gamma }{\sqrt{m e^{\gamma t}}}P_{p}(t) - 2\sqrt{m e^{\gamma t}}\left[\omega_{0}^{2} X_{p}(t) - f_{\mathrm{d}} \cos(\omega t)\right]\right), \quad (26) $$

$$ g_{5}(t) = \frac{P_{p}^{2}(t)}{\hbar m e^{\gamma t}}+\frac{\gamma^{2}}{8\Omega}\left[2|\nu|^{2}-\left(\mu\nu+\mu^{*}\nu^{*}\right) +1\right] +\frac{i\gamma}{4(\mu-\nu)\left(\mu^{*}-\nu^{*}\right)} \left[|\mu|^{2}\left(\nu^{2}-\nu^{*2}\right)-|\nu|^{2}\left(\mu^{2}-\mu^{*2}\right)\right. \left.+ (2|\nu|^{2}+1)\left(\mu\nu^{*}-\mu^{*}\nu\right) +(\mu-\mu^{*})(\nu-\nu^{*})\right]. \quad (27) $$

g₅の最後の項（μ−μ*）（ν−ν*）を含む部分は純虚数であり、位相として不適切です。したがって、μとνの少なくとも一方を実数値として選ぶことでこの項を除去します。この処置は、μとνの絶対位相ではなく相対位相のみが物理的意味を持つため、一般性を損なわずに常に可能です。

式(22)の積分を実行すると：

$$ \gamma_{G}(t) = A_{0}^{2}\left[\left(\frac{\gamma^{2}}{4\Omega}+\Omega\right)t + \frac{g_{1}}{\Omega}\sin(\Omega t+2\varphi) \sin(\Omega t) +\frac{g_{2}}{\Omega} \cos(\Omega t+2\varphi) \sin(\Omega t) \right] -A_{0}\left[ \left(\frac{2m\Omega}{\hbar} \right)^{1/2}\omega X_{p,0} \bar{g}_{3}(t) +\sqrt{\frac{2m}{\hbar\Omega}}\frac{1}{4\omega^{2}+\gamma^{2}}\bar{g}_{4}(t) \right] + \bar{g}_{5}(t) +\gamma_{G}(0), \quad (28) $$

ここで、\(\bar {g}_{i}(t)~(i=3,4,5)\)は：

$$ \bar{g}_{i}(t) = G_{i}(t) -G_{i}(0), \quad (29) $$

と定義され、

$$ G_{3}(\tau) = e^{\gamma \tau/2}\left(\frac{1}{4(\Omega+\omega)^{2}+\gamma^{2}} \left\{2(\Omega+\omega)\sin[(\Omega+\omega)\tau+\varphi-\delta] +\gamma \cos[(\Omega+\omega)\tau+\varphi-\delta] \right\} - \frac{1}{4(\Omega-\omega)^{2}+\gamma^{2}} \left\{ 2(\Omega-\omega) \sin[(\Omega-\omega)\tau+\varphi+\delta]+\gamma \cos[(\Omega-\omega)\tau+\varphi+\delta]\right\} \right), \quad (30) $$

$$ G_{4}(\tau) = e^{\gamma \tau/2} \left\{X_{p,0} \left\{ \gamma\omega[ 2\omega\cos(\omega \tau-\delta)-\gamma \sin(\omega \tau-\delta)] -2\omega_{0}^{2} \left[2\omega\sin(\omega \tau-\delta)+\gamma \cos(\omega \tau-\delta)\right] \right\}+2f_{\mathrm{d}} [ 2\omega\sin(\omega \tau) +\gamma \cos(\omega \tau)] \right\}, \quad (31) $$

$$ G_{5}(\tau) = \frac{m\omega^{2}}{2\hbar}X_{p,0}^{2} \frac{e^{\gamma \tau}}{\gamma \left(4\omega^{2}+\gamma^{2}\right)} \left\{ \gamma^{2}+4\omega^{2} -\gamma^{2} \cos[2(\omega\tau -\delta)] -2\gamma\omega \sin[2(\omega \tau -\delta)] \right\} +\frac{\gamma^{2} \tau}{8\Omega}\left[2|\nu|^{2}-\left(\mu\nu+\mu^{*}\nu^{*}\right)+1\right] +\frac{i\gamma \tau}{4(\mu-\nu)\left(\mu^{*}-\nu^{*}\right)} \left[|\mu|^{2}\left(\nu^{2}-\nu^{*2}\right)-|\nu|^{2}\left(\mu^{2}-\mu^{*2}\right)\right. \left.+\left(2|\nu|^{2}+1\right)\left(\mu\nu^{*}-\mu^{*}\nu\right)\right]. \quad (32) $$

以上で、圧縮状態における完全な幾何学的位相が式(28)および式(23)、(24)、(29)〜(32)により得られました。

幾何学的位相の時間発展を図1〜4に示します。図1から、幾何学的位相は振動し、その包絡線が時間とともに増大することがわかります。包絡線の増大はA₀が大きいほど顕著です。振動パターンはμとνの値が増加するにつれて徐々に不規則になり、振動振幅も時間経過とともに大きくなります。

図2は減衰因子γの影響を示します。ω=0.3（a）、0.99（b）、5（c）の場合で、圧縮パラメータμ=√2、ν=1（q圧縮状態）としています。共鳴に近い場合（b）では幾何学的位相が時間に対して非常に急速に増大します。図3はμ=√2、ν=−1（p圧縮状態）の場合で、図2と全体的に大きな差はなく、μとνの絶対値が変わらない限り圧縮の種類（qまたはp）によらず幾何学的位相の発展はほぼ無関係であることが確認されます。

図2a、3aと図2c、3cを比較すると、ωが共鳴角周波数より大きい場合に幾何学的位相がやや急速に変化することがわかります。系の共鳴状態における幾何学的位相の時間振る舞いは大きな関心事です[22, 23]。図2b、3bは、ωが共鳴角周波数に近いときに幾何学的位相が非常に急速に増大することを示しています。これは、幾何学的位相の大きさが波動関数の時間変化に関連するため、この状況では波動関数が時間とともに大きく変化することを意味します。実際、ワイヤ振動の振幅は共鳴状態で著しく増大します。浮遊CNTナノワイヤ共振器の共鳴角周波数は高いだけでなく、非常に高い品質因子で広範囲にチューニング可能です[3]。したがって、系の振動モードは完全に減衰するまで長時間維持されます[11]。

図4は駆動力振幅f_dの影響を示します。f_dが増加すると、幾何学的位相の時間増大率が急速になります。

結論

シュレディンガー方程式に基づく量子力学の枠組みで、CNTベースナノワイヤ共振器における圧縮状態の幾何学的位相を調査しました。系を記述するハミルトニアンの時間依存性に対し、時間で記述されるハミルトニアンの場合に量子解を導出する有力なツールである不変演算子法を導入しました。この手法により、CNTナノワイヤ振動における幾何学的位相の解析的公式を取得しました。

力学振動の理論的理解に必要な位相効果の詳細な解析を行いました。本研究での幾何学的位相の展開は、厳密な数学的評価に基づく完全な量子論的アプローチです。幾何学的位相は機械的パラメータの変化に敏感で、多くの場合振動を示します。圧縮パラメータが発展に及ぼす影響も分析しました。共鳴角周波数付近で幾何学的位相の蓄積が時間とともに強く増大することを確認しました。

本成果は、CNTナノワイヤ振動に現れる幾何学的位相の時間振る舞いを示しています。本研究で与えられた幾何学的位相の解析は、系のトポロジカルな特徴だけでなく、他のナノワイヤベース力学振動子の動力学的振動を理解する上でも重要です。特に、量子情報技術やその他の量子ベース産業へのシステム応用に必要な共鳴状態の位相特性を明らかにしました[40]。本研究で用いた類似の手法とフレームワークは、超伝導ファブリ・ペロー共振器[41]、ナノカンチレバー[42]、量子ビット-共振器-原子ハイブリッドシステム[43]など、他のナノシステムにも拡張可能です。