OpenMXについて

OpenMX（Open source package for Material eXplorer）は密度汎関数理論（DFT）[1]、 ノルム保存型擬ポテンシャル[32,33,34,35,36]および擬原子基底関数[41]に基づき、 原子レベルから物質の第一原理シミュレーションを実行するためのソフトウェア・パッケージです。 OpenMXで用いられている計算手法、アルゴリズム、またそのプログラム上での実装は並列コンピュータ上でのMPI並列や OpenMP/MPIハイブリッド並列による大規模第一原理電子状態計算を実現するために、入念に設計されています。 OpenMXにおいて実現されたDFTの効率的な実装により、固体、表面、界面、液体、そして低次元物質などの幅広い物質の 電子構造、磁気構造、また幾何構造を第一原理的に計算することが可能です。 数百コアを有する並列コンピュータを使用することにより、1000個の原子で構成される系を通常の対角化手法を用いて扱うことができます。 さらに数千CPUのコアを有する並列コンピュータを使用すれば、一万個以上の原子から構成される系の第一原理による電子状態計算さえも、OpenMXに実装されているO($N$)法を用いて実行可能です。 多くの元素に対して最適化された擬ポテンシャルおよび基底関数がデータベースとして整備されており、またその精度はベンチマーク計算に よって検証されています。そのためユーザは自分でこれらのデータを準備する必要がなく、速やかに各自のシミュレーションを開始することができます。 OpenMXには磁気特性、誘電特性、電気伝導特性、反応エネルギー障壁などの様々な物理的・化学的特性を計算するための機能が実装されており、 量子力学の第一原理からナノスケールの物質群を深く理解するための強力なツールとして、幅広い活用が 期待されます。 これまでの応用計算の事例はOpenMXのWebサイト(http://www.openmx.-square.org/)に記載されていますので、参考にして下さい。 OpenMXの開発は、2000年に尾崎グループにより開始され、その後、本マニュアルの冒頭のリストに挙げられている数多くの開発者が 本オープンソースパッケージの開発に貢献してきました。 プログラム・パッケージとソースコードはGNU General Public License version 3 (GPLv3) [102]に準じて配布されており、 OpenMXのWebサイトからダウンロードすることができます。

OpenMX Ver. 3.9 の特徴と機能は以下のようになります。

• クラスタ計算法、バンド計算法、O($N$)法、低次スケーリング法による全エネルギーと力の計算
• 交換相関エネルギーに対する局所密度近似（LDA, LSDA）[2,3,4]および一般化勾配近似（GGA）[5]
• DFT+$U$法[20]
• ノルム保存型擬ポテンシャル [2,33,34,36]
• 変分的に最適化された擬原子基底関数 [41]
• 擬ポテンシャル法の枠組みでの完全および半相対論的取り扱い [12,32,16]
• ノンコリニアDFT法[8,9,10,11]
• ノンコリニアDFT法におけるスピン磁気モーメントおよび軌道磁気モーメントの方位制御 [13]
• Berry位相法による巨視的分極率の計算[15]
• O($N$)分割統治法（DC法）[50]
• 局在化自然軌道によるO($N$)分割統治法（DC-LNO法）[51]
• O($N$)Krylov部分空間法
• ELPA [39]による並列固有値ソルバ
• 単純混合法、RMM-DIIS法 [58]、GR-Pulay法 [57], Kerker法 [59]、Kerker因子による重み付きRMM-DIIS法 [58]による電荷密度の混合、 ハミルトニアン行列に対するRMM-DIIS法による電荷密度の混合 [58]
• 交換結合パラメータ [17,18]の計算
• 有効遮蔽媒質法(ESM) [125,128]
• 走査型トンネル顕微鏡（STM）シミュレーション [71]
• DFT-D2、DFT-D3法によるVan der Waals相互作用 [135,136,137]
• バンドアンフォールディング法 [142]
• Nudged Elastic Band法（NEB法）[72]
• 固体中の内殻軌道の絶対束縛エネルギーの計算[88]
• 光伝導率と誘電関数[98]
• 電荷ドーピング
• のこぎり波による均一電界の導入
• 完全および制限付き構造最適化
• 非平衡グリーン関数（NEGF）法による電気伝導計算 [73]
• 最大局所化ワニエ関数の計算
• NVEアンサンブル分子動力学
• 速度スケーリング法 [30]およびNose-Hoover法 [31]によるNVTアンサンブル分子動力学
• マリケン(Mulliken)法、ボロノイ(Voronoi)法、ESP法による電荷・スピン数の解析
• 自然軌道解析 [7]
• 波動関数および電子・スピン密度分布の解析
• バンド計算によるバンド分散の解析
• 状態密度（DOS）および射影DOSの計算
• 柔軟なデータ書式(flexible data format)による入力ファイル
• BoltzTrapへのインタフェース [100,101]
• Wannier90へのインタフェース [145]
• 電荷密度などを可視化するためのXCrySDenへのインタフェース [105]
• 完全な動的メモリ割当て
• メッセージ・パッシング・インタフェース (MPI)による並列計算
• OpenMPによる並列計算
• 開発者向けの利便性の高いユーザーインタフェース

コリニアDFT法ではスカラー相対論的擬ポテンシャルが、またノンコリニアDFT法では完全相対論的擬ポテンシャルが 利用可能です。 またスピンおよび軌道磁気モーメントを制御するために制限ノンコリニアDFT法が実装されています。 これらの手法は複雑なノンコリニア磁気構造やスピン軌道相互作用を調べるのに役に立ちます。 通常の対角化による計算は、数千コアまでの並列化が実現できるELPAに基づく並列固有値ソルバ [39]及びScaLAPACKによって行われます。 この機能により、通常の対角化法を用いて1000個の原子からなる系の計算が可能になります。 この高並列化対角化手法によってクラスタ、分子、スラブおよび固体に対する大規模計算が可能ですが、線形スケーリング法および 低次スケーリング法も固有値ソルバとして利用可能です。これらの低次スケーリング法に対しては計算精度と計算効率に対して 注意深い検討を行うことで、1万原子を越える系をも取り扱うことが可能です。 結晶格子可変の構造最適化やバンドアンフォールディング法も有用な機能です。 またOpenMX Ver. 3.9の重要な新機能として、分割統治法と局在自然軌道に基づく新しいO($N$)法と X線光電子分光法(XPS)で観測される内殻電子束縛エネルギーの絶対値計算が挙げられます。 OpenMXの開発は継続して行われています。本オープンソースコードの発展に寄与して頂ける開発者の方のご参加を歓迎致します。