유사퍼텐셜

양자화학 계산물리 계산화학 마지막 수정: 2026-07-07 06:27:51

유사퍼텐셜 Pseudopotential
분야	계산물리 × 양자화학
목적	코어 전자를 유효 퍼텐셜로 대체
쓰이는 곳	밀도범함수이론 평면파 계산
대표 종류	노름보존, 울트라소프트, PAW
대표 소프트웨어	VASP, Quantum ESPRESSO

1. 개요[편집]

원자핵 근처의 전자는 미친 듯이 진동한다. 그걸 다 그리느니, 안 보이게 덮어버리자.

유사퍼텐셜(pseudopotential, 擬似-)은 원자핵에 강하게 묶인 코어(core) 전자와 핵의 강한 쿨롱 퍼텐셜을, 화학적으로 중요한 원자가(valence) 전자 입장에서 본 부드러운 유효 퍼텐셜로 대체하는 근사 기법이다. 첫 원리(ab initio) 밀도범함수이론 계산, 특히 평면파 기저를 쓰는 계산에서 계산량을 극적으로 줄이기 위해 거의 필수적으로 동원된다.

동기는 단순하다. 결합·반응·전도 같은 화학과 물성은 거의 전부 바깥쪽 원자가 전자가 담당한다. 반면 안쪽 코어 전자는 원자가 어떤 환경에 놓이든 거의 변하지 않는 “얼어붙은” 상태다. 그런데 이 코어 전자와 핵 근처의 원자가 파동함수는 마디(node)를 만들며 격렬하게 진동해서, 이를 정직하게 표현하려면 어마어마하게 많은 평면파가 필요하다.¹ 안 변하는 코어를 위해 계산 자원을 갈아넣는 건 낭비다. 그래서 코어를 통째로 유효 퍼텐셜로 봉인하고, 원자가 파동함수를 핵 근처에서 매끄럽게(마디 없이) 다시 그린다.

2. 아이디어[편집]

핵심은 실제(all-electron) 파동함수를 두 영역으로 나누는 것이다.

컷오프 반지름 $r_c$ 바깥: 유사파동함수(pseudo wavefunction)가 실제 파동함수와 정확히 일치해야 한다. 화학이 벌어지는 영역이므로 손대면 안 된다.
$r_c$ 안쪽: 마디를 없애고 부드럽게 만든다. 어차피 코어 영역이라 화학에 직접 관여하지 않으므로 형태를 바꿔도 된다.

이렇게 만든 유사파동함수는 핵 근처에서 진동하지 않으므로 적은 수의 기저로 표현된다. 계산은 실제 퍼텐셜 $V(r) \sim -Z/r$ 대신, 발산하지 않는 부드러운 유사퍼텐셜 $V_{\text{ps}}(r)$ 을 쓴다.

\hat{H}_{\text{ps}} = -\frac{\hbar^2}{2m}\nabla^2 + V_{\text{ps}}(\mathbf{r})

유사퍼텐셜은 보통 각운동량 $l$ 마다 다른 성분을 갖는 비국소(non-local) 형태다. s, p, d 전자가 핵을 각각 다르게 느끼기 때문이다.

3. 종류[편집]

시대순으로 세 세대가 있다.

종류	특징	장단점
노름보존 (norm-conserving)	$r_c$ 안쪽 전하량을 실제와 일치시킴	정확하고 이론이 깔끔, 대신 여전히 무거움
울트라소프트 (ultrasoft, USPP)	노름보존 조건을 완화해 더 부드럽게	평면파 수 급감, 대신 구현 복잡
PAW	코어 정보를 복원 가능하게 보존	정확도·효율 균형, 현대 표준

노름보존 유사퍼텐셜은 유사·실제 파동함수가 $r_c$ 안에서 같은 전하량(norm)을 갖도록 강제한다. 산란 성질을 잘 보존해 신뢰성이 높다. 울트라소프트(밴더빌트, Vanderbilt)는 이 노름보존 조건을 버려 파동함수를 훨씬 부드럽게 만들고, 대신 잃어버린 전하를 보정항으로 되살린다. 이후 블뢰흘(Blöchl)의 PAW(Projector Augmented-Wave) 방법이 코어 근처 정보를 투영자로 복원 가능하게 유지하면서 효율까지 챙겨, 오늘날 VASP를 비롯한 대부분의 코드에서 사실상 표준으로 자리 잡았다.²

4. 생성과 전이성[편집]

유사퍼텐셜은 하늘에서 떨어지지 않는다. 고립된 원자에 대해 하트리-폭 방법이나 DFT로 실제 전(全)전자 계산을 먼저 수행한 뒤, 그 결과를 재현하도록 유사퍼텐셜을 역설계한다. 이때 가장 중요한 덕목이 전이성(transferability)이다. 원자 하나로 만든 유사퍼텐셜이 결정·분자·표면 등 전혀 다른 환경에서도 실제 계산을 잘 흉내 내야 한다는 뜻이다.

전이성과 부드러움(적은 평면파)은 대체로 상충한다. 너무 부드럽게 만들면 컷오프 에너지는 낮아지지만 낯선 환경에서 답이 틀어진다. 그래서 유사퍼텐셜 제작은 이 둘 사이의 줄타기이며, 검증된 라이브러리(pseudopotential library)를 가져다 쓰는 게 국룰이다. 직접 만드는 건 전문가 영역이고, 잘못 만든 유사퍼텐셜은 밴드 구조 계산 결과를 조용히 오염시킨다.³

5. 계산에서의 위치[편집]

유사퍼텐셜은 평면파 DFT 파이프라인의 입력 데이터다. 전형적인 흐름은 다음과 같다.

각 원소마다 검증된 유사퍼텐셜(또는 PAW 데이터셋)을 준비한다.
평면파 기저의 컷오프 에너지 $E_{\text{cut}}$ 를 수렴할 때까지 올린다. 이 값은 쓰는 유사퍼텐셜이 얼마나 부드러운지에 직접 좌우된다.
슈뢰딩거 방정식에 해당하는 콘-샴(Kohn-Sham) 방정식을 자기무장(self-consistent)으로 반복해 풀어 전자밀도를 얻는다.
총에너지, 밴드 구조 계산, 상태밀도, 힘·응력 등을 뽑는다.

여기서 얻은 원자 간 힘은 분자동역학의 첫 원리 시뮬레이션으로 연결되기도 한다. 즉 유사퍼텐셜은 물질의 전자구조부터 동역학까지 이어지는 계산 사슬의 첫 단추다. 이 단추를 잘못 끼우면 아래 계산이 전부 어긋난다.

6. 관련 문서[편집]

핵 근처 원자가 파동함수는 코어 전자와 직교해야 해서 마디가 생기고, 그 마디를 표현하려면 짧은 파장 성분, 즉 고에너지 평면파가 잔뜩 필요하다. 컷오프 에너지가 하늘로 치솟는 주범이 바로 이 마디들이다. ↩
엄밀히 말하면 PAW는 유사퍼텐셜과 전전자 방법의 중간쯤에 있다. 코어 근처 정보를 버리지 않고 투영자로 되살릴 수 있게 남겨두기 때문이다. 그래서 “PAW는 유사퍼텐셜이 아니다”라고 깐깐하게 지적하는 사람도 있는데, 실무에선 대개 같은 서랍에 넣어둔다. ↩
유사퍼텐셜이 틀려도 계산은 멀쩡히 수렴하고 그럴듯한 숫자를 뱉는다는 게 함정이다. 그래서 논문에 쓴 유사퍼텐셜의 출처와 컷오프 수렴을 명시하는 게 계산 재현성의 기본 예의다. ↩