ホーム > 研究活動 > 口頭発表(2013) > 大規模スーパーセルを対象とした遮蔽KKR 法に基づく第一原理計算コードの開発/Development of a first-principles code for large super-cells based on the screened KKR method

研究活動

大規模スーパーセルを対象とした遮蔽KKR 法に基づく第一原理計算コードの開発/Development of a
first-principles code for large super-cells based on the screened KKR method

文部科学省科研費新学術領域「コンピューティクスによる物質デザイン:複合相関と非平衡ダイナミクス」

2013年7月8日(月)

S. Doi(ISSP), H. Akai(ISSP) ( Tokyo Univ )

Abstract

The KKR Green's function method is a method for first-principles electronic structure calculations.
In this method, instead of solving eigenvalue problems, Green's function of the Kohn-Sham equation is directly calculated.
By taking the imaginary part of the Green's function, we can obtain the electron density very accurately and efficiently.
The Green's function is calculated from that of an arbitrary reference system by solving the Dysontype equation.
In the process of solving the Dyson equation, inversions of matrices of a dimension proportional to N (N: the number of atoms in a unit-cell) is required.
In the usual KKR method, the free space is chosen as the reference system.
When N becomes large, the computational cost increases drastically.
This is because O(N3) operations are required to invert a N×N matrix.
To overcome this difficulty, the concept of screened KKR was proposed.
In the screened KKR method, a system with repulsive muffin-tin potentials of a constant height is chosen as the reference system.
In this reference system, the Green's function decays exponentially in real-space if the energy range is sufficiently lower than the repulsive potential height.
Then we can ignore long range propagations without sacrificing the accuracy.
In this case, the matrix that we must invert becomes sparse.
Therefore, one can expect large improvement of the calculational efficiency by use of an efficient sparse matrix solver.
This technique has been mainly applied to the multilayered systems which extend to only one dimensional direction assuming two-dimensional periodicity.
For such systems, the matrix becomes block tridiagonal and the computational cost is exactly proportional to the layer thickness.
It is desirable to apply the screened KKR technique to more general and complex systems which extend to all three dimensional directions.
For such general large super-cells, we must treat general non-Hermitian block sparse matrices.
As a sparse matrix solver, the iterative method, which is often used to solve linear simultaneous equations with sparse coefficient matrices, is more efficient than the direct method such as the LU-factorization.
This algorithm is suitable for the parallel computing.
We developed a full-potential based screened KKR code which can treat such large super-cells by combining the iterative method for solving the Dyson equation with parallel computing techniques.
Through the test calculations, it is shown that our code is reliable and well parallelized w.r.t the number of atoms in the super-cell.
By the use of our code, it is possible to execute the all-electron level firstprinciples electronic structure calculations for super-cells which include thousands of atoms within a reasonable time.


研究活動

元素戦略拠点

触媒・電池元素戦略拠点
触媒・電池元素戦略研究拠点 (京都大学)
東工大元素戦略拠点
東工大元素戦略拠点 (東京工業大学)
構造材料元素戦略研究拠点
構造材料元素戦略研究拠点 (京都大学)
高効率モーター用磁性材料技術研究組合
高効率モーター用 磁性材料技術研究組合