Quantum ESPRESSO

소프트웨어 양자화학 계산물리 마지막 수정: 2026-07-16 04:31:57

1. 개요[편집]

Quantum ESPRESSO opEn-Source Package for Research in Electronic Structure, Simulation, and Optimization
약칭	QE
기반 이론	밀도범함수이론 (DFT)
기저 함수	평면파 (plane wave)
언어	Fortran 90 (MPI + OpenMP)
개발	DEMOCRITOS / SISSA (트리에스테) 및 국제 커뮤니티
라이선스	GPL (무료·오픈소스)

이름에 에스프레소가 들어가지만, 계산이 끝나기를 기다리는 동안 마시는 커피는 대개 에스프레소가 아니라 대용량 아메리카노다.

Quantum ESPRESSO(약칭 QE)는 평면파 기저와 유사퍼텐셜을 사용하는 오픈소스 밀도범함수이론 계산 패키지다. 이탈리아 트리에스테의 DEMOCRITOS/SISSA를 중심으로 파올로 잔노치(Paolo Giannozzi) 등이 여러 갈래의 기존 코드(PWscf, CP90, FPMD)를 하나의 배포판으로 통합해 2009년에 정식 발표했고, GPL 라이선스로 배포된다.¹ 이름은 억지스러운 두문자 조합이지만, 개발진이 이탈리아 사람이라 아무도 토를 달지 않았다.

VASP와 동일한 이론적 토대(평면파 + 유사퍼텐셜 + 주기경계조건)에 서 있으면서 값이 0원이라는 점, 그리고 소스가 열려 있어 새 기능을 얹기 좋다는 점 때문에 학계에서 특히 강한 존재감을 가진다. 응집물질물리·촉매·재료 계산에서 VASP의 가장 현실적인 대안으로 꼽힌다.

2. 구성 요소[편집]

QE는 하나의 실행 파일이 아니라 여러 실행 파일의 연합체다. 각 단계마다 다른 바이너리를 호출하고 그 결과 디렉토리를 다음 단계에 물려주는 구조라, 사용자는 자연스레 셸 스크립트 장인이 된다.

pw.x (PWscf) — 심장. 자기무장(SCF) 계산, 구조 최적화, 가변셀 최적화, 밴드 계산, 제일원리 분자동역학.
ph.x (PHonon) — 밀도범함수 섭동이론(DFPT) 기반 포논 분산, 유전 텐서, 보른 유효 전하, 전자-포논 결합.
cp.x (CP) — 카-파리넬로 분자동역학.
dos.x / projwfc.x / bands.x — 상태밀도와 투영 상태밀도(PDOS), 밴드 구조 계산 후처리.
pp.x — 전하밀도·퍼텐셜·STM 상 등 필드량 추출.
neb.x — 넛지 탄성 밴드법으로 반응 경로와 활성화 장벽.

여기에 와니에 함수 변환(Wannier90 연계), 전자-포논 계산(EPW), TDDFT 모듈 등이 위성처럼 붙는다.

3. 입력의 문법[편집]

QE의 입력은 포트란 네임리스트 + 카드 조합이다. VASP가 파일 네 개(INCAR/POSCAR/POTCAR/KPOINTS)로 역할을 쪼갠 것과 달리, QE는 텍스트 파일 하나에 전부 몰아넣는다.

&CONTROL — 계산 종류(scf, relax, vc-relax, nscf, bands, md), 출력 경로.
&SYSTEM — 셀과 원자 종류 수, 평면파 컷오프, 스미어링, 스핀.
&ELECTRONS — SCF 수렴 기준, 전하밀도 믹싱.
&IONS, &CELL — 이온·셀 이완 알고리즘.
카드: ATOMIC_SPECIES, ATOMIC_POSITIONS, K_POINTS, CELL_PARAMETERS.

주의할 함정이 두 개 있다. 첫째, **단위가 리드베리(Ry)**다. eV에 익숙한 채로 ecutwfc = 500을 적으면 500 Ry = 약 6800 eV라는 초현실적 컷오프로 계산이 영원히 끝나지 않는다.² 둘째, 파동함수 컷오프 ecutwfc와 전하밀도 컷오프 ecutrho가 따로 있다. 노름보존 유사퍼텐셜이면 기본값인 4배가 맞지만, 울트라소프트나 PAW를 쓰면 증강 전하가 날카로워서 8~12배를 줘야 한다. 이 값을 방치한 채 나온 힘과 응력은 조용히 틀린다.

\frac{\hbar^2}{2m}|\mathbf{k}+\mathbf{G}|^2 \leq E_{\text{cut}}

평면파의 미덕은 여전하다. 완전성이 단 하나의 스칼라로 제어되므로 수렴 검사가 단조롭고 깔끔하다. 물론 컷오프· $k$ 점·셀 크기 3종 수렴 검사를 안 하고 나온 숫자는 검증 및 확인의 관점에서 그냥 난수다.

4. 유사퍼텐셜과 k점[편집]

QE는 유사퍼텐셜을 코드에 동봉하지 않고 UPF 형식 파일을 사용자가 골라 오는 방식이다. 이게 자유이자 저주다. 같은 원소에 대해 수십 종의 UPF가 인터넷에 떠다니고, 어느 것을 고르냐에 따라 격자 상수가 갈린다.

종류	컷오프 요구	정확도	비고
노름보존(NC)	큼	높음	이론적으로 깔끔, 하이브리드·GW에 선호
울트라소프트(US)	작음	양호	노름 조건을 버리고 증강 전하 도입
PAW	작음	높음	전전자 정보 복원, 실무 기본값

이 혼돈을 정리하려는 시도가 검증된 유사퍼텐셜 라이브러리다. SSSP(Standard Solid State Pseudopotentials), PSlibrary, PseudoDojo, GBRV 등이 원소별로 “이걸 쓰고 컷오프는 이만큼 주라”는 표까지 제공한다. 고민할 시간에 SSSP 권장값부터 쓰는 것이 국룰.³

브릴루앙 영역 적분은 몬크호스트-팩(Monkhorst-Pack) 자동 격자로 지정한다. 금속은 페르미면에서 점유수가 불연속이라 성긴 격자로는 에너지가 요동치므로, 스미어링(마르차리-반더빌트 콜드 스미어링이 금속의 단골)이나 사면체법으로 완충해야 한다. 절연체는 스미어링 없이 고정 점유수로 두면 된다.

5. 병렬화[편집]

QE는 병렬화 축이 여러 겹이라는 것이 특징이다. 실행 시 플래그로 계층을 직접 쌓는다.

이미지(-ni) — NEB의 각 이미지, 포논의 각 $q$ 점처럼 서로 독립인 계산을 통째로 분배. 확장성이 거의 선형.
풀(-nk) — $k$ 점을 그룹에 나눠 준다. 통신량이 적어 가성비가 가장 좋은 축.
밴드(-nb) / 평면파 — 풀 내부에서 밴드와 $\mathbf{G}$ 벡터를 쪼갠다. 평면파 분산은 고속 푸리에 변환의 전역 통신을 유발해, 노드가 늘수록 효율이 꺾인다.
여기에 노드 안에서는 OpenMP 스레드가 붙는다.

$k$ 점 8개짜리 계산에 128코어를 -nk 없이 던지면 FFT 통신에 다 잡아먹힌다. 병렬 컴퓨팅의 교훈이 그대로 재현되는 현장. 근래 버전은 GPU 오프로딩도 지원한다.

6. VASP와의 대비[편집]

비용 — QE는 GPL, VASP는 유료 그룹 라이선스. 학생 여럿에게 계정을 뿌려야 하는 연구실에서 이 차이는 작지 않다.
소스 접근 — QE는 열려 있어 새 범함수나 물성 계산 루틴을 직접 얹을 수 있다. 방법론 개발자가 QE를 선호하는 이유.
DFPT — QE의 ph.x는 선형 응답으로 임의의 $q$ 점 포논을 초격자 없이 계산한다. 유한 변위(초격자) 방식에 의존하는 워크플로 대비 구조적 강점이다.
생태계 — QE는 AiiDA 워크플로 엔진, Materials Cloud와 밀착해 재현 가능한 자동화가 잘 갖춰져 있다. 반대로 VASP는 상용 지원과 “리뷰어가 익숙해한다”는 관성이 강점.
안정성 체감 — QE는 SCF가 잘 안 붙을 때 믹싱 파라미터(mixing_beta)를 손으로 낮춰가며 달래야 하는 경우가 상대적으로 잦다는 평이 있다.

정확도 자체는 둘 다 같은 이론(콘-샴 방정식 + 교환-상관 범함수)을 푸는 코드이므로, 동일한 유사퍼텐셜 품질과 동일한 수렴 조건에서는 유의미한 차이가 나지 않는 것이 정상이다.⁴

7. 관련 문서[편집]

Giannozzi et al., J. Phys.: Condens. Matter 21, 395502 (2009). 2017년에 후속 논문(J. Phys.: Condens. Matter 29, 465901)이 나왔다. QE로 논문을 쓰면 이 둘을 인용하는 것이 예의다. ↩
1 Ry ≈ 13.6057 eV. 실무에서 흔히 쓰는 ecutwfc는 30~80 Ry 범위이며, 산소나 전이금속이 끼면 위쪽이다. eV 감각으로 숫자를 적었다가 큐 대기시간까지 통째로 날려본 사람이 한둘이 아니다. ↩
SSSP는 원소마다 여러 유사퍼텐셜을 델타 인자(Δ-factor)와 포논·응력 벤치마크로 걸러 “efficiency”와 “precision” 두 세트를 제공한다. 남이 다 해놓은 수렴 검사를 안 쓸 이유가 없다. ↩
실제로 여러 전전자·평면파 코드를 같은 문제에 물려 비교한 대규모 재현성 연구가 있었고, 잘 만든 코드들끼리는 상태방정식이 거의 겹친다는 결론이 나왔다. 코드 탓하기 전에 대개 입력이 범인이다. ↩