Differences

This shows you the differences between two versions of the page.

--- 물리:bbgky_계층 [2022/04/21 17:24] – [푸아송 괄호 계산] jiwon
+++ 물리:bbgky_계층 [2022/04/22 17:04] – [BBGKY 계층] jiwon
@@ Line 26: / Line 26: @@
 이 계의 해밀토니안이 다음과 같이 주어진다고 생각해보자.
 $$H(\mathbf p, \mathbf q) = \sum_{i=1}^N \left[\frac{\vec p_i^2}{2m}+U(\vec q_i)\right]+\frac12\sum_{i=1}^N\sum_{j=1}^NV(\vec q_i-\vec q_j)$$
-$U(\vec q_i)$는 외부 퍼텐셜, $V(\vec q_i-\vec q_j)$는 두 입자 상호작용을 나타낸다. 현재 우리의 관심사는 $\rho_s$의 시간변화이기 때문에 위 해밀토니안을 세 부분으로 나누어보도록 하자.
+$U(\vec q_i)$는 외부 퍼텐셜, $V(\vec q_i-\vec q_j)$는 두 입자 상호작용을 나타낸다. 지금은 전체 $N$개의 입자 중 $s$개를 부분계로 선택해서 이들의 시간에 따른 분포를 볼 것이다. 이를 위해 전체 해밀토니안을 세 부분으로 나누어보자.
 $$H_s = \sum_{n=1}^s \left[\frac{\vec p_n^2}{2m}+U(\vec q_n)\right]+\frac12\sum_{n=1}^s\sum_{m=1}^sV(\vec q_n-\vec q_m)$$
 $$H_{N-s} = \sum_{i=s+1}^N \left[\frac{\vec p_i^2}{2m}+U(\vec q_i)\right]+\frac12\sum_{i=s+1}^N\sum_{j=s+1}^NV(\vec q_i-\vec q_j)$$
 $$H' = \sum_{n=1}^s\sum_{i=s+1}^NV(\vec q_n-\vec q_j)$$
-여기에 [[물리:리우빌 정리]]를 적용하면, $\rho_s$의 시간변화는
+$\rho_s$의 시간변화는 [[물리:리우빌 정리]]를 적용하면 다음과 같이 주어진다.
 $$\frac{\partial\rho_s}{\partial t} = \int\prod_{i=s+1}^Nd^3p_id^3q_i\frac{\partial\rho}{\partial t} = -\int\prod_{i=s+1}^Nd^3p_id^3q_i \{\rho,H_s+H_{N-s}+H'\}$$
-로 주어진다.
 ==== 푸아송 괄호 계산 ====
@@ Line 57: / Line 55: @@
 마찬가지로 푸아송 괄호를 풀어서 적어보면
 \begin{align*}
-	\int \prod_{i=s+1}^N d^3p_id^3q_i\sum_{j=1}^N\left[\frac{\partial\rho}{\partial\vec q_j}\cdot\frac{\partial H'}{\partial\vec p_j}-\frac{\partial\rho}{\partial\vec p_j}\cdot\frac{\partial H'}{\partial\vec q_j}\right]
+	-\int \prod_{i=s+1}^N d^3p_id^3q_i\sum_{j=1}^N\left[\frac{\partial\rho}{\partial\vec q_j}\cdot\frac{\partial H'}{\partial\vec p_j}-\frac{\partial\rho}{\partial\vec p_j}\cdot\frac{\partial H'}{\partial\vec q_j}\right]
 \end{align*}
 이다. $\partial H'/\partial\vec q_j=0$이 되고, $j\le s$일 때는
@@ Line 69: / Line 67: @@
 이므로, 이를 모두 모아서 적어보면
 \begin{align*}
-	&\int \prod_{i=s+1}^N d^3p_id^3q_i\sum_{j=1}^N\left[\frac{\partial\rho}{\partial\vec q_j}\cdot\frac{\partial H'}{\partial\vec p_j}-\frac{\partial\rho}{\partial\vec p_j}\cdot\frac{\partial H'}{\partial\vec q_j}\right]\\
+	&-\int \prod_{i=s+1}^N d^3p_id^3q_i\sum_{j=1}^N\left[\frac{\partial\rho}{\partial\vec q_j}\cdot\frac{\partial H'}{\partial\vec p_j}-\frac{\partial\rho}{\partial\vec p_j}\cdot\frac{\partial H'}{\partial\vec q_j}\right]\\
 	=&\int \prod_{i=s+1}^N d^3p_id^3q_i\left[\sum_{n=1}^s\frac{\partial\rho}{\partial\vec p_n}\cdot\sum_{j=s+1}^N\frac{\partial V(\vec q_n-\vec q_j)}{\partial\vec q_n} + \sum_{j=s+1}^N\frac{\partial\rho}{\partial\vec p_j}\cdot\sum_{n=1}^s\frac{\partial V(\vec q_j-\vec q_n)}{\partial\vec  q_j}\right]
 \end{align*}
-이 된다.
+이 된다. 두 번째 항 중에서 $j$인덱스를 하나 골라 계산해보면
+\begin{align*}
+	&\int d^3p_jd^3q_j\frac{\partial\rho}{\partial\vec p_j}\sum_{n=1}^s\cdot\frac{\partial V(\vec q_j-\vec q_n)}{\partial\vec q_j}\\
+	=&\int d^3q_j\left[\rho\sum_{n=1}^s\cdot\frac{\partial V(\vec q_j-\vec q_n)}{\partial\vec q_j}\right]_{\vec p_j\text{ at }\infty}-\int d^3q_jd^3p_j\rho\sum_{n=1}^s\cdot\frac{\partial^2 V(\vec q_j-\vec q_n)}{\partial\vec p_j\partial\vec q_j}\\
+	=&0
+\end{align*}
+가 된다. 그리고 모든 입자는 같다고 생각하기 때문에 첫 항의 $j$에 대한 합은 $s+1$	번째 입자를 $N-s$번 고려한 것으로 생각할 수 있다. 정리해보면,
+\begin{align*}
+	&(N-s)\int\prod_{i=s+1}^Nd^3q_id^3p_i \sum_{i=1}^s \frac{\partial V(\vec q_n-\vec q_{s+1})}{\partial\vec q_n}\cdot\frac{\partial\rho}{\partial\vec p_n}\\
+	=&(N-s)\sum_{n=1}^s\int d^3q_{s+1}d^3p_{s+1}\frac{\partial V(\vec q_n-\vec q_{s+1})}{\partial\vec q_n}\cdot\frac{\partial\rho_{s+1}}{\partial\vec p_n}
+\end{align*}
+이다. 이제 계산해두었던 푸아송 괄호들을 다 모아보면
+$$
+	\frac{\partial\rho_s}{\partial t}-\{H_s,\rho_s\}=(N-s)\sum_{n=1}^s\int d^3q_{s+1}d^3p_{s+1}\frac{\partial V(\vec q_n-\vec q_{s+1})}{\partial\vec q_n}\cdot\frac{\partial\rho_{s+1}}{\partial\vec p_n}
+$$
+를 얻는다.