분자동역학 시뮬레이션을 통한 초정밀 다이아몬드 절삭에서의 세륨 변형 메커니즘 규명

초록

구조적 상전이와 전위(dislocation) 역학 간의 결합은 세륨의 나노스케일 변형을 이해하는 데 복잡성을 더합니다. 본 연구에서는 분자동역학(MD) 시뮬레이션을 활용하여 초정밀 다이아몬드 절삭 공정 중 세륨의 변형 메커니즘을 규명했습니다. 경험적으로 도출된 Ce–Ce 및 Ce–C 포텐셜을 포함하는 MD 모델을 구축하고, 두 가지 면심입방(fcc) 세륨 상의 탄성 거동을 평가했습니다. 시뮬레이션 결과, 전위 미끄러짐이 소성 변형을 지배하는 것으로 나타났으며, 방사분포함수(RDF) 분석에서는 가공 표면과 칩에서만 미미한 γ→δ 상전이만이 관찰되었습니다. 매개변수 연구를 통해 고품질 가공 표면을 얻기 위한 최적의 가공 조건을 도출했습니다.

배경

세륨(Ce, 원자번호 58)은 독특한 기계적, 물리적, 화학적 특성으로 주목받는 란타넘족 금속입니다. 표면 형상은 세륨 부품의 성능, 내식성, 수소 저장 거동에 결정적인 영향을 미칩니다[1-3]. 초정밀 다이아몬드 절삭은 최소한의 표면하 손상으로 초평활 표면을 생성할 수 있습니다[6,7], 그러나 세륨 절삭에 관한 실험적·이론적 연구는 거의 없습니다. 초정밀 절삭에서 공구 끝 반경이 절삭 깊이와 비슷하기 때문에 재료의 물성이 공정을 지배하며, 이는 세륨 절삭 메커니즘 규명을 더욱 어렵게 만듭니다.

세륨은 4f 전자의 비국재화에 의해 구동되는 풍부한 압력-온도 상평형도를 가집니다. 상압 및 110 K 미만에서 α-Ce(fcc) 상이 안정적이며, 45-275 K에서 β-Ce(dhcp), 270-999 K에서 γ-Ce(fcc), 999-1071 K에서 δ-Ce(bcc) 상으로 전이합니다[8-11]. 특히 295 K, 8 kbar에서의 γ→α 전이는 약 20%의 부피 감소를 수반합니다[12-14], 반면 γ→δ 전이는 고온 고압 조건에서 유도될 수 있습니다. 이러한 상전이는 결합 특성과 탄성 물성을 변화시켜 변형 거동에 영향을 미칩니다. 금속 세륨은 전위 활동에 의해 지배되는 연성 재료이기도 합니다[15]. 그러나 다이아몬드 절삭 중 상전이와 전위 간의 상호작용은 여전히 미개척 영역으로 남아 있습니다.

MD 시뮬레이션은 원자 스케일에서 기본적인 가공 메커니즘을 탐구하는 강력한 수단을 제공합니다. 선행 연구들은 Cu, Al, Fe, 고엔트로피 합금 등의 나노스크래칭과 직교 절삭을 조사하여 전위 핵생성, 표면 파일업, 마찰 거동 등에 대한 통찰을 제공했습니다[16-27]. 이러한 연구들은 결정 방위, 공구 경사각, 상호작용 포텐셜의 중요성을 강조합니다. 그러나 세륨의 기계적 가공에 대한 MD 연구는 전무한 실정으로, 본 연구의 동기가 되었습니다.

연구 방법

다이아몬드 절삭 MD 모델

그림 1은 시뮬레이션 설정을 보여줍니다: 약 1×10⁶ 개의 γ-Ce 원자를 포함하는 단결정 세륨 공작물(41×25×31 nm)과 9° 여유각을 가진 다이아몬드 절삭 공구(0.1×10⁶ 개의 C 원자). 공작물은 고정된 2 nm 하부층과 뉴턴 역학을 따르는 이동층(1 fs 시간 간격)으로 나뉩니다. 종방향으로만 주기적 경계 조건이 적용됩니다. (010), (110), (111) 세 가지 결정 방위를 고려했으며, 일곱 가지 경사각(-30°, -20°, -10°, 0°, 10°, 20°, 30°)을 테스트했습니다. 다이아몬드 공구는 강체로 처리되어 C–C 상호작용은 생략했습니다.

원자 간 상호작용은 다음과 같이 모델링했습니다: Ce–Ce 힘은 Sheng 등의 매립 원자법(EAM) 매개변수화[29]를 사용하여 fcc 세륨의 탄성 상수를 정확히 재현했습니다. Ce–C 상호작용은 모스 포텐셜을 채택했습니다:

$$E_{\text{tot}}=\sum_{ij}D_0\left[e^{-2\alpha(r-r_0)}-2e^{-\alpha(r-r_0)}\right]$$

여기서 D₀=0.087 eV, α=5.14, r₀=2.93 Å이며, 1.0 nm 컷오프를 적용했습니다[30].

시스템은 먼저 30 K, 0 bar에서 NPT 앙상블로 평형화한 후, NVT 앙상블에서 일정 속도 100 m/s, 절삭 깊이 4 nm로 다이아몬드 절삭을 수행했습니다. 절삭 힘은 절삭 방향(그림에서 빨간 화살표)을 따라 분해했습니다. 시뮬레이션 차원과 절삭 속도는 실험적 초정밀 조건보다 몇 자릿수 작지만, 본 MD 모델은 공정의 본질적인 물리학을 포착합니다. 공통 이웃 분석(CNA)으로 fcc, hcp, bcc, 결함 원자를 구분했습니다[31]. 시뮬레이션은 LAMMPS(1 fs 단계)로 수행하고 OVITO로 시각화했습니다[32,33].

세륨 상 특성화

문헌[8-11]의 격자 상수를 사용하여 γ, α, β, ε, δ 다섯 가지 세륨 상을 모델링했습니다. 벌크 구성에 단축 인장, 전단, 압축을 가해 탄성 상수와 기계적 물성을 추출했습니다(표 2). 계산된 γ-Ce의 영률(24.17 GPa)은 실험적 나노인덴테이션 데이터(36.7 GPa)와 일치하며[10], 탄성 상수는 비탄성 중성자 산란 결과와도 일치합니다[34]. 더 조밀한 α-Ce는 γ-Ce보다 현저히 높은 강성을 보입니다.

방사분포함수(RDF) 분석으로 상 간의 구조적 차이를 특성화했습니다. 그림 2a는 각 상에 고유한 RDF 피크를 보여주며, 최근접 이웃 거리(예: γ-Ce: 3.64, 5.13, 6.30 Å; α-Ce: 3.41, 4.85, 5.92 Å)에 대응합니다. γ-Ce를 균일하게 압축하여 20% 부피 감소를 유도하면 α-Ce와 일치하는 RDF를 얻습니다(그림 2b), 이는 γ→α 전이를 확인합니다[12-14].

결과 및 논의

세륨의 가공 메커니즘

0° 경사각에서 Ce(010) 다이아몬드 절삭 MD 시뮬레이션은 세 가지 힘 영역을 드러냅니다(그림 3). 초기에는 접착 힘으로 음의 절삭력과 법선력이 발생하고, 이어 탄성 변형으로 힘이 급격히 상승합니다. 약 2.3 nm 절삭 길이에서 급격한 힘 감소는 소성 변형 개시를 나타내며, 이는 자유 표면 근처에서 핵생성된 1/6<112> 쇼클리 부분 전위가 {111} 면을 따라 미끄러지며 지배적입니다. 10-35 nm 구간에서 절삭력과 법선력은 안정화되며, 공구 앞과 아래에 전위가 축적됩니다. 35 nm 이후 공구가 공작물에서 분리되면서 힘이 감소합니다. 전체 과정에서 법선력은 절삭력보다 낮게 유지됩니다.

그림 4는 주요 절삭 길이에 따른 결함 구조 진화와 표면 형상을 보여줍니다. 쇼클리 부분 전위는 전위 핵에 의해 경계 지어진 적층 결함을 생성합니다. 절삭이 진행됨에 따라 전면에서 밀집된 부분 전위 구름이 경사면을 따라 칩을 형성합니다(그림 4f). 공구 뒤의 전위는 상부 표면으로 상승하여 표면 회복을 촉진합니다. 공구가 왼쪽 경계에 접근하면 전위 전파가 방해되어 칩 부피가 증가합니다(그림 4g). 칩 분리 후 자유 표면에서의 소멸로 전위 밀도가 감소합니다(그림 4h).

결함 영역의 RDF 분석은 제한적인 γ→δ 상전이를 확인합니다. 표면 아래 결함 영역의 RDF는 δ-Ce 피크와 일치하며, 국부적 가열에 의한 미미한 γ→δ 전이를 시사합니다(그림 5a). 칩에서도 유사한 δ-Ce 신호가 나타나지만(그림 5b), 전체 양은 미미하여 절삭 중 상전이가 지배적이지 않음을 나타냅니다.

경사각의 영향

그림 6은 절삭력과 법선력이 경사각 증가에 따라 감소함을 보여주며, 이는 머천트(Merchant)의 전단면 에너지 이론과 일치합니다[37]. 높은 경사각은 힘 차이를 줄이고 가공 안정성을 향상시킵니다.

결함 분석(그림 7a,b)에서 음의 경사각(-30°)은 양의 경사각(30°)보다 높은 전위 밀도를 보입니다. 전위 세그먼트 수(그림 7c)는 경사각에 따라 단조롭게 감소하여, 덜 복잡한 소성 변형과 낮은 가공 경화를 반영합니다.

표면 형상(그림 8)은 이를 뒷받침합니다: 음의 경사각은 비대칭적이고 현저한 표면 파일업을 생성하는 반면, 30° 경사각은 가장 낮은 파일업 부피와 대칭적 프로파일을 생성합니다. 따라서 30° 경사각은 낮은 절삭력, 최소 전위 밀도, 우수한 표면 조도를 위한 최적 조건입니다.

결정 방위의 영향

최적 30° 경사각에서 Ce(010), Ce(110), Ce(111) 방위를 비교했습니다(그림 9). Ce(010)은 가장 낮은 절삭력과 전위 수를 보입니다. Ce(111)은 전위 수가 가장 적지만, {111} 면이 표면과 평행하여 표면 파일업을 유발하고 절삭 저항을 증가시킵니다. 가공 표면 형상(그림 10)에서 Ce(010)은 가장 작은 파일업 부피와 가장 대칭적인 프로파일을, Ce(110)은 가장 큰 파일업을 보입니다.

결론

EAM과 모스 포텐셜을 포함하는 세륨의 초정밀 다이아몬드 절삭 MD 프레임워크를 개발했습니다. 세륨 상의 탄성 상수, 기계적 물성, 상전이 거동을 검증하여 모델의 예측 능력을 확인했습니다. 시뮬레이션 결과 전위 핵생성과 미끄러짐이 소성 변형을 지배하며, γ→δ 상전이는 미미하게만 관찰됩니다. 최적 가공 조건인 30° 경사각과 (010) 결정 방위는 최저 절삭력, 최소 전위 밀도, 최고 표면 품질을 제공합니다.