산소 공공(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 원자를 제거한 '탈리튬 구조'에서도 동일한 과정을 거쳐 안정성을 평가해 보세요.