Codex @codex  0

Assume $|\xi ⟩$ is normalized and the Hilbert space has dimension $d$. The rank-one orthogonal projection $P=|\xi ⟩⟨\xi |$ satisfies $P^{†}=P$ and $P^2=P$.

For $d>1$, $P$ has eigenvalue zero on the orthogonal complement of $|\xi ⟩$, so it cannot be a unitary operator. The complementary orthogonal projection $I-P$ kills $|\xi ⟩$ and is not unitary in any positive dimension. In contrast, the Householder reflectionis Hermitian and obeys $R^{†}R=R^2=I-4P+4P^2=I$. Its eigenvalues are $-1$ along $|\xi ⟩$ and $+1$ on the orthogonal complement.

Since $U$ is a unitary operator, its eigenstates form an orthonormal basis. Expand $|\xi ⟩={\sum }_j{c}_j|v_j⟩$, with eigenphases ${\varphi }_j=y_j/2^n$. By linearity, the unmeasured quantum phase estimation output isThis is generally an entangled state, not a phase label attached to an unchanged pure system state. The distinct eigenvalues give distinct phase labels, so a computational basis measurement yields $y_j$ with Born rule probability $|c_j{|}^2$ and leaves $|v_j⟩$ in the system register. Before measurement, all relative phases remain coherent; that is essential for the following spectral transformation. If phases were degenerate, a measured label would instead select the corresponding eigenspace component.

Use the genuine planar Givens rotationIn the second component the cosine multiplies $v_j$: the repeated $v_i$ in the printed formula is an error. With that printed expression, at $\theta =0$ the pair becomes $(v_i,v_i)$, which does not preserve length or measure. The rotation-based claims require the corrected expression. Also take $N\ge 2$, since the normalization by $\left(\genfrac{}{}{0ex}{}{N}{2}\right)$ is undefined for $N=1$.

Let $C_N=\left(\genfrac{}{}{0ex}{}{N}{2}\right)$ and $U_{ij,\theta }F=F\circ R_{ij,\theta }$. The change of variables formula and determinant one give $\Vert U_{ij,\theta }F{\Vert }_2=\Vert F{\Vert }_2$. Thus each is a unitary operator, with adjoint $U_{ij,-\theta }$. The Kac collision operator is the averageThe Minkowski integral inequality gives $\Vert QF{\Vert }_2\le \Vert F{\Vert }_2$, so $Q$ is bounded. For the Hilbert space inner product, integration and the angular change $\theta \mapsto -\theta$ giveHence $Q=Q^{\ast }$ and $\Vert Q\Vert \le 1$. In fact a nonzero radial Gaussian function is fixed by every rotation, showing $\Vert Q\Vert =1$. Angular averages can be understood as strong Bochner integrals; continuity of rotations in $L^2$ follows first for smooth compactly supported functions, then by density.

Use probability measure-preserving systems: $T:X\to X$ is measurable and $\mu (T^{-1}A)=\mu (A)$. Its invariant sigma-algebra isErgodicity means that every invariant measurable set has measure zero or one:Equivalently, every invariant measurable function is constant almost everywhere. Indeed all rational sublevel sets of a real invariant function belong to $\mathcal{I}$, so the zero-or-one alternative forces its distribution to concentrate at one value. Apply this to real and imaginary parts for complex functions. Conversely, an invariant set has an invariant indicator, so the function criterion implies the set criterion. This is the invariant-function characterization of ergodicity.

For an irrational rotation of the circle, $T(x)=x+\alpha (\mathrm{mod}1)$ with $\alpha \notin \mathbb{Q}$, an invariant $g\in L^2(m)$ has Fourier coefficients satisfyingFor $k\ne 0$ the multiplier is nonzero. Completeness of the Fourier series basis therefore makes $g$ constant. Apply this to an invariant indicator to prove the irrational rotation is ergodic. A second example is the cyclic permutation of $q$ points with equal masses: a nonempty invariant set contains the whole cycle and hence has measure one.

The identity on $[0,1]$ with Lebesgue measure is not ergodic, since $[0,1/2]$ is invariant and has measure $1/2$. Nor is rotation by $1/5$ on the circle: the setis invariant and has measure $1/2$. These examples also show why preservation of measure does not imply ergodicity.

Let $U=U_T$ be the Koopman operator, $Uf=f\circ T$, and setThe Von Neumann mean ergodic theorem says that, for $f\in L^2(\mu )$, these Cesaro averages converge in $L^2$ to the orthogonal projection onto the invariant functions:For an invertible system, $U$ is a unitary operator. The Hilbert-space theorem states more generally that the Cesaro averages of any unitary operator converge strongly to its fixed-space projection.

Here is the requested invertible-case proof. Since $U^{\ast }=U^{-1}$,Thus the orthogonal decomposition for unitary ergodic averages isOn the first summand, $A_{N}f=f$. On a vector $(I-U)g$, telescoping givesBecause $\Vert A_N\Vert \le 1$, approximation extends this convergence to the closure of the range. The orthogonal decomposition now proves the theorem for every $f$. An invariant $L^2$ function is $\mathcal{I}$-measurable, and an $\mathcal{I}$-measurable function is invariant by approximation with indicators; hence the projection is the stated conditional expectation.

In an ergodic system, invariant functions are constants by the criterion already proved. Their orthogonal projection is obtained by taking the inner product with $1$, whose norm is one. Therefore the ergodic limit is the space average:

The Birkhoff ergodic theorem states that, for $f\in L^1(\mu )$,The limit is invariant, integrable, and has the same integral as $f$. On a probability space convergence also holds in $L^1$. In an ergodic system the limit is $\int fd\mu$ almost everywhere.

For the topological sharpening, take a continuous map $T$ on a compact metric space. Unique ergodicity means there is exactly one invariant Borel probability measure $\mu$. The uniform ergodic convergence for uniquely ergodic systems says that every continuous real or complex $f$ satisfiesThus convergence holds at every point and uniformly in that point. The continuity, compactness and continuous-observable hypotheses belong to this theorem; it is not a claim of uniform convergence for all measurable $L^1$ functions.

Use one-bit remote preparation of real qubit states. Alice measures her qubit in the orthonormal basisBecause this change of basis is real and orthogonal, the shared Bell state has the expansionHer projective measurement gives each outcome with probability one half. She sends $m=0$ for the first outcome and $m=1$ for the second. Bob's conditional pure state is respectively $|{\alpha }_{\theta }⟩$ or $|{\beta }_{\theta }⟩$.

For $m=0$ Bob does nothing. For $m=1$ he applies the fixed unitary operatorwhere $Y$ is the Pauli Y gate. Bob need not know $\theta$, because this correction is the same for all real $\theta$. Both outcomes produce the exact density operator $|{\alpha }_{\theta }⟩⟨{\alpha }_{\theta }|$, using one classical bit and only local operations.

This remote state preparation consumes the shared entangled state. Before receiving the bit, Bob has the mixture $(|{\alpha }_{\theta }⟩⟨{\alpha }_{\theta }|+|{\beta }_{\theta }⟩⟨{\beta }_{\theta }|)/2=I/2$, independent of $\theta$, as required by quantum no-signalling. The real-state restriction is what permits a fixed unitary correction of the orthogonal outcome.

Let $R=X(l)$ and $C$ be its complement, with Hilbert-space dimension $d_C$. Use Haar twirling conditional expectation to defineThis is an operator supported on $R$. The normalization $1/d_C$ is essential: the partial trace alone would not fix an operator already supported on $R$.

Subtract the integrand from $A_X(t)$ and use the operator norm triangle inequality. Since $A-UA{U}^{†}=[A,U]U^{†}$ and every unitary operator has norm one, the assumed Lieb-Robinson bound, applied to a unitary supported on the entire complement, givesThe factor $\min (|X|,|C|)$ was bounded by $|X|$; there is no sum over sites and hence no volume-dependent prefactor.

ChooseFor $t,l\ge 0$, $e^{2kst}-1\le 2kst{e}^{2kst}$, and thereforeThus $A_{X(l)}(t)=A_R(t)$ is the desired Lieb-Robinson localization by Haar twirling. The displayed choice assumes the usual positive constants $\mu ,s$; if the interaction bound has $s=0$, the error is identically zero and any positive $v$ works. For negative times the same argument uses $|t|$.

Use an $n$-qubit state as the witness. Uniformly choose a term $h_j$ and perform the two-outcome positive operator-valued measure $\{I-h_j,h_j\}$. Regard the second outcome as an energy flag. For witness $\rho$, its probability isThis is an efficient local-energy measurement verification: a unitary operator on the $k$ relevant qubits and one ancilla qubit can implement the measurement using the isometryThe local dimension is constant; the term matrices and the measurement can be approximated to inverse-polynomial accuracy. Uniform choice of $j$ and the classical threshold comparison also have polynomial overhead.

SetIn the nontrivial case $0\le a<b\le M$, $\delta$ is inverse-polynomial. Cases outside this energy range are decided directly from $0\le H\le MI$.

Ask for $r$ witness registers, measure an independently sampled local term on each, and accept if the mean number of energy flags is at most $\theta$. An honest product ground state on a YES instance has flag probability at most $p_{\mathrm{y}}$. The Hoeffding inequality gives failure at most $e^{-r{\delta }^2/2}$. For completeness, the required Bernoulli concentration estimate follows by bounding the centered log moment-generating function $f(\lambda )=\ln \mathbb{E}{e}^{\lambda (Y-p)}$: $f(0)=f^{\prime }(0)=0$ and $f^{\prime \prime }(\lambda )$ is a tilted Bernoulli variance, at most $1/4$. Hence $f(\lambda )\le {\lambda }^2/8$; exponential Markov's inequality optimized at $\lambda =4u$ yields $\mathrm{Pr}(\overline{Y}-p\ge u)\le e^{-2r{u}^2}$, and the lower-tail version follows in the same way.

Soundness must also cover entangled witness registers. The averaged flag effect is $R=H/M$. All its eigenvalues are at least $p_{\mathrm{n}}$ on a NO instance. The repeated verifier's acceptance effect is a polynomial in the commuting tensor-factor operators $R_j$ and $I-R_j$. In their product eigenbasis its eigenvalues are lower-tail probabilities of independent Bernoulli variables whose parameters are each at least $p_{\mathrm{n}}$. Their lower tail is bounded by the identical-parameter tail at $p_{\mathrm{n}}$, using the same uniform-variable coupling as in QMA parallel repetition with entangled witnesses. Therefore its operator norm is at most $e^{-r{\delta }^2/2}$, which bounds acceptance for any entangled witness.

Choosinggives completeness at least $2/3$ and soundness at most $1/3$. A slightly larger constant absorbs implementation errors; for example, approximate each flag effect in norm to at most $\delta /8$, retain a gap at least $3\delta /4$, and increase $r$ by a constant factor. The total witness length $rn$ and circuit size are polynomial. Thus

The absolute value of an operator, trace norm, and operator norm are respectivelyThe positive square root of an operator makes $|L|$ well defined. Its eigenvalues $s_j\ge 0$ are the singular values of $L$, so $\Vert L{\Vert }_1={\sum }_j{s}_j$ and $\Vert L{\Vert }_{\mathrm{op}}={\max }_j{s}_j$.

Use the polar decomposition of a bounded operator $L=U|L|$, with the given unitary operator $U$. In an orthonormal eigenbasis $e_j$ for $|L|$,The operator norm assumption and unitarity imply $|⟨e_j|ZU|e_j⟩|\le \Vert Z{\Vert }_{\mathrm{op}}\le 1$. The triangle inequality therefore yields the required trace estimate:Choosing $Z=U^{†}$ attains equality, which also gives the finite-dimensional trace duality formula $\Vert L{\Vert }_1={\max }_{\Vert Z{\Vert }_{\mathrm{op}}\le 1}|\mathrm{Tr}(ZL)|$. A singular $L$ causes no difficulty: in a finite-dimensional square space the polar factor can be extended to a unitary operator on the complementary subspaces.

Take $|\psi ⟩$ to be normalized and let ${\Pi }_{\mathcal{G}}$ be the orthogonal projection onto the good subspace. The two reflection operators areThe first fixes the hyperplane orthogonal to $|\psi ⟩$ and negates its normal direction. The second fixes ${\mathcal{G}}^{\perp }$ and negates $\mathcal{G}$. Both are Hermitian operators and unitary operators, since their defining projectors square to themselves. For a non-normalized nonzero vector, divide $|\psi ⟩⟨\psi |$ by $⟨\psi |\psi ⟩$.

Put $p=⟨\psi |{\Pi }_{\mathcal{G}}|\psi ⟩={\sin }^2\theta$, with $0<\theta <\pi /2$, and define normalized orthogonal good and bad vectorsThen $|\psi ⟩=\sin \theta |g⟩+\cos \theta |b⟩$. The amplitude amplification theorem states that the iterate $Q=-I_{\psi }{I}_{\mathcal{G}}$ satisfies, for every nonnegative integer $m$,To prove it, the span of $|g⟩,|b⟩$ is invariant under both reflection operators. In that ordered orthonormal basis,Multiplying the second matrix into $(\sin \theta ,\cos \theta {)}^T$ advances the angle by $2\theta$. Induction proves the formula, and projecting onto $\mathcal{G}$ gives $p_m$. The minus sign in $Q$ is just a global phase choice, but it makes the displayed state formula exact.

If $p$ is known, choose $m$ nearest to $\pi /(4\theta )-1/2$, with $m\ge 0$. The final angle is within $\theta$ of $\pi /2$, so $p_m\ge {\cos }^2\theta =1-p$. For $p\le 1/2$, this is at least $1/2$ after $O(1/\sqrt{p})$ iterations, and is close to one when $p$ is small. For $p>1/2$, measuring the initial state already has success above $1/2$. Repeated independent attempts with a success test give any fixed desired success probability. This is the quadratic improvement of amplitude amplification over repeated sampling, which takes $O(1/p)$ attempts. At $p=0$ there is no good component to amplify, and at $p=1$ the initial state is already entirely good; these endpoint cases do not require the undefined normalized component vectors.

For first-order two-local Hamiltonian simulation, regard each term $H_k$ as a Hermitian operator, as is standard for a local-Hamiltonian decomposition. The Lie-Trotter product formula for bounded operators statesin operator norm; the reversed factor ordering is just as valid as the displayed one. Take $A_k=-i{H}_k$ and use the circuitEach factor is a unitary operator on at most two qubits and hence is one allowed two-qubit gate, with a spectator identity for a one-qubit term. The factors in the rightmost part of each written product act first. Merely quoting the limit would not establish the requested rate, so we give a quantitative error bound.

For two Hermitian operators $A,B$ and $t\ge 0$, differentiatinggives $F^{\prime }(s)=i{e}^{-i(t-s)(A+B)}[B,e^{-isA}]e^{-isB}$. The commutator identity obtained by differentiating $e^{iuA}B{e}^{-iuA}$ implies $\Vert [B,e^{-isA}]\Vert \le s\Vert [A,B]\Vert$. Since the surrounding factors are unitary, integration yields the first-order unitary product-formula error boundSplitting the terms successively and using the triangle inequality therefore gives the one-step estimatefor $M>1$. Here submultiplicativity of the operator norm gives $\Vert [H_k,H_{\ell }]\Vert \le 2\Vert H_k\Vert \Vert H_{\ell }\Vert <2$. The ordering chosen for the products changes the sign of local commutators, not this norm bound. Applying part (i) to the $r$ steps givesChoose $r=\max (1,\lceil M(M-1)/(2ϵ)\rceil )$. For $M>1$ the error is strictly below $ϵ$; $M=1$ is exact already with $r=1$. The quantum circuit contains $Mr$ two-qubit gates. A sufficient gate count isThus the worst-case polynomial degree in this first-order bound is six. Terms with disjoint supports commute and may improve the count in particular decompositions, but the stated assumptions alone suffice for this bound, without a further sparsity or bounded-degree condition. The requested gate model allows arbitrary two-qubit unitaries, so no additional compilation error is incurred here.

The spectral norm is the induced Euclidean operator norm:It is the largest singular value. From the definition, multiplying on either side by a unitary operator leaves the norm unchanged: right multiplication permutes the unit sphere of possible inputs, and left multiplication preserves output lengths. In particular every unitary operator has norm one.

For unitary product telescoping, replace the factors one at a time. With empty products interpreted as the identity,All cross terms cancel. By the triangle inequality and unitary invariance of the spectral norm,This bound is independent of the dimension and does not assume that any of the factors commute.

For $H=H_0+H_{\mathrm{int}}$, the interaction picture has $H_I(t)=e^{i{H}_0(t-t_0)}{H}_{\mathrm{int}}(t)e^{-i{H}_0(t-t_0)}$ andThe trial solution $|\psi (t){⟩}_I=U(t,t_0)|\psi (t_0){⟩}_I$ works for arbitrary initial states precisely whenFor $t\ge t_0$, the Dyson series writesIt is also ${\sum }_{n\ge 0}(-i{)}^n/n!{\int }_{[t_0,t{]}^n}{d}^{n}sT\{H_I(s_1)\cdots H_I(s_n)\}$: the cube consists of $n!$ ordered simplexes with the same time-ordered integrand. Differentiating the nested integrals makes the latest operator $H_I(t)$ the leftmost factor, leaving the $(n-1)$-fold integral. Thus ${\partial }_{t}U=-i{H}_I(t)U$ and at $t=t_0$ every nonconstant term vanishes. This proves the differentiation of the Dyson time-ordered exponential identity without incorrectly treating the generally noncommuting Hamiltonians as scalars.

For a Hermitian Hamiltonian operator, ${\partial }_t(U^{†}U)=i{U}^{†}{H}_{I}U-i{U}^{†}{H}_{I}U=0$, so the initial condition makes $U$ a unitary operator. Evolution backwards in time is its inverse, with anti-time ordering. For bounded norm-continuous $H_I$ the series converges in operator norm and the differentiation is justified by the factorial estimates. In quantum field theory, unbounded operators and products of fields require domain or regulator assumptions; this computation is the formal perturbative identity unless those analytic assumptions are supplied.

Let $V$ denote the coherent exact quantum phase estimation unitary operator from the preceding part, with the data quantum register written first. For every eigenvector of $A$,because $e^{2\pi iA}|u_j⟩=e^{2\pi i{\lambda }_j}|u_j⟩$ and ${\lambda }_j=c_j/2^n$ lies in $[0,1)$, with no phase-aliasing ambiguity. Apply $V$ coherently to the input quantum state and a clean quantum register; do not measure the eigenvalue label. By linearity it produces ${\sum }_j{\beta }_j|u_j⟩|c_j⟩$.

Append a quantum ancilla in $|0⟩$, and use the provided quantum variable rotation with ${\theta }_{c_j}=\arcsin {\lambda }_j$, choosing the nonnegative cosine. The resulting normalized quantum state isThe sum covers both terms; this fixes the potentially ambiguous sum placement in the abbreviated display. The orthonormal basis of eigenvectors and the unit length of each rotated quantum ancilla show that its norm is one. This is quantum spectral filtering by multiplication, not the reciprocal rotation used in the HHL algorithm.

Use the Minkowski metric $g^{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$ throughout. Taking the Hermitian conjugation of the massive Dirac equation and using the Dirac adjoint givesTaking the matrix transpose and multiplying by the charge conjugation matrix $C$ yieldsThus the charge conjugation of a Dirac field preserves the free massive Dirac equation:The unitarity of $C$ and ${\gamma }^0{C}^{-1}{\gamma }^0=-C^{-1}{\gamma }^{0T}$ also giveConsequentlyHere $\hat{C}$ is the unitary operator acting on the relativistic quantum field, whereas $C$ acts on Dirac spinor indices.

The bosonic annihilation operator $a^{\lambda }(p)$ removes a vector particle of four-momentum $p$ and polarization vector ${\epsilon }^{\mu ,\lambda }(p)$. The bosonic creation operator $c^{\lambda †}(p)$ creates its antiparticle with the same four-momentum and polarization label. The polarization vector describes the corresponding spin-one mode expansion of a free field.

Because charge conjugation acts through a unitary operator, it leaves the numerical polarization vector coefficients unchanged. Applying the specified charge conjugation to the operators givesComparison with the Hermitian conjugation of the original expansion therefore givesThe factor ${\eta }_C$ is the intrinsic charge-conjugation phase.

The flavour-diagonal weak neutral current coupled to the real Z boson has real vector and axial current couplings. With ${\eta }_C(Z)=-1$, part (e) shows that this tree-level interaction is invariant under CP symmetry, although the simultaneous presence of the Dirac current and axial current generally violates charge conjugation and parity symmetry in quantum field theory separately. This conclusion concerns the tree-level weak neutral current, rather than every interaction or radiative effect in the Standard Model.

The weak charged current couples distinct quark flavours and involves the Cabibbo-Kobayashi-Maskawa matrix:Here $u_{Li}$ and $d_{Lj}$ are left-handed quark fields. A CP symmetry transformation exchanges the two conjugate-field operators; because it acts through a unitary operator, it does not complex-conjugate their numerical coefficients. Equality with the original interaction therefore requires a Cabibbo-Kobayashi-Maskawa matrix that can be made real by quantum field rephasings. With three generations, a physical CP-violating phase remains in general, measured by the Jarlskog invariant. Therefore the W boson interactions can violate CP symmetry. Merely seeing a complex coefficient is insufficient: the phase must survive quantum field rephasing, unlike the real two-generation Cabibbo angle rotation.

In time reversal acoustics, the array records the incoming signal, reverses each recorded time trace, and re-emits it through the same medium. This is phase conjugation in the frequency domain. For the time-harmonic wave convention $e^{-i\omega t}$, reversing a real time trace replaces its positive-frequency wave amplitude by its complex conjugate. The medium must remain unchanged between recording and re-emission.

Let $f(\mathbf{z})={\Psi }_0(0,\mathbf{z})$ and $G_L(\mathbf{\zeta },\mathbf{z})=G(L,\mathbf{\zeta };0,\mathbf{z})$, with the chosen source and array normalizations incorporated into the Green function. Let $a(\mathbf{\zeta })\ge 0$ be the array's aperture weight, equal to the indicator of its receiving region for an ideal uniform array. The recorded field isBy wave reciprocity, back-propagation has the same Green function with the source and receiver exchanged. Therefore the physically re-emitted, back-propagated wave amplitude isIf $Hf=R$ and $P_A$ multiplies by $a$, then ${\Psi }_B=H^T{P}_A\overline{Hf}$, where $T$ denotes the transpose without conjugation. Taking a final complex conjugate instead defines the adjoint reconstruction $\overline{{\Psi }_B}=H^{\ast }{P}_{A}Hf$. This distinction prevents an erroneous conjugation in the time reversal operator.

For a localized Gaussian beam or acoustic point source in a homogeneous medium, a finite aperture admits a limited range of angles. The focal width is of order $\lambda L/A$ when $A$ denotes the aperture diameter. In a random medium, multiple scattering creates paths with a larger angular spread. Each reversed path retraces its route, and the paths interfere constructively at the source. This can produce a larger effective aperture in time reversal and a narrower focus, even though the unreversed field has a complicated speckle pattern.

This comparison concerns a homogeneous reference medium; a deterministic heterogeneous medium can also provide useful multipath propagation. Suitable scale limits or frequency and spatial averaging can make refocusing self-averaging. Such self-averaging is not automatic for every monochromatic source and every random realization. Wave absorption, changing medium parameters, unresolved paths or poor array coverage can spoil refocusing. With complete capture of the propagating modes and a lossless unitary operator $H$, ideal adjoint reconstruction is already exact in either medium. Random scattering can improve finite-aperture wave focusing through angular diversity. See the regime-dependent analysis in Statistical stability in time reversal.