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

Resumen Técnico: Dirradicales basados en Trianguleno como Plataformas Cuánticas Moleculares

Planteamiento del Problema
El desarrollo de la tecnología cuántica molecular requiere la identificación de moléculas que posean simultáneamente tiempos de coherencia de espín prolongados e interfaces espín-ópticas eficientes, idealmente operables a temperatura ambiente. Si bien los defectos de estado sólido, como los centros nitrógeno-vacante (NV) en el diamante, ofrecen un punto de referencia para estas propiedades —exhibiendo direccionamiento de espín, lectura mediante resonancia magnética detectada ópticamente (ODMR) y tiempos de coherencia prolongados—, replicar esta funcionalidad en sistemas moleculares sigue siendo un desafío significativo. Los candidatos moleculares actuales suelen tener dificultades para combinar un estado fundamental triplete bien aislado, una alta estabilidad química y un cruce entre sistemas (ISC) altamente selectivo al espín, necesario para el control óptico. Además, se ha carecido de un marco teórico unificado de primeros principios que describa tanto la decoherencia de espín del estado fundamental como la dinámica de espín selectiva del estado excitado en sistemas orgánicos de capa abierta.

Metodología
Los autores emplean métodos avanzados de química cuántica ab initio de primeros principios para investigar la estructura electrónica y las propiedades dinámicas del trianguleno (1), su derivado de azacatión (2) y el cristal de 2,6,10-tri-terc-butil-4,8,12-trimesil-trianguleno (3).

• Estructura Electrónica: Se utilizan simulaciones de química cuántica multirreferencial (CASSCF con correcciones de correlación dinámica NEVPT2) para determinar la escala de energía de los estados electrónicos libres de espín. El estudio utiliza el resolvedor ICE-CI para manejar grandes espacios activos (hasta 22 electrones en 19 orbitales), reduciéndose a un espacio activo de NEVPT2(10,10) para un análisis detallado.
• Dinámica Espín-Fonón: La relajación espín-red ($T_1$) y el desfase puro ($T_2$) se modelan mediante un marco de acoplamiento vibrónico relativista. El Hamiltoniano incluye términos no relativistas, acoplamiento espín-órbita (SOC), acoplamiento espín-espín (SSC) e interacciones Zeeman. El acoplamiento vibrónico se trata mediante la expansión de Taylor del Hamiltoniano electrónico con respecto a los desplazamientos nucleares (términos lineales y cuadráticos). Las tasas de transición se calculan utilizando ecuaciones maestras cuánticas, contabilizando procesos de uno y dos fonones (absorción, emisión y combinación) dentro de una aproximación de baño de Markov.
• Desfase Espín-Espín: La decoherencia debida a los espines nucleares circundantes se simula utilizando el método de Expansión de Correlación de Clúster (CCE) (vía pyCCE), modelando el espín electrónico central interactuando con un baño de espines nucleares de hidrógeno o deuterio.
• Dinámica del Estado Excitado: Las tasas de ISC se computan mediante contribuciones de dos fonones mediados por transiciones virtuales a estados electrónicos intermedios, incorporando acoplamiento vibrónico tanto lineal como cuadrático.

Contribuciones Clave y Resultados

1. Estructura Electrónica y Analogía con NV:
   El estudio confirma que el trianguleno (1) y su derivado de azacatión (2) poseen un estado fundamental triplete ($T_0$) bien separado del primer estado excitado singlete ($S_0$) por aproximadamente quite 0.5 eV. Esta separación de energía se asemeja estrechamente a la de los centros NV. El azacatión (2) exhibe una brecha de energía triplete-triplete reducida en comparación con 1, con estados triplete excitados ($T_1, T_2$) posicionados de manera similar a los estados excitados de los centros NV. Se muestra que la estructura electrónica del derivado cristalino (3) es similar a la de 1 y 2, aunque el empaquetamiento cristalino reduce ligeramente la simetría y las brechas de energía.

2. Relajación Espín-Fonón ($T_1$):
   Los cálculos predicen tiempos de relajación espín-fonón excepcionalmente largos para las moléculas aisladas 1 y 2 a temperatura ambiente (300 K), con $T_1 \approx 1$ ms para 1 y $T_1 \approx 27$ ms para 2. El $T_1$ extendido en 2 se atribuye a un acoplamiento vibrónico reducido. Sin embargo, la funcionalización y la inclusión en una red cristalina (compuesto 3) reducen drásticamente el $T_1$ a $\sim 0.1$ $\mu$s a 300 K. Esta reducción es impulsada por fonones de baja energía asociados con los grupos sustituyentes voluminosos (terc-butilo y mesitilo) y las interacciones de la red.

3. Desfase Espín-Espín ($T_2$):
   Las simulaciones de decoherencia inducida por espín nuclear revelan que, en un entorno hidrogenado, el $T_2$ se limita a la escala de microsegundos (por ejemplo, $T_2 \approx 2.42$ $\mu$s para 3a-H). Sin embargo, el estudio demuestra que la ingeniería química puede mitigar esto significativamente. En un escenario hipotético donde la molécula central está deuterada y colocada en un entorno libre de espín nuclear (por ejemplo, un solvente o matriz deuterada), el $T_2$ aumenta en más de un orden de magnitud a 0.21 ms a 10 K.

4. Cruce entre Sistemas Selectivo al Espín (ISC):
   El cálculo de las tasas de ISC de dos fonones a 300 K para la molécula 2 revela procesos altamente selectivos al espín. El estudio identifica que la selectividad de espín de la ISC es sensible al alineamiento de energía de los estados y la geometría molecular.

• En la geometría del estado fundamental, la ISC desde $T_1$ hacia $S_2$ es altamente eficiente para el componente $M_S=0$, mientras que la ISC inversa ($S_0 \to T_0$) favorece $M_S=\pm 1$.
• Tras la relajación geométrica al estado triplete excitado, la selectividad puede cambiar, favoreciendo potencialmente $M_S=\pm 1$ para la ISC $T_1 \to S_2$.
• Crucialmente, los autores encuentran que el "sabor" de la selectividad de espín (favorecer $M_S=0$ o $M_S=\pm 1$) puede modularse mediante la ingeniería de los niveles de energía de los estados involucrados. Esto sugiere una vía para lograr la dinámica de ISC específica requerida para el contraste de ODMR y la inicialización de espín, de manera similar a los centros NV.

5. Mecanismos de Acoplamiento Vibrónico:
   El análisis descompone el acoplamiento vibrónico en contribuciones de interacciones de Coulomb ($J$) e interacciones relativistas ($K$). Se encuentra que $J$ domina las transiciones que preservan el espín, mientras que $K$ (impulsada por SOC) impulsa las transiciones de mezcla de espín ($\Delta M_S = \pm 1$). El estudio destaca que la selectividad de espín está intrínsecamente ligada a la simetría de los modos vibracionales: los modos fuera del plano (baja energía) favorecen las transiciones $\Delta M_S = \pm 1$, mientras que los modos en el plano (alta energía) favorecen $\Delta M_S = 0$.

Significancia y Reivindicaciones
El artículo afirma proporcionar una base teórica ab initio unificada para comprender la decoherencia de espín y los procesos de estado excitado selectivos al espín en sistemas orgánicos de capa abierta. Postula que los dirradicales basados en trianguleno representan una arquitectura químicamente sintonizable y flexible para desarrollar qubits moleculares con direccionabilidad óptica y tiempos de coherencia prolongados.

Los autores enfatizan que, si bien los derivados actuales (como 3) sufren de una coherencia reducida debido a los fonones de la red y de los sustituyentes de baja energía, el marco teórico identifica estrategias de diseño químico específicas para superar estas limitaciones. Estas incluyen reemplazar los sustituyentes flexibles por grupos rígidos o integrar el núcleo en marcos orgánicos covalentes (COFs) rígidos para suprimir las vibraciones de baja energía. Además, el estudio sugiere que la selectividad de espín de la ISC es una propiedad general de los sistemas $\pi$-conjugados que puede diseñarse mediante el ajuste de las energías de los estados, ofreciendo una hoja de ruta para diseñar la próxima generación de tecnologías cuánticas ópticamente activas. El trabajo marca un paso hacia la simulación cuantitativa de tasas de ISC resueltas por espín en moléculas utilizando enfoques totalmente ab initio, yendo más allá de las consideraciones de equilibrio estático hacia vías dinámicas de múltiples fonones.