개념완성

슬레이터 행렬식 (Slater Determinant)

1. 정의

NN개의 동일 페르미온의 반대칭 파동함수를 체계적으로 구성하는 방법이다.

정의8.5슬레이터 행렬식

NN개의 페르미온이 각각 1입자 상태 ψ1,ψ2,,ψN\psi_1, \psi_2, \ldots, \psi_N을 점유할 때, 반대칭 NN-입자 파동함수:

Ψ(x1,x2,,xN)=1N!ψ1(x1)ψ2(x1)ψN(x1)ψ1(x2)ψ2(x2)ψN(x2)ψ1(xN)ψ2(xN)ψN(xN)\Psi(\mathbf{x}_1, \mathbf{x}_2, \ldots, \mathbf{x}_N) = \frac{1}{\sqrt{N!}} \begin{vmatrix} \psi_1(\mathbf{x}_1) & \psi_2(\mathbf{x}_1) & \cdots & \psi_N(\mathbf{x}_1) \\ \psi_1(\mathbf{x}_2) & \psi_2(\mathbf{x}_2) & \cdots & \psi_N(\mathbf{x}_2) \\ \vdots & \vdots & \ddots & \vdots \\ \psi_1(\mathbf{x}_N) & \psi_2(\mathbf{x}_N) & \cdots & \psi_N(\mathbf{x}_N) \end{vmatrix}

여기서 xi=(ri,si)\mathbf{x}_i = (\mathbf{r}_i, s_i)는 위치와 스핀을 합친 좌표이며, 1/N!1/\sqrt{N!}은 규격화 인자이다.

2. 슬레이터 행렬식의 성질

참고핵심 성질

행렬식의 성질로부터 다음이 자동으로 보장된다:

  1. 반대칭성: 두 행(두 입자)을 교환하면 행렬식의 부호가 바뀐다. P^ijΨ=Ψ\hat{P}_{ij}\Psi = -\Psi 자동 만족.

  2. 파울리 배타 원리: 두 열(두 상태)이 같으면 행렬식이 0이 된다. 즉, 두 페르미온이 같은 양자 상태를 점유할 수 없다.

  3. 비분리성: 슬레이터 행렬식은 일반적으로 단순 곱 상태로 분해할 수 없다. 이는 페르미온의 양자 상관(correlation)을 나타낸다.

3. 명시적 예: 2-전자 시스템

예제헬륨 원자 기저 상태의 슬레이터 행렬식

헬륨 원자의 기저 상태에서 두 전자는 같은 궤도 ψ1s\psi_{1s}를 점유하되, 스핀이 반대이다.

ψ1(x)=ψ1s(r)α(s)\psi_1(\mathbf{x}) = \psi_{1s}(\mathbf{r})\alpha(s), ψ2(x)=ψ1s(r)β(s)\psi_2(\mathbf{x}) = \psi_{1s}(\mathbf{r})\beta(s)

슬레이터 행렬식:

Ψ(x1,x2)=12ψ1s(r1)α(s1)ψ1s(r1)β(s1)ψ1s(r2)α(s2)ψ1s(r2)β(s2)\Psi(\mathbf{x}_1, \mathbf{x}_2) = \frac{1}{\sqrt{2}} \begin{vmatrix} \psi_{1s}(\mathbf{r}_1)\alpha(s_1) & \psi_{1s}(\mathbf{r}_1)\beta(s_1) \\ \psi_{1s}(\mathbf{r}_2)\alpha(s_2) & \psi_{1s}(\mathbf{r}_2)\beta(s_2) \end{vmatrix}=12ψ1s(r1)ψ1s(r2)[α(s1)β(s2)α(s2)β(s1)]= \frac{1}{\sqrt{2}}\psi_{1s}(\mathbf{r}_1)\psi_{1s}(\mathbf{r}_2)[\alpha(s_1)\beta(s_2) - \alpha(s_2)\beta(s_1)]

공간 부분: 대칭 ψ1s(r1)ψ1s(r2)\psi_{1s}(\mathbf{r}_1)\psi_{1s}(\mathbf{r}_2)

스핀 부분: 반대칭 단일항 12(αββα)\frac{1}{\sqrt{2}}(\alpha\beta - \beta\alpha)

4. 리튬 원자의 예

예제리튬 원자의 기저 상태

리튬 (Z=3Z = 3)의 기저 상태 전자 배치: 1s22s11s^2 2s^1

세 개의 스핀-궤도: ψ1=ψ1sα\psi_1 = \psi_{1s}\alpha, ψ2=ψ1sβ\psi_2 = \psi_{1s}\beta, ψ3=ψ2sα\psi_3 = \psi_{2s}\alpha

Ψ=13!ψ1s(r1)α(1)ψ1s(r1)β(1)ψ2s(r1)α(1)ψ1s(r2)α(2)ψ1s(r2)β(2)ψ2s(r2)α(2)ψ1s(r3)α(3)ψ1s(r3)β(3)ψ2s(r3)α(3)\Psi = \frac{1}{\sqrt{3!}} \begin{vmatrix} \psi_{1s}(\mathbf{r}_1)\alpha(1) & \psi_{1s}(\mathbf{r}_1)\beta(1) & \psi_{2s}(\mathbf{r}_1)\alpha(1) \\ \psi_{1s}(\mathbf{r}_2)\alpha(2) & \psi_{1s}(\mathbf{r}_2)\beta(2) & \psi_{2s}(\mathbf{r}_2)\alpha(2) \\ \psi_{1s}(\mathbf{r}_3)\alpha(3) & \psi_{1s}(\mathbf{r}_3)\beta(3) & \psi_{2s}(\mathbf{r}_3)\alpha(3) \end{vmatrix}

이 행렬식을 전개하면 6개의 항이 나온다. 세 번째 전자의 스핀을 β\beta로 선택하면 동등한 에너지의 다른 슬레이터 행렬식을 얻는다.

5. 하트리-폭 방법

정의8.6하트리-폭 방법

하트리-폭(Hartree-Fock) 방법은 다전자 파동함수를 단일 슬레이터 행렬식으로 근사하되, 변분 원리를 적용하여 최적의 1입자 궤도를 결정하는 방법이다.

변분 조건 δE[{ψi}]=0\delta E[\{\psi_i\}] = 0 (규격화 조건 ψiψj=δij\langle\psi_i|\psi_j\rangle = \delta_{ij} 하에서)을 적용하면 하트리-폭 방정식을 얻는다:

f^ψi(x)=ϵiψi(x)\hat{f}\psi_i(\mathbf{x}) = \epsilon_i\psi_i(\mathbf{x})

여기서 폭 연산자 f^\hat{f}는:

f^=h^+j(J^jK^j)\hat{f} = \hat{h} + \sum_j(\hat{J}_j - \hat{K}_j)
  • h^\hat{h}: 1전자 해밀토니안 (운동 에너지 + 핵-전자 인력)
  • J^j\hat{J}_j: 쿨롱 연산자 (직접 전자-전자 반발)
  • K^j\hat{K}_j: 교환 연산자 (교환 효과)

6. 슬레이터 행렬식의 한계와 전자 상관

참고전자 상관 (Electron Correlation)

단일 슬레이터 행렬식은 교환 상관(exchange correlation)만 포함하고, 쿨롱 상관(Coulomb correlation, 또는 동적 상관)은 포함하지 않는다. 정확한 파동함수를 얻으려면 여러 슬레이터 행렬식의 선형 결합이 필요하다:

Ψexact=IcIΦI\Psi_{\text{exact}} = \sum_I c_I \Phi_I

여기서 ΦI\Phi_I는 서로 다른 전자 배치에 해당하는 슬레이터 행렬식이다. 이것이 배치 상호작용(configuration interaction, CI) 방법의 기본 아이디어이다.

상관 에너지:

Ecorr=EexactEHFE_{\text{corr}} = E_{\text{exact}} - E_{\text{HF}}

원자의 경우 상관 에너지는 전체 에너지의 1%\sim 1\% 정도이지만, 화학 결합 에너지와 비슷한 크기이므로 화학적 정확도를 위해서는 반드시 고려해야 한다.