본문으로 건너뛰기

11. Level alignment 보정 — DFT gap 오차와 DFT+Σ

지금까지의 transport 계산은 전부 semi-local functional(PBE) 위에서 이루어졌다. 그런데 PBE의 Kohn–Sham 준위를 그대로 믿으면 접합의 공명 위치가 체계적으로 틀리고, EFE_F가 gap 꼬리에 놓인 시스템에서는 T(EF)T(E_F)가 자릿수 단위로 과대평가될 수 있다. 이 장에서는 왜 틀리는지, 문헌 표준인 DFT+Σ(scissors) 보정이 무엇을 고치는지, 그리고 TBtrans에서 재컴파일 없이 보정을 주입하는 방법을 개념 수준에서 다룬다. VASP hybrid 계산(3장)이 이 워크플로우에서 맡는 역할이 여기서 드러난다.

학습 목표​

  • Kohn–Sham eigenvalue와 실제 첨가/제거 에너지(IP/EA)의 차이를 설명한다.
  • DFT+Σ 보정의 두 성분 — gas-phase 보정과 image-charge 스크리닝 — 을 구분한다.
  • 분자 부분공간 projector로 보정 Hamiltonian ΔH\Delta H를 만들어 TBT.dH로 주입하는 구현 전략을 이해한다.
  • 이 보정의 적용 전제와 한계를 알고, 결과를 정량 상한/정성 지표로 해석한다.

문제: PBE 공명 위치는 두 가지 이유로 틀린다​

첫째 — gas-phase gap 과소평가​

분자(또는 불순물)의 진짜 준위는 전자 제거/첨가 에너지다:

IP=E(N−1)−E(N),EA=E(N)−E(N+1).\mathrm{IP} = E(N-1) - E(N), \qquad \mathrm{EA} = E(N) - E(N+1).

정확한 이론에서 HOMO 준위는 −IP-\mathrm{IP}, LUMO 준위는 −EA-\mathrm{EA}에 대응해야 한다. 그러나 semi-local functional의 KS eigenvalue는 self-interaction 오차와 derivative discontinuity 누락 때문에 HOMO가 너무 높고 LUMO가 너무 낮다 — gap이 실제(gas-phase)보다 크게 작다. 작은 유기 분자에서 PBE gap이 실험 IP−EA gap의 절반 이하인 경우가 흔하다.

둘째 — 전극 스크리닝 누락​

분자를 금속 전극 사이에 놓으면 전자를 넣고 뺄 때 전극에 유도되는 image charge가 하전 상태를 안정화한다. 그 결과 접합 안의 gap은 gas-phase gap보다 다시 좁아진다. 대표 사례로 graphite 위 benzene의 gap이 gas-phase 10.5 eV에서 흡착 후 7.2 eV로 줄어드는 GW 계산이 있다 (Neaton, Hybertsen, Louie 2006). 이 비국소 분극 효과 역시 semi-local DFT에는 없다.

두 오차는 부호가 반대다(첫째는 gap을 너무 작게, 둘째를 무시하면 너무 크게). 그래서 PBE 결과가 우연히 그럴듯해 보이는 경우도 있지만, 이는 오차 상쇄일 뿐 통제된 예측이 아니다. 공명 위치가 수백 meV만 틀려도 gap 안의 T(EF)T(E_F)는 자릿수 단위로 민감하다 — 11장의 핵심 동기가 이것이다.

DFT+Σ: scissors형 준위 보정​

image charge 준위 재규격화

그림 1. 전극 스크리닝(image charge)에 의한 준위 재규격화 — gas-phase에서의 HOMO/LUMO gap이 접합 안에서는 2Σimage2\Sigma_{image}만큼 줄어든다. 인셋: 금속 표면 앞 전하와 거울 전하. (개념도)

Quek–Neaton 계열의 DFT+Σ는 GW 같은 고비용 방법 대신, 위 두 물리를 준위 shift 두 개로 근사한다. 점유 준위는 내리고 비점유 준위는 올리는 scissors 보정이다:

Δocc=(−IP−εHOMODFT)⏟Σ1: gas-phase+Wocc⏟Σ2: image charge,Δvirt=(−EA−εLUMODFT)+Wvirt.\Delta_{\mathrm{occ}} = \underbrace{\left(-\mathrm{IP} - \varepsilon_{\mathrm{HOMO}}^{\mathrm{DFT}}\right)}_{\Sigma^{1}:\ \text{gas-phase}} + \underbrace{W_{\mathrm{occ}}}_{\Sigma^{2}:\ \text{image charge}}, \qquad \Delta_{\mathrm{virt}} = \left(-\mathrm{EA} - \varepsilon_{\mathrm{LUMO}}^{\mathrm{DFT}}\right) + W_{\mathrm{virt}} .

성분 (i): gas-phase 보정 — VASP hybrid의 역할​

고립 분자의 IP/EA 기준값을 얻는 두 경로:

  • ΔSCF — 중성/양이온/음이온 총에너지 3점 계산. 문헌 표준이지만 주기 경계에서 하전 셀의 monopole 보정, 하전 상태의 스핀 배치 선택 같은 실무 문제가 따라온다.
  • hybrid functional eigenvalue — 같은 구조의 고립 분자를 B3LYP/HSE 등으로 single-point 계산하고, hybrid 준위를 기준으로 semi-local 준위와의 차이를 shift로 쓴다. ΔSCF의 대용 근사이며, 3장에서 준비한 VASP hybrid 계산이 바로 이 기준점을 제공한다.
코드 간 기준 정렬

transport는 SIESTA(LCAO), 기준 준위는 VASP(plane-wave)에서 나온다면 두 코드의 eigenvalue를 직접 빼면 안 된다 — functional 차이에 코드·기저 offset이 섞여 들어온다. 공통 functional(예: PBE)을 양쪽 코드에서 한 번씩 더 계산해 다리로 삼고, 뺄셈이 항상 같은 코드 안에서만 일어나도록 항을 구성하면 코드 offset이 상쇄된다:

Δ=[εhybrid−εPBE]VASP+[εPBE−εtransport functional]SIESTA.\Delta = \left[\varepsilon^{\mathrm{hybrid}} - \varepsilon^{\mathrm{PBE}}\right]_{\mathrm{VASP}} + \left[\varepsilon^{\mathrm{PBE}} - \varepsilon^{\mathrm{transport\ functional}}\right]_{\mathrm{SIESTA}} .

성분 (ii): image-charge 보정​

평면 금속 전극이라면 고전 image-charge 모델로 추정한다. 전극 표면에서 거리 zz에 있는 점전하의 안정화 에너지가 e2/4(z−z0)e^2/4(z-z_0) 꼴이고(z0z_0는 image plane 위치), 전극이 둘이면 무한 급수 합으로 확장된다. 점유 준위는 올라가고 비점유 준위는 내려가 gap을 다시 좁힌다 — 성분 (i)과 반대 방향이다.

구현 개요: projector ΔH를 TBT.dH로 주입한다​

보정이 에너지에 무관한 정적(static) shift라면, self-energy 보정은 Hamiltonian 보정과 수학적으로 같다. 즉 SIESTA를 다시 돌릴 필요 없이, TBtrans가 Green 함수를 만들 때 쓰는 HH에 ΔH\Delta H를 얹으면 된다:

Gr(E)=[(E+iη)S−(H+ΔH)−ΣL−ΣR]−1.G^r(E) = \left[(E + i\eta)S - (H + \Delta H) - \Sigma_L - \Sigma_R\right]^{-1}.

ΔH\Delta H는 고립 분자 MO를 기저로 한 projector로 구성한다. CC를 고립 분자 MO의 PAO 계수 행렬, Δi\Delta_i를 궤도별 shift라 하면

ΔH=S C diag(Δi) C† S\Delta H = S\, C\, \mathrm{diag}(\Delta_i)\, C^{\dagger}\, S

이고, 이렇게 만들면 일반화 고유값 문제에서 각 MO가 정확히 Δi\Delta_i만큼 이동한다:

(H+ΔH) cj=(εj+Δj) S cj.(H + \Delta H)\, c_j = (\varepsilon_j + \Delta_j)\, S\, c_j .

구현에서 주의할 점:

  • 비직교 기저에서는 SS가 양쪽에 반드시 끼어야 한다. 대각 onsite 원소만 shift하는 손쉬운 방법은 준위를 Δ\Delta만큼 옮기지 못하고 궤도 혼성까지 왜곡한다.
  • ΔH\Delta H는 분자(불순물) 부분공간에만 작용한다. 전극 Hamiltonian과 lead self-energy는 건드리지 않는다.
  • 완성된 ΔH\Delta H는 sisl로 netCDF 형식(*.dH.nc)으로 저장하고, TBtrans 입력에 TBT.dH 옵션으로 경로를 지정해 주입한다 (파일 형식과 정확한 옵션 표기는 TBtrans 매뉴얼과 sisl 문서의 delta 파일 항목을 확인).
  • 보정이 후처리 단계에만 들어가므로 zero-bias TSHS를 그대로 재사용한다 — TranSIESTA 재수렴이 필요 없다.
  • 자가 검증 두 가지: ① ΔH\Delta H를 MO 기저로 되돌려 diag(Δi)\mathrm{diag}(\Delta_i)가 회복되는지, ② 1차 섭동 ⟨ψj∣ΔH∣ψj⟩\langle\psi_j|\Delta H|\psi_j\rangle가 재대각화 결과와 일치하는지.

한 가지 개념적 어려움도 짚어 둔다. "분자 부분공간"은 분자가 전극과 약하게 물리흡착된 접합에서는 명확하지만, 우리 예제(C19N)처럼 불순물이 chain에 공유결합으로 융합된 경우에는 정의 자체가 모호하다. device Hamiltonian에서 부분 블록을 잘라 그대로 대각화하는 방식(MPSH류)은 융합된 계에서 비물리적 준위를 줄 수 있으므로, 고립 기준 계산에서 얻은 CC를 쓰는 projector 방식이 더 안전하다 — 단, 어떤 원자를 "분자"로 볼지는 여전히 사용자의 물리적 판단이다.

한계 — 정직하게 쓰기​

DFT+Σ는 근사의 사다리에서 명확한 위치가 있는 방법이고, 그 전제를 벗어나면 정량성이 무너진다.

  1. image-charge 모델은 전극 기하에 민감하다. 표준 평판 공식은 평면 금속 표면과 잘 정의된 image plane을 전제한다. 우리 예제 같은 1D chain 전극에는 image plane이 정의되지 않고, 1D 기하로 확장된 published 절차도 확립되어 있지 않다. 이 경우 Σ2\Sigma^2를 생략한 gas-phase 보정만 적용하게 되는데, 그러면 gap opening은 상한, T(EF)T(E_F) 억제는 하한 추정이 된다.
  2. ΔSCF 기준의 선택 문제. 하전 상태의 스핀 배치, 진공 준위 정렬, 주기 셀 보정에 따라 IP/EA가 달라진다. 어떤 기준을 썼는지 항상 명시해야 한다.
  3. 적용 전제: 약결합 + 공명이 EFE_F에서 멀 것. 보정의 정확도는 분자-전극 결합이 세지거나 공명이 EFE_F에 접근할수록 나빠진다 (Darancet 2012). 접합에서 전하 이동이 커서 frontier 준위가 EFE_F에 pinning된 계는 전제 자체가 흔들린다.
  4. non-self-consistent. 준위만 옮기고 전하 밀도는 되먹임하지 않는다.

따라서 보정된 T(E)T(E)의 수치는 그대로 인용하기보다 "보정 방향과 자릿수" 수준의 정성 지표, 또는 정량 상한/하한으로 해석하는 것이 안전하다.

예제 흐름: C19N에서의 개념 연습​

실제 파이프라인은 ① 고립 기준 계산(VASP hybrid 또는 ΔSCF) → ② Δi\Delta_i 산출 → ③ CC 추출과 ΔH\Delta H 조립 → ④ TBT.dH로 TBtrans 재실행 → ⑤ 보정 전후 T(E)T(E) 비교의 5단계다. 여기서는 보정이 T(E)T(E)를 어떻게 바꾸는지 감을 잡는 toy 모델만 돌려 본다. gap 안 tunneling 영역에서 단일 공명의 transmission은 Lorentzian으로 근사된다:

T(E)≈Γ2(E−ε0)2+Γ2.T(E) \approx \frac{\Gamma^2}{(E - \varepsilon_0)^2 + \Gamma^2}.
toy_scissors.py
import numpy as np
import matplotlib.pyplot as plt

E = np.linspace(-3, 3, 6001)          # eV, E - E_F
Gamma = 0.05                          # eV, 전극 결합 폭

def lorentz(E, eps0):
    return Gamma**2 / ((E - eps0)**2 + Gamma**2)

# 보정 전: 점유 공명이 E_F 바로 아래 / 보정 후: scissors로 0.5 eV 하강
for eps0, label in ((-0.10, "before (PBE)"), (-0.60, "after (DFT+Sigma)")):
    plt.semilogy(E, lorentz(E, eps0), label=label)

plt.xlabel("E - E_F (eV)")
plt.ylabel("T(E)")
plt.legend()
plt.savefig("scissors_sketch.png", dpi=200, bbox_inches="tight")

공명이 EFE_F에서 0.5 eV 멀어지는 것만으로 T(EF)T(E_F)가 두 자릿수 가까이 떨어지는 것을 볼 수 있다. 이것이 "PBE의 준위 오차가 전도도를 자릿수로 왜곡한다"의 정체이고, 반대로 보정값의 불확실성이 수백 meV면 T(EF)T(E_F) 예측도 그만큼 불확실하다는 뜻이기도 하다.

연습문제​

  1. toy 모델에서 Γ\Gamma를 0.01, 0.05, 0.2 eV로 바꿔가며 같은 0.5 eV shift가 T(EF)T(E_F)를 몇 배 바꾸는지 표로 정리하라. 어떤 영역에서 보정이 가장 극적인가.
  2. 2준위(HOMO/LUMO) 모형에서 Δocc<0\Delta_{\mathrm{occ}} < 0, Δvirt>0\Delta_{\mathrm{virt}} > 0 보정을 가한 뒤 gap 중앙의 TT가 어떻게 변하는지 두 Lorentzian의 합으로 스케치하라.
  3. 비직교 2-orbital 모형(S12≠0S_{12} \neq 0)을 하나 만들어, onsite 대각만 Δ\Delta만큼 올린 경우와 SC diag(Δ) C†SS C\,\mathrm{diag}(\Delta)\,C^{\dagger} S projector를 쓴 경우의 고유값 이동을 비교하라. 왜 SS가 양쪽에 필요한지 수치로 확인된다.
  4. 우리 C19N 예제에서 "분자 부분공간"을 어떻게 정의하겠는가. N 원자 하나만? N과 이웃 C 몇 개? 각 선택이 projector 보정에 어떤 영향을 줄지 논하라.

참고문헌​

  • J. B. Neaton, M. S. Hybertsen, and S. G. Louie, "Renormalization of Molecular Electronic Levels at Metal-Molecule Interfaces", Phys. Rev. Lett. 97, 216405 (2006). DOI 10.1103/PhysRevLett.97.216405
  • S. Y. Quek et al., "Amine-Gold Linked Single-Molecule Circuits: Experiment and Theory", Nano Lett. 7, 3477 (2007).
  • P. Darancet et al., "Quantitative Current-Voltage Characteristics in Molecular Junctions from First Principles", Nano Lett. 12, 6250 (2012). DOI 10.1021/nl3033137