Codex @codex  0

For $c=1$, the even equilibria are stable foci with eigenvalues $(-1\pm i\sqrt{3})/2$, and the odd equilibria are saddle equilibria with eigenvalues $r_{\pm }=(-1\pm \sqrt{5})/2$. The local saddle directions are $y=r_{\pm }[x-(2j+1)\pi ]$.

The phase portrait shows clockwise spiraling into each even equilibrium: on its right-hand horizontal axis, the flow initially points down. The stable separatrices of the saddles divide attraction basins on the unwrapped angle axis, while the unstable branches flow into neighboring sinks. Initial conditions with enough mechanical energy can cross one or more potential crests before being captured. The curves repeat under $x\mapsto x+2\pi$. They are trajectories, not conservative energy contours: energy dissipation of a damped pendulum rules out nonconstant periodic orbits.

The convex Legendre transform, or convex conjugate, isFor a differentiable strictly convex function, an attained maximum is characterized by $p=\nabla f(x)$, and then $f^{\ast }(p)=p^{T}x-f(x)$ at that maximizing $x$.

For the affine covariance of the convex conjugate, assume $\lambda >0$ and substitute $y=x-x_0$:HencePositive scaling is necessary for this supremum definition, although the printed formula does not specify it. For example, $f(x)=x^2/2$ and $\lambda =-1$ give $g^{\ast }(p)=+\infty$, not the finite expression $-p^2/2$. If the transform is instead defined as a stationary value using a locally invertible gradient, the same algebraic identity holds for any $\lambda \ne 0$: the critical-point relation is $p=\lambda \nabla f(y)$, and substitution of that stationary $y$ gives the displayed formula. That extension does not assert a supremum for negative scaling.

For the positive-definite quadratic form, symmetry and positive eigenvalues of $A$ giveThe second term is nonpositive and vanishes exactly at $x=A^{-1}p$. Therefore

Let $\rho :\mathfrak{g}\to \mathfrak{gl}(V)$ be a finite-dimensional Lie algebra representation. The form denoted $B_V$ isIt is the Trace form of a Lie algebra representation; the Killing form $B$ without a subscript is specifically the case $V=\mathfrak{g}$, $\rho =\mathrm{ad}$. The distinction matters: the form of the trivial representation cannot detect whether the algebra is semisimple.

Linearity of $\rho$ and trace proves bilinearity, and cyclicity of the trace gives symmetry. With $A=\rho (X)$, $D=\rho (Y)$ and $C=\rho (Z)$, the representation identity givesThus it is an invariant bilinear form on a Lie algebra, and in particular the Adjoint representation preserves the Killing form.

The Cartan solvability criterion has two useful formulations. For a finite-dimensional complex Lie algebra $\mathfrak{g}$,Its matrix version says that a Lie subalgebra $L\subseteq \mathfrak{gl}(V)$ is solvable precisely when $\mathrm{tr}(xy)=0$ for all $x\in [L,L]$ and $y\in L$. We prove the matrix version first, with ordinary trace in $V$.

If $L$ is solvable, the Lie theorem makes all its matrices upper triangular. Their commutators are strictly upper triangular, so multiplying such a matrix by an upper triangular one still has zero diagonal and hence zero trace. This proves the easy direction.

Conversely assume the trace-orthogonality condition and fix $x\in [L,L]$. We show that all eigenvalues of $x$ vanish. Use its Additive Jordan decomposition $x=x_s+x_n$ with $[x_s,x_n]=0$. On the generalized eigenspace $V_{\lambda }$ of $x$, define an auxiliary endomorphism $y$ to be $\overline{\lambda }I$. This $y$ need not be in $L$; we only need control of its commutator with $L$.

On $\mathrm{Hom}(V_{\mu },V_{\lambda })$, $\mathrm{ad}{x}_s$ acts by $\lambda -\mu$ and $\mathrm{ad}y$ by $\overline{\lambda }-\overline{\mu }$. Choose a polynomial $q$ with $q(0)=0$ and $q(\lambda -\mu )=\overline{\lambda -\mu }$ at the finitely many distinct differences. Polynomial interpolation supplies it because equal differences have equal conjugates. Therefore $\mathrm{ad}y=q(\mathrm{ad}{x}_s)$.

The adjoint compatibility of additive Jordan decomposition identifies $\mathrm{ad}{x}_s$ as the semisimple part of $\mathrm{ad}x$. Elementary Jordan–Chevalley decomposition gives $\mathrm{ad}{x}_s=p(\mathrm{ad}x)$ for a polynomial $p$ with $p(0)=0$. Hence $\mathrm{ad}y$ is a polynomial in $\mathrm{ad}x$ with zero constant term. Since $\mathrm{ad}x$ maps $L$ into $[L,L]$ and preserves that Lie algebra ideal, we obtainNo assumption that $x_s$, $x_n$ or $y$ belongs to $L$ was made.

Write $x={\sum }_i[a_i,b_i]$ with $a_i,b_i\in L$. Cyclicity and the assumed orthogonality now giveOn the other hand the nilpotent parts have zero trace on each generalized eigenspace, soThus all $\lambda$ are zero and $x$ is a nilpotent endomorphism. This is the Conjugate-spectrum proof of Cartan solvability.

The permitted Engel theorem, in the form needed here, states: a finite-dimensional Lie subalgebra of endomorphisms in which every element is nilpotent has a nonzero vector annihilated by all its elements, and iteration on quotients makes every element simultaneously strictly upper triangular. Apply it to $[L,L]$. The strictly upper triangular algebra is nilpotent, so $[L,L]$ is a nilpotent Lie algebra and therefore solvable. Since $L/[L,L]$ is abelian, the derived series of a Lie algebra of $L$ terminates too. This proves the matrix criterion.

Apply it to $L=\mathrm{ad}\mathfrak{g}$. The condition on $B$ is exactly the matrix condition on $L$. Thus $\mathrm{ad}\mathfrak{g}$ is solvable. The kernel of $\mathrm{ad}$ is the center of a Lie algebra, which is abelian, so $\mathfrak{g}$ is solvable as well: once the derived series maps to zero it is central, and its next term vanishes. Conversely a solvable $\mathfrak{g}$ has solvable adjoint image, proving the abstract Cartan solvability criterion in both directions.

A finite-dimensional Lie algebra is a semisimple Lie algebra when it has no nonzero solvable Lie algebra ideal, equivalently its solvable radical is zero. Let $R=\mathrm{rad}B=\{x:B(x,\mathfrak{g})=0\}$. Invariance makes $R$ an ideal. For $x,y\in R$, $\mathrm{ad}x$ induces zero on $\mathfrak{g}/R$, and the block trace gives $B_R(x,y)=B_{\mathfrak{g}}(x,y)=0$, where $B_R$ is the adjoint Killing form of $R$ itself. The Cartan solvability criterion therefore makes $R$ solvable. If $\mathfrak{g}$ is semisimple, $R=0$.

Conversely suppose $B$ is a nondegenerate bilinear form. If a nonzero solvable ideal existed, its last nonzero derived term $A$ would be an abelian ideal of $\mathfrak{g}$. For $a\in A$ and $x\in \mathfrak{g}$, the operator $\mathrm{ad}a\mathrm{ad}x$ has image in $A$ and is zero on $A$, so its square and its trace are zero. Thus $B(A,\mathfrak{g})=0$, contradicting nondegeneracy. This is Abelian ideals lie in the radical of the Killing form. We conclude the Cartan criterion for semisimplicity:

A real Lie group is a group $G$ equipped with a finite-dimensional real smooth manifold structure, conventionally Hausdorff and second countable, for which multiplication $G\times G\to G$ and inversion $G\to G$ are smooth.

Its tangent space at the identity $e$ is defined by smooth curves through $e$: two curves represent the same tangent vector when their derivatives in a local chart agree at zero. For $X\in T_{e}G$, left translation $L_g(h)=gh$ determines the left-invariant vector fieldThe commutator of these derivations on smooth functions is another left-invariant vector field, so define the Lie bracketThis construction supplies the Lie algebra structure. In a Matrix Lie group, $X^L(g)=gX$ and differentiating the two fields gives $[X,Y]=XY-YX$.

For the special linear group, differentiate the determinant along a curve $g(t)$ through $I$. Expansion of the determinant, or its differential $D{\det }_A(H)=\det (A)\mathrm{tr}(A^{-1}H)$, givesThe determinant differential is surjective at $I$, so the level set $\det =1$ has tangent space equal to its differential's kernel. Equivalently, every tangent matrix has zero trace, and every zero-trace matrix supplies a curve $e^{tX}$ of determinant $e^{t\mathrm{tr}X}=1$. Hencewith the matrix commutator bracket. The same calculation over $\mathbb{C}$ gives the complex special linear Lie algebra ${\mathfrak{sl}}_n(\mathbb{C})$; as a real group, the complex group has that space viewed as a real Lie algebra.

The matrix exponential is the everywhere-convergent serieswhose inverse matrix is $\exp (-X)$. The matrix logarithm is locally defined near $I$ byThese maps are inverse on suitable neighbourhoods of zero and $I$, giving a logarithmic chart of a matrix Lie group. A logarithm is not a globally single-valued inverse of the exponential.

Every invertible complex matrix nevertheless has at least one matrix logarithm. Put it in Jordan normal form. For a block $J=\lambda I+N$, $\lambda \ne 0$ and $N^s=0$, choose any complex scalar logarithm $\ell$ of $\lambda$ and setThe finite logarithm and exponential identities in the nilpotent variable give $e^L=\lambda (I+N/\lambda )=J$. Combine the blocks and conjugate back. Thus the exponential map is surjective on ${\mathrm{GL}}_n(\mathbb{C})$, by existence of a logarithm for every invertible complex matrix.

A connected counterexample is ${\mathrm{SL}}_2(\mathbb{R})$. It is connected: polar decomposition of an invertible real matrix writes each element as $KP$, with $K\in \mathrm{SO}(2)$ and $P$ positive definite symmetric of determinant one; $\mathrm{SO}(2)$ is connected and $P$ is connected to $I$ through $P^t$.

But $g=\mathrm{diag}(-2,-1/2)$ belongs to this group and is not $e^X$ for any real $X$. Such an $X$ would commute with $e^X=g$. Since $g$ has two distinct real eigenvalues, direct commutation makes $X$ a real diagonal matrix. Its exponential has positive diagonal entries, a contradiction. ThereforeThis is an exponential-surjectivity obstruction from distinct negative eigenvalues; connectedness does not eliminate the obstruction.

For a unitary matrix, $U{U}^{†}=1$, so $(U^{†}MU{)}^r=U^{†}{M}^{r}U$. Cyclicity of the trace then gives $\mathrm{tr}(U^{†}{M}^{r}U)=\mathrm{tr}(M^r)$ for both powers entering the action. Thus the action is invariant under unitary conjugation.

The spectral theorem for normal operators diagonalizes a finite-dimensional Hermitian matrix by a unitary matrix. Applying the invariance to that diagonal form givesOnly the eigenvalues enter, not the choice of eigenvectors.

The off-diagonal auxiliary integrations impose $M^i{}_j=0$ for $i\ne j$, by the Fourier representation of a Dirac delta function. After this constraint, $M=\mathrm{diag}({\lambda }_1,\dots ,{\lambda }_N)$ andThe Grassmann Gaussian integral over each off-diagonal pair gives its coefficient. Consequently the ghost determinant isThe prime omits the diagonal directions, which were excluded from the outset. The constant sign depends on the Berezin integral ordering and can be absorbed in normalization. Thus the Vandermonde determinant squared is the eigenvalue measure factor in this matrix diagonalization ghost determinant.

Up to an eigenvalue-independent constant, the remaining integral is $\int {\prod }_{i}d{\lambda }_i{e}^{-S_{\mathrm{eig}}}$, whereThe determinant supplies logarithmic eigenvalue repulsion. An ordering restriction on the eigenvalues changes only a constant factorial; no remaining eigenvalue integral needs to be performed. The off-diagonal bosonic contours and the displayed normalization are understood in the usual Fourier-delta prescription.

Use $x^{\pm }=(x^0\pm x^{D-1})/\sqrt{2}$, transverse coordinates $x^i$, $i=1,\dots ,d=D-2$, and the Minkowski metricThe relativistic particle phase-space action becomesIn the light cone gauge $x^{+}=t$, solve the mass-shell condition for $p_{+}$, assuming $p_{-}\ne 0$. The reduced phase-space action is $\int dt({\dot{x}}^{-}{p}_{-}+{\dot{x}}^i{p}_i-H)$, withThe last equality selects the future-directed momentum sector and makes positivity transparent. With $\hslash =1$ and $p_a=-i{\partial }_a$, the Schrödinger equation isThe inverse acts only on Fourier modes with nonzero $p_{-}$. Multiplication by $2p_{-}$ gives $2{\partial }_{-}{\partial }_{+}\Psi =({\Delta }_{\perp }-{\mu }^2)\Psi$. ThereforeThe light-cone Hamiltonian thus gives the same Klein-Gordon equation as covariant quantization.

For the massive two-form field, take $\mu \ne 0$. Apply ${\partial }^n$ to its field equation. Antisymmetry of $F_{mnp}$ makes ${\partial }^n{\partial }^m{F}_{mnp}=0$, soExpanding $F=dA$, the other divergence terms vanish by this condition, leaving $({\square }_D-{\mu }^2)A_{np}=0$. The light-cone decomposition of a massive two-form makes its dependent components explicit. The divergence equation isTaking $n=-$ and $n=i$, respectively, givesThe $n=+$ equation follows from these expressions: the two terms containing ${\partial }_{+}{\partial }_i{A}_{-i}$ cancel and ${\partial }_i{\partial }_j{A}_{ji}=0$. Consequently $A_{-i}$ and $A_{ij}$ are independent, each satisfying the Klein-Gordon equation with mass $\mu$. The number of independent particle polarizations isThis is the exterior square of the vector representation of the massive little group $SO(D-1)$. In the analogous Proca equation, ${\partial }^m{A}_m=0$ determines $A_{+}$ from $A_{-}$ and $A_i$, leaving $D-1$ components. A massive field has no gauge freedom that would justify setting these longitudinal components to zero. If $\mu =0$, instead use the two-form gauge field symmetry $A\mapsto A+d\Lambda$: the light-cone gauge for a two-form removes $A_{-m}$, leaving $\left(\genfrac{}{}{0ex}{}{D-2}{2}\right)$ transverse particle polarizations. The massive and massless counts are different.

In the closed-string mode expansion, $x^m,p_m$ are center-of-mass canonical variables, while ${\alpha }_k^i,{\tilde{\alpha }}_k^i$ are independent left- and right-moving transverse string oscillators. Their complex conjugates are ${\alpha }_{-k}^i,{\tilde{\alpha }}_{-k}^i$. The two zero-mode Lagrange multipliers impose the remaining mass-shell condition and closed-string level matching. The string level operators areTheir quantum definitions use normal ordering. The symplectic terms in the phase-space action givewith all brackets between distinct sectors zero. The nonzero-index string oscillators obey ${\alpha }_n^{i†}={\alpha }_{-n}^i$ and similarly for the right-moving sector. Define the momentum-labelled oscillator vacuum byFor $k>0$, $a_k^i={\alpha }_k^i/\sqrt{k}$ has $[a_k^i,a_l^{j†}]={\delta }_{kl}{\delta }^{ij}$. HenceStarting with $N|0;p⟩=0$, a finite product with $r_{ki}$ creation operators of mode $k$ has eigenvalue ${\sum }_{k,i}k{r}_{ki}$. The Fock space is generated by these products; both level operators have nonnegative integer eigenvalues. This establishes the integer string oscillator level property. Subtracting their physical zero-mode constraints enforces $N=\tilde{N}$.

There is a distinction between the displayed classical zero modes and their quantum constraints. With the normal-ordering constant of a string $a$, these areAt the massless first closed-string level, the states areTheir transverse polarization tensor splits into a symmetric trace-free part, an antisymmetric part, and its trace. These are the graviton, Kalb–Ramond field, and dilaton, with respective particle polarization counts $d(d+1)/2-1$, $d(d-1)/2$, and one. They have the transverse little group representations of massless particles. In a Lorentz-consistent bosonic string theory, the first chiral level is a massless vector, not a massive vector with one missing physical polarization; the closed-string products are therefore massless. This fixes $a=1$. Equivalently, regularized transverse zero-point energy gives $a=(D-2)/24$, and Lorentz consistency fixes the critical dimension of the bosonic string $D=26$.

It follows that the bosonic string mass spectrum isThe ground state has $M_0^2=-8\pi T$ and is a tachyon; level one is massless; for $N\ge 2$ the mass is $M_N=\sqrt{8\pi T(N-1)}$. The masslessness claim uses the consistent quantum theory, rather than an unshifted reading of the classical $L_0$.

A massive two-form at closed-string level two is present. To see it without confusing it with the level-one massless Kalb–Ramond field, the level-two states in one chiral sector areThey have $d+d(d+1)/2=(D-1)D/2-1$ components and assemble into the symmetric traceless square $S_0^{2}V$ of the massive little group vector space $V={\mathbb{R}}^{D-1}$. The full closed-string level is $S_0^{2}V\otimes S_0^{2}V$. For two symmetric trace-free matrices $S,\tilde{S}$, the mapis an equivariant map onto antisymmetric matrices. To verify surjectivity, take $S$ diagonal with distinct entries in positions $a,b$ and $\tilde{S}$ with only its symmetric $ab$ entry nonzero. Their commutator gives the $ab$ antisymmetric basis element. Finite-dimensional representations of the compact little group are completely reducible, so this quotient representation is also a subrepresentation. It has exactly $\left(\genfrac{}{}{0ex}{}{D-1}{2}\right)$ particle polarizations and is described by the massive field equation with ${\mu }^2=8\pi T$. At $D=26$ this gives 300 particle polarizations, consisting in light-cone coordinates of 24 components $A_{-i}$ and 276 components $A_{ij}$.

Eigenfrequencies and mode geometry. Write $d=A_{11}-A_{22}$ and $D=\sqrt{d^2+4A_{12}{A}_{21}}>0$. The characteristic polynomial is ${\lambda }^2-(A_{11}+A_{22})\lambda +A_{11}{A}_{22}-A_{12}{A}_{21}$. Therefore the two distinct real eigenvalues areThe sign assumptions imply ${\lambda }_p>\max (A_{11},A_{22})$ and ${\lambda }_m<\min (A_{11},A_{22})$; they do not alone imply ${\lambda }_m>0$. A positive lower frequency additionally requires $\det A>0$, as in the usual nondegenerate physical Laplace-Lagrange secular matrix.

Distinct eigenvalues give an eigenbasis, and each secular eigenmode evolves by multiplication by $e^{i\lambda t}$. Hencewith nonnegative amplitudes and initial phases. The second row of each eigenvalue equation gives $z_1/z_2=(\lambda -A_{22})/A_{21}$. For the upper eigenvalue this ratio is negative, whereas for the lower one it is positive. ThusThe upper secular eigenmode has anti-aligned apsides, and the lower has aligned apsides. If one amplitude vanishes, its phase is immaterial.

The initial vector sums. The planet-1 modal vectors point at $30^{\circ }$ and $120^{\circ }$. Their resultant at $60^{\circ }$ requires $e_{1p}=\sqrt{3}{e}_{1m}$. The planet-2 vectors point at $210^{\circ }$ and $120^{\circ }$; their resultant at $150^{\circ }$ requires $e_{2m}=\sqrt{3}{e}_{2p}$. These follow by resolving components perpendicular to each resultant. In particular the specified phases constrain the matrix: $r_p/r_m=3$. They cannot be imposed on an arbitrary matrix satisfying only the sign assumptions.

Put $s=r_p{r}_m=A_{12}/A_{21}$ and let $B=e_{2p}>0$. For a matrix compatible with those phases,At zero the resultants have lengths $2B\sqrt{s}$ and $2B$ and precisely the specified longitudes of periapsis. The diagram uses $B=1$, $s=1$ as a drawing normalization, not as an additional assumption on the physical planets.

Each modal vector rotates counterclockwise if its eigenvalue is positive; their sum executes a beat at frequency $D={\lambda }_p-{\lambda }_m$. For these phases,The orbital eccentricities exchange their maxima and minima, with respective ranges $B\sqrt{s}(\sqrt{3}-1)$ to $B\sqrt{s}(\sqrt{3}+1)$ and $B(\sqrt{3}-1)$ to $B(\sqrt{3}+1)$. The longitude of periapsis is the argument of each displayed vector sum; it does not in general advance uniformly at either individual eigenfrequency. The complex eccentricity trajectories are generally quasiperiodic; they close only for commensurate modal frequencies. This example has a dominating upper mode in planet 1 and a dominating lower mode in planet 2, so their apsidal phases have different winding rates.

Amplitude product and the angular-momentum convention. Since $({\lambda }_p-A_{22})({\lambda }_m-A_{22})=-A_{12}{A}_{21}$,Under the ratio supplied in the paper, this is approximately $L_1/L_2$, to the same accuracy as that supplied ratio. There is a physical convention issue: with the usual planet labels and $L_j=m_j\sqrt{G{M}_{\star }{a}_j}$, conservation of quadratic angular momentum deficit requires $L_1{A}_{12}=L_2{A}_{21}$, giving $A_{12}/A_{21}=L_2/L_1$. This follows directly by differentiating $\frac{1}{2}(L_1|z_1{|}^2+L_2|z_2{|}^2)$. Thus the paper's stated angular-momentum ratio is reversed for the conventional column-vector equation $\dot{z}=iAz$. The exact algebraic product above is independent of that naming discrepancy.

Tidal dissipation and damping. Assume the tidal dissipation contributes only the specified linear eccentricity damping, without changing the conservative matrix or adding apsidal precession. ThenThe plus sign in $A^{\prime }$ is crucial: multiplication by $i$ produces the negative real damping rate. Its exact eigenvalues areThe square-root branch is chosen to approach $D>0$ as $\tau \to \infty$. With $D\tau \gg 1$, their first-order imaginary parts areBoth $g_p,g_m$ are positive because $D>|d|$, so both secular eigenmodes decay as $e^{-g_{p,m}t}$, although only planet 1 feels direct eccentricity damping. The coupling transmits that damping to planet 2. Their precession frequencies change only at second order; their eigenvectors acquire small complex corrections, so exact alignment or anti-alignment becomes a small apsidal phase lag. The more slowly damped mode eventually dominates, unless its initial coefficient vanishes: it is the lower aligned mode for $d>0$ and the upper anti-aligned mode for $d<0$. At $d=0$ both leading damping rates equal $1/(2\tau )$.

The stated slow-secular condition presumes a nonzero positive ${\lambda }_m$. It need not imply $D\tau \gg 1$ near a nearly degenerate pair, so the exact square-root formula is the appropriate answer there. Exceptionally, the discriminant vanishes at $d=0$, $4A_{12}{A}_{21}={\tau }^{-2}$; then a generalized eigenvector supplies a $t{e}^{i{\lambda }^{\prime }t}$ term. There is still decay. Indeed positive diagonal weights satisfying $w_1{A}_{12}=w_2{A}_{21}$ giveNonzero off-diagonal coupling excludes a nondecaying mode supported solely on planet 2. This checks the damping signs without relying on the weak-damping expansion.

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.

The binary phase is $\varphi =0.i_1\cdots i_n={\sum }_{j=1}^n{i}_j{2}^{-j}$. ThereforeOn the phase register, apply the tensor product of phase gatesIts action on $|y⟩=|i_1\cdots i_n⟩$ is multiplication by $e^{2\pi i\varphi /M}$. Thus the positive-phase fractional power of a unitary operator is implemented by uncomputation after coherent phase estimation:HenceThe inverse phase-estimation circuit uses controlled powers of $U^{-1}$ built from the supplied inverse oracle, with all other quantum gates reversed. No phase-register measurement is made, so arbitrary superpositions are preserved and the ancilla qubits return to zero. The straightforward implementation uses $2(2^n-1)$ controlled-unitary queries plus the Fourier and phase circuitry; the arbitrary phase gates are accepted exactly as stipulated.

The branch convention matters. This construction uses the phase representative $0<\varphi <1$ specified here, corresponding to argument in $(0,2\pi )$. It implements that explicitly defined root, even for $\varphi >1/2$. The usual complex principal branch with argument in $(-\pi ,\pi ]$ would choose a different root on some eigenvalues. No substitution of that alternative branch is implicit.

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.

Let $L=2^n$ and use an $n$-qubit phase register. Begin with $|0^n⟩|v⟩$ and apply $H^{\otimes n}$ to the phase register, creating $L^{-1/2}{\sum }_{m=0}^{L-1}|m⟩|v⟩$. Controlled powers implementFor phase qubit $j$, counted from the most significant bit, the controlled power is $U^{2^{n-j}}$. It can be built from $2^{n-j}$ calls to the supplied controlled-$U$ gate. By quantum phase kickback, the phase-register state isApply the inverse quantum Fourier transform to obtain the exact quantum phase estimation mappingA computational basis measurement of the first register determines $y$ with certainty, hence $\varphi =y/L$ and the eigenvalue $e^{2\pi i\varphi }$. Exactness follows from the promised dyadic phase; no approximation or continued-fraction reconstruction is needed.

With only controlled-$U$ available as a query, the repeated-power construction uses $L-1=2^n-1$ oracle calls. The other Fourier-transform circuitry has polynomial size in $n$ in the ideal phase-gate model. The cost of exact phase estimation on a dyadic spectrum is therefore not polynomial in $n$ in this primitive-query model unless powered queries have additional implementations. The task does not require such a polynomial bound.

Use the unitary change of basis from part (b). It reduces the Feynman-Kitaev Hamiltonian to $A+B$, whereThe positive eigenvalues of $Q$ are positive integers, since its commuting ancilla projectors act on different qubits. The positive spectral gap of $B$ is $1-\cos (\pi /(T+1))\ge 2/(T+1{)}^2$. Thus both positive spectra are bounded below by $\mu \ge 2/(T+1{)}^2$, for $T\ge 1$.

Let $|u⟩=(T+1{)}^{-1/2}{\sum }_t|t⟩$, and split the work space into $G=\ker Q$ and $G^{\perp }$. The common ground space is $K=G\otimes |u⟩$. A unit vector in $\ker B$ orthogonal to $K$ has the form $|z⟩|u⟩$, where $z\in G^{\perp }$. Its orthogonal projection onto $\ker A$ simply removes its time-zero component, so the projected norm is $\sqrt{T/(T+1)}$. Consequently the smallest angle between two subspaces, after removing their common intersection, satisfiesThe Kitaev geometrical lemma now givesHere $1-\sqrt{1-x}\ge x/2$ supplies the penultimate step. HenceIf there are no input constraints, $Q=0$ and the propagation gap is already $\Omega (T^{-2})$, which is stronger.

The printed geometric-lemma notation needs a correction: the maximum overlap defines $\cos \vartheta$, not $\vartheta$, and is taken over normalized vectors in the two kernels, with the common ground space removed. The ground space restriction is essential when many quantum witnesses are allowed.

QMA is the class of promise problems for which a uniform polynomial-size quantum circuit checks a polynomial-size quantum witness. A YES instance has some quantum witness accepted with probability at least $2/3$; on a NO instance every quantum witness is accepted with probability at most $1/3$. The verifier initializes its own ancilla qubits to zero. Mixed quantum witnesses do not improve the maximum because acceptance is linear in their density operator.

For a verifier $U_x$ and quantum witness embedding $J|w⟩=|w⟩|0\cdots 0⟩$, define the witness acceptance operatorThe Rayleigh quotient gives ${\max }_{\Vert w\Vert =1}⟨w|F|w⟩={\lambda }_{\max }(F)$. Every entry of $F$ can be recomputed by the polynomial-storage path summation of part (b); no full quantum witness matrix needs to be stored.

Let the quantum witness have $m$ qubits, so its dimension is $D=2^m$. For this positive operator, trace-power witness optimization usesThis is the useful exponentiated form of the supplied logarithmic inequality. Positivity is required; such a logarithmic statement is not true for an arbitrary operator. No logarithms or roots need to be calculated.

Choose $d=2m+2$. On a NO instance,since $2^m(2/3{)}^{2m+2}=(4/9)(8/9{)}^m<1$. On a YES instance, $\mathrm{Tr}(F^d)\ge (2/3{)}^d>(1/2{)}^d$. Thus comparison of $2^d\mathrm{Tr}(F^d)$ with one distinguishes the cases using only multiplication and a final comparison.