状態密度のソースを表示

{{翻訳直後|[[:en:Special:Redirect/revision/767410454|en:Density of states (19:29, 25 February 2017 UTC)]]|date=2017年5月}}

[[固体物理学]]および[[物性物理学]]において、系の'''状態密度'''（じょうたいみつど、{{Lang-en-short|density of states}}, '''DOS'''）とは、微小なエネルギー区間内に存在する、系の占有しうる状態数を各エネルギーごとに記述する物理量である。[[気相]]中の[[原子]]や[[分子]]のような{{仮リンク|孤立系|en|Isolated_system}}とは異なり、[[確率密度関数|密度分布]]は[[スペクトル密度]]のような[[離散]]分布ではなく[[連続]]分布となる。あるエネルギー準位において DOS が高いことは、そこに占有しうる状態が多いことを意味する。DOS がゼロとなることは、系がそのエネルギー準位を占有しえないことを意味する。一般的に DOS とは、空間的および時間的に平均されたものを言う。局所的な変動は[[局所状態密度]] ({{Lang|en|local density of states}}, LDOS) と呼ばれ区別される。

== 概要 ==
量子系において、波もしくは波動的粒子は系によって定まる波長と伝播方向をもつモードもしくは状態を占める。特定の状態のみが許容されることも多く、系によっては物質の原子間距離と原子核電荷により特定の波長の電子のみが存在を許容される場合もある。また、物質の結晶構造により波が一方向にのみ伝播を許容され、別の方向への伝播が抑制されるような系もある。したがって、特定の波長においては多くの状態が許容され、別のエネルギー準位には全く許容される状態が存在しないということがありうる。量子系により、[[電子]]や[[光子]]、[[フォノン]]の状態密度をエネルギーもしくは[[波数ベクトル]] {{mvar|k}} の関数として計算することができる。DOS を表わす記号としては、{{mvar|g, ρ, D, n, N}} などが用いられる。エネルギーの関数としての DOS と波数ベクトルの関数としての DOS の間の変換は、系ごとに決まる {{mvar|E}} と {{mvar|k}} との間の[[分散関係]]が分かっていれば行うことができる。

たとえば、半導体中の電子の状態密度は伝導帯端においては低く、電子の占有できる状態は少ない。電子のエネルギーが増えるにつれて状態密度も増加し、占有できる状態が増える。しかし、バンドギャップ中には電子の占有できる状態は存在しないため、伝導帯端の電子は別のモードへと遷移するために少くとも {{math|''E''<sub>g</sub>}} だけのエネルギーを失う必要がある。

一般的に、系の[[位相幾何学]]的性質が状態密度の主な性質を決定する。[[中性子星]]中の[[ニュートロニウム]]や[[金属]]中の自由電子ガス（[[フェルミ縮退|縮退物質]]と[[フェルミ気体]]の例）のような最も良く知られた系では、3次元[[ユークリッド空間]]の位相構造を持つ。より知られていない系としては、[[グラファイト]]層中の[[二次元電子ガス]]や [[MOSFET]] 型素子中の[[量子ホール効果]]系は二次元ユークリッド空間の位相構造を持つ。さらに知られていないものとしては、[[カーボンナノチューブ]]や{{仮リンク|量子ワイヤ|en|Quantum_wire}}、[[朝永–ラッティンジャー液体|朝永・ラッティンジャー液体]]などは1次元位相幾何を持つ。1D および 2D 位相幾何を持つ系は、[[ナノテクノロジー]]と[[物質科学]]の進展につれてよりよく知られるようになると考えられる。

== 定義 ==
系が状態 i をとるときの系のエネルギーが <math>\mathcal{E}_i</math> で与えられるとき、状態密度は
:<math>D(E) =\sum_i \delta(E-\mathcal{E}_i)</math>
で定義される。

系の状態が連続パラメータ {{Mvar|λ}} で指定される場合の定義は
:<math>D(E) =\int \delta(E-\mathcal{E}(\lambda))\, \mu(\lambda)\, \mathrm d^f\lambda</math>

である。ここで {{Mvar|μ}} はパラメータ {{Mvar|λ}} の張る状態空間の体積を与える[[測度]]である。

古典系の状態は正準変数の組 {{Math|(''p'', ''q'')}} で指定され、系のエネルギーは[[ハミルトン関数]]で与えられる。正準変数の張る[[位相空間 (物理学)|位相空間]]の体積は、一対の {{Math|d''p''d''q''}} ごとに[[プランク定数]] {{Mvar|h}} で割る約束で、自由度が {{Math|2''f''}} での状態密度は
:<math>D(E) =\frac{1}{h^f} \int \delta(E-\mathcal{H}(p,q))\, \mathrm d^fp\, \mathrm d^fq</math>
となる。

状態密度をエネルギー {{Mvar|E}} まで積分すれば
:<math>N(E) =\int^E_{-\infty} D(E')\, \mathrm dE' =\sum_i \int^E_{-\infty} \delta(E'-\mathcal{E}_i)\, \mathrm dE'</math>
として、系のエネルギーが <math>\mathcal{E}_i < E</math> である[[状態数|状態の数]]を与える。状態数が微分可能である場合には、状態密度は
:<math>D(E) =\frac{\mathrm dN}{\mathrm dE}</math>
で与えられる。エネルギーが E である状態の[[縮退度]]は
:<math> g(E) = \lim _{\Delta E\to 0} \int^E_{E-\Delta E} D(E')\, \mathrm dE' = \lim _{\Delta E\to 0} D(E)\, \Delta E</math>
で与えられる。ここで最後の等式は、積分の[[平均値の定理]]が妥当なときのみ成り立つ。

== 対称性と状態密度 ==
[[ファイル:Brillouin_Zone_(1st,_FCC).svg|サムネイル|260x260ピクセル|[[面心立方格子]]における第一[[ブリルアンゾーン|ブリュアンゾーン]]。[[切頂八面体]]となっている。高対称性線および点に対称性ラベルを付記してある。]]
DOS 計算の行える系は多種多様である。凝縮系において重要な性質は系の[[結晶構造|微視的構造]]のもつ[[対称性]]である。 [[流体]]、[[ガラス]]、[[アモルファス]]固体は回転対称性のある分散関係を持つ。球対称な系では、たとえば関数の積分などは一次元となる。なぜなら、計算が分散関係の動径パラメータにのみ依存するからである。

例えば[[単結晶]]からなる系などの[[等方的と異方的|非等方]]な系においては、状態密度がある結晶学的方位と別の方位とでは異るので、角度に依存する計算および計測が必要となる。非等方な問題は計算が難しくなり、また非等方な状態密度は可視化するのも難しくなる。そのため、ある特定の点のみを計算したり、射影状態密度 ({{Lang|en|projected density of states}}, PDOS) を計算したりといった手法がよく用いられる。
[[ファイル:Dual_Cube-Octahedron.svg|左|サムネイル|151x151ピクセル|八面体]]
粉末試料や多結晶試料に対する測定には、系の分散関係の定義域（一つの[[ブリュアンゾーン]]とすることが多い）全体にわたる積分が必要となる。系の対称性が高い場合、系の分散関係を表わす関数の形は分散関係の定義域全体にわたって何度も繰り返し表われる。このような場合、DOS の計算は還元ゾーンのみについての計算に帰着し、相当に省力化できる<ref name="Harrison">{{cite book|author=Walter Ashley Harrison|title=Electronic Structure and the Properties of Solids|url=https://books.google.com/books?id=R2VqQgAACAAJ|year=1989|publisher=Dover Publications|isbn=0-486-66021-4}}</ref>。面心立方格子のブリュアンゾーンは[[点群]] ''O<sub>h</sub>'' の、完全{{仮リンク|八面体対称性|en|Octahedral_symmetry}}をもつ48重対称性を持つ。よって、ブリュアンゾーン全体にわたる積分をその48分の1の部分領域にわたる積分に帰着することができる。{{仮リンク|結晶構造の周期表|en|Periodic_table_(crystal_structure)}}に示される通り、面心立方格子をとる元素は[[ダイヤモンド]]、[[ケイ素|シリコン]]、[[白金]]など多く、これらのブリュアンゾーンおよび分散関係は48重対称性を持つ。
[[ファイル:Rhombicdodecahedron.gif|サムネイル|256x256ピクセル|[[立方晶系|体心立方格子]]の第一ブリュアンゾーン。 [[菱形十二面体]]となる。]]
良く知られている結晶構造として、体心立方格子と六方最密充填格子の二つが挙げられる。体心立方格子は点群 ''T<sub>h</sub>'' の24重{{仮リンク|label=黄鉄鉱型対称性|四面体対称性|en|Tetrahedral_symmetry}}を持つ。六方最密充填格子は点群 ''D<sub>3h</sub>'' の12重{{仮リンク|label=プリズム二面体対称性|二面体対称性|en|Dihedral_symmetry_in_three_dimensions}}を持つ。点群の対称性の特性の網羅的リストについては、{{仮リンク|点群指標表|en|List_of_character_tables_for_chemically_important_3D_point_groups#Character_tables}}を参照のこと。

一般に、対称性が高く位相的次元の低い系ほど DOS の計算は容易である。回転対称性のある分散関係の状態密度は解析的に計算可能であることが多い。鋼やシリコンなど、実用上の興味の対象となる物質は高い対称性を持っていることが多いので、このことは幸運である。

== 波数空間の位相幾何 ==
[[ファイル:K-space.JPG|右|サムネイル|250x250ピクセル|図1: 三次元自由度を持つ電子の波数空間における球面]]
状態密度は対象の次元に依存する。次元の果たす役割は、DOS の単位 (Energy<sup>−1</sup>Volume<sup>−1</sup>) からも明らかである。系が二次元的になる極限において体積は面積となり、一次元的となる極限においては長さとなる。ここでいう体積とは波数空間上の、分散関係から導かれる{{仮リンク|等エネルギー面|en|Constant_energy_surface}}で囲われる領域の体積であることに注意が必要である。固体中の電子の分散関係は[[バンド構造]]を成している。三次元的波数空間の例を図1に示す。系の次元そのものが系内の粒子の運動量を規定することが見てとれる。

=== 波数ベクトル状態密度（球） ===
DOS を計算するにはまずある {{mvar|k}} に対して波数空間上の領域 {{math|[''k'', ''k''+d''k'']}} 内に含まれる状態数 {{Mvar|N}} を数える必要がある。これは、ある {{Mvar|k}} に対する {{Mvar|n}} 次元波数空間全体の体積 {{Math|&Omega;<sub>''n'', ''k''</sub>}} を {{Mvar|k}} で微分することで得られる。三次元、二次元、一次元波数空間の体積、面積、長さは次のように表わされる。
: <math>\Omega_n(k) = c_n k^n</math>
ここで、 {{Math|''c''<sub>''n''</sub>}} は波数空間の次元 {{Mvar|n}} に依存して位相幾何学的に定まる定数で、一次元、二次元、三次元ユークリッド波数空間に対してはそれぞれ以下のように定まる。
: <math>c_1 = 2,\quad c_2 = \pi,\quad c_3 = \frac{4 \pi}{3}</math>

この式によれば、波数ベクトル状態密度 {{Mvar|N}} は {{Math|&Omega;<sub>''n'', ''k''</sub>}} を {{Mvar|k}} で微分することにより次のように得られる。
: <math>N_n(k) = \frac{\mathrm d\Omega_n(k)}{\mathrm dk} = n\, c_n\, k^{(n - 1)} </math>
これを一次元、二次元、三次元の場合に明示的に書き下すと次のようになる。
: <math>N_1(k) = 2</math>
: <math> N_2(k) = 2 \pi k </math>
: <math> N_3(k) = 4 \pi k^2 </math>
一つの状態は波長 {{Mvar|λ}}J の粒子を含むことができる程度に大きい。波長と波数 {{Mvar|k}} との間の関係式は以下のようになる。
: <math>k = \frac{2\pi}{\lambda}</math>
長さ {{Mvar|λ}} の量子系は粒子を閉じ込める系の大きさ {{Mvar|L}} に依存する。最後に、状態密度 {{Mvar|N}} に係数 {{Math|''s''/''V''<sub>''k''</sub>}} をかける。ここで、{{Mvar|s}} はスピンや偏極などの物理現象に起因する内部自由度である。このような物理現象が無い場合は {{Math|1=''s''=1}} となる。{{Math|''V''<sub>''k''</sub>}} は波数空間上の、ある {{mvar|k}} よりも小さい波数ベクトルを全て含む体積である。

=== エネルギー状態密度 ===
DOS の計算の最後として、あるエネルギー <math>E</math> に対して 定まる区間 {{Math|[''E'', ''E''+d''E'']}} に含まれる体積あたりの状態数を計算する。一般的な系の DOS は次のような形式となる。
: <math>D_n\left(E\right) = \frac{\mathrm d\Omega_n(E)}{\mathrm dE}</math>
ここまでの式は、分散関係が[[単調増加]]する球対称な系に対してのみ成り立つ。 一般に、分散関係 {{Math|''E''(''k'')}} は球対称ではなく、単調増加でもないことが多い。{{Mvar|D}} を {{Mvar|E}} の関数として分散関係 {{Math|''E''(''k'')}} の[[逆関数]]を用いてここまでの式中に現われていた {{Mvar|k}} の関数 {{Math|&Omega;<sub>''n''</sub>(''k'')}} をエネルギーの関数 {{Math|&Omega;<sub>''n''</sub>(''E'')}} に置き換える必要がある。これは分散関係が球対称でなかったり単調増加しなかったりする場合は容易ではなく、ほとんどの場合において DOS は数値的に計算される。より詳細な導出もある<ref>[http://britneyspears.ac/physics/dos/dos.htm Sample density of states calculation]</ref><ref>[http://ece-www.colorado.edu/~bart/book/book/chapter2/ch2_4.htm Another density of states calculation]</ref>。

== 分散関係 ==
粒子の[[運動エネルギー]]は波数ベクトル {{Mvar|k}} の大きさと向きに依存する。たとえば[[フェルミ気体]]中の[[電子]]の運動エネルギーは以下のように得られる。
: <math>E = E_0 + \frac{(\hbar k)^2}{2m}</math>
ここで {{Mvar|m}} は電子質量である。この分散関係は球対称かつ単調増加であるから、DOS を容易に計算することができる。

[[ファイル:Monoatomic_chain_phonon_dispersion.svg|右|サムネイル|300x300ピクセル|図2:単原子鎖フォノンの分散関係]]
原子鎖の縦モード[[フォノン]]の分散関係は、図2に示すような 1 次元 {{Mvar|k}} 空間上の運動エネルギーについての関数となり、数式で表わすと以下のようになる。
: <math>E =  2 \hbar \omega_0 |\sin(ka/2)| </math>
ここで <math>\omega_0 = \sqrt{k_F / m}</math> は振動子周波数、{{Mvar|m}} は原子の質量、{{Math|''k''<sub>F</sub>}} は原子間に働く力の力定数、{{Mvar|a}} は原子間距離である。力定数が小さく、{{Math|''k'' ≪ ''&pi;'' / ''a''}} が満たされるような値である場合は分散関係は線形となる。
: <math>E = \hbar \omega_0 ka </math>
{{Math|''k'' &asymp; ''&pi;'' / ''a''}} の場合は以下のようになる。
: <math>E = 2 \hbar \omega_0 |\cos(\pi/2 - ka/2)| </math>
変数変換 {{Math|1=''q'' = ''k'' &minus; ''&pi;''/''a''}} を施して {{Mvar|q}} が小さくなるとき、分散関係は以下のように書ける。
: <math>E = 2 \hbar \omega_0 [1 - (qa/2)^2] </math>

=== 等方的分散関係 ===
ここで言及した二つの例は次のように書ける。
: <math>E = E_0 + c_k k^p </math>
この種の[[分散関係]]はエネルギーが波数ベクトルの長さのみに依存し、向きに依存しないため[[等方性媒質|等方]]的な分散関係といえる。このとき逆に、波数ベクトルの大きさはエネルギーを用いて以下のように書ける。
: <math> k = \left(\frac{E-E_0}{c_k}\right)^{1/p} </math>
また、{{Mvar|k}} よりも小さい波数ベクトルを含む {{Mvar|n}} 次元 {{Mvar|k}} 空間上の体積は次のように書ける。
: <math>\Omega_n(k) = c_n k^n</math>
したがって、等方的分散関係から、被占有状態の体積は以下のように書ける。
: <math>\Omega_n(E) = \frac{c_n}{c_k^{n/p}}\left(E-E_0\right)^{n/p}\ ,</math>
この体積をエネルギーで微分すれば等方的分散関係に対する DOS を得ることができる。
: <math>D_n\left(E\right) = \frac {\mathrm d}{\mathrm dE}\Omega_n(E) = \frac{n c_n}{p c_k^{n/p}}\left(E-E_0\right)^{(n/p - 1)} </math>

=== 放物線分散関係 ===
[[ファイル:Free-electron_DOS.svg|右|サムネイル|300x300ピクセル|Figure 3: Free-electron DOS in 3-dimensional k-space]]
フェルミ気体中の自由電子などのように分散関係が放物線を描く ({{Math|1=''p'' = 2}}) 場合、{{Mvar|n}} 次元系における状態密度 <math>D_n\left(E\right)</math> は以下のようになる。
: <math>D_1\left(E\right)  = \frac {1}{ \sqrt{c_k(E-E_0)}}</math>

: <math>D_2\left(E\right) = \frac {\pi}{c_k}</math>

: <math>D_3\left(E\right) = 2 \pi \sqrt{\frac{E-E_0}{c_k^3}} </math>
ここで <math>E > E_0</math> とし、<math>D(E) = 0</math> と <math>E < E_0</math> の場合はする。

1 次元系では DOS は {{Mvar|E}} が {{Math|''E''<sub>0</sub>}} に落ちる際に発散する。2 次元系では {{Mvar|E}} に依存しなくなる。3 次元系では状態密度はエネルギーの平方根に比例して増加する<ref name="Kittel">{{cite book|author=Charles Kittel|title=Introduction to Solid State Physics|edition=7th|year=1996|publisher=Wiley|page=Equation (37), p. 216|nopp=true|isbn=0-471-11181-3}}</ref>。

係数部分を全て書き下すと、3 次元系における DOS は以下のように書ける。
: <math>N(E) = \frac {V}{2\pi^2} \left(\frac {2m}{\hbar^2}\right)^{3/2}\sqrt{E-E_0}</math>
ここで {{Mvar|V}} は総体積であり、{{Math|''N''(''E''&minus;''E''<sub>0</sub>)}} には2重のスピン縮退を含む。

=== 線形分散関係 ===
[[フォノン|光子]]や[[フォノン|音響フォノン]]、特定の固体中の電子バンドのように分散関係が線形 ({{Math|1=''p'' = 1}}) のとき、1、2、3 次元系におけるエネルギーに対する DOS はそれぞれ以下のようになる。
: <math>D_1\left(E\right) = \frac{1}{c_k} </math>
: <math>D_2\left(E\right) = \frac{2 \pi}{c_k^2}\left(E-E_0\right) </math>
: <math>D_3\left(E\right) = \frac{4 \pi}{c_k^3}\left(E-E_0\right)^2 </math>

== 状態密度と分布関数 ==
状態密度は[[フォノン|固体中の運動エネルギー理論]]において重要な役割を果たす。状態密度と[[確率密度関数|確率密度分布]]との積は熱平衡状態にある系について、あるエネルギーにおける単位体積あたりの被占有状態数を与える。この値は物質の様々な物性を調べる際に広く用いられている。ここで、確率密度分布と状態密度からどのように物性を得るかの例をいくつか挙げる。
[[ファイル:Fermi-Dirac_product_with_Density_of_States.svg|サムネイル|320x320ピクセル|図4: ある半導体におけるフェルミ・ディラック分布（青色）、状態密度（橙色）、それらの積（緑色）。緑の線のふくらみのうち下の方は[[正孔]]のエネルギーであるから、確率分布関数としては {{Math|1 &minus; ''f''(''x'')}} を用いてある。]]
[[フェルミ分布関数|フェルミ・ディラック統計]]: 図4に示すフェルミ・ディラック分布は、熱平衡状態においてフェルミオンが特定の量子状態を占有する確率を与える。[[フェルミ粒子|フェルミオン]]は[[パウリの排他原理|パウリの排他律]]に従う粒子であり、例えば電子、陽子、中性子などが挙げられる。この分布関数は次のように書ける。
: <math>f_{\mathrm{FD}}(E) = \frac{1}{\exp\left(\frac{E-\mu}{k_\mathrm{B} T}\right)+1}</math>
{{Mvar|&mu;}} は[[化学ポテンシャル]]（{{Math|1=''T'' = 0}} の場合[[フェルミ準位]]と呼び {{Math|''E''<sub>F</sub>}}  と書く）、{{Math|''k''<sub>B</sub>}} は[[ボルツマン定数]]、{{Mvar|T}} は[[温度]]である。図4に示す、フェルミ・ディラック分布関数と3次元半導体の状態密度の積がキャリア密度やエネルギーバンドギャップなどの物性についての知識を得るために用いられる。

[[ボース分布関数|ボース・アインシュタイン統計]]: ボース・アインシュタイン分布関数は熱平衡にある系においてボソンがある量子状態を占有する確率を表わす。[[ボース粒子|ボソン]]はパウリの排他律に従わない粒子で、例えばフォノンや光子が挙げられる。この分布関数は以下のように書ける。
: <math>f_{\mathrm{BE}}(E) = \frac{1}{\exp\left(\frac{E-\mu}{k_BT}\right)-1}</math>
これら二つの分布関数から、[[内部エネルギー]] {{Mvar|U}}、粒子数 {{Mvar|n}}、[[熱容量|比熱容量]] {{Mvar|C}}、[[熱伝導率]]  {{Mvar|k}} を計算することができる。 これらの物性値と、密度関数と分布関数との関係式は、状態密度を {{Math|''D''(''E'')}} ではなく {{Math|''g''(''E'')}} と書くと、以下のようになる。
: <math>U = \int E\, f(E)\, g(E)\,\mathrm dE</math>
: <math>n = \int f(E)\, g(E)\,\mathrm dE</math>
: <math>C = \frac{\partial}{\partial T}\int E\, f(E)\, g(E) \,\mathrm dE</math>
: <math>k = \frac{1}{d}\frac{\partial}{\partial T}\int E f(E)\, g(E)\, \nu(E)\, \Lambda(E)\,\mathrm dE</math>
{{Mvar|d}} は次元数、{{Mvar|&nu;}} は音速、{{Math|&Lambda;}} は[[平均自由行程]]である。

== 応用 ==
状態密度は物理学の多くの分野で登場し、量子力学的現象の説明の助けとなる。

=== 量子化 ===
微視的構造に対して状態密度を計算すると、次元が減るにつれて電子の分布が変化することがわかる。特定のエネルギー領域において{{仮リンク|量子ワイヤー|en|Quantum_wires}}の DOS は、バルク半導体の DOS に比べて実際に高くなり、[[量子ドット]]の DOS は特定のエネルギーに量子化される。

=== フォトニック結晶 ===
光の波長スケールの繰り返し構造を用いると光子の状態密度を操作することができる。構造によっては特定の色（エネルギー）の光を完全に禁止し、DOS がゼロとなるエネルギー領域、フォトニックバンドギャップを作り出すことができる。また、別種の構造ではある方向にのみ光の伝搬を抑制し、鏡、導波管、発振器を構成することもできる。このような繰り返し構造を[[フォトニック結晶]]と呼ぶ。ナノ構造を施した媒質では状態密度よりも場所ごとに異なる[[局所状態密度]] (LDOS) の考え方のほうがより適している。

== 状態密度の計算 ==
[[化合物]]や[[生体物質|生体分子]]、[[高分子]]など、興味の対象となる系は一般的に複雑である。これらの系は解析的に状態密度を計算するには複雑すぎ、ほとんどの場合それは不可能である。高精度の状態密度をコンピュータシミュレーションにより計算するアルゴリズムがいくつか知られている。その一つとして[[ワン・ランダウ法|ワン・ランダウのアルゴリズム]]が挙げられる<ref>{{cite journal|date=Mar 2001|title=Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of States|journal=Phys. Rev. Lett.|volume=86|issue=10|pages=2050–2053|publisher=American Physical Society|arxiv=cond-mat/0011174|bibcode=2001PhRvL..86.2050W|author1=Wang, Fugao|author2=Landau, D. P.|pmid=11289852|doi=10.1103/PhysRevLett.86.2050}}</ref>。

ワン・ランダウ法の枠組みの中では、状態密度に関する事前知識は一切必要がない。まず系のコスト関数（たとえばエネルギー）を離散化し、階級 {{Mvar|i}} に到達するごとに状態密度の[[ヒストグラム]] {{Math|''g''(''i'')}} を次のように更新する。
: <math>  g(i)  \rightarrow   g(i)  +   f  </math>
ここで {{Mvar|f}} は修正因子である。この階級に特定の回数 ({{Val|10|-|15}}) だけ到達するごとに修正因子は何らかの基準により減少させる。例えば以下のように処理する。
: <math> f_{n+1} \rightarrow (1/2)f_{n}</math>
ここで {{Mvar|n}} は更新が {{Mvar|n}} 回目であることを示す。特定の閾値に修正因子が到達したとき、たとえば {{Math|''f''<sub>n</sub> < {{1e-|8}}}} となったときにシミュレーションを終了する。

ワン・ランダウ法は[[マルチカノニカル法]]や[[レプリカ交換法]]などに比べていくつかの利点を持っている。たとえば、状態密度がシミュレーションの主目的として算出される。また、ワン・ランダウ法は完全に温度非依存である。この性質により、[[タンパク質]]のような非常にでこぼこした[[エネルギー地形]]を持つ系に対しても状態密度を計算することができる<ref name="ojeda">{{cite journal|year=2010|title=Electric Field-Driven Disruption of a Native beta-Sheet Protein Conformation and Generation of a Helix-Structure|journal=Biophysical Journal|volume=99|issue=2|pages=595–599|bibcode=2010BpJ....99..595O|author1=Ojeda, P.|author2=Garcia, M.|doi=10.1016/j.bpj.2010.04.040|pmid=20643079|pmc=2905109}}</ref>。

数学的には、状態密度は[[被覆写像]]を用いて形式化することができる<ref>{{cite journal|author=Adachi T. and [[Toshikazu Sunada|Sunada. T]]|date=1993|title=Density of states in spectral geometry of states in spectral geometry|journal=Comment. Math. Helvetici|volume=68|pages=480–493}}</ref>。

== 関連項目 ==
{{Columns-list|colwidth=20em|
* [[有効質量]]
* [[バンド構造]]
* [[:en:k·p perturbation theory|k·p perturbation theory]]
* [[半導体]]
* [[電気伝導体]]
* [[価電子帯]]
* [[クローニッヒ・ペニーのモデル]]
* [[強結合近似]]
* [[マフィンティンポテンシャル]]
}}

== 出典 ==
{{Reflist|2}}

== 関連文献 ==
* Chen, Gang. Nanoscale Energy Transport and Conversion. New York: Oxford, 2005
* Streetman, Ben G. and Sanjay Banerjee. Solid State Electronic Devices. Upper Saddle River, NJ: Prentice Hall, 2000.
* Muller, Richard S. and Theodore I. Kamins. Device Electronics for Integrated Circuits. New York: John Wiley and Sons, 2003.
* Kittel, Charles and Herbert Kroemer. Thermal Physics. New York: W.H. Freeman and Company, 1980
* Sze, Simon M. Physics of Semiconductor Devices. New York: John Wiley and Sons, 1981

== 外部リンク ==
* [http://nanohub.org/resources/5784 Online lecture:ECE 606 Lecture 8: Density of States] by M. Alam
* [https://web.archive.org/web/20121208070536/http://sapienzalab.org:80/scientists-shed-light-on-glowing-materials/ Scientists shed light on glowing materials] How to measure the Photonic LDOS

{{Normdaten}}
{{DEFAULTSORT:しようたいみつと}}
[[Category:密度]]
[[Category:固体物理学]]
[[Category:量子力学]]
[[Category:物理学の概念]]
[[Category:物性物理学]]
[[Category:バンド計算]]
[[Category:統計力学]]