From Quantum Mechanics to Materials

A material theory is not obtained by making a one-electron Schrödinger equation larger. The microscopic problem already contains many electrons, nuclei, Fermi statistics, Coulomb interactions, boundaries, and a state. One-particle orbitals, bands, phonons, screened interactions, and lattice Hamiltonians appear only after specific reorganizations or reductions of that problem.

This page teaches how to audit one such chain. For every arrow from microscopic ingredients to a measured quantity, record what was retained, what was eliminated, where each parameter came from, why the step is justified, and what observation could make it fail. The result is a provenance ledger: a bounded material claim rather than a slogan about deriving everything from first principles.

Computational Quantum Matter takes the declared model and observable from that ledger into method routing, numerical validation, probe comparison, and a bounded computational conclusion; this page retains the physical reduction and parameter provenance.

When a reduction targets paired or coherent-flow matter, Superfluidity and Superconductivity owns the choice of phase, response, microscopic, defect, or device route; this page continues to own the provenance of the reduction feeding that route.

When a reduction targets a topological material claim, carry its retained subspace, symmetry, gap or node, interaction, boundary, and probe assumptions into Topological Quantum Matter, which selects the local specialist branch; this page continues to own the provenance of the reduction.

Helpful background. The Time-Independent Schrödinger Equation supplies the eigenvalue language used below, while Identical Particles and Exchange Symmetry explains fermionic antisymmetry. Readers wanting the derivational machinery can follow Many-Particle Hamiltonians and Effective Hamiltonians and Scale Separation. None is required merely to use the ledger.

The Electron–Ion Problem Comes First

Within a declared nonrelativistic Coulomb model, organize the microscopic Hamiltonian as

$$H_{\mathrm{micro}} =H_{\mathrm{Coul}}+H_{\mathrm{rel}}+H_{\mathrm{ext}},$$

where

$$\begin{aligned} H_{\mathrm{Coul}} &=T_{\mathrm e}+T_{\mathrm I}+V_{\mathrm{eI}}\\ &\quad+V_{\mathrm{ee}}+V_{\mathrm{II}}. \end{aligned}$$

Here $H_{\mathrm{Coul}}$ collects the nonrelativistic Coulomb terms: $T_{\mathrm e}$ and $T_{\mathrm I}$ are electronic and ionic kinetic energies, while the three $V$ terms are electron–ion, electron–electron, and ion–ion interactions. The term $H_{\mathrm{rel}}$ collects retained relativistic corrections such as spin–orbit coupling, and $H_{\mathrm{ext}}$ describes declared fields or environments. Composition, charge state, geometry, boundary conditions, and thermodynamic preparation are also part of the problem.

Even this is not the ultimate exact theory. Nuclear structure, radiation, pair creation, and fully relativistic quantum electrodynamics have been excluded. “Exact” below therefore always means exact within a stated model, Hilbert space, and numerical representation.

The electronic state must be antisymmetric under exchange of identical electrons. That statement is exact within ordinary fermionic quantum mechanics, but it does not account for Coulomb correlation. A Slater determinant builds antisymmetry into a trial state; it is not generally the exact state of an interacting material. Likewise, translation symmetry may label a many-body state by total crystal momentum without factorizing it into independent one-electron Bloch orbitals.

The useful question is therefore not “How does one particle become a material?” Define

$$\begin{aligned} X_{\mathrm{micro}} &\equiv(\mathcal C,H_{\mathrm{micro}},\rho,O),\\ X_{\mathrm{work}} &\equiv(H_{\mathrm{work}},\rho_{\mathrm{work}},\\ &\qquad O_{\mathrm{work}};\epsilon_{\mathcal R}), \end{aligned}$$

and audit the chain

$$\begin{aligned} X_{\mathrm{micro}} &\xrightarrow{\mathcal R}X_{\mathrm{work}},\\ X_{\mathrm{work}} &\xrightarrow{\mathcal S}\langle O_{\mathrm{work}}\rangle \xrightarrow{\mathcal M}d_{\mathrm{pred}}. \end{aligned}$$

$\mathcal C$ declares composition and conditions; $\rho$ is the prepared microscopic state or ensemble; and $O$ is the microscopic observable. The reduction $\mathcal R$ produces a working Hamiltonian, state, and usually a dressed or projected observable. The solver $\mathcal S$ evaluates that working problem, while $\mathcal M$ maps the result through matrix elements, geometry, resolution, and instrument response to the predicted data record $d_{\mathrm{pred}}$. The reduction remainder $\epsilon_{\mathcal R}$ is meaningful only after its norm, observable, scale, and regime have been named.

The Material-Model Provenance Ledger

Emergence and Effective Degrees of Freedom owns the general ledger for retained variables, matching, errors, and breakdown. A material calculation adds four provenance questions that cannot be inferred from the final Hamiltonian alone.

Import three generic fields from that owner rather than redefining them: retained versus eliminated degrees of freedom; map status and matching; and controls, errors, breakdown, and stopping. The map taxonomy is exact reorganization, controlled approximation, uncontrolled but testable modeling approximation, or phenomenological/inferential step. Numerical convergence is an orthogonal solver test, not another type of physical arrow.

Then add four material-specific entries.

1. Material target and endpoints. Name the input composition, structure, environment, state, scale, and boundary data. Then name the output observable, probe, resolution, and desired accuracy. “Understand silicon” is not a target; “predict low-field conductivity over a declared temperature and carrier-density window” is.
2. Structural and active-space provenance. Record the crystal or sample idealization, core and valence partition, active bands, orbitals, spins, valleys, lattice modes, boundaries, and sectors. Identify remote bands, high-energy charge configurations, defects, baths, or fast modes projected out, averaged over, frozen, or ignored.
3. Parameter and screening provenance. Mark every parameter as bare, symmetry-fixed, calculated, screened, downfolded, renormalized, or fitted, including its basis, frequency window, and scale. Track which polarization channels the retained solver will treat.
4. Operator-to-probe map. State how the microscopic observable is projected or dressed and how matrix elements, sample geometry, contacts, selection rules, backgrounds, and instrumental resolution enter $d_{\mathrm{pred}}$.

Finally apply the generic stopping rule: separate model-reduction error, solver or discretization error, finite-size and boundary error, and measurement or calibration uncertainty. Give a gap ratio, mass ratio, convergence sequence, sum rule, cross-method benchmark, or falsifying observation for every non-exact arrow. Stop only when the combined uncertainty meets the target and every omitted effect is bounded within the residual budget, varied in a sensitivity analysis, or explicitly excluded from the licensed domain. Otherwise the omission blocks the quantitative claim.

The labels are not a prestige ranking. An exact basis transformation may be useless for the chosen observable, while a simple phenomenological model may be the most reliable description in a carefully calibrated window.

First Reductions: Nuclei, Periodicity, and Effective Electrons

Electronic–nuclear separation

For fixed nuclear coordinates $\mathbf R$, define a conditional electronic problem,

$$\begin{aligned} H_{\mathrm e}(\mathbf R)\Phi_\alpha &=E_\alpha(\mathbf R)\Phi_\alpha,\\ U_\alpha(\mathbf R) &=E_\alpha(\mathbf R)+V_{\mathrm{II}}(\mathbf R). \end{aligned}$$

The Born–Oppenheimer approximation neglects or organizes derivative couplings between electronic surfaces by exploiting nuclear–electronic scale separation. It does not declare nuclei permanently classical or motionless. Small electronic gaps, fast driving, light nuclei, and strong nonadiabatic coupling can reopen the eliminated electronic channels.

An equilibrium structure is a minimum of an appropriate energy or free-energy surface under specified conditions. Replacing the actual finite, defective, thermally fluctuating sample by an infinite periodic crystal is another modeling step. Periodicity enables crystal-momentum organization, but omits surfaces, disorder, strain gradients, and finite geometry unless they are restored explicitly.

Why one-particle equations reappear

Mean-field, density-functional, quasiparticle, and fitted band descriptions often use an auxiliary eigenproblem of the form

$$h_{\mathrm{eff}}(\mathbf k;\varepsilon_{n\mathbf k}) \phi_{n\mathbf k} =\varepsilon_{n\mathbf k}\phi_{n\mathbf k}.$$

This is a reduced description of a many-electron problem, not a return to literal independent electrons. Depending on the method, $h_{\mathrm{eff}}$ can be a local density-dependent Kohn–Sham operator, a nonlocal Fock operator, a matrix or spinor Hamiltonian, or an energy-dependent quasiparticle operator. The meaning of $\varepsilon_{n\mathbf k}$ changes with that choice. In particular, a Kohn–Sham eigenvalue is not generically a charged excitation energy, and a self-consistent orbital does not by itself certify the neglected correlation.

For a periodic effective one-body Hamiltonian, Bloch’s theorem organizes solutions into bands. The symmetry decomposition is exact for that Hamiltonian; the construction of the Hamiltonian and the interpretation of its bands are separate ledger entries. Carry the description provenance, state and filling, gap or crossing claim, target observable, and validity window into the Band Theory and Electronic Structure gateway for downstream interpretation. The Band Theory Overview owns the detailed band-to-many-body dictionary.

Screening and parameter matching

Eliminated electronic polarization changes the interaction seen by retained carriers. Use the convention $P=\delta n/\delta U_{\mathrm{tot}}$, for which stable static electronic polarization is negative and the Dyson equation is $W_r=v+vP_rW_r$. For a target subspace $\mathcal T$, separate the irreducible polarization schematically as

$$\begin{aligned} P&=P_{\mathcal T}+P_r,\\ W_r(\omega) &=\left(1-vP_r(\omega)\right)^{-1}v. \end{aligned}$$

Here $v$ is the bare Coulomb kernel, $P_r$ contains polarization outside the retained target processes, and $W_r$ is the interaction screened only by that eliminated polarization. Products and inverses may be matrices in position, orbital, reciprocal-lattice, and frequency labels. Projected interaction parameters are matrix elements such as $U_{abcd}(\omega)=\langle ab|W_r(\omega)|cd\rangle$. The fully screened interaction would also resum $P_{\mathcal T}$; if the low-energy solver treats that polarization explicitly, it must be excluded from $W_r$ to avoid counting it twice. Replacing the resulting nonlocal, frequency-dependent interaction by a dielectric constant or an onsite $U$ is a further scale- and basis-dependent approximation. The Random Phase Approximation gives one organizational scheme in appropriate regimes; it is not a universal definition of screening.

Combining a screened interaction calculated for one retained subspace with hoppings or self-energies derived from another can double count or omit polarization. Parameter provenance is therefore part of the physics, not administrative metadata.

Lattice Motion and the Return of Nuclear Degrees of Freedom

Expanding an adiabatic potential surface about equilibrium positions $\mathbf R^0$, write $\Delta U_0(\mathbf u)\equiv U_0(\mathbf R^0+\mathbf u)-U_0(\mathbf R^0)$. Then

$$\Delta U_0(\mathbf u)=U_{\mathrm{harm}}(\mathbf u)+\cdots,$$

where

$$\begin{aligned} U_{\mathrm{harm}}(\mathbf u) &\equiv\frac12\sum_{I\alpha,J\beta}\\ &\quad\Phi_{I\alpha,J\beta}u_{I\alpha}u_{J\beta}. \end{aligned}$$

The force-constant matrix $\Phi$, together with ionic masses, defines harmonic normal modes; quantizing those modes gives phonons. Carry the retained lattice variables, approximation order, state, target observable, and failure test into the Lattice Vibrations and Collective Modes gateway when the pipeline continues to anharmonic, electron–phonon, carrier-dressing, charge-mode, or bound-pair physics. The harmonic truncation, any finite-range force-constant truncation, and the chosen treatment of anharmonicity are approximations. The normal-mode diagonalization of the declared quadratic Hamiltonian is exact.

Electronic parameters also depend on displacement. For a normal coordinate $Q_\lambda$,

$$\begin{aligned} t_{ab}(\mathbf u) &=t_{ab}^{(0)}\\ &\quad+\sum_\lambda \left.\frac{\partial t_{ab}}{\partial Q_\lambda}\right|_0 Q_\lambda+\cdots. \end{aligned}$$

The derivative supplies an electron–phonon vertex after the coordinate is quantized. This makes clear why “fixed ions” and “phonons” are not competing theories: the first is a conditional electronic step, while the second restores controlled fluctuations about the resulting structure.

Uniform translation cannot change the energy of a free bulk crystal. With $l$ labeling cells and $\kappa$ a basis ion, the harmonic force constants obey

$$\sum_{l\kappa'} \Phi_{0\kappa\alpha,l\kappa'\beta}=0.$$

The sum covers every cell and basis ion. A substrate, pinning potential, or other external constraint changes this statement. In an otherwise free-crystal calculation, violation is a useful diagnostic of inconsistent forces, truncation, or numerical error.

From Bands to a Retained Material Hamiltonian

Suppose a composite set of bands is isolated over the energy and momentum window of interest. A unitary transformation within that complete retained set is an exact representation change. Wannier Functions owns the exact Bloch-to-localized-basis construction and its gauge, localization, and obstruction tests. Exponentially localized Wannier orbitals require a sufficiently regular, topologically unobstructed frame; a nonzero Chern bundle forbids such a basis, and disentangling an overlapping energy window introduces another modeling choice. Selecting bands from the full space, truncating matrix-element range, and replacing dynamic interactions by static ones are approximations.

Before that reduction, Band Structure Workflows owns the crystal-to-band execution record and the validated full-zone source states, projectors, and provenance passed downstream. This page retains the end-to-end distinction between exact representation changes and material-model approximations.

A common working Hamiltonian has the schematic decomposition

$$\begin{aligned} H_{\mathrm{work}} &=H_{\mathrm{el}}+H_{\mathrm{ph}}\\ &\quad+H_{\mathrm{e-ph}}+\cdots, \end{aligned}$$

with

$$H_{\mathrm{el}}=H_t+H_U,$$ $$H_t=\sum_{ab}t_{ab}c_a^\dagger c_b,$$ $$H_U=\frac12\sum_{abcd}U_{abcd} c_a^\dagger c_b^\dagger c_d c_c,$$

and

$$H_{\mathrm{ph}} =\sum_\lambda\hbar\omega_\lambda \left(b_\lambda^\dagger b_\lambda+\frac12\right).$$

The fermionic operators $c_a^\dagger,c_a$ create and remove an electron in retained state $a$, which must specify orbital, cell, spin, and any other active label. The bosonic operators $b_\lambda^\dagger,b_\lambda$ act on retained phonon mode $\lambda$. The tight-binding and effective-Hamiltonian pages own the detailed constructions.

The same retained manifold can support different reductions. Weak residual interactions may justify a quasiparticle or screened-band description. Near an integer filling with a narrow band, a Hubbard model may expose competition between hopping and local repulsion. For the repulsive single-band Hubbard model at one electron per site with $U\gg|t|$, a two-site or bipartite nearest-neighbor reduction gives an antiferromagnetic exchange $J=4|t|^2/U+O(|t|^4/U^3)$ when direct exchange, multiorbital terms, and additional hopping paths are absent. That is a controlled statement about a sector and energy window, not an exact identity between the material and a spin model. Hubbard Physics in Materials owns the material-specific continuation.

When a projection produces magnetic degrees of freedom, Magnetism and Spin Systems routes the resulting moment, exchange, ordered-state, excitation, and probe claims. This page retains the provenance and validity audit for the projection itself.

Worked Ledger: Low-Field Electrons in Crystalline Silicon

Choose one illustrative bounded target: predict the longitudinal four-probe dc conductivity of a homogeneous, unstrained, bulk phosphorus-doped n-type silicon bar from $250$ to $350\,\mathrm K$, with an independently measured electron density from $10^{15}$ to $10^{17}\,\mathrm{cm}^{-3}$, compensation below $1\%$, and independently characterized ionized-impurity density, to within $10\%$. Require an Ohmic field window with negligible Joule heating and $eE\ell_{\mathrm{mfp}}\ll k_{\mathrm B}T$, where $\ell_{\mathrm{mfp}}$ is the carrier mean free path. These values define a ledger exercise, not a universal silicon benchmark.

1. Composition and microscopic scope. Begin with silicon nuclei and electrons in the nonrelativistic Coulomb model, plus any required spin–orbit and external-field terms. Charge neutrality, isotopic composition, dopants, sample geometry, and temperature belong in the declaration. This is the microscopic model, not yet a computational method.

2. Equilibrium crystal. Use electronic–nuclear scale separation to obtain an adiabatic surface and a diamond-lattice equilibrium structure. The Born–Oppenheimer step is approximate; minimizing the chosen surface is a numerical task. Compare lattice constants, forces, phonon stability, and energy differences against convergence tests and experiment. A defect-rich, strained, nanoscale, or strongly driven sample can invalidate the perfect-bulk replacement.

3. Effective electronic structure. Replace core and high-energy electronic structure by a declared effective description and solve for periodic valence and conduction states. The Bloch decomposition is exact for the resulting periodic one-body operator; the core treatment, exchange-correlation approximation, basis cutoff, and interpretation of eigenvalues are not. Converge basis and sampling, and benchmark the band gap and valley structure relevant to the target.

4. Low-energy projection. Unstrained bulk silicon has six symmetry-related conduction minima along the $\Delta$ directions, the equivalent $\langle100\rangle$ valley axes. For sufficiently weak fields and low carrier excess energy, retain a valley label $\nu$ and expand about each minimum $\mathbf k_{0\nu}$. Define $\Delta\varepsilon_\nu(\mathbf q)\equiv\varepsilon_\nu(\mathbf k_{0\nu}+\mathbf q)-\varepsilon_c$; then

$$\Delta\varepsilon_\nu(\mathbf q) =\frac{\hbar^2}{2} q_i(m_\nu^{*-1})_{ij}q_j +\cdots.$$

Each effective-mass tensor has longitudinal and transverse principal axes and is matched to local band curvature. Conductivity sums the populated valley contributions. This truncation eliminates remote bands and nonparabolic terms. Its failure tests include occupation far from the minima, strong confinement, large fields, or target resolution comparable with neglected splittings; intervalley phonon scattering is instead a standard retained process whenever the temperature and carrier distribution activate it.

5. Scattering and lattice motion. Restore phonon-mediated intravalley and, when active, intervalley scattering, along with impurities and other relevant channels. Electron–phonon vertices come from displacement derivatives of the electronic Hamiltonian; impurity parameters require their own structural and screening model. A relaxation time $\tau$ is not a universal material constant: it depends on energy, temperature, carrier density, valley, scattering channel, and approximation. Converge band and phonon interpolation meshes, check positivity and detailed balance of the collision kernel, compare a relaxation-time closure with the full linearized collision operator, and test the predicted temperature and density dependence at points not used for parameter matching.

6. Transport solver. Use the Transport, Response, and Optics gateway to decide whether the observable calls for a bulk coefficient, kinetic closure, correlation kernel, or terminal framework. In a semiclassical regime, the Boltzmann transport equation may connect the dispersion and collision kernel to conductivity. First verify the nondegenerate condition $n\ll N_c(T)$, where $N_c$ is the effective conduction-band density of states. For each populated valley $\nu$ and principal direction $i$, define a thermal wave number $k_{\mathrm{th},\nu i}=\sqrt{2m^*_{\nu i}k_{\mathrm B}T}/\hbar$ and a directional mean free path $\ell_{\mathrm{mfp},\nu i}=v_{\mathrm{th},\nu i}\tau_{\mathrm{tr},\nu i}$. Check $k_{\mathrm{th},\nu i}\ell_{\mathrm{mfp},\nu i}\gg1$, quasiparticle width small compared with the resolved band scale, sample dimensions large compared with every relevant mean free path, the weak-field condition above, and stability against the collision-memory approximation. Failure of these tests routes to coherent, quantum-kinetic, or finite-geometry transport. Numerical convergence of the solver does not validate those physical assumptions.

7. Measurement map. Convert the calculated conductivity tensor into the voltage-to-current ratio for the declared four-probe bar geometry. Verify the geometric factor by varying voltage-lead spacing or sample dimensions, confirm linearity under current reversal, bound contact injection and inhomogeneity, calibrate thermometry and carrier density independently, and demonstrate negligible self-heating. How Quantum Matter Is Measured and Data Interpretation and Pitfalls own that final inference layer.

The completed ledger therefore records the electron–ion model and prepared sample as input; the multivalley carrier distribution as the working state; valleys, anisotropic masses, phonons, impurities, and geometry as retained data; core states and remote bands as eliminated data; calculated band curvatures and matched collision parameters as distinct provenance classes; model, discretization, scattering, and measurement uncertainties separately; and four-probe low-field bulk conductivity as the licensed output. It stops only after the combined uncertainty is below $10\%$ across the declared window, including held-out temperature and density points, and plausible omitted channels have been varied or bounded.

Changing the target reopens the ledger. Indirect optical absorption requires photon and phonon matrix elements, excitonic effects, and optical selection rules that a low-field mobility model can omit. Ultrafast excitation may invalidate adiabatic separation and equilibrium occupations. A nanostructure may require explicit surfaces, valleys, confinement, and contacts. “The silicon Hamiltonian” is therefore not one permanent object.

Why Simplified Hamiltonians Are Powerful

For $N_{\mathrm e}$ electrons in $M$ spin-orbitals, the fixed-particle sector has dimension

$$\dim\mathcal H_{N_{\mathrm e}}=\binom{M}{N_{\mathrm e}}.$$

At half filling with $M=40$, $\binom{40}{20}=137{,}846{,}528{,}820$. This count explains why brute-force state storage rapidly fails, but it does not prove that every material problem is computationally hopeless. Symmetry sectors, locality, sparsity, perturbative scales, tensor structure, stochastic estimators, and observable-specific reductions can make structured questions tractable.

A simplified Hamiltonian earns trust by preserving the structure needed for the target while making the discarded structure testable. The strongest ledger therefore separates three axes:

• model validity: does the retained Hamiltonian represent the material and regime?
• physical convergence: are basis, size, boundary, cutoff, and approximation sequences controlled?
• solver convergence: has the declared finite problem actually been solved to the claimed tolerance?

Agreement on one benchmark validates only that benchmark contract. A ground-state energy match does not automatically certify spectral weight, transport, topology, or a phase assignment.

When the Ledger Must Be Reopened

Revisit an earlier arrow whenever the target crosses its validity boundary. Common triggers include:

• electronic excitation energies approaching eliminated-band gaps;
• interaction matrix elements comparable with the separation to remote states;
• strong nonadiabatic nuclear motion or ultrafast driving;
• disorder, surfaces, interfaces, or strain that destroy the assumed translation symmetry;
• temperatures or fields that populate new valleys, bands, phonon branches, or phases;
• probe resolution that exposes operator dressing or matrix elements omitted from the working model;
• a fitted parameter being transferred to a different basis, environment, or frequency window;
• disagreement among observables that should be described by the same retained theory.

The response is not automatically to restore every microscopic variable. Often the correct repair is a new effective theory with a larger retained space, frequency-dependent interactions, dressed observables, or a different state and solver.

Canonical Boundaries

This page owns the end-to-end provenance audit. It does not rederive the pieces it connects.

• Quantum Matter owns the broad scientific synthesis, and the Quantum Matter Map owns the reusable system–Hamiltonian–state–observable–resolution graph.
• Choosing a Model for Quantum Matter owns adequacy comparisons among candidate retained descriptions; this page audits the provenance and validity of a selected material pipeline.
• Born–Oppenheimer Scale Separation owns the electronic–nuclear derivation and nonadiabatic corrections.
• Bloch’s Theorem and the Band Theory Overview own periodic-state and band-interpretation details.
• Random Phase Approximation, Phonons, and Effective Hamiltonians in Quantum Matter own their respective derivations.
• Emergence and Effective Degrees of Freedom owns the model-independent matching theory. Later Quantum Matter chapters own phases, quasiparticles, topology, materials phenomenology, and experiments.

Common Pitfalls

• “Pauli exclusion is the electron interaction.” Antisymmetry changes allowed states and exchange structure; Coulomb interaction is a distinct term.
• “Born–Oppenheimer means frozen classical nuclei.” It organizes coupled electronic and nuclear motion; vibrations and nonadiabatic corrections remain.
• “Bloch bands require noninteracting electrons.” Translation symmetry also labels interacting many-body sectors. Independent-particle bands require additional reduction.
• “Self-consistent bands are measured excitation energies.” Their interpretation depends on the auxiliary theory and observable.
• “Screening is one dielectric constant.” It can depend on momentum, frequency, environment, matrix channel, and retained subspace.
• “First principles means assumption-free.” Composition, boundary conditions, relativistic scope, effective cores, functionals, cutoffs, and solvers remain declared choices.
• “A converged model calculation proves the material model.” Solver convergence, physical convergence, and model validity are different questions.
• “One Hamiltonian predicts every experiment.” Effective theories are state-, scale-, and observable-dependent.

Exercises

1. Classify the arrows

Classify each operation: (a) a unitary rotation within a complete retained band subspace; (b) Born–Oppenheimer separation; (c) an exact energy-dependent projection retaining every induced operator; (d) dropping all but onsite terms after that projection; (e) fitting a static $U$ to one measured spectrum; and (f) convolving a calculated response with a calibrated resolution function.

Solution

(a) is an exact representation change, although localization can be topologically obstructed. (b) is a controlled approximation only in a declared adiabatic regime. (c) is a formal exact reorganization within the original model; it becomes approximate when energy dependence or induced operators are truncated. (d) is a model reduction whose control requires range and larger-space tests. (e) is phenomenological matching and licenses transfer only after independent validation. (f) is part of the probe map: exact for the declared calibrated response model, but subject to calibration and background uncertainty.

2. Audit the fixed-ion step

Starting from $H_{\mathrm{micro}}$, identify which terms belong to the conditional electronic Hamiltonian at fixed $\mathbf R$. Why must $V_{\mathrm{II}}$ be restored in the adiabatic surface used to find equilibrium positions?

Solution

$T_{\mathrm I}$ acts on nuclear motion and is omitted from the conditional electronic eigenproblem. The electronic Hamiltonian retains $T_{\mathrm e}+V_{\mathrm{eI}}+V_{\mathrm{ee}}$ together with declared relativistic and external electronic terms. The ion–ion energy is constant with respect to electronic coordinates, but varies with $\mathbf R$; it must be added to $E_\alpha(\mathbf R)$ to obtain correct forces and equilibrium geometry. Any retained relativistic correction or external potential with explicit nuclear-coordinate dependence must likewise be included on the adiabatic surface.

3. Translation symmetry without independent electrons

Explain why a translation-invariant electron–electron interaction preserves total crystal momentum but does not imply a product of one-electron Bloch states.

Solution

A simultaneous lattice translation of every electron leaves coordinate differences and a periodic ionic background invariant, so the many-body Hamiltonian commutes with the total translation operator. Its eigenstates can therefore be labeled by a total crystal momentum defined modulo a reciprocal-lattice vector. The interaction still couples individual momenta subject to total-momentum conservation and generally entangles the electrons; symmetry labeling is not factorization.

4. Count a fermionic sector

Derive the dimension of the $N$-electron sector built from $M$ spin-orbitals and evaluate it for $N=20$, $M=40$. What does this count fail to establish?

Solution

An antisymmetric occupation basis is specified by choosing which $N$ of the $M$ spin-orbitals are occupied, giving $\binom{M}{N}$. Thus $\binom{40}{20}=137{,}846{,}528{,}820$. The count establishes rapid generic growth, but not the cost of every observable or algorithm. Symmetry, sparsity, scale separation, tensor structure, and approximations can substantially reduce a structured problem.

5. Change the target observable

The silicon ledger was built for low-field conductivity. What must be restored or re-matched to predict absorption near an indirect band edge?

Solution

The target now needs electromagnetic coupling and optical matrix elements, electron–hole interactions when excitonic effects matter, phonon modes and electron–phonon vertices to supply the required crystal momentum, finite-temperature occupations, linewidth mechanisms, and the spectrometer’s polarization and resolution. A relaxation time fitted to dc mobility is not automatically transferable to the optical process.

6. Audit a single-band truncation

Suppose $E_{\mathrm{target}}\ll\Delta_{\mathrm{remote}}$, but an interaction scale satisfies $U\sim\Delta_{\mathrm{remote}}$. Is a naive single-band model controlled?

Solution

Not from the target energy alone. The interaction can mix the retained band with remote states, renormalize parameters and observables, or generate additional operators. One must compare coupling matrix elements with the separation, perform a larger-space benchmark or controlled downfolding, and test whether the single-band predictions remain stable.

7. Audit a correlated-material spectrum

A workflow chooses a narrow $d$-orbital window, imports a static interaction calculated in a different orbital basis, adds a self-energy correction, solves the finite model to tight numerical tolerance, and identifies its raw spectral function with an angle-resolved photoemission spectrum. List at least four ledger defects before accepting the material claim.

Solution

First, the active-space and orbital gauges differ, so the imported interaction must be transformed or recomputed. Second, the polarization excluded from that interaction must be reconciled with the self-energy and solver to avoid gaps or double counting. Third, tight solver convergence certifies only the declared finite model; basis, frequency, size, analytic-continuation, and model uncertainties remain. Fourth, photoemission intensity contains photon-energy- and polarization-dependent matrix elements, occupations, escape-depth effects, backgrounds, and instrumental resolution; it is not the bare spectral function. Competing orbital windows, interaction frequencies, surface reconstructions, and alternative self-energies should be varied before a unique material interpretation is claimed.

References

• F. Aryasetiawan, M. Imada, A. Georges, G. Kotliar, S. Biermann, and A. I. Lichtenstein, “Frequency-Dependent Local Interactions and Low-Energy Effective Models from Electronic Structure Calculations,” Physical Review B 70, 195104 (2004), doi:10.1103/PhysRevB.70.195104.
• N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976.
• M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Clarendon Press, 1954.
• M. Born and R. Oppenheimer, “Zur Quantentheorie der Molekeln,” Annalen der Physik 389, 457–484 (1927), doi:10.1002/andp.19273892002.
• M. L. Cohen and J. R. Chelikowsky, Electronic Structure and Optical Properties of Semiconductors, 2nd ed., Springer, 1989.
• P. Coleman, Introduction to Many-Body Physics, Cambridge University Press, 2015.
• G. F. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, Cambridge University Press, 2005, doi:10.1017/CBO9780511619915.
• P. Hohenberg and W. Kohn, “Inhomogeneous Electron Gas,” Physical Review 136, B864–B871 (1964), doi:10.1103/PhysRev.136.B864.
• W. Kohn and L. J. Sham, “Self-Consistent Equations Including Exchange and Correlation Effects,” Physical Review 140, A1133–A1138 (1965), doi:10.1103/PhysRev.140.A1133.
• R. M. Martin, Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press, 2004.
• P. Y. Yu and M. Cardona, Fundamentals of Semiconductors, 4th ed., Springer, 2010, doi:10.1007/978-3-642-00710-1.
• J. M. Ziman, Electrons and Phonons, Oxford University Press, 1960.