局所密度近似のソースを表示

'''局所密度近似'''（きょくしょみつどきんじ、{{lang-en-short|Local Density Approximation}}、略称'''LDA'''）は、[[密度汎関数理論]]（DFT）における理論に現れる[[交換相関エネルギー|交換相関]]（XC）エネルギー[[汎関数]]に対する近似の一部類である。空間中の各点での{{仮リンク|電子密度|en| Electron density|preserve=1}}の値だけに依っている（密度の導関数や[[コーン–シャム方程式|コーン–シャム軌道]]には依存しない）。多くのアプローチによってXCエネルギーに対する局所近似を得ることができる。しかしながら、圧倒的に成功を収めている局所近似は[[ジェリウムモデル|均一電子ガス]]（HEG）モデルから導かれたものである。この点に関しては、LDAはHEG近似に基づく汎関数と一般的に同義である。

一般に、スピン非偏極系について、交換相関エネルギーに対する局所密度近似は次のような関数系を仮定する。
:<math>E_{\mathrm{xc}}^{\mathrm{LDA}}[\rho] = \int \rho(\boldsymbol{r})\epsilon_{\mathrm{xc}}(\rho(\boldsymbol{r}))\ \mathrm{d}\boldsymbol{r}\ </math>
上式において、''ρ''は[[電子密度]]、''ε''<sub>xc</sub>は電荷密度''ρ''を持つ均一電子ガスの粒子毎の交換相関エネルギーである。この仮定では空間の各点で（つまり局所的に）電子の'''交換・相関エネルギー密度'''<math>\epsilon_{\rm xc}</math>が決まっており、<math>\epsilon_{\rm xc}</math>はその場所の電子密度<math>n(\boldsymbol{r})</math>だけの関数になっている。この交換相関エネルギーは交換項と相関項に線形に分解される。
:<math>E_{\mathrm{xc}} = E_{\mathrm{x}} + E_{\mathrm{c}}</math>
こうすることで、''E''<sub>x</sub>と''E''<sub>c</sub>について別々の式を探すことができる。交換項はHEGに対して単純な解析形を取る。相関密度については限定的な式しか厳密に知られておらず、''ε''<sub>c</sub>に対する膨大な数の異なる近似が生み出された。

[[密度汎関数理論#ホーヘンベルク・コーンの定理|ホーヘンベルグ・コーンの定理]]によれば、この<math>E_{\rm xc}</math>は取り扱う系に依存しない普遍的な関数である。よって、もし局所密度近似が妥当であれば、<math>\epsilon_{\rm xc}</math>は（計算しやすい）一様電子系について求めた値でも、実際に計算したい系の値でも同じはずである。このようにして、一様電子系についてもとめた<math>\epsilon_{\rm xc}</math>を用いることが正当化され、実際の計算に用いることができる。

実際に用いられる<math>\epsilon_{\rm xc}</math>の関数形は、厳密に求められる低密度、高密度の極限からの[[外挿]]によるもの<ref name=wigner/><ref name="von_Barth"/><ref name="JanakMoruzzi1975">{{cite journal|last1=Janak|first1=J. F.|last2=Moruzzi|first2=V. L.|last3=Williams|first3=A. R.|title=Ground-state thermomechanical properties of some cubic elements in the local-density formalism|journal=Physical Review B|volume=12|issue=4|year=1975|pages=1257–1261|doi=10.1103/PhysRevB.12.1257}}</ref><ref name="GunnarssonLundqvist1976">{{cite journal|last1=Gunnarsson|first1=O.|last2=Lundqvist|first2=B. I.|title=Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism|journal=Physical Review B|volume=13|issue=10|year=1976|pages=4274–4298|doi=10.1103/PhysRevB.13.4274}}</ref><ref name="MacDonaldVosko1979">{{cite journal|last1=MacDonald|first1=A H|last2=Vosko|first2=S H|title=A relativistic density functional formalism|journal=Journal of Physics C: Solid State Physics|volume=12|issue=15|year=1979|pages=2977–2990|doi=10.1088/0022-3719/12/15/007}}</ref>や、[[モンテカルロ法]]を使ったもの<ref name="Ceperley1978">{{cite journal|last1=Ceperley|first1=D.|title=Ground state of the fermion one-component plasma: A Monte Carlo study in two and three dimensions|journal=Physical Review B|volume=18|issue=7|year=1978|pages=3126–3138|doi=10.1103/PhysRevB.18.3126}}</ref><ref name="CeperleyAlder1980">{{cite journal|last1=Ceperley|first1=D. M.|last2=Alder|first2=B. J.|title=Ground State of the Electron Gas by a Stochastic Method|journal=Physical Review Letters|volume=45|issue=7|year=1980|pages=566–569|doi=10.1103/PhysRevLett.45.566}}</ref><ref name="PerdewZunger1981">{{cite journal|last1=Perdew|first1=J. P.|last2=Zunger|first2=Alex|title=Self-interaction correction to density-functional approximations for many-electron systems|journal=Physical Review B|volume=23|issue=10|year=1981|pages=5048–5079|doi=10.1103/PhysRevB.23.5048}}</ref>などがある。

局所密度近似は、[[一般化勾配近似]]（GGA）や[[混成汎関数]]といった交換相関エネルギーに対するより洗練された近似の構築において重要である。これは、いかなる近似交換相関汎関数も均一電子ガスの厳密な結果を再現することが望まれるためである。こういったものとして、LDAはこういった汎関数の陽な混成要素としてしばしば取り入れられている。

== 応用 ==
局所密度近似はGGAと同様に[[固体物理学|固体物理学者]]によって半導体酸化物や[[スピントロニクス]]を含む半導体素材中の電子および磁気相互作用を解釈するためのDFT研究において広範に利用されている。これらの計算研究の重要性は、第一原理に基づく解析を必要とする合成パラメータに対する高い感受性を引き起す系の複雑さに由来する。ドープされた半導体酸化物中の[[フェルミ準位]]とバンド構造の予測はCASTEPやDMol3といったシミュレーションパッケージに取り入れられたLDAを使ってしばしば行われる<ref>{{cite journal| last1=Segall| first1=M.D.| last2=Lindan| first2=P.J | title= First-principles simulation: ideas, illustrations and the CASTEP code | journal= Journal of Physics: Condensed Matter | year= 2002| volume=14| issue=11| pages=2717|bibcode = 2002JPCM...14.2717S |doi = 10.1088/0953-8984/14/11/301 }}</ref>。しかしながら、LDAおよび[[一般化勾配近似|GGA]]としばしば関係している[[バンドギャップ]]値の過小評価は、こういった系における不純物媒介伝導性とキャリア媒介磁性の両方またはいずれか一方の誤った予測をもたらしうる<ref>{{cite journal| last1=Assadi| first1=M.H.N| title= Theoretical study on copper's energetics and magnetism in TiO<sub>2</sub> polymorphs| journal= Journal of Applied Physics | year=2013| volume=113| issue=23| pages= 233913–233913–5| doi=10.1063/1.4811539|arxiv = 1304.1854 |bibcode = 2013JAP...113w3913A | s2cid=94599250|display-authors=etal}}</ref>。1998年に始まった{{仮リンク|固有値についてのレイリーの定理|en|Rayleigh theorem for eigenvalue}}の応用によって、LDAポテンシャルを使って、材料のほとんど正確なバンドギャップの計算が可能となっている<ref>{{Cite journal|last1=Zhao|first1=G. L.|last2=Bagayoko|first2=D.|last3=Williams|first3=T. D.|date=1999-07-15|title=Local-density-approximation prediction of electronic properties of GaN, Si, C, andRuO2|journal=Physical Review B|volume=60|issue=3|pages=1563–1572|doi=10.1103/physrevb.60.1563|issn=0163-1829}}</ref><ref>{{Cite journal|last=Bagayoko|first=Diola|date=December 2014|title=Understanding density functional theory (DFT) and completing it in practice|journal=AIP Advances|volume=4|issue=12|pages=127104|doi=10.1063/1.4903408|issn=2158-3226|doi-access=free}}</ref>。[[密度汎関数理論#ホーヘンベルク・コーンの定理|DFTの第2定理]]に対する誤解は、LDAおよびGGA計算によるバンドギャップの過小評価の大半を説明するように思われる。

== 均一電子ガス ==
電子密度にのみ依存した''ε''<sub>xc</sub>に対する近似は数多くのやり方で開発することができる。最も成功を収めているやり方は均一電子ガス（HEG）に基づく。これは、相互作用のある''N''個の電子を、系を中性に保つ正の背景電荷を有する体積''V''に置くことによって構築される。''N''および''V''は次に、電子密度 (''ρ''&nbsp;=&nbsp;''N''&nbsp;/&nbsp;''V'') を有限に保つようなやり方で無限大まで持っていかれる。これは、全エネルギーが運動エネルギーおよび交換-相関エネルギーのみからの寄与によって構成され、波動関数が平面波の観点から表現できるため有用な近似である。具体的には、一定密度''ρ''に対して、交換エネルギー密度は''ρ''<sup>⅓</sup>に比例する。

== 交換汎関数 ==
HEGの交換エネルギー密度は解析的に知られている。交換に対するLDAは、密度が均一でない系における交換エネルギーがHEGの結果を各点に適用することによって得られるという近似の下でこの式を使用して、以下の式を得る<ref name="parryang">{{cite book|last=Parr|first=Robert G|author2=Yang, Weitao |title=Density-Functional Theory of Atoms and Molecules|publisher=Oxford University Press|location=Oxford |year=1994|isbn=978-0-19-509276-9}}</ref><ref>{{cite journal|last=Dirac|first=P. A. M.|year=1930|title=Note on exchange phenomena in the Thomas-Fermi atom|journal=Proc. Camb. Phil. Soc.|volume=26|pages=376–385|doi=10.1017/S0305004100016108|issue=3|bibcode = 1930PCPS...26..376D |doi-access=free}}</ref>。
:<math>E_{\mathrm{x}}^{\mathrm{LDA}}[\rho] = - \frac{3}{4}\left( \frac{3}{\pi} \right)^{1/3}\int\rho(\boldsymbol{r})^{4/3}\ \mathrm{d}\boldsymbol{r}</math>

== 相関汎関数 ==
HEGの相関エネルギーに対する解析表式は、それぞれ無限に弱い相関と無限に強い相関に対応する高密度および低密度限界で利用可能である。電子密度''ρ''を持つHEGについて、相関エネルギー密度の高密度限界<ref name="parryang"/>は
:<math>\epsilon_{\mathrm{c}} = A\ln(r_{s}) + B + r_{s}(C\ln(r_{s}) + D)</math>
であり、低密度限界は
:<math>\epsilon_{\mathrm{c}} = \frac{1}{2}\left(\frac{g_{0}}{r_{s}} + \frac{g_{1}}{r_{s}^{3/2}} + \dotsb\right)</math>
である。上式において、Wigner-Seitzパラメータ<math>r_{s}</math>は無次元である<ref name="Murray Gell-Mann and Keith A. Brueckner 1957 364">{{cite journal | title = Correlation Energy of an Electron Gas at High Density | author = Murray Gell-Mann and Keith A. Brueckner | journal = Phys. Rev. | volume = 106 | pages = 364–368 | year = 1957 | doi = 10.1103/PhysRev.106.364 | issue = 2| bibcode = 1957PhRv..106..364G | url = https://authors.library.caltech.edu/3713/1/GELpr57b.pdf }}</ref>。これは、厳密に1つの電子を包含する球の半径をボーア半径で割った値として定義される。Wigner-Seitzパラメータ<math>r_s</math>は密度と以下の式で結び付けられる。
:<math>\frac{4}{3}\pi r_{s}^{3} = \frac{1}{\rho}</math>

密度の全領域に対する解析表式は多体摂動論に基づいて提案されてきた。計算された相関エネルギーは2ミリハートリー以内で[[量子モンテカルロ法|量子モンテカルロ]]シミュレーションの結果と一致する。

HEGのエネルギーに対する精密な量子モンテカルロシミュレーションは複数の中間的値の密度について実行され、次々に相関エネルギー密度の精密な値を与えてきた<ref>{{cite journal | title = Ground State of the Electron Gas by a Stochastic Method | author = D. M. Ceperley and B. J. Alder | journal = Phys. Rev. Lett. | volume = 45 | pages = 566–569 | year = 1980 | doi = 10.1103/PhysRevLett.45.566 | bibcode=1980PhRvL..45..566C | issue = 7| url = https://digital.library.unt.edu/ark:/67531/metadc1059358/ }}</ref>。相関エネルギー密度に対する最も人気のあるLDAは、厳密に知られている漸近挙動を再現しながら、シミュレーションから得られたこれらの正確な値を内挿する。''ε''<sub>c</sub>に対する異なる解析形式を使った様々なアプローチによって相関汎関数に対する複数のLDAが生み出されてきた。
* Vosko-Wilk-Nusair (VWN)<ref name="vwn">{{cite journal | title = Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis | author = S. H. Vosko, L. Wilk and M. Nusair | journal = Can. J. Phys. | volume = 58 | pages = 1200–1211 |  year = 1980 | doi = 10.1139/p80-159 |bibcode = 1980CaJPh..58.1200V | issue = 8 | url = https://escholarship.org/content/qt23j4q7zm/qt23j4q7zm.pdf?t=obc5l4 }}</ref>
* Perdew-Zunger (PZ81)<ref name="pz81">{{cite journal | title = Self-interaction correction to density-functional approximations for many-electron systems | author = J. P. Perdew and A. Zunger | journal = Phys. Rev. B | volume = 23 | pages =  5048–5079 | year = 1981 | doi = 10.1103/PhysRevB.23.5048 |bibcode = 1981PhRvB..23.5048P | issue = 10 | doi-access = free }}</ref>
* Cole-Perdew (CP)<ref>{{cite journal | title = Calculated electron affinities of the elements | author = L. A. Cole and J. P. Perdew | journal =  Phys. Rev. A | volume = 25 | pages =  1265–1271  | year = 1982 | doi = 10.1103/PhysRevA.25.1265 |bibcode = 1982PhRvA..25.1265C | issue = 3 }}</ref>
* Perdew-Wang (PW92)<ref name=pw92>{{cite journal | title = Accurate and simple analytic representation of the electron-gas correlation energy | author = John P. Perdew and Yue Wang | journal = Phys. Rev. B | volume = 45 | pages = 13244–13249 | year = 1992 | doi = 10.1103/PhysRevB.45.13244 | pmid = 10001404 |bibcode = 1992PhRvB..4513244P | issue = 23 }}</ref>

これらや、DFTそれ自身の形式的樹立よりさえも前から存在するのがHEGモデルから[[メラー＝プレセット法|摂動論的]]に得られるWigner相関汎関数である<ref name=wigner>{{cite journal | title = On the Interaction of Electrons in Metals | author = E. Wigner | journal = Phys. Rev. | volume = 46 | pages = 1002–1011 | year = 1934 | doi = 10.1103/PhysRev.46.1002 |bibcode = 1934PhRv...46.1002W | issue = 11 }}</ref>。

== スピン偏極 ==
[[スピン偏極]]系への密度汎関数の拡張は、厳密なスピンスケーリングが知られている交換については明快であるが、相関についてはさらなる近似が用いられなければならない。DFTにおけるスピン偏極系は2つのスピン密度''ρ''<sub>α</sub>および''ρ''<sub>β</sub>（ ''ρ''&nbsp;=&nbsp;''ρ''<sub>α</sub>&nbsp;+&nbsp;''ρ''<sub>β</sub>）を用い、局所スピン密度近似（Local Spin Density Approximation, LSDA）の形式は
:<math>E_{\mathrm{xc}}^{\mathrm{LSDA}}[\rho_{\alpha},\rho_{\beta}] = \int\mathrm{d}\boldsymbol{r}\ \rho(\boldsymbol{r})\epsilon_{\mathrm{xc}}(\rho_{\alpha},\rho_{\beta})</math>
である。LSDAは[[バンド計算]]において[[磁性]]（[[強磁性]]、[[反強磁性]]、[[ハーフメタリック]]など）やスピンの問題（スピン分極）を扱う時に使用される。

交換エネルギーについては、（局所密度近似に対してのみではない）厳密な結果がスピン非偏極汎関数の観点から知られている<ref>{{cite journal|last=Oliver|first=G. L.|author2=Perdew, J. P.  |year=1979|title=Spin-density gradient expansion for the kinetic energy|journal=Phys. Rev. A|volume=20|pages=397–403|doi=10.1103/PhysRevA.20.397|bibcode = 1979PhRvA..20..397O|issue=2 }}</ref>。
:<math>E_{\mathrm{x}}[\rho_{\alpha},\rho_{\beta}] = \frac{1}{2}\bigg( E_{\mathrm{x}}[2\rho_{\alpha}] + E_{\mathrm{x}}[2\rho_{\beta}] \bigg)\ .</math>

相関エネルギー密度のスピン依存性は相対スピン偏極度
:<math>\zeta(\boldsymbol{r}) = \frac{\rho_{\alpha}(\boldsymbol{r})-\rho_{\beta}(\boldsymbol{r})}{\rho_{\alpha}(\boldsymbol{r})+\rho_{\beta}(\boldsymbol{r})}</math>
を導入することによってアプローチする。<math>\zeta = 0\,</math>は等しい<math>\alpha\,</math>および<math>\beta\,</math>スピン密度を持つ常磁性スピン非偏極状況に対応しするが、<math>\zeta = \pm 1</math>は一方のスピン密度が消滅する強磁性状況に対応する。全密度および相対偏極度の所与の値に対するスピン相関エネルギー密度''ε''<sub>c</sub>(''ρ'',''ς'') は極値を内挿するように構築される。いくつかの形式がLDA相関汎関数と共に開発されてきた<ref name="vwn"/><ref name="von_Barth">{{cite journal|last=von Barth|first=U.|author2=Hedin, L. |year=1972|title=A local exchange-correlation potential for the spin polarized case|journal=J. Phys. C: Solid State Phys.|volume=5|pages=1629–1642|doi=10.1088/0022-3719/5/13/012|bibcode = 1972JPhC....5.1629V|issue=13 }}</ref>。

== 実例となる計算 ==
LDA計算は実験値とまあまあの一致を示す。
{| class="wikitable"
|+ イオン化ポテンシャル (eV) <ref name="parryang"/>
|-
! 
! LSDA
! LDA
! [[ハートリー＝フォック方程式|HF]]
! 実験値
|- align="center"
|  H || 13.4 || 12.0 ||    13.6  || 13.6
|- align="center"
| He || 24.5 || 26.4 ||      || 24.6
|- align="center"
| Li || 5.7  || 5.4  || 5.3  || 5.4
|- align="center"
| Be || 9.1  ||      || 8.0  || 9.3
|- align="center"
| B  || 8.8  ||      || 7.9  || 8.3
|- align="center"
| C  || 12.1 ||      || 10.8 || 11.3
|- align="center"
| N  || 15.3 ||      || 14.0 || 14.5
|- align="center"
| O  || 14.2 || 16.5 || 11.9 || 13.6
|- align="center"
| F || 18.4 ||       || 16.2 || 17.4
|- align="center"
| Ne || 22.6 || 22.5 || 19.8 || 21.6
|}

{| class="wikitable"
|+ 結合長（オングストローム）<ref name="parryang"/>
|-
!
! 実験値
! LSDA
! 誤差
|- align="center"
| H<sub>2</sub> || 0.74 || 0.77 || 0.03
|- align="center"
| Li<sub>2</sub>  || 2.67 || 2.71 || 0.04
|- align="center"
| B<sub>2</sub> || 1.59 || 1.60 || 0.02
|- align="center"
| C<sub>2</sub> || 1.24 || 1.24 || 0.00
|- align="center"
| N<sub>2</sub> || 1.10	|| 1.10 || 0.00
|- align="center"
| O<sub>2</sub> || 1.21 || 1.20 || 0.01
|- align="center"
| F<sub>2</sub> || 1.42 || 1.38 || 0.04
|- align="center"
| Na<sub>2</sub> || 3.08 || 3.00 || 0.08
|- align="center"
| Al<sub>2</sub> || 2.47 || 2.46 || 0.01
|- align="center"
| Si<sub>2</sub> || 2.24 || 2.27 || 0.03
|- align="center"
| P<sub>2</sub> || 1.89	|| 1.89 || 0.01
|- align="center"
| S<sub>2</sub> ||  1.89 || 1.89 || 0.00
|- align="center"
| Cl<sub>2</sub> || 1.99 || 1.98 || 0.01
|- align="center"
| Average ||  ||  || 0.02
|}

== 交換-相関ポテンシャル ==
局所密度近似に対する交換-相関エネルギーに対応する交換-相関ポテンシャルは以下の式で与えられる<ref name="parryang"/>。
:<math>v_{\mathrm{xc}}^{\mathrm{LDA}}(\boldsymbol{r}) = \frac{\delta E^{\mathrm{LDA}}}{\delta\rho(\boldsymbol{r})} = \epsilon_{\mathrm{xc}}(\rho(\boldsymbol{r})) + \rho(\boldsymbol{r})\frac{\partial \epsilon_{\mathrm{xc}}(\rho(\boldsymbol{r}))}{\partial\rho(\boldsymbol{r})}</math>

有限の系においては、LDAポテンシャルは指数関数的な形で漸近的に減衰する。これは誤りである。真の交換-相関ポテンシャルはクーロン的によりゆっくりと減衰する。人為的に急速な減衰は、ポテンシャルが束縛できるコーン・シャム軌道の数（つまり、ゼロ未満のエネルギーを持つ軌道の数）に現れる。LDAポテンシャルはリュードベリ系列を支持できず、ポテンシャルが束縛するそれらの状態はエネルギーが高過ぎる。これはエネルギー的に高過ぎる[[HOMO]]エネルギーをもたらし、[[クープマンズの定理]]に基づく[[イオン化ポテンシャル]]に対する予測は精度が低い。そのうえ、LDAは[[陰イオン]]といった電子豊富種のまずい描写を与える。こういった場合、LDAはしばしば追加の電子を束縛することができず、陰イオン種が不安定であると誤って予測する<ref name="pz81"/><ref>{{cite book|last=Fiolhais|first=Carlos|author2=Nogueira, Fernando |author3=Marques Miguel |title=A Primer in Density Functional Theory|publisher=Springer|year=2003|isbn=978-3-540-03083-6|page=60}}</ref>。

== LDAを越える試み==
LDAを越える試みとは、局所密度近似 (LDA) の問題点を解消する新たな手法を見出す試みの総称である。

局所密度近似は大変成功した近似であるが、実際の系に対する様々な計算の結果、その限界もまた露わになってきた。代表的な問題点とその克服に向けたアプローチについて記述する。

=== 代表的な問題点（限界） ===
#[[半導体]]、[[絶縁体]]において[[バンドギャップ]]が実験値より過小な値となる。
#[[鉄]]の[[強磁性]]結晶構造（体心立方構造：BCC）が安定とならない。（他にも安定構造や電子状態がLDAが原因で、実際のものと一致しない場合がある。GGA近似を行うことで修正される場合がある。）。
#[[活性化エネルギー]]の過小評価。
#鏡像ポテンシャルが記述できない（[[表面]]）。
#[[自己相互作用補正]]の問題。
#絶対零度（基底状態）での計算が前提。←[[密度汎関数理論]]
#励起状態に対する計算の正しさの保証がない（これは、むしろ密度汎関数理論の問題）。など

=== 問題を克服する手段・手法 ===
以下のようなものが提案、試行されている。
*[[一般化勾配近似|GGA]]（Generalized Gradient Approximation, 一般化された密度勾配近似）
*SIC（Self-Interaction Correction, [[自己相互作用補正]]）
*[[GW近似]]
*[[LDA plus U|LDA+U]] (LSDA+U)
*[[TDDFT]] ([[TDLDA]])（Time-Dependent DFT, 時間発展を考慮した密度汎関数理論）
さらに、交換項を（[[ハートリー-フォック法]]での交換項として）厳密に取り扱うアプローチ (Exact Exchange)、[[密度汎関数理論]]の[[有限温度への拡張]]や、[[電子]]の[[多体問題]]をより直接的に扱う方法（[[量子モンテカルロ法]]による）、また動的平均場法などの強相関電子系でのモデル計算で開発された手法と組み合わせ、電子相関の効果を導入する研究がされているが、まだ汎用的な計算手法とは言い難く、簡単な系でのテスト計算どまりである。

== 脚注 ==
{{reflist|2}}

==参考文献==
* {{cite journal | title=Self-Consistent Equations Including Exchange and Correlation Effects | journal=Physical Review | year=1965 | author=W. Kohn |coauthors=L. J. Sham | volume=140 | issue=4A | pages=A1133-1138| doi=10.1103/PhysRev.140.A1133}}

==関連項目==
*[[電子相関]]
*[[Xα法]]
*[[密度汎関数法]]
*[[第一原理バンド計算]]

{{DEFAULTSORT:きよくしよみつときんし}}
[[Category:密度汎関数理論]]