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

Resumo Técnico: Dirradicais baseados em Trianguleno como Plataformas Quânticas Moleculares

Definição do Problema
O desenvolvimento da tecnologia quântica molecular requer a identificação de moléculas que possuam simultaneamente tempos de coerência de spin longos e interfaces spin-ópticas eficientes, idealmente operáveis à temperatura ambiente. Embora defeitos em estado sólido, como os centros de vacância de nitrogênio (NV) em diamante, ofereçam um parâmetro de referência para essas propriedades — exibindo inicialização de spin endereçável opticamente, leitura via ressonância magnética detectada opticamente (ODMR) e longos tempos de coerência — replicar essa funcionalidade em sistemas moleculares permanece um desafio significativo. Os atuais candidatos moleculares frequentemente lutam para combinar um estado fundamental de tripleto bem isolado, alta estabilidade química e uma interconversão entre sistemas (ISC) altamente seletiva ao spin, necessária para o controle óptico. Além disso, carecia-se de um arcabouço teórico unificado de primeira ordem (ab initio) que descrevesse tanto a decoerência de spin do estado fundamental quanto a dinâmica de spin seletiva do estado excitado em sistemas orgânicos de camada aberta.

Metodologia
Os autores empregam métodos avançados de primeira ordem (ab initio) para investigar a estrutura eletrônica e as propriedades dinâmicas do trianguleno (1), de seu derivado aza-cátion (2) e do cristal de 2,6,10-tri-tert-butil-4,8,12-trimesil-trianguleno (3).

• Estrutura Eletrônica: Simulações de química quântica multirreferencial (CASSCF com correções de correlação dinâmica NEVPT2) são utilizadas para determinar a escada de energia dos estados eletrônicos livres de spin. O estudo utiliza o solver ICE-CI para lidar com grandes espaços ativos (até 22 elétrons em 19 orbitais), reduzindo para um espaço ativo NEVPT2(10,10) para análise detalhada.
• Dinâmica Spin-Fônon: O relaxamento spin-rede ($T_1$) e o desfasamento puro ($T_2$) são modelados usando um arcabouço de acoplamento vibrônico relativístico. O Hamiltoniano inclui termos não relativísticos, acoplamento spin-órbita (SOC), acoplamento spin-spin (SSC) e interações Zeeman. O acoplamento vibrônico é tratado via expansão de Taylor do Hamiltoniano eletrônico em relação aos deslocamentos nucleares (termos lineares e quadráticos). As taxas de transição são calculadas usando equações mestras quânticas, contabilizando processos de um e dois fônons (absorção, emissão e combinação) dentro de uma aproximação de banho Markoviano.
• Desfasamento Spin-Spin: A decoerência causada por spins nucleares é simulada usando o método de Expansão de Correlação de Cluster (CCE) (via pyCCE), modelando o spin eletrônico central interagindo com um banho de spins nucleares de hidrogênio ou deutério.
• Dinâmica do Estado Excitado: As taxas de ISC são computadas via contribuições de dois fônons mediados por transições virtuais para estados intermediários, incorporando acoplamento vibrônico linear e quadrático.

Principais Contribuições e Resultados

1. Estrutura Eletrônica e Analogia com NV:
   O estudo confirma que o trianguleno (1) e seu derivado aza-cátion (2) possuem um estado fundamental de tripleto ($T_0$) bem separado do primeiro estado excitado singleto ($S_0$) por aproximadamente 0,5 eV. Essa separação de energia assemelha-se de perto à dos centros NV. O aza-cátion (2) exibe um gap tripleto-tripleto reduzido em comparação ao (1), com estados de tripleto excitados ($T_1, T_2$) posicionados de forma semelhante aos estados excitados do centro NV. A estrutura eletrônica do derivado cristalino (3) mostra-se semelhante a 1 e 2, embora o empacotamento cristalino reduza ligeiramente a simetria e os gaps de energia.

2. Relaxamento Spin-Fônon ($T_1$):
   Os cálculos preveem tempos de relaxamento spin-fônon excepcionalmente longos para as moléculas isoladas 1 e 2 à temperatura ambiente (300 K), com $T_1 \approx 1$ ms para 1 e $T_1 \approx 27$ ms para 2. O $T_1$ estendido em 2 é atribuído a uma força de acoplamento vibrônico reduzida. No entanto, a funcionalização e a inclusão em uma rede cristalina (composto 3) reduzem drasticamente o $T_1$ para $\sim 0,1$ $\mu$s a 300 K. Essa redução é impulsionada por fônons de baixa energia associados aos grupos substituintes volumosos (tert-butil e mesil) e interações de rede.

3. Desfasamento Spin-Spin ($T_2$):
   Simulações de decoerência induzida por spin nuclear revelam que, em um ambiente hidrogenado, o $T_2$ é limitado à escala de microssegundos (ex: $T_2 \approx 2,42$ $\mu$s para 3a-H). Contudo, o estudo demonstra que a engenharia química pode mitigar significativamente isso. Em um cenário hipotético onde a molécula central é deuterada e colocada em um ambiente livre de spin nuclear (ex: um solvente ou matriz deuterada), o $T_2$ aumenta em mais de uma ordem de magnitude para 0,21 ms a 10 K.

4. Interconversão entre Sistemas (ISC) Seletiva ao Spin:
   O cálculo das taxas de ISC de dois fônons a 300 K para a molécula 2 revela processos altamente seletivos ao spin. O estudo identifica que a seletividade de spin da ISC é sensível ao alinhamento de energia dos estados e à geometria molecular.

• Na geometria do estado fundamental, a ISC de $T_1$ para $S_2$ é altamente eficiente para o componente $M_S=0$, enquanto a ISC reversa ($S_0 \to T_0$) favorece $M_S=\pm 1$.
• Ao relaxar a geometria para o estado de tripleto excitado, a seletividade pode mudar, potencialmente favorecendo $M_S=\pm 1$ para a ISC $T_1 \to S_2$.
• Crucialmente, os autores descobrem que o "sabor" da seletividade de spin (favorecendo $M_S=0$ ou $M_S=\pm 1$) pode ser modulado pela engenharia dos níveis de energia dos estados envolvidos. Isso sugere um caminho para alcançar a dinâmica de ISC específica necessária para o contraste de ODMR e inicialização de spin, similar aos centros NV.

5. Mecanismos de Acoplamento Vibrônico:
   A análise decompõe o acoplamento vibrônico em contribuições de interações Coulombianas ($J$) e interações relativísticas ($K$). Observou-se que $J$ domina as transições que preservam o spin, enquanto $K$ (impulsionado por SOC) conduz as transições de mistura de spin ($\Delta M_S = \pm 1$). O estudo destaca que a seletividade de spin está intrinsecamente ligada à simetria dos modos vibracionais: modos fora do plano (baixa energia) favorecem transições $\Delta M_S = \pm 1$, enquanto modos no plano (alta energia) favorecem $\Delta M_S = 0$.

Significância e Alegações
O artigo afirma fornecer uma base teórica ab initio unificada para compreender a decoerência de spin e os processos de estado excitado seletivos ao spin em sistemas orgânicos de camada aberta. Ele postula que os dirradicais baseados em trianguleno representam uma arquitetura quimicamente ajustável e flexível para o desenvolvimento de qubits moleculares com endereçabilidade óptica e longos tempos de coerência.

Os autores enfatizam que, embora os derivados atuais (como o 3) sofram com a redução da coerência devido aos fônons de substituintes e de rede de baixa energia, o arcabouço teórico identifica estratégias específicas de design químico para superar essas limitações. Estas incluem substituir substituintes flexíveis por grupos rígidos ou integrar o núcleo em estruturas orgânicas covalentes (COFs) rígidas para suprimir as vibrações de baixa energia. Além disso, o estudo sugere que a seletividade de spin da ISC é uma propriedade geral de sistemas $\pi$-conjugados que pode ser engenheirada através do ajuste dos níveis de energia dos estados, oferecendo um roteiro para o design de tecnologias quânticas moleculares opticamente ativas de próxima geração. O trabalho marca um passo em direção à simulação quantitativa de taxas de ISC resolvidas por spin em moléculas usando abordagens totalmente ab initio, indo além das considerações de equilíbrio estático para caminhos dinâmicos de múltiplos fônons.