Microstructural Insights into Fast Ion Transport in Solid Electrolytes via Multiscale Modeling

기술 요약: 다중 스케일 모델링을 통한 고체 전해질 내 빠른 이온 수송의 미세구조적 통찰

문제 정의
고성능 전고체 배터리 개발은 고체 전해질 내에서 빠른 이온 수송을 가능하게 하는 미세구조적 특징에 대한 불완전한 이해로 인해 저해되고 있다. 역사적인 노력들은 벌크 초이온 전도도를 최적화하는 데 집중해 왔으나, 실제 고체 전해질은 다결정 형태이며 밀집된 결정립계(GB) 네트워크를 포함하고 있다. 산화물 전해질의 경우 결정립계 효과가 부분적으로 이해되어 있으나, 황화물 전해질에서의 결정립계 역할은 형태학적 특성을 규명하고 결정립계의 기여도를 벌크 전도와 분리하기 어려운 실험적 어려움으로 인해 여전히 불분명하다. 또한, 다결정 내 이온 수송을 모델링하는 것은 계산적으로 매우 까다롭다. 이는 원자 단위의 호핑(hopping) 메커니즘과 연속체 규모의 퍼콜레이션(percolation) 사이를 연결해야 하며, 복잡한 결정립계 상호작용을 정확하게 묘사할 수 있는 상호작용 포텐셜과 구조적 무질서를 포착할 수 있을 만큼 충분히 큰 시뮬레이션 셀을 요구하는데, 전통적인 ab initio 방식으로는 적절한 스케일에서 이를 달iling할 수 없다.

방법론
저자들은 밀도범함수이론(DFT), 머신러닝 포텐셜(MLP), 분자 동역학(MD), 그리고 유한요소(FE) 시뮬레이션을 통합하는 계층적 다중 스케일 모델링 프레임워크를 제안한다.

1. 머신러닝 포텐셜 (MLPs): 정확도와 계산 비용 사이의 균형을 맞추기 위해, 저자들은 폐쇄 루프 능동 학습(active learning) 스킴을 사용하여 경량 모멘트 텐서 포텐셜(MTP)을 개발하였다. 이 과정은 결정립계를 포함하는 전체 미세구조로부터 국부적인 원자 환경을 추출하고, 이를 바탕으로 DFT 계산을 수행하며, MTP가 근사 ab initio 정확도에 도달할 때까지 복잡성을 높여가며 반복적으로 학습시키는 과정을 포함한다.
2. 분자 동역학 (MD): 훈련된 MLP는 바이크리스탈(bicrystal), 컬럼형 결정(columnar crystal), 그리고 무작위 결정립계 네트워크를 가진 다결정체를 포함하는 대규모 대리 결정립계 모델에 적용되어 나노초 단위의 MD 시뮬레이션을 수행하는 데 사용되었다. 이 시뮬레이션은 300 K에서 600 K 사이의 온도 범위에서 Li-이온의 무질서, 상관관계 및 확산 역학을 포착하였다. 주요 지표로는 자기 확산도(self-diffusivities), 확산 장벽, 그리고 구조적 무질서를 정량화하기 위한 과잉 엔트로피(excess entropy)가 포함되었다.
3. 유한요소 (FE) 시뮬레이션: MD에서 도출된 자기 확산도와 결정립계 폭은 연속체 규모의 FE 시뮬레이션을 위한 입력값으로 사용되었다. 스케일 분리를 처리하기 위해 두 가지 접근법이 활용되었다:
   • 해상된 결정립계 (Resolved GBs): 나노 크기의 결정립(~25 nm)에 대해, 푸리에 가속 노달 솔버(FANS)를 사용하여 결정립계를 부피 기반으로 해상하였다.
   • 축소된 결정립계 (Collapsed GBs): 마이크로 크기의 결정립(~1 µm)에 대해, 결정립계를 두께가 축소된 유한 요소로 표현하였다.
     이 접근법을 통해 소자 관련 결정립 크기를 가진 다결정 시스템의 거시적 확산도 및 전도도를 예측할 수 있었다.

주요 기여 및 결과
본 연구는 아르지로다이트(argyrodite) 고체 전해질인 $\text{Li}_6\text{PS}_5\text{X}$ ($\text{X} \in \{\text{Cl, Br, I}\}$)에 적용되어, 음이온 정렬, 할로겐 치환 및 결정립계 형태의 영향을 체계적으로 조사하였다.

• 결정립계 에너지 및 구조: 음이온이 정렬된 $\text{Li}_6\text{PS}_5\text{Cl}$에서 결정립계 형성 에너지는 일반적인 결정립계의 경우 12–22 $\text{meV \AA}^{-2}$, $\Sigma3$의 경우 약 5 $\text{meV \AA}^{-2}$로 낮게 나타났으며, 이는 합성된 샘플 내에 풍부하고 무작위적인 결정립계가 존재함을 시사한다. 결정립계는 과잉 엔트로피 증가와 서브 나노미터 크기의 공극(void)을 통해 구조적 무질서를 유도하지만, 이러한 공극이 확산도에 미치는 영향은 미미했다.
• 음이온 정렬 및 치환 효과:
  • 음이온 무질서: 음이온이 무질서한 $\text{Li}_6\text{PS}_5\text{Cl}$에서 결정립계는 장벽으로 작용하여 Li-이온의 자기 확산도를 벌크 대비 약 60% 감소시킨다. 그러나 특징적인 호핑 빈도는 여전히 높은 수준($\sim10^8 \text{ Hz}$)을 유지하며, 이는 실험적으로 관찰되는 낮은 결정립계 저항과 일치한다.
  • 할로겐 치환: 정렬된 상(phase)에서 Li-이온 확산 장벽과 음이온 반경 사이의 선형 스케일링 관계가 확립되었다. 음이온 반경이 커질수록 (Cl $\rightarrow$ Br $\rightarrow$ I) 확산 장벽은 0.1 Å 증가당 약 80 meV씩 증가한다.
  • 상반된 결정립계 효과: 결정립계의 영향은 벌크의 특성에 따라 달라진다. 빠른 확산을 보이는 상(음이온 무질서 $\text{Li}_6\text{PS}_5\text{Cl}$)에서는 결정립계가 확산을 억제한다. 반면, 느린 확산을 보이는 상(음이온 정렬 $\text{Li}_6\text{PS}_5\text{I}$)에서는 결정립계가 확산을 최대 60배까지 향상시키며 연결 채널 역할을 한다.
• 비-아레니우스(Non-Arrhenius) 수송: 시뮬레이션은 음이온 정렬된 $\text{Li}_6\text{PS}_5\text{I}$에서 관찰되는 비-아레니우스 수송 거동을 성공적으로 재현하였으며, 수송 메커니즘이 벌크 지배(고온)에서 결정립계 지배(저온)로 전환되는 전이 온도(~250 K)를 식별하였다. 이러한 거동은 이 특정 시스템에서 결정립계의 활성화 에너지가 벌크보다 낮기 때문인 것으로 분석된다.
• 결정립 크기 효과: 연속체 시뮬레이션은 결정립 크기를 나노 규모로 줄이면 비-아레니우스 거동이 완화됨을 보여주었다. $\text{Li}_6\text{PS}_5\text{Cl}$의 경우, 결정립 미세화는 초이온 전도도를 저해하지 않으면서 입계 간 접촉을 개선한다. 반면 $\text{Li}_6\text{PS}_5\text{I}$의 경우, 나노 크기화는 고유한 초이온 특성이 부족한 시스템에서 수송을 활성화하는 전략이 될 수 있다.

의의 및 주장
본 논문은 자사의 다중 스케일 모델링 프레임워크가 원자 단위의 메커니즘과 거시적 소자 성능 사이의 엄격하고 정량적인 연결 고리를 제공하며, 다결정 고체 전해질 모델링의 오랜 과제를 해결한다고 주장한다.

• 메커니즘적 통찰: 본 연구는 결정립계가 보편적으로 이온 수송을 방해하는 것이 아니라, 결정립계의 효과(강화 또는 억제)가 벌크 상의 고유한 확산도에 의해 결정된다는 점을 명확히 하였다.
• 설계 가이드라인: 연구 결과는 미세구조 설계에 대한 사전(a priori) 지침을 제공한다. $\text{Li}_6\text{PS}_5\text{Cl}$과 같은 빠른 확산형 황화물의 경우, 전도도를 희생하지 않으면서 접촉을 개선하기 위해 결정립 미세화가 유리하다. $\text{Li}_6\text{PS}_5\text{I}$와 같은 느린 확산형 시스템의 경우, 나노 크기화가 수송을 활성화하는 전략이 될 수 있다.
• 조성 파라미터: 본 연구는 음이온 반경이 정렬된 아르지로다이트의 확산 장벽을 결정하는 핵심 조성 파라미터임을 밝혀냈으며, 이는 격자 편극성(lattice polarizability)만이 유일한 지배적 요인이라는 기존의 관점에 도전한다.
• 실험적 검증: $\text{Li}_6\text{PS}_5\text{I}$에서 예측된 비-아레니우스 거동과 전이 온도는 실험적 관찰 결과와 일치하며, 이는 결정립계가 이 물질들의 이러한 현상의 기원임을 시사한다.

저자들은 결론적으로, 이 접근법이 이상적인 벌크 모델을 넘어 실제적인 다결정 표현으로 확장함으로써 미세구조적 복잡성을 명시적으로 고려한 차세대 고체 전해질 설계를 가능하게 한다고 밝히고 있다.