[논문 리뷰] The Exact Form of the Green's Function of the H\"uckel (Tight Binding) Model

Ramis Movassagh, Yuta Tsuji외 1명|arXiv (Cornell University)|2014. 07. 17.

Graph theory and applications인용 수 3

한 줄 요약

이 논문은 선형 체인과 순환 시스템에서 힉켈(H"uckel) 모델의 정확한 그린 함수의 형태를 해석적으로 유도하며, 역행렬이 실수 대칭 행렬이며, 짝수성과 거리에 기반한 폐쇄형 표현식을 따른다는 것을 증명한다. 고리 구조에서는 역행렬이 토플리츠 행렬임을 보이며, d차원 격자에서의 가역성 조건을 수론적 조건으로 설정한다. 이는 전자 이동도 및 전도도 계산에 응용된다.

ABSTRACT

The applications of the H\uckel (tight binding) model are ubiquitous in quantum chemistry and solid state physics. The matrix representation is isomorphic to an unoriented vertex adjacency matrix of a bipartite graph, which is also the Laplacian matrix plus twice the identity. In this paper, we analytically calculate the determinant and, when it exists, the inverse of this matrix in connection with the Green's function, $\mathbf{G}$, of the $N imes N$ H\uckel matrix for linear chains and cyclic systems. For an open linear chain we prove that $\mathbf{G}$ is a real symmetric matrix whose entries are $G\left(r,s ight)=\left(-1 ight)^{\frac{r+s-1}{2}}$ when $ $$r$ is even and $s<r$ is odd; $G\left(r,s ight)=0$ otherwise. A corollary is a closed form expression for a Harmonic sum (Eq. 9). For a ring we calculate the inverse, give formulas for the entries and find that it is always a Toeplitz matrix. We then extend the results to $d-$dimensional lattices, whose linear size is $N$. The existence of the inverse becomes a question of number theory. We prove that the inverse exists if and only if $N+1$ and $d$ are odd and $d$ is smaller than the smallest divisor of $N+1$. We corroborate our results by numerical demonstrations of the entry patterns of the Green's function and discuss applications related to transport and conductivity.

연구 동기 및 목표

• 일차원 선형 체인과 순환 시스템에서 힉켈 행렬의 정확한 해석적 형태를 유도하는 것.
• 특히 고차원 격자에서 힉켈 행렬의 역행렬이 존재하는 조건을 규명하는 것.
• 힌켈 행렬의 가역성과 수학적 이론, 특히 시스템 크기에 대한 나눗셈 조건 간의 연관성을 설정하는 것.
• 전자적 성질, 예를 들어 전도도와 이동도를 직접 계산할 수 있도록 매트릭스 원소에 대한 폐쇄형 표현식을 제공하는 것.

제안 방법

• 힌켈 행렬은 이분 그래프의 인접행렬로 표현되며, 라플라스 행렬에 두 배의 항등행렬을 더한 것과 동형이다.
• 선형 체인과 고리 구조에 특화된 조합론적 및 대수적 기법을 사용하여 행렬식과 역행렬을 해석적으로 계산한다.
• 선형 체인의 경우, 짝수성 기반의 부호 규칙을 사용하여 그린 함수 원소를 유도한다: $ G(r,s) = (-1)^{(r+s-1)/2} $ 이며, $ r $ 가 짝수이고 $ s < r $ 가 홀수일 때 성립하며, 그 외의 경우는 0이다.
• 고환 시스템의 경우, 역행렬이 항상 토플리츠 행렬임을 보이며, 삼각함수 합과 대칭성 고려를 통해 원소를 표현할 수 있다.
• d차원 격자로의 분석을 확장하기 위해, 가역성 조건을 $ N+1 $ 과 차원 $ d $ 를 포함하는 수론 문제로 환원한다.
• 수치적 검증을 통해 그린 함수의 원소 패턴이 분석적 예측과 일치함을 확인한다.

실험 결과

연구 질문

• RQ1 열린 경계 조건을 가진 일차원 힉켈 모델에서 그린 함수의 정확한 해석적 형태는 무엇인가?
• RQ2 힌켈 모델에서 선형 체인과 고환 시스템 간의 역행렬의 구조는 어떻게 다를까?
• RQ3 d차원 격자에서 크기 $ N $ 의 힉켈 행렬의 역행렬이 존재하는 수론적 조건은 무엇인가?
• RQ4 고환 시스템에서 그린 함수 원소에 대한 폐쇄형 표현식을 도출할 수 있으며, 항상 토플리츠 행렬인가?
• RQ5 유도된 그린 함수는 이러한 시스템에서 전자 이동도 및 전도도에 어떤 영향을 미치는가?

주요 결과

• 열린 선형 체인의 경우, 그린 함수는 실수 대칭 행렬이며, 원소가 $ G(r,s) = (-1)^{(r+s-1)/2} $ 로 주어지며, $ r $ 가 짝수이고 $ s < r $ 가 홀수일 때 성립하며, 그 외의 경우는 0이다.
• 이 결과의 따름정리로, 논문의 식 (9)에 제시된 특정 조화합에 대한 폐쇄형 표현식이 도출된다.
• 고환 시스템에서는 힉켈 행렬의 역행렬이 항상 토플리츠 행렬이며, 원소는 사이트 색인의 거리와 짝수성에 따라 달라진다.
• d차원 격자에서, 역행렬이 존재하는 것은 $ N+1 $ 과 $ d $ 가 모두 홀수여야 하며, $ d $ 가 $ N+1 $ 의 최소 소인수보다 작아야 한다는 조건과 동치이다.
• 수치 시뮬레이션을 통해 그린 함수의 예측된 원소 패턴이 분석적 유도와 일치함을 확인하였다.
• 유도된 그린 함수를 통해 힉켈 모델 시스템에서 이동도 및 전도도 성질을 직접 계산할 수 있다.