Modeling oxygen-void interactions in uranium nitride

기술 요약: 우라늄 질화물 내 산소-공공 상호작용 모델링

문제 정의
우라늄 질화물(UN)은 높은 우라늄 밀도와 열전도율 덕분에 유망한 핵연료 후보 물질이다. 그러나 UN은 핵분열 가스 기포의 성장으로 인한 현저한 스웰링(swelling, 부풀음) 현상으로 인해 성능이 저하된다. 실험적 증거에 따르면, 산소 불순물은 특히 농도가 1000–1500 ppm을 초과할 때 스웰링을 크게 악화시킨다. UN의 스웰링을 설명하기 위한 메커니즘적 다중 스케일 모델(예: SIFGRS 모델)이 존재하지만, 이들은 일반적으로 산소가 없는 UN에 대해 도출된 일정한 표면 에너지 값을 사용한다. 산소 불순물이 보이드(void) 핵 생성 및 핵분열 가스 기포 성장에 영향을 미치는 근본적인 메커니즘은 여로 알려져 있지 않다. 구체적으로, 산소 편석이 보이드 및 핵분열 가스 기포의 표면 에너지를 어떻게 변화시키는지에 대한 정량적 이해가 부족하며, 이는 스웰링 거동을 예측하는 데 결정적인 매개변수이다.

방법론
본 연구는 UN 내 산소 불순물과 보이드 표면 사이의 상호작용을 정량화하기 위해 제1원리 기반의 열역학 모델을 개발한다. 이 접근 방식은 밀도 범함수 이론(DFT) 계산을 결합한 결함 화학 및 운동학적 모델링을 통합한다:

1. DFT 계산: Vienna Ab-initio Simulation Package(VASP)와 PBE 범함수를 사용하여, 벌크 및 (001) 표면에서의 다양한 산소 결함(치환형 $O_N$, 침입형 $O_i$, 그리고 $\{O_N:v_N\}$과 같은 클러스터)에 대한 형성, 삽입, 결합 및 이동 에너지를 계산한다. 곡선 형태의 보이드 표면을 근사하기 위해 평면 슬래브의 표면 에너지를 계산한다.
2. 결함 화학: 질량 작용 법칙에 기반한 열역학 모델을 구축하여 벌크 및 보이드 표면에서의 산소 결함 상대 농도를 결정한다. 이는 결함 반응, 화학 퍼텐셜(산화물 형성 에너지에 참조됨), 그리고 화학 양론적 상태(U-rich, 중간, N-rich)의 영향을 고려한다.
3. 운동학적 모델링: 서로 다른 이동 메커니즘($O_i$ 및 $\{O_N:v_N\}$)에 대한 산소 확산도를 계산한다. 핵 생성 단계 동안 보이드 표면에 도달하기 위한 산소 확산의 유한한 시간을 고려하기 위해 운동학적 보정 계수($\alpha$)를 도입한다.
4. 표면 에너지 감소 모델: 산소 흡착에 따른 표면 에너지 변화($\Delta\sigma$)를 도출한다. 이 모델은 불활성 가스 원자가 표면 에너지에 영향을 미치지 않는다고 가정하며, 보이드와 핵분열 가스 기포의 표면을 동일하게 취급한다. 분석은 이차 $UO_2$ 상의 형성을 배제하고 열역학적 요인만을 분리하기 위해 산소를 유일한 불순물로 간주하며 산소에 집중한다.

주요 결과
분석 결과, 치환형 산소가 표면 질소 자리에 편석되는 것($O_N^{(s)}$)이 표면 에너지 감소($|\Delta\sigma|$)의 주요 동력인 반면, 표면 홀(hollow) 자리에 있는 산소($O_i^{(s)}$)는 미미하거나 때로는 상쇄되는 역할을 한다는 것이 밝혀졌다.

• 온도 의존성: $|\Delta\sigma|$의 크기는 세 가지 뚜렷한 영역을 보인다:
  • 저온 영역: 산소 확산이 느려 보이드 표면으로 충분히 편석되지 못하므로 에너지 감소가 미미하다.
  • 중간 온도 영역: 산소의 삽입이 최적화되는 지점에서 최대 감소가 발생한다. 기준 사례(산소 1500 ppm, 기공도 5%, 보이드 반경 1 nm)의 경우, 최대 감소량($\approx -0.44$ J/m²)은 약 1524 K에서 나타난다. 이 온도 범위는 실험적으로 보고된 UN의 "브레이크어웨이 스웰링(breakaway swelling)" 시작 온도(1473–1573 K)와 상당 부분 겹친다.
  • 고온 영역: 산소의 벌크 용해도가 증가하여 표면 편석을 저해함에 따라 감소 효과가 감소한다.
• 매개변수 민감도:
  • 산소 농도: $|\Delta\sigma|$는 산소 함량에 따라 대략 선형적으로 비례한다.
  • 보이드 반경: 이 효과는 작은 공동(cavity, $R_v \lesssim$ 수 nm)에서 가장 두드하게 나타난다. 더 큰 보이드의 경우, 산소가 표면 에너지를 유의미하게 낮추기 위해서는 더 높은 온도가 필요하다.
  • 기공도: 기하학적 요인이 강한 의존성을 시사하지만, 기하학적 프리팩터(prefactor)가 표면 자격 농도 항에 의해 일부 상쇄되기 때문에 기공도 의존성은 제한적이고 완만하다.
• 전자 구조: 전하 밀도 차이 분석은 $O_N^{(s)}$가 음이온 격자 자리에 있는 배위가 부족한 우라늄 이웃들의 배위를 만족시키기 때문에 에너지적으로 선호됨을 확인시켜 준다. 반면 $O_i^{(s)}$는 좌절된 홀 자리 구성을 나타낸다.

의의 및 주장
본 논문은 산소 농도와 보이드 핵 생성 열역학 사이의 최초의 정량적 연결 고리를 제공한다고 주장한다. 본 연구의 주요 의의는 산소 오염된 UN에서 관찰되는 가속화된 스웰링을 설명하기 위해 UN의 메커니즘적 스웰링 모델을 재조정하는 데 있다:

1. 모델과 실험의 일치: 본 연구는 산소에 의한 표면 에너지 감소가 산소로 오염된 UN에서 관찰되는 가속화된 스웰링을 설명하는 데 필수적임을 입증한다. 계산된 표면 에너지 감소량(전형적인 기포 크기에 대해 약 0.30–0.37 J/m²)은 저자의 DFT 계산값(1.59 J/m²)과 기존의 성공적인 스웰링 모델에서 사용된 낮은 값(1.11 J/m²) 사이의 차이와 같은 차수(order)에 해당한다.
2. 메커니즘적 통찰: 결과는 산소가 보이드 엠브리오(embryo)와 핵분열 가스 기포를 안정화하여 핵 생성의 열역학적 장벽을 낮춘다는 것을 시사한다. 작은 보이드의 경우, 이 감소는 핵 생성 장벽을 약 62% 낮추어 실험적 스웰링 데이터와 일치하는 급격한 핵 생성율 증가를 초래한다.
3. 예측 능력: 본 프레임워크는 산소 함량과 가속화된 스웰링 사이의 연관성에 대한 물리적으로 근거 있는 설명을 제공하며, 이는 다중 스케일 연료 성능 코드에 통합될 수 있는 메커니즘을 제공한다.

저자들은 평면 (001) 표면을 곡선형 보이드에 적용한 가정, 고농도에서의 불순물-불순물 간 상호작용 무시, 탄소나 핵분열 생성물과 같은 다른 불순물의 제외 등 한계점을 인정한다. 그러나 그들은 산소 편석(질소 자리로의 산소 편석)이라는 지배적인 메커니즘이 이러한 단순화에도 불구하고 견고하게 유지된다고 주장한다.