[논문 리뷰] Numerical analysis of the planewave discretization of orbital-free and Kohn-Sham models Part I: The Thomas-Fermi-von Weizacker model

Éric Cancès, Rachida Chakir외 1명|ArXiv.org|2009. 09. 08.

Acoustic Wave Resonator Technologies참고 문헌 7인용 수 4

한 줄 요약

이 논문은 밀도-functional 이론의 엄밀한 비 궤도함수 모델인 주기적 톰슨-페르미-폰 바이츠러크(TFW) 모델에 대한 평면파 이산화의 사전 오차 추정을 제공한다. 스펙트럴 및 허수 스펙트럴 푸리에 방법에 대해 색소프 노름에서 최적 수렴 속도를 확립하며, 에너지 커파트와 격자 해상도가 증가함에 따라 이산 최소화자가 연속 해로 최적 순서로 수렴함을 증명한다.

ABSTRACT

We provide {\it a priori} error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the condensed phase. The TFW model is stricly convex with respect to the electronic density, and allows for a comprehensive analysis (Part I). This is not the case for the Kohn-Sham LDA model, for which the uniqueness of the ground state electronic density is not guaranteed. Under a coercivity assumption on the second order optimality condition, we prove in Part II that for large enough energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in the vicinity of any Kohn-Sham ground state, and that this minimizer is unique up to unitary transform. We then derive optimal {\it a priori} error estimates for both the spectral and the pseudospectral discretization methods.

연구 동기 및 목표

• 주기적 톰슨-페르미-폰 바이츠러크(TFW) 모델에 대한 평면파 이산화의 엄밀한 사전 오차 추정을 제공함: 이는 궤도함수 밀도-functional 이론의 기초 모델이다.
• 색소프 노름 하에서 스펙트럴 및 허수 스펙트럴(평면파 기반) 이산화 방법의 수렴 성질을 분석함.
• 에너지 커파트와 격자 해상도가 증가함에 따라 TFW 에너지 함수의 이산 최소화자에 대한 최적 수렴 속도를 확립함.
• 지역 밀도 근사(LDA) 하에서 더 복잡한 코른-샤무 모델의 분석을 위한 수학적 기초를 마련함: 이는 제2부에서 다루어진다.

제안 방법

• 단위 세포 Γ = [0, L)³ 과 쌍대 격자 R* = (2π/L)ℤ³ 를 사용하여 TFW 모델을 주기적 설정으로 기술함.
• 파동수 k ∈ R* 에 대해 e_k(x) = |Γ|⁻¹/² e^{ik·x} 인 평면파 기저 함수를 사용하며, 이는 L²_#(Γ, ℂ) 에서 정규직교 기저를 이룸.
• 에너지 커파트 N_c 까지의 푸리에 모드를 포함하는 실수 값 함수의 유한차원 부분공간 V_{N_c} ⊂ H^s_#(Γ) 를 정의함.
• 최적 근사의 정의를 위해 투영 연산자 Π_{N_c} 를 사용하고, 색소프 임베딩과 역부등식을 통해 오차 한계를 유도함.
• 직각 격자 G_{N_g} = (L/N_g)ℤ³ 에서 이산 푸리에 변환(FFT)을 사용하여 적분을 근사하고 TFW 함수에서의 비선형성을 다룸.
• 총 오차를 截斷(N_c) 과 적분(quadature, N_g) 성분으로 분해하여 오차 추정을 유도함. 연속 해와 이산 해의 차이에 대한 L² 및 H¹ 노름의 경계를 사용함.

실험 결과

연구 질문

• RQ1 TFW 모델에 대한 평면파 이산화의 색소프 노름에서 최적 수렴 속도는 무엇인가?
• RQ2 에너지 커파트와 격자 해상도가 증가함에 따라 스펙트럴 및 허수 스펙트럴(FFT 기반) 평면파 방법의 오차 감소는 어떻게 비교되는가?
• RQ3 TFW 에너지 함수의 이산 최소화자가 H²_#(Γ) 에서 진정한 최소화자로 수렴하는 조건은 무엇인가?
• RQ4 해의 정(regularity)과 이산화 매개변수 N_c 및 N_g 에 따라 수렴을 명시적 상수로 정량화할 수 있는가?
• RQ5 헤시안의 강제성(coercivity)이 이산 해의 유일성과 수렴성을 보장하는 데 어떤 역할을 하는가?

주요 결과

• 적절한 정규성 조건 하에서 연속 해와 그 평면파 근사 사이의 H¹_#(Γ) 노름 오차는 O(N_c^{-1/2} N_g^{-7/2}) 의 속도로 감소한다.
• N_c → ∞ 이고 N_g ≥ 4N_c + 1 이면 이산 최소화자 u_{N_c,N_g} 가 진정한 해 u 에 대해 H²_#(Γ) 에서 수렴하며, 이는 ℝ³ 에서 L^∞ 수렴을 의미한다.
• 충분히 큰 N_c 에 대해 이산 해 u_{N_c,N_g} 는 0 에서 일정 거리 이상 떨어져 있어, H^{m+1/2−ε}_#(Γ) 에서 비선형 항의 균일 유계성을 보장한다.
• 에너지 함수 E^{TFW}_{N_g} 의 수렴 속도는 최적이며, E^{TFW}_{N_g}(v_{N_c}) - E^{TFW}_{N_g}(u_{N_c,N_g}) 의 차이는 H¹_# 노름에서 양의 정부호 이차형식에 의해 아래로 유계이다.
• 충분히 큰 N_c 와 N_g ≥ 4N_c + 1 에서 이산 최소화자 u_{N_c,N_g} 는 이산 에너지 함수의 헤시안의 강제성으로 인해 유일하다.
• 이산 고유값 λ_{N_c,N_g} 가 연속 고유값 λ 로 수렴하는 것은 최적 순서 O(N_c^{-1/2} N_g^{-m}) 로 이루어지며, m ≥ 4 이다.