Condensed Matter: Graduate-Level Formalism & Experiment

Condensed Matter Physics » Graduate-Level Formalism

Graduate-Level Formalism & Experiment

This is the theory-first companion to the Condensed Matter hub. It collects the graduate-level machinery — second quantization, Green’s functions, topological band theory, and the renormalization-group view of quantum criticality — that the narrative pages use informally. Theory comes first; the experimental probes that test it are gathered at the end and treated in full on the dedicated Experimental Techniques page.

Current Research Areas

2D Materials

  • Graphene: Dirac fermions
  • Transition metal dichalcogenides
  • van der Waals heterostructures
  • Moiré superlattices

Quantum Materials

  • Weyl and Dirac semimetals
  • Axion insulators
  • Quantum spin liquids
  • Majorana fermions

Non-equilibrium Physics

  • Floquet engineering
  • Many-body localization
  • Time crystals
  • Driven-dissipative systems

Graduate-Level Mathematical Formalism

Second Quantization in Condensed Matter

Field operators for fermions:

\[\psi(\mathbf{r}) = \sum_k \phi_k(\mathbf{r}) c_k\] \[\psi^\dagger(\mathbf{r}) = \sum_k \phi_k^*(\mathbf{r}) c_k^\dagger\]

Anticommutation relations:

\[\{\psi(\mathbf{r}), \psi^\dagger(\mathbf{r}')\} = \delta(\mathbf{r} - \mathbf{r}')\] \[\{\psi(\mathbf{r}), \psi(\mathbf{r}')\} = \{\psi^\dagger(\mathbf{r}), \psi^\dagger(\mathbf{r}')\} = 0\]

General Hamiltonian:

\[H = \int d\mathbf{r} \, \psi^\dagger(\mathbf{r})\left[-\frac{\hbar^2\nabla^2}{2m} + V(\mathbf{r})\right]\psi(\mathbf{r}) + \frac{1}{2}\int d\mathbf{r} \, d\mathbf{r}' \, \psi^\dagger(\mathbf{r})\psi^\dagger(\mathbf{r}')U(\mathbf{r}-\mathbf{r}')\psi(\mathbf{r}')\psi(\mathbf{r})\]

Many-Body Green’s Functions

Definitions

Single-particle Green’s function:

\[G(\mathbf{r},t;\mathbf{r}',t') = -i\langle T[\psi(\mathbf{r},t)\psi^\dagger(\mathbf{r}',t')]\rangle\]

Spectral function:

\[A(\mathbf{k},\omega) = -2\text{Im}[G^R(\mathbf{k},\omega)]\]

Dyson equation:

\[G = G_0 + G_0 \Sigma G\]

Where $\Sigma$ is the self-energy.

Matsubara Formalism

Imaginary time:

\[G(\mathbf{r},\tau;\mathbf{r}',\tau') = -\langle T_\tau[\psi(\mathbf{r},\tau)\psi^\dagger(\mathbf{r}',\tau')]\rangle\]

Matsubara frequencies:

  • Fermions: $\omega_n = (2n+1)\pi/\beta$
  • Bosons: $\omega_n = 2n\pi/\beta$

Analytic continuation: $i\omega_n \to \omega + i\delta$

Advanced Band Theory

k·p Method

Near band extrema:

\[H = E_0 + \frac{\hbar^2 k^2}{2m^*} + \frac{\hbar}{m_0}\sum_i k_i p_i + O(k^2)\]

Kane model for narrow gap semiconductors:

\[H = \begin{pmatrix} E_c + \frac{\hbar^2 k^2}{2m_c} & Pk \\ Pk & E_v - \frac{\hbar^2 k^2}{2m_v} \end{pmatrix}\]

Wannier Functions

Construction from Bloch states:

\[w_n(\mathbf{r} - \mathbf{R}) = \frac{V}{(2\pi)^3} \int_{BZ} d\mathbf{k} \, e^{-i\mathbf{k} \cdot \mathbf{R}} \psi_{n\mathbf{k}}(\mathbf{r})\]

Maximally localized Wannier functions: Minimize spread

\[\Omega = \sum_n\left[\langle w_n|\mathbf{r}^2|w_n\rangle - \langle w_n|\mathbf{r}|w_n\rangle^2\right]\]

Topological Band Theory

Berry connection:

\[\mathbf{A}_n(\mathbf{k}) = i\langle u_{n\mathbf{k}}|\nabla_\mathbf{k}|u_{n\mathbf{k}}\rangle\]

Berry curvature:

\[\boldsymbol{\Omega}_n(\mathbf{k}) = \nabla_\mathbf{k} \times \mathbf{A}_n(\mathbf{k}) = i\sum_{m\neq n} \frac{\langle u_{n\mathbf{k}}|\nabla_\mathbf{k} H|u_{m\mathbf{k}}\rangle \times \langle u_{m\mathbf{k}}|\nabla_\mathbf{k} H|u_{n\mathbf{k}}\rangle}{(E_n - E_m)^2}\]

Z₂ invariant:

\[(-1)^\nu = \prod_{i=1}^4 \frac{\text{Pf}[w(\mathbf{k}_i)]}{\sqrt{\det[w(\mathbf{k}_i)]}}\]
Where $w_{mn} = \langle u_{m\mathbf{k}} -i\partial_{k_\mu} u_{n\mathbf{k}}\rangle$

Superconductivity: Advanced Theory

Bogoliubov-de Gennes Formalism

BdG Hamiltonian:

\[H_{BdG} = \begin{pmatrix} H_0(\mathbf{k}) & \Delta(\mathbf{k}) \\ \Delta^*(\mathbf{k}) & -H_0^*(-\mathbf{k}) \end{pmatrix}\]

Nambu spinor: $\Psi = (c_{\mathbf{k}\uparrow}, c_{-\mathbf{k}\downarrow}^\dagger)^T$

Quasiparticle spectrum:

\[E_\mathbf{k} = \pm\sqrt{\xi_\mathbf{k}^2 + |\Delta_\mathbf{k}|^2}\]

Ginzburg-Landau Theory

GL functional:

\[F = \int d^3\mathbf{r} \left[\alpha|\psi|^2 + \frac{\beta}{2}|\psi|^4 + \frac{1}{2m^*}|(-i\hbar\nabla - e^*\mathbf{A})\psi|^2 + \frac{B^2}{2\mu_0}\right]\]

GL equations:

\[\alpha\psi + \beta|\psi|^2\psi + \frac{1}{2m^*}(-i\hbar\nabla - e^*\mathbf{A})^2\psi = 0\] \[\mathbf{j} = \frac{e\hbar}{2m^*i}(\psi^*\nabla\psi - \psi\nabla\psi^*) - \frac{4e^2}{m^*}|\psi|^2\mathbf{A}\]
Coherence length: $\xi = \frac{\hbar}{\sqrt{2m^* \alpha }}$

Penetration depth: $\lambda = \sqrt{\frac{m^*}{\mu_0 4e^2 n_s}}$

Josephson Effects

Josephson relations:

\[I = I_c \sin(\phi)\] \[\frac{\partial\phi}{\partial t} = \frac{2eV}{\hbar}\]

RCSJ model:

\[C \frac{d^2\phi}{dt^2} + \frac{1}{R}\frac{d\phi}{dt} + I_c \sin(\phi) = I\]

Shapiro steps: $V_n = \frac{n\hbar\omega}{2e}$

Quantum Hall Physics

Landau Levels

Single particle states:

\[\psi_{n,m}(z) = (z - z_m)^n e^{-|z - z_m|^2/(4l_B^2)}\]

Where $l_B = \sqrt{\hbar/(eB)}$ is magnetic length.

Projected density operators:

\[\rho_q = \sum_k c_{k+q}^\dagger c_k e^{iq \times k l_B^2/2}\]

Composite Fermion Theory

CF transformation:

\[\Psi_{CF} = P_{LLL} \prod_{i<j}(z_i - z_j)^2 \Phi_{fermions}\]

Effective magnetic field:

\[B_{eff} = B - 2\phi_0\rho\]

Where $\phi_0 = h/e$ is flux quantum.

Chern-Simons Theory

Effective action:

\[S = \int d^3x \left[\frac{\epsilon^{\mu\nu\lambda}}{4\pi} a_\mu\partial_\nu a_\lambda + j^\mu a_\mu\right]\]

Statistical transmutation: Fermions $\leftrightarrow$ Bosons + flux

Strongly Correlated Electrons

Hubbard Model Extensions

t-J model (large U limit):

\[H = -t\sum_{\langle ij\rangle,\sigma} P(c_{i\sigma}^\dagger c_{j\sigma} + \text{h.c.})P + J\sum_{\langle ij\rangle}\left(\mathbf{S}_i \cdot \mathbf{S}_j - \frac{n_i n_j}{4}\right)\]

Where $P$ projects out double occupancy.

Anderson model (impurity):

\[H = \sum_{k\sigma}\epsilon_k c_{k\sigma}^\dagger c_{k\sigma} + \sum_\sigma \epsilon_d d_\sigma^\dagger d_\sigma + Un_{d\uparrow}n_{d\downarrow} + V\sum_{k\sigma}(c_{k\sigma}^\dagger d_\sigma + \text{h.c.})\]

Dynamical Mean-Field Theory (DMFT)

Self-consistency equations:

\[G_{loc}(\omega) = \sum_k G(\mathbf{k},\omega)\] \[G^{-1}(\mathbf{k},\omega) = \omega + \mu - \epsilon_\mathbf{k} - \Sigma(\omega)\] \[\Gamma(\omega) = G_0^{-1}(\omega) - G_{loc}^{-1}(\omega)\]

Anderson impurity problem:

\[H_{imp} = \epsilon_d d^\dagger d + Un_{d\uparrow}n_{d\downarrow} + \sum_k V_k(c_k^\dagger d + \text{h.c.}) + \sum_k \epsilon_k c_k^\dagger c_k\]

Slave Particle Methods

Slave boson representation:

\[c_{i\sigma} = b_i^\dagger f_{i\sigma}\]

Constraint: $b_i^\dagger b_i + \sum_\sigma f_{i\sigma}^\dagger f_{i\sigma} = 1$

Mean-field decoupling: $\langle b_i\rangle \neq 0$ describes coherent quasiparticles

Topological Phases: Advanced Topics

Topological Field Theory

Chern-Simons term:

\[S_{CS} = \frac{k}{4\pi} \int d^3x \, \epsilon^{\mu\nu\lambda} A_\mu\partial_\nu A_\lambda\]

BF theory:

\[S_{BF} = \frac{K_{IJ}}{2\pi} \int d^3x \, \epsilon^{\mu\nu\lambda} a_\mu^I\partial_\nu a_\lambda^J\]

Topological Order

Ground state degeneracy on torus: Depends on topology

Modular matrices: $S$ and $T$ characterize anyon statistics

\[S_{ab} = \frac{1}{\mathcal{M}} \sum_c \frac{N_{ab}^c d_c}{d_a d_b}\]

Topological entanglement entropy:

\[S = \alpha L - \gamma\]

Where $\gamma = \ln(\mathcal{M})$ is universal.

Symmetry-Protected Topological Phases

Classification by cohomology: $H^{d+1}(G, U(1))$

Matrix product state representation:

\[|\psi\rangle = \sum_{s_1...s_N} \text{Tr}[A^{s_1}...A^{s_N}]|s_1...s_N\rangle\]

Symmetry: $u(g)A^s u^\dagger(g) = \sum_{s’} U(g)_{ss’} A^{s’}$

Quantum Phase Transitions & Scaling

A quantum phase transition (QPT) is a transition between distinct ground states of a many-body system, driven not by temperature but by a non-thermal coupling $g$ (pressure, magnetic field, doping, chemical potential) at strictly $T = 0$. At the critical coupling $g_c$ the ground-state energy is non-analytic and a characteristic energy scale — the gap $\Delta$ — vanishes. Although the transition itself sits at $T=0$ and is unreachable in experiment, the quantum critical point (QCP) controls a wedge-shaped region of the finite-$T$ phase diagram, the quantum critical fan, where the only relevant scale is temperature itself. This is why a $T=0$ singularity leaves fingerprints (anomalous transport, $T$-linear resistivity, non-Fermi-liquid specific heat) at experimentally accessible temperatures.

Diverging length and time

Approaching $g_c$ a correlation length and a correlation time both diverge:

\[\xi \sim |g - g_c|^{-\nu}, \qquad \xi_\tau \sim \xi^{z} \sim |g - g_c|^{-\nu z}\]

The exponent $\nu$ is the same correlation-length exponent as in a classical transition, but space and time scale differently, encoded in the dynamic critical exponent $z$. Equivalently the characteristic energy (the gap, or $\hbar/\xi_\tau$) softens as

\[\Delta \sim |g - g_c|^{\nu z}, \qquad \omega \sim k^{z}.\]

For a Lorentz-invariant critical theory $z = 1$ (energy scales like momentum); for a metallic QCP with overdamped order-parameter dynamics one typically finds $z = 2$ or $z = 3$, reflecting Landau damping by the particle-hole continuum.

Quantum-to-classical mapping

The Feynman path integral for a $d$-dimensional quantum system at $T=0$ is a statistical-mechanics problem in $d + z$ effective dimensions: imaginary time $\tau \in [0, \beta\hbar)$ acts as $z$ extra spatial directions. A $T=0$ quantum transition in $d$ space dimensions therefore maps onto a classical transition in

\[d_{\text{eff}} = d + z\]

dimensions. The mapping is exact for many models (e.g. the transverse-field Ising chain $\leftrightarrow$ the 2D classical Ising model, with $z=1$). Finite temperature plays the role of a finite system size $L_\tau = \beta\hbar$ in the imaginary-time direction, so

\[\xi_\tau \lesssim \beta\hbar \;\Longrightarrow\; T \gtrsim \frac{\hbar}{\xi_\tau} \sim |g - g_c|^{\nu z},\]

and the boundary $T \sim |g - g_c|^{\nu z}$ traces the edges of the quantum critical fan.

Scaling hypothesis and exponents

Near the QCP the singular part of the free-energy density obeys a scaling form. With reduced coupling $t = (g - g_c)/g_c$ and a symmetry-breaking field $h$,

\[f_s(t, h, T) = b^{-(d+z)}\, f_s\!\left(b^{1/\nu}\,t,\; b^{y_h}\,h,\; b^{z}\,T\right),\]

for an arbitrary rescaling factor $b$. Choosing $b$ to remove one argument generates the standard exponents. Choosing $b = |t|^{-\nu}$ gives the order parameter, gap, and correlation functions; choosing $b = T^{-1/z}$ collapses finite-temperature data. The familiar critical exponents are then fixed by $\nu$, $z$, and the field dimension $y_h$:

Exponent Definition Scaling relation    
$\nu$ $\xi \sim t ^{-\nu}$ correlation length
$z$ $\xi_\tau \sim \xi^{z}$ dynamic / time    
$\beta$ $m \sim t ^{\beta}$ order parameter
$\gamma$ $\chi \sim t ^{-\gamma}$ susceptibility
$\eta$ $G(k) \sim k^{-2+\eta}$ anomalous dimension    
$\delta$ $m \sim h^{1/\delta}$ critical isotherm    

Only two are independent; the rest follow from hyperscaling, which for a QPT uses the effective dimension $d + z$:

\[2 - \alpha = \nu(d + z), \qquad \gamma = \nu(2 - \eta), \qquad \beta = \tfrac{1}{2}\nu(d + z - 2 + \eta).\]

Finite-size / finite-$T$ scaling. On a system of linear size $L$ the magnetization obeys

\[M(t, h, L) = L^{-\beta/\nu}\, f\!\left(t\,L^{1/\nu},\; h\,L^{y_h}\right),\]

and at finite temperature the same form holds with $L_\tau = \beta\hbar$ in the time direction — the practical route to extracting exponents from quantum Monte Carlo.

Renormalization group and universality

The exponents follow from RG flow near the QCP. Linearizing the RG transformation about the fixed point, the eigenvalues $\lambda_i = b^{y_i}$ of the relevant couplings fix $\nu = 1/y_t$ and $y_h$. The upper critical dimension is $d_c^{+} = 4$ for the effective dimension, i.e. $d + z = 4$ separates the mean-field regime ($d + z > 4$, Gaussian fixed point, classical exponents with logarithmic corrections) from the fluctuation-dominated regime ($d + z < 4$, Wilson–Fisher fixed point with non-trivial exponents). For a 2D antiferromagnet ($d=2$, $z=1$) one has $d + z = 3 < 4$, so fluctuations matter and the transition is in the 3D classical Heisenberg/O(3) universality class.

Universality is the central payoff: $\nu$, $z$, $\eta$ depend only on the spatial dimension, the order-parameter symmetry, and the range of interactions — not on microscopic details. Systems as different as a transverse-field magnet and a superfluid-insulator transition share exponents if they share a fixed point.

Representative quantum critical points

  • Transverse-field Ising chain $H = -J\sum_i \sigma_i^z\sigma_{i+1}^z - h\sum_i \sigma_i^x$. Exactly solvable by Jordan–Wigner; the QCP at $h = J$ has $z = 1$, $\nu = 1$, $\eta = 1/4$, mapping to the 2D classical Ising model.
  • Bose–Hubbard / superfluid–Mott insulator. At integer filling the tip of the Mott lobe is a relativistic $z=1$ transition in the 3D XY class; away from the tip, density fluctuations give $z = 2$ (mean-field, $d + z = 4$).
  • Heavy-fermion antiferromagnetic QCP. Field- or pressure-tuned (e.g. CeCu$_6$, YbRh$_2$Si$_2$). Spin-density-wave (Hertz–Millis) treatment gives $z = 2$ (AFM) or $z = 3$ (FM); Kondo-breakdown scenarios predict local quantum criticality with $\omega/T$ scaling.
  • Dilute Bose gas / magnetization onset. A field-tuned transition with $z = 2$, $\nu = 1/2$ — the canonical $d + z = d + 2$ example.

Beyond Landau: deconfined quantum criticality

The Landau–Ginzburg–Wilson paradigm assumes a single order parameter. Some quantum transitions evade it. The Néel–VBS transition of a 2D quantum antiferromagnet is between two different broken symmetries (antiferromagnetic Néel order and valence-bond-solid order) that share no group-subgroup relation, so no single order-parameter field describes both phases. Instead the critical theory is written in terms of fractionalized spinons $z_\alpha$ ($CP^{1}$ field) coupled to an emergent U(1) gauge field $a_\mu$:

\[S = \int d^2x \, d\tau \left[|(\partial_\mu - i a_\mu)z|^2 + s\,|z|^2 + u\,(|z|^2)^2\right] + \frac{1}{2e^2}\int d^2x\, d\tau\, f_{\mu\nu}^2.\]

The spinons are confined on both sides of the transition (recombining into magnons or VBS order) but become deconfined exactly at the critical point — a genuinely non-Landau, “beyond order-parameter” continuous transition.

Effective Field Theory in Condensed Matter

The renormalization group reframes condensed-matter problems as effective field theories (EFTs): rather than tracking every electron, one writes the most general local action consistent with the system’s symmetries for the relevant low-energy degrees of freedom (the order parameter, a Goldstone mode, a gauge field, a Dirac cone), then organizes terms by their importance under coarse-graining. The same logic that gives Ginzburg–Landau theory its $|\psi|^4$ form underlies the entire low-energy description of quantum matter.

Symmetries fix the action

The construction is symmetry-first. One identifies the slow fields and the exact and emergent symmetries (translations, rotations, time reversal, particle–hole, a global U(1) or O(N), gauge invariance), and writes every local term compatible with them. Symmetry decides which terms may appear; the RG then decides which ones matter. Two consequences are immediate:

  • Goldstone modes. Spontaneously broken continuous symmetries guarantee gapless modes whose action is fixed by the broken symmetry — e.g. magnons in a Heisenberg magnet, the phase mode of a superfluid, $\mathcal{L} = \tfrac{1}{2}\rho_s(\nabla\theta)^2 + \ldots$
  • Forbidden couplings. Time-reversal forbids a Chern–Simons term; particle–hole symmetry constrains the BdG Hamiltonian; lattice symmetry can protect band touchings (Dirac/Weyl points) against gapping.

Scaling dimensions and relevance

Under a coarse-graining step $x \to b\,x$, $\tau \to b^{z}\tau$ each field $\phi$ acquires a scaling dimension $[\phi] = \Delta_\phi$ fixed by demanding the leading (Gaussian) part of the action be scale-invariant. For a relativistic scalar in $d+z$ effective dimensions, requiring $\int d^d x\, d\tau\,(\partial\phi)^2$ to be dimensionless gives

\[\Delta_\phi = \frac{d + z - 2}{2}.\]

A coupling $g_n$ multiplying an operator $\mathcal{O}_n$ of dimension $\Delta_n$ then has dimension $[g_n] = (d + z) - \Delta_n$, and flows as

\[g_n(b) = b^{(d+z) - \Delta_n}\, g_n.\]

This single counting classifies every term:

  Condition Behavior under RG Example    
Relevant $\Delta_n < d+z$ grows; drives system off fixed point mass term $s \psi ^2$
Marginal $\Delta_n = d+z$ logarithmic; needs higher order $ \psi ^4$ at $d+z=4$
Irrelevant $\Delta_n > d+z$ decays; drops out at low energy gradient$^4$, $ \psi ^6$ in $d{+}z{<}3$

The handful of relevant and marginal couplings constitute the universal data; the infinite tower of irrelevant operators only supplies non-universal corrections to scaling. This is why universality holds — microscopic complexity is encoded in irrelevant operators that the flow discards.

RG flow and fixed points

Combining the linear scaling with loop corrections gives the beta functions $\beta_n = dg_n/d\ell$ (with $\ell = \ln b$). Their zeros are fixed points, where the theory is scale-invariant:

  • The Gaussian fixed point $g_n = 0$ governs the mean-field regime above the upper critical dimension $d + z > 4$.
  • A non-trivial Wilson–Fisher fixed point appears below it. The classic $\phi^4$ example, organized in $\epsilon = 4 - (d+z)$, has at one loop
\[\beta_u = -\epsilon\, u + \frac{(N+8)}{8\pi^2}\, u^2,\]

with an infrared-stable zero at $u^{*} = 8\pi^2\epsilon/(N+8)$. Linearizing about a fixed point, the eigenvalues $y_i$ of the relevant directions fix the universal exponents ($\nu = 1/y_t$, etc.), tying this EFT section directly back to the scaling exponents above.

Emergent IR symmetries are common: a lattice model with only discrete rotation symmetry can flow to a Lorentz- and conformally invariant fixed point, and gauge fields ($CP^{1}$, Chern–Simons, $\mathbb{Z}_2$) routinely emerge as low-energy descriptions of fractionalized phases even though the microscopic Hamiltonian has no gauge structure at all.

Computational Methods

Density Functional Theory for Solids

Kohn-Sham equations:

\[\left[-\frac{\hbar^2\nabla^2}{2m} + v_{eff}(\mathbf{r})\right]\phi_i(\mathbf{r}) = \epsilon_i\phi_i(\mathbf{r})\]

Exchange-correlation functionals:

  • LDA: $\epsilon_{xc}[n] = \epsilon_{xc}(n)$
  • GGA: $\epsilon_{xc}[n,\nabla n]$
  • Hybrid: Mix exact exchange

Band structure calculations: Plane wave basis, pseudopotentials

Quantum Monte Carlo

Variational QMC:

\[E = \frac{\langle\Psi_T|H|\Psi_T\rangle}{\langle\Psi_T|\Psi_T\rangle}\]

Diffusion QMC: Project out ground state

\[|\Psi_0\rangle = \lim_{t\to\infty} e^{-Ht}|\Psi_T\rangle\]

Sign problem: Constrains fermionic/frustrated systems

Tensor Network Methods

iPEPS for 2D systems:

\[|\Psi\rangle = \sum_s \text{tTr}[A^{s_{1,1}}...A^{s_{N,N}}]|s\rangle\]

Corner transfer matrix: Compute observables

Time evolution: TEBD, MPO methods

import numpy as np
from scipy.linalg import expm

def tebd_step(psi, U_bonds, chi_max):
    """Time-evolving block decimation step"""
    for bond in range(0, len(psi)-1, 2):  # Even bonds
        psi = apply_two_site_gate(psi, U_bonds[bond], bond, chi_max)
    for bond in range(1, len(psi)-1, 2):  # Odd bonds  
        psi = apply_two_site_gate(psi, U_bonds[bond], bond, chi_max)
    return psi

def apply_two_site_gate(psi, U, bond, chi_max):
    """Apply two-site gate with truncation"""
    # Contract tensors
    theta = np.tensordot(psi[bond], psi[bond+1], axes=([2],[0]))
    theta = np.tensordot(U, theta, axes=([2,3],[0,2]))
    
    # SVD and truncate
    theta = theta.transpose(0,2,1,3).reshape(d*chi_l, d*chi_r)
    u, s, vh = np.linalg.svd(theta, full_matrices=False)
    
    # Truncate to chi_max
    chi_new = min(len(s), chi_max)
    u = u[:, :chi_new]
    s = s[:chi_new]
    vh = vh[:chi_new, :]
    
    # Update MPS tensors
    psi[bond] = u.reshape(chi_l, d, chi_new)
    psi[bond+1] = (np.diag(s) @ vh).reshape(chi_new, d, chi_r)
    
    return psi

Research Frontiers

Quantum Materials Design

Materials informatics: Machine learning for materials discovery

Heterostructure engineering: Designer quantum phases

Moiré systems: Tunable strongly correlated physics

Non-equilibrium Phenomena

Floquet engineering: Light-induced topological phases

\[H_F = H_0 + V \cos(\omega t)\]

Ultrafast spectroscopy: Pump-probe dynamics

Many-body localization: Breakdown of thermalization

Quantum Technologies

Topological quantum computing: Anyonic braiding

Quantum sensors: NV centers, SQUIDs

Coherent quantum devices: Josephson junctions, quantum dots

Unconventional Superconductivity

Iron-based superconductors: Multi-orbital physics

Heavy fermion superconductors: Quantum criticality

Organic superconductors: Low dimensionality

Interface superconductivity: STO/LAO, FeSe/STO

Correlated Topology

Twisted bilayer graphene: Flat bands and superconductivity

Magnetic topological insulators: Quantum anomalous Hall effect

Weyl-Kondo semimetals: Topology meets strong correlations

Experimental Probes of the Formalism

Theory is only as good as the spectra it predicts. The probes below close the loop — each measures a quantity that the formalism above computes directly (a spectral function, a Fermi-surface cross-section, a local density of states). They are summarized here for completeness; the full methodology, geometry, and worked data-analysis examples live on the dedicated Experimental Techniques page, which also covers transport, neutron and X-ray scattering, and thermodynamic measurements.

Transport, scattering, and thermodynamics at a glance

  • Transport. Four-probe resistivity, the Hall effect (carrier density and mobility), and quantum oscillations to map the Fermi surface.
  • Spectroscopy. ARPES, STM/STS, neutron scattering (magnetic structure and excitations), and X-ray scattering (crystal structure).
  • Thermodynamics. Specific heat (locating phase transitions), magnetization, and thermal expansion (lattice coupling).

ARPES (Angle-Resolved Photoemission)

Directly images the single-particle spectral function $A(\mathbf{k},\omega)$ of the Green’s-function section above.

Intensity:

\[I(\mathbf{k},\omega) \propto |M_{fi}|^2 f(\omega) A(\mathbf{k},\omega)\]

Where $M_{fi}$ is matrix element, $f(\omega)$ is Fermi function.

Self-energy extraction:

\[\Sigma'(\mathbf{k},\omega) = \omega - \epsilon_\mathbf{k}^0 - \text{Re}[\Sigma(\mathbf{k},\omega)]\] \[\Sigma''(\mathbf{k},\omega) = \text{Im}[\Sigma(\mathbf{k},\omega)]\]

Quantum Oscillations

Lifshitz-Kosevich formula:

\[M \propto \left(\frac{T}{B}\right)^{1/2} R_T R_D R_S \sin\left(\frac{2\pi F}{B} + \phi\right)\]

Where:

  • $R_T$ = thermal damping
  • $R_D$ = Dingle factor
  • $R_S$ = spin factor
  • $F$ = oscillation frequency

Fermiology: Extract Fermi surface, effective mass, scattering rate

STM/STS

Tunneling current:

\[I \propto \int_{-eV}^0 d\omega \, \rho_s(\omega)\rho_t(\mathbf{r},\omega+eV)T(\omega,eV)\]

Differential conductance:

\[\frac{dI}{dV} \propto \rho_s(E_F)\rho_t(\mathbf{r},eV)\]

Quasiparticle interference: Fourier transform reveals $\mathbf{q} = \mathbf{k}_f - \mathbf{k}_i$

References and Further Reading

Classic Textbooks

  1. Ashcroft & Mermin - Solid State Physics
  2. Kittel - Introduction to Solid State Physics
  3. Mahan - Many-Particle Physics
  4. Abrikosov, Gorkov & Dzyaloshinski - Methods of Quantum Field Theory in Statistical Physics

Advanced Monographs

  1. Coleman - Introduction to Many-Body Physics
  2. Wen - Quantum Field Theory of Many-Body Systems
  3. Bernevig & Hughes - Topological Insulators and Topological Superconductors
  4. Tinkham - Introduction to Superconductivity

Specialized Topics

  1. Giamarchi - Quantum Physics in One Dimension
  2. Sachdev - Quantum Phase Transitions
  3. Girvin & Yang - Modern Condensed Matter Physics
  4. Phillips - Advanced Solid State Physics

Recent Reviews

  1. Keimer et al. - From quantum matter to high-temperature superconductivity in copper oxides (2015)
  2. Armitage, Mele & Vishwanath - Weyl and Dirac semimetals in three-dimensional solids (2018)
  3. Balents et al. - Superconductivity and strong correlations in moiré flat bands (2020)
  4. Khajetoorians et al. - Creating designer quantum states of matter atom-by-atom (2019)
← Superconductivity, Quantum Hall & Topological Phases Condensed Matter Physics (Hub) →

See Also