수소가 주도하는 Zr(0001) 면 전위 결함 에너지 변화: 첫 원리 계산 연구

초록

본 연구는 밀도범함수이론(DFT)을 이용해 수소가 hcp Zr의 basal 면 {0001}에서 <10‑10> 및 <11‑20> 전위 전이(Generalized Stacking Fault, GSF) 에너지에 미치는 영향을 조사하였다. 수소 원자가 전위면에 가까울수록 <10‑10> 방향의 GSF 에너지가 크게 감소하고, 수소 농도가 증가하면 두 전위 방향 모두에서 에너지 장벽이 추가로 낮아진다. Bader 전하 분석을 통해 Zr 원자와 수소 사이의 전하 이동이 Zr 원자 간의 정전기적 반발을 강화시켜 전위 에너지 감소를 초래한다는 메커니즘을 제시한다.

상세 분석

본 논문은 hcp Zr의 basal 면 {0001}에서 두 가지 전위 전이 경로, 즉 1/3 <10‑10>와 1/3 <11‑20>에 대한 GSF 에너지 곡선을 DFT(PBE‑GGA)와 PAW 방법으로 정밀하게 계산하였다. 3 × 3 × 2 초셀(36 Zr 원자) 내에 수소를 텐타헤드(T) 인터스티셜 자리에 배치하고, H 원자를 슬립면에서 1층, 3층, 5층에 각각 위치시켜 H1, H3, H5 모델을 만든 뒤, H5에 추가로 두 번째 H 원자를 삽입한 H57 모델까지 고려하였다.

1. 전위 에너지 변화

   • 순수 Zr의 <10‑10> GSF 에너지 최대값(γ_US)은 0.276 J m⁻²이며, 전위 위치는 0.7–0.8 배 변위에서 발생한다.

   • H 원자가 슬립면에서 멀리(H1, H3) 있을 경우 γ_US는 약 5 % 정도 감소하지만, 슬립면에 가장 가까운 H5에서는 γ_US가 0.152 J m⁻²(≈45 % 감소)까지 급격히 낮아진다. γ_SF(전위 고정 에너지) 역시 0.039 J m⁻²(≈83 % 감소)로 크게 감소한다.

   • H57(수소 농도 증가)에서는 γ_US가 0.070 J m⁻², γ_SF가 -0.115 J m⁻²까지 내려가, 전위가 실제로 안정적인 결함 상태가 됨을 보여준다.

   • <11‑20> 방향에서는 전반적인 감소폭이 <10‑10>보다 작다. H5에서 γ_US는 0.407 J m⁻²(≈9 % 감소) 수준이며, H57에서도 0.285 J m⁻²로 약간 낮아진다. 이는 두 전위면 사이의 원자 간 거리 차이(≈0.15 Å)로 인한 정전기 반발 감소가 원인으로 해석된다.

2. 구조적 요인

   • H 원자는 T‑site에 안정적으로 자리잡으며, H 삽입 시 주변 Zr 원자와의 거리(특히 슬립면 인접 Zr층) 가축소된다. 그러나 단순히 격자 팽창만으로는 GSF 감소를 설명하기 어려워, 전하 재분배가 핵심 메커니즘임을 제시한다.

3. Bader 전하 분석

   • H 원자는 약 0.79 e의 전자를 받아 전기음성도가 높아지며, 인접 Zr 원자는 각각 ~0.2 e(단일 H)에서 최대 0.6 e(두 개의 H, H57)까지 전자를 잃는다.
   • 전하가 동일한 부호를 가진 Zr 원자들 사이에 축적되면 정전기적 반발력이 증가하고, 이는 전위면을 이동시키는 에너지 장벽을 낮추는 효과로 작용한다.
   • C와 N을 H 대신 삽입한 경우 전하 이동량이 더 커서 GSF 에너지 감소폭도 더 크게 나타났으며, 이는 전하 이동량과 GSF 감소 사이의 정량적 연관성을 뒷받침한다.

4. 전위 방향 의존성

   • <10‑10> 방향은 원자 배열이 촘촘해 전하 이동에 따른 정전기 반발이 크게 작용한다. 반면 <11‑20> 방향은 원자 간 간격이 더 넓어 동일한 전하 이동이 미치는 효과가 약해, 에너지 감소가 제한적이다.

5. 의의 및 적용

   • 수소가 Zr 내부에 고용될 경우, 특히 전위면 근처에 축적될 경우 전위 이동이 용이해져 전위 기반 연성(플라스틱 변형)이 촉진될 수 있다. 이는 경수로 연료봉 피복재에서 수소 침투에 따른 기계적 거동 변화를 이해하는 데 중요한 통찰을 제공한다.