Low-energy model for doped graphene nanoribbons

技術要約：ドープされたグラフェンナノリボンに対する低エネルギーモデル

問題提起
本論文は、化学ポテンシャルがバルクの単一粒子バンド内に存在する、自由空間におけるドープされたグラフェンナノリボン（GNR）の多体電子挙動を扱っている。GNRはトランジスタや量子ドットのようなナノ電子デバイスの有望な候補であるが、境界条件によって$\pm k$波数が混合するため、バレー量子数が無効となり、その挙動は複雑化する。従来の理論研究は、主に密度汎関数理論（DFT）または平均場ハバードモデルに依存してきた。しかし、DFTは一体理論であるため、クーロンブロック現象や近藤物理学といった相関電子現象を正しく記述できず、しばしばクーロン・ダイヤモンドの数や位置、および励起状態の線に関する誤った予測を与える。一方、実空間における標準的なハバードモデルは、輸送実験に関連する狭いエネルギー窓の多体状態を厳密に解くには計算量的に困難である。著者らは、ドープされたGNRにおける運動エネルギーと電子間相互作用の競合を捉えるための低エネルギー有効モデルを開発することで、このギャップを埋めることを目的としている。具体的には、開殻（open-shell）および高スピン状態の出現を理解することを目指している。

手法
著者らは、以下の多段階の理論的アプローチを用いている：

1. 厳密なマッピング： 実空間のサイト基底で定義された拡張ハバードモデル（オンサイト相互作用$U^0$およびサイト間相互作用$U^n$を含む）から出発する。彼らは、このモデルをリボンの単一粒子バルク固有状態の基底で書かれたカナモリ・ハミルトニアンへと厳密にマッピングする。この変換は、有限長のリボンに対するタイトバインディングモデルの解析解を利用している。
2. 解析的導出： 著者らは、マッピングの結果として生じるクーロン行列要素（直接積分および交換積分）の解析的な式を導出する。そして、これらの要素の、リボン幅（$M_y$）および長さ（$M_x$）の関数としてのスケーリング挙動を分析する。
3. パラメータ決定： モデルを簡略化するため、拡張ハバードパラメータを繰り込まれたオンサイト相互作用$U^0$へと集約する。著者らは、様々な幅（$N=7, 9, 17$）および長さのアーミチェア型リボンに対するDFT計算に、オンサイト・ハバードモデルの平均場解をフィッティングさせることで、$U^0$の値を決定する。その結果、$U^0 \approx 4.9$ eVがDFTのギャップおよび励起エネルギーに対して良好な一致を示すことを見出した。
4. 低エネルギー・ハミルトニアンの構築： 著者らは、化学ポテンシャル付近の少数の「開殻（open shells）」状態のみを保持し、より深い「閉殻（closed shells）」を不活性として扱う、低エネルギー多体ハミルトニアンを構築する。このハミルトニアンには、オンサイト反発、直接クーロン相互作用、強磁性的交換相互作用、およびペア・ホッピング項が含まれる。
5. 厳密対角化： 低エネルギー・ハミルトニアンを（2〜3個のシェルについては解析的に、より大きな系については数値的に）厳密に解き、グランド状態の特性、特に開殻および高スピン構成となる条件を決定する。
6. 輸送モデリング： 転送ハミルトニアンを導入して、逐次電子トンネリング（SET）領域における遷移速度を計算し、導出された多体状態が電気輸送にどのように影響するかを分析する。

主要な貢献と結果

• カナモリ・マッピング： 本論文は、拡張ハバードモデルからGNR用のカナモリ・モデルへの厳密なマッピングを確立している。これにより、固有状態基底におけるハバード相互作用が、直接クーロン積分だけでなく、強磁性的交換相互作用やペア・ホッピング相互作用をも生成することが明らかになった。
• スケーリング則： 著者らは、クーロン行列要素がリボンのサイズに応じて3つの異なるタイプ、すなわち同一状態積分（$U$）、同一バンド内積分（$V$）、およびバンド間積分（$W$）に集約されることを決定した。彼らは、$V \approx \frac{2}{3}U$ および $W \approx \frac{2}{3}V$ という極めて高い精度でのスケーリング則を見出した。さらに、リボンの次元の関数としての$U$に関する解析的なサイズ・スケーリング関係も提供している。
• 高スピン・グランド状態： 低エネルギー・ハミルトニアンの解析により、強磁性的交換相互作用がドープされたGNRにおける開殻・高スピン多体状態を安定化させ得ることが示された。
  • 2シェルモデルの場合、有効交換定数が正であるとき（これは、シェル間のエネルギー差が相互作用強度に対して小さいときに起こる）、グランド状態は高スピン・トリプレット（$S=1$）となる。
  • 3シェルモデルの場合、高スピン・クァルテット（$S=3/2$）が存在するための条件が導出されている。
• 定量的推定： 様々なサイズの リボンにこれらの条件を適用することで、著者らはこのような高スピン状態が見出される確率を推定している。自由空間のリボンの場合、シェル・ペアの約29%、およびシェル・トリプル（三つ組）の約17%が、$S=1$および$S=3/2$の状態を満たす条件を満たしている。これらの確率はリボンのサイズとともに増加するが、飽和する。
• 遮蔽効果： 本研究では、基板効果をシミュレートするために有効誘電率（$\epsilon$）を組み込んでいる。彼らは、誘電体遮蔽が有効なクーロン相互作用を減少させ、その結果、開殻・高スピン挙動を示すリボンの割合を減少させることを発見した。
• 輸送ブロック現象： モデルは、輸送実験における「シェル・ブロック（shell-blockade）」および「スピン・ブロック（spin-blockade）」のメカニズムを予測している。シェル・ブロックは、2つの単一粒子状態の同時変化を必要とする遷移が禁止されることから生じる。スピン・ブロックは、単一電子トンネリングが、互換性のない全スピンを持つ状態間（例：シングレットから高スピン状態への遷移）のギャップを埋めることができないことから生じる。

意義
本論文は、提案された低エネルギー・カナモリ・ハミルトニアンが、GNRの固有状態の狭いエネルギー窓に対処するための、クリーンかつ解析的に扱いやすい枠組みを提供すると主張している。これは、単一または数個のリボン接合に関する近年の輸送実験に直接関連している。DFTとは異なり、このアプローチは多体相関を正しく捉えており、クーロン・ブロックや安定性図における励起状態の線に関する理論的基礎を提供する。高スピン・開殻状態の堅牢なメカニズムの特定は、GNRがスピン・ブロックのような複雑な磁気的および輸送現象を示し得ることを示唆しており、これらはナノ電子応用において利用できる可能性がある。本研究は、DFTが欠いている、特に多体準粒子ポールのエネルギー位置と数に関する必要な多体補正を提供することで、DFTを補完するものである。