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

幾何構造最適化における初期ヘシアン推定の数理的基礎と拡張:Schlegel Hessianについて
2026-01-25
幾何構造最適化の収束効率を決定づける初期ヘシアン推定法(Schlegel Hessian)について、H. Bernhard Schlegelによる1984年の提唱から1997年の重元素への拡張に至るまでの理論的背景、数理的アルゴリズム、および経験的パラメータ決定のプロセスを詳述する。原子価座標系を用いた力場構築と座標変換の数学的定式化に焦点を当てる。
冗長内部座標系における自動化された鞍点探索アルゴリズムの数理と実装:Sellaについて
2026-01-25
ポテンシャルエネルギー曲面(PES)上の一次鞍点(遷移状態)を探索するための新規アルゴリズム『Sella』について、その数学的基礎、座標変換の幾何学的処理、およびNull-Space SQPを用いた制約付き最適化手法を詳述する。特に、自動化を阻害する直線結合角問題へのダミー原子を用いた対処法と、反復的ヘシアン対角化による計算コスト削減効果に焦点を当てる。
振動解析における幾何学的基礎: 質量荷重ヘシアンの導出と配置空間の平坦性に関する考察
2026-01-24
分子振動解析において中心的役割を果たす質量荷重座標系への変換について、その数理的導出、計量テンソルの定義、および空間の平坦性(曲率の不在)に関する幾何学的正当性を詳細に解説する。
深層学習による交換相関汎関数の構築:Skalaアーキテクチャの数理と物理的意義
2026-01-24
Microsoft Researchによって提案された深層学習ベースの交換相関汎関数「Skala」について、その数理的構造、物理的制約条件の充足、および既存の「ヤコブの梯子」に対する位置づけを論じる。局所特徴量からの非局所相互作用の学習と計算コストの並立に関する理論的背景を詳述する。
電子相関の二分法:静的相関と動的相関の定義、弊害、および解決手法
2026-01-21
Hartree-Fock法で欠落する「電子相関」を、近距離相互作用(動的)と多配置性(静的)の観点から分類・定義する。単一参照理論が破綻するメカニズム(静的相関の無視)と、化学精度が得られない理由(動的相関の無視)を対比し、CASSCFやCCSD(T)などの対応手法をマッピングする。
基底関数系の収束性と誤差要因:BSSEとCBS極限
2026-01-21
量子化学計算における離散化誤差の根源である「基底関数系」の選択指針。Pople系とDunning系の設計思想の違い、基底関数重なり誤差(BSSE)の補正法、および完全基底系(CBS)極限への外挿について記述する。