고압·고온에서 광물·수용액 리튬 동위 원소 분별 예측: 효율적인 ab initio 방법

초록

본 연구는 고압·고온 환경(지각·맨틀)에서 리튬 동위 원소(⁶Li/⁷Li) 평형 분별 계수를 효율적으로 계산하는 새로운 ab initio 방법을 제시한다. 전체 진동 스펙트럼 대신 분별 원소에 작용하는 힘 상수만을 이용한 ‘단일 원자 근사’를 도입해 계산 비용을 크게 절감하면서도 실험값과 1‰ 이내의 정밀도를 달성하였다. 스투롤라이트, 스포듐, 다양한 종류의 마이카와 물‑리튬 용액 사이의 분별을 검증했으며, 원자 배위, Li‑O 결합 길이, 압축·열팽창이 분별에 미치는 메커니즘을 원자 수준에서 해석한다. 이 방법은 결정과 액체를 동일한 주기적 경계조건으로 다룰 수 있어 복잡한 시스템에도 적용 가능하다.

상세 분석

이 논문은 기존의 전통적인 정상모드 분석(NMA)이나 클러스터 모델링이 갖는 계산량 문제를 극복하고자, ‘단일 원자 근사(single‑atom approximation)’라는 새로운 접근법을 도입하였다. Bigeleisen‑Mayer 이론에 따르면, 온도 T가 충분히 높고(보통 T > 600 K) 진동수 ω가 2kBT/ħ 보다 작을 경우(즉 u = ħω/kBT < 2) 전체 진동 스펙트럼 대신 해당 원소에 작용하는 3개의 힘 상수(Ax,Ay,Az)만으로 β‑factor를 근사할 수 있다. 이는 β ≈ 1 + (Δm · ħ²)/(24 kB² T²) · ΣAi 형태로, 계산 복잡도가 N‑atom 시스템에 대해 O(N)에서 O(1)로 감소한다는 점에서 혁신적이다.

연구팀은 이 근사를 검증하기 위해 DFT 기반 CPMD 코드를 사용해 스투롤라이트, 스포듐, 그리고 여러 마이카 폴리타입(1M, 2M1, 2M2, 3T) 구조와 Li⁺이 용해된 물‑리튬 용액을 각각 5 Å 이상 간격의 주기적 셀로 모델링하였다. 고체는 최소 40원자, 액체는 64개의 물 분자와 1개의 Li⁺, 1개의 F⁻(전하 중화)로 구성하였다. 구조 최적화와 AIMD 시뮬레이션은 BLYP 교환‑상관 함수와 Goedecker 의사퍼텐셜을 사용했으며, 진동수와 힘 상수 계산을 위해 140 Ryd의 높은 에너지 컷오프를 적용하였다.

β‑factor는 두 가지 방식으로 산출되었다. 첫 번째는 전통적인 전체 진동 스펙트럼을 이용한 Eq. 1(완전 모드 분석)이며, 두 번째는 Eq. 3(단일 원자 근사)이다. 두 방법 모두 T > 600 K 구간에서 거의 일치했으며, 특히 1000 K에서 차이는 0.01 ‰ 이하로 미미했다. 이렇게 얻은 β‑factor를 이용해 α = β_A/β_B를 계산하면, 실험적으로 보고된 Li 동위 원소 분별(예: 스투롤라이트‑용액, 스포듐‑용액, 마이카‑용액)과 1 ‰ 이내의 차이만을 보였다.

또한, 힘 상수와 배위 환경을 분석함으로써 다음과 같은 메커니즘적 통찰을 얻었다. (1) Li⁺이 4‑배위(짧은 Li‑O 결합)일수록 힘 상수가 크게 나타나며, 무거운 ⁷Li가 선호되는 경향이 있다. (2) 압력이 증가하면 Li‑O 결합이 단축되고 힘 상수가 상승해 분별계수가 더욱 양(⁷Li‑rich)으로 변한다. (3) 고온에서는 물 분자와의 수소 결합 네트워크가 동적으로 재배열되면서 배위 수가 4→5→6으로 변동하고, 이는 평균 힘 상수의 감소와 함께 ⁶Li‑선호 효과를 강화한다. 이러한 원자‑수준의 설명은 기존 실험만으로는 파악하기 어려웠던 ‘압축·열팽창·배위 변화가 동위 원소 분별에 미치는 복합 효과’를 명확히 해준다.

마지막으로, 이 방법은 주기적 경계조건을 사용해 고체와 액체를 동일한 양자역학적 프레임워크로 다루므로, 복잡한 광물·용액 시스템(예: 다중 원소, 다상 반응)에도 확장 가능하다. 계산 비용이 크게 절감되면서도 정확도가 유지되므로, 향후 지구화학·지구물리학에서 고압·고온 동위 원소 모델링의 표준 도구로 자리 잡을 전망이다.