Average density of Bloch electrons in a homogeneous magnetic field: A second-order response

Resumen Técnico: Densidad Promedio de Electrones de Bloch en un Campo Magnético Homogéneo

Planteamiento del Problema
Este trabajo investiga la modificación de la densidad electrónica promedio en un sólido cristalino tridimensional (3D) limpio (ya sea un metal o un aislante) sometido a un campo magnético homogéneo, estático y débil. Mientras que las respuestas de campo fuerte, como el efecto Hall cuántico, están bien estudiadas, el régimen de campo débil —específicamente la respuesta no lineal (de segundo orden) de la densidad de carga— ha recibido comparativamente menos atención. Los autores pretenden calcular la expansión de la densidad promedio hasta el segundo orden en el campo magnético ($B$) para un sistema de múltiples bandas de simetría arbitraria, abordando tanto los regímenes lineales como los no lineales.

Metodología
Los autores desarrollan un marco teórico basado en electrones no interactuantes en una red 3D con un Hamiltoniano de partícula única $h = \frac{1}{2m}(\mathbf{p} - e\mathbf{A})^2 + U$. Las características metodológicas clave incluyen:

• Invariancia de Calibre (Gauge Invariance): El cálculo es explícitamente invariante de calibre, lo que asegura que los observables físicos sean independientes de la elección de fase de los estados de Bloch ($|nk\rangle \to e^{i\alpha(\mathbf{k})}|nk\rangle$).
• Regularización del Operador de Posición: Para manejar la naturaleza singular del operador de posición en la base de Bloch, el campo magnético se introduce inicialmente con una modulación espacial mediante un vector de onda $\mathbf{g} = q\hat{x}$. El límite de campo homogéneo ($q \to 0$) se toma solo al final del cálculo.
• Expansión de la Densidad: El operador de densidad se expande en modos de Fourier espaciales. La densidad promedio se deriva sumando estos modos en el límite de modulación evanescente.
• Teoría de Perturbaciones: La matriz de densidad se expande a primer y segundo orden en el campo magnético. Los autores tratan los procesos intrabanda e interbanda en igualdad de condiciones e incluyen procesos de relajación (tomando el límite $\tau \to \infty$ para sistemas limpios).
• Magnitudes Geométricas: El formalismo incorpora naturalmente las propiedades geométricas de las funciones de onda de Bloch, específicamente la curvatura de Berry y el tensor métrico cuántico de Fubini-Study.

Resultados Clave

1. Respuesta Lineal (Primer Orden):

   • Para aislantes (donde el potencial químico se encuentra dentro de una brecha de energía), la respuesta de densidad lineal recupera la conocida fórmula de Streda: $\frac{\partial n}{\partial B} = -\frac{1}{e}\sigma_H$, que relaciona el cambio de densidad con la conductividad Hall anómala ($\sigma_H$).
   • Para metales (donde el potencial químico cruza una banda de energía), los autores derivan una corrección a la fórmula de Streda. La respuesta lineal incluye un término adicional proporcional a los momentos magnéticos orbitales ($m_n$) en la superficie de Fermi:
     $$\frac{\partial n}{\partial B} = -\frac{1}{e}\sigma_H + \frac{1}{V} \sum_{n\mathbf{k}} f'_n m_n$$
     Aquí, $f'_n$ es la derivada de la función de Fermi, localizando la contribución en la superficie de Fermi.

2. Respuesta No Lineal (Segundo Orden):

   • La respuesta de densidad de segundo orden ($\varsigma^{(2)}$) se deriva como una expresión general que involucra funciones de onda de Bloch. Esta expresión es compleja, comprendiendo contribuciones de estados de la superficie de Fermi (FS) y de no-FS, así como procesos de una banda, dos bandas y tres bandas.
   • Momentos Magnéticos Geométricos: Un hallazgo significativo es la identificación de una contribución puramente geométrica a la respuesta de segundo orden. El tensor métrico cuántico de Fubini-Study ($g_{ab}^n$), que caracteriza la distancia entre estados de Bloch en el plano proyectivo complejo, induce un "momento magnético cuadrático". Este momento surge de la rotación de las funciones de onda de Bloch en el espacio de parámetros más que del momento angular de espín u orbital.
   • Para procesos intrabanda en un modelo simple de dos bandas, el tensor métrico cuántico proporciona la contribución dominante a la respuesta de densidad de segundo orden.

3. Implicaciones Termodinámicas (Efectos Magnetovolumétricos):

   • Los autores argumentan que un cambio en la densidad electrónica promedio requiere un cambio en el volumen del cristal o en la presión para mantener la periodicidad.
   • Derivan expresiones que vinculan los cambios fraccionales en el volumen ($V$) y la presión ($P$) con la respuesta de densidad. En presencia de ruptura de la simetría de inversión temporal, el cambio de volumen lineal es proporcional a la conductividad Hall anómala, lo que sugiere posibles saltos discontinuos en el volumen en el régimen Hall cuántico.

4. Estimaciones Cuantitativas:

   • En un modelo específico de dos bandas cerca del centro de la zona, el cambio fraccional de densidad se estima en el orden de $n^{(2)}/n^{(0)} \sim 10^{-4}\%$.
   • La magnitud de este efecto escala como $\sim 1/k_F^4$, lo que implica que los sistemas con superficies de Fermi más pequeñas son más sensibles a estos cambios de densidad inducidos por el campo magnético.

Significancia y Reivindicaciones
El artículo afirma proporcionar un método riguroso e invariante de calibre para calcular la densidad promedio de electrones de Bloch en campos magnéticos débiles, extendiéndose más allá de la dinámica semiclásica.

• Avance Teórico: Unifica la descripción de las respuestas lineales y no lineales, aclarando el papel de los momentos magnéticos orbitales en metales e introduciendo un nuevo mecanismo geométrico (momentos magnéticos inducidos por el métrico cuántico) para los cambios de densidad no lineales.
• Aplicabilidad: Las expresiones analípticas derivadas están diseñadas para ser entradas de simulaciones numéricas a gran escala (por ejemplo, códigos de primeros principios) que puedan dar cuenta de estructuras cristalinas detalladas y simetrías arbitrarias.
• Contexto Experimental: Los autores sugieren que estas correcciones cuadráticas pueden ser relevantes para interpretar la respuesta no lineal de la resistividad transversal observada en experimentos a campos magnéticos más altos ($\sim 1$ T), donde los modelos semiclásicos estándar de dos portadores podrían ser insuficientes.
• Modestia: Los autores señalan que, si bien el marco analítico está establecido, las predicciones cuantitativas precisas para materiales específicos requieren simulaciones numéricas que incorporen todos los términos de la expresión general derivada. No pretenden haber verificado experimentalmente estos cambios de densidad específicos, sino que proporcionan las herramientas teóricas para hacerlo.