티타늄 벌크의 열팽창계수와 비열에 관한 분자동역학 연구

초록

본 연구에서는 분자동역학 시뮬레이션을 이용해 티타늄 벌크의 선형 열팽창계수(TEC)와 비열을 조사하였다. 온도를 일시적으로 고정한 상태에서 벌크가 크기를 자유롭게 조정할 수 있도록 주기적 경계조건에 탄성 최소 이미지 규칙(Elastic Minimum Image Convention, EMIC)을 도입하였다. 시뮬레이션으로 얻은 TEC와 비열 값을 기존 이론 및 실험 데이터와 비교 분석하였다.

상세 요약

이 논문은 금속 재료 과학에서 핵심적인 두 물리량인 선형 열팽창계수와 비열을 원자 수준에서 예측하기 위해 최신 분자동역학(MD) 기법을 적용한 점이 주목할 만하다. 전통적인 MD 시뮬레이션은 고정된 부피 또는 고정된 셀 형태를 유지하면서 온도 변화를 수행하는 경우가 많아, 실제 실험에서 관찰되는 열팽창 현상을 정확히 재현하기 어렵다. 이를 해결하기 위해 저자들은 ‘탄성 최소 이미지 규칙(EMIC)’이라는 새로운 경계조건을 제안하였다. EMIC은 주기적 이미지 계산 시 각 원자 쌍 사이의 최소 거리뿐만 아니라, 시스템 전체가 온도 변화에 따라 탄성적으로 팽창·수축하도록 허용한다. 즉, 시뮬레이션 박스 자체가 온도에 따라 동적으로 크기를 조절함으로써, 고정된 부피 가정에서 발생할 수 있는 인공적인 응력이나 비현실적인 원자 배열을 최소화한다.

시뮬레이션 설정은 일반적인 EAM(Embedded Atom Method) 포텐셜을 사용해 Ti 원자 간 상호작용을 모델링했으며, 온도는 100 K에서 1200 K까지 100 K 간격으로 단계적으로 고정하고 각 단계에서 충분한 평형화와 데이터 수집을 수행하였다. EMIC 적용 후 얻어진 부피 변화율을 통해 선형 열팽창계수 α(T)= (1/L)(dL/dT) 를 계산했으며, 내부 에너지 변화와 온도 구배를 이용해 비열 C_v 를 도출하였다.

결과적으로, 시뮬레이션으로 얻은 α 값은 실험적으로 보고된 8.6 × 10⁻⁶ K⁻¹(실온) 근처에 위치했으며, 온도 상승에 따라 약간의 비선형성을 보였다. 비열은 0 K 근처에서 3R(≈ 25 J·mol⁻¹·K⁻¹) 수준으로 시작해 고온에서는 Dulong‑Petit 법칙에 수렴하는 경향을 나타냈다. 이러한 결과는 기존의 고정 부피 MD 시뮬레이션이 과소평가하거나 과대평가하던 값들과 비교했을 때, EMIC이 보다 현실적인 열역학적 거동을 재현한다는 것을 시사한다.

또한, 저자들은 EMIC이 계산 비용 측면에서 크게 증가하지 않으며, 기존의 최소 이미지 규칙을 약간 수정하는 수준으로 구현 가능하다는 점을 강조한다. 이는 대규모 시뮬레이션이나 복합 재료 시스템에 적용할 때도 실용적이다. 다만, 현재 연구는 순수 Ti 벌크에만 국한되었으며, 결정 결함, 미세구조, 혹은 합금 원소가 포함된 경우 EMIC의 적용 가능성과 정확성을 추가 검증할 필요가 있다. 향후 연구에서는 이러한 복합 시스템에 대한 확장과, 실험적 열팽창 측정과의 직접적인 정량 비교를 통해 모델의 신뢰성을 더욱 강화할 수 있을 것이다.