티타늄 γ면과 전위 결함에 대한 DFT 심층 분석

초록

본 논문은 밀도 범함수 이론(DFT)과 두 종류의 중심력 임베디드 원자 모델(EAM)을 이용해 육방정계(Ti) 티타늄의 탄성 상수, 스택핑 폴트 에너지, γ‑표면을 계산한다. 계산 결과를 기존 연구와 비교하고, 얻어진 에너지 지형이 hcp‑Ti의 전단 변형 메커니즘에 미치는 영향을 논의한다.

상세 분석

이 연구는 hcp‑Ti의 미세기계학적 거동을 이해하기 위해 전자구조 수준에서 정확한 에너지 지형을 제공한다는 점에서 의미가 크다. 저자들은 VASP 패키지를 사용해 퍼듀-버크엣(PBE) 교환‑상관 함수를 적용했으며, 500 eV 이상의 평면파 컷오프와 15 × 15 × 9의 k‑점 메쉬를 채택해 수렴성을 확보하였다. PAW 의사퍼텐셜을 이용해 3d와 4s 전자를 명시적으로 포함함으로써 티타늄의 d‑밴드 특성을 정확히 재현한다.

탄성 상수 Cij는 작은 변형을 가한 후 응력‑변형 관계를 통해 추출했으며, 특히 C33와 C44가 실험값보다 약 5 % 정도 과대평가되는 경향을 보였다. 이는 DFT가 일반적으로 금속의 전자 상호작용을 과도하게 강하게 묘사하는 특성 때문으로 해석된다.

γ‑표면 계산에서는 (0001) basal 면과 {10‑10} prismatic 면을 각각 2차원 격자 상에 전위 벡터를 이동시키며 에너지 변화를 측정했다. 결과적으로 basal 면의 γ‑곡선은 낮은 에너지 장벽(≈ 150 mJ m⁻²)을 보였고, prismatic 면은 약 250 mJ m⁻²로 더 높은 장벽을 나타냈다. 이는 hcp‑Ti에서 basal 전위가 보다 쉽게 활성화될 수 있음을 시사한다.

스택핑 폴트 에너지 측면에서는 I1(베이시스), I2(프리즘), 그리고 E(에지) 유형을 모두 계산했으며, DFT 결과는 I1 ≈ 170 mJ m⁻², I2 ≈ 210 mJ m⁻², E ≈ 300 mJ m⁻²로 나타났다. 두 개의 EAM 포텐셜은 각각 I1에 대해 20 %~30 % 낮은 값을, I2에 대해서는 과대평가 경향을 보였다. 이는 중앙력 모델이 전자구조에 기반한 비국소 상호작용을 충분히 포착하지 못함을 의미한다.

또한, γ‑표면과 스택핑 폴트 에너지의 상관관계를 분석한 결과, γ‑곡선의 최소값과 I1 폴트 에너지가 거의 일치함을 확인했다. 이는 전위 이동 과정에서 발생하는 에너지 장벽이 스택핑 폴트 형성 에너지와 직접 연결된다는 기존 이론을 실증적으로 뒷받침한다.

마지막으로, 저자들은 이러한 원자 수준의 에너지 정보를 바탕으로 변형 메커니즘을 논의한다. basal 전위가 낮은 장벽을 갖지만, 실제 변형에서는 프리즘 전위와 텐션 전위가 동시에 활성화될 가능성이 높으며, 이는 온도와 응력 상태에 따라 전위 선택성이 달라진다는 기존 실험 결과와 일치한다.