[논문 리뷰] Two-temperature equations of state for d-band metals irradiated by femtosecond laser pulses

K. P. Migdal, N. A. Inogamov외 2명|arXiv (Cornell University)|2017. 02. 02.

Laser Material Processing Techniques인용 수 3

한 줄 요약

이 연구는 밀도함수이론(DFT)과 볼츠만 운동이론을 사용하여 펌터세컨드 레이저 조사 하에서 d-대역 금속(Cu, Fe, Ta)에 대한 이온도 온도 방정정식을 개발한다. 비평형 상태에서 전자 온도가 최대 55,000 K일 때, 격자 불안정성, 전자 열전도도 및 저항률을 예측하며, 특히 Cu와 Ni에서 전자-이온 결합 및 전도도에 강한 전자 온도 의존성이 나타나며, 고온에서 QMD 시뮬레이션과 양호한 일치를 보인다.

ABSTRACT

The cold curves for energy and pressure of Copper, Iron, and Tantalum were obtained using methods of the density functional theory. We consider hydrostatic and uniaxial deformations in the range from double compression of the initial volume per atom to double stretching. The presence of allotropic transformation from $α$ - phase of Iron to the hexaferrum with the growth of pressure is observed. In the case of hydrostatic deformations we also have obtained analogous cold curves, but with non-zero electronic temperatures in the range up to 5 eV. The similar volume and electronic temperature ranges have been considered recently. The behavior of electronic internal energy, pressure, and density of states was investigated in the volume and temperature ranges called above. The maximum hydrostatic strains and the types of lattice instabilities were theoretically predicted for the considered metals. The influence of high electronic temperature on the electronic heat conductivity and electric resistivity has been provided for d-band metals by the approach based on the solution of Boltzmann kinetic equation in $τ$-approximation. This data is compared with the results of quantum molecular dynamics for Gold.

연구 동기 및 목표

• 펌터세컨드 레이저 조사 하에서 비평형 전자 및 이온 부분계를 고려하여 d-대역 금속(Cu, Fe, Ta)에 대한 이온도 온도 방정정식을 개발한다.
• 고온 전자 온도(최대 55,000 K) 및 정적/이축 응력이 전자 내부 에너지, 압력 및 격자 안정성에 미치는 영향을 조사한다.
• 이완 시간(τ) 근사법을 사용한 볼츠만 운동방정식을 통해 전자 열전도도(κ₂ₜ) 및 저항률(σ)을 계산하며, 온도 의존성 전자 상태 밀도를 포함한다.
• DFT로 계산된 탄성 모odulus를 이용해 극한 조건에서 격자 불안정성(부피, 비틀림, 테트라고널) 및 동소변환(예: Fe의 bcc에서 hcp로의 전이)을 예측한다.
• 실험 데이터(예: Fe의 램프 웨이브 압축) 및 금속 Au에 대한 양자 분자역학(QMD) 시뮬레이션과의 비교를 통해 결과를 검증한다.

제안 방법

• Cu(fcc), Fe(bcc, hcp), Ta(bcc)에 대해 VASP와 PAW 허위전자함수(LDA, PBE, PBEpv)를 사용하여 영온도 상태(격자 온도 0 K)의 냉각 곡선을 계산하였으며, 이는 이중 압축에서 이중 팽창까지의 부피 범위를 포함한다.
• 페르미-디랙 스메어링과 DFT로부터 유도된 전자 상태 밀도(DOS)를 사용하여 최대 55,000 K까지의 비제로 전자 온도를 포함하여 이온도(2T) 열역학 계산을 수행하였다.
• 이완 시간(τ) 근사법을 사용한 볼츠만 운동방정식을 통해 전자 열전도도(κ₂ₜ) 및 저항률(σ)을 계산하였으며, 전자-격자 상호작용은 수정된 디바이 방법으로 모델링하였다.
• 탄성 모듈러스 조합(c₁₁+2c₁₂, c₁₁−c₁₂, c₄₄)을 계산하여 격자 불안정성을 식별하였으며, 음수 값은 각각 부피, 비틀림, 테트라고널 불안정성을 나타낸다.
• Fe의 램프 웨이브 압축 실험 데이터(560 GPa) 및 Au에 대한 QMD 시뮬레이션과의 비교를 통해 결과를 검증하였으며, 특히 전도도 및 부피에 따른 압력 예측에서 양호한 일치를 보였다.
• 전자 전도도에서 d-전자 기여를 고려하기 위해 s-전자와 d-전자 기여를 비교하였으며, Cu와 Au에서 s-전자 기여가 지배적임을 확인하였다.

실험 결과

연구 질문

• RQ1 고온 전자 온도(최대 55,000 K)는 Cu, Fe, Ta와 같은 d-대역 금속의 상태방정식, 압력 및 내부 에너지에 어떻게 영향을 미치는가?
• RQ2 극한의 정적 및 이축 응력 하에서 어떤 종류의 격자 불안정성(부피, 비틀림, 테트라고널)이 발생하며, 어떤 부피/온도 임계값에서 나타나는가?
• RQ3 전자 온도 의존성 전자 상태 밀도가 전자-이온 결합 및 전자 열전도도에 미치는 영향은 무엇인가?
• RQ4 s-전자와 d-전자 기여는 전자 전도도에서 얼마나 다를까? 비평형 조건에서 저항률과 전도도에 어떤 영향을 미치는가?
• RQ5 DFT 기반 이온도 온도 상태방정식 및 전도도 특성은 실험적 램프 웨이브 압축 데이터 및 Au와 Fe에 대한 QMD 시뮬레이션과 얼마나 잘 일치하는가?

주요 결과

• PBEpv 허위전자함수는 Cu와 Fe의 압력-부피 실험 데이터와 가장 양호한 일치를 보였으며, LDA 및 PBE보다 부압축 영역에서 뛰어난 성능을 보였다.
• Fe의 경우, 전자 온도가 최대 55,000 K일 때도 bcc에서 hcp로의 동소변환이 명확히 관측되지 않았으며, 이는 고온 전자 에너지가 격자 기반 상전이를 억제함을 시사한다.
• Cu의 전자-이온 결합(α)은 전자 온도 상승에 따라 뚜렷한 감소를 보였으며, 이는 온도 의존성 DOS의 영향 때문이다. κ₂ₜ의 피크(약 5,000 K)는 d-밴드 에지 이하에서 약화된 d-밴드 영향(ε₂,Au = −1.7 eV)에 기인한다.
• 모든 전자 온도 범위에서 d-전자 전도도 대비 s-전자 전도도 비율(κ₂ₜ,ₔ/κ₂ₜ,ₛ)은 0.02 이하로 유지되어 s-전자 기여가 전자 전도도에서 지배적임을 확인한다.
• Au의 경우, 전자 온도(Tₑ)가 300 K 및 1,200 K일 때 저항률(σ)은 거의 변화하지 않았지만, Ni는 융해 온도 근처에서 σ가 급격히 감소하여 전이 금속에서 전자 온도 의존성이 뚜렷함을 보였다.
• 고온 전자 온도(>30,000 K)에서 DFT 기반 κ₂ₜ 결과는 QMD 시뮬레이션과 양호한 일치를 보였으며, 이는 이온도 모델의 정확성이 극한 조건에서 검증됨을 의미한다.