장거리 전자 요동의 비국소 상관을 pairwise·many-body·비국소 범함수로 보완하기

DFT · PART 5 · 숫자의 정확도와 이론의 적용 범위를 시험하는 동안02 / 08

층상재료, 분자결정, 유기분자 흡착은 전하가 크게 겹치지 않아도 서로 끌립니다. 밀도와 기울기만 보는 범함수는 멀리 떨어진 두 영역의 상관된 요동을 충분히 알기 어렵습니다.

장거리 분산력은 왜 준국소 범함수에서 사라지며 어떤 방식으로 되살릴까요?

이번 장의 답London dispersion은 순간 dipole 요동의 장거리 상관으로 $-C_6/R^6$ 꼬리를 갖습니다. pairwise DFT-D, density-dependent 비국소 범함수, many-body 모델 중 계에 맞는 방법을 선택하고 중거리 damping·screening·구조를 검증합니다.

01분산력의 양자적 기원

평균 dipole이 0인 두 중성 원자도 순간적인 전자 요동으로 dipole을 만들고 서로의 요동을 상관시켜 에너지를 낮춥니다. 멀리서는 대표적으로 $-C_6/R^6$으로 감소합니다.

LDA는 우연히 결합을 만들 수 있지만 올바른 장거리 꼬리를 갖지 않고, GGA도 두 밀도 영역이 겹치지 않으면 서로의 존재를 알기 어렵습니다. 따라서 물리적 분산 처리를 별도로 포함해야 합니다.

02pairwise에서 many-body까지

접근핵심주의
DFT-D원자쌍 $C_n/R^n$ 보정범함수별 damping·환경 의존
vdW-DF류밀도 이중적분 비국소 상관교환 선택·중거리 구조
TS/MBD류밀도 기반 polarizability와 many-body 응답금속 screening·매개변수 범위

작은 분자·유기결정에서는 pairwise가 효율적이지만 큰 공액계·층상물질·금속 표면에서는 many-body screening과 비가산성이 중요할 수 있습니다. 같은 ‘vdW 포함’ 결과도 방법의 물리가 다릅니다.

03결합에너지와 구조를 함께 검증하기

분산 보정은 결합거리·층간거리·흡착높이·에너지를 동시에 바꿉니다. 에너지 하나만 실험과 맞아도 구조가 틀리면 힘과 진동·반응경로가 잘못될 수 있습니다.

vdW 계산 점검
  • 보정과 기본 범함수의 권장 조합 사용
  • 진공·셀 크기와 장거리 cutoff 수렴
  • 흡착높이·층간거리·결합에너지 동시 benchmark
  • 금속·이온계에서 screening과 전하 상태 적합성 확인

04이 글의 용어 사전

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

London dispersion
순간 dipole 요동 사이의 상관에서 생기는 장거리 인력입니다.
polarizability
외부 전기장에 의해 전자분포와 dipole이 얼마나 쉽게 변하는지 나타내는 응답량입니다.
C6 계수
두 종 사이 장거리 $-C_6/R^6$ 분산 세기를 정하는 계수입니다.
damping function
짧은 거리에서 분산 보정과 기본 DFT 상관의 이중계수를 줄이는 함수입니다.
many-body dispersion
원자쌍 합을 넘어 여러 dipole 요동의 집단적 결합과 screening을 포함한 분산 처리입니다.
비가산성nonadditivity
전체 상호작용이 독립적인 원자쌍 에너지 합으로 표현되지 않는 성질입니다.

05참고자료와 더 읽을거리

  1. Grimme, Semiempirical GGA-type density functional with long-range dispersion correction (2006)장거리 분산력을 원자쌍 보정으로 추가하는 DFT-D 접근의 대표 논문입니다.
  2. Sholl & Steckel, Density Functional Theory: A Practical Introduction평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
  3. Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996)재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.
NEXT DFT NOTE스핀 분극과 자성 계산