분자 동역학으로 본 <100> 인장 하중 탄탈륨 단결정의 탄성 거동: 온도, 압력, 상전이 통찰
초록
원자론적 시뮬레이션은 탄탈륨(Ta)의 기계적 응답을 지배하는 메커니즘 규명에 기여합니다. 본 연구에서는 분자 동역학(MD)을 이용해 온도와 압력이 <100> 인장 하중 하 탄탈륨 단결정 탄성 특성에 미치는 영향을 조사했습니다. 먼저 Zhou-EAM과 Ravelo-EAM 두 EAM 퍼텐셜을 넓은 압력 범위에서 탄탈륨 탄성 상수 재현 능력 비교. Ravelo-EAM 퍼텐셜이 140 GPa까지 탄성 계수 정확하게 포착해 우수. 이 퍼텐셜 기반 MD 시뮬레이션에서 1 K에서 체심입방(BCC)→면심입방(FCC) 상전이가 파단에 앞서 발생, 모델 크기나 변형률 속도 무관. 1-1500 K에서 <100> 탄성 계수 E₁₀₀은 온도에 따라 선형 감소, 항복 응력은 이차 함수 의존. 압력 의존 시뮬레이션(0-140 GPa)에서 E₁₀₀ 단조 증가. 이 결과들은 극한 조건 하 탄탈륨 탄성 거동 예측을 위한 포괄적 원자론적 프레임워크 제공.
배경
탄탈륨은 상온에서 체심입방(BCC) 격자 결정화, 뛰어난 상 안정성, 높은 융점(3269 K), 내산화성으로 유명. 이 특성들로 미세전자 확산 방지막, 내마모 코팅, 고온 초합금 후보로 이상적. 다이아몬드 앤빌 셀, 밀도 범함수 이론(DFT), MD 등 실험·이론 연구로 고압·고온 하 Ta 거동 탐구, BCC 상이 135 GPa까지 지속되고 탄성 계수가 압력에 따라 증가함 일관되게 확인. 이러한 발전에도 불구하고 Ta의 인장 하중 동적 응답에 대한 체계적 원자론적 연구는 부족. 본 연구는 Ta의 <100> 인장 시 크기, 변형률 속도, 온도, 압력 효과를 탄성 특성 평가와 인장 유도 상전이 가능성 탐구로 다룹니다.
방법 및 실험 세부사항
물리 모델링
BCC 단위 격자를 <001>, <010>, <100> 방향으로 복제해 입방 Ta 셀 구성, 격자 상수 a=b=c=3.301 Å. 네 가지 셀 크기(12a=3.96 nm, 18a=5.94 nm, 24a=7.92 nm, 30a=9.90 nm), 3456~54000 원자 포함. 그림 1은 주요 구조로 쓰인 3.96 nm 모델.
MD 시뮬레이션 프로토콜
모든 시뮬레이션 LAMMPS로 수행. Zhou-EAM과 Ravelo-EAM 두 EAM 퍼텐셜을 0 K와 정수압 하 탄성 상수 재현 능력 평가. Ravelo-EAM 퍼텐셜(Ravelo et al., Phys. Rev. B 2013, 88, 134101)이 관심 압력 범위에서 Zhou-EAM 능가. 직접법으로 C₁₁, C₁₂ 계산 후 식 (1)로 <100> 탄성 계수 유도. 후속 시뮬레이션: 1 fs 타임스텝, 전 방향 주기적 경계 조건, 50 ps NVT 완화 후 NPT 평형으로 정수압(0-140 GPa) 부과. x 방향 인장 하중을 5×10⁸ s⁻¹~7.5×10⁹ s⁻¹ 변형률 속도로 15% 신장까지 적용.
결과 및 토의
인장 과정
그림 3은 영압에서 <100> 인장 시 응력-변형률 곡선과 해당 원자 구성. 탄성 변형 ~7% 변형률까지 선형, 이후 ~7.4%에서 BCC→FCC 변환 개시, ~9.8% 완료. FCC 상 지속되다 ~13% 근처 표면 파단. 변환 및 파단 메커니즘은 모든 모델 크기와 변형률 속도에서 일관, 탄성 계수와 항복 응력에 크기/속도 영향 미미.
크기 및 변형률 속도 효과
표 2: 모델 크기(3.96-9.90 nm)가 탄성 계수(~139 GPa)나 항복 응력 영향 없음. 표 3: 변형률 속도 5×10⁸~7.5×10⁹ s⁻¹ 사이에서도 두 특성 변화 없음. 따라서 3.96 nm 셀과 5×10⁸ s⁻¹ 변형률 속도 후속 분석 채택.
온도 의존성
그림 4: 1 K~1500 K 응력-변형률 곡선 변화. 탄성 계수 E₁₀₀과 항복 응력 모두 온도 증가에 따라 단조 감소. 계수 ~42% 감소(136.5→76.7 GPa), 항복 응력 ~51% 감소(~8→4 GPa). 선형 피팅: E₁₀₀ = 138.07 – 0.04094 T (GPa), 1500 K까지 유효. 항복 응력은 이차 함수 피팅.
압력 의존성
그림 6: 여러 온도에서 정수압(0-140 GPa) 대비 E₁₀₀. 1 K에서 동적 결과가 식 (1) 기반 정적 계산과 잘 일치. ≤600 K에서 20-60 GPa 사이 오목 영역 출현, 선행 MD 연구와 일치. 고온 곡선은 단조. 모든 온도에서 ~40 GPa 이상 압력에서 E₁₀₀ 증가, 탄탈륨 강성이 압축 하 강화 확인.
결론
광범위한 온도(1-1500 K)와 압력(0-140 GPa)에서 <100> 인장 하 탄탈륨 단결정 MD 종합 조사 수행. Ravelo-EAM 퍼텐셜이 고압 시뮬레이션에 신뢰성 입증. 주요 발견: (1) 저온에서 파단 전 BCC→FCC 상전이 선행; (2) 탄성 계수와 항복 응력 각각 선형, 이차 함수로 온도 증가 시 감소; (3) ~40 GPa 이상 압력에서 탄성 계수 증가; (4) 모델 크기와 변형률 속도는 측정된 탄성 특성에 미미한 영향.