MLIP 기반 산소 공공 형성 에너지 분석
양극재(LiCoO2) 안정성 예측 시뮬레이션 가이드
1. 계산 원리 및 공식
산소 공공(Oxygen Vacancy) 형성 에너지는 격자에서 산소가 이탈할 때 드는 에너지를 의미하며, 배터리의 열적 안정성과 전압 강하를 예측하는 핵심 지표입니다.
Evac = Edefective + ½EO₂ - Epristine
- Epristine: 결함이 없는 완전한 슈퍼셀의 에너지
- Edefective: 산소 원자 1개를 제거한 후 이완(Relax)된 구조의 에너지
- ½EO₂: 기체 상태 산소 원자 1개의 기준 화학 퍼텐셜
2. 완벽한 격자구조(Pristine) 확보 소스
| 데이터 소스 | 특징 및 장점 |
|---|---|
| Materials Project | 무기 화합물 DB의 표준. 최적화된 CIF/POSCAR 즉시 획득 가능. |
| ICSD / COD | 실험적으로 검증된 결정학 구조 데이터 확보에 유리. |
| ASE / Pymatgen | Python 라이브러리를 통해 단위 격자를 대형 슈퍼셀로 확장 생성. |
3. 예제 데이터: LiCoO2 슈퍼셀 (POSCAR)
아래는 산소 공공 계산을 위한 기본적인 슈퍼셀 구조 파일 예시입니다.
LiCoO2_Pristine_Supercell_Example
1.0
8.4312 0.0000 0.0000 # Lattice Vector X
-4.2156 7.3016 0.0000 # Lattice Vector Y
0.0000 0.0000 14.0500 # Lattice Vector Z
Li Co O
12 12 24 # 3x3x1 Supercell 구성 원자 수
Direct
0.00000 0.00000 0.00000 # Li 좌표들...
0.33333 0.66667 0.12500 # Co 좌표들...
0.33333 0.66667 0.24000 # O 좌표 (이 중 하나를 선택해 제거)
4. MLIP를 이용한 분석 프로세스
Step 1 구조 최적화 (Pristine Relaxation)
가져온 구조를 선택한 MLIP 모델(CHGNet, MACE 등)로 다시 한 번 최적화하여 Epristine을 구합니다.
Step 2 산소 결함 생성 (Defect Generation)
슈퍼셀 내부의 산소 원자 중 하나를 제거합니다. 층상 구조인 경우 다양한 위치의 산소를 제거해 보며 에너지를 비교합니다.
Step 3 결함 구조 최적화 (Defective Relaxation)
산소가 빠진 자리 주변 원자들이 재배치(Relaxation)되도록 MLIP로 계산하여 Edefective를 얻습니다.
Step 4 최종 에너지 도출
준비된 공식에 대입하여 형성 에너지를 계산합니다. MLIP는 이 과정을 DFT보다 1,000배 이상 빠르게 수행합니다.
💡 핵심 팁: 배터리 양극재는 충전/방전 상태에 따라 산소 공공 에너지가 달라집니다. Li 원자를 제거한 '탈리튬 구조'에서도 동일한 과정을 거쳐 안정성을 평가해 보세요.
