計算材料科学

LDAやGGAの交換相関汎関数を使った標準的なDFT計算の結果の傾向として正しいものはどれか。

  1. A半導体や絶縁体のバンドギャップを実験値よりおおむね大きく見積もる傾向があり、シリコンではGGAで約2 eVになる
  2. B半導体や絶縁体のバンドギャップを実験値より小さく見積もる傾向があり、シリコンではGGAで約0.6 eVになる
  3. C0 Kの基底状態だけでなく、有限温度の自由エネルギーや融点も、構造を1回最適化するだけで直接得られる
答えと解説を見る

正解:B. 半導体や絶縁体のバンドギャップを実験値より小さく見積もる傾向があり、シリコンではGGAで約0.6 eVになる

LDAやGGAでは、コーン・シャム方程式の固有値の差からバンドギャップを読むと、実験値より小さくなることがよく知られています。シリコンではGGAで約0.6 eVと、実験値1.12 eVのおよそ半分で、大きく見積もるとする説明は傾向が逆です。正確な値が必要なときは、ハイブリッド汎関数やGW近似を使います。一方、格子定数、弾性率、形成エネルギーのような基底状態の全エネルギーに関わる量は、多くの物質でよい精度で求まります。通常のDFTは0 Kの基底状態の計算なので、有限温度の性質を得るには、格子振動(フォノン)の計算や分子動力学を組み合わせる必要があります。また、標準的な汎関数はファンデルワールス力をうまく表せないので、層状物質や吸着では補正を加えます。

誤りを報告

この分野を解く 分野ページへ

計算材料科学のほかの問題

  1. ある金属の108原子のスーパーセル(周期的な境界条件)をDFTで計算すると、全エネルギーは −405.00 eV だった。原子を1個取り除いて空孔を1つつくり、構造を緩和すると −400.10 eV になった。空孔形成エネルギーを E_v = E(N−1) − [(N−1)/N]E(N) で求めるとどれに最も近いか。
  2. 分子動力学法で、時間刻み2 fsで運動方程式を積分して、1 nsの間の原子の動きを計算したい。必要な時間ステップの数はどれか。
  3. 分子動力学法で使う原子間ポテンシャルについての説明として正しいものはどれか。
  4. レナード・ジョーンズ・ポテンシャル V(r) = 4ε[(σ/r)¹² − (σ/r)⁶] で、アルゴンを σ = 0.34 nm として表す。エネルギーが最小になる原子間距離 r₀ はどれに最も近いか。
  5. 分子動力学法の計算で、単原子の系の原子1個あたりの運動エネルギーの平均が0.0517 eVだった。エネルギー等分配則から求めた系の温度はどれに最も近いか。ボルツマン定数を k = 8.617×10⁻⁵ eV/K とする。