분자동역학 시뮬레이션에서 접촉각 자동 추정 알고리즘

초록

본 논문은 분자동역학(MD) 시뮬레이션에서 액체 방울의 접촉각을 자동으로 추정하는 새로운 알고리즘을 제안한다. 마할라노비스 거리(MND)로 기체와 액체 분자를 구분하고, 4차 B‑스플라인 보간으로 2차원 밀도 그리드를 생성한다. 밀도 필터링을 통해 기-액 인터페이스를 추출한 뒤, 란다우 방법으로 원을 피팅하고 원의 방정식을 풀어 접촉각을 계산한다. 이 과정을 방울을 수직축을 중심으로 회전시키며 반복함으로써 전 시간에 걸친 접촉각 변화를 효율적으로 얻는다. 다양한 물 모델과 온도조절 방법에 적용해 검증하였다.

상세 요약

이 연구는 MD 기반 접촉각 측정의 두 가지 핵심 문제, 즉 (1) 기체와 액체 분자를 정확히 구분하는 방법과 (2) 인터페이스를 정밀하게 재구성하는 방법을 동시에 해결한다는 점에서 의미가 크다. 기존에는 밀도 임계값을 단순히 적용하거나, 시각화 이미지에서 수동으로 경계선을 추출하는 방식이 주를 이루었다. 그러나 이러한 방법은 잡음에 민감하고, 대규모 시뮬레이션에서 시간 효율성이 떨어진다. 저자들은 마할라노비스 거리(MND)를 이용해 각 입자의 위치와 전체 분포의 공분산 행렬을 고려함으로써, 비정상적인 분포를 가진 기체 분자를 효과적으로 배제한다. MND는 다변량 정규분포 가정 하에 이상치 탐지에 널리 쓰이며, 여기서는 액체 내부의 고밀도 영역과 저밀도 기체 영역을 통계적으로 구분하는 데 적합하다.

그 다음 단계에서는 4차 B‑스플라인 보간을 사용해 입자 데이터를 2D 격자에 매핑한다. 4차 스플라인은 2차 혹은 3차 보간보다 부드러운 연속성을 제공해, 특히 인터페이스 근처의 급격한 밀도 변화에서도 인공적인 진동을 최소화한다. 이때 격자 해상도와 스플라인 커널 폭을 적절히 선택하면, 계산 비용을 크게 늘리지 않으면서도 고해상도 밀도 프로파일을 얻을 수 있다.

밀도 필터링 단계에서는 사전 정의된 임계값을 적용해 액체와 기체 사이의 등밀도 면을 추출한다. 여기서 중요한 점은 필터링 파라미터가 물리적 의미(예: 평형 라플라스 압력)와 연결될 수 있다는 점이다. 추출된 인터페이스 점들을 란다우(Landau) 방법으로 원에 피팅하는데, 란다우 방법은 최소제곱 문제를 선형화하여 수치적 안정성을 확보한다. 원의 중심과 반지름을 구한 뒤, 고정된 고체 표면과의 접점에서 접선의 기울기를 계산해 접촉각을 도출한다.

특히 저자들은 방울을 수직축을 중심으로 다각도(예: 0°~360°)로 회전시키며 동일한 절차를 반복함으로써, 비대칭성이나 표면 이질성으로 인한 접촉각 변동을 포착한다. 이는 기존에 단일 평면에서만 측정하던 방식과 달리, 3차원적인 접촉각 분포를 시간에 따라 추적할 수 있게 한다.

알고리즘의 효율성은 두 가지 측면에서 검증된다. 첫째, MND와 B‑스플라인 보간은 O(N) 수준의 연산 복잡도를 유지한다(여기서 N은 입자 수). 둘째, 원 피팅과 회전 반복은 GPU 가속을 통해 실시간 수준으로 처리 가능하다. 결과적으로 수십 나노초에서 마이크로초까지의 시뮬레이션 구간을 자동으로 분석할 수 있다.

한계점으로는 (1) 마할라노비스 거리 계산 시 공분산 행렬이 충분히 샘플링되지 않으면 이상치 판정이 부정확해질 수 있다. (2) 인터페이스가 강하게 비구면(예: 타원형, 불규칙형)인 경우 원 피팅이 오차를 유발한다. (3) 온도와 압력에 따라 밀도 임계값을 동적으로 조정해야 하는데, 현재는 사용자 경험에 의존한다는 점이다. 이러한 점들은 향후 연구에서 다변량 커널 밀도 추정(KDE)이나 비선형 피팅(예: 타원, 베지어 곡선)으로 보완될 수 있다.

전반적으로 이 논문은 MD 시뮬레이션에서 접촉각을 자동, 정량적으로 추출하는 새로운 워크플로우를 제시함으로써, 습윤 현상, 캡릴러리 현상, 그리고 표면 에너지 계산에 필요한 데이터 획득을 크게 간소화한다는 점에서 높은 실용적 가치를 가진다.