Microwave-free vector magnetometry and crystal orientation determination with Nitrogen-Vacancy centers using Bayesian inference

Resumen Técnico: Magnetometría de Vectores sin Microondas y Determinación de la Orientación Cristalina con Centros Nitrógeno-Vacante mediante Inferencia Bayesiana

Planteamiento del Problema
Los centros nitrógeno-vacante (NV) en diamante se han establecido como sensores cuánticos de alta sensibilidad para magnetometría, termometría y detección de campos eléctricos. La magnetometría convencional de NV suele depender de la Resonancia Magnética Detectada Ópticamente (ODMR), la cual requiere campos de microondas (MW) para manipular los estados de espín. Aunque es robusta, la entrega de MW introduce limitaciones prácticas, incluyendo la perturbación ambiental, el calentamiento y la necesidad de conexiones galvánicas que complican la integración, particularmente bajo condiciones criogénicas.

Existen alternativas sin microondas, como aquellas que explotan el anticruzamiento de niveles en el estado fundamental o la mezcla de subniveles de espín inducida por campos magnéticos. Sin embargo, estos métodos a menudo requieren una alineación estricta del campo magnético externo con ejes cristalográficos específicos o instrumentación óptica de alta resolución para mapear patrones de fotoluminiscencia (PL). Un cuello de botella crítico en la magnetometría vectorial es la necesidad de conocer la orientación cristalográfica precisa de los centros NV. La calibración tradicional implica espectroscopía ODMR bajo campos controlados, mientras que los métodos de orientación sin MW suelen estar limitados por la precisión de la imagen y el modelado de dipolos. Trabajos recientes que utilizan la relajación cruzada NV-NV han demostrado sensibilidad, pero han estado limitados por protocolos que requieren conocimiento previo de la cristalografía y una alineación precisa del campo para asignar resonancias de manera inequívoca.

Metodología
Los autores proponen un marco general para la magnetometría vectorial sin microondas y la determinación de la orientación cristalina que aprovecha la inferencia bayesiana para analizar mapas de fotoluminiscencia (PL) en 2D. El núcleo del método reside en el fenómeno de relajación cruzada: cuando las transiciones de espín de orientaciones NV inequivalentes se vuelven degeneradas, los procesos de flip-flop se ven potenciados, lo que provoca una caída medible en la intensidad de la PL.

1. Modelo Físico: El sistema utiliza un campo de sesgo sintonizable ($B_{bias}$) aplicado a lo largo de un eje fijo (eje $z$ del laboratorio) y un campo externo estático ($B_{ext}$). La muestra de diamante se rota alrededor de este eje (o equivalentemente, el campo transversal se rota). El campo magnético total en el marco de la muestra se determina mediante el ángulo de rotación $\phi$ y la fuerza del sesgo.
2. Condiciones de Resonancia: La relajación cruzada ocurre en superficies planas específicas en el espacio del campo magnético definidas por relaciones lineales entre los componentes del campo (por ejemplo, $B^s_x = \pm B^s_y$). Estas condiciones se manifiestan como caídas agudas en la señal de PL. El modelo contabiliza las degeneraciones discretas inherentes a la simetría tetraédrica ($O$ grupo) de los NV, donde múltiples orientaciones físicamente distintas producen patrones de PL idénticos.
3. Modelo Analítico Directo (Forward Model): En lugar de diagonalizar repetidamente el Hamiltoniano, los autores derivan una expresión analítica de forma cerrada para la señal de PL esperada. Este modelo incorpora las condiciones de resonancia como caídas tipo Lorentziano, parametrizadas por el vector de campo externo, la matriz de orientación de la muestra, el ancho de línea y los pesos de contraste.
4. Inferencia Bayesiana: El problema inverso —extraer el vector de campo magnético y la orientación del cristal a partir del mapa de PL— se resuelve mediante estimación bayesiana.
   • Distribución Posterior: El método computa la distribución de probabilidad posterior completa $P(\theta|y)$, donde $\theta$ representa los parámetros desconocidos (componentes de campo, ángulos de orientación) y $y$ son los datos experimentales de PL.
   • Gestión de la Degeneración: A diferencia de los métodos de ajuste de picos que fuerzan una única solución, el enfoque bayesiano captura naturalmente las degeneraciones discretas como picos distintos en la distribución posterior, cuantificando la incertidumbre y la ambigüedad.
   • Reducción de Simetría: Para asegurar la eficiencia computacional, el espacio de parámetros se restringe a un dominio mínimo e inequivalente por simetría, explotando la simetría tetraédrica de los centros NV, lo que reduce el espacio de búsqueda de 24 configuraciones equivalentes a un subconjunto manejable.

Contribuciones Clave

• Marco de Alineación Libre: El trabajo introduce un protocolo que no requiere conocimiento previo de la orientación del cristal ni una alineación precisa del campo externo con los ejes cristalográficos. Tanto el vector de campo magnético como la orientación de los NV se extraen simultáneamente de los datos.
• Eficiencia Analítica: Al derivar un modelo analítico para las resonancias de relajación cruzada, los autores evitan el cuello de botella computacional de la diagonalización repetida del Hamiltoniano, permitiendo una inferencia bayesiana rápida y escalable.
• Robustez ante la Degeneración: El marco gestiona explícitamente las simetrías discretas del sistema NV, proporcionando distribuciones posteriores completas que revelan múltiples soluciones válidas (por ejemplo, orientaciones relacionadas por inversión) en lugar de descartarlas o seleccionar un ajuste óptimo arbitrario.
• Demostración Experimental: Los autores validan experimentalmente el método utilizando un microscopio de fluorescencia de campo ancho y un conjunto de diamantes. Reconstruyeron con éxito la orientación de un cristal de diamante con una alineación desconocida y posteriormente utilizaron esta orientación calibrada para reconstruir un vector de campo magnético externo desconocido.

Resultados

• Determinación de la Orientación: En un experimento de prueba de concepto, el equipo determinó la orientación de un cristal de diamante respecto al marco del laboratorio. El análisis bayesiano del mapa de PL experimental produjo dos regiones de alta probabilidad correspondientes a orientaciones relacionadas por inversión, consistentes con la simetría del sistema. Los parámetros reconstruidos ($\alpha, \beta, \zeta$) reprodujeron el patrón de PL experimental con alta fidelidad.
• Magnetometría Vectorial: Utilizando la orientación determinada previamente, los autores reconstruyeron un campo magnético externo desconocido. Las distribuciones posteriores marginales para los componentes del campo ($b_z, b_\perp, \phi_0$) mostraron picos bien localizados. Los valores reconstruidos fueron compatibles con mediciones independientes de ODMR.
• Análisis de Simetría: Las simulaciones demostraron que cuando el eje de rotación se alinea con direcciones cristalinas de alta simetría (por ejemplo, [100], [110], [111]), los mapas de PL exhiben simetrías de rotación y espejo que conducen a ambigüedades en la reconstrucción del campo. Sin embargo, al elegir un eje de rotación arbitrario (rompiendo la simetría), se logra una reconstrucción única.

Significado
El artículo establece una ruta general hacia magnetómetros de NV compactos y libres de alineación, adecuados para aplicaciones de detección práctica donde las técnicas de ODMR convencionales son impracticables debido a restricciones de calentamiento o integración. Al aprovechar las características de relajación cruzada mediante inferencia bayesiana, el método elimina la necesidad de una alineación de campo estricta y del conocimiento previo de la cristalografía, ofreciendo una solución robusta para la magnetometría vectorial y la determinación de la orientación en configuraciones de campo arbitrarias. El trabajo amplía el alcance de la detección basada en NV al proporcionar un marco que cuantifica naturalmente las incertidumbres y gestiona las degeneraciones discretas intrínsecas de la simetría de los NV.