상태밀도

계산물리 양자역학 마지막 수정: 2026-07-13 04:30:11

1. 개요[편집]

상태밀도(Density of States, DOS)는 단위 에너지 구간당 전자가 차지할 수 있는 양자 상태가 몇 개나 있는지를 나타내는 함수 $g(E)$ 다. 쉽게 말해 “에너지 $E$ 근처에 자리(좌석)가 얼마나 촘촘히 깔려 있는가”를 세는 지표다. 고체 안 전자는 밴드 구조 계산이 보여주듯 연속적인 에너지 밴드를 이루는데, 밴드가 평평한 곳에는 상태가 빽빽하게 몰려 DOS가 크고, 밴드가 가파른 곳에는 성기게 퍼져 DOS가 작다. 물질이 금속이냐 절연체냐, 어떤 빛을 흡수하느냐, 열을 얼마나 저장하느냐가 전부 이 $g(E)$ 의 생김새에 새겨져 있다.¹

수식으로는 밴드 에너지 $\varepsilon_n(\mathbf{k})$ 를 브릴루앙 영역 전체에 걸쳐 에너지별로 합산한 것으로 정의한다.

g(E) = \sum_n \int_{\text{BZ}} \frac{d^3k}{(2\pi)^3}\, \delta\!\left(E - \varepsilon_n(\mathbf{k})\right)

여기서 $n$ 은 밴드 지수, 적분은 파수공간(k-space)의 브릴루앙 영역(BZ)에 대한 것이며, 디랙 델타 함수 $\delta$ 가 “에너지가 정확히 $E$ 인 상태만 골라 세는” 필터 역할을 한다. DFT 계산 결과물 중 밴드 구조와 함께 가장 먼저 뽑아 보는 그림이 바로 이 DOS다.

2. 자유전자 모형: √E 의 세계[편집]

가장 단순한 출발점은 상호작용 없는 자유전자 기체다. 이 경우 에너지가 운동량의 제곱에 비례( $\varepsilon = \hbar^2 k^2 / 2m$ )하므로, 3차원에서 상태밀도는 에너지의 제곱근에 비례한다.

g(E) \propto \sqrt{E}

이 $\sqrt{E}$ 꼴은 금속의 전도 전자를 이해하는 교과서적 기준선이다. 흥미롭게도 차원이 낮아지면 모양이 확 바뀐다. 2차원 전자계에서는 $g(E)$ 가 계단 함수(상수)가 되고, 1차원 나노선에서는 $1/\sqrt{E}$ 로 발산하는 뾰족한 봉우리(반 호브 특이점, van Hove singularity)가 나타난다. 저차원 소재의 독특한 광·전자 물성이 상당 부분 이 DOS 모양 차이에서 비롯된다.²

3. 페르미 준위와 금속/절연체 판별[편집]

전자는 페르미온이라 파울리 배타 원리를 따르므로, 절대영도에서 낮은 에너지 상태부터 차곡차곡 채워 올라간다. 이 채워진 맨 꼭대기 에너지가 페르미 준위(Fermi level) $E_F$ 다. DOS와 페르미 준위의 관계만 보면 물질의 전기적 성격이 대번에 갈린다.

금속: $E_F$ 가 밴드 한가운데를 지나 그 자리에 상태가 있다. $g(E_F) \neq 0$ 이므로 전자가 조금만 에너지를 얻어도 곧장 빈 상태로 이동할 수 있어 전류가 흐른다.
절연체·반도체: $E_F$ 가 상태가 텅 빈 구간, 즉 밴드갭(band gap) 안에 놓인다. $g(E)=0$ 인 이 금지 영역을 뛰어넘어야 전도가 가능하니, 갭의 크기가 곧 절연 성능이다.

그래서 계산된 DOS에서 페르미 준위 근처에 상태가 있는지 없는지를 확인하는 것이, 새 물질이 금속인지 반도체인지 판정하는 첫 관문이다. 다만 표준 밀도범함수이론은 밴드갭을 과소평가하는 고질병이 있어, DOS로 읽은 갭 값을 그대로 믿으면 안 된다는 점은 늘 따라붙는 경고다.³

4. 부분 상태밀도(PDOS)[편집]

전체 DOS는 “총 좌석 수”만 알려줄 뿐, 그 상태가 어느 원자의 어느 궤도(s, p, d, f)에서 왔는지는 말해주지 않는다. 이걸 분해한 것이 부분 상태밀도(Projected/Partial DOS, PDOS)다. 각 콘-샴 상태를 특정 원자의 특정 궤도에 사영(projection)해 기여도를 나눈다.

PDOS는 화학 결합과 전자 구조를 해부하는 강력한 현미경이다. 예컨대 전이금속 산화물에서 페르미 준위 근처 상태가 금속 $d$ 궤도에서 오는지 산소 $p$ 궤도에서 오는지를 보면 자성·촉매 활성의 기원을 짚을 수 있다. 두 원자의 PDOS 봉우리가 같은 에너지에서 겹치면 그 원자들끼리 공유 결합(혼성)을 이룬다는 신호이기도 하다. 고유값 문제를 풀어 얻은 궤도 계수가 있어야 계산되므로, PDOS는 밴드 계산의 부산물로 자연스럽게 딸려 나온다.

5. 어떻게 계산하나: 스미어링과 테트라헤드론[편집]

정의식의 델타 함수는 수학적으로는 우아하지만 컴퓨터가 직접 다루기엔 곤란하다. 무한히 뾰족한 봉우리를 유한 개의 k-점 샘플로 적분해야 하기 때문이다. 실무에서는 두 가지 처방을 쓴다.

가우시안 스미어링(Gaussian smearing): 델타 함수를 폭 $\sigma$ 짜리 가우스 봉우리로 뭉갠다. 구현이 간단하지만 $\sigma$ 를 크게 잡으면 날카로운 특징이 뭉개지고, 작게 잡으면 k-점이 부족해 DOS가 톱니처럼 우글거린다. 이 폭 조절이 은근히 손이 많이 간다.
테트라헤드론 방법(tetrahedron method): 브릴루앙 영역을 사면체로 잘게 쪼갠 뒤 각 사면체 안에서 밴드를 선형 보간해 델타 적분을 해석적으로 수행한다. 스미어링 폭이라는 인위적 변수가 없어 정확한 DOS를 주므로, 정밀한 계산에서는 이쪽이 표준이다.⁴

DOS는 여기서 그치지 않고 비열, 전기 전도도, 광흡수 스펙트럼, 초전도 전이온도 추정 등 온갖 물성 계산의 출발 재료로 재활용된다. 전자 구조 계산에서 밴드 구조가 “지도”라면, 상태밀도는 그 지도를 에너지 축 하나로 압축한 “인구 통계”인 셈이다.⁵

6. 관련 문서[편집]

반 호브 특이점(van Hove singularity)이라 불리는 DOS의 뾰족한 봉우리들은 밴드가 극값(평평)을 이루는 지점에서 생긴다. 이 특이점이 페르미 준위 근처에 오면 초전도나 자성 같은 극적인 현상이 잘 일어나서, 물성 물리학자들이 눈에 불을 켜고 찾는 지형지물이다. ↩
그래서 나노선·양자우물·2차원 소재가 벌크와 전혀 다른 광학·수송 특성을 보인다. “차원을 낮췄더니 물성이 바뀌더라”의 배후에는 거의 항상 DOS 모양의 변화가 있다. ↩
이른바 “밴드갭 문제”. LDA/GGA 범함수는 실제 갭을 30~100% 과소평가하기도 한다. 진지하게 갭을 논하려면 하이브리드 범함수나 GW 근사로 보정하는 것이 정석이다. DOS 그림에서 갭이 좁아 보인다고 “이거 금속이네” 하고 성급하게 논문 쓰면 곤란하다. ↩
특히 금속의 페르미 면 적분처럼 페르미 준위 근처가 예민한 계산에서는 블뢰흘(Blöchl)이 보정한 테트라헤드론 방법이 사실상 국룰이다. 스미어링은 빠른 대신 정밀도에서 밀린다. ↩
밴드 구조 그림은 방향( $\Gamma$ -X-M 같은 k-경로)별로 에너지를 보여주는 반면, DOS는 방향 정보를 다 적분해 버려 에너지 축 하나로 뭉갠다. 둘을 나란히 붙여 그리는 이유는 서로가 서로의 부족한 정보를 메워주기 때문이다. ↩