カーボンナノチューブ・グラフェン・量子ドット - 高性能化の設計原理
この章を完了すると、以下を説明できるようになります：
ナノ材料（Nanomaterials）とは、少なくとも1つの次元が概ね1〜100 nmの範囲にある材料を指します。この寸法領域では、バルク材料とは質的に異なる物性が現れます。その源泉は主に次の2つです：
直径10 nmの金ナノ粒子では全原子の約16%が表面に露出しています（本章3.5節で計算します）。バルクの金（表面原子割合はほぼ0%）が化学的に不活性なのに対し、金ナノ粒子は室温でCO酸化を触媒します。表面原子は配位不飽和で反応性が高く、この違いがナノ材料の触媒・光学・電子物性を生み出します。
ナノ材料は「ナノスケールに閉じ込められている方向の数」によって分類されます。閉じ込め方向の数だけ電子の自由度が失われ、状態密度（Density of States, DOS）の形が変わります。
flowchart LR
A[ナノ材料の次元分類] --> B[0次元 0D
量子ドット]
A --> C[1次元 1D
ナノチューブ・ナノワイヤ]
A --> D[2次元 2D
グラフェン・ナノシート]
B --> B1[3方向を閉じ込め
離散準位]
C --> C1[2方向を閉じ込め
1方向に伝導]
D --> D1[1方向を閉じ込め
面内に伝導]
style A fill:#f093fb
style B fill:#e3f2fd
style C fill:#e8f5e9
style D fill:#fff3e0
| 次元 | 閉じ込め方向 | 代表例 | 状態密度の特徴 |
|---|---|---|---|
| 0D（0次元） | 3方向すべて | 量子ドット、フラーレン | 離散的なデルタ関数状 |
| 1D（1次元） | 2方向 | カーボンナノチューブ、ナノワイヤ | 1/√Eの発散（van Hove特異点） |
| 2D（2次元） | 1方向 | グラフェン、遷移金属ダイカルコゲナイド | 階段関数状（一定） |
| 3D（バルク） | なし | 結晶・多結晶固体 | √Eに比例 |
半径 r の球について、表面積 S = 4πr² と体積 V = (4/3)πr³ の比は次式で与えられます：
すなわち比表面積は半径に反比例して増大します。粒径を1/10にすると単位体積あたりの表面積は10倍になり、表面に露出する原子の割合が飛躍的に高まります。この幾何学的事実が、ナノ材料の高い触媒活性・溶解性・反応性の根拠です。定量的な計算は3.5節のPython実践で扱います。
カーボンナノチューブ（Carbon Nanotube, CNT）は、グラフェンシートを円筒状に丸めた1次元構造です。どの方向に丸めるかを表すのがカイラルベクトル（chiral vector）で、グラフェンの基本並進ベクトル a₁, a₂ を用いて整数対 (n, m) で指定されます：
この (n, m) がCNTの幾何と電子物性をすべて決めます。丸め方によって3つの型に分類されます：
直径 d とカイラル角 θ は (n, m) から解析的に求まります：
CNTの最も重要な性質は、幾何構造だけで金属になるか半導体になるかが決まることです。判定則は次のように表せます：
この規則はグラフェンのバンド構造（3.3節）から導かれます。グラフェンはK点で価電子帯と伝導帯が接する半金属で、CNTでは円周方向の波数が量子化されるため、許される波数線がK点を通るかどうかで金属性が決まります。統計的には、ランダムに合成したCNTの約1/3が金属的、2/3が半導体的になります（3.4節のPython実践で確認します）。
半導体型CNTのバンドギャップ E_g は直径 d にほぼ反比例し、E_g ≈ 0.8 eV / d[nm] 程度になります。直径1 nmのCNTなら約0.8 eVで、シリコン（1.1 eV）に近く、トランジスタチャネル材料として有望です。直径を選ぶことでバンドギャップを設計できる点がCNTの魅力です。
CNTは強いsp²炭素-炭素結合に由来する卓越した物性を示します：
工業的に主流の合成法は化学気相成長（Chemical Vapor Deposition, CVD）です。触媒金属（Fe、Co、Ni）のナノ粒子上に炭化水素ガス（メタン、エチレンなど）を供給し、600〜1000°Cで分解・析出させてCNTを成長させます。
CVDでは直径はある程度制御できますが、カイラリティ（(n,m)）を単一に揃えることは依然として困難です。半導体型のみを高純度で得るには、合成後の分離精製（密度勾配超遠心分離、ゲルクロマトグラフィーなど）が必要になります。単一カイラリティCNTの直接合成は現在も活発な研究テーマです。
グラフェン（Graphene）は、炭素原子がsp²結合で蜂の巣格子（honeycomb lattice）を組んだ、厚さ原子1個分の2次元材料です。2004年にGeimとNovoselovが粘着テープによる機械的剥離で単離し、2010年のノーベル物理学賞につながりました。
蜂の巣格子は2つの炭素原子（A副格子・B副格子）を単位胞に持ちます。強束縛近似（tight-binding）で解くと、価電子帯と伝導帯がブリルアンゾーンの6つのK点（およびK'点）で接触することがわかります。
グラフェンの決定的な特徴は、K点近傍でエネルギーが波数に対して線形に分散することです：
この接触点をディラック点（Dirac point）と呼びます。通常の半導体では E ∝ k² の放物線分散ですが、グラフェンでは線形分散のため、電子はあたかも質量ゼロの相対論的粒子（ディラックフェルミオン）のように振る舞います。フェルミ速度 v_F は光速の約1/300です。
グラフェンは価電子帯と伝導帯が1点で接するため、バンドギャップが厳密にゼロの半金属（semimetal）です。これは高移動度をもたらす一方、トランジスタとしてはオフ状態を作れない弱点になります。ナノリボン化・二層化・基板との相互作用などでギャップを開く研究が進められています。
| 製造法 | 品質 | 面積・量 | 主な用途 |
|---|---|---|---|
| 機械的剥離（スコッチテープ法） | 最高品質・欠陥極小 | 微小片のみ | 基礎研究・物性測定 |
| CVD（Cu箔上成長） | 高品質・大面積 | ロール・ツー・ロール可 | 透明電極・フレキシブル素子 |
| SiC熱分解（エピタキシャル成長） | 高品質 | ウェハスケール | 高周波電子デバイス |
| 酸化グラフェン還元（rGO） | 欠陥多め・低コスト | 大量生産・溶液プロセス | 導電インク・複合材・電極 |
量子ドット（Quantum Dot, QD）は、直径がおよそ2〜10 nmの半導体ナノ結晶です。電子と正孔（ホール）の運動が3方向すべてにわたって結晶サイズに閉じ込められるため、エネルギー準位が離散化し、バンドギャップが粒径に依存して変化します。この現象を量子閉じ込め効果（Quantum Confinement Effect）と呼びます。
閉じ込めが顕著になるのは、結晶半径が励起子ボーア半径（exciton Bohr radius）と同程度以下になったときです。CdSeでは励起子ボーア半径が約5.6 nmで、粒径がこれより小さくなるとバンドギャップが顕著に拡大します。
最も基本的な描像は量子力学の「箱の中の粒子（particle in a box）」です。1辺 L の3次元無限井戸に閉じ込められた粒子の基底エネルギーは E ∝ 1/L² で増大します。すなわち小さい箱ほど基底準位が高くなるため、粒径が小さいほどバンドギャップが広がります。
この直感を半導体ナノ結晶に定量化したのがBrus方程式（Brus equation）です。半径 R の球状ナノ結晶の実効バンドギャップは次式で表されます：
右辺は3項からなります：
粒径が小さくなると第2項（1/R²）が第3項（1/R）より速く増大するため、正味ではバンドギャップが拡大し、発光波長が短波長側（青色側）へシフトします。この関係を3.5節でPythonで計算します。
バンドギャップが粒径で決まるため、量子ドットは同一組成のまま粒径を変えるだけで発光色を連続的に調整できるという際立った特徴を持ちます。CdSe量子ドットでは、粒径を約2 nmから6 nmへ変えるだけで青から赤まで発光色が変化します。
高性能な量子ドットの多くはCdSeやCdTeなどカドミウムを含み、毒性と環境規制（RoHS等）が課題です。代替としてInP系やカーボンドット、ペロブスカイト量子ドットなどのカドミウムフリー材料が研究されていますが、発光効率や安定性でCd系に及ばない場合があり、材料開発が続いています。
本節では、これまで学んだ3つのテーマを実際に計算します。すべてのコードはNumPyのみで自己完結して動作し、実行結果を併記しています。
(n, m) 指数から直径・カイラル角・幾何型・金属/半導体を判定します。判定則 (n − m) mod 3 = 0 が金属、そして 0 ≤ m ≤ n ≤ 20 の全列挙で金属比率が約1/3になることを確認します。
# ===================================
# Example 1: CNTのカイラリティ分類
# ===================================
import numpy as np
a = 0.246 # グラフェンの格子定数 [nm]
def cnt_diameter(n, m):
"""カイラル指数 (n, m) からCNT直径を計算 [nm]"""
return a * np.sqrt(n**2 + n*m + m**2) / np.pi
def chiral_angle(n, m):
"""カイラル角を計算 [degree]"""
return np.degrees(np.arctan(np.sqrt(3)*m / (2*n + m)))
def cnt_type(n, m):
"""金属性・半導体性を判定"""
if (n - m) % 3 == 0:
return "metallic"
return "semiconducting"
def cnt_class(n, m):
"""幾何学的分類（アームチェア/ジグザグ/カイラル）"""
if m == 0:
return "zigzag"
if n == m:
return "armchair"
return "chiral"
# 代表的な (n, m) の分類
examples = [(5,5), (9,0), (10,0), (7,3), (6,4), (10,10), (8,4), (11,7)]
print("(n, m) type geometry d [nm] theta [deg]")
print("-" * 58)
for n, m in examples:
print(f"({n:2d},{m:2d}) {cnt_type(n,m):14s} {cnt_class(n,m):9s} "
f"{cnt_diameter(n,m):5.3f} {chiral_angle(n,m):5.2f}")
# 0 <= m <= n <= 20 の全 (n, m) で金属比率を集計
total = 0
metal = 0
for n in range(1, 21):
for m in range(0, n + 1):
total += 1
if (n - m) % 3 == 0:
metal += 1
print()
print(f"Enumerated (n,m) with 0<=m<=n<=20: {total}")
print(f"Metallic count: {metal} ({100*metal/total:.1f}%)")
# 実行結果:
# (n, m) type geometry d [nm] theta [deg]
# ----------------------------------------------------------
# ( 5, 5) metallic armchair 0.678 30.00
# ( 9, 0) metallic zigzag 0.705 0.00
# (10, 0) semiconducting zigzag 0.783 0.00
# ( 7, 3) semiconducting chiral 0.696 17.00
# ( 6, 4) semiconducting chiral 0.683 23.41
# (10,10) metallic armchair 1.356 30.00
# ( 8, 4) semiconducting chiral 0.829 19.11
# (11, 7) semiconducting chiral 1.231 22.69
#
# Enumerated (n,m) with 0<=m<=n<=20: 230
# Metallic count: 83 (36.1%)
アームチェア型 (5,5)・(10,10) は常に金属的（n − m = 0）で、直径は (10,10) が (5,5) の2倍になっています。全列挙では金属が230個中83個（36.1%）で、理論値の「約1/3が金属」とよく一致します。
CdSe量子ドットについて、半径を1.5〜4.0 nmまで変化させてバンドギャップと発光波長を計算します。粒径が小さいほどギャップが広く、発光が短波長（青側）になることを確認します。
# ===================================
# Example 2: Brus方程式（量子ドット）
# ===================================
import numpy as np
# 物理定数（SI単位）
h = 6.62607015e-34 # プランク定数 [J s]
me0 = 9.1093837015e-31 # 電子の静止質量 [kg]
e = 1.602176634e-19 # 電気素量 [C]
eps0 = 8.8541878128e-12 # 真空の誘電率 [F/m]
c = 2.99792458e8 # 光速 [m/s]
# CdSeのパラメータ
Eg_bulk = 1.74 # バルクバンドギャップ [eV]
m_e = 0.13 * me0 # 電子の有効質量
m_h = 0.45 * me0 # 正孔の有効質量
eps_r = 10.6 # 比誘電率
def brus_gap_eV(R_nm):
"""Brus方程式で実効バンドギャップを計算 [eV]"""
R = R_nm * 1e-9
confinement = (h**2 / (8 * R**2)) * (1/m_e + 1/m_h) / e # 量子閉じ込め項 [eV]
coulomb = 1.8 * e / (4 * np.pi * eps0 * eps_r * R) # クーロン項 [eV]
return Eg_bulk + confinement - coulomb
def emission_nm(Eg_eV):
"""バンドギャップから発光波長を計算 [nm]"""
return h * c / (Eg_eV * e) * 1e9
print("CdSe quantum dot: Brus equation")
print(f"{'radius [nm]':>11} {'diameter [nm]':>13} {'E_gap [eV]':>11} {'lambda [nm]':>12}")
print("-" * 50)
for R in [1.5, 2.0, 2.5, 3.0, 3.5, 4.0]:
Eg = brus_gap_eV(R)
lam = emission_nm(Eg)
print(f"{R:11.1f} {2*R:13.1f} {Eg:11.3f} {lam:12.1f}")
print()
print(f"Bulk CdSe gap {Eg_bulk} eV -> lambda {emission_nm(Eg_bulk):.1f} nm")
# 実行結果:
# CdSe quantum dot: Brus equation
# radius [nm] diameter [nm] E_gap [eV] lambda [nm]
# --------------------------------------------------
# 1.5 3.0 3.234 383.4
# 2.0 4.0 2.550 486.3
# 2.5 5.0 2.239 553.8
# 3.0 6.0 2.073 598.2
# 3.5 7.0 1.974 627.9
# 4.0 8.0 1.912 648.5
#
# Bulk CdSe gap 1.74 eV -> lambda 712.6 nm
半径1.5 nmでは発光波長383 nm（紫）、4.0 nmでは649 nm（赤）と、粒径だけで可視域全体をカバーできます。バルクCdSe（712 nm、近赤外）と比べ、ナノ結晶化によりバンドギャップが大きく拡大していることがわかります。この単調な粒径-発光関係がディスプレイの色設計の基盤です。
Brus方程式は有効質量近似に基づくため、半径が1 nmを下回る強閉じ込め領域ではバンドギャップを過大評価する傾向があります。上の計算で半径1.5 nm以上から始めているのはこのためです。より小さな粒子では強束縛法や第一原理計算が必要になります。
球状ナノ粒子について、表面から原子1層分の殻に含まれる原子を「表面原子」とみなし、その割合を粒径の関数として計算します。粒径が小さいほど表面原子割合が急増することを確認します。
# ===================================
# Example 3: 表面原子割合の計算
# ===================================
import numpy as np
# 殻モデル: 表面から原子1個分（原子直径）の厚さの殻を「表面」とみなす
# F_surface = 1 - ((R - t)/R)^3, t = 表面殻の厚さ [nm]
r_atom = 0.144 # 金の金属半径 [nm]
t = 2 * r_atom # 表面殻の厚さ（原子直径1個分）[nm]
def surface_fraction(D_nm):
"""直径 D の球ナノ粒子の表面原子割合"""
R = D_nm / 2.0
if R <= t:
return 1.0
return 1.0 - ((R - t) / R)**3
def n_total(D_nm):
"""おおよその全原子数（充填率0.74のfccを仮定）"""
R = D_nm / 2.0
return 0.74 * (R / r_atom)**3
print("Gold nanoparticle: surface-atom fraction (r_atom = 0.144 nm)")
print(f"{'diameter [nm]':>13} {'~N_total':>10} {'surface fraction':>17}")
print("-" * 44)
for D in [1, 2, 5, 10, 20, 50, 100]:
print(f"{D:13d} {n_total(D):10.0f} {surface_fraction(D)*100:16.1f}%")
# 実行結果:
# Gold nanoparticle: surface-atom fraction (r_atom = 0.144 nm)
# diameter [nm] ~N_total surface fraction
# --------------------------------------------
# 1 31 92.4%
# 2 248 63.9%
# 5 3872 30.7%
# 10 30978 16.3%
# 20 247825 8.4%
# 50 3872258 3.4%
# 100 30978063 1.7%
直径1 nmでは全原子の9割以上が表面に露出しているのに対し、100 nmではわずか1.7%です。この表面原子割合の急増が、金ナノ粒子の触媒活性・融点降下・プラズモン発色といったサイズ依存物性の根本にあります。3.1.3節の S/V = 3/r の関係が、原子数レベルでも同様に効いていることがわかります。
この章を完了すると、以下を説明できるようになります：
カイラル指数 (12, 6) のカーボンナノチューブは金属的でしょうか、半導体的でしょうか。判定則を用いて答えてください。
正解: 金属的
解説:
判定則は「(n − m) が3の倍数なら金属的」です。
n − m = 12 − 6 = 6 = 3 × 2 なので3の倍数です。したがって (12, 6) は金属的と判定されます。
なお n = 2m の場合でも、判定に効くのは n − m の値だけである点に注意してください。
次のナノ材料を 0D・1D・2D に分類してください： (a) グラフェン、(b) CdSe量子ドット、(c) 単層カーボンナノチューブ。
正解:
閉じ込め方向の数（3・1・2）に対応して、自由に動ける方向の数が 0・2・1 になります。
同じCdSeで作った2つの量子ドットA（半径2.0 nm）とB（半径3.5 nm）があります。より短波長（青色側）で発光するのはどちらですか。理由も述べてください。
正解: A（半径2.0 nm）
解説:
量子閉じ込め項は 1/R² で効くため、半径が小さいほどバンドギャップが広くなります。バンドギャップが広いほど発光光子のエネルギーが高く、波長は短くなります。
本文3.5.2の計算では、半径2.0 nmで発光波長486 nm、半径3.5 nmで628 nmでした。したがって半径の小さいAの方が短波長（青緑）、Bの方が長波長（赤）で発光します。
アームチェア型 (10, 10) CNTの直径を計算してください。格子定数 a = 0.246 nm を用います。式 d = (a/π)·√(n² + nm + m²) を使ってください。
解答:
n = m = 10 を代入します。
n² + nm + m² = 100 + 100 + 100 = 300
√300 ≈ 17.32
d = (0.246 / π) × 17.32 = 0.0783 × 17.32 ≈ 1.356 nm
本文3.5.1のコード出力（(10,10) → 1.356 nm）と一致します。参考までに (5,5) では √75 ≈ 8.66 となり d ≈ 0.678 nm で、ちょうど半分の直径になります。
直径10 nmの金ナノ粒子では表面原子割合が約16%、直径2 nmでは約64%でした。触媒として使う場合、なぜ小さい粒子ほど有利になりやすいのか、表面効果の観点から説明してください。また小さくしすぎる欠点も1つ挙げてください。
解答例:
小さい粒子が触媒に有利な理由:
触媒反応は表面に露出した原子上で進行します。表面原子割合が高いほど、同じ質量の金でも反応に寄与できる原子が多くなり、単位質量あたりの活性が高まります。また表面原子は配位数が少なく（配位不飽和）、反応分子を吸着・活性化しやすいため、活性点として働きやすくなります。S/V = 3/r の関係により、粒径を小さくするほど比表面積が増大します。
小さくしすぎる欠点（いずれか1つ）:
青色バックライト（450 nm）を用いるQDディスプレイで、緑色（約530 nm）と赤色（約630 nm）を発するCdSe量子ドットを設計します。本文3.5.2の計算結果を参考に、それぞれ必要なおおよその半径を推定し、あわせて量子ドットをディスプレイに用いる利点を有機蛍光体と比較して2つ述べてください。
解答:
半径の推定（3.5.2の表から内挿）:
| 目標発光波長 | 表の近い点 | 推定半径 |
|---|---|---|
| 緑 約530 nm | 2.0 nm→486 nm, 2.5 nm→554 nm | 約2.3〜2.4 nm |
| 赤 約630 nm | 3.5 nm→628 nm | 約3.5 nm |
緑色は半径約2.3 nm、赤色は半径約3.5 nmの量子ドットで得られます。青色バックライトの光を、緑・赤の量子ドットが吸収してそれぞれの色で再発光する方式です。
有機蛍光体に対する利点（いずれか2つ）:
補足: 実用ではカドミウムフリー化（InP系など）や、粒径分布を狭める合成制御が課題になります。
第3章では、ナノ材料の基礎（サイズ効果・次元性）と、カーボンナノチューブ・グラフェン・量子ドットの構造と物性、そしてPythonによる特性計算を学びました。次の第4章では、これらのナノ材料や先進材料を実デバイス・システムへ統合する設計原理を扱います。