Triangulene-based diradicals as a blueprint for molecular quantum platforms with optical addressability and long spin coherence times

기술 요약: 분자 양자 플랫폼으로서의 트리앙굴렌 기반 디라디칼

문제 제념
분자 양자 기술의 발전을 위해서는 긴 스핀 결맞이 시간(coherence time)과 효율적인 스핀-광학 인터페이스를 동시에 보유하며, 이상적으로는 상온에서 작동 가능한 분자를 식별하는 것이 필요하다. 다이아몬드의 질소-공석(NV) 센터와 같은 고체 상태 결함은 광학적으로 주소 지정 가능한 스핀 초기화, 광 검출 자기 공명(ODMR)을 통한 판독, 그리고 긴 결맞이 시간을 나타내며 이러한 특성의 기준점을 제공하지만, 분자 시스템에서 이러한 기능을 재현하는 것은 여전히 큰 과제로 남아 있다. 현재의 분자 후보군들은 잘 격리된 삼중항 바닥 상태, 높은 화학적 안정성, 그리고 광학적 제어에 필수적인 고도로 스핀 선택적인 계간 전이(intersystem crossing, ISC)를 결합하는 데 어려움을 겪는 경우가 많다. 또한, 개방형 껍질(open-shell) 유기 시스템에서 바닥 상태 스핀 결맞음 해제(decoherence)와 들뜬 상태의 스핀 선택적 역학을 모두 설명할 수 있는 통합된 제일원리 이론 프레임워크가 그동안 부족했다.

방법론
저자들은 트리앙굴렌(1), 그 아자-양이온 유도체(2), 그리고 2,6,10-tri-tert-butyl-4,8,12-trimesyl-triangulene 결정(3)의 전자 구조와 역학적 특성을 조사하기 위해 고급 ab initio 제일원리 방법을 사용한다.

• 전자 구조: 다중 참조 양자 화학 시뮬레이션(NEVPT2 동적 상관 보정이 포함된 CASSCF)을 사용하여 스핀 자유 전자 상태의 에너지 층계를 결정한다. 연구에서는 최대 22개의 전자와 19개의 궤도를 갖는 거대 활성 공간을 처리하기 위해 ICE-CI 솔버를 사용하며, 상세 분석을 위해 NEVPT2(10,10) 활성 공간으로 축소하여 사용한다.
• 스핀-포논 역학: 스핀-격자 완화($T_1$) 및 순수 탈위상($T_2$)은 상대론적 진동 결합(vibronic coupling) 프레임워크를 사용하여 모델링된다. 해밀토니안은 비상대론적 항, 스핀-궤도 결합(SOC), 스핀-스핀 결합(SSC), 그리고 제만(Zeeman) 상호작용을 포함한다. 진동 결합은 핵 변위에 대한 전자 해밀토니안의 테일러 전개(선형 및 이차 항)를 통해 처리된다. 전이율은 마르코프 배스 근사(Markovian bath approximation) 내에서 일포논 및 이포논 과정(흡수, 방출, 결합 과정)을 고려한 양자 마스터 방정식을 통해 계산된다.
• 스핀-스핀 탈위상: 주변 핵 스핀에 의한 결맞음 해제는 클러스터 상관 확장(Cluster Correlation Expansion, CCE) 방법(pyCCE 이용)을 사용하여, 중심 전자 스핀이 수소 또는 중수소 핵 스핀의 배스와 상호작용하는 모델을 통해 시뮬레이션된다.
• 들뜬 상태 역학: ISC 비율은 중간 전자 상태로의 가상 전이를 매개로 하는 이차 포논 기여를 포함하며, 선형 및 이차 진동 결합을 통합하여 계산된다.

주요 기여 및 결과

1. 전자 구조 및 NV 센터와의 유사성:
   본 연구는 트리앙굴렌(1)과 그 아자-양이온 유도체(2)가 약 0.5 eV 정도의 간격으로 첫 번째 단일항 들 excited state($S_0$)와 잘 분리된 삼중항 바닥 상태($T_0$)를 보유하고 있음을 확인한다. 이 에너지 간격은 NV 센터와 매우 유사하다. 아자-양이온(2)은 1에 비해 감소된 삼중항-삼중항 에너지 간격을 나타내며, 들뜬 삼중항 상태($T_1, T_2$)는 NV 센터의 들뜬 상태와 유사하게 위치한다. 결정 유도체(3)의 전자 구조는 1 및 2와 유사하지만, 결정 패킹이 대칭성을 약간 낮추고 에너지 간격을 줄이는 것으로 나타났다.

2. 스핀-포논 완화 ($T_1$):
   계산 결과, 고립된 분자 1과 2의 상온(300 K)에서의 스핀-포논 완화 시간은 매우 길 것으로 예측되었다 ($T_1 \approx 1$ ms for 1, $T_1 \approx 27$ ms for 2). 2에서 $T_1$이 더 긴 이유는 진동 결합 강도가 감소했기 때문이다. 그러나 기능화 및 결정 격자 내 포함(화합물 3)은 300 K에서 $T_1$을 $\sim 0.1$ $\mu$s로 급격히 감소시킨다. 이러한 감소는 부피가 큰 치환기(tert-butyl 및 mesityl)와 관련된 저에너지 포논 및 격자 상호작용에 의해 구동된다.

3. 스핀-스핀 탈위상 ($T_2$):
   핵 스핀에 의한 결맞음 해제를 시뮬레이션한 결과, 수소화된 환경에서 $T_2$는 마이크로초 단위로 제한된다 (예: 3a-H의 경우 $T_2 \approx 2.42$ $\mu$s). 그러나 본 연구는 화학적 엔지니어링을 통해 이를 크게 완화할 수 있음을 보여준다. 중심 분자가 중수소화되고 핵 스핀이 없는 환경(예: 중수소화된 용매 또는 매트릭스)에 놓이는 가상의 시나리오에서, $T_2$는 10 K에서 0.21 ms로 10배 이상 증가한다.

4. 스핀 선택적 계간 전이 (ISC):
   분자 2에 대한 300 K에서의 이차 포논 ISC 비율 계산은 고도로 스핀 선택적인 과정을 드러낸다. 연구는 ISC의 스핀 선택성이 상태의 에너지 정렬과 분자 기하 구조에 민감함을 확인했다.

   • 바닥 상태 기하 구조에서, $T_1$에서 $S_2$로의 ISC는 $M_S=0$ 성분에 대해 매우 효율적인 반면, 역방향 ISC($S_0 \to T_0$)는 $M_S=\pm 1$을 선호한다.
   • 들뜬 삼중항 상태로의 기하 구조 완화가 일어나면, 선택성이 변하여 $T_1 \to S_2$ ISC에서 $M_S=\pm 1$을 선호할 수 있다.
   • 결정적으로, 저자들은 스핀 선택성의 "맛(flavor)"(즉, $M_S=0$ 또는 $M_S=\pm 1$을 선호하는 성질)이 관련 상태들의 에너지 준위를 조절함으로써 제어될 수 있음을 발견했다. 이는 NV 센터와 유사하게 ODMR 대비(contrast) 및 스핀 초기화에 필요한 특정 ISC 역학을 달성하기 위한 경로를 제시한다.

5. 진동 결합 메커니즘:
   분석을 통해 진동 결합을 쿨롱 상호작용($J$)과 상대론적 상호작용($K$)의 기여로 분해하였다. $J$는 스핀 보존 전이를 지배하는 반면, $K$(SOC에 의해 구동됨)는 스핀 혼합 전이($\Delta M_S = \pm 1$)를 구동한다는 것을 발견했다. 연구는 스핀 선택성이 진동 모드의 대칭성과 본질적으로 연결되어 있음을 강조한다: 저에너지의 면외(out-of-plane) 모드는 $\Delta M_S = \pm 1$ 전이를 선호하는 반면, 고에너지의 면내(in-plane) 모드는 $\Delta M_S = 0$ 전이를 선호한다.

의의 및 주장
본 논문은 개방형 껍질 유기 시스템에서 스핀 결맞음 해제 및 스핀 선택적 들뜬 상태 과정을 이해하기 위한 통합된 ab initio 이론적 토대를 제공한다고 주장한다. 저자들은 트리앙굴렌 기반 디라디칼이 광학적으로 주소 지정 가능하고 긴 결맞이 시간을 가진 분자 큐비트를 개발하기 위한 화학적으로 튜닝 가능하며 유연한 구조라고 제안한다.

저자들은 현재의 유도체(3과 같은)가 저에너지 격자 및 치환기 포논으로 인해 감소된 결맞음을 겪지만, 이론적 프레임워크가 이러한 한계를 극克服하기 위한 구체적인 화학적 설계 전략을 식별해 냈음을 강조한다. 여기에는 유연한 치환기를 단단한 그룹으로 교체하거나, 저에너지 진동을 억제하기 위해 핵심 구조를 단단한 공유 결합 유기 골격(COF)에 통합하는 것이 포함된다. 또한, 연구는 ISC의 스핀 선택성이 상태 에너지를 튜닝함으로써 엔지니어링할 수 있는 $\pi$-공액 시스템의 일반적인 속성임을 시사하며, 이는 차세대 광학 활성 분자 양자 기술을 설계하기 위한 로드맵을 제공한다. 이 작업은 정적인 평형 고려를 넘어 동적인 다중 포논 경로를 사용하여 분자의 스핀 분해 ISC 비율을 완전히 ab initio 방식으로 정량적으로 시뮬레이션하는 단계로 나아가는 이정표가 된다.