Home
974 words
5 minutes
原子単位系 (Hartree Atomic Units):定義、物理的意義、および換算係数

last_modified: 2026-01-21

1. 概要 (Overview)#

量子力学的な微視的スケールにおいて、SI単位系(kg, m, s, J)は数値の桁が極端に小さくなり(例:電子質量 1031\approx 10^{-31} kg)、取り扱いに不便である。 また、数式中に物理定数(,me,e,4πϵ0\hbar, m_e, e, 4\pi\epsilon_0)が頻出し、記述が煩雑となる。 これを解決するため、計算化学ではこれらの基本定数をすべて 「1」 と置く Hartree原子単位系 (a.u.) を採用する。これにより、方程式は無次元化され、数値計算上の誤差(アンダーフロー等)も抑制される。

2. 基本単位の定義 (Base Units)#

以下の4つの物理定数を、原子単位系における「1」と定義する。

物理量記号定義 (a.u.)SI単位系での値 (近似値)
質量mem_e11電子静止質量 (9.109×10319.109 \times 10^{-31} kg)
電荷ee11素電荷 (1.602×10191.602 \times 10^{-19} C)
作用 (角運動量)\hbar11ディラック定数 (1.055×10341.055 \times 10^{-34} J s)
静電定数1/4πϵ01/4\pi\epsilon_011クーロン定数 (8.987×1098.987 \times 10^9 N m2^2 C2^{-2})

この定義により、非相対論的シュレーディンガー方程式(水素原子)は以下のように極限まで単純化される。

[22me2e24πϵ0r]Ψ=EΨa.u.[1221r]Ψ=EΨ\left[ -\frac{\hbar^2}{2m_e}\nabla^2 - \frac{e^2}{4\pi\epsilon_0 r} \right] \Psi = E \Psi \quad \xrightarrow{\text{a.u.}} \quad \left[ -\frac{1}{2}\nabla^2 - \frac{1}{r} \right] \Psi = E \Psi

3. 組立単位と換算係数 (Derived Units and Conversions)#

基本単位の組み合わせにより導出される、実用上重要な単位群である。 ※以下の換算値は2022年CODATA推奨値に基づく。

3.1 長さ (Length): Bohr#

  • 単位名: Bohr (または a0a_0)
  • 定義式: a0=4πϵ02mee2a_0 = \frac{4\pi\epsilon_0 \hbar^2}{m_e e^2}
  • 物理的意味: 水素原子の基底状態におけるボーア半径。
  • 換算:
    • 1 Bohr0.529177 A˚1 \text{ Bohr} \approx \mathbf{0.529177} \text{ \AA}
    • 1 Bohr5.29177×1011 m1 \text{ Bohr} \approx 5.29177 \times 10^{-11} \text{ m}

3.2 エネルギー (Energy): Hartree#

  • 単位名: Hartree (または EhE_h, Ha)
  • 定義式: Eh=mee4(4πϵ0)22=2mea02E_h = \frac{m_e e^4}{(4\pi\epsilon_0)^2 \hbar^2} = \frac{\hbar^2}{m_e a_0^2}
  • 物理的意味: 水素原子の基底状態のポテンシャルエネルギーの絶対値(結合エネルギーの2倍)。
  • 換算 (重要):
    • 1 Hartree27.2114 eV1 \text{ Hartree} \approx \mathbf{27.2114} \text{ eV}
    • 1 Hartree627.509 kcal/mol1 \text{ Hartree} \approx \mathbf{627.509} \text{ kcal/mol}
    • 1 Hartree2625.50 kJ/mol1 \text{ Hartree} \approx 2625.50 \text{ kJ/mol}
    • 1 Hartree219474.6 cm11 \text{ Hartree} \approx 219474.6 \text{ cm}^{-1} (波数)

3.3 時間 (Time)#

  • 単位名: atomic unit of time
  • 定義式: t0=Eht_0 = \frac{\hbar}{E_h}
  • 物理的意味: 電子がボーア軌道を1ラジアン回転するのに要する時間。極めて短時間であり、アト秒 (101810^{-18} s) 領域のダイナミクスで重要となる。
  • 換算:
    • 1 a.u.2.41888×1017 s1 \text{ a.u.} \approx \mathbf{2.41888 \times 10^{-17}} \text{ s} (24.19\approx 24.19 as)

3.4 速度 (Velocity)#

  • 単位名: atomic unit of velocity
  • 定義式: v0=a0t0=e24πϵ0v_0 = \frac{a_0}{t_0} = \frac{e^2}{4\pi\epsilon_0 \hbar}
  • 物理的意味: 水素原子の基底状態における古典的な電子速度。光速 cc との関係は微細構造定数 α1/137\alpha \approx 1/137 で表される(v0=αcv_0 = \alpha c)。
  • 換算:
    • 1 a.u.2.18769×106 m/s1 \text{ a.u.} \approx 2.18769 \times 10^6 \text{ m/s}

3.5 双極子モーメント (Electric Dipole Moment)#

  • 単位名: atomic unit of dipole moment
  • 定義式: ea0ea_0
  • 物理的意味: 単位電荷 +1+11-1 が 1 Bohr 離れた時の双極子モーメント。
  • 換算:
    • 1 a.u.2.54174 Debye (D)1 \text{ a.u.} \approx \mathbf{2.54174} \text{ Debye (D)}

4. 注意点 (Notes)#

  1. Rydberg単位との混同: 分光学ではRydberg単位 (1 Ry=1/2 Hartree13.6 eV1 \text{ Ry} = 1/2 \text{ Hartree} \approx 13.6 \text{ eV}) が用いられることがあるが、量子化学プログラム(Gaussian, GAMESS等)の出力は原則として Hartree である。係数 1/21/2 の取り違えに注意すること。
  2. 分子構造の出力: 座標データは入力がAngstromであっても、内部計算および出力ログ(アーカイブ)ではBohrに変換されている場合が多い。可視化ソフトなどで読み込む際に単位系の指定を誤ると、分子が約半分 (0.5290.529倍) に潰れて表示される。

結論 (Conclusion)#

計算化学プログラムの出力ログを解析する際は、以下の換算係数を定数として実装することが推奨される。

  • Length: Bohr * 0.529177 -> Angstrom
  • Energy: Hartree * 627.509 -> kcal/mol (反応熱議論用)
  • Energy: Hartree * 27.2114 -> eV (軌道エネルギー議論用)
原子単位系 (Hartree Atomic Units):定義、物理的意義、および換算係数
https://ss0832.github.io/posts/20260121_compchem_unit/
Author
ss0832
Published at
2026-01-21
License
CC BY-NC-SA 4.0

Related Posts

Kabsch アルゴリズムによる最小二乗重ね合わせと RMSD 算出:数学的導出・歴史的背景・実利的応用
2026-03-16
分子構造の比較において中心的な役割を果たす Kabsch アルゴリズムについて、特異値分解(SVD)に基づく数学的導出、歴史的成立過程、実装上の注意点(キラリティ処理を含む)、および計算化学・構造生物学における実利的成果を包括的に解説する。
O1NumHess: O(1) 勾配評価によるセミ数値的ヘシアン行列算出法の数理的構造と実装
2026-02-15
分子系のヘシアン行列算出において、原子数に依存しない定数オーダー O(1) の勾配計算回数で高精度な近似を実現する O1NumHess アルゴリズムについて、その数理的背景(Off-Diagonal Low-Rank 性質)、実装の詳細、および計算精度とコストに関する実証的評価を包括的に解説する。
密度汎関数理論における交換相関汎関数の系統的最適化:B97汎関数の数学的構造と有限級数展開による拡張性
2026-02-04
Axel D. Beckeによる1997年の論文 (J. Chem. Phys. 107, 8554) に基づき、B97汎関数の設計思想とその数学的定式化について詳述する。特に、勾配補正項の有限級数展開(冪級数展開)を用いた系統的な最適化手法、半無限区間の有限変数への写像、およびG2テストセットを用いたパラメータ決定プロセスに焦点を当てる。また、この手法が示唆するGGA/正確交換混合形式の精度の限界と、物理的妥当性と過剰適合(オーバーフィッティング)の境界に関する考察を行う。
デカルト座標系における角度ポテンシャル勾配の特異点除去:数値的正則化手法の数理的導出と実装
2026-02-02
分子動力学および構造最適化計算において、3点間角度ポテンシャルの評価時に発生する特異点(ゼロ除算および桁落ち)の問題を数学的に解析し、テイラー展開を用いた正則化手法(Regularization)を詳細に論じる。Szabo-OstlundのBoys関数処理との類似性、誤差解析、およびPythonによる参照実装を含め、学術的・技術的観点から包括的に解説する。
冗長内部座標系における測地線アプローチを用いた幾何構造最適化の加速
2026-01-29
Hermesら(2021)によって提案された、冗長内部座標多様体上の測地線を利用した構造最適化手法の数理的定式化、アルゴリズムの実装、およびその有効性について詳細に解説する。従来のニュートン法との幾何学的差異と、微分幾何学的アプローチによる収束性の改善に焦点を当てる。
有効曲率の幾何学的導出と実装オプション
2026-01-29
論文(The Journal of Chemical Physics 1997, 107 (2), 375–384.)の「多項式フィッティング」と、ライプニッツ則に基づく「解析的分解」の等価性を証明し、実務的な実装分岐を考察する。