Defect assignment of the clock site in $^{229}\text{Th:CaF}_2$

문제 정의
고체 상태 $^{229}$Th 핵시계의 성능은 호스트 결정 내 토륨 핵의 미시적 환경에 결정적으로 의존한다. $^{229}$Th:CaF$_2$는 높은 VUV 투과도와 충분한 도핑 밀도 덕분에 유망한 후보 물질이지만, 사중극자 분할(quadrupole-split)이 일어나는 지배적인 시계 활성 스펙트럼 사이트의 구체적인 결함 구성은 논쟁의 대상이 되어 왔다. 최근의 분광학적 연구는 이 사이트를 국부적인 전하 보상이 없는 토륨 이량체(thorium dimer, 인접한 두 개의 Th$^{4+}$ 이온)로 해석된 "D 센터"로 할당했다. 그러나 이 사이트의 조성과 구조는 시계 개발에 있어 매우 중요하다. 만약 전이가 Th-Th 응집체에 의해 호스트된다면, 선폭과 주파수 이동은 농도 의존적인 다중 도펀트 복합체에 묶이게 될 것이다. 반대로, 만약 스펙트럼이 국부적으로 전하가 보상된 단일 Th 사이트에서 기인하는 것이라면, 최적화 노력은 인접한 전하 보상 결함과 변형(strain)을 제어하는 데 집중되어야 한다.

방법론
저자들은 열역학적 추정, 밀도 범함수 이론(DFT) 계산, 그리고 전기장 구배(EFG) 비교를 결합하여 지배적인 결함 구성을 재평가하였다.

• 열역학: 본 연구는 제안된 Th-Th 이량체와 다양한 국부 전하 보상 단일-Th 모티프의 형성 자유 에너지를 평가한다. 여기에는 결합 에너지 계산과 관련 온도(300 K 및 CaF$_2$ 용융점 근처)에서의 구성 엔트로피 고려가 포함된다.
• DFT 계산: 저자들은 PBE 범함수와 PAW(projector-augmented-wave) 방법을 사용하는 VASP 코드를 사용하여 DFT 계산을 수행하였다. 이들은 치환된 Th$^{4+}$ 사이트 주변의 중성, 동일 화학양론적 2불소 침입형(2F$_i$) 모티프를 모델링하기 위해 5 $\times$ 5 $\times$ 5 슈퍼셀(377개 원자)을 사용하였다. 모든 구조에 대해 스핀-궤도 결합(SOC)이 포함되었다.
• 스펙트럼 비교: 연구는 릴랙스된 다양한 결함 기하 구조에 대한 주 전기장 구배 성분($V_{ZZ}$)과 비대칭 매개변수($\eta$)를 계산한다. 계산된 값들은 최근 문헌(예: Zhang et al., Morawetz et al., Hiraki et al.)에 보고된 실험적 전이 주파수 및 EFG 매개변수와 비교된다.

주요 기여 및 결과

1. Bare Dimer 가설의 기각: 열역학적 분석 결과, bare Th-Th 이량체의 형성은 에너지적으로 불리하다는 것이 나타났다. 분리된 결함으로부터 이량체를 형성하는 데 필요한 계산된 에너지는 +0.56 eV이며, 이는 인접한 Th$^{4+}$ 치환 사이의 반발력을 반영한다. 또한, 구성 엔트로피는 시계용으로 사용되는 희박 결정(농도 $\sim$0.02 at.%)에서 이량체 형성을 강력하게 저해한다(농도가 높을수록 자유 에너지 비용이 증가함). 이전 연구들에서 관찰된 Th 클러스터링은 시계 결정에서 사용되는 농도(0.02 at.%)보다 훨씬 높은 도핑 농도(1.1 at.%)에 기인한 것으로 간서된다.
2. 지배적 사이트 식별: 계산 결과, 국부적으로 전하가 보상된 단일-Th 결함이 열역학적으로 선호되는 구성임이 확인되었다. 구체적으로, 고립된 Th$^{4+}$가 Ca$^{2+}$ 자리에 위치하고 두 개의 인접한 불소 침입형에 의해 보상되는 90$^\circ$ 모티프(Th$^{\bullet\bullet}_{Ca} + 2F'_{i}$)가 두 침입형 가족 중 가장 낮은 에너지를 가진다.
3. 스펙트럼 일관성: 90$^\circ$ 2-F$_i$ 모티프는 107.650 V Å$^{-2}$의 계산된 $V_{ZZ}$ 값을 생성하며, 이는 106.3에서 109.1 V Å$^{-2}$ 사이의 실험값과 매우 잘 일치한다. 이는 지배적인 사중극자 분할 사이트가 이량체가 아닌, 국부적으로 보상된 단일-Th 중심과 일치함을 확인시켜 준다.
4. 부차적 사이트 할당: 본 연구는 "1st-2nd obtuse" 혼합 쉘 침입형 모티프를 다음으로 가능성 높은 구성으로 식별하였다. 이 구조는 더 높은 에너지를 갖는 약한 스펙트럼 사이트(참고문헌 [8]의 Site 3)에 대한 타당한 할당을 제공하며, 계산된 $V_{ZZ}$ 값은 -309.049 V Å$^{-2}$로 실험적 관찰과 일치한다.
5. 이전 불일치의 해소: 저자들은 90$^\circ$ 모티프에 대해 높은 $V_{ZZ}$(223 V Å$^{-2}$)를 예측했던 이전 계산들이 격자 벡터가 릴랙스된 작은 슈퍼셀을 사용한 결과로 인한 아티팩트(artifact)였을 가능성이 높다고 언급하였다. 현재의 더 큰 고정 격자 계산은 관측된 분기 규모(branch scale)를 회복하였으며, 이는 불일치가 불소 보상에 대한 증거가 아니라 수렴 문제였음을 나타낸다.

의의 및 주장
본 논문은 $^{229}$Th:CaF$_2$의 지배적인 시계 활성 사중극자 분할 사이트가 피할 수 없는 토륨 응집이 아니라 국부적인 불소 보상에 의해 제어된다고 주장한다. 이러한 재할당은 고체 상태 핵시계 개발에 직접적인 함의를 갖는다:

• 선폭 제어: 관찰된 농도 의존적 선폭은 중성인 국부 보상 결함 복합체가 탄성원(Kanzaki 힘)으로 작용하여 불균일 확장을 유도한다는 그림과 일치한다. 이는 쿨롱-유사 장거리 상호작용을 생성할 것으로 예상되는 uncompensated dimer의 경우와 대조적이다.
• 최적화 전략: 유용한 사중극자 분할 사이트의 비율을 높이고 선폭을 줄이기 위해서는 국부 불소 보상의 결함 화학을 엔지니어링하고, 변형을 관리하며, 방사선 유도 결함 전환을 완화하는 데 집중해야 한다. 이는 토륨 응집을 촉진하는 것이 올바른 방향이 아님을 시사한다.
• 미시적 기초: 본 연구 결과는 넓어짐(broadening) 메커니즘에 대한 미시적 기초를 제공하며, $^{229}$Th 삽입에 결함 보상이 필요하지 않은 대체 호스트(예: ThF$_4$)가 추가적인 개선을 제공할 수 있음을 시사한다.

저자들은 "O center"(zero-EFG 사이트)에 대해서는 신중한 태도를 유지하며, 그 기원이 원격 보상, 동적 평균화 또는 더 큰 결함 복합체와 관련된 별개의 문제임을 언급하였다. 또한 Site 4의 할당은 현재의 정적 모델 세트 내에서 아직 해결되지 않은 상태로 남아 있다.