混合基底法による効率的なバンド計算
第一原理
バンド計算は、物質の電子状態を量子力学に基づいて計算する手法として広く用いられています。この計算において基底関数の選択は計算精度と効率に大きな影響を与えます。
平面波基底は周期境界条件下で扱いやすく、多くの計算パッケージで採用されていますが、原子核のポテンシャルが深い場合、計算精度を確保するために非常に多くの
平面波が必要となり、計算コストが急激に増加します。特に、炭素や3d遷移金属など、原子核のポテンシャルが深い元素を含む系では、この問題が顕著になります。
そこで、計算効率を向上させるために考案された手法が混合基底法です。混合基底法は、
平面波基底に加えて、
局在基底を組み合わせて用いることで、計算コストを削減します。
局在基底としては、ガウス型基底がしばしば用いられます。
平面波基底は、空間全体に広がる波動関数で、周期的なポテンシャルに対して効率的に計算できます。一方、
局在基底は原子核の近傍に局在した波動関数であり、原子核ポテンシャルの急峻な変化を精度良く記述できます。混合基底法では、原子核近傍の電子状態を
局在基底で、それ以外の領域を
平面波基底で記述することで、計算精度を維持しつつ必要な基底関数の数を削減します。
ガウス型基底を用いる利点は、行列要素の計算を解析的に解ける点にあります。
平面波基底同士の積分は数値積分を用いる必要があるのに対し、ガウス型基底では解析解が得られるため、計算速度の大幅な向上に繋がります。これは計算時間の大幅な削減に貢献します。
しかし、混合基底法には課題も存在します。
平面波基底と
局在基底は異なる性質を持つため、それらを組み合わせた基底関数の集合は
正規直交系ではなくなります。この非直交性は、計算の複雑化を招きます。具体的には、行列の対角化や、力やストレスといった物理量の計算が困難になります。これらの物理量は、物質の構造や物性を理解するために重要であり、その計算の複雑化は、プログラムのコード化や、実際の計算において大きな障壁となります。
混合基底法は、計算効率の向上と精度の確保という相反する要求のバランスを取るための高度な手法です。計算コストの削減効果は、対象とする系や計算精度によって大きく変化します。そのため、計算対象や目的に応じて、適切な基底関数を選ぶことが重要です。混合基底法を用いた計算を行う際には、非直交性による計算の複雑化を考慮し、適切な計算手法を選択する必要があります。また、Pulay補正などの手法を用いて計算精度を向上させることも重要になります。