[논문 리뷰] Anharmonic phonon frequency shift in MgB2

Michele Lazzeri, Matteo Calandra외 1명|arXiv (Cornell University)|2003. 06. 25.

Superconductivity in MgB2 and Alloys참고 문헌 29인용 수 63

한 줄 요약

이 연구는 고체 상태 이론에서 밀도함수이론 내에서 다체 섭동이론을 사용하여 MgB2의 비조화 격자파 주파수 이완을 계산하며, 브릴루아 존 전역에서 삼중 및 사중 격자파 산란을 명시적으로 고려한다. 주요 결과로는 Γ 점에서 E2g 모드의 +5.4% 이완이 도출되었으며, 이는 이전의 추정치보다 훨씬 낮은 값으로, 이론적 예측가 X선 산란 실험 결과와 일치시킨다.

ABSTRACT

We compute the anharmonic shift of the phonon frequencies in MgB2, using density functional theory. We explicitly take into account the scattering between different phonon modes at different q-points in the Brillouin zone. The shift of the E2g mode at the Gamma point is +5 % of the harmonic frequency. This result comes from the cancellation between the contributions of the four- and three-phonon scattering, respectively +10 % and -5 %. A similar shift is predicted at the A point, in agreement with inelastic X-ray scattering phonon-dispersion measurements. A smaller shift is observed at the M point.

연구 동기 및 목표

• MgB2의 E2g 격자파 모드에서 이론적 예측의 큰 비조화 이완과 비탄성 X선 산란 실험 데이터 사이의 오랫동안 지속된 괴리 문제를 해결하기 위해.
• 특히 브릴루아 존의 Γ 및 A 점에서 비조화 주파수 이완을 결정하는 데 기여하는 삼중 및 사중 격자파 산란 과정의 기여도를 정량화하기 위해.
• 모든 브릴루아 존에서 운동량 및 모드 의존성 산란을 명시적으로 포함함으로써 이전의 고정된 격자파 근사 모델을 개선하기 위해.
• 큰 비조화 이완을 암시하는 Raman 스펙트로스코피의 모순되는 실험적 해석과, 최소한의 이완을 보여주는 더 직접적인 X선 측정 결과 사이의 모순을 해소하기 위해.
• 다체 섭동이론을 사용하여 MgB2에서의 격자파 비조화성을 정량적으로 정확하게 기술하는 제1원리적 기술을 제공하기 위해.

제안 방법

• 총 에너지의 도함수를 사용하여 '2n+1' 정리에 기반해 제3 및 제4계력 상수를 계산하기 위해 밀도함수 섭동이론(DFPT)을 적용하였다.
• 자기에너지의 실수부가 주파수 이완을 제공하는 다항 섭동 전개를 사용하여 격자파 자기에너지(자기에너지의 실수부)를 계산하였으며, 이는 T(타플), L(루프), B(버블) 도형을 포함하였다.
• 다른 파장벡터(q)와 다른 모드(j)를 가진 격자파 간의 산란을 명시적으로 포함하여 고정된 격자파 모델의 한계를 극복하였다.
• 힘 상수 계산을 위해 1×1×1 초세포와 4×4×2 격자 메쉬를 사용하였으며, 브릴루아 존에 대한 합산을 위해 20×3 격자를 사용하여 수렴성을 확보하였다.
• 온도 의존적 보즈-아인슈타인 점유 인자와 온도에 따른 격자 열팽창을 기하학적 최적화를 통해 고려하였다.
• 수렴성이 잘 확보된 16×16×16 k-점 메쉬를 사용하여 0 K 및 300 K에서 계산를 수행하여 온도 의존성을 평가하였다.

실험 결과

연구 질문

• RQ1 이론적 보고와 실험적 보고가 상반된 바를 감안할 때, MgB2의 E2g 격자파 모드에서 Γ 점에서의 진짜 비조화 주파수 이완의 크기는 얼마인가?
• RQ2 서로 반대 부호를 가진 삼중 및 사중 격자파 산란 기여도가 어떻게 합쳐져 네트워크 비조화 이완을 결정하는가?
• RQ3 왜 Raman 스펙트로스코피 측정은 큰 비조화 이완(~18%)을 암시하는 반면, 비탄성 X선 산란 측정은 주로 조화 모델과 일치하는가?
• RQ4 브릴루아 존 전역에서 운동량 및 모드 의존성 산란을 포함함으로써 고정된 격자파 근사 모델에 비해 예측된 비조화 이완이 얼마나 변화하는가?
• RQ5 일부 모델이 예측한 바와 같이, E2g 모드의 비조화 이완은 Γ–A 경로를 따라 균일한가?

주요 결과

• Γ 점에서 E2g 모드의 비조화 주파수 이완은 +3.50 meV이며, 이는 300 K에서의 조화 주파수(64.51 meV)보다 5.4% 증가한 것이다.
• 이 이완은 사중 격자파 산란으로부터 기인한 +10.68 meV 기여와 삼중 격자파 산란으로부터 기인한 −3.03 meV 기여 간의 상쇄로 발생한다.
• A 점에서의 총 이완은 4.66 meV로 Γ 점의 이완과 거의 동일하여, Γ–A 경로를 따라 균일한 비조화 이완이 일어남을 시사한다.
• M 점에서의 이완은 더 작으며, E2g 모드의 비조화성이 Γ–A 경로를 따라 가장 두드러짐을 확인한다.
• E2g 주파수 이완의 온도 의존성은 극히 미미하며, 0 K에서의 68.26 meV에서 300 K에서의 68.01 meV로 약간 감소한다.
• 결과는 비탄성 X선 산란 측정과 뛰어난 일치를 보이며, 실험 오차 범위 내에서 E2g 분산이 조화 모델과 일치함을 확인한다.