PAW 방법

계산물리 양자화학 마지막 수정: 2026-07-20 04:19:22

1. 개요[편집]

PAW 방법 Projector Augmented-Wave method
제안	P. E. Blöchl (1994)
재정식화	Kresse & Joubert (1999)
정체	전전자 ↔ 유사 파동함수 선형변환
정확도 / 비용	전전자급 정확도, 유사퍼텐셜급 비용
대표 코드	VASP, Quantum ESPRESSO, ABINIT

PAW 방법(Projector Augmented-Wave method, 투영자 증강파동 방법)은 계산하기 쉬운 매끄러운 유사 파동함수와, 물리적으로 진짜인 전전자(all-electron) 파동함수 사이를 정확한 선형변환으로 이어 주는 전자구조 계산 기법이다. 1994년 페터 블뢰흘이 제안했고, 이름 그대로 “투영자(projector)“로 코어 영역의 정보를 되살려(augment) 붙인다. 핵심 셀링 포인트는 딱 하나 — 유사퍼텐셜만큼 싸게 계산하면서 전전자만큼 정확하다. 그래서 오늘날 VASP를 비롯한 평면파 밀도범함수이론 코드의 사실상 기본값이다.¹

PAW가 왜 필요한지 이해하려면, 이웃한 두 문서와의 역할 분담을 먼저 정리하는 게 빠르다. 이 셋은 자주 한 서랍에 뭉뚱그려지지만 층위가 다르다.

2. 유사퍼텐셜·평면파 기저와 무엇이 다른가[편집]

평면파 기저 — 파동함수를 무엇으로 전개하느냐의 문제다. 역격자 벡터로 지수화한 평면파의 합. 즉 *기저(basis)*의 선택이다.
유사퍼텐셜 — 코어 근처를 어떻게 처리하느냐의 문제다. 골치 아픈 코어 전자와 핵의 강한 퍼텐셜을 부드러운 유효 퍼텐셜로 **대체(replace)**해 버린다. 이 과정에서 코어 정보는 원리적으로 버려진다.
PAW 방법 — 역시 코어 처리 기법이지만, 유사퍼텐셜과 결정적으로 다르다. 코어 근처의 격렬한 진동을 매끄럽게 다루되, 그 정보를 버리지 않고 변환으로 되살릴 수 있게 보존한다.

한 줄 요약: 평면파는 기저, 유사퍼텐셜과 PAW는 코어 처리 방식이고, 그중 유사퍼텐셜은 코어를 봉인해 버리는 반면 PAW는 봉인하되 열쇠를 남겨둔다. 그래서 PAW는 종종 “유사퍼텐셜과 전전자 방법의 중간”이라고 불린다.² 셋은 경쟁 관계가 아니라 평면파 기저 위에서 PAW(또는 유사퍼텐셜)로 코어를 처리하는 식으로 함께 쓰인다.

3. 변환 연산자[편집]

PAW의 수학적 심장은 매끄러운 유사 파동함수 $|\tilde\psi\rangle$ 를 진짜 전전자 파동함수 $|\psi\rangle$ 로 보내는 선형 변환 $\mathcal{T}$ 이다.

|\psi\rangle = \mathcal{T}\,|\tilde\psi\rangle = |\tilde\psi\rangle + \sum_i \left( |\phi_i\rangle - |\tilde\phi_i\rangle \right)\langle \tilde p_i | \tilde\psi \rangle

이 식이 PAW의 전부라 해도 과언이 아니다. 뜯어보면 이렇다.

$|\tilde\psi\rangle$ : 평면파로 값싸게 전개되는 매끄러운 유사 파동함수. 실제로 코드가 변수로 들고 다니는 것.
$|\phi_i\rangle$ : 각 원자 주위에서 정확한 슈뢰딩거 방정식을 푼 전전자 부분파(partial wave). 코어 근처의 진짜 마디까지 담고 있다.
$|\tilde\phi_i\rangle$ : 그 부분파를 코어 영역에서 매끄럽게 다듬은 유사 부분파.
$\langle\tilde p_i|$ : 투영자(projector). 유사 파동함수에서 각 부분파 성분을 뽑아내는 역할.

직관은 명쾌하다. 원자를 둘러싼 증강 영역(augmentation sphere) 안에서만 매끄러운 버전 $|\tilde\phi_i\rangle$ 을 빼고 진짜 버전 $|\phi_i\rangle$ 을 도로 끼워 넣는 것이다. 영역 밖에서는 두 부분파가 같아 보정항이 0이 되므로, 결과적으로 코어 근처의 급격한 진동과 마디가 정확히 복원된다. 유사퍼텐셜이 근사로 버린 코어 정보를, PAW는 이 변환식으로 원할 때마다 재구성한다. 전자밀도, 힘, 초미세 상호작용, NMR 화학 이동처럼 핵 근처 파동함수의 정확한 값이 필요한 물리량에서 PAW가 유사퍼텐셜을 압도하는 이유다.

4. 에너지와 세 조각 분해[편집]

PAW에서 총에너지를 비롯한 모든 물리량은 세 항으로 깔끔하게 쪼개진다. 매끄러운 부분은 평면파 격자에서, 증강 영역 안의 두 항(진짜 부분파 기여, 유사 부분파 기여)은 각 원자 중심의 방사형 격자(radial grid)에서 따로 계산한다.

E = \tilde{E} + \sum_a \left( E^1_a - \tilde{E}^1_a \right)

$\tilde E$ 는 평면파에서 값싸게 얻는 매끄러운 항, $E^1_a$ 는 원자 $a$ 증강구 안의 전전자 항, $\tilde E^1_a$ 는 같은 구 안의 유사 항이다. 뒤의 두 항은 원자별로 국소화되어 있어 한 번 만들어 두면 재활용되는 데이터셋으로 미리 계산해 배포한다. 이것이 원소마다 하나씩 있는 “PAW 데이터셋”의 정체이며, 사용자는 평면파 기저의 컷오프 에너지처럼 이 데이터셋을 골라 넣기만 하면 된다. 계산은 여전히 콘-샴 방정식을 자기무장(self-consistent)으로 반복해 푸는 표준 DFT 흐름 위에서 돌아간다.

5. 실무에서[편집]

정확도의 표준이 되다 — 여러 코드의 상태방정식을 대규모로 대조한 국제 벤치마크에서, 잘 만든 PAW 데이터셋은 순수 전전자 코드와 실용적으로 구분이 안 될 만큼 일치한다는 결론이 났다. PAW가 “값싼 근사”가 아니라 “정확한 방법”으로 대접받게 된 결정적 계기다.
데이터셋 품질이 전부 — 유사퍼텐셜과 마찬가지로 PAW도 데이터셋이 틀리면 계산은 멀쩡히 수렴하면서 조용히 틀린 숫자를 낸다. 어떤 원자가 상태를 코어로 얼렸는지(예: 전이금속의 semicore p 상태 포함 여부)가 밴드 구조 계산과 자기 모멘트를 좌우한다.
증강 전하와 이중 격자 — 유사 밀도만으로는 부족해 증강 전하(compensation charge)를 더해야 하므로, 밀도용 격자를 파동함수용보다 훨씬 촘촘히(VASP의 ENAUG) 잡아야 하는 경우가 있다. 이 점은 울트라소프트 유사퍼텐셜과 사정이 비슷하다.³
비교하자면 — 하트리-폭 방법이나 국소 기저를 쓰는 분자 양자화학 진영과 달리, PAW는 주기계·평면파 문명의 산물이다. 고체·표면·결함 계산이 주무대다.

6. 관련 문서[편집]

블뢰흘의 1994년 원논문은 처음엔 즉각적인 대세가 아니었다. 판을 뒤집은 건 1999년 크레세와 주베르가 PAW를 울트라소프트 유사퍼텐셜과 같은 틀로 재정식화해 VASP에 이식하면서다. 이론은 블뢰흘이, 대중화는 VASP가 했다는 게 이 바닥의 공통된 평가. ↩
그래서 “PAW는 유사퍼텐셜이냐 아니냐”로 티키타카가 자주 벌어진다. 형식상으로는 유사퍼텐셜 계열의 방법론을 포함하지만, 코어 정보를 복원 가능하게 남겨둔다는 점에서 전전자 방법에 더 가깝다. 실무에선 “PAW 데이터셋”이라 부르며 유사퍼텐셜 라이브러리와 같은 서랍에 넣어두지만, 엄밀한 사람 앞에서 “PAW = 유사퍼텐셜”이라고 하면 정정당한다. ↩
계산은 잘 되는데 힘이나 응력이 미묘하게 안 맞으면 ENAUG(증강 격자 컷오프)를 의심하라는 게 VASP 사용자들 사이의 정설이다. 파동함수 컷오프만 올리고 밀도 격자를 방치하면 생기는 전형적 증상. ↩