🔬 condensed matter

Voids in liquids: peculiarities of molecular dynamics simulation of fluid systems

이 논문은 분자 동역학 시뮬레이션에서 관찰되는 액체 내의 거대 공동이 임계 온도 이상이거나 액체 - 기체 2 상 영역에 있을 때 발생하는 현상임을 다양한 예시를 통해 규명하고 있습니다.

원저자: Yu. D. Fomin

게시일 2026-04-20

📖 3 분 읽기☕ 가벼운 읽기

원저자: Yu. D. Fomin

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

이 논문은 **"액체 안에 커다란 구멍 (공기 주머니) 이 생기는 이유"**에 대한 과학적 수수께끼를 푸는 이야기입니다.

일반적으로 우리는 액체 (물, 기름 등) 를 생각할 때 "분자들이 빽빽하게 모여 있어 구멍이 있을 수 없다"고 생각합니다. 하지만 컴퓨터 시뮬레이션 (분자 동역학) 을 해보면, 액체처럼 보이는 시스템 안에 **엄청나게 큰 구멍 (Void)**이 생기는 것을 종종 발견합니다. 마치 물속에 거품이 아닌, 완전히 빈 공간이 생기는 것처럼 보이죠.

저자 (유리 D. 포민) 는 이 현상이 "시뮬레이션 오류"가 아니라, 시스템이 실제로 어떤 상태에 있기 때문이라고 설명합니다. 핵심은 두 가지 상황입니다.

1. "초critical (Supercritical) 상태"의 액체: 뜨거운 안개와 얼음의 중간

비유: "뜨거운 찜질방 안의 안개"

보통 물은 차가우면 얼고, 따뜻하면 액체가 되며, 뜨거우면 기체 (수증기) 가 됩니다. 하지만 **임계점 (Critical Point)**이라는 아주 뜨거운 온도를 넘어서면, 액체와 기체의 구분이 사라집니다. 이를 '초임계 유체'라고 합니다.

상황: 이 상태에서는 액체처럼 보이지만, 사실은 액체와 기체가 뒤섞인 '혼돈' 상태입니다.
현상: 분자들이 뭉쳐서 작은 덩어리 (액체 같은 부분) 를 만들기도 하고, 그 사이사이로 거대한 빈 공간 (기체 같은 부분) 이 생깁니다.
해석: 시뮬레이션에서 큰 구멍이 보인다면, 그 액체가 너무 뜨거워서 더 이상 '진짜 액체'가 아니라, 액체와 기체가 섞인 '초임계 유체' 상태일 가능성이 큽니다. 마치 뜨거운 찜질방 안의 안개처럼, 물방울 (분자 덩어리) 과 수증기 (빈 공간) 가 공존하는 상태인 것입니다.

2. "두 상 (Phase) 공존" 상태: 물과 공기가 섞인 병

비유: "반쯤 채운 물병을 흔들면?"

액체와 기체가 공존하는 상태 (예: 끓는 물 위쪽의 수증기) 에서 컴퓨터 시뮬레이션을 할 때, 컴퓨터의 '벽' (주기적 경계 조건) 때문에 이상한 일이 일어납니다.

상황: 컴퓨터는 유한한 공간 (상자) 안에서 분자들을 움직입니다. 만약 액체와 기체가 섞인 상태를 시뮬레이션하면, 분자들이 뭉쳐서 '액체 방울'을 만들거나, 반대로 '액체 판'이 생기고 그 안에 '구멍'이 뚫리는 형태를 취합니다.
현상: 실제로는 액체와 기체가 층을 이루고 있어야 하는데, 컴퓨터 상자의 크기 제한 때문에 액체가 구멍을 뚫고 기체가 그 안으로 들어가는 듯한 모양이 됩니다.
해석: 큰 구멍이 보인다면, 그 시스템이 **진짜 액체가 아니라, 액체와 기체가 섞인 '혼합 상태'**에 있을 확률이 높습니다. 마치 물병에 물과 공기가 반반 섞여 있을 때, 물이 한쪽으로 몰리고 공기가 다른 쪽으로 몰리는 것처럼, 시뮬레이션 상자 안에서 액체가 뭉치면서 빈 공간이 생긴 것입니다.

구체적인 사례: 텔루륨 (Tellurium) 의 비밀

논문에서는 **텔루륨 (Te)**이라는 금속 원소를 액체 상태로 시뮬레이션한 연구들을 분석합니다. 이전 연구들에서는 텔루륨 액체 안에 전체 부피의 1/3 에 달하는 거대한 구멍들이 있다고 보고했습니다.

과거의 오해: "텔루륨 액체는 특이해서 구멍이 많구나!"라고 생각했습니다.
이 논문의 결론: "아니야, 텔루륨이 특이한 게 아니라, 시뮬레이션 조건이 잘못 설정되었을 뿐이야."
- 연구자들이 설정한 온도와 압력이, 텔루륨이 액체와 기체가 섞인 상태 (두 상 공존 영역) 에 해당했습니다.
- 즉, 텔루륨 액체 자체가 구멍을 만든 게 아니라, 시뮬레이션이 '끓는 물' 상태를 모방하다가 액체와 기체가 섞인 모습을 보여준 것입니다.

요약: 왜 이 논문이 중요한가?

이 논문은 과학자들에게 중요한 경고를 줍니다.

"컴퓨터 시뮬레이션에서 액체 안에 거대한 구멍이 보인다면, 그것은 액체의 성질이 아니라, 시뮬레이션이 '임계점 이상의 뜨거운 상태'이거나 '액체와 기체가 섞인 상태'를 보여주고 있다는 신호입니다."

한 줄 요약:
액체 속에 커다란 구멍이 보인다면, 그 액체가 '변태'한 게 아니라 너무 뜨겁거나 (초임계), 액체와 공기가 섞여서 (두 상 공존) 그렇게 보이는 것입니다. 시뮬레이션 결과를 해석할 때는 이 점을 꼭 기억해야 합니다.

논문 요약: 유체 시스템의 분자 동역학 시뮬레이션에서 관찰되는 공동 (Void) 의 기원

1. 연구 배경 및 문제 제기 (Problem)

배경: 분자 동역학 (Molecular Dynamics, MD) 시뮬레이션은 유체 시스템의 특성을 연구하는 강력한 도구로, 실험적으로 측정하기 어려운 구조적 특성이나 확산 계수 등을 계산할 수 있게 해줍니다.
문제: 일부 MD 시뮬레이션 연구에서 액체 내부에 거대한 공동 (voids/cavities) 이 형성되는 현상이 보고되었습니다. 이는 일반적인 액체의 상식 (밀집된 상태) 과 모순되며, 예를 들어 액체 텔루륨 (Tellurium) 의 경우 시스템 전체 부피의 26~35% 가 공동으로 차지한다는 결과가 나왔습니다.
가설: 이러한 거대한 공동이 실험적 한계 때문인지, 아니면 시뮬레이션 방법론의 인공적 산물 (artifact) 인지 명확한 설명이 부족했습니다. 본 연구는 이러한 거대 공동의 발생 원인을 규명하고, 시뮬레이션에서 공동이 관찰되는 열역학적 조건을 예측하는 것을 목표로 합니다.

2. 연구 방법론 (Methodology)

저자는 다양한 유체 시스템을 대상으로 MD 시뮬레이션을 수행하여 공동 형성 메커니즘을 분석했습니다.

시뮬레이션 시스템 및 모델:
- 2 차원 및 3 차원 Lennard-Jones (LJ) 유체: 임계점 부근 및 상변화 영역에서의 거동을 분석.
- 초임계 탄소 (Supercritical Carbon): GAP-20 포텐셜을 사용하여 고온 (4000~7000 K) 및 다양한 압력 조건에서 시뮬레이션.
- 물 (Water): TIP4P/2005 모델을 사용하여 NVT(정적적) 및 NPT(정압적) 앙상블에서 상변화 영역 시뮬레이션.
- 액체 텔루륨 (Liquid Tellurium): Vienna Ab initio Simulation Package (VASP) 를 이용한 ab-initio MD 시뮬레이션 수행 (기존 연구 [2] 의 결과 재현 및 검증).
분석 도구:
- LAMMPS (고전적 MD) 및 VASP (양자 역학적 MD) 사용.
- OVITO 소프트웨어를 통한 시각화 및 구조 분석.
- 클러스터링 분석, 퍼콜레이션 임계값 확인, 상태 방정식 (EOS) 비교.

3. 주요 결과 (Key Results)

연구 결과는 액체 내부에 거대한 공동이 관찰되는 경우, 해당 시스템이 진정한 액체 상태가 아님을 의미하며, 다음 두 가지 열역학적 상태 중 하나에 해당함을 규명했습니다.

A. 초임계 유체 (Supercritical Fluid, SCF) 영역

조건: 기체 - 액체 상전이의 임계 온도 ( $T_c$ ) 이상인 경우.
발견: $T > T_c$ $T > T_{c}$ 인 초임계 영역에서는 밀도가 높아도 (액체 밀도에 근접할지라도) 시스템이 불균질한 구조를 가집니다.
- 저밀도: 기체 매질 내 작은 입자 클러스터.
- 고밀도: 액체와 유사한 구조지만 내부에 거대한 공동 (voids) 이 존재.
예시: 2D LJ 유체 및 GAP-20 탄소 모델에서 $T > T_c$ 조건에서 밀도가 높을지라도 공동이 명확히 관찰됨을 확인했습니다. 이는 액체가 아닌 초임계 유체의 고유한 불균질성 (inhomogeneity) 때문입니다.

B. 기체 - 액체 2 상 공존 영역 (Two-phase Liquid-Gas Region)

조건: 임계 온도 이하 ( $T < T_c$ ) 이면서 기체와 액체가 공존하는 영역.
발견: 주기적 경계 조건 (Periodic Boundary Conditions, PBC) 하에서 시스템 크기가 유한할 때, 2 상 영역은 기체와 액체가 분리된 형태 (구, 원통, 판 등) 로 존재합니다.
- NVT 앙상블 (부피 고정): 밀도가 액체 영역에 가깝더라도 시스템 전체가 액체처럼 보이지만, 내부에 거대한 기체 공동 (cavity) 이 형성됩니다. 이는 실제 액체가 아니라 기체 - 액체 공존 상태의 한 형태입니다.
- NPT 앙상블 (압력 고정): 시스템이 자발적으로 기체 또는 액체 단일 상으로 수렴합니다.
시스템 크기 효과: 작은 시스템에서는 복잡한 구조 (구, 원통, 판) 가 관찰되지만, 열역학적 극한 (무한대) 에서는 액체 방울이 기체로 둘러싸인 형태가 되어야 합니다. 시뮬레이션에서 관찰되는 '액체 내 공동'은 2 상 공존 영역에서의 PBC 효과로 인한 인공적 구조일 가능성이 높습니다.

C. 액체 텔루륨 (Liquid Tellurium) 사례 분석

기존 연구 [2] 에서 액체 텔루륨 내부에 거대한 공동이 보고된 이유를 재분석했습니다.
해당 연구의 온도 조건 (560~970 K) 은 텔루륨의 끓는점 ( $T_b \approx 1260$ K) 에 가깝거나 그 이하이며, 밀도 값이 2 상 영역에 해당할 가능성이 높습니다.
저자가 수행한 NPT 시뮬레이션 결과, 해당 조건에서 시스템이 2 상 영역에 위치함을 확인했습니다. 따라서 기존 연구에서 관찰된 '공동'은 액체 내부의 결함이 아니라, 액체 - 기체 2 상 공존 상태에서의 기체 영역으로 해석되어야 합니다.

4. 주요 기여 및 결론 (Key Contributions & Conclusions)

공동의 기원 규명: MD 시뮬레이션에서 관찰되는 거대 공동은 액체의 고유한 특성이 아니라, ① 초임계 온도 ( $T > T_c$ ) 이상의 초임계 유체 상태이거나 ② 기체 - 액체 2 상 공존 영역에 있을 때 발생하는 현상임을 증명했습니다.
앙상블의 중요성 강조: 2 상 영역에서 NVT(부피 고정) 시뮬레이션은 인공적인 공동 구조를 생성할 수 있으나, NPT(압력 고정) 시뮬레이션은 시스템이 단일 상 (기체 또는 액체) 으로 수렴하게 하여 물리적 상태를 올바르게 반영함을 보였습니다.
기존 연구 재해석: 액체 텔루륨 등 복잡한 유체 시스템에 대한 기존 ab-initio MD 연구 결과 (거대 공동 존재) 가 잘못 해석되었을 가능성을 지적했습니다. 이는 실험적 밀도 값이 시뮬레이션 조건에서 2 상 영역에 해당할 수 있으며, 계산된 압력이 음수 (negative pressure) 를 보이는 것과 일치합니다.
실무적 시사점: 유체 시스템의 MD 시뮬레이션 결과에서 거대 공동이 관찰된다면, 해당 시스템이 액체 상태가 아니거나 (초임계 상태), 2 상 영역에 있음을 의미하므로 열역학적 조건 (온도, 압력, 밀도) 을 임계점 및 공존 곡선과 비교하여 신중하게 해석해야 함을 강조합니다.

5. 의의 (Significance)

본 논문은 분자 동역학 시뮬레이션 데이터 해석 시 발생할 수 있는 오해를 방지하고, 액체와 초임계 유체, 2 상 영역의 구조적 특성을 올바르게 구분하는 기준을 제시합니다. 특히, 실험적으로 접근하기 어려운 고온/고압 조건이나 ab-initio 시뮬레이션 결과에서 관찰되는 비정상적인 구조 (공동) 가 실제 물리적 현상인지 시뮬레이션의 인공적 산물인지를 판단하는 데 중요한 지침을 제공합니다.