미스핏 각도에 따른 경계 윤활 특성 연구

초록

본 연구는 좌표·속도 의존 감쇠를 포함한 랭게빈 분자동역학을 이용해, 기판-윤활제 상호작용이 강하고 윤활제 내부 결합이 약한 ‘소프트’ 윤활제 층이 두 고체 사이에 존재할 때의 마찰 거동을 조사한다. 저속 구동에서는 스틱-슬립이 나타나며, 슬라이딩 중 윤활제가 완전히 녹을 수도, 층간 미끄럼(LoLS) 형태를 유지할 수도 있다. 특히, 두 고체 표면 사이의 미스핏 각도가 정지 마찰 한계와 동적 마찰에 미치는 영향을 체계적으로 분석하고, 다결정 접촉을 가정한 정지 마찰 임계값 분포를 제시한다.

상세 분석

이 논문은 전통적인 마찰 시뮬레이션에서 흔히 가정되는 일정한 감쇠 계수를 넘어, 입자 위치와 속도에 따라 변하는 비선형 감쇠를 적용한 랭게빈 방정식 기반 분자동역학(MD) 모델을 구축하였다. ‘소프트’ 윤활제라 함은 윤활제-기판 간 포텐셜 깊이가 윤활제-윤활제 간 포텐셜보다 크게 설정된 경우를 의미한다. 이러한 설정은 실제 실험에서 관찰되는, 기판에 강하게 흡착되면서도 내부는 비교적 유동적인 유기막층을 모사한다. 시뮬레이션은 두 평면 기판 사이에 2~3층 두께의 윤활제 원자를 배치하고, 한쪽 기판을 일정 속도로 이동시키는 방식으로 진행되었다. 저속 구동(≈0.01 σ/τ 이하)에서는 시스템이 주기적인 정지-이동(stick‑slip) 현상을 보였으며, 이는 정지 마찰력(F_s)과 동적 마찰력(F_k)의 차이에 의해 구동력이 급격히 축적된 뒤 급방출되는 메커니즘으로 해석된다.

슬라이딩 과정에서 윤활제의 구조적 변화를 두 가지 모드로 구분한다. 첫 번째는 ‘액체 슬라이딩(Liquid Sliding, LS)’으로, 구동에 의해 윤활제 층이 완전히 녹아내려 비정질 액체 상태가 되면서 기판과의 마찰이 점성에 의해 지배된다. 두 번째는 ‘층 위 층 슬라이딩(Layer‑over‑Layer Sliding, LoLS)’으로, 윤활제 층이 고체‑유사한 구조를 유지하면서 서로 겹쳐진 층 사이에서 미끄러지는 형태이다. LoLS는 특히 저속 구간에서 우세하며, 이때 마찰계수는 미스핏 각도에 민감하게 변한다.

미스핏 각도(θ)는 두 기판 격자 사이의 회전각을 의미한다. θ=0°(완전 정렬)에서는 기판 격자와 윤활제 격자가 일치해 강한 포텐셜 ‘잠금’이 발생, 정지 마찰이 최대치에 도달한다. 반면 θ≈30°~45° 근처에서는 격자 간 비정합이 커져 잠금 효과가 약화되고, 정지 마찰이 급격히 감소한다. 특히 θ≈30°에서 LoLS가 가장 안정적으로 유지되며, LS 전이 없이도 낮은 구동력으로 지속적인 미끄럼이 가능함을 확인했다. 감쇠 계수의 위치·속도 의존성은 이러한 전이 현상을 더욱 뚜렷하게 만든다; 높은 상대 속도 구간에서는 감쇠가 크게 증가해 에너지 소산이 급격히 늘어나 LS로 전이되는 경향이 있다.

다결정 접촉을 모사하기 위해, 다양한 θ 값을 갖는 미시적 접촉점들의 정지 마찰 임계값을 통계적으로 수집하고, 이를 확률 분포 P(F_s)로 변환하였다. 결과는 θ가 무작위적으로 분포된 경우, P(F_s)는 평균값 주변에 좁은 피크를 보이지만, 큰 θ(≈0°)와 작은 θ(≈45°) 사이에 두 개의 부수 피크가 나타나며, 이는 특정 결정학적 정렬이 마찰을 크게 억제하거나 강화함을 의미한다. 이러한 분포는 실제 다결정 마찰면에서 관측되는 ‘마찰 강도 변동성’과 일맥상통한다.

요약하면, 본 연구는 (1) 비선형 감쇠를 포함한 랭게빈 MD가 경계 윤활 현상을 정밀히 재현함을, (2) 미스핏 각도가 정지·동적 마찰, LS↔LoLS 전이, 그리고 마찰 임계값 분포에 결정적 영향을 미친다는 점을, (3) 다결정 접촉에서 마찰 강도의 통계적 변동성을 정량적으로 설명할 수 있는 모델을 제시한다는 점에서 의미가 크다.