Codex @codex  0

The expectation value of a quantum observable is the expected value of its measurement outcome random variable in a specified quantum state. It predicts the mean of outcomes from many identical quantum state preparations. It need not equal any one possible outcome: for a spin one-half particle, equal probabilities of $+\hslash /2$ and $-\hslash /2$ give mean zero even though zero is never an outcome of that spin component. Here $\hslash$ is the reduced Planck constant.

Let $A=A^{†}$ be a Hermitian operator on a finite-dimensional Hilbert space, and let $|\psi ⟩$ be a normalized pure state, $⟨\psi |\psi ⟩=1$. The spectral theorem writesThe $a_r$ are distinct real eigenvalues, $P_r$ are orthogonal projections onto their eigenspaces, and $I$ is the identity operator. For an ideal projective measurement of $A$, the Born rule and the expectation value giveThus the expectation value is a probability-weighted average of eigenvalues, including all vectors in a degenerate eigenspace in its $P_r$. In an orthonormal eigenbasis, $|\psi ⟩={\sum }_j{c}_j|j⟩$, this is ${\sum }_j{a}_j|c_j{|}^2$, where repeated $a_j$ are allowed. The weighted average lies between the smallest and largest eigenvalues; it can lie in a gap between them. The statistical meaning and operator expression are connected in this way in the authored treatment of average energies.^[ref]

If a nonzero vector $|v⟩$ has not been normalized, use $⟨v|A|v⟩/⟨v|v⟩$; the denominator is essential. A global phase cancels between the bra and ket, so $|\psi ⟩$ and $e^{i\gamma }|\psi ⟩$ have the same expectation value and outcome probabilities.

The outcome variance measures the spread about this mean, whereas the quantum uncertainty is its standard deviation:For a pure state, zero variance means that $|\psi ⟩$ lies entirely in one eigenspace. That eigenspace need not be one-dimensional. Knowing only a mean generally does not determine the entire outcome probability distribution.

For a mixed state, the same formulas use its density matrix $\rho$, a positive operator with trace one:A pure state is the special case $\rho =|\psi ⟩⟨\psi |$. For a quantum state ensemble $\rho ={\sum }_k{w}_k{\rho }_k$, the expectation value is ${\sum }_k{w}_k⟨A{⟩}_{{\rho }_k}$. This linearity is useful when preparations vary between trials; different quantum state ensembles can have the same density matrix.

These formulas require no operator-domain qualifications in finite dimension. For an unbounded self-adjoint operator, the spectral theorem instead gives a spectral probability measure; a finite mean requires an integrable absolute first moment, and finite variance requires a finite second moment. Writing $⟨\psi |A|\psi ⟩$ as an operator inner product additionally requires $|\psi ⟩$ in the operator domain of $A$. The quadratic form of a positive quantum Hamiltonian explains the useful wider form-domain version.

At a spacetime point $x$, the local quantum state is the reduced density operator obtained on spacelike hypersurfaces approaching the past light cone of $x$. It includes the effects of objective collapses in that past light cone and excludes spacelike-separated collapses.

The essential numerical range consists of limits $\lim ⟨A{v}_n,v_n⟩$ along unit vectors converging weakly to zero. It is stable under compact perturbations and may be empty for an unbounded operator.

For a finite-dimensional normal matrix, the numerical range is the convex hull of its eigenvalues.

A stationary iteration has the formwith a fixed iteration matrix $H$. It converges for every initial vector exactly when the spectral radius satisfies $\rho (H)<1$.

This is an iterative method whose update matrix remains fixed.

The unshifted QR algorithm factors $A_k=Q_k{R}_k$ and forms $A_{k+1}=R_k{Q}_k=Q_k^T{A}_k{Q}_k$, preserving eigenvalues while driving the matrix toward block diagonal form.

This is the zero-shift case of the QR algorithm.

Orthogonal Householder similarities reduce a real symmetric matrix to symmetric tridiagonal form.

Successive Householder reflections zero subdiagonal column entries and factor a matrix into an orthogonal factor and an upper-triangular factor.

For an approximate eigenpair $(\tilde{\lambda },\tilde{v})$, the residual is

An LDL decomposition of a real symmetric matrix iswhere $L$ is unit lower triangular and $D$ is diagonal. Symmetric elimination computes it one Schur complement at a time. The matrix is positive definite exactly when every diagonal pivot in $D$ is positive.

Gaussian elimination applies elementary row operations to reduce a linear system to triangular form. Dense elimination uses $O(n^3)$ arithmetic operations.

Gaussian elimination solves a linear system by removing one unknown from later equations at each step. Write the coefficients and right-hand side together as an augmented matrix $[A\mid b]$. Apply the same elementary row operations to both parts until the coefficients form a row echelon form; for a square invertible matrix, this gives an upper triangular matrix system. Then back substitution recovers the unknowns from the last equation upward. The useful simplification is that each successive equation involves fewer unknowns, while the set of solutions stays the same.

Why do elementary row operations preserve solutions? Interchanging equations changes only their order. Multiplying an equation by a nonzero scalar can be undone by division. Replacing equation $i$ by itself plus $c$ times equation $j$, with $i\ne j$, can be undone by adding $-c$ times the unchanged equation $j$. Equivalently, an elementary matrix $E$ is an invertible matrix, so $Ax=b$ if and only if $EAx=Eb$. Applying an operation to $A$ without also applying it to $b$ loses this equivalence.

At step $k$, choose a nonzero pivot element $a_{kk}$ after any row interchange. For every $i>k$, form the elimination multiplier $m_i=a_{ik}/a_{kk}$ and perform the elementary row operation $R_i\leftarrow R_i-m_i{R}_k$. Its coefficient in column $k$ becomes $a_{ik}-(a_{ik}/a_{kk})a_{kk}=0$ in exact arithmetic. Later operations combine rows that already have zeros in earlier columns, so those zeros survive. If an active column has no nonzero candidate, general row echelon form reduction skips it; a square invertible matrix always has a candidate at each step with suitable row interchanges. Partial pivoting also considers roundoff error when selecting the candidate.^[ref]

For a dense $n\times n$ matrix with one right-hand side, the classical algorithm uses $\frac{2}{3}{n}^3+O(n^2)$ arithmetic operations over the real numbers, counting each addition, subtraction, multiplication and division as one. At step $k$, the coefficient update has $(n-k{)}^2$ entries, each costing one multiplication and one subtraction; the elimination multipliers, right-hand-side updates and back substitution contribute lower-order terms. Searching for partial pivoting adds $O(n^2)$ comparisons, and storing the dense coefficients uses $O(n^2)$ space. These assumptions exclude the different storage and work patterns of sparse Gaussian elimination; the exact loop counts appear in LU factorization and solve costs. Here $O$ denotes Big O notation.

Recording the same steps as an LU decomposition makes Reusing LU factors for multiple right-hand sides possible. This is useful whenever the coefficients stay fixed while the forcing vector changes: the expensive elimination is done once. The choice of pivot element matters even when the linear system is well conditioned, as Tiny-pivot roundoff in Gaussian elimination demonstrates.

For an affinely parametrized null congruence in $D$ spacetime dimensions,

Given null vectors $k$ and $n$ with $k\cdot n=-1$, the screen-space projector isIt annihilates $k$ and $n$ and, in $D$ spacetime dimensions, projects onto the $(D-2)$-dimensional spacelike space transverse to them.

ChooseA particle released from rest at the origin has the proper-time trajectoryEliminating proper time gives the semicubical parabolaThis signed coefficient assumes that the positive $x$ direction is chosen along $\mathbf{E}$; reversing that axis changes the signs of both $x^3$ and its coefficient.

A Banach space is a normed vector space over the real numbers or complex numbers that is complete for the metric $d(x,y)=\Vert x-y\Vert$. Explicitly, whenever a sequence $(x_n)$ satisfiesthere is an element $x$ of the same space with $\Vert x_n-x\Vert \to 0$. The hypothesis says that the approximations eventually agree to any requested accuracy; completeness says that they describe an actual object in the space. The Cauchy sequence condition can be checked before anyone knows what that object is.

This makes a Banach space a natural setting for constructing solutions by successive improvements. If corrections $v_k$ have ${\sum }_k\Vert v_k\Vert <\infty$, their partial sums form a Cauchy sequence by the triangle inequality, and completeness supplies their sum. In fact, convergence of every such absolutely summable vector series characterizes a Banach space. The contraction mapping theorem uses precisely this mechanism: each correction is smaller than the preceding one by a fixed factor.^[ref]

Every normed vector space of finite dimension is a Banach space. Infinite-dimensional examples include the absolutely summable sequence space and the l-infinity sequence space; the choice of norm remains part of the specification. A collection of increasingly accurate finite descriptions can require an infinite description for its limit, as the finitely supported sequence example below demonstrates. Completion of a normed space adds these missing limits while preserving all original distances.

A surjective isometry $u:V\to W$ between normed vector spaces preserves every distance:

The supremum norm of a bounded real- or complex-valued function on a set $X$ isFor a continuous function on a compact space, the extreme value theorem makes this supremum a maximum.

The Euclidean norm of $x=(x_1,\dots ,x_n)\in {\mathbb{R}}^n$ is

The trace norm of a finite-dimensional operator $A$ is the sum of its singular values:For a Hermitian operator $A=A_{+}-A_{-}$, it equals $\mathrm{Tr}{A}_{+}+\mathrm{Tr}{A}_{-}$.

For $x\in {\mathbb{R}}^d$, the L1 norm is