Electronic Structure Methods Map

Electronic-structure methods approximate different parts of the same interacting-electron problem. Hartree–Fock optimizes one determinant, configuration interaction mixes determinants, coupled cluster organizes excitations exponentially, and density-functional theory replaces the many-electron wavefunction by a density-based formulation. Choosing among them requires knowing which features of the state and which observable matter.

A method that gives a useful equilibrium geometry may be inadequate for bond breaking, a diffuse excited state, or a small spin splitting. This guide connects the assumptions of each method with concrete comparisons that reveal its limitations.

Required background. Electronic structure introduces the ansätze and orbital bases; the molecular Hamiltonian defines the fixed-nuclei approximation. Helpful background. The Hartree–Fock example shows a self-consistency calculation.

Electronic states and the property being calculated

Start with a physical question: an equilibrium bond length, a vertical excitation at fixed nuclei, an adiabatic energy after geometry relaxation, or a response to a pulse. These require different calculations even for the same molecule.

Specify the charge, electronic symmetry, geometry, Hamiltonian, orbital basis, and property definition. For example, “a gas-phase vertical singlet excitation at the ground-state equilibrium geometry” distinguishes the electronic gap from a vibrational band maximum or an excitation after structural relaxation.

For an energy difference $\Delta E=E_B-E_A$, use balanced approximations for both states. A larger basis or more flexible correlation treatment for only one side can worsen the difference even if it improves that state’s absolute energy. In quantum chemistry, a method combined with its basis and auxiliary approximations is often called a model chemistry.

Method-Family Overview

The table summarizes characteristic uses and limitations. Helgaker, Jørgensen, and Olsen develop the molecular methods in detail.

Family	Basic object	Often useful for	Characteristic risk	Essential check
HF	optimized determinant	reference orbitals, qualitative states, baseline energies	missing correlation; symmetry breaking	stability, spin, and orbital character
Kohn–Sham DFT	density represented by auxiliary orbitals	structures, densities, large systems, routine screening	functional-dependent bias and qualitative failures	functional class, grid, basis, and benchmark
CI	linear determinant expansion	finite-basis benchmarks, selected states, spectroscopy	combinatorial growth; truncated CI lacks size extensivity	excitation-space and basis convergence
CC	exponential excitation ansatz	high-accuracy single-reference energies and properties	breakdown near degeneracy; nonvariational truncation	reference diagnostics and hierarchy comparison
multireference	optimized active-space expansion	bond breaking, near-degenerate states, intersections	active-space and state-averaging dependence	occupations, active-space enlargement, root tracking
QMC	sampled many-electron distribution	correlated ground-state benchmarks and large parallel workloads	nodal, time-step, population, and statistical bias	trial-state and sampling convergence
time-dependent or response	response vectors or propagated state/density	spectra, polarizabilities, driven dynamics	state-character and kernel/truncation failures	sum rules, time/frequency convergence, cross-method tests

The phrase “often useful” is intentionally weaker than “reliable.” Each row contains several distinct methods, implementations, and approximation levels.

Hartree–Fock

Choose HF when reference orbitals, a qualitative determinant description, or a baseline for correlation methods are needed. The Hartree–Fock construction gives the ansatz and variational bound; the helium calculation implements its self-consistency loop.

Converge the SCF equations, then test stability, spin, and orbital character. A stationary solution need not be the lowest determinant in its class. Restricted, unrestricted, and generalized determinants impose different constraints; a lower unrestricted energy may reflect broken spin symmetry rather than a satisfactory description of the target state. The total-energy upper bound concerns the same Hamiltonian and does not certify energy differences or other properties.

HF treats exchange within one determinant but omits correlation beyond that reference. Exchange and correlation explains the energy distinction; correlation energy is not itself an observable. For bond breaking, follow the changing configuration character rather than judging adequacy from the equilibrium energy alone. The H$_2$ decision case below illustrates this failure of a single reference.

Density-Functional Theory

Choose a named functional for a declared property and physical regime. The density-functional framework explains the density and orbital constructions. A reproducible calculation also specifies its basis or grid, spin treatment, integration settings, dispersion model if used, and self-consistent solution.

Local, semilocal, meta-GGA, hybrid, range-separated hybrid, and double-hybrid approximations make different tradeoffs. Their conceptual hierarchy is not an accuracy ordering for every observable. Compare materially different functional classes after numerical convergence; agreement among closely related functionals can reflect shared bias. Cohen, Mori-Sánchez, and Yang discuss these limitations.

Match the comparison to the likely failure: delocalization and fractional-charge errors, static-correlation or fractional-spin errors, spin-state sensitivity, diffuse or charge-transfer excitations, and missing long-range dispersion require different checks. Also test density sensitivity, meta-GGA integration grids, and unintended spin or orbital solutions. A dispersion correction can repair a missing asymptotic interaction without repairing density or static-correlation errors.

Configuration Interaction

Use CI for transparent state analysis, spectroscopic models, active-space solvers, or small-system benchmarks. Configuration interaction develops the determinant expansion, finite-space exactness, and excitation truncations. FCI is exact only for the declared orbital space, Hamiltonian, active-electron or frozen-core choice, and symmetry sector; it does not remove basis or physical-model errors.

Excitation-rank-truncated CI is variational in its retained linear space but generally lacks the disconnected products needed for size consistency and extensivity. The many-fragment exercise quantifies how the defect grows in a CISD-like calculation. Check separated-fragment limits before transferring a small-system error estimate to a larger system.

Selected CI adds determinants adaptively rather than by a fixed excitation rank. Record convergence against selection thresholds, perturbative corrections, and extrapolation choices; one selected expansion is not evidence that the FCI limit has been reached.

Coupled Cluster

The coupled-cluster construction explains the exponential ansatz, disconnected products, and projected equations. For a stable single-reference state, compare an excitation hierarchy: CCD retains doubles, CCSD adds singles, CCSD(T) estimates connected triples perturbatively, and iterative triples or higher corrections provide more demanding benchmarks. CCSD(T)‘s success near many equilibrium geometries does not establish its accuracy for every molecule or property.

A small amplitude residual establishes numerical convergence of the selected equations. It does not validate the truncated state model, and the projected energy has no general variational upper bound. Near-degenerate occupations, small denominators, strongly stretched bonds, double-excitation character, and interacting electronic surfaces call for the state diagnostics below and a higher-rank or multireference comparison. The H$_2$ decision case explains the special two-electron CCSD limit.

CC energies, gradients, response properties, and transition moments use different left and right states or Lagrangian machinery. Specify whether the reported property uses relaxed response, a finite-field derivative, or a transition matrix element. Converging an energy does not validate a polarizability, gradient, or EOM-CC transition property; Bartlett and Musiał review these distinctions.

Multireference Methods

Choose a multireference treatment when several configurations are essential across the states and geometries being compared. The active-space construction defines inactive, active, and external orbitals and the joint optimization in CASSCF. CAS$(n,m)$ specifies $n$ active electrons in $m$ active spatial orbitals; reproducibility additionally requires the orbital identities, included states, and state-averaging protocol.

Build the space from the process: breaking bonds and their antibonding partners, metal $d$ or radical orbitals, participating lone pairs and acceptors, and valence or Rydberg partners needed to balance several states. Use occupations and orbital character across the geometry path. An energy window alone can lose the intended orbitals when their ordering changes.

CASSCF captures the configurational rearrangements inside its active space while usually leaving substantial dynamic correlation outside it. Compare enlarged spaces and appropriate subsequent treatments, such as CASPT2, NEVPT2, multireference CI or CC, larger selected-CI or DMRG solvers, and tailored or embedding approaches. These differ in intruder-state behavior, contraction, size consistency, Hamiltonian partition, and state treatment; they are not interchangeable corrections.

State-specific orbitals may favor one root. State averaging instead optimizes a weighted mean of the included energies, with nonnegative weights summing to one; both the weights and root selection affect the orbitals. Along a path, track symmetry, state overlaps, transition densities, dominant configurations, occupations, and relevant observables. Root number alone is not a state identity.

Quantum Monte Carlo

Choose variational Monte Carlo to optimize Jastrow or multideterminant trial states, explore explicit electronic correlation, and estimate their energies or other observables. That treatment derives local-energy sampling and its variational assumptions. Trial-state flexibility and optimization control accuracy; autocorrelation and effective sample size control sampling precision. The variational bound applies to the exact trial-state expectation, so a finite noisy estimate can fall below the exact ground energy.

Choose diffusion Monte Carlo for imaginary-time projection while separately testing the fermionic nodal constraint. Under the usual fixed-node assumptions for a local real Hamiltonian, the projected energy is an upper bound; nonlocal pseudopotentials and their localization treatments need their own variational analysis. Validate the estimator for each non-energy property, because operators that do not commute with the Hamiltonian can have mixed-estimator bias.

Vary production length, trial-state optimization, time step, walker population, nodal quality, finite size and boundary conditions, and pseudopotential treatment independently. Population control can correlate walkers and bias results. Use blocking or another autocorrelation-aware analysis of the effective sample count: halving a standard error requires roughly four times as many effective samples. Neither a small error bar nor longer sampling removes fixed-node, localization, or other systematic bias. Foulkes et al. review these error channels.

Time-Dependent and Excited-State Methods

Choose among frequency-domain linear response, equation-of-motion or state-interaction calculations, and real-time propagation according to the observable. These are related constructions, not interchangeable ways of calculating every spectrum. The excited-state overview distinguishes energy definitions, state character, and transition intensities.

Linear response and state character

Response poles locate excitation energies and their residues encode transition strengths. TDHF, random-phase approximations, TDDFT, and coupled-cluster response differ in their reference and response kernel. In practical molecular TDDFT, common adiabatic kernels can describe many valence single excitations but lack the frequency structure required for generic double-excitation poles. Test long-range charge transfer, Rydberg diffuseness, functional-dependent state ordering, and behavior near conical intersections explicitly. The Runge–Gross density–potential mapping for a fixed initial state, under its regularity assumptions and up to a purely time-dependent scalar potential, does not establish the accuracy of an approximate kernel.

Equation-of-motion sectors

EOM-CC diagonalizes the similarity-transformed Hamiltonian in a retained excitation-operator manifold. The unprojected similarity transform preserves the spectrum; projection and truncation introduce the approximation. Specify the sector: neutral excitation, ionization, electron attachment, or spin flip. An eigenvalue of the projected transformed Hamiltonian is a total electronic energy in its sector; a neutral excitation energy is its difference from the CC ground-state energy. An explicitly energy-shifted implementation reports that difference directly.

Reference quality and the retained manifold determine which state characters are balanced. EOM-CCSD contains singles and doubles, but a predominantly double-excited state generally needs higher ranks for balanced relaxation and correlation when more than two electrons are correlated. Check transition moments with the corresponding left/right-state treatment, as developed by Stanton and Bartlett, and track states by character rather than energy order.

Real-time propagation

Real-time propagation can treat broadband response and nonlinear driving. Converge the time step, total propagation time, spatial representation, absorbing boundaries, field envelope and gauge, and conservation or norm drift. Fourier resolution scales inversely with total time, while the highest reliably represented frequency depends on the time step and propagator. Converged ground-state energies alone do not test these propagation and spectral limits.

Computational cost and orbital count

Let $N_o$ and $N_v$ denote occupied and virtual orbital counts and let $N$ denote a rough overall orbital scale. Representative dense canonical algorithms have the following leading operations:

Method	Representative leading work	Common shorthand	Important qualifier
HF	four-index integral handling plus SCF diagonalization	up to $N^4$ conventionally	direct, density-fitted, local, and grid methods change cost
Kohn–Sham DFT	grid evaluation plus orbital solution	often $N^3$–$N^4$	grid, exact exchange, periodicity, and sparsity dominate
MP2	$N_o^2N_v^3$	$N^5$	integral transformation and storage matter
CISD	roughly $N_o^2N_v^4$	$N^6$	determinant-space dimensions and symmetry alter prefactors
CCSD	$N_o^2N_v^4$	$N^6$	memory and communication can be limiting
CCSD(T) triples step	$N_o^3N_v^4$	$N^7$	variants and implementations differ
FCI	combinatorial determinant space	exponential	feasible size depends strongly on symmetry and solver
CASSCF	active-space CI plus orbital optimization	exponential in active space	inactive/external orbital work also matters
QMC	samples times per-sample evaluation	no single universal power	variance, nodes, autocorrelation, and target error control cost

These estimates apply to the stated algorithmic forms; sparsity and local approximations can change them.

Why the same scaling can behave differently

Wall time depends on:

• prefactors and tensor sparsity;
• number of electrons and basis functions;
• occupied-to-virtual ratio;
• molecular symmetry;
• integral approximation;
• memory, storage, and communication;
• parallel efficiency;
• convergence iterations;
• requested number of states or properties;
• and software implementation.

A nominally lower-scaling method can be slower for a particular system, and a local approximation can change both cost and error.

Method Choice by Physical Regime

The following entries are starting hypotheses. Each requires the final column.

Target regime	Plausible starting families	Main danger	Useful comparison
closed-shell equilibrium structure	DFT or HF-based correlation	hidden conformer, basis, or functional bias	basis/grid test and higher-level points
single-reference thermochemistry	CC hierarchy or calibrated composite protocol	basis and higher-excitation imbalance	basis extrapolation and hierarchy increments
bond breaking	CASSCF plus dynamic correlation; selected CI	changing active orbitals and unbalanced fragments	occupations and dissociation limits
weak intermolecular binding	dispersion-aware DFT or correlated wavefunction method	BSSE and geometry sensitivity	counterpoise or CBS analysis and benchmark
transition-metal spin gap	several DFT classes and multireference checks	functional, spin, active-space, and relativistic sensitivity	multiple references and state diagnostics
valence excitation	TDDFT, algebraic-diagrammatic construction (ADC), EOM-CC, or multireference response	state character and root switching	basis, functional/rank, and benchmark states
Rydberg excitation	diffuse-basis response or EOM method	missing diffuse space and continuum mixing	augmentation sequence and orbital extent
charge-transfer excitation	range-aware TDDFT or wavefunction method	asymptotic-potential error	donor–acceptor separation trend and benchmark
double excitation	multireference or higher EOM/response treatment	inadequate singles-doubles response space	configuration analysis and higher method
heavy-element property	relativistic DFT or wavefunction method	Hamiltonian and picture-change error	relativistic-level and core-space tests
large periodic ground state	periodic DFT or QMC	finite size, $k$ mesh, functional, pseudopotential	cell/$k$ convergence and independent benchmark

The table does not imply that every listed implementation is available or equally mature for gradients, excited states, spin–orbit coupling, or periodic boundary conditions.

For a periodic crystal whose chosen family is an independent-particle or Kohn–Sham method, Band Structure Workflows explains crystal structure, Brillouin-zone sampling, magnetic states, and convergence.

Interpreting state diagnostics

No scalar diagnostic universally partitions electronic states into “single-reference” and “multireference” classes.

Orbital gaps

A small HF or Kohn–Sham gap can warn of near-degeneracy or SCF instability, but:

• the value is method and orbital dependent;
• a small gap need not imply strong static correlation;
• a large gap does not exclude a poorly described excited state;
• and a Kohn–Sham gap is not generally a many-body excitation gap.

Natural occupations

For a one-particle reduced density matrix,

$$\gamma\phi_p = n_p\phi_p,$$

the eigenvalues $n_p$ are natural occupations. For the spin-summed spatial density matrix they lie between $0$ and $2$. Departures from $2$ and $0$ can reveal correlation relative to a closed-shell reference, but singly occupied orbitals are expected in ordinary open-shell states. For the full spin-orbital density matrix, a pure Slater determinant has only $0$ and $1$ occupations; this is the sharper determinant test. Spin summation can hide that distinction.

Occupations depend on the approximate state used to compute $\gamma$. They are more informative when followed along the relevant geometry path and compared across methods or active spaces.

Spin contamination

With spin operators measured in units of $\hbar$, for a nominal spin $S$,

$$\langle \hat S^2\rangle = S(S+1)$$

for a pure spin eigenstate. A larger unrestricted value signals contamination by other spin sectors. A small deviation does not by itself validate the energy, and a broken-symmetry state can sometimes be a useful intermediate model if interpreted explicitly.

Coupled-cluster amplitude diagnostics

$T_1$, $D_1$, largest amplitudes, and perturbative-triples fractions probe different features of a chosen CC solution. Their thresholds are empirical and system dependent.

Use them as prompts for:

• orbital stability checks;
• natural-occupation analysis;
• comparison to higher excitation rank;
• a multireference calculation;
• or a smaller-basis FCI/selected-CI benchmark.

DFT functional spread

Spread across functionals is useful only when the set spans meaningfully different approximations. It is not a calibrated uncertainty distribution, and cherry-picking a functional after seeing the benchmark leaks reference information into the prediction.

State overlap

For geometry steps $R_k$ and $R_{k+1}$, an overlap measure

$$\left| \langle \Psi_I(R_k) | \Psi_J(R_{k+1}) \rangle \right|$$

can assist root tracking. In nonorthogonal orbital spaces, the overlap must be evaluated with the proper metric and orbital transformation.

State overlap should be combined with symmetry, transition properties, and configuration character near avoided or true crossings.

Convergence Is Multidimensional

For a scalar result $Q$, consider a two-dimensional method–basis grid:

	basis $B_1$	basis $B_2$	basis $B_3$
method $M_1$	$Q_{11}$	$Q_{12}$	$Q_{13}$
method $M_2$	$Q_{21}$	$Q_{22}$	$Q_{23}$
method $M_3$	$Q_{31}$	$Q_{32}$	$Q_{33}$

Define local increments

$$\delta_B Q = Q_{i,j+1} - Q_{ij}$$

and

$$\delta_M Q = Q_{i+1,j} - Q_{ij}.$$

If these increments change substantially across the grid, basis and method errors are coupled. A correction computed in a very small basis may not transfer cleanly to a larger basis.

A cross increment

A simple nonadditivity diagnostic is

$$\Delta_{\times}Q = Q_{22} - Q_{21} - Q_{12} + Q_{11}.$$

If $\Delta_{\times}Q$ is small relative to the target tolerance, an additive composite estimate may be plausible. If it is large, method and basis increments cannot be treated as independent without further evidence.

Energy differences require balanced convergence

For

$$\Delta E = E_B-E_A,$$

the relevant basis error is

$$\delta_B\Delta E = \delta_B E_B - \delta_B E_A.$$

Large absolute basis errors can cancel, or small-looking component errors can reinforce. Converge the reported difference directly.

Geometry-dependent balance

A method adequate near equilibrium may fail at dissociation or along a reaction coordinate. Convergence should be sampled across:

• minima;
• transition regions;
• asymptotic fragments;
• avoided crossings;
• and any geometry where state character changes.

One converged stationary point does not validate an entire potential-energy surface.

Sources of uncertainty

A conceptual decomposition is

$$\Delta Q = \Delta Q_{\mathrm{model}} + \Delta Q_{\mathrm{repr}} + \Delta Q_{\mathrm{method}} + \Delta Q_{\mathrm{solver}} + \Delta Q_{\mathrm{sampling}} + \Delta Q_{\mathrm{property}}.$$

The terms denote:

Term	Examples
$\Delta Q_{\mathrm{model}}$	relativity, nuclear motion, environment, finite cell, pseudopotential
$\Delta Q_{\mathrm{repr}}$	orbital basis, grid, plane-wave cutoff, auxiliary basis
$\Delta Q_{\mathrm{method}}$	functional approximation, excitation truncation, active space, fixed nodes
$\Delta Q_{\mathrm{solver}}$	incomplete SCF, eigensolver, propagation, optimizer convergence
$\Delta Q_{\mathrm{sampling}}$	Monte Carlo variance, autocorrelation, finite population
$\Delta Q_{\mathrm{property}}$	finite field, numerical derivative, response truncation, state assignment

This is a schematic decomposition, not a theorem of additivity. Error channels interact. For example, basis incompleteness can change CC amplitudes, active-space occupations, DFT density error, or QMC nodes.

Numerical precision versus model accuracy

Suppose an SCF energy is converged to $10^{-10}E_{\mathrm h}$. That number sets a solver tolerance. It does not establish:

• basis convergence to $10^{-10}E_{\mathrm h}$;
• correlation accuracy to that level;
• physical-model completeness;
• or agreement with experiment to ten decimal places.

Report enough digits to reproduce arithmetic, but round interpreted claims to the supported uncertainty.

Statistical and systematic uncertainty

For QMC,

$$Q = \bar Q \pm \sigma_{\bar Q}$$

reports sampling precision under the estimator assumptions. Fixed-node, time-step, finite-size, and pseudopotential biases are systematic and require separate tests.

For deterministic electronic structure, a spread across basis sets or methods is not automatically a confidence interval. It measures sensitivity to those particular changes.

Worked Decision Cases

Case 1: H₂ dissociation

At large separation, bonding and antibonding configurations must both contribute to the singlet. Restricted HF keeps unwanted ionic components; unrestricted HF can lower the energy by breaking spin symmetry. A CAS$(2,2)$ state includes both occupations and gives the correct qualitative neutral-fragment limit, although larger orbital spaces are needed for quantitative energies away from that limit.

H$_2$ is an important exception to generic warnings about truncated CC: with only two electrons, singles and doubles exhaust the excitation ranks. CCSD can therefore reproduce FCI in the same orbital space when the desired root is found and the reference has nonzero overlap with it. Perturbative triples add nothing for two electrons. Numerical conditioning or root selection can still cause difficulty. To study genuine excitation-rank failure during bond breaking, use a system with more correlated electrons, such as simultaneous breaking of several bonds.

Case 2: Weak dimer binding

Claim: Determine a gas-phase interaction energy at a declared geometry.

Risks: Dispersion, basis-set superposition, geometry sensitivity, and small differences of large total energies.

Define the interaction energy as

$$\Delta E_{\mathrm{int}} = E_{AB} - E_A - E_B$$

with geometry and monomer deformation conventions stated.

Reasonable map:

1. Use a dispersion-capable DFT or correlated wavefunction method.
2. Add diffuse and polarization functions where needed.
3. Test basis convergence and Boys–Bernardi counterpoise sensitivity;
4. compare a smaller model or selected geometries to a high-level correlation benchmark;
5. report deformation, zero-point, and thermal corrections separately.

A tiny SCF residual is irrelevant if the interaction energy changes by more than the target tolerance under basis augmentation.

Case 3: Lowest bright valence excitation

Claim: Predict a vertical excitation and oscillator strength near a closed-shell equilibrium geometry.

Risks: State ordering, diffuse contamination, double-excitation character, and method-dependent transition moments.

Reasonable map:

1. Diagnose orbital and configuration character.
2. Test diffuse augmentation even for a nominally valence state.
3. Compare TDDFT functional classes or use EOM-CC/ADC where feasible.
4. inspect transition density and oscillator strength, not energy alone;
5. track the same state across basis and method changes;
6. benchmark a representative subset against a higher excitation rank or multireference treatment.

An energy match with a measured band maximum is not enough if vibronic structure, solvent shift, and vertical-to-band assignments are unresolved.

Case 4: Transition-metal spin gap

Claim: Determine the energy difference between two low-lying spin states.

Risks: Several near-degenerate $d$ occupations, functional sensitivity, state-specific orbital relaxation, scalar relativity, spin–orbit coupling, and environmental effects.

Reasonable map:

1. Search multiple orbital occupations and spin solutions.
2. inspect natural occupations and metal–ligand active orbitals;
3. compare materially different DFT approximations without treating their spread as a confidence interval;
4. construct and enlarge a balanced active space;
5. add dynamic correlation and relativistic corrections consistently;
6. test ligand geometry, solvation, and thermal terms if the comparison is experimental.

No single diagnostic threshold can certify this problem class.

Exercises

Diagnose size inconsistency

Consider $M$ identical noninteracting two-electron fragments in disjoint orthogonal orbital spaces. Each fragment retains a reference configuration $|0\rangle$ and a double excitation $|d\rangle$, with Hamiltonian

$$h=\begin{pmatrix}0&v\\v&\Delta\end{pmatrix}, \qquad \Delta>0,\quad v\in\mathbb R.$$

The two-fragment calculation shows why a global singles-and-doubles truncation loses simultaneous local double excitations. Here retain only the all-reference state $|\Phi_0\rangle$ and the $M$ states $|j\rangle$ with fragment $j$ doubly excited. No local single excitations occur in this model.

One isolated fragment has energy $e=(\Delta-\sqrt{\Delta^2+4v^2})/2$. Find the truncated energy for $M$ fragments and compare it with $Me$ as $M\to\infty$ at fixed $v\ne0$. Why does agreement through second order in $v$ at fixed $M$ fail to establish size extensivity?

Solution

The projected Hamiltonian has $\langle\Phi_0|H|j\rangle=v$ and $\langle i|H|j\rangle=\Delta\delta_{ij}$. Only the normalized symmetric state $|S\rangle=M^{-1/2}\sum_j|j\rangle$ couples to the reference. In the span of $|\Phi_0\rangle,|S\rangle$,

$$H_{\mathrm{bright}} =\begin{pmatrix}0&\sqrt M\,v\\\sqrt M\,v&\Delta\end{pmatrix}.$$

The $M-1$ orthogonal combinations remain at $\Delta$. Thus

$$E_{\mathrm{tr}}(M) =\frac{\Delta-\sqrt{\Delta^2+4Mv^2}}2, \qquad E_{\mathrm{exact}}(M)=Me.$$

At fixed nonzero $v$ and fixed $\Delta$,

$$E_{\mathrm{tr}}(M) =-|v|\sqrt M+\frac{\Delta}{2}+O(M^{-1/2}), \qquad \frac{E_{\mathrm{tr}}(M)}{M}\longrightarrow0.$$

The exact energy per fragment instead remains $e<0$. Although the retained space grows with $M$, it allows only one fragment to be excited at a time; it cannot describe independent correlation on every fragment.

For fixed $M$ and $Mv^2/\Delta^2\ll1$, expansion gives

$$E_{\mathrm{tr}}(M)-Me =M(M-1)\frac{v^4}{\Delta^3} +O\!\left(M^3\frac{v^6}{\Delta^5}\right).$$

The second-order terms agree, but the expansion is not uniform as $M$ grows. One cannot extrapolate a small-fragment defect linearly or use that weak-coupling expansion at arbitrarily large $M$. Exponential factorization $e^{\sum_jT_j}=\prod_j e^{T_j}$ retains the missing products; it is exact in this two-state-per-fragment model, without making doubles-only coupled cluster exact for general fragments.

Test an additive composite estimate

For one energy difference in eV, a method–basis grid is

	$B_1$	$B_2$
$M_1$	2.10	2.24
$M_2$	2.31	2.52

Compute the cross increment

$$\Delta_{\times}Q = Q_{22}-Q_{21}-Q_{12}+Q_{11}.$$

Would a target tolerance of $0.02\,\mathrm{eV}$ support treating method and basis corrections as independent?

Solution

Substitution gives

$$\begin{aligned} \Delta_{\times}Q &= 2.52-2.31-2.24+2.10 \\ &= 0.07\,\mathrm{eV}. \end{aligned}$$

This exceeds the $0.02\,\mathrm{eV}$ target. The method change depends materially on the basis, so an additive correction assembled from the smallest calculations is not supported at the requested tolerance.

Interpret a coupled-cluster diagnostic

A calculation has a small CCSD residual, a $T_1$ diagnostic of $0.031$, a large perturbative triples contribution, and natural occupations $1.73$ and $0.27$ for two frontier spatial orbitals. What can and cannot be concluded?

Solution

The small residual shows that the chosen CCSD equations were solved numerically. It does not validate the single-reference truncation.

The elevated $T_1$, large triples correction, and noninteger frontier occupations are mutually consistent warnings that more than one configuration may matter. None is an individual proof that CCSD(T) fails, and no universal numerical cutoff settles the question.

A responsible response is to:

• check HF stability and orbital alternatives;
• compare CCSD, CCSD(T), and a higher excitation treatment in a smaller basis;
• construct a physically motivated active space;
• inspect behavior across the relevant geometry region;
• and compare the target energy difference, not only total energies.

Budget QMC sampling

A QMC estimate has a standard error of $0.40\,\mathrm{m}E_{\mathrm h}$ after a fixed effective sample count. If all else is unchanged, by what factor must the effective sample count increase to reduce the standard error to $0.10\,\mathrm{m}E_{\mathrm h}$?

Solution

The desired standard error is four times smaller. Since

$$\sigma_{\bar Q} \propto N_{\mathrm{eff}}^{-1/2},$$

the effective sample count must increase by

$$\left( \frac{0.40}{0.10} \right)^2 = 16.$$

This improves sampling precision only. It does not reduce fixed-node, time-step, finite-size, or pseudopotential bias.

Choose methods for two excited states

At one geometry, state $A$ is dominated by one valence single excitation. State $B$ is dominated by a simultaneous two-electron promotion. Explain why one method may not describe both states with balanced accuracy.

Solution

Common adiabatic TDDFT kernels describe single particle–hole excitations naturally but lack the frequency structure needed for generic double-excitation poles. EOM-CCSD includes explicit singles and doubles, so it can contain both states; however, predominantly double-excited states generally need triple and higher configurations for balanced relaxation and correlation.

For a system with more than two correlated electrons, balanced options for state $B$ can therefore require:

• higher EOM excitation rank;
• a multireference state-interaction method;
• selected CI;
• or another approach whose response space explicitly contains the relevant configurations.

The method should be selected for the hardest state needed by the claim, and both roots should be tracked by character rather than energy order.

References

• Bartlett, R. J., and M. Musiał. “Coupled-cluster theory in quantum chemistry.” Reviews of Modern Physics 79, 291–352 (2007). doi:10.1103/RevModPhys.79.291.
• Boys, S. F., and F. Bernardi. “The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors.” Molecular Physics 19, 553–566 (1970). doi:10.1080/00268977000101561.
• Cohen, A. J., P. Mori-Sánchez, and W. Yang. “Challenges for density functional theory.” Chemical Reviews 112, 289–320 (2012). doi:10.1021/cr200107z.
• Foulkes, W. M. C., L. Mitas, R. J. Needs, and G. Rajagopal. “Quantum Monte Carlo simulations of solids.” Reviews of Modern Physics 73, 33–83 (2001). doi:10.1103/RevModPhys.73.33.
• Helgaker, T., P. Jørgensen, and J. Olsen. Molecular Electronic-Structure Theory. Wiley (2000). doi:10.1002/9781119019572.
• Hohenberg, P., and W. Kohn. “Inhomogeneous electron gas.” Physical Review 136, B864–B871 (1964). doi:10.1103/PhysRev.136.B864.
• Kohn, W., and L. J. Sham. “Self-consistent equations including exchange and correlation effects.” Physical Review 140, A1133–A1138 (1965). doi:10.1103/PhysRev.140.A1133.
• Roos, B. O., P. R. Taylor, and P. E. M. Siegbahn. “A complete active space SCF method (CASSCF) using a density matrix formulated super-CI approach.” Chemical Physics 48, 157–173 (1980). doi:10.1016/0301-0104(80)80045-0.
• Runge, E., and E. K. U. Gross. “Density-functional theory for time-dependent systems.” Physical Review Letters 52, 997–1000 (1984). doi:10.1103/PhysRevLett.52.997.
• Stanton, J. F., and R. J. Bartlett. “The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties.” Journal of Chemical Physics 98, 7029–7039 (1993). doi:10.1063/1.464746.

Further Study

Electronic structure develops the ansätze in more detail. Potential-energy surfaces connects electronic energies to nuclear motion, while AMO benchmarks provides small calculations for checking implementations.