대규모 DFT 프로그램 CONQUEST의 다중 사이트 로컬 오비탈 최적화

초록

본 논문은 CONQUEST 코드에서 다중 사이트 서포트 함수(multi‑site support functions)를 수치적으로 최적화하는 방법을 제시한다. 기존에는 로컬 필터 대각화(LFD) 기법으로 계수를 얻었지만, 최적화 과정을 추가함으로써 작은 절단 반경에서도 정확도가 크게 향상되고, 에너지 최소화가 변분 원칙을 만족함을 확인하였다. 실험은 실리콘 결정, 벤젠 분자, 그리고 수화 DNA 시스템을 대상으로 수행되었다.

상세 분석

CONQUEST는 선형 스케일 DFT를 구현하기 위해 지역화된 지원 함수(support functions)를 사용한다. 전통적인 지원 함수는 각 원자에 할당된 의사 원자 궤도(pseudo‑atomic orbitals, PAOs)를 그대로 사용하거나, 다중 원자에 걸쳐 확장된 다중 사이트 지원 함수로 구성한다. 다중 사이트 지원 함수는 목표 원자와 일정 거리 이내에 있는 이웃 원자들의 PAO를 선형 결합해 만든다. 이 접근법은 최소한의 기반 함수 수로도 큰 기반의 정확도를 재현할 수 있다는 장점이 있다. 그러나 초기 계수 결정에 LFD 방법을 적용하면, 절단 반경이 작을 경우 계수가 최적이 아니어서 전자 밀도와 총 에너지에 오차가 누적된다.

본 연구는 이러한 한계를 극복하기 위해 LFD로 얻은 초기 계수를 출발점으로, 변분 원칙에 기반한 수치 최적화 절차를 도입한다. 구체적으로는, 각 다중 사이트 지원 함수에 대해 계수 벡터를 변수로 두고, 전체 에너지에 대한 그라디언트를 계산한 뒤, 제한된 라인 서치와 BFGS와 같은 준뉴턴 방법을 이용해 최소화한다. 최적화 과정에서 사용되는 에너지 함수는 Kohn‑Sham 에너지와 동일하며, 전자 밀도와 전위가 자체적으로 업데이트된다.

실험 결과는 세 가지 시스템에 대해 상세히 제시된다. 첫째, 실리콘 결정에서는 절단 반경을 5.0 Å에서 3.0 Å로 감소시켜도, 최적화된 다중 사이트 지원 함수는 기존 전통적 PAO 기반 계산과 비교해 격자 상수와 밴드 구조에서 차이가 10⁻⁴ eV 수준으로 감소한다. 둘째, 벤젠 분자에서는 결합 길이와 전자 밀도 차이가 LFD만 사용했을 때보다 약 5배 이상 개선되었다. 셋째, 수화 DNA 시스템(수천 원자 규모)에서는 최적화 전후의 총 에너지 차이가 0.02 eV 이하로 수렴했으며, 메모리 사용량과 계산 시간도 크게 증가하지 않아 실용성을 입증했다.

또한, 최적화가 변분 원칙을 만족한다는 점을 확인하기 위해, 최적화 전후의 에너지 곡선을 동일한 구조 변형에 대해 비교하였다. 최적화 후에는 에너지 감소가 항상 구조 변형에 대한 미분(힘)과 일치했으며, 이는 수치 최적화가 물리적으로 일관된 결과를 제공함을 의미한다.

이러한 결과는 다중 사이트 지원 함수가 작은 절단 반경에서도 높은 정확도를 유지하도록 하는 새로운 패러다임을 제시한다. 특히, 대규모 시스템에서 메모리와 계산 비용을 최소화하면서도 변분적 안정성을 확보할 수 있다는 점은 CONQUEST와 같은 선형 스케일 DFT 코드의 적용 범위를 크게 확대한다.