[논문 리뷰] On the Concept of Static Structure Factor

Kai Zhang|arXiv (Cornell University)|2016. 06. 11.

Mechanical Engineering and Vibrations Research참고 문헌 6인용 수 7

한 줄 요약

이 논문은 고체물리학에서 정적 구조 인자 $S(\mathbf{k})$에 둘러싸인 개념적이고 계산적 혼동을 명확히 하며, 표준 표현 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} (g(\mathbf{r}) - 1) e^{i\mathbf{k} \cdot \mathbf{r}}$가 올바르고 물리적으로 의미 있는 반면, $g(\mathbf{r})$의 비수렴 적분으로 인해 수학적으로 잘못 정의된 대안 표현 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} g(\mathbf{r}) e^{i\mathbf{k} \cdot \mathbf{r}}$는 잘못되었다고 주장한다. 시뮬레이션에서 직접 방법 $S(\mathbf{k}) = \frac{1}{N} \left\langle \left| \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i} \right|^2 \right\rangle$를 검증하여, 이 방법이 주기적 시스템에서 $S(\mathbf{k}) \sim N$ 스케일링을 정확히 반영하고 액체에서는 푸리에 변환 방법과 일치하지만, 결정 또는 조절된 상에서는 일치하지 않는다는 것을 보여준다.

ABSTRACT

We clarify the confusion in the expression of the static structure factor S(k) in the study of condensed matters and discuss its explicit form that can be directly used in calculations and computer simulations.

연구 동기 및 목표

• 고체물리계에서 정적 구조 인자 $S(\mathbf{k})$에 대한 표준 표현에 오랫동안 지속된 문헌 내 혼동을 해결하기 위해.
• 표준 표현 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} (g(\mathbf{r}) - 1) e^{i\mathbf{k} \cdot \mathbf{r}}$가 수렴 가능하고 올바르며, $g(\mathbf{r})$의 비수렴 적분으로 인해 수학적으로 잘못 정의된 표현 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} g(\mathbf{r}) e^{i\mathbf{k} \cdot \mathbf{r}}$는 비수렴함을 보여주기 위해.
• 모oleculer 시뮬레이션에서 신뢰할 수 있고 정확한 계산 방법으로서 직접 방법 $S(\mathbf{k}) = \frac{1}{N} \left\langle \left| \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i} \right|^2 \right\rangle$를 검증하기 위해.
• 등방성이고 균일한 시스템에 대해서만 유효한 푸리에 변환 방법 $S(k) = 1 + 4\pi\rho \int_0^\infty (g(r)-1) r^2 \frac{\sin(kr)}{kr} dr$가 조절된 또는 결정상에서 실패하는 이유를 밝히기 위해.
• $S(\mathbf{k}=0) = N$의 물리적 의미와 압축성 모듈러스 $\kappa$에 대한 $\lim_{\mathbf{k} \to 0} S(\mathbf{k}) = \rho k_B T \kappa$의 정확한 물리적 해석을 명확히 하기 위해.

제안 방법

• 밀도 변동 $\hat{\rho}_\mathbf{k} = \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i}$의 푸리에 변환으로부터 정적 구조 인자를 유도하여 직접 표현 $S(\mathbf{k}) = \frac{1}{N} \left\langle \left| \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i} \right|^2 \right\rangle$를 도출한다.
• 이 직접 방법이 주기적 구조(예: 결정, 레이어 구조)에서 $S(\mathbf{k}) \sim N$ 스케일링을 자연스럽게 반영함을 보여주며, 물리적 기대와 일치함을 확인한다.
• 무한대에서 $\lim_{r \to \infty} g(\mathbf{r}) = 1$이므로 $g(\mathbf{r})$의 푸리에 변환 $\int g(\mathbf{r}) e^{i\mathbf{k} \cdot \mathbf{r}} d\mathbf{r}$가 수렴하지 않아, $g(\mathbf{r})$ 기반의 표준 $S(\mathbf{k})$ 표현이 수학적으로 잘못 정의됨을 분석한다.
• 쌍체 상관 함수 $h(\mathbf{r}) = g(\mathbf{r}) - 1$을 도입하여, 그 푸리에 변환 $\hat{h}_\mathbf{k}$가 잘 수렴하며 $S(\mathbf{k}) = 1 + \rho \hat{h}_\mathbf{k} + N\delta_{\mathbf{k},0}$의 정확한 표현을 이끌어내는 방법을 제시한다.
• 역 푸리에 변환 관계 $g(\mathbf{r}) = 1 + \frac{1}{(2\pi)^3} \int d\mathbf{k} \frac{S(\mathbf{k}) - 1}{\rho} e^{-i\mathbf{k} \cdot \mathbf{r}}$를 유도하여, $S(\mathbf{k}) - 1$이 $\mathbf{k} \to 0$ 포함 모든 $\mathbf{k}$에 대해 통합되어야만 $g(\mathbf{r})$를 복원할 수 있음을 보여준다.
• Lennard-Jones 시스템을 액체, 레이어 구조, 결정 상에서 수치적으로 시험하여, 직접 방법(2)과 푸리에 변환 방법(6)을 비교함으로써, 액체에서는 일치하지만 주기적 상에서는 일치하지 않음을 확인한다.

실험 결과

연구 질문

• RQ1 장거리 상관이 있는 시스템에서 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} g(\mathbf{r}) e^{i\mathbf{k} \cdot \mathbf{r}}$ 표현이 수학적으로 문제가 되는 이유는 무엇인가요?
• RQ2 모의 실험 및 이론 계산에 사용하기 위한, 반경분포함수 $g(\mathbf{r})$를 기반으로 한 정적 구조 인자 $S(\mathbf{k})$의 정확하고 수렴 가능한 표현은 무엇인가요?
• RQ3 다른 상에서 $S(\mathbf{k})$를 캡처하는 데 있어 직접 방법 $S(\mathbf{k}) = \frac{1}{N} \left\langle \left| \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i} \right|^2 \right\rangle$와 $g(\mathbf{r})$의 푸리에 변환 간의 비교는 어떻게 되나요?
• RQ4 표준 푸리에 변환 방법 $S(k) = 1 + 4\pi\rho \int_0^\infty (g(r)-1) r^2 \frac{\sin(kr)}{kr} dr$가 조절된 또는 결정상에서 실패하는 이유는 무엇인가요?
• RQ5 $S(\mathbf{k}=0) = N$의 물리적 해석은 무엇이며, 이는 압축성 모듈러스 $\kappa$와 어떻게 관련이 있나요?

주요 결과

• $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} (g(\mathbf{r}) - 1) e^{i\mathbf{k} \cdot \mathbf{r}}$는 수학적으로 잘 정의되어 있으며 올바르며, 반면 $S(\mathbf{k}) = 1 + \rho \int d\mathbf{r} g(\mathbf{r}) e^{i\mathbf{k} \cdot \mathbf{r}}$는 무한대에서 $g(\mathbf{r})$의 비적분 가능성으로 인해 발산한다.
• 직접 방법 $S(\mathbf{k}) = \frac{1}{N} \left\langle \left| \sum_i e^{i\mathbf{k} \cdot \mathbf{r}_i} \right|^2 \right\rangle$는 주기적 시스템(예: 결정, 레이어)에서 $S(\mathbf{k}) \sim N$ 스케일링을 정확히 반영하며, 모의 실험에 매우 강력하고 정확하다.
• 밀도 $\rho\sigma^3 = 0.8$인 Lennard-Jones 액체에서, 직접 방법과 $g(r)$의 푸리에 변환 방법이 모든 $k$ 값에서 완벽하게 일치한다.
• 밀도 $\rho\sigma^3 = 0.4$인 레이어 상에서, 직접 방법은 $k^* = 0.36744$에서 $S(k^*) = 270.55$의 주요 피크를 드러내며, 이는 $N = 2000$에 비례하여 스케일링됨을 보여주며 장거리 질서를 나타낸다.
• 밀도 $\rho\sigma^3 = 1.2$인 면심정육면계 결정에서, $g(r)$의 푸리에 변환은 $S(k)$의 피크 높이를 정확히 재현하지 못하며, 이는 $N$에 비례해야 하지만 직접 방법은 정확히 이러한 스케일링을 반영한다.
• $S(\mathbf{k}) = 1 + \rho \hat{h}_\mathbf{k} + N\delta_{\mathbf{k},0}$ 표현은 $\mathbf{k} = 0$에서 $N$에 의존하는 발산을 정확히 반영하며, $\lim_{\mathbf{k} \to 0} S(\mathbf{k}) = \rho k_B T \kappa$는 물리적으로 의미 있는 탄성 모듈러스 한계이다.