Exact Breit-Rabi formulae for the nitrogen-vacancy center in diamond

기술 요약: 다이아몬드 내 질소-공석(NV) 센터를 위한 정확한 브레이트-라비(Breit-Rabi) 공식

문제 정의
다이아몬드 내 질소-공석(NV) 센터의 초미세 구조, 특히 질소 핵과의 상호작용은 양자 센싱 및 컴퓨팅 응용 분야에서 매우 중요하다. 브레이트-라비 공식은 알칼리 원자(전자 스핀 $S=1/2$)의 초미세 구조에 대한 정확한 해석적 해를 제공하지만, $S=1$인 전자 스핀을 가진 NV 센터에는 적용할 수 없다. 결과적으로, 축 방향 자기장 하에서의 NV 센터 초미세 제만 에너지 고유값(eigenvalues)과 고유상태(eigenstates)를 결정하는 작업은 역사적으로 수치적 방법이나 해석적 근사법에 의존해 왔다. 이러한 근사법들은 준위 반교차(level anticrossing, LAC) 부근이나, 단순 선형 표현이 무너지는 제로 및 초저자기장(ZLF) 응용 분야와 같은 영역에서는 불충분한 경우가 많다.

방법론
저자들은 $^{14}\text{N}$ (핵 스핀 $I=1$) 및 $^{15}\text{N}$ (핵 스핀 $I=1/2$) 동위원소 모두에 대해 NV 센터의 초미계 에너지 고유값 및 고유상태에 대한 정확한 폐쇄형(closed-form) 식을 유도한다.

1. 미세 구조 유도: 저자들은 임의의 자기장 하에서 NV 센터의 미세 구조를 먼저 다룬다. 좌표계를 회전시켜 가로 방향 자기장을 $x$축에 정렬함으로써, $3 \times 3$ 미세 구조 해밀토니안을 실수 대칭 행렬로 축소한다. 이 행렬의 특성 다항식이 3차임을 입증한다. 3차 방정식의 삼각 함수 해법(특히 실수 근에 맞게 조정된 카르다노의 방법)을 사용하여, 자기장 성분들에 대한 고유값 및 고 고유벡터의 정확한 식을 유도한다.
2. 초미세 유도 ($^{15}\text{N}$): $^{15}\text{N}$의 경우 ($S=1, I=1/2$), 전체 해밀토니안은 $6 \times 6$ 행렬이다. 저자들은 총 자기 양자수 $m_F = m_S + m_I$의 보존을 활용하여 해밀토니안을 블록 대각(block-diagonal) 형태로 분해한다. 이를 통해 문제를 독립적인 $1 \times 1$ 및 $2 \times 2$ 블록을 푸는 문제로 축소하여 표준적인 해석적 해법을 적용할 수 있게 한다.
3. 초미세 유도 ($^{14}\text{N}$): $^{14}\text{N}$의 경우 ($S=1, I=1$), 해밀토니안은 $9 \times 9$ 행렬이다. $^{15}\text{N}$의 경우와 달리 초기 행렬은 블록 대각 형태가 아니다. 저자들은 기저 상태(basis states)를 재정렬하기 위해 특정 치환 행렬을 적용하며, 이를 통해 $1 \times 1$, $2 \times 2$, 그리고 $3 \times 3$ 블록으로 구성된 블록 대각 구조를 드러낸다. 이 $3 \times 3$ 블록은 앞서 유도한 미세 구조 행렬과 수학적으로 동일하며, 이를 통해 앞선 삼각 함수 기반의 3차 해법을 재사용할 수 있다.
4. 근사식 회복: 저자들은 준위 반교차 지점 근처에서 자신들의 정확한 식에 대해 테일러 전개(Taylor expansion) 및 푸아소 전개(Puiseux expansion)를 수행한다. 이를 통해 Auzinsh 등이 제안했던 이차 근사식을 해석적으로 회복함으로써, 새로운 정확한 공식과 기존의 확립된 근사 모델 사이의 엄밀한 수학적 연결 고리를 제공한다.

주요 기여 및 결과

• 정확한 폐쇄형 해법: 본 논문은 두 종류의 $^{14}\text{N}$ 및 $^{15}\text{N}$ 동위원소에 대해 NV 센터의 초미세 에너지 고유값 및 고유상태에 대한 최초의 정확한 폐쇄형 해석적 식을 제공한다. 이 식들은 저자기장 제만 효과부터 고자기장 파셴-백(Paschen-Back) 영역까지 모든 자기장 영역에서 유효하다.
• 준위 반교차(LAC) 예측: 저자들은 바닥 상태(GSLAC) 및 들뜬 상태(ESLAC)의 준위 반교차에 대한 정밀한 조건을 유도한다.
  • $^{15}\text{N}$의 경우, GSLAC는 $B \approx 102.318$ mT에서 예측된다.
  • $^{14}\text{N}$의 경우, 서로 다른 하위 준위에 대한 별도의 GSLAC를 식별한다: $m_F = \pm 1$ 준위의 경우 $B \approx 102.522$ mT에서, $m_F = 0$ 준위의 경우 $B \approx 102.247$ mT에서 발생한다.
  • 유도된 공식들은 LAC 위치가 가로 방향 초미세 결합($A_\perp$)에 의존하는 항들에 의해 이동됨을 보여주며, 이는 제만 분리가 제로 자기장 분리(zero-field splitting)와 일치하는 지점에서 LAC가 발생한다는 기존의 가정을 교정한다.
• 검증: 유도된 해석적 곡선들을 점 단위의 수치 평가와 비교하였다. 저자들은 자신들의 정확한 공식이 float64 정밀도로 수치 결과와 일치함을 확인하였다.
• 근사법의 개선: 정확한 해를 전개함으로써, 자기장 및 에너지 준위의 2차 변화를 고려한 수정된 이차 근사식을 도출하였으며, 이는 LAC 근처에서 기존 근사법보다 우수한 정확도를 보여준다.

의의
본 논문은 이러한 정확한 공식들이 ZULF NMR 및 자력계와 같은 저자기장 응용 분야를 포함하여, NV 기반 양자 기술의 설계 및 해석에 필수적이라고 주장한다. 정확한 해법은 시뮬레이션에서 수치적 대각화를 수행할 필요성을 제거하며, 준위 반교차에서 결합 파라미터(예: 들뜬 상태 초미세 텐서 및 바닥 상태 유사 슈타르크(pseudo-Stark) 텐서)를 측정하기 위한 엄밀한 토대를 제공한다. 나아가, 본 연구는 9개 준위 시스템의 복잡성에도 불구하고, NV 센터 특유의 대칭성을 통해 알칼리 원자에서 발견되는 것과 유사한 정확한 해석적 해를 구할 수 있음을 입증하며 NV 초미세 시스템의 수학적 구조를 명확히 한다.