なぜ量子化学なのか — 古典化学の限界
なぜ量子化学なのか — 古典化学の限界
古典電磁気学をそのまま当てはめれば原子は 秒で潰れる — 量子力学とはその崩壊を止める理論である。
はじめに
本書は 14 章で量子化学をひととおり通覧する。最初の章(0章)は、その出発点 — なぜ化学に量子力学が必要なのか — を一行で答えるための章である。読み終えたとき、読者は「古典化学では、なぜ原子が潰れないかすら説明できない」という事実を、次元解析と一つの積分で示せるようになり、1 章以降に出てくる波動関数・演算子・量子数がどんな欠落を埋めにくる道具なのかを頭の中に絵として持てるようになる。
本論 1 — ラザフォードのパズル
1911 年、ラザフォード(E. Rutherford) はアルファ粒子散乱の実験から原子の核モデルを確立した。極めて小さな正電荷の核のまわりを電子が回っている、という絵である。ところがこの絵を古典電磁気学にそのまま放り込んだ瞬間、問題が始まる。
加速する電荷は電磁波を放射しエネルギーを失う。その放射率は ラーモアの公式(Larmor formula) により
となる。ここで は電子の電荷、 は加速度、(イプシロン・ゼロ、真空の誘電率)、 は光速である。ボーア半径 で円運動する電子の速さはおよそ 。ここで (アルファ、微細構造定数)は
で定義される無次元量である。この速さで円を描く電子はラーモアの公式が命じるだけのエネルギーを失い続け、次元見積もりでは原子サイズから核へ吸い込まれるまで長くても 秒しかかからない。つまり古典電磁気学にしたがえば すべての原子は 1 ナノ秒も保たない。しかし実験では原子はほぼ永遠に安定している。古典化学は「なぜ原子は潰れないのか」というもっとも基本的な問いに答えを持たないのである。
本論 2 — 量子力学による一行の答
量子力学の答えは短い。原子のエネルギーは離散化(量子化)されており、その最も低い基底状態(ground state)はそれ以上下がる先がないので放射できない。 電子が核に吸い込まれないのは摩擦がないからではなく、「より低い状態」がそもそも存在しないからだ。
この絵の中で、ボーア半径の値は偶然ではない。電子が半径 の領域に閉じ込められているとすれば、量子力学的な運動エネルギーは概略
だけ持ち上がり(領域が狭いほど運動量の不確定性が大きい)、クーロン位置エネルギーは
だけ下がる。両者の和を について最小化すると、最適半径はちょうど
に落ちる。 は「これ以上圧縮されたくない運動エネルギー」と「もっと引き寄せたいクーロン引力」とのつり合い点なのである。この一行が本書全体のスローガンだ — 量子力学とは、原子が潰れないようにする理論である。
本論 3 — この本の 14 章でやること
本書は次のように進む。1 章で量子力学の言語(波動関数、演算子、測定)を定義し、2–3 章でシュレディンガー方程式を 1 次元箱・調和振動子・水素原子に適用する。4–6 章は多電子原子(スピン、パウリの排他原理、スレーター行列式)と周期表の量子的根拠。7–9 章は化学結合 — イオン性・共有結合・分子軌道法(MO)と原子軌道(AO)。10–11 章は分子の振動・回転・電子スペクトル。12–14 章は固体とバンド構造、そして最後に量子化学が現代の材料・触媒・創薬にどうつながっているかを見る。
前提はだいたい次のとおりである。高校化学(周期表、価電子)、1 変数・多変数の微積分を少し、複素数 — そして好奇心。量子力学の予備知識は仮定しない。波動関数、演算子、量子数、スレーター行列式、MO/AO といった用語はすべて初出のときに二文ほどの平易な解説とともに定義する。
Pythonで確かめる
# 古典 vs 量子: ボーア半径と古典的崩壊時間を一度に計算する。
import numpy as np
hbar = 1.0545718e-34 # J·s
me = 9.1093837e-31 # kg
e = 1.602176634e-19 # C
eps0 = 8.8541878128e-12 # F/m
c = 2.99792458e8 # m/s
# 量子の予測: ボーア半径
a0 = 4 * np.pi * eps0 * hbar**2 / (me * e**2)
print(f"ボーア半径 a0 = {a0*1e12:.2f} pm")
# 古典的崩壊時間 (ラーモアの公式 + ビリアル定理、オーダー見積もり)
# E_K = e^2 / (8 pi eps0 a0), 速度 v = alpha * c
alpha = e**2 / (4 * np.pi * eps0 * hbar * c)
v = alpha * c
a_acc = v**2 / a0 # 向心加速度
P = (1 / (6 * np.pi * eps0)) * e**2 * a_acc**2 / c**3
E_K = e**2 / (8 * np.pi * eps0 * a0)
tau = E_K / P # エネルギー/放射率
print(f"古典的崩壊時間 tau ~ {tau:.2e} s")
print("実験: 原子は事実上永久に安定 — このギャップこそ量子力学の動機。")
出力では がおよそ 53 pm、 が 秒の近くに落ちる。古典理論は「原子は一瞬も保たない」と言い、実験は「永遠に保つ」と言う。このギャップを 14 章かけて埋めていくのが本書である。
次章へ
1章: 量子化学の言語では、波動関数 、演算子 、測定値と固有値方程式 の意味を定義する。この章で出た「基底状態にはそれより下がない」というスローガンが、次章ではエルミート演算子の離散スペクトルという数学的事実として再登場することになる。