[논문 리뷰] Semi-classical limit of the Levy-Lieb functional in Density Functional Theory

Mathieu Lewin|arXiv (Cornell University)|2017. 06. 01.

Random Matrices and Applications참고 문헌 16인용 수 50

한 줄 요약

이 논문은 Bindini-De Pascale 정규화를 혼합 양자 페르미 상태로 확장하여 DFT에서 Levy-Lieb 함수의 준고전적 한계를 증명하고, 이를 쿨롱 비용을 갖는 다중-마르건 최적 수송과 연결한다.

ABSTRACT

In a recent work, Bindini and De Pascale have introduced a regularization of $N$-particle symmetric probabilities which preserves their one-particle marginals. In this short note, we extend their construction to mixed quantum fermionic states. This enables us to prove the convergence of the Levy-Lieb functional in Density Functional Theory , to the corresponding multi-marginal optimal transport in the semi-classical limit. Our result holds for mixed states of any particle number $N$, with or without spin.

연구 동기 및 목표

• DFT에서 1-입자 밀도를 보존하면서 N-입자 페르미 상태를 정규화할 필요성을 제시합니다.
• 기존의 보손 정규화를 혼합 페르미 상태로 확장하여 엄밀한 준고전적 분석을 가능하게 합니다.
• 준고전적 한계에서 Levy-Lieb 함수가 다중-마르건 쿨롱 수송으로 수렴하는 것을 확립합니다.

제안 방법

• 대칭 N-입자 밀도를 정규화하되 1-입자 밀도 rho_P를 보존하는 양자 확장 Gamma_epsilon를 도입합니다.
• χ_epsilon와 Slater 행렬식을 이용해 반대대칭을 강제하는 혼합 페르미 상태로 Gamma_epsilon를 구성합니다.
• Tr(-Δ) Gamma_epsilon가 N배의 H1 운동에너지의 제곱근 rho_P의 제곱의 적분과 χ의 ∇의 제곱의 (1/ε^2) 항의 합을 합산한 식과 같다는 것을 보임으로써 (1.6)을 확립합니다.
• Φ에 대한 기대값이 P에 따라 고전적 기대값으로 수렴함을 보이고 명시적 오차 한계를 제시합니다(1.7).
• 정의된 구성을 이용해 Levy-Lieb 함수의 DFT에 대한 준고전적 한계 결과를 도출합니다.

실험 결과

연구 질문

• RQ1 Bindini-De Pascale 유형의 정규화를 보손 상태에서 혼합 페르미 상태로 확장해도 같은 1-입자 밀도를 유지할 수 있는가?
• RQ2 확장된 정규화가 운동에너지의 제어를 가능하게 하고 준고전(η -> 0) 한계에서 Levy-Lieb 함수의 수렴을 보장하는가?
• RQ3 Levy-Lieb 함수의 준고전적 한계를 쿨롱 다중-마르건 수송 문제로 특징지을 수 있는가?
• RQ4 양자 혼합 상태에서 준고전적 수송 한계로의 전이를 할 때의 정량적 오차 추정은 무엇인가?

주요 결과

• 저자는 1-입자 밀도 P를 보존하는 페르미 혼합 상태 확장 Gamma_epsilon를 정의합니다.
• 그들은 Tr(-Δ) Gamma_epsilon가 N(∫|∇√ρ_P|^2 + (1/ε^2) ∫|∇χ|^2)임을 보이는 정확한 운동에너지 항등식을 증명합니다.
• Gamma_epsilon가 축소 밀도 및 기대값의 관점에서 고전적 N-입자 밀도 P로의 수렴을 확립합니다(축소 밀도에 의한 수렴 및 (1.7)을 통해).
• 준고전적 스케일링에서 Levy-Lieb 함수가 η로 나누어진 값이 다중-마르건 쿨롱 수송 에너지보다 크지 않으면서 그 에너지에 O(√η) 및 O(η) 항이 더해지는 구간에 끼여 있음을 보입니다.
• 주된 준고전적 결과는 η→0일 때 E_OT(ρ) ≤ E(η^3 ρ(η⋅))/η ≤ E_OT(ρ) + C(√η + η)로 연결되어 DFT를 다중-마르건 최적 수송에 연결합니다.