Dissociation Line and Driving Force for Nucleation of the Multiple Occupied Hydrogen Hydrate from Computer Simulation

문제 정의
수소(H₂) 하이드레이트는 가스 운송 및 저장의 유망한 대안으로, 금속 수소화물과 비교할 만한 부피 저장 용량을 유지하면서도 원자재 비용을 절감할 수 있는 잠재력을 가지고 있다. 그러나 수소 하이드레이트의 실질적인 적용은 상평형과 케이지 점유율을 결정하는 구체적인 요인에 대한 포괄적인 이해가 부족하다는 점에 의해 저해되고 있다. 하이드레이트 형성에 필요한 극고압 및 저온 조건 때문에 실험적으로 하이드레이트의 특성을 결정하는 것은 매우 까다로우나, 이론 및 시뮬레이션 접근 방식은 분자적 관점을 제공한다. 구조 II (sII) H₂ 하이드레이트에서 작은 케이지(D)와 큰 케이지(H)의 최적 점유율에 대해서는 여전히 논쟁이 지속되고 있다. 기존 연구들은 D 케이지의 이중 점유가 가능하다고 제안하는 반면, 다른 연구들은 그것이 구조를 불안정하게 만든다고 지적하며 상충되는 결과를 보여주었다. 또한, 기존 시뮬레이션 모델들은 다양한 압력 범위에서 물에 대한 H₂ 용해도와 3상 공존 온도($T_3$)를 동시에 정확하게 예측하는 데 어려움을 겪는 경우가 많다.

방법론
저자들은 H₂ 하이드레이트의 해리선을 결정하기 위해 GROMACS 소프트웨어 패키지를 이용한 분자 역학(MD) 시뮬레이션을 사용하였다. 물 분자는 TIP4P/Ice 포텐셜을 사용하여 모델링하였고, H₂ 분자는 수정된 Silvera 및 Goldman (SG) 모델을 사용하여 기술하였다. 방법론의 핵심적인 측면은 물의 산소와 H₂ 분자 사이의 교차 상호작용에 대한 베르텔롯(Berthelot) 결합 규칙의 수정이었다. 저자들은 온도 의존적 보정 계수 $\chi$를 도입하였으며, 상호작용 포텐셜을 정밀하게 조정하기 위해 새로운 파라미터 $b_0$를 제안하였다.

하이드레이트(H), 수용액 액체(Lw), 순수 H₂ 액체(LH₂) 상이 공존하는 3상 공존 온도($T_3$)를 결정하기 위해 두 가지 별도의 시뮬레이션 방법이 활용되었다:

1. 용해도 방법: 이 접근 방식은 두 개의 분리된 평형 상태를 시뮬레이션하였다: 순수 H₂(LH₂)와 접하고 있는 수용액 상(Lw–LH₂)과 H₂ 하이드레이트(H–Lw)와 접하고 있는 수용액 상(H–Lw). 일정한 압력(185 MPa)에서 두 시스템에 대한 수용액 내 H₂ 용해도를 온도의 함수로 도식화함으로써, 두 용해도 곡선의 교차점을 통해 $T_3$ 값을 얻었다. 이 방법은 1-1, 1-2, 2-2, 1-3, 1-4 구성(단일 및 다중 점유 포함)을 포함한 다양한 하이드레이트 점유율에 적용되었다.
2. 직접 공존 기법: 이 방법에서는 하나의 시뮬레이션 박스 안에 모든 세 상(H, Lw, LH₂)이 서로 접촉하도록 하였다. 고정된 압력(100, 185, 300 MPa)에서 온도를 변화시키며 계의 위치 에너지를 모니터링함으로써, 저자들은 하이드레이트 상이 성장하는지 또는 녹는지를 결정하였다. $T_3$는 하이드레이트가 성장하는 최고 온도와 녹는 최저 온도 사이의 범위로 식별되었다.

또한, 고전 핵생성 이론(CNT) 근사를 사용하여 과냉각 및 하이드레이트 점유율의 함수로서 핵생성의 구동력($\Delta\mu_N^{EC}$)을 계산하였다.

주요 결과

• 점유율이 $T_3$에 미치는 영향: 185 MPa에서의 용해도 방법을 사용했을 때, 본 연구는 하이드레이트 점유율이 해리 온도에 미치는 영향이 미미하다는 것을 발견하였다. 테스트된 모든 점유율(1-1, 1-2, 2-2, 1-3, 1-4)은 약 264–266 K의 오차 범위 내에서 유사한 $T_3$ 값을 나타냈다.
• 열역학적 안정성 및 핵생성: 유사한 $T_3$ 값에도 불구하고, 핵생성 구동력 분석은 상당한 열역학적 차이를 드러냈다. 작은 D 케이지에는 1개의 H₂ 분자가 있고 큰 H 케이지에는 3개의 H₂ 분자가 있는 구성(1-3 점유)이 가장 음(-)의 $\Delta\mu_N^{EC}$를 나타내어, 가장 열역학적으로 유리한 구성임을 보여주었다. 반대로, 작은 D 케이지의 이중 점유(2-2)는 강한 반발력으로 인해 가장 불리한 것으로 나타났으며, 가장 양(+)의 $\Delta\mu_N^{EC}$ 값을 가졌다.
• 모델 개선: Michalis 등이 제안한 원래의 베르텔롯 수정 방식($b_0=0$)은 H₂ 용해도는 정확하게 예측하였으나, 실험적 $T_3$ 값을 5–7 K 정도 과소평가하였다. 저자들은 $\chi$ 인자 식에 $b_0 = 0.1$이라는 새로운 파라미터를 도입하여 새로운 수정을 제안하였다. 이러한 조정을 통해 $T_3$ 예측값이 크게 개선되어 100, 185, 300 MPa 전 범위에서 실험 데이터와 매우 잘 일치하게 되었으나, 이는 H₂ 용해도의 약간의 과대평가를 초래하였다.
• 격자 상수: 본 연구는 sII 하이드레이트의 격자 상수($a$)가 온도와 점유율 모두에 따라 증가함을 관찰하였다. 특히, 2-2 점유는 작은 D 케이지에 갇힌 두 H₂ 분자 사이의 강한 반발력으로 인해 1-4 점유보다 더 큰 격자 상수를 나타냈다.

의의 및 주장
본 논문은 H₂ 하이드레이트의 점유율이 해리 온도와 핵생성 구동력 모두에 미치는 영향을 컴퓨터 시뮬레이션으로 결정한 첫 사례라고 주장한다. 저자들은 자신들의 연구가 가장 안정적인 화학 양론적 구조에 대한 모호성을 해결하였으며, 연구된 조건 하에서 1-3 점유가 열역학적으로 선호되는 상태임을 확인하였다고 강조한다. 이 결과는 H₂ 하이드레이트의 이론적 저장 용량이 3.2 wt.%임을 시사하는데, 이는 모든 케이지의 이중 점유와 관련된 5.0 wt.%의 용량보다는 낮지만 가장 안정적인 구성을 나타낸다.

또한 본 연구는 힘장(force field) 파라미터화의 중요성을 강조하며, 베르텔롯 결합 규칙에 대한 간단한 수정이 시뮬레이션된 용해도와 상평형 온도 사이의 불일치를 해소할 수 있음을 입증하였다. 저자들은 1-3 점유가 가장 안정적이지만, 시뮬레이션 파라미터(특히 $b_0$ 파라미터)의 선택을 통해 연구자의 특정 요구 사항에 따라 용해도의 정확한 예측 또는 해리선의 정확한 예측 중 하나를 우선시할 수 있다고 결론지었다.