Capítulo 10 — Conexión con Mecánica Cuántica

La mecánica estadística clásica se construye sobre el espacio de fases y la distribución de probabilidad de Gibbs. En el mundo cuántico, la noción de trayectoria en el espacio de fases pierde sentido debido al principio de incertidumbre, y el concepto de probabilidad adquiere una dualidad profunda: la probabilidad intrínseca de la mecánica cuántica (superposición) se combina con la probabilidad estadística del desconocimiento (ensambles). El operador matriz densidad unifica ambas fuentes de incertidumbre en un formalismo elegante que generaliza la mecánica estadística al dominio cuántico y establece el lenguaje natural para la teoría de la información cuántica.

1. Operador matriz densidad

Consideremos un sistema cuántico del cual solo conocemos que se encuentra en el estado puro $|\psi_i\rangle$ con probabilidad clásica $p_i$, donde $p_i \geq 0$ y $\sum_i p_i = 1$. El operador matriz densidad (o simplemente operador densidad) se define como:

$$\hat{\rho} = \sum_i p_i \, |\psi_i\rangle \langle\psi_i|$$

El operador $\hat{\rho}$ actúa sobre el espacio de Hilbert $\mathcal{H}$ del sistema y satisface tres propiedades fundamentales que lo caracterizan como un estado físico legítimo:

Hermiticidad: $\hat{\rho}^\dagger = \hat{\rho}$. Esto garantiza que todos sus autovalores son reales y que los valores medios de observables son reales.
Traza unitaria: $\text{Tr}(\hat{\rho}) = \sum_i p_i = 1$. Es la condición de normalización: la probabilidad total de encontrar al sistema en algún estado es uno.
Semidefinido positivo: $\langle \phi | \hat{\rho} | \phi \rangle \geq 0$ para todo $|\phi\rangle \in \mathcal{H}$. Todos los autovalores de $\hat{\rho}$ son no negativos y representan probabilidades.

Un estado se dice puro si y solo si $\hat{\rho}^2 = \hat{\rho}$ (proyector de rango 1), equivalentemente $\text{Tr}(\hat{\rho}^2) = 1$. En ese caso, $\hat{\rho} = |\psi\rangle\langle\psi|$ para algún $|\psi\rangle$. Un estado es mezcla cuando $\text{Tr}(\hat{\rho}^2) < 1$, reflejando ignorancia clásica adicional a la incertidumbre cuántica intrínseca.

Valor medio de un observable: El valor de expectación de cualquier observable representado por un operador hermítico $\hat{O}$ se obtiene mediante la traza $\langle\hat{O}\rangle = \text{Tr}(\hat{\rho}\,\hat{O})$. Esta expresión contiene simultáneamente el promedio cuántico sobre cada estado puro y el promedio estadístico sobre el ensamble.

2. Ecuación de von Neumann

Así como la distribución de probabilidad clásica $\rho(q, p, t)$ evoluciona según la ecuación de Liouville, el operador densidad cuántico evoluciona según la ecuación de von Neumann. Partiendo de la ecuación de Schrödinger $i\hbar \frac{d}{dt}|\psi_i\rangle = \hat{H}|\psi_i\rangle$ para cada estado puro de la mezcla, se obtiene:

$$i\hbar \frac{\partial \hat{\rho}}{\partial t} = [\hat{H}, \hat{\rho}] = \hat{H}\hat{\rho} - \hat{\rho}\hat{H}$$

Esta ecuación es la contraparte cuántica exacta de la ecuación de Liouville clásica $\partial\rho/\partial t = \{H, \rho\}$, con la sustitución del corchete de Poisson por el conmutador (multiplicado por $i\hbar$):

$$\{A, B\} \;\longrightarrow\; \frac{1}{i\hbar} [\hat{A}, \hat{B}]$$

Una consecuencia inmediata es que si $\hat{\rho}$ depende solo del hamiltoniano, $[\hat{H}, \hat{\rho}] = 0$, entonces $\partial\hat{\rho}/\partial t = 0$ y el estado es estacionario. Los ensambles de equilibrio (microcanónico, canónico, gran canónico) son casos particulares de estados estacionarios.

La evolución temporal para un sistema aislado se escribe formalmente mediante el operador de evolución $\hat{U}(t) = e^{-i\hat{H}t/\hbar}$:

$$\hat{\rho}(t) = \hat{U}(t)\,\hat{\rho}(0)\,\hat{U}^\dagger(t) = e^{-i\hat{H}t/\hbar}\,\hat{\rho}(0)\,e^{i\hat{H}t/\hbar}$$

Reversibilidad: La ecuación de von Neumann, al igual que la de Liouville clásica, es reversible en el tiempo. La irreversibilidad macroscópica (crecimiento de la entropía) no está en la dinámica microscópica sino en el coarse-graining y en la elección de condiciones iniciales. Este es el corazón del problema de la flecha del tiempo en mecánica estadística.

3. Entropía de von Neumann

La entropía de von Neumann generaliza la entropía de Gibbs-Shannon al dominio cuántico y constituye la medida fundamental de desorden e información para un estado cuántico:

$$S(\hat{\rho}) = -k_B \, \text{Tr}\left(\hat{\rho} \ln \hat{\rho}\right)$$

En la base que diagonaliza $\hat{\rho}$, con autovalores $\lambda_i$, la entropía adopta la forma familiar de Shannon:

$$S(\hat{\rho}) = -k_B \sum_i \lambda_i \ln \lambda_i$$

Propiedades fundamentales de la entropía de von Neumann:

No negatividad: $S(\hat{\rho}) \geq 0$, con $S = 0$ si y solo si $\hat{\rho}$ es un estado puro.
Cota superior: Para un espacio de Hilbert de dimensión $D$, $S(\hat{\rho}) \leq k_B \ln D$, alcanzando el máximo en el estado máximamente mezclado $\hat{\rho} = \hat{I}/D$.
Concavidad: $S(\sum_i p_i \hat{\rho}_i) \geq \sum_i p_i S(\hat{\rho}_i)$, reflejando que mezclar incrementa la entropía.
Subaditividad: $S(\hat{\rho}_{AB}) \leq S(\hat{\rho}_A) + S(\hat{\rho}_B)$, con igualdad solo para estados producto $\hat{\rho}_{AB} = \hat{\rho}_A \otimes \hat{\rho}_B$.
Invariancia unitaria: $S(\hat{U}\hat{\rho}\,\hat{U}^\dagger) = S(\hat{\rho})$, la evolución unitaria aislada conserva la entropía.

Teorema H cuántico: La ecuación de von Neumann para un sistema aislado implica $dS/dt = 0$. Para obtener crecimiento de entropía ($dS/dt \geq 0$) es necesario introducir interacción con un entorno o realizar una traza parcial sobre grados de libertad inaccesibles (operador densidad reducido). Esta es la esencia de la decoherencia cuántica.

4. Ensamble canónico cuántico

En completo paralelismo con el ensamble canónico clásico, el ensamble canónico cuántico describe un sistema en equilibrio térmico con un reservorio a temperatura $T$. Maximizando la entropía de von Neumann sujeta a la condición de energía media fija $\langle\hat{H}\rangle = E$, el operador densidad resulta:

$$\hat{\rho} = \frac{e^{-\beta \hat{H}}}{Z}, \qquad Z = \text{Tr}\left(e^{-\beta \hat{H}}\right)$$

donde $\beta = 1/(k_B T)$ y $Z$ es la función de partición canónica cuántica. La traza puede evaluarse en cualquier base ortonormal completa; la elección más natural es la base de autoestados de energía $\hat{H}|n\rangle = E_n|n\rangle$:

$$Z = \sum_n e^{-\beta E_n}$$

La energía libre de Helmholtz retiene su relación fundamental:

$$F = -k_B T \ln Z, \qquad E = \langle\hat{H}\rangle = -\frac{\partial \ln Z}{\partial \beta}, \qquad S = \frac{E - F}{T}$$

En el límite de altas temperaturas ($\beta \to 0$), $e^{-\beta\hat{H}} \approx \hat{I} - \beta\hat{H}$ y el operador densidad tiende al estado máximamente mezclado $\hat{\rho} \approx \hat{I}/D$, recuperando el postulado de igual probabilidad a priori de la mecánica estadística clásica.

Conexión con la integral de camino: La traza $\text{Tr}(e^{-\beta\hat{H}})$ admite una representación como integral funcional (Feynman-Kac) en la que la temperatura inversa $\beta$ juega el papel de un tiempo imaginario $\tau = i\beta\hbar$. Esta analogía, fundamental en teoría cuántica de campos y materia condensada, unifica la mecánica estadística cuántica con la formulación de integral de camino de la mecánica cuántica.

5. Entrelazamiento y su medida

El entrelazamiento cuántico es la correlación no clásica más profunda que puede existir entre subsistemas. Un estado puro bipartito $|\psi\rangle \in \mathcal{H}_A \otimes \mathcal{H}_B$ se dice separable si puede escribirse como producto $|\psi\rangle = |\phi\rangle_A \otimes |\chi\rangle_B$; de lo contrario, está entrelazado. La descomposición de Schmidt proporciona una representación canónica:

$$|\psi\rangle = \sum_{i=1}^{r} \sqrt{\lambda_i} \, |i_A\rangle \otimes |i_B\rangle$$

donde $\lambda_i \geq 0$, $\sum_i \lambda_i = 1$, $r \leq \min(d_A, d_B)$ es el número de Schmidt, y $\{|i_A\rangle\}, \{|i_B\rangle\}$ son bases ortonormales en cada subsistema. Los coeficientes de Schmidt $\lambda_i$ son los autovalores de la matriz densidad reducida de cualquiera de los subsistemas:

$$\hat{\rho}_A = \text{Tr}_B(|\psi\rangle\langle\psi|) = \sum_i \lambda_i \, |i_A\rangle\langle i_A|$$

La entropía de entrelazamiento (o entropía de von Neumann del subsistema reducido) cuantifica el grado de entrelazamiento de un estado puro bipartito:

$$S_{\text{ent}}(|\psi\rangle) = S(\hat{\rho}_A) = -\sum_i \lambda_i \ln \lambda_i = S(\hat{\rho}_B)$$

Para un estado producto ($r = 1$), $S_{\text{ent}} = 0$. Para un estado máximamente entrelazado en dimensión $d$, $\lambda_i = 1/d$ para todo $i$, y $S_{\text{ent}} = \ln d$ (en unidades de $k_B$).

Monogamia del entrelazamiento: A diferencia de las correlaciones clásicas, el entrelazamiento cuántico es monógamo: si dos sistemas $A$ y $B$ están máximamente entrelazados, ninguno de ellos puede tener entrelazamiento con un tercer sistema $C$. Esta propiedad es fundamental para la seguridad de protocolos de criptografía cuántica como BB84.

6. Aplicación: información cuántica

La teoría de la información cuántica extiende la teoría clásica de Shannon al dominio cuántico, donde el bit clásico $\{0, 1\}$ es reemplazado por el qubit (quantum bit), un sistema de dos niveles cuyo estado general es una superposición $|\psi\rangle = \alpha|0\rangle + \beta|1\rangle$ con $|\alpha|^2 + |\beta|^2 = 1$. La esfera de Bloch parametriza geométricamente el espacio de estados puros de un qubit:

$$|\psi\rangle = \cos\frac{\theta}{2}\,|0\rangle + e^{i\phi}\sin\frac{\theta}{2}\,|1\rangle, \qquad \hat{\rho} = \frac{1}{2}\left(\hat{I} + \mathbf{r} \cdot \boldsymbol{\hat{\sigma}}\right)$$

donde $\boldsymbol{\hat{\sigma}} = (\hat{\sigma}_x, \hat{\sigma}_y, \hat{\sigma}_z)$ son las matrices de Pauli y $|\mathbf{r}| \leq 1$ es el vector de Bloch. Los estados puros corresponden a la superficie de la esfera ($|\mathbf{r}| = 1$); los estados mezcla al interior ($|\mathbf{r}| < 1$).

Las puertas lógicas cuánticas son operaciones unitarias sobre qubits: la puerta de Hadamard $\hat{H}|0\rangle = (|0\rangle + |1\rangle)/\sqrt{2}$ crea superposición, la puerta CNOT entrelaza dos qubits, y el conjunto $\{\hat{H}, \hat{T}, \text{CNOT}\}$ es universal para computación cuántica.

El operador densidad es el lenguaje natural para describir canales cuánticos ruidosos, representados por mapas completamente positivos que preservan la traza (operadores de Kraus):

$$\hat{\rho} \;\longrightarrow\; \mathcal{E}(\hat{\rho}) = \sum_k \hat{E}_k \,\hat{\rho}\, \hat{E}_k^\dagger, \qquad \sum_k \hat{E}_k^\dagger \hat{E}_k = \hat{I}$$

Conexión conceptual: La mecánica estadística y la información cuántica comparten un núcleo matemático común centrado en el operador densidad y la entropía de von Neumann. La segunda ley de la termodinámica, el teorema de no-clonación cuántica, los límites de Landauer sobre el costo energético del borrado de información y la corrección de errores cuánticos son facetas de una misma estructura profunda que conecta física, información y computación.

7. Puente con los cursos de Mecánica Cuántica y Álgebra Lineal

Este capítulo final establece deliberadamente conexiones con los conocimientos adquiridos en los cursos previos de la carrera. Los conceptos aquí presentados no son nuevos en su esencia matemática, sino una reorganización y extensión de herramientas ya familiares:

De Mecánica Cuántica: La matriz densidad es precisamente el objeto que permite describir estados que no son vectores de estado puros. La ecuación de von Neumann es la ecuación de Schrödinger para operadores. La entropía de entrelazamiento cuantifica el grado de correlación cuántica, generalizando la noción de estados entrelazados estudiados en ejemplos como el par EPR.
De Álgebra Lineal: El operador densidad es un operador lineal en el espacio de Hilbert con propiedades bien definidas (hermítico, traza, positividad). La diagonalización de $\hat{\rho}$ no es más que el teorema espectral. La descomposición de Schmidt es la descomposición en valores singulares (SVD) aplicada al estado bipartito. Las puertas lógicas cuánticas son matrices unitarias actuando sobre el producto tensorial de espacios.
De Mecánica Estadística: La función de partición canónica $Z = \text{Tr}(e^{-\beta\hat{H}})$ es el análogo cuántico de la integral sobre el espacio de fases $Z = \int e^{-\beta H(q,p)} dq dp$. La entropía de von Neumann reemplaza la entropía de Gibbs $S = -k_B \int \rho \ln \rho \, dq dp$.

La unificación final es elegante: en el formalismo de la matriz densidad, la distinción entre probabilidad clásica (ignorancia) y probabilidad cuántica (superposición) desaparece. Ambas quedan codificadas en un único objeto matemático —el operador $\hat{\rho}$— cuyas propiedades algebraicas contienen toda la información física accesible sobre el sistema.

Lecturas recomendadas para profundizar: Quantum Computation and Quantum Information de Nielsen & Chuang (Cambridge, 2010) desarrolla exhaustivamente la teoría de la información cuántica. Statistical Mechanics: Entropy, Order Parameters, and Complexity de Sethna (Oxford, 2021) explora las conexiones entre mecánica estadística, información y sistemas complejos.