slab 모델과 진공 기준으로 facet 안정성·결합·전자 방출을 계산하기
벌크 내부 원자는 사방에서 결합하지만 표면 원자는 배위가 끊기고 전자 밀도가 진공으로 새어 나옵니다. 같은 물질도 facet과 coverage에 따라 안정성과 반응성이 달라집니다.
주기 경계조건 안에서 열린 표면과 진공을 어떻게 만들고 에너지 기준을 맞출까요?
이번 장의 답충분히 두꺼운 slab와 진공을 반복하고 두 표면·면적을 고려해 표면에너지를 구합니다. 흡착에너지는 깨끗한 slab와 gas reference를 같은 방법으로 비교하며, 일함수는 진공준위와 Fermi 준위 차이입니다.
01slab 모델의 네 가지 길이
표면 법선 방향으로 유한한 원자층과 충분한 진공층을 두고 이를 주기적으로 반복합니다. 수렴시킬 길이는 slab 두께, 진공 두께, 평면 내 셀 크기와 k-point 밀도입니다. 아래층을 고정하면 bulk constraint를 흉내 내지만 층 수와 고정 범위를 검사해야 합니다.
대칭 slab는 두 표면이 같아 순 dipole이 작지만 흡착종이 양쪽에 필요할 수 있습니다. 비대칭 slab는 한쪽 표면만 연구하기 쉽지만 인공 전기장과 dipole correction을 점검해야 합니다.
02표면에너지와 흡착에너지의 기준
분자 기준은 충분한 진공·적절한 spin·동일한 pseudopotential와 범함수로 계산합니다. coverage가 달라지면 흡착종 사이 측면 상호작용과 표면 재구성이 달라지므로 단일 분자 에너지를 고유한 상수처럼 사용하면 안 됩니다.
03진공준위와 일함수
slab의 평면평균 electrostatic potential이 진공 중앙에서 평탄해지는 값을 $V_{vac}$라 두면 일함수는 전자를 Fermi 준위에서 진공으로 꺼내는 에너지 $\Phi=V_{vac}-E_F$입니다.
흡착으로 표면 dipole이 바뀌면 work function도 변합니다. 비대칭 slab의 양쪽 진공준위가 다를 수 있으므로 관심 면의 plateau를 사용하고 dipole correction·진공 두께를 수렴시킵니다.
- slab 층 수와 중앙층의 bulk 회복
- 진공 두께와 potential plateau
- 평면 k-grid와 lateral cell/coverage
- 고정층·완화층·대칭 및 dipole correction
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- slab model
- 유한 원자층과 진공을 주기적으로 반복해 표면을 표현하는 계산 모델입니다.
- 표면에너지surface energy
- 벌크를 절단해 단위면적의 표면을 만드는 초과 에너지입니다.
- 흡착에너지adsorption energy
- 흡착종과 깨끗한 표면이 결합할 때의 기준 상태 대비 에너지 변화입니다.
- 일함수work function
- 고체 Fermi 준위의 전자를 진공준위로 꺼내는 최소 에너지입니다.
- dipole correction
- 비대칭 slab의 주기 반복에서 생기는 인공 전기장을 제거하는 보정입니다.
- 표면 재구성surface reconstruction
- 끊어진 결합과 전하를 낮추기 위해 표면 원자 배열이 bulk 절단면과 다르게 바뀌는 현상입니다.
05참고자료와 더 읽을거리
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.
- Grimme, Semiempirical GGA-type density functional with long-range dispersion correction (2006) — 장거리 분산력을 원자쌍 보정으로 추가하는 DFT-D 접근의 대표 논문입니다.