[논문 리뷰] Wavelets for Density-Functional Theory and Post-Density-Functional-Theory Calculations

Bhaarathi Natarajan, Mark E. Casida외 2명|arXiv (Cornell University)|2011. 10. 21.

NMR spectroscopy and applications참고 문헌 82인용 수 5

한 줄 요약

이 논문은 BIGDFT 코드에서 다우비치 웨이브렛의 다중해상도 및 적응적 성질을 활용하여 밀도-functional 이론(DFT)과 선형응답 시간의존 DFT(LR-TD-DFT)의 웨이브렛 기반 구현을 제시한다. 이는 기저함수 한계 수렴을 효율적으로 달성한다. 결과적으로 웨이브렛 기반 기초함수는 궤도 및 흡수 에너지 계산에서 가우시안형 궤도함수(GTOs)보다 수렴 속도가 빠르며, CO 분자의 계산 결과는 고정밀도 GTO 코드와 일치하며, 대규모 순서-N 전자구조 계산에 유망하다.

ABSTRACT

We give a fairly comprehensive review of wavelets and of their application to density-functional theory (DFT) and to our recent application of a wavelet-based version of linear-response time-dependent DFT (LR-TD-DFT). Our intended audience is quantum chemists and theoretical solid-state and chemical physicists. Wavelets are a Fourier-transform-like approach which developed primarily in the latter half of the last century and which was rapidly adapted by engineers in the 1990s because of its advantages compared to standard Fourier transform techniques for multiresolution problems with complicated boundary conditions. High performance computing wavelet codes now also exist for DFT applications in quantum chemistry and solid-state physics, notably the BigDFT code described in this chapter. After briefly describing the basic equations of DFT and LR-TD-DFT, we discuss how they are solved in BigDFT and present new results on the small test molecule carbon monoxide to show how BigDFT results compare against those obtained with the quantum chemistry gaussian-type orbital (GTO) based code deMon2k. In general, the two programs give essentially the same orbital energies, but the wavelet basis of BigDFT converges to the basis set limit much more rapidly than does the GTO basis set of deMon2k. Wavelet-based LR-TD-DFT is still in its infancy, but our calculations confirm the feasibility of implementing LR-TD-DFT in a wavelet-based code.

연구 동기 및 목표

• BIGDFT 코드에서 고체 상태 DFT 및 LR-TD-DFT를 위한 웨이브렛 기반 접근법을 개발하고 검증한다.
• 웨이브렛 기반 기초함수는 전통적인 가우시안형 궤도함수(GTO) 기초함수 집합보다 기초함수 한계 수렴 속도가 더 빠르다는 것을 입증한다.
• 웨이브렛 이산화를 사용한 LR-TD-DFT의 구현 가능성과 정확도를 평가하여 전자 흡수 에너지 및 진동자 강도를 계산한다.
• 탄소 산화물과 같은 기준 시스템에서 웨이브렛 기반 결과를 고정밀도 GTO 기반 계산(DEMON2K)과 비교한다.
• 웨이브렛을 대규모 순서-N 전자구조 시뮬레이션의 실용적이고 적응 가능한 대안으로 정립한다.

제안 방법

• Kohn-Sham 궤도함수와 밀도를 다중해상도 분석(MRA) 프레임워크 내에서 다우비치 웨이브렛을 사용하여 적응형 공간 해상도로 표현한다.
• 웨이브렛 격자 위에서 유사 유한차분 이산화를 적용하여 Kohn-Sham 방정식을 해결함으로써 운동에너지 및 위치에너지 항을 효율적으로 계산한다.
• 비국소적 효과함수 및 포아송 연산자를 웨이브렛 적응형 수치 기법을 사용하여 적용하며, 계산 효율성을 높이기 위해 커널 분리 기법을 활용한다.
• 선형응답 시간의존 DFT(LR-TD-DFT)를 Tamm-Dancoff 근사(TDA)를 통해 구현하여 흡수 에너지 및 진동자 강도를 계산한다.
• TD-DFT 계산에서 발생하는 대규모 희박 고유값 문제를 해결하기 위해 Krylov 부분공간 방법을 사용한다.
• 탄소 산화물 시스템을 시험 시스템으로 삼아 고정밀도 GTO 기반 계산(DEMON2K)과의 결과를 비교 검증한다.

실험 결과

연구 질문

• RQ1 분자의 시스템에서 웨이브렛 기반 DFT는 기존의 GTO 기초함수 집합보다 기초함수 한계 수렴 속도가 더 빠를 수 있는가?
• RQ2 웨이브렛 기반 LR-TD-DFT 계산은 기존의 확립된 GTO 기반 코드와 비교해 흡수 에너지 및 진동자 강도를 얼마나 정확하게 재현하는가?
• RQ3 대규모 시스템에 대해 웨이브렛 기반 DFT 및 TD-DFT의 계산 효율성과 확장성은 어떠한가?
• RQ4 웨이브렛의 적응형 다중해상도 성질은 복잡한 전자 밀도 특성 또는 경계 조건을 가진 시스템을 다루는 데 자연스럽게 유리한가?
• RQ5 정확도를 손상시키지 않고 고성능 순서-N 프레임워크 내에서 웨이브렛 기반 TD-DFT를 구현할 수 있는가?

주요 결과

• BIGDFT의 웨이브렛 기반 DFT는 DEMON2K의 GTO 기반 DFT보다 기초함수 한계 수렴 속도가 현저히 빠르며, 궤도 에너지 수준은 유사하다.
• 탄소 산화물의 경우, BIGDFT에서의 HOMO 에너지(−12.85 eV)는 DEMON2K 결과(−12.84 eV)와 매우 유사하여 안정적인 기본 상태 기술이 가능함을 시사한다.
• BIGDFT에서 계산한 첫 번째 단일 상태 흡수 에너지(S1, 8.95 eV)는 DEMON2K 결과(8.43 eV)와 양호한 일치를 보이지만, 실험값(8.4 eV)에 더 가까운 것은 DEMON2K 결과이다.
• BIGDFT의 진동자 강도(0.853)는 DEMON2K 결과(0.232)보다 뚜렷이 높은데, 이는 TD-LDA에서 이온화 연속체가 너무 낮게 설정된 데 기인한 것으로 보이며, 연속체 효과에 민감함을 시사한다.
• 비록 초기 단계에 있지만, BIGDFT 내 LR-TD-DFT의 구현은 실현 가능하며, 장애 요소가 없음을 확인하여 향후 고성능 순서-N 웨이브렛 기반 TD-DFT 코드 개발에 기여할 수 있다.
• 웨이브렛 접근법은 자연스럽게 적응형 해상도를 지원하여 복잡한 경계 조건을 가진 다중해상도 문제를 효율적으로 다룰 수 있으며, 이는 표준 푸리에 및 GTO 방법에 비해 핵심적인 장점이다.