플라스마 가림효과
플라스마 가림효과(Plasma screening)란 이온화 된 전자 즉, 플라스마가 들뜨게 되면 전자 간의 상호작용했던 쿨롱 퍼텐셜이 서로 가림을 받는 현상이 일어나게 되는데 이러한 현상을 가리켜 플라즈마 가림효과라 일컫는다. 플라즈마 가림현상으로 나타나는 대표적인 특징은 물질의 유전상수 변화이다. 주어진 퍼텐셜에서 플라즈마 가림현상으로 일어나게 되는 유전상수의 변화를 나타내는 식이 린드하드 식이다. 1954년 린드하드가 알아낸 이 유전함수는 흔히 RPA (Random Phase Approximation)이라고 많이 불린다. 이는 동적상태인 formula_1와 정적상태인 formula_2로 나누어 모델링 한 것인데 특히, 전자 가스와 같은 플라즈몬에 대해서 잘 예측하게 해준다. 그 린드하드식은 아래와 같다.
여기서, formula_3은 동적 유전상수(Dynamic Dielectric Funtion)이고 formula_4 는 페르미-디락 분포를 따르는 캐리어 분포 함수를 formula_5는 formula_6을 푸리에 변환한 함수이다.
어떤 시스템의 Dynamic Electronic Response를 보기 위한 방법으로 주로 사용되는 것이 Random Phase Approximation이다. RPA에서 전자는 외부 퍼텐셜 formula_7 과 가림효과를 받은 가림 퍼텐셜 formula_8을 합친 formula_6에 의해서만 반응을 하고 외부 퍼텐셜formula_7이 특정한 단일 주파수 formula_11로 진동한다고 가정을 한 방법이다. 따라서, 전체 퍼텐셜 formula_6이 동적 유전상수에 대한 기여는 평균화 되어 특정한 Wave Vector formula_13 만이 기여를 하게 된다.
전자 플라즈자에서 다체계 상호작용의 대표적인 특징이 플라즈마 가림효과이다. 플라즈마 가림효과는 양자역학적으로 접근을 해야한다. 단일 입자에 관한 해밀토니안을 아래와 같이 표현 할 수 있다. 여기에서 formula_14 이다. formula_15 은 쿨롱 퍼텐셜이고, formula_16 는 가림을 받는 입자들에 의해 유도된 퍼텐셜이다.
위에 있는 해밀토니안을 푸리에 변환을 하면 아래와 같은 식이 된다. (일반적으로 전자 플라즈마에 대해서 분석하기 위해서 Charge Density Fluctuation을 계산해야 한다. 그런데 Second Quantization Formalism을 이용하면 Charge Density Operatior의 기대값을 구할 수 있다. 따라서 formula_18 나 formula_19와 같은 Operator가 수식에 나오게 되는 것이다.)
이제, 하이젠버그 운동방정식을 이용하면 formula_21" "의 운동에 대해 기술 할 수 있게 된다.
앞의 결과와 Random Phase Approximation을 한다면, formula_21" "에 대한 기대값을 구할 수 있게 되고,
또한 전자가 " "formula_26" "에 비례해서 움직인다고 가정을 하면 즉, " "formula_27" " 이라고 생각 할 수 있다. 따라서 앞의 식을 바탕으로 다음과 같은 식을 유도할 수 있다.
결국, formula_29에 대해서 정리하면 아래와 같게 된다.
더 나아가기 이전에, 두 가지 함수를 정의해 보도록 하겠다.
먼저 Polarization Function으로 formula_31 으로 표현을 하며, 그 꼴을 아래와 같다.
다음으로 Electron Charge Denisty Operator이다.
따라서 위에서 정의한 두 식과 앞의 유도한 식을 이용하면 Electron Charge Denisty Operator은 아래와 같이 새롭게 정의가 된다.
가림효과를 받은 입자의 퍼텐셜은 포아송 방정식을 만족하므로 다음과 같이 쓰여 질 수 있다.
위 식을 푸리에 변환을 하고, 앞서 구한 Electron Charge Denisty Operator을 대입하면,
맨 처음에 formula_37이라고 했으므로, 위의 결과 들을 이용하면 formula_38를 가림을 받게되는 실질적인 퍼텐셜(Dynamically Screened Coulomb Potential)인 formula_39으로 새로 정의할 수 있다.
이 때,
formula_41
으로 formula_42 이 동적 유전상수(Dynamic Dielectric Funtion)로 정의 된다.
그러므로, 우리는 Polarization Function formula_31의 꼴을 알고 있으므로 린드하드식을 다음과 같이 유도할 수 있게 되는 것이다.
여기에 언급한 린드하드 식은 3차원과 2차원에 한 해서 적용이 가능하다. 그리고 formula_45" "와 formula_46" "에 대한 의존성이 있기 때문에 공간과 시간에 대한 분포를 포함한 식이라고 볼 수 있다. 식에서 사용된 입자의 밀도 함수 " "formula_47" "는 플라즈마의 경우에는 페르미-디락 함수" "formula_48" "와 같다. 바로 아래에서 차원에 따라서, 또 공간과 시간의 극한영역에서는 린드하드식이 어떻게 되는지 알아보도록 하겠다.
긴파장 영역은 formula_49으로 가능 영역이다. 왜냐하면, formula_50는 파수이기 때문에 formula_51가 커지면 0으로 가기 때문이다.
린드하드 식의 분모를 보면,
린드하드 식의 분자를 보면,
위 식을 린드하드 식에 대입을 하고 formula_54으로 가는 경우를 취하면,
여기서 formula_56, formula_57, formula_58 이다. 이는 고전적인 유전상수 모델(즉, Drude 모델)과 같은 결과이다.
정적인 상태는 formula_59 일 때 이므로, 린드하드 식은 다음과 같이 된다.
평형상태의 페르미-디락 분포상태를 가정하면,
여기서 formula_62, formula_63이다.
그러므로, 정적 상태인 3차원 린드하드식은 아래와 같이 유도된다.
여기서 formula_65를 새롭게 정의 했는데 이를 3D Screening Wave Number라고 한다.
formula_66
또한, 앞의 결과를 이용하여 3차원 일 때 가림효과를 고려한 새로운 쿨롱 포텐션을 구할 수 있다.
이 퍼텐셜은 주파수 도메인으로 표현되었으므로, 푸리에 변환을 해 공간 상에서 어떤 함수를 갖는지 알아보도록 하면,
이 되어 우리가 흔히 보았던 유카와 퍼텐셜이 됨을 알 수 있다. 두가지 퍼텐셜을 비교하면, 아래 그래프와 같다.
3D Screening Wave Number를 축퇴 상태와 비축퇴 상태 두 가지 경우에 대해서 생각해 볼 수 있다. 먼저 축퇴상태를 생각해보면, 온도가 0이 상태이므로 종종 이러한 상태의 가림을 Thomas-Fermi Screening이라고 한다. 축퇴상태의 전자기체는 다음과 같은 밀도함수를 가진다.
이 상태에서 페르미 에너지는 formula_70와 같아 지므로,
을 얻을 수 있다.
따라서 축퇴상태의 3D Screening Wave Number는 아래와 같이 된다. 이와 같은 상태를 Thomas-Fermi Screening이라고 하므로 3D Thomas-Fermi Screening Wave Number라고 부른다.
다음으로 비축퇴 상태를 생각해 보면, 페르미 분포로부터 다음 식을 유도할 수 있다.
따라서 아래와 같이 비축퇴 상태의 3D Screening Wave Number를 구할 수 있고, 이를 Debye-Hückel Screening Wave Number라고 부른다.
긴파장 영역은 formula_75으로 가능 영역이다. 왜냐하면, q는 파수이기 때문에 formula_76가 커지면 0으로 가기 때문이다.
린드하드 식의 분모를 보면,
린드하드 식의 분자를 보면,
위 식을 린드하드 식에 대입을 하고 formula_54으로 가는 경우를 취하면 원하는 결과를 얻을 수 있다.
여기서 formula_56, formula_82이고, 특히 2D plasma frequency를 formula_83 와 같이 정의하였다.
정적인 상태는 formula_59 일 때 이므로, 린드하드 식은 다음과 같이 된다.
평형상태의 페르미-디락 분포상태를 가정하면,
여기서 formula_62, formula_63이다.
따라서, 2D 유전상수를 구할 수 있다.
여기서 formula_90 로 정의를 했다. 3차원일 때와 마찬가지로 formula_65 를 알면 가림효과를 고려한 쿨롱 퍼텐셜을 쉽게 유도 할 수 있다. 그 결과는 아래와 같다.
2차원 페르미 기체의 Chemical Potenital은 formula_93 이므로 formula_94 가 된다.
따라서 formula_95 가 됨을 알 수 있다. 수식을 보면 2차원 Screening Wave Number는 3차원이 밀도에 의존하는데 반해, 그렇지 않다는 것을 확인할 수 있다.