Matriz S

En mecánica cuántica, la matriz S o matriz de scattering es un tipo de formalismo usado para calcular el resultado de un problema de scattering (dispersión) de partículas cuánticas interactuantes. Matemáticamente viene dado por un operador S entre dos espacios de Hilbert, cuando el conjunto de estados cuánticos admisibles inicial y final es finito entonces el operador S se reduce a una matriz y de ahí su nombre (ya que en sus inicios fue concebido como una matriz entre un conjunto posible de estados).

En el formalismo habitual, muchos problemas de dispersión con partículas subatómicas que se influyen mutuamente mediante interacción electromagnética, interacción fuerte o interacción débil no pueden ser resueltos de manera exacta. El formalismo de la matriz S permite realizar los cálculos numéricos para muchos casos que no admiten un tratamiento exacto. La idea básica del método consiste en suponer que los estados inicial y final de un sistema de partículas interectuantes son autoestados del hamiltoniano libre (sin interacción). El estado inicial se considera un estado del pasado remoto (que físicamente se concibe como el estado de las partículas cuando están muy lejos entre sí, antes de empezar a interacutar, y por tanto son autoestados "libres"), mientras que el estado final se considera también como un estado libre de interacción en el futuro remoto, como el que alcanzarán las partículas cuando se hayan separado definitivamente y no se ejerzan influencias mutuas. Estas consideraciones permiten construir los espacios de Hilbert de las partículas entrantes y salientes, a partir del espacio de Hilbert de una partícula libre aislada.

En la física de partículas de altas energías, resulta importante calcular las probabilidades de los diferentes posibles resultados (estados finales) de diversos experimentos de dispersión. Estos experimentos pueden dividirse en tres etapas:

Para que el formalismo de la matriz S sea aplicable es necesario poseer una teoría física del proceso de interacción que permita calcular las probabilidades de los diferentes resultados posibles. El proceso de colisión no es enteramente determinista, además el producto o estado final depende de la energía entrante. La energía inicial afecta críticamente a cual será rl resultado más probable, es decir, según la energía de las partículas entrantes el resultado más probable puede ser uno u otro. Cuando las condiciones de densidad de partículas son las suficientemente bajas se considera que puede usarse con bastante buena aproximación el formalismo de la matriz S, para aproximar la solución exacta de la teoría cuántica de campos para el problema de dispersión.

La matriz S esá estrechamente relacionada con las amplitudes de probabilidad de transición de la mecánica cuántica y con las secciones eficaces de las diversas interacciones. Las componentes o entradas numéricas de la matriz S se conocen como amplitudes de scattering. La amplitud de probabilidad de observar un proceso de scattering con un estado inicial ${displaystyle |a angle }$ y con un estado final ${displaystyle |b angle }$ viene dada por definición por:

${displaystyle mathrm {Prob} (a o b)equiv ,_{out}langle b|a angle _{in}=,_{out}langle b|S|a angle _{out}=,_{in}langle b|S|a angle _{in}equiv S_{ba}}$

Cada elemento de matriz S representa, por tanto, una amplitud de probabilidad de un proceso físico. Donde debe señalarse que los estados a y b son estados idealizados o asintóticos definidos por la ausencia de interacción entre ambos (debido a que ambos describan partículas alejadas de la localización donde se produce la interacción).

La matriz S así definida depende de la energía de la colisión, si se extiende al plano complejo la función que da la matriz S en términos de la energía, resulta que los polos de dicha función pueden identificarse con los estados ligados los estados virtuales o las resonancias. Los punto de ramifición de la matriz S en el plano complejo están asociados a la abertura de un canal de scattering.

En el enfoque hamiltoniano de la teoría cuántica de campos, la matriz S puede calcularse como la exponencial temporalmente ordenada de la integral del hamiltoniano en la "imagen" de interacción. El cálculo de dicha exponencial puede expresarse también como integral de camino. En cualquiera de esas dos representaciones puede usarse un cálculo perturbativo de la matriz S mediante diagramas de Feynman.

En teoría de la dispersión, la matriz S es un operador que aplica al estado de las partículas entrantes, el estado de una partícula saliente en la "imagen" de Heisenberg. Esto resulta muy útil ya que frecuentemente no es posible describir exactamente todos los detalles de la interacción (o al menos algunos de los más interesantes).

Formalmente la matriz S se define como un operador unitario entre los espacios de Hilbert.

${displaystyle S:{mathcal {H}}_{in} o {mathcal {H}}_{out}}$

El primer espacio de Hilbert contiene vectores que representan los estados asintóticamente posibles de una partícula o conjunto de partículas "entrante" y el segundo representa al conjunto de estados asintóticos de las partículas "salientes". La idea básica es representar el efecto de una colisión o interacción compleja de partículas observando las partículas entrantes y observando las partículas salientes. La matriz S puede predecir la probabilidad de cada posible estado final en función del estado inicial, de hecho el estado inicial y final están relacionados por:

${displaystyle |psi angle _{in}equiv S|psi angle _{out};;quad |psi angle _{out}equiv S^{dagger }|psi angle _{in}}$

Técnicamente la matriz se puede definir para cualquier espacio-tiempo asintóticamente soluble sin horizonte, aunque comúnmente se plantea sobre un espacio-tiempo plano de Minkowski, ya que matemáticamente es la situación más simple. Para este caso particularmente simple, el espacio de Hilbert entrante y saliente es el espacio donde actúa una representación irreducible unitaria del grupo de Lorentz inhomogéneo. Y además la matriz S admite una representación como producto de operadores unitarios de evolución.

Usando la notación de Dirac, el estado del vacío cuántico se representa sencillamente como ${displaystyle |0 angle }$. Si ${displaystyle a^{dagger }(k)}$ es un operador de creación, su conjugado hermítico, es el operador de aniquilación o destrucción, cuya acción sobre el vector asociado al estado de vacío cuántico es:

${displaystyle a(k)left|0 ight angle =0.}$

En estas condiciones se definen dos tipos de operadores de creación y destrucción que actúan sobre diferentes espacios de Hilbert (el espacio de estados inicial i, y el espacio de estados final f), ${displaystyle a_{i}^{dagger }(k)}$ and ${displaystyle a_{f}^{dagger }(k)}$. De tal manera que:

Es posible asumir que ${displaystyle left|I,0 ight angle }$ y ${displaystyle left|F,0 ight angle }$ son ambos invariantes bajo traslaciones y que los estados ${displaystyle left|I,k_{1}ldots k_{n} ight angle }$ y ${displaystyle left|F,p_{1}ldots p_{n} ight angle }$ son autoestados del operador de momento lineal ${displaystyle {mathcal {P}}^{mu }}$, "conectando" y "desconectando" adiabáticamente la interacción.

En la "imagen" de Heisenberg los estados son independientes del tiempo, de tal manera que los estados iniciales pueden expresarse como combinación lineal de una base para los estados finales y viceversa, tal como sigue:

donde ${displaystyle left|C_{m} ight|^{2}}$ es la probabilidad de que la interacción transforme ${displaystyle left|I,k_{1}ldots k_{n} ight angle }$ en ${displaystyle left|F,p_{1}ldots p_{m} ight angle }$

De acuerdo con el teorema de Wigner, ${displaystyle S}$ debe ser un operador unitario tal que ${displaystyle leftlangle I,eta ight|Sleft|I,alpha ight angle =S_{alpha eta }=leftlangle F,eta |I,alpha ight angle }$. Es más, ${displaystyle S}$ deja el estado de vacío cuántico invariante y las transformaciones de campos del espacio de estados iniciales en estados del espacio final:

${displaystyle Sleft|0 ight angle =left|0 ight angle ,qquad phi _{f}=S^{-1}phi _{i}S}$

Si la S describe la interacción correctamente debe cumplirse que dos estados inicial y final, para los cuales su amplitud de probabilidad no sea nula deban tener la misma energía y momento (y eventualmente otras leyes de conservación también deben cumplirse).

Usando la imagen de Heisenberg para relacionarlos operadores relevantes en la interacción, se tiene que:

Por tanto,

porque

Substituyendo la expresión explícita para U se obtiene:

Esta ecuación sin embargo no es explícitamente covariante.

La matriz S fue introducida por primera vez por John Archibald Wheeler en un artículo de 1937 titulado "'On the Mathematical Description of Light Nuclei by the Method of Resonating Group Structure'".^[1]​ En ese artículo Wheeler introdujo una matriz de scattering que era una matriz unitaria de coeficientes que conectaban "el comportamiento asintótico de una solución particular arbitraria [de un sistema de ecuaciones integrales] con las soluciones de una forma estándar".^[2]​

En la década de 1940 Werner Heisenberg desarrolló, independientemente, la idea de la matriz S. Debido a las divergencias que plagaban la teoría cuántica de campos tal como se usaba en ese tiempo, Heisenberg trató de aislar las características esenciales de la teoría que pudieran no verse afectadas por cambios futuros en la teoría. Al hacer eso introdujo una matriz S "característica" unitaria.^[2]​

El hamiltoniano cuántico que describe un sistema puede ser dividido en una parte que representa la evolución libre sin interacción (hamiltoniano libre o no perturbado) ${displaystyle scriptstyle {hat {H}}^{0}}$, y que por tanto no contiene términos de interacción y una segunda parte que incluye todos los términos de interacción (hamiltoniano de intercción o perturbación) ${displaystyle scriptstyle {hat {H}}^{I}}$.

En la imagen de evolución temporal de Schrödinger es válida la ecuación de Schrödinger:

${displaystyle ihbar {frac {partial left|psi _{S}(t) ight angle }{partial t}}=H_{S}left|psi _{S}(t) ight angle }$

Definiendo la llamada imagen de interacción como:

${displaystyle left|psi _{I}(t) ight angle =e^{{frac {i}{hbar }}H_{S}^{0}t}left|psi _{S}(t) ight angle }$

de aquí puede escribirse:

es decir, en la imagen de interacción, los estados evolucionan de acuerdo con una dinámica dada excplusivamente por el hamiltoniano de interacción.

Para obtener la matriz S, se definen sus elementos como:

${displaystyle S_{f,i}=leftlangle psi _{f}|psi (t=+infty ) ight angle =leftlangle psi _{f}|S|psi _{i} ight angle }$

Para calcular explícitamente ${displaystyle S_{f,i}}$ se reescribe en forma integral la ecuación de Schrödinger para la función de donda en la imagen de interacción:

${displaystyle left|psi _{I}(t) ight angle =left|psi _{I}(t=-infty ) ight angle -ihbar int _{-infty }^{t}d au H_{S}^{I}left( au ight)left|psi _{I}( au ) ight angle }$

Dado que ${displaystyle left|psi _{I}( au ) ight angle }$ satisface una ecuación análoga es posible iterar una y otra vez esta ecuación hasta obtener un desarrollo perturbativo:

Y a partir de ese desarrollo se obtiene:

${displaystyle S=1+sum _{n=1}^{infty }int _{-infty }^{infty }dt_{1}int _{-infty }^{t_{1}}dt_{2}int _{-infty }^{t_{2}}dt_{3}ldots int _{-infty }^{t_{n-1}}dt_{n}left(H_{S}^{I}(t_{1})cdot H_{S}^{I}(t_{2})cdot H_{S}^{I}(t_{3})cdot ldots cdot H_{S}^{I}(t_{n}) ight)}$. ${displaystyle =1+sum _{n=1}^{infty }{frac {1}{n!}}int _{-infty }^{infty }dt_{1}int _{-infty }^{+infty }dt_{2}int _{-infty }^{+infty }dt_{3}ldots int _{-infty }^{+infty }dt_{n}Tleft(H_{S}^{I}(t_{1})cdot H_{S}^{I}(t_{2})cdot H_{S}^{I}(t_{3})cdot ldots cdot H_{S}^{I}(t_{n}) ight)}$,

donde ${displaystyle T}$ indica el producto cronológicamente ordenado de los operadores entre paréntesis:

${displaystyle Tphi (x)psi (y)= heta (x^{0}-y^{0})left(phi (x)psi (y) ight)+ heta (y^{0}-x^{0})left(psi (y)phi (x) ight)}$

donde ${displaystyle heta left(x^{0}-y^{0} ight)}$ es la función escalón unitario de Heaviside.

El desarrollo anterior es precisamente la serie de Dyson para la matriz ${displaystyle S}$. Y por tanto dicha matriz puede ser calculada hasta cualquier orden de aproximación mediante la serie anterior, usando el teorema de Wick.

De manera alternativa puede usarse la fórmula de reducción de LSZ para calcular los elementos de la matriz S. Este procedimiento es más sofisticado y es preferible en teoría cuántica de campos, no hace uso de la serie de Dyson, aunque sí usa la funciones de Green proporcionandas por la formulación mediante integrales funcionales.

Para ilustrar este tipo de cálculo se considera un caso particular, el de un campo escalar ${displaystyle phi }$ con masa asociada m y con una acción dada por_

${displaystyle {mathcal {S}}[phi ]=int { ext{d}}^{4}x{frac {1}{2}}left(partial _{mu }phi partial ^{mu }phi -m^{2}phi ^{2} ight)-V(phi )}$

donde ${displaystyle V(phi )}$ puede ser, por ejemplo, un término de interacción del tipo ${displaystyle lambda phi ^{4}}$, que por el momento no es necesario especificar.

Las funciones de Green de n puntos están definidas como los valores esperados sobre el vacío del producto cronológico ordenado de n campos:

${displaystyle G^{(n)}(x_{1},x_{2},dots ,x_{n})equiv langle 0|Tleft(phi (x_{1})phi (x_{2})dots phi (x_{n}) ight)|0 angle ={frac {int {mathcal {D}}phi ,phi (x_{1})phi (x_{2})dots phi (x_{n}),e^{i{mathcal {S}}}}{int {mathcal {D}}phi ,e^{i{mathcal {S}}}}}}$

Estas funciones son calculables perturbativamente a través del citado teorema de Wick. Puede demostrarse que la transformada de Fourier de las funciones de Green tienen polos que se corresponden con las masas físicas de las partículas, debido a que cuando ${displaystyle p_{mu }p^{mu }=m^{2}}$ dichas funciones tienen un polo. A estos polos le corresponden propiamente los estados asintóticos de la teoría: de hecho, estos estados son creados o destruidos por los campos "in" y "out", que satisfacen la ecuación de Klein-Gordon

que difieren de las ecuaciones de movimiento correctas solo por la ausencia del potencial de interacción. en consecuencia, y de modo intutivo, es necesario extraer la contribución de los polos de las funciones de Green para obtener las funcioone de Green construidas con los campos asintóticos, que generan propiamente los elementos de la matriz S deseados. Si en el estado inicial están presentes m partículas con momentos lineales q₁,...,q_m y en el estado final hay n partículas con momentos p₁,...,p_n, la fórmula que describe el procedimiento viene dada por:

El proceso de extracción del polo es más evidente si la fórmula se escrib en términos de la transformada de Fourier de la función de Green. Aparte de la multiplicación por algunas constantes (entre las cuales están la constante del normalización de los campos Z) la fórmula muestra que basta multiplicar las funciones de Green por los factores ${displaystyle p^{2}-m^{2}}$, que eliminan los polos, y después hacer el límite on-shell de los momentos lineales, es decir, ${displaystyle p^{2} o m^{2}}$ correspondiente al las partículas físicas:

---