We apply the well known Rayleigh-Ritz method (RRM) to the projection of a Hamiltonian operator chosen recently for the extension of the Rayleigh-Ritz variational principle to ensemble states. By means of a toy model we show that the RRM eigenvalues approach to those of the projected Hamiltonian from below in most cases but a few ones. We also discuss the effect of an energy shift and the projection of the identity operator.
논문 ID : 2411.14490제목 : 투영 해밀토니안에 대한 레일리-리츠 방법의 적용저자 : Francisco M. Fernández (INIFTA, DQT, La Plata, Argentina)분류 : quant-ph (양자물리학)발표 시간 : 2024년 12월 31일 (arXiv v5 버전)논문 링크 : https://arxiv.org/abs/2411.14490 본 논문은 고전적인 레일리-리츠 방법(RRM)을 투영 해밀토니안 연산자에 적용한다. 이 연산자는 최근 레일리-리츠 변분 원리를 앙상블 상태로 확장하기 위해 사용되었다. 장난감 모델을 통해 대부분의 경우 RRM 고유값이 투영 해밀토니안의 고유값에 아래에서 수렴함을 보여주지만, 소수의 예외도 존재함을 보인다. 본 논문은 또한 에너지 이동과 항등 연산자 투영의 영향을 논의한다.
전통적 RRM의 한계 : 고전적 레일리-리츠 방법은 원자 및 분자 전자 구조 연구에서 가장 널리 사용되는 방법 중 하나이며, 주요 장점은 RRM 고유값이 물리 시스템의 정확한 에너지에 위에서 수렴한다는 것이다.새로운 이론 발전 : Ding 등은 최근 레일리-리츠 변분 원리를 앙상블 상태로 확장하는 방법을 제안했으며, D차원 힐베르트 공간 위의 특수한 해밀토니안 연산자를 도입했다.이론적 공백 : 이 연산자가 실제 물리 문제의 처리에는 비현실적이지만, 그 수학적 성질, 특히 RRM이 투영 해밀토니안에 적용될 때의 거동은 추가 연구의 가치가 있다.투영 해밀토니안에서 RRM의 수렴 특성을 탐구하고, 특히 위에서 수렴하는 고전적 특성이 유지되는지 여부를 파악 투영 연산자가 변분 방법의 수렴 거동에 미치는 영향을 이해 확장된 변분 원리에 대한 이론적 기초 제공 투영 해밀토니안에 대한 RRM 적용을 처음으로 체계적으로 연구 하여 고전 이론과 다른 수렴 거동을 발견대부분의 경우 RRM 고유값이 아래에서 수렴함을 증명 하며, 이는 전통적인 위에서의 수렴과 대조됨구체적인 수치 예제 제공 하여 1차원 상자 속 입자 모델을 통해 이론적 예측을 검증에너지 이동과 항등 연산자 투영의 영향을 분석 하여 이론 체계를 완성N-D개의 영 고유값이 정확히 주어지는 흥미로운 현상 발견 레일리-리츠 방법이 투영 해밀토니안 연산자에 적용될 때의 수렴 특성을 연구하며, 투영 해밀토니안은 다음과 같이 정의된다:
H D = ∑ k = 0 D − 1 E k ∣ ψ k ⟩ ⟨ ψ k ∣ H_D = \sum_{k=0}^{D-1} E_k |\psi_k\rangle \langle\psi_k| H D = ∑ k = 0 D − 1 E k ∣ ψ k ⟩ ⟨ ψ k ∣
입력: 비직교 기저 함수 집합 { ∣ u j ⟩ } \{|u_j\rangle\} { ∣ u j ⟩}
출력: 근사 고유값 및 고유벡터
제약: 기저 함수가 투영 부분공간을 포함하는 더 큰 힐베르트 공간을 생성
무한 차원 힐베르트 공간 H \mathcal{H} H 위의 해밀토니안 연산자 H H H 를 고려하고, D차원 부분공간 H D \mathcal{H}_D H D 위의 투영을 다음과 같이 정의한다:
H D = ∑ k = 0 D − 1 E k ∣ ψ k ⟩ ⟨ ψ k ∣ H_D = \sum_{k=0}^{D-1} E_k |\psi_k\rangle \langle\psi_k| H D = ∑ k = 0 D − 1 E k ∣ ψ k ⟩ ⟨ ψ k ∣
여기서 { ∣ ψ k ⟩ , k = 0 , 1 , . . . , D − 1 } \{|\psi_k\rangle, k=0,1,...,D-1\} { ∣ ψ k ⟩ , k = 0 , 1 , ... , D − 1 } 는 부분공간 H D ⊂ H \mathcal{H}_D \subset \mathcal{H} H D ⊂ H 를 생성한다.
RRM은 행렬 원소를 계산해야 한다:
( H D ) i j = ∑ k = 0 D − 1 E k ⟨ u i ∣ ψ k ⟩ ⟨ ψ k ∣ u j ⟩ (H_D)_{ij} = \sum_{k=0}^{D-1} E_k \langle u_i|\psi_k\rangle \langle\psi_k|u_j\rangle ( H D ) ij = ∑ k = 0 D − 1 E k ⟨ u i ∣ ψ k ⟩ ⟨ ψ k ∣ u j ⟩
세속 방정식은: ∣ H − W S ∣ = 0 |H-WS| = 0 ∣ H − W S ∣ = 0
여기서 H i j = ⟨ u i ∣ H D ∣ u j ⟩ H_{ij} = \langle u_i|H_D|u_j\rangle H ij = ⟨ u i ∣ H D ∣ u j ⟩ , S i j = ⟨ u i ∣ u j ⟩ S_{ij} = \langle u_i|u_j\rangle S ij = ⟨ u i ∣ u j ⟩
수렴성 분석 : 고전적 RRM과 다른 수렴 패턴, 특히 아래에서의 수렴 현상 발견경계 조건 분석 : 코시-슈바르츠 부등식을 통해 특수한 경우의 상한과 하한 분석영 고유값의 정확성 : N-D개의 영 고유값이 정확히 주어짐을 증명1차원 상자 속 입자 모델 사용:
H = − 1 2 d 2 d x 2 H = -\frac{1}{2}\frac{d^2}{dx^2} H = − 2 1 d x 2 d 2
경계 조건: ψ ( 0 ) = ψ ( 1 ) = 0 \psi(0) = \psi(1) = 0 ψ ( 0 ) = ψ ( 1 ) = 0
정확한 해:
고유값: E k = k 2 π 2 2 , k = 1 , 2 , . . . E_k = \frac{k^2\pi^2}{2}, k = 1,2,... E k = 2 k 2 π 2 , k = 1 , 2 , ... 고유 함수: ψ k ( x ) = 2 sin ( k π x ) \psi_k(x) = \sqrt{2}\sin(k\pi x) ψ k ( x ) = 2 sin ( kπ x ) 비직교 기저 함수 집합 채택:
u i ( x ) = x i ( 1 − x ) , i = 1 , 2 , . . . u_i(x) = x^i(1-x), i = 1,2,... u i ( x ) = x i ( 1 − x ) , i = 1 , 2 , ...
투영 차원: D = 1, 2, 3, 4 기저 함수 개수: N = 1 ~ 19 에너지 이동: c = -5 표 1은 고전적 RRM의 상한 수렴 특성을 보여준다:
N=1일 때: W 1 = 5 > E 1 = 4.934802200 W_1 = 5 > E_1 = 4.934802200 W 1 = 5 > E 1 = 4.934802200 N이 증가함에 따라 모든 고유값이 위에서 정확한 값으로 수렴 D=1인 경우 (표 2) :
N=1: W 1 = 4.927671482 < E 1 = 4.934802200 W_1 = 4.927671482 < E_1 = 4.934802200 W 1 = 4.927671482 < E 1 = 4.934802200 아래에서 정확한 값으로 수렴 D=2인 경우 (표 3) :
두 개의 영이 아닌 고유값이 모두 아래에서 수렴 나머지 고유값은 정확히 영 D=3인 경우 (표 4) :
N=1일 때: W 1 = 4.988506932 > E 1 W_1 = 4.988506932 > E_1 W 1 = 4.988506932 > E 1 (상한) N>1일 때: 모든 고유값이 아래에서 수렴 표 6은 E k = k 2 π 2 / 2 − 5 E_k = k^2\pi^2/2 - 5 E k = k 2 π 2 /2 − 5 인 경우를 보여준다:
W 1 W_1 W 1 이 상한이 됨W 2 W_2 W 2 와 W 3 W_3 W 3 는 하한 특성 유지표 5는 기댓값 ⟨ H D ⟩ \langle H_D\rangle ⟨ H D ⟩ 이 D에 따라 변하는 것을 보여준다:
D가 증가함에 따라 기댓값이 ⟨ H ⟩ = 5 \langle H\rangle = 5 ⟨ H ⟩ = 5 로 수렴 해석적 증명: lim D → ∞ ⟨ H D ⟩ = 5 \lim_{D\to\infty}\langle H_D\rangle = 5 lim D → ∞ ⟨ H D ⟩ = 5 수렴 방향의 변화 : 고전적 RRM과 달리 투영 해밀토니안의 RRM 고유값은 주로 아래에서 수렴영 고유값의 정확성 : N-D개의 영 고유값이 정확히 주어짐차원 의존성 : 수렴 거동은 투영 차원 D와 기저 함수 개수 N의 관계에 따라 달라짐MacDonald (1933): RRM 상한 수렴의 이론적 기초 확립 Pilar, Szabo & Ostlund: 양자 화학에서의 응용 Ding 등 (2024): 앙상블 상태의 변분 원리 확장 제안 본 저자의 이전 연구: RRM의 일반적 논의 기존 연구와 비교하여 본 논문이 처음으로:
투영 해밀토니안의 RRM 거동을 체계적으로 연구 아래에서의 수렴이라는 새로운 현상 발견 구체적인 수치 검증 제공 수렴 특성의 변화 : RRM이 투영 해밀토니안에 적용될 때, 고유값은 주로 아래에서 수렴하며, 이는 고전적인 상한 특성과 대조된다.정확한 영 고유값 : RRM은 N-D개의 정확한 영 고유값을 제공하며, 이는 근사 해가 N→∞일 때만 정확한 해로 수렴하더라도 성립한다.조건 의존성 : 수렴 방향은 구체적인 투영 차원과 기저 함수 선택에 따라 달라진다.모델의 비물리성 : 연구된 투영 해밀토니안은 실제 물리 문제에서 비현실적이다.기저 함수 제한 : 결과는 특정 기저 함수 선택에 의존한다.이론적 설명의 불완전성 : 아래에서의 수렴 현상에 대한 엄격한 수학적 증명이 부족하다.투영 해밀토니안 RRM 수렴성의 엄격한 이론 개발 실제 물리 시스템에서의 응용 가능성 탐구 다른 유형의 투영 연산자의 특성 연구 이론적 혁신 : RRM의 투영 해밀토니안에서의 새로운 특성을 밝혀내어 전통적 인식에 도전충분한 수치 검증 : 여러 차원과 매개변수의 체계적 계산을 통해 이론적 예측 검증심층적 분석 : 수치 결과뿐만 아니라 해석적 분석과 이론적 설명 제공명확한 작성 : 논리 구조가 명확하고 수학적 유도가 엄밀함제한된 실용성 : 연구된 투영 해밀토니안은 실제 물리적 의미가 부족하다.불완전한 이론 : 아래에서의 수렴 현상에 대한 완전한 수학적 증명이 부족하다.좁은 적용 범위 : 결과는 특정 모델과 기저 함수 선택에 매우 의존적이다.이론적 기여 : 변분 방법 이론에 새로운 관점과 이해 제공방법론적 의의 : 새로운 수치 방법 개발에 영감을 줄 수 있음재현성 : 상세한 수치 결과를 제공하여 검증 및 확장 용이이론 양자역학 연구 변분 방법의 수학적 분석 투영 연산자 이론의 발전 수치 방법의 벤치마크 테스트 MacDonald, J. K. L. (1933). Successive approximations by the Rayleigh-Ritz variation method. Physical Review, 43, 830. Ding, L., Hong, C-L., & Schilling, C. (2024). Quantum 8, 1525. Fernández, F. M. (2024). On the Rayleigh-Ritz variational method. arXiv:2206.05122. 본 논문은 상대적으로 추상적인 수학 문제를 다루고 있지만, 투영 해밀토니안에서 RRM의 새로운 특성을 밝혀낸 것은 중요한 이론적 가치를 가지며, 변분 방법의 수렴 거동을 이해하기 위한 새로운 관점을 제공한다.