Excited-state Properties Beyond the Excitation Energy from Orbital-Optimized Density Functional Calculations II: Absorption Spectra

本研究は、非直交行列式に対するLöwdinの形式論を、射影型拡張波（PAW）法を用いた軌道最適化密度汎関数計算へと拡張するものであり、この手法が単一行列式の励起状態については吸収スペクトルやピーク強度を定性的に再現する一方で、多配置状態には苦慮し、また完全交換や自己相互作用補正による系統的な改善が見られないことを示している。

要約

全体像：踊る分子のスナップショットを撮る

分子を、小さくて複雑なダンス・グループだと想像してみてください。光を当てると、ダンサー（電子）はより高いエネルギーレベルへと跳ね上がり、ダンスの動きを変えます。科学者たちは、この「ジャンプ」について2つのことを予測したいと考えています。

1. エネルギー： ジャンプするためにどれくらいのエネルギーが必要か？（ジャンプの高さのようなもの）
2. 明るさ： ジャンプした時の光のフラッシュはどれくらい明るいか？（これは「振動子強度」または「吸収強度」と呼ばれます）

長い間、コンピュータプログラムはジャンプの「高さ」を予測することには長けていましたが、「明るさ」については間違えることがよくありました。この論文では、これらのジャンプを計算するための新しい手法である軌道最適化（Orbital-Optimized: OO）密度汎関数理論を紹介しています。

問題点：「非直交」なダンスフロア

標準的なコンピュータモデルでは、通常、「基底状態」（休息中のダンス）と「励起状態」（ジャンプ中のダンス）は、別々の部屋に立っている二人のように、完全に分離しており、重なり合わないと仮定します。これにより、数学的な計算が簡単になります。

しかし、この新しい「軌道最適化」手法は、より現実的です。これは、電子がジャンプに合わせて特別に再配置されることを許容します。問題は、この再配置によって、「休息」状態と「ジャンプ」状態がもはや別々の部屋ではなくなり、同じ部屋に入り、わずかに重なり合うようになることです。

比喩： ハグをしている二人の間の距離を測ろうとしている場面を想像してください。もし二人が離れて立っていると仮定してしまえば、測定結果は間違ったものになります。これらの状態が「ハグ（重なり）」しているため、光のフラッシュの明るさを計算するのは非常に難しくなります。これまでの研究は主にジャンプの「高さ」が正しいかどうかをチェックしてきましたが、状態が重なり合っている時に「明るさ」の計算が正確かどうかまではチェックしていませんでした。

彼らがしたこと：新手法のテスト

研究者たちはこの新手法を用い、少数の分子（水、アンモニア、ホルムアルデヒド、メタノール、エチレン）を用いてテストを行いました。彼らは以下の点を確認しようとしました。

1. この手法は、光のフラッシュの「明るさ」を正しく予測できるか？
2. 計算を行う際に使う「道具」（数学的基底関数系）は影響するのか？
3. 「ゲームのルール」（電子の振る舞いを記述する数学的公式）を変えることで、明るさのエラーを修正できるのか？

結果：一筋縄ではいかない結果

得られた結果を簡単に説明すると、以下の通りです。

1. 「道具箱」が重要（基底関数系）
基底関数系を、カメラの「解像度」だと考えてください。

• 低解像度（単純な道具）： 単純なカメラを使うと、「拡散」した電子（リュードベリ状態のように、分子から遠くに漂っている電子）の詳細を見逃してしまいます。そのため、明るさの計算が少しずれてしまいます。
• 高解像度（複雑な道具）： 非常に詳細な道具（プレーンウェーブやダブル増強基底関数系など）を使用したとき、明るさの計算結果は非常に一貫したものになりました。
• 教訓： 明るさを正しく捉えるには、特に遠くに漂う電子に対して、高解像度のカメラが必要になります。

2. 「ルール」だけでは解決しない（交換相関汎関数）
研究者たちは、シミュレーションの「ルール」（PBE、PBE0、あるいは自己相互作用補正を加えたものなどの数学的公式）を変えて試みました。

• 結果： ルールを変えることは、ジャンプの「高さ（エネルギー）」を修正する助けにはなりましたが、明るさのエラーを一貫して修正することにはなりませんでした。
• 比喩： ぼやけた写真を修正しようとしている場面を想像してください。レンズを変えたり、照明を変えたり、フィルターをかけたりしました。時には写真が鮮明になることもありましたが、多くの場合、明るさは依然として間違ったままでした。あらゆる分子に対して明るさを修正できるような「魔法のルール」は存在しませんでした。

3. 真の犯人は「ソロ」か「グループ」か（単一構成 vs 多構成）
これが最も重要な発見です。この手法は、励起状態が**「ソロダンサー」**（単一で明確な構成）である場合には非常にうまく機能しました。

• ソロダンサー： アンモニアのような単純なジャンプの場合、この手法は明るさを完璧に予測しました。
• グループダンサー： 電子のダンスが多くの異なる可能性が混ざり合った複雑な状態（多構成状態）である場合、この手法は明るさを正しく予測できませんでした。
• 具体的な失敗： 水、ホルムアルデヒド、エチレンにおいて大きなエラーが見られました。これらのケースでは、励起状態はさまざまなダンスの動きが入り混じった複雑なものです。コンピュータモデルは、その状態を一つのクリーンな「ソロ」の動きとして強制してしまうため、明るさを誤ってしまうのです。
• 重なりの問題： 彼らは、基底状態と励起状態の「ハグ（重なり）」がエラーの原因であるかどうかを調べました。しかし、重なり具合を変えても明るさのエラーは変わらないことが分かりました。つまり、重なりが主因ではなく、「ソロかグループか」という性質が原因だったのです。

旧来の手法（LR-TDDFT）との比較

彼らはこの新手法を標準的な手法（LR-TDDFT）と比較しました。

• 標準的な手法： 複雑な「グループ」のダンス（エチレンで見られるような明るいフラッシュ）の明るさを予測することには長けていますが、遠くに漂う「ソロ」ダンサー（リュードベリ状態）のエネルギーを予測するのは苦手です。
• 新しい手法（OO）： 「ソロ」ダンサー（リュードベリ状態）のエネルギーを予測することには優れていますが、複雑な「グループ」のダンスの明るさを予測することには苦戦します。

結論

この論文は、新しい「軌道最適化」手法が分子の光吸収を予測するための強力なツールであることを示していますが、「ただし、条件があります」。

• 単純な「ソロの動き」（リュードベリ状態）には非常によく機能します。
• しかし、ダンスが多くの動きが混ざり合った複雑なもの（多構成状態）になると、苦戦します。
• 数学的な公式や「カメラの解像度」を変えるだけでは、これらの複雑なダンスによるエラーを修正することはできません。これらを解決するには、単なるソロではなく、複雑なグループダンスを扱える手法が必要なのです。

要するに、この手法は電子のジャンプの「高さ」と単純なジャンプの「明るさ」を理解するための素晴らしい一歩ですが、最も複雑で混沌とした電子のダンスを正確に予測するためには、まださらなる助けが必要です。

技術概要

技術要約：励起エネルギーを超えた励起状態の性質：軌道最適化密度汎関数計算 II：吸収スペクトル

問題提起
軌道最適化（OO）密度汎関数計算は、励起状態を基底状態の線形応答として扱うのではなく、各状態に対して変分的に最適化された軌道を用いることで、励起状態を記述する状態特異的なアプローチを提供する。この手法は、標準的な線形応答時間依存密度汎関数理論（LR-TDDFT）が、不正確な漸近ポテンシャルの挙動により Rydberg 状態の記述に失敗することが多い中で、励起エネルギーの改善に寄与することが示されている。しかし、これは遷移特性の計算において重大な困難をもたらす。各電子状態は異なる分子軌道のセットによって記述されるため、OO 状態は一般に非直交である。この非直交性は、遷移双極子モーメント（TDM）および振動強度の評価を複雑にする。これまでのベンチマークは、主に励起エネルギー、あるいは最初の励起状態（HOMO–LUMO 遷移）に限定された振動強度の計算に焦点を当ててきた。その結果、より広い範囲の原子価および Rydberg 励起における OO 法の振動強度の性能、および非直交性がこれらの特性に与える影響については、十分に理解されていない。

手法
本研究では、一連の小分子（H₂O, CH₂O, NH₃, CH₃OH, C₂H₄）を対象とし、価電子状態および Rydberg 状態を 10 eV までカバーする OO 密度汎関数計算の吸収スペクトルへの適用を拡張する。

• 理論的枠組み: 著者らは、変分崩壊を回避するために、エネルギー汎関数の停留点（鞍点）として励起状態を特定する直接的な軌道最適化スキームを採用している。結果として得られる状態の非直交性を扱うために、Löwdin の非直交行列式に関する定式化を射影増加平面波（PAW）フレームワークへと拡張している。これにより、非直交な Slater 行列式間の行列要素を厳密に計算することが可能となる。
• 実装: この手法は GPAW ソフトウェア内に実装されている。遷移行列要素は、滑らかな擬波動関数と全電子波動関数の間の変換を考慮しつつ、PAW フレームワーク内で導出される。この定式化には、増加球内における重なり行列および一電子演算子（具体的には双極子演算子）に対する補正が含まれる。
• 計算の詳細: 計算は、様々な交換相関（xc）汎関数（一般化勾配近似である PBE、ハイブリッド汎関数である PBE0、およびグローバルスケーリング因子（$\alpha = 1$ および $\alpha = 1/2$）を用いた明示的な Perdew–Zunger 自己相互作用補正（SIC）を伴う PBE）を用いて行われた。平面波（PW）および原子軌道線形結合（LCAO）（aug-cc-pVDZ および d-aug-cc-pVDZ）の両方の基底関数が PAW フレームワーク内で使用された。結果は、高レベルの多参照リファレンスデータ（MS-CASPT2, CCSD, CCSD(T), および CCSDT）および LR-TDDFT 計算と比較された。
• スピン純化: 開殻一重項状態については、混合スピンおよび三重項計算から導出されたスピン純化公式を用いて、エネルギーおよび TDM を取得した。

主要な貢献

1. Löwdin の規則の PAW への拡張: 本論文は、平面波および格子ベースの表現を用いて OO 状態の遷移特性を計算することを可能にするために、非直交行列式に対する Löwdin の規則を PAW フレームワーク内で特別に導出し、実装している。
2. 系統的なベンチマーク: 本研究は、多様な分子に対する複数の励起状態（HOMO–LUMO 遷移を超えて）の OO 振動強度を包括的に評価し、異なる基底関数および xc 汎関数の影響を検証している。
3. 非直交性の分析: 本研究は、OO 計算における基底状態と励起状態の間のゼロでない重なりが、振動強度の系統的な誤差を導入するかどうかを明示的に調査している。これは、以前の文献で提起された懸念であるが、完全な吸収スペクトルについてはまだ解決されていない。

結果

• 基底関数依存性: 拡散的な Rydberg 状態の振動強度は、基底関数の表現に敏感である。単一に増強された LCAO 基底関数（aug）は、二重に増強された（d-aug）および平面波（PW）表現と比較して、励起エネルギーを過大評価し、高度に拡散した状態に対して不正確な振動強度を与えることが多い。しかし、明るい価電子遷移は基底関数に対してほとんど敏感ではない。
• 汎関数依存性: xc 汎関数（PBE, PBE0, PBE-SIC/2, PBE-SIC）を変更すると、励起エネルギーおよび特定の振動強度は変化するが、全データセットにわたって精度の系統的な階層を生み出すことはない。SIC は漸近ポテンシャルを改善するが、振動強度の誤差を一貫して修正することはない。
• 単一構成 vs 多構成的性質: 励起状態の電子的な性質に基づいて明確な区別が現れる。
  • 単一構成的状態: 単一の行列式によって支配される状態（例：アンモニアの多くの Rydberg 状態）については、OO 計算はリファレンスの振動強度およびピーク位置と極めて良好な一致を示す。
  • 多構成的状態: 強力な多構成的性質または顕著な構成混合を持つ状態（例：H₂O S3/S4, CH₂O S5, C₂H₄ S1）については、OO 計算は大幅な不一致を示す。具体的には、単一行列式による記述は正しい強度分布を捉えることができず、テストされたすべての汎関数にわたって、振動強度の過大評価または過小評価が生じる。
• 非直交性と重なり: 基底状態と励起状態の重みの大きさと振動強度の誤差との間に、明確な相関は見られない。例えば、ハイブリッド汎関数や SIC によって重みを減少させても、H₂O S3 のような問題のある状態における振動強度が系統的に改善されることはなかった。したがって、非直交性は、最大の偏差を生じさせる誤差の主要な原因とは特定されなかった。
• LR-TDDFT との比較: OO 法は、Rydberg 状態のエネルギーおよび強度を記述する上で、一般に LR-TDDFT よりも優れた性能を示す。逆に、強い構成混合を伴う明るい価電子遷移については、LR-TDDFT がより優れた記述を提供することが多い。これは、単一行列式による OO アプローチが苦戦する領域である。

意義および結論
本論文は、OO 密度汎関数計算が、単一構成的な Rydberg 状態の励起エネルギーおよび吸収スペクトルの予測において堅牢なツールである一方で、その振動強度の定量的信頼性は、コーン・シャム記述の単一行列式という性質によって制限されると結論付けている。この手法は、励起状態が支配的な単一構成によって適切に記述されている場合には良好に機能するが、強力な多構成的性質を持つ状態に対しては、汎関数や基底関数の選択にかかわらず失敗する。

著者らは、特定の遷移（例：水 S3/S4, ホルムアルデヒド S5, エチレン S1）で見られる大きな振動強度の誤差は、非直交性や基底関数の不完全性のアーティファクトではなく、単一行列式近似に固有のものであると主張している。したがって、多構成的性質が顕著な系において、OO-KS スペクトルは強度の観点からは定性的（例：ピーク位置）に解釈されるべきである。本研究は、今後の改善には、単に汎関数や基底関数を修正することではなく、明示的に多構成的励起状態を扱うように OO フレームワークを拡張する必要があることを示唆している。