Accurate Solvation Properties in supercritical CO$_2$ with Molecular Density Functional Theory

기술 요약: 분자 밀도 범함수 이론을 이용한 초임계 이산화탄소 내 정확한 용매화 특성 연구

문제 제기
초임계 이산화탄소(scCO2)는 추출에서 나노입자 합성까지 다양한 산업적 응용 분야에서 매우 효율적이고 환경 친화적인 용매이다. 이러한 공정을 최적화하려면 다양한 열역학적 조건에 걸쳐 용매화 자유 에너지 및 용매화 구조와 같은 용매화 특성을 정확하게 예측하는 것이 필요하다. 분자 역학(MD) 시뮬레이션은 높은 예측 정확도와 미시적 통찰력을 제공하지만, 고급 샘플링 전략(예: Widom 삽입법, 열역학적 적분)을 통한 열역학적 성질 계산 시 과도한 계산 비용이 발생한다는 단점이 있다. 반대로, 기존의 파라미터 기반 방법(상태 방정식, 밀도 기반 또는 용매화 모델)은 계산 비용은 낮지만, 새로운 용질이나 시스템 변경에 따라 재파라미터화가 필요할 정도로 전이성(transferability)이 제한적이다. 따라서 MD의 정확성과 대규모 파라미터 연구에 필요한 계산 효율성을 결합한 방법론이 필요하다.

방법론
본 연구는 분자 유체에 적응된 고전 밀도 범함수 이론(cDFT) 형식인 분자 밀도 범함수 이론(MDFT)을 채택한다. 이 접근 방식은 용매 밀도 $\rho(\mathbf{r}, \omega)$와 용매화 자유 에너지를 결정하기 위해 용매화 범함수 $\Delta\Omega_{solv}[\rho]$를 최소화한다. 범함수는 세 가지 성분으로 분해된다:

1. 이상적 기여 ($\Delta\Omega_{id}$): 용매 밀도의 엔트로피를 설명한다.
2. 외부 기여 ($F_{ext}$): 표준 포스 필드(CO2를 위한 EPM2, 분자 용질을 위한 OPLS-AA, 물을 위한 SPC/E)로부터 유도된 용질-용매 상호작용을 나타낸다.
3. 과잉 기여 ($\Delta\Omega_{exc}$): 용매-용매 상호작용을 나타낸다. 이 항은 균일 참조 유체(Homogeneous Reference Fluid, HRF) 근사를 사용하여 근사화되는데, 이는 용질 근처의 밀도 요동이 균일한 벌크 상태의 요동과 유사하다고 가정한다. HRF 항은 사전 MD 시뮬레이션 및 적분 방정식 이론으로부터 얻은 순수 벌크 용매의 직접 쌍 상관 함수(DCF)에 의존하는 이차 전개식이다.

HRF 근사의 한계(고차 요동을 무시함)를 해결하기 위해, 저자들은 브리지 범함수($\Delta\Omega_{bridge}$, 정확한 과잉 범함수와 HRF 항 사이의 차이)에 대해 세 가지 수준의 근사법을 조사하였다:

• 브리지 없음 (No Bridge): $\Delta\Omega_{bridge} = 0$.
• 하드 스피어 브리지 (Hard-Sphere Bridge, HSB): 벌크 자유 에너지에 맞춰 조정된 하드 스피어 과잉 범함수(FMT)를 사용하여 브리지를 근사한다.
• 조립형 브리지 (Coarse-Grained Bridge, CGB): 벌크 상태 방정식의 다항식 전개로부터 유도된 조립형 밀도장(coarse-grained density field)의 함수로 브리지를 모델링한다.

MDFT 계산은 공간 격자 $0.3$Å, 최대 각도 지수$n_{max}=4$를 가진 입방체 박스($L=30$ Å)에서 수행되었다. 이 결과들은 4000개의 CO2 분자와 베넷 수락 비(Bennett acceptance ratio) 방법을 사용한 MD 시뮬레이션과 벤치마킹되었다.

주요 기여

• 형식론 적용: 본 연구는 공간적 위치와 분자 방향 모두에 의존하는 분자 cDFT 형식을 scCO2에 성공적으로 적용하였다.
• 브리지 범함수 파라미터화: 저자들은 scCO2를 위한 두 가지 구별된 브리지 범함수(HSB 및 CGB)를 구축하고 파라미터화하였다. 이들은 조립형 격리 길이 $\sigma_g = 1.45$ Å를 가진 CGB 범함수가 상태 방정식으로부터 유도된 균일 브리지 기여를 정확하게 재현함을 입증하였다.
• 계산 효율성: 본 연구는 MDFT의 계산적 이점을 정량화하였는데, 단일 계산에 노트북에서 1분 미만(약 0.01 CPU·시간)이 소요되는 반면, 동일한 MD 시뮬레이션에는 약 640 CPU·시간이 소요됨을 명시하였다.

결과
MDFT 프레임워크는 구형(아르곤 유사) 용질과 9종의 분자 용질(물, 에탄올, 톨루엔, CO2 자신 포함)에 대해 MD 벤치마크 테스트를 거쳤으며, 임계점 근처의 상태 지점($n_0 = 0.8 n_c$, $T = 1.05 T_c$)에서 수행되었다.

• 구조적 특성: 모든 용질에 대해, MDFT가 예측한 방사 분포 함수(RDF)는 MD 결과와 거의 완벽하게 일치하는 중첩을 보였다. 이러한 일치는 어떤 브리지 근사(HRF, HSB, CGB)를 사용하더라도 유지되었다.
• 용매화 자유 에너지:
  • 구형 용질의 경우, HRF 근사(브리지 없음)는 큰 용질($\sigma > 7$ Å)에 대해 용매화 자유 에너지를 체계적으로 과대평가하였다. 그러나 HSB 및 CGB 보정은 넓은 크기와 상호작용 에너지 범위에 걸쳐 MD 벤치마크와 구별할 수 없을 정도로 일치하는 결과를 냈다.
  • 분자 용질의 경우, 세 가지 MDFT 변형(HRF, HSB, CGB) 모두 평균 상대 오차 약 1% 내외로 MD 용매화 자유 에너지를 재현하였다.
  • 특히, scCO2 내의 분자 용질에 대해 브리지 범함수의 기여는 무시할 수 있는 수준인 것으로 나타났다. 즉, HRF 근사만으로도 열역학적 성질을 포착하기에 충분했으며, 이는 액체 상태의 물에 대한 이전의 MDFT 연구에서 브리지 항이 필수적이었던 것과 대조적이다. 이는 scCO2에 강력한 방향성 수소 결합 네트워크가 부재하기 때문인 것으로 분석된다.
• 정확도 대비 비용: MDFT 접근법은 구조적 및 열역학적 특성 모두에서 정량적 수준의 정확도를 달와는 동시에, MD 대비 약 60,000배의 계산 효율성을 달성하였다.

의의 및 주장
본 논문은 MDFT가 scCO2의 용매화를 모델링하는 데 있어 견고하고 계산 효율적인 대안임을 주장한다. 주요 의의는 벌크 용매의 상관 관계가 정확히 알려져 있다면, HRF 근사가 충분히 MD 벤치마크 결과를 회복할 수 있음을 입증했다는 점에 있다. 이 발견은 물과 같은 다른 유체에서 흔히 요구되는 복잡한 브리지 보정이 scCO2의 특정 유체 특성 덕분에 불필요할 수 있음을 시사한다.

저자들은 현재의 연구가 전체 상평면에 대한 전이성을 확립하는 것이 아니라, 대표적인 근임계 상태 지점에서 방법을 검증하는 것임을 명시하였다. 그들은 이 프레임워크가 계면 및 구속 조건과 같은 더 복리한 시나리오로 MDFT를 확장할 수 있는 토대를 제공하며, 향후 개발에서는 서로 다른 열역학적 조건에 대한 전이성을 달성하기 위해 머신러닝 기반 범함수를 활용할 수 있다고 결론지었다. 식별된 즉각적인 효용성은 빠르고 정확한 용매화 특성 스크리닝을 통해 scCO2 기반 산업 공정의 합리적 설계 및 최적화를 촉진할 수 있는 잠재력이다.