에탄올 내 메탄·질소·산소·이산화탄소 헨리 상수의 분자 시뮬레이션 예측

초록

본 연구는 새롭게 개발한 강체 비극성 유니티드 원자 모델을 이용해 에탄올의 순수 물성(포화액밀도·증기압)을 재현하고, 윈도우 테스트 입자 삽입법으로 273–498 K 구간에서 메탄, 질소, 산소, 이산화탄소의 무한 희석 잔여 화학퍼텐셜을 계산해 헨리 상수를 예측한다. 전통적인 Lorentz‑Berthelot 결합 규칙(ξ = 1)으로는 실험값과 평균 20 % 오차(이산화탄소는 70 %)를 보였으며, 하나의 실험 혼합 데이터에 맞춘 수정된 결합 규칙(ξ ≠ 1)을 적용하면 오차를 크게 감소시킨다.

상세 분석

본 논문은 두 가지 주요 과업을 수행한다. 첫 번째는 에탄올 분자를 위한 새로운 강체 비극성(United‑Atom) 모델을 설계·검증하는 것이며, 두 번째는 이 모델과 기존에 검증된 2CLJQ(두‑중심 Lennard‑Jones + 점‑쿼드러플) 용질 모델들을 결합해 다양한 온도에서 헨리 상수를 예측하는 것이다.

에탄올 모델은 기존의 복잡한 다중‑점 전하·극성 모델과 달리, 메틸·메틸렌·히드록실 그룹을 각각 하나의 Lennard‑Jones 사이트와 전하(히드록실·메틸렌에 배치)로 단순화하였다. 양자화학(GAMESS, HF/6‑31G)으로 얻은 최적 기하구조를 기반으로, Ungerer 등(1995)의 AUA4 파라미터를 차용해 메틸·메틸렌의 σ·ε 값을 설정하고, 히드록실 사이트의 σ·ε와 전하, 그리고 히드록실 LJ 중심의 미세 오프셋을 실험 포화액밀도·증기압에 맞게 최적화하였다. 결과적으로 모델은 실험 증기압(오차 3.7 %), 포화액밀도(0.3 %), 기화열(0.9 %)에서 기존 최고 수준인 TraPPE‑UA(오차 8.4 %, 1.4 %, 1.6 %)보다 우수한 정확도를 보였다. 또한 임계점(Tc = 514 K, ρc = 5.94 mol L⁻¹, pc = 5.89 MPa)도 실험값에 근접하였다.

용질 모델은 이전 연구에서 순수 물성(포화액밀도·증기압)만을 이용해 파라미터화된 2CLJQ 형태이며, 메탄·질소·산소·CO₂ 모두에 대해 동일한 모델링 프레임워크를 적용했다. 이는 용질‑용질 상호작용을 정확히 재현함과 동시에, 용매‑용질 상호작용을 결합 규칙에 의존하도록 만든다.

용매‑용질 비결합 파라미터는 기본적으로 Lorentz‑Berthelot 규칙(σab = (σaa+σbb)/2, εab = √εaaεbb)을 사용한다. 이 경우 ξ = 1으로 설정해 ‘예측 모드’를 수행했으며, 273–498 K 구간에서 메탄·질소·산소는 평균 20 % 정도의 오차를 보였지만, CO₂는 70 %에 달하는 큰 편차를 나타냈다. 이는 비극성·극성 차이가 큰 CO₂와 에탄올 사이의 LJ ε 파라미터가 Lorentz‑Berthelot 규칙으로는 충분히 반영되지 않음을 시사한다.

이를 보완하기 위해 ξ를 하나의 실험 헨리 상수에 맞춰 최적화한 ‘조정 모드’를 도입하였다. ξ > 1(즉, εab를 증가)으로 설정하면 CO₂에 대한 오차가 크게 감소하고, 다른 세 용질도 실험값에 거의 일치한다. 이는 단일 상태 독립 바이너리 파라미터 하나만으로도 복잡한 용매‑용질 상호작용을 효과적으로 보정할 수 있음을 보여준다.

또한 논문은 다른 결합 규칙(Kohler, Hudson‑McCoubrey 등)이 주로 비극성 기체 혼합에 기반을 두고 있어, 현재 연구 대상인 극성 용매‑극성·비극성 용질 시스템에 적용하면 오히려 εab를 과소평가하게 됨을 지적한다. 따라서 수정된 Lorentz‑Berthelot 규칙이 가장 실용적이며, 추가적인 실험 데이터가 최소한 하나만 있으면 광범위한 온도·조성 영역에 적용 가능함을 입증한다.

전반적으로 이 연구는 (1) 계산 효율성이 높은 강체 비극성 에탄올 모델을 성공적으로 구축, (2) 기존 용질 모델과 결합해 무한 희석 상태에서의 헨리 상수를 정확히 예측, (3) 단일 바이너리 파라미터 조정을 통해 실험 데이터와의 일치를 크게 향상시킨다는 세 가지 핵심 성과를 제공한다. 이러한 접근은 산업용 혼합물 설계, 공정 시뮬레이션, 그리고 AIChE 2004 물성 예측 챌린지와 같은 대규모 베enchmark에 직접 활용될 수 있다.