En espectroscopia de RMN , el formalismo del operador de producto es un método utilizado para determinar el resultado de secuencias de pulsos de una manera rigurosa pero sencilla. Con este método es posible predecir cómo evoluciona la magnetización en masa con el tiempo bajo la acción de pulsos aplicados en diferentes direcciones. Es una mejora neta con respecto al modelo vectorial semiclásico que no es capaz de predecir muchos de los resultados en espectroscopia de RMN y es una simplificación del formalismo de matriz de densidad completa.
En este modelo, para un solo espín , existen cuatro operadores base: , , y que representan respectivamente la polarización (diferencia de población entre los dos estados de espín), la coherencia cuántica simple (magnetización en el plano xy) y el operador unitario. Existen muchos otros operadores no clásicos para sistemas acoplados . Utilizando este enfoque, la evolución de la magnetización bajo precesión libre se representa mediante y corresponde a una rotación sobre el eje z con un ángulo de fase proporcional al desplazamiento químico del espín en cuestión:
Los pulsos en torno a los ejes x e y se pueden representar mediante y respectivamente; estos permiten interconvertir la magnetización entre planos y, en última instancia, observarla al final de una secuencia. Dado que cada espín evolucionará de manera diferente según su desplazamiento, con este formalismo es posible calcular exactamente dónde terminará la magnetización y, por lo tanto, idear secuencias de pulsos para medir la señal deseada excluyendo otras.
El formalismo del operador de producto es particularmente útil para describir experimentos en dos dimensiones como COSY, HSQC y HMBC.
Motivación para conjuntos de partículas de espín 1/2
A lo largo de esta sección, la constante de Planck se reduce para mayor comodidad.
El formalismo del operador producto se aplica habitualmente a conjuntos de partículas de espín 1/2 , ya que el hecho de que los operadores individuales satisfagan , donde es el operador identidad , hace que las relaciones de conmutación de los operadores producto sean particularmente simples. En principio, el formalismo podría extenderse a espines superiores, pero en la práctica se utiliza con más frecuencia el tratamiento general del tensor esférico irreducible . Por ello, a continuación solo consideramos el caso de espín 1/2.
La idea principal del formalismo es facilitar el seguimiento del operador de densidad del sistema , que evoluciona bajo un hamiltoniano según la ecuación de Liouville-von Neumann como
Para un hamiltoniano independiente del tiempo, el operador de densidad hereda sus soluciones del operador de evolución temporal de Schrödinger como
Dualidad operador-estado de densidad
Supongamos que un único espín-1/2 está en el estado , que es un estado propio del operador de espín z , es decir . De manera similar . Al hacer uso de la expansión de un operador hermítico en términos de proyecciones sobre sus mercados propios con valores propios como , el operador de densidad asociado es
donde es el operador de identidad. De manera similar, el operador de densidad para el estado es
Dado que los operadores de espín no tienen traza y el valor esperado de un operador para un sistema con operador de densidad es , los términos proporcionales al operador unitario no afectan las expectativas de los operadores de espín. Además, esas partes no evolucionan en el tiempo, ya que conmutan trivialmente con el hamiltoniano. Por lo tanto, esos términos pueden ignorarse, y el estado corresponde a un operador de densidad , mientras que el estado corresponde a un operador de densidad . Exactamente de la misma manera, la polarización a lo largo del eje x positivo, es decir, un estado , corresponde a un operador de densidad . Esta idea se extiende naturalmente a múltiples espines, donde los estados y operadores son productos directos de estados y operadores de espín único. Por lo tanto, los términos del operador en el operador de densidad tienen una dualidad directa con los estados.
En el caso de dos espines , los términos en el operador de densidad (ignorando la identidad por sí sola) pueden interpretarse como que representan
- - magnetización longitudinal
- - magnetización transversal en fase, que es la cantidad observable en RMN.
- - magnetización longitudinal antifase
- - orden longitudinal de dos espines
- - otras coherencias, que son más difíciles de interpretar, pero que pueden evolucionar hacia otros términos
donde eg es una abreviatura del producto Kronecker , donde es el operador de identidad en el giro y, de manera similar, es una abreviatura de .
Los factores de dos en los operadores de dos espines "verdaderos" sirven para permitir relaciones de conmutación convenientes en este caso específico de espín 1/2 (véase más abajo). Nótese también que podríamos optar por expandir el operador de densidad en la base, etc., donde los operadores transversales han sido reemplazados por operadores de elevación y descenso . Con la detección de cuadratura, el observable asociado con un espín individual es efectivamente el no hermítico , por lo que a veces esto es más conveniente.
Evolución del operador de densidad
Consideremos operadores que obedecen las relaciones de conmutación cíclica.
De hecho, solo las dos primeras relaciones son necesarias para la siguiente derivación, pero como normalmente trabajamos con operadores asociados a direcciones cartesianas, como los operadores de momento angular individuales, el tercer conmutador sigue un argumento de simetría.
Introduzca también el superoperador de conmutación de un operador (en nuestro caso, esto está relacionado más formalmente con la representación adjunta del álgebra de Lie cuyos elementos son ), que actúa como
En particular, para los operadores cíclicos, tenemos
y en consecuencia para enteros
Una identidad para dos operadores es
que se puede derivar poniendo donde es un parámetro escalar, diferenciando ambos lados con respecto a , y notando que ambos lados satisfacen la misma ecuación diferencial en ese parámetro, con la misma condición inicial en . En particular, para algún parámetro escalar , tenemos
donde la igualdad final se sigue de reconocer la serie de Taylor para el seno y el coseno. Ahora supongamos que el operador de densidad en el momento cero es , y se le permite evolucionar libremente bajo el hamiltoniano donde es un escalar. Usando los resultados anteriores, el operador de densidad en algún momento posterior estará dado por
La interpretación de esto es que, aunque el momento angular de espín nuclear en sí no está conectado a las rotaciones en el espacio tridimensional de la misma manera que lo está el momento angular, la evolución del operador de densidad puede verse como rotaciones en un espacio abstracto, en el que los operadores son los generadores de rotaciones sobre los ejes. Un ejemplo de un conjunto de generadores de este tipo son los propios operadores de espín.
Ahora también introducimos la "notación de flecha" que se utiliza típicamente en RMN, que escribe la evolución general dada anteriormente como una abreviatura.
- .
Con referencia más específica a los pulsos de radiofrecuencia aplicados durante los experimentos de RMN, un pulso duro con un ángulo de punta alrededor de una dirección se escribe como sobre la flecha y corresponde a tomar como generador de rotación en la Ecuación 1. Cuando no hay ambigüedad, la etiqueta de la flecha puede omitirse o, por ejemplo, puede ser texto en su lugar.
Obsérvese que un cálculo más complicado ahora se ha reducido a un procedimiento más simple que no requiere conocimiento de la mecánica cuántica subyacente, especialmente porque los subespacios de los operadores cíclicos se pueden tabular de antemano.
Ejemplos
El pulso de reenfoque de 180°
El hamiltoniano para un solo espín que evoluciona bajo un desplazamiento químico de frecuencia angular es
lo que significa que en un conjunto de muchos de estos espines con desplazamientos químicos ligeramente diferentes, hay un desfase de la magnetización en el plano - . Considere la secuencia de pulsos
- — — —
donde es un intervalo de tiempo. Partiendo de un estado de equilibrio con toda la polarización a lo largo del eje , la evolución de un espín individual en el conjunto es
Por tanto, esta secuencia reorienta la magnetización transversal producida por el primer pulso, independientemente del valor del desplazamiento químico.
Como indicación de la utilidad del formalismo, supongamos en cambio que tratamos de alcanzar el mismo resultado utilizando solo estados y, por lo tanto, los operadores de evolución temporal de Schrödinger. Esto equivale a intentar simplificar el propagador unitario que lleva el estado inicial al estado final como , donde explícitamente
En esencia, queremos encontrar el propagador en la forma , es decir, como un exponente único de una combinación de operadores, porque eso da el hamiltoniano efectivo que actúa durante la secuencia. Dado que los argumentos de los exponenciales en la forma original del propagador no conmutan, esto equivale a resolver un ejemplo específico del problema de Baker–Campbell–Hausdorff (BCH). En este caso relativamente simple, podemos resolver el problema de BCH utilizando el hecho de que para el operador unitario , el operador y la función , así como la similitud matemática de los operadores de espín con los generadores de rotación físicos, que nos permiten escribir
Por lo tanto , solo queda el efecto del pulso de 180°, que concuerda con el tratamiento del operador de producto. Para secuencias más grandes de pulsos, este tratamiento de estado rápidamente se vuelve aún más difícil de manejar, a menos que se utilicen métodos más avanzados como la teoría hamiltoniana exacta y efectiva (que proporciona expresiones de forma cerrada para los propagadores entrelazados a través del teorema de Cayley-Hamilton y las descomposiciones propias).
La amplitud de un eco de Hahn en un campo magnético no homogéneo
Como una extensión del pulso de reenfoque tratado anteriormente, considere un conjunto de dos pulsos con ángulos de giro arbitrarios y , es decir, una secuencia
- — — —
donde nuevamente es un intervalo de tiempo. Si se eliminan libremente los términos irrelevantes, la evolución para un solo giro con un desfase hasta justo después del segundo pulso es
Consideremos ahora un conjunto de espines en un campo magnético que es lo suficientemente heterogéneo como para desfasar por completo los espines en el intervalo entre los pulsos. Después del segundo pulso, podemos descomponer los términos restantes en una suma de dos poblaciones de espines que difieren solo en el signo del término, en el sentido de que para un espín individual tenemos
donde utilizamos las identidades y .
Son los espines de la nueva población que se han generado por el segundo pulso, es decir, el que tiene , los que conducirán a la formación de un eco después de la evolución para el siguiente intervalo. Por lo tanto, recordando incluir el introducido por el primer pulso, la amplitud del eco de Hahn resultante en relación con el producido por una secuencia de pulsos de reenfoque ideal de 90°—180° es aproximadamente
Tenga en cuenta que este no es un resultado exacto, ya que solo considera el reenfoque de la polarización que era transversal inmediatamente antes del segundo pulso. En realidad, habrá más componentes transversales originados por la inclinación de la magnetización longitudinal que permaneció después del primer pulso. Sin embargo, para muchos ángulos de inclinación, esta es una buena regla general.
Para llegar a este resultado utilizando el formalismo de estado, tendríamos que haber evaluado de manera no trivial el propagador de rotación como
y luego evaluar una probabilidad de transición considerando el resultado de aplicar esto a un estado que representa polarización en el plano transversal.
DEPT (Mejora sin distorsión por transferencia de polarización)
DEPT (Distortionless Enhancement by Polarisation Transfer) es una secuencia de pulsos que se utiliza para distinguir entre la multiplicidad de enlaces de hidrógeno al carbono, es decir, puede separar los grupos C, CH, CH 2 y CH 3 . Esto se logra aprovechando el acoplamiento carbono - hidrógeno heteronuclear y variando el ángulo de la punta del pulso final en la secuencia. La secuencia de pulsos básica se muestra a continuación.
Bajo el supuesto de acoplamiento débil, los términos de desplazamiento químico conmutan con el término de acoplamiento en el hamiltoniano. Por lo tanto, podemos ignorar el desplazamiento químico reenfocado (ver § El pulso de reenfoque de 180°) en los dos intervalos que contienen pulsos , a saber y , y además abstenernos de evaluar la evolución del desplazamiento químico en el último período . El tiempo de separación de los pulsos se ajusta a la fuerza de acoplamiento (con el coeficiente hamiltoniano asociado ) de modo que satisfaga
- ,
porque entonces el primer término en el operador de densidad evolucionado en la Ecuación 2 desaparece bajo la evolución de acoplamiento pura entre los pulsos.
es
Etiquete el espín del hidrógeno como , y el espín del carbono como . Para fines ilustrativos, suponemos que el estado de equilibrio solo tiene polarización en el espín - (en realidad, también habrá polarización en el espín, con las poblaciones relativas determinadas por los factores térmicos de Boltzmann ). El hamiltoniano de acoplamiento - es
lo que da la siguiente evolución
Los conmutadores no triviales utilizados para identificar el subespacio cíclico son
y en consecuencia la siguiente rotación cíclica
donde utilizamos la 'identidad de producto mixto' , que relaciona los productos de matriz y Kronecker para dimensiones compatibles de , y también el hecho de que dado que los dos valores propios de cualquiera de los operadores de espín 1/2 son , cualquiera de sus cuadrados está dado por el teorema de Cayley-Hamilton .
Obsérvese también que el término es invariante bajo la evolución del acoplamiento β. Es decir, el término conmuta con el hamiltoniano y, en este caso, esto se puede confirmar manualmente evaluando el conmutador utilizando las representaciones matriciales de los operadores de espín.
es2
Ahora etiquete los dos espines del hidrógeno como y el espín del carbono como . El hamiltoniano de acoplamiento es ahora
lo que da la siguiente evolución
donde 'otros' denota varios términos que pueden ignorarse con seguridad porque no evolucionarán hacia una polarización transversal observable en el espín objetivo . Los conmutadores cíclicos necesarios para tratar la evolución del acoplamiento son los siguientes tres conjuntos (y sus versiones, si es necesario)
es3
Un tratamiento similar (pero más extenso) da como resultado el término observable final .
APT (Prueba de protones unidos)
Consulte § DEPT (Mejora sin distorsión por transferencia de polarización) para conocer la notación utilizada en este ejemplo.
APT es similar a DEPT en el sentido de que detecta la multiplicidad de carbono. Sin embargo, tiene degeneraciones adicionales: proporciona señales positivas idénticas para C y CH 2 , y señales negativas idénticas para CH y CH 3 . A continuación se muestra una variación de la secuencia de pulsos básica.
La observación clave es que, dado que podemos ignorar nuevamente el desplazamiento químico reenfocado, la única dinámica relevante ocurre en el intervalo sin desacoplamiento de hidrógeno, donde podemos considerar únicamente el acoplamiento . Al usar un intervalo dos veces más largo que en el caso DEPT, nos aseguramos de que un operador de densidad de al comienzo del intervalo simplemente tenga su signo invertido después del acoplamiento (ya que esto corresponde a en el tratamiento general, y ). Los hamiltonianos para los acoplamientos a cada uno de los átomos de hidrógeno vecinos separados conmutan, por lo que el efecto general es multiplicarse por un factor . Esto motiva el signo alterno de la señal mencionada anteriormente.
Referencias
- James Keeler. "Understanding NMR Spectroscopy" (reimpreso en la Universidad de Cambridge ) . Universidad de California, Irvine . Consultado el 5 de agosto de 2012 .
- David Donne; David Gorenstein. "Una representación gráfica del formalismo del operador de producto" (PDF) . Universidad de Texas . Archivado desde el original (PDF) el 2013-05-14 . Consultado el 2012-08-05 .
- Hennig, Jürgen (1991). "Ecos: cómo generarlos, reconocerlos, utilizarlos o evitarlos en secuencias de imágenes por RM. Parte I: Propiedades fundamentales y no tan fundamentales de los ecos de espín". Conceptos en resonancia magnética . 3 (3): 125–143. doi :10.1002/cmr.1820030302.
- Mateescu, Gheorghe D; Valeriu, Adrián (1993). "Matriz de densidad de RMN 2D y tratamiento del operador del producto". Revista de Educación Química . 70 (6): A172. Código bibliográfico :1993JChEd..70S.172.. doi :10.1021/ed070pA172.3.
- Siminovitch, David; Untidt, Thomas; Nielsen, Niels Chr. (2004). "Teoría hamiltoniana exacta y efectiva. II. Expansión polinómica de funciones matriciales y operadores exponenciales unitarios entrelazados". The Journal of Chemical Physics . 120 (1): 51–66. Bibcode :2004JChPh.120...51S. doi :10.1063/1.1628216. PMID 15267261.
- Zhang, Yuning; Han, Fei; Jerschow, Alexej (2013). "Formalismo del operador de producto". EMag Res. doi : 10.1002/9780470034590.emrstm1310. ISBN 978-0470034590.