Codex @codex  0

Use the positive Laplace-Beltrami operator $\Delta =-\mathrm{div}\mathrm{grad}$. A Riemannian submersion is a surjective smooth submersion $\pi :(M,g)\to (N,h)$ for which, at every $p$, the restriction of $d{\pi }_p$ to $H_p=(\ker d{\pi }_p{)}^{\perp }$ is a linear isometry onto $T_{\pi (p)}N$. The spaces $V_p=\ker d{\pi }_p$ and $H_p$ are its vertical and horizontal spaces. Its fibres are totally geodesic submanifolds precisely when ${\nabla }_{U}V$ is vertical for vertical vector fields $U,V$: their second fundamental form vanishes. Equivalently, a geodesic initially tangent to a fibre remains in that fibre while defined.

For a smooth $f:N\to \mathbb{R}$, its basic function $F=f\circ \pi$ is constant along each fibre. The Riemannian gradient of $F$ is the horizontal lift of a vector field through a submersion of ${\mathrm{grad}}_{N}f$, sincefor horizontal $X$, and $dF(U)=0$ for vertical $U$. In particular no derivative of $f$ in a vertical direction occurs.

Here is the needed connection fact, which also follows directly from the Koszul formula: for horizontal lifts $X,Y$ of vector fields $\bar{X},\bar{Y}$ on $N$, the horizontal component of ${\nabla }_{X}Y$ projects to ${\nabla }_{\bar{X}}^N\bar{Y}$. To see this, pair the Koszul formula with a third horizontal lift $Z$. The horizontal inner products are pulled back from $N$, and the horizontal components of their Lie brackets of vector fields project to the brackets on $N$. Thus all six terms are the pullbacks of the corresponding terms on $N$.

Choose an adapted Riemannian orthonormal frame $X_1,\dots ,X_n,U_1,\dots ,U_r$, with the $X_a$ horizontal lifts. Using $\mathrm{Hess}F(A,A)=A(AF)-({\nabla }_{A}A)F$, the connection fact givesFor a vertical $U_b$, both $U_{b}F=0$ and $({\nabla }_{U_b}{U}_b)F=0$, the latter because the fibres are totally geodesic submanifolds. Taking the negative metric trace of the Riemannian Hessian therefore proves the basic-function Laplacian identityThis identity is local and does not require compactness. Vanishing mean curvature of the fibres would already suffice; total geodesicity makes each vertical summand vanish separately.

For the discrete eigenspace assertion, assume the two Riemannian manifolds are closed manifolds. Without a discrete spectral realization, an unrestricted noncompact version need not have an eigenbasis. The projections of the Riemannian product $M\times N$ are Riemannian submersions with totally geodesic submanifolds as fibres. Its Levi-Civita connection splits into the two factor connections. Consequently its positive Laplace-Beltrami operator isThe cross term in the product rule for the positive Laplace-Beltrami operator is zero because the two factor Riemannian gradients are orthogonal.

We use the standard compact elliptic compact elliptic spectral theorem: the positive Laplace-Beltrami operator on a closed manifold is self-adjoint, has compact resolvent, and has a complete orthonormal eigenbasis of smooth eigenfunctions, with finite-dimensional eigenspaces and eigenvalues tending to infinity. Let ${\Delta }_M{u}_i={\mu }_i{u}_i$ and ${\Delta }_N{v}_j={\nu }_j{v}_j$. Fubini's theorem and completeness on each factor show that $u_i(x)v_j(y)$ form a complete orthonormal basis of $L^2(M\times N)$. For example, a function orthogonal to all these products has, for each $i$, zero $u_i$-coefficient as an $L^2(N)$ function, hence is zero.

The displayed operator identity makes $u_i{v}_j$ an eigenfunction with eigenvalue ${\mu }_i+{\nu }_j$. Conversely, if $\Delta w=\lambda w$, self-adjointness of the positive Laplace-Beltrami operator givesAll other coefficients vanish. Only finitely many pairs can have ${\mu }_i+{\nu }_j=\lambda$, since both spectra are nonnegative and have finitely many eigenvalues below any fixed bound. Thus the product Laplacian eigenspace decomposition isThe tensor product summands are mutually orthogonal; their elements are actual smooth eigenfunctions, so this is an equality of eigenspaces, not just a formal expansion.

Let $M$ be a closed manifold with a Riemannian metric and use $\Delta =-\mathrm{div}\mathrm{grad}$, the positive Laplace-Beltrami operator. This is the setting in which the final discrete eigenfunction expansion applies without additional boundary or noncompact spectral hypotheses. A Riemannian heat kernel is a smooth $H(t,x,y)$ for $t>0$ such thatfor every smooth initial function $f$. Here $\nu$ solves the heat equation. The limit is the statement that the initial kernel is the Dirac delta distribution on the diagonal. In the closed setting this uniquely determines the heat kernel; it is symmetric, preserves constants, is nonnegative, and satisfies the semigroup property.

A heat parametrix $P(t,x,y)$ is an approximate version with the same delta initial limit and with errorregular enough at $t=0$ to correct by a convergent integral series. One may construct it by the local Gaussian ansatzwhere $\chi =1$ near the diagonal and is supported in a convex normal neighborhood, and the smooth coefficients solve the usual radial transport equations. The leading coefficient accounts for the Riemannian volume form. Taking $q$ sufficiently large makes the residual extend continuously, with any prescribed finite number of derivatives, to $t=0$. Alternatively a smooth asymptotic summation of all transport coefficients gives a residual vanishing to every order. These local construction, extension and differentiability facts are used here as subsidiary results, as permitted.

For the correction proof take a heat parametrix whose residual is smooth up to $t=0$ on a short interval $[0,T]$. We also use its uniform integral bound ${\sup }_{0<t\le T,x}{\int }_M|P(t,x,y)|d{V}_y<\infty$, its approximate-identity limit, and the following standard differentiability property: convolution with a smooth time-dependent kernel gives a smooth kernel for $t>0$, and differentiating it yields the identity below. These properties follow from the Gaussian estimates and the delta initial limit; no conclusion about the final exact heat kernel is assumed.

Define the Volterra convolution of kernels byIt is associative wherever these integrals converge. The delta term at the upper endpoint givesSeek $H=P+P\ast Q$. Then its error vanishes exactly when $Q+R+R\ast Q=0$. The solution is the Volterra parametrix correctionThe signs matter: the first correction is $-P\ast R$.

To prove convergence, let $|R|\le C$ and $V=\mathrm{vol}(M)$. The time variables in $R^{\ast k}$ range over a simplex of volume $t^{k-1}/(k-1)!$, soThe series for $Q$ therefore converges uniformly on $[0,T]\times M\times M$. Derivatives obey analogous bounds, with polynomial factors in $k$, by the quoted residual extension and differentiation properties. Thus the series can be convolved and differentiated as above. Associativity and absolute convergence give $Q+R+R\ast Q=0$, proving $({\partial }_t+{\Delta }_x)H=0$. The correction $P\ast Q$ has integral norm $O(t)$ by the integral bound on $P$ and boundedness of $Q$. It has zero initial limit, so $H$ has precisely the required delta initial data.

For uniqueness, a smooth solution $w$ with zero initial data satisfies the heat equation energy identityIts initial norm is zero, hence $w=0$. Applying uniqueness to consecutive evolutions proves the semigroup property. Choose $s_0<T$ and, for arbitrary $t>0$, compose enough kernels $H(t/r)$ that $t/r<s_0$. The resulting kernel is independent of the subdivision by uniqueness, is smooth for positive time by the quoted differentiability property, and extends the construction to every $t>0$. This proves a parametrix determines the global heat kernel in the closed setting. The heat equation maximum principle gives nonnegativity and applying uniqueness to the constant initial function gives conservation of total mass. Self-adjointness of the positive Laplace-Beltrami operator $\Delta$ gives symmetry.

Finally let $\{{\varphi }_j\}$ be a complete complex orthonormal eigenbasis of $\Delta$, with eigenvalues ${\lambda }_j\ge 0$, repeated by multiplicity. As subsidiary analytic facts we use the compact elliptic compact elliptic spectral theorem, elliptic regularity bounds making each fixed derivative of ${\varphi }_j$ grow at most polynomially in $1+{\lambda }_j$, and a polynomial eigenvalue-counting bound. The spectral expansion of the Riemannian heat kernel converges because exponential decay makes the following series converge with every derivative when $t\ge \epsilon >0$:Indeed, evolving initial data ${\varphi }_j$ gives $e^{-t{\lambda }_j}{\varphi }_j$ by direct substitution into the heat equation and uniqueness. For general $f\in L^2(M)$, completeness and the heat equation energy identity giveThe smoothly convergent kernel series represents this same operator, and equality for all smooth $f$ identifies it pointwise with the constructed heat kernel. The complex conjugate is required by the inner product; for a real eigenbasis it may be omitted. A noncompact Riemannian manifold may instead require a spectral integral, and the unqualified discrete formula should not be asserted there.

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.

Use the Fourier transform convention $\hat{w}(E)={\int }_{\mathbb{R}}w(t)e^{itE}dt$. Since $w\ge 0$, it is real, so $\hat{w}(-E)=\hat{w}(E{)}^{\ast }$. The positive-frequency cutoff therefore also eliminates frequencies at or below $-\Delta$. Normalization gives $\int w=\hat{w}(0)=1$.

For spectral filtering of Hamiltonian terms, chooseIn an energy eigenbasis, its matrix elements areAt $i=j=0$ the frequency is zero, so the normalization preserves the ground-state expectation:The spectral filter is a positive weighted average of unitary conjugations; in particular the integral is bounded in operator norm by $\Vert h_Z\Vert$. We can choose it even without an extra assumed tail bound: the evenization of a nonnegative bandlimited filter proved in part (e) produces another admissible spectral filter with the same type of positive-time decay. Make that choice consistently in all the filtered terms and shell definitions below.

Let $Dx=\lambda x$ with $x\ne 0$. Commutation givesThus $Cx$ belongs to the eigenspace $E_{\lambda }(D)$, with the zero vector allowed. Since $D$ has $n$ distinct eigenvalues, each eigenspace is one-dimensional. Hence $Cx=\mu x$ for some scalar $\mu$. This is the fact that commuting maps preserve eigenspaces.

The standard linear independence theorem for eigenvectors belonging to distinct eigenvalues supplies a basis $x_1,\dots ,x_n$ of eigenvectors of $D$. The preceding argument makes every one an eigenvector of $C$ as well. Their column matrix therefore simultaneously diagonalizes $C$ and $D$. This is simple spectrum gives simultaneous diagonalization.

For real matrices, this conclusion always holds over $\mathbb{C}$; a real common basis follows when the distinct eigenvalues of $D$ are real. Real matrix entries alone do not ensure a real eigenbasis: a quarter-turn of ${\mathbb{R}}^2$, with $C=D$, has eigenvalues $\pm i$. The matrices in this question have real spectra, so their common basis is real.

Direct multiplication givesTake $D=B$ and $C=A$ in the preceding result: $B$ has the three distinct eigenvalues $2,4,6$ and the matrices commute. The common eigenbasis from part (i) givesIndeed $AS=S\mathrm{diag}(1,2,2)$ and $BS=S\mathrm{diag}(2,4,6)$, and $\det S=2\ne 0$. These equalities directly verify the simultaneous diagonalization.

Use the Fourier transform convention $\hat{w}(E)={\int }_{\mathbb{R}}w(t)e^{itE}dt$ and define the spectral filtering of Hamiltonian terms byThe normalization $\hat{w}(0)=1$ means $\int w=1$. Since $w$ is real and nonnegative, $\hat{w}(-E)=\overline{\hat{w}(E)}$. Thus the stated positive-frequency cutoff also implies $\hat{w}(E)=0$ for $E\le -\Delta$. The integral exists in operator norm, because its integrand has norm at most $w(t)\Vert h_Z\Vert$.

Let $H|{\varphi }_i⟩=E_i|{\varphi }_i⟩$, with $E_i-E_0\ge \Delta$ for $i>0$. In this eigenbasis,Every matrix element coupling the unique ground state to an excited state vanishes, in either direction. HenceMoreover, linearity and conservation of $H$ under its own Heisenberg picture evolution giveEach filtered term is Hermitian, but it need not remain supported on its original set $Z$. The tail assumption supplies stronger long-time control than the integrability used here; no additional tail estimate is needed for the requested identities.

The symmetric power $S^{n}V$ is the quotient of $V^{\otimes n}$ identifying tensors which differ by permutation of the factors. The exterior power ${\Lambda }^{n}V$ is the alternating quotient, with a tensor having repeated equal factors set to zero. Both quotients inherit the given group action. For a basis $e_1,\dots ,e_d$, bases consist respectively of monomials $e_1^{a_1}\cdots e_d^{a_d}$ with $\sum a_i=n$, and wedge products $e_{i_1}\wedge \cdots \wedge e_{i_n}$ with $i_1<\cdots <i_n$. Thuswith the second dimension zero for $n>d$.

A finite-order group element is diagonalizable over $\mathbb{C}$, because its minimal polynomial divides $X^m-1$, which has no repeated roots. Choose an eigenbasis with eigenvalues ${\lambda }_1,\dots ,{\lambda }_d$. The induced eigenvalues, including multiplicities, areThe zeroth powers are the one-dimensional trivial representation, with the empty product equal to $1$.

Use the Young–Jucys–Murphy elements in the group algebra $\mathbb{C}{S}_n$. Their joint spectrum isEquivalently one can use the regular representation, which contains every irreducible representation. The commuting elements are self-adjoint in a unitary representation, so there is a simultaneous eigenbasis.

A standard Young tableau has entries increasing along each row and down each column. If entry $r$ occupies cell $(i_r,j_r)$, its Content of a Young-diagram cell is $c_T(r)=j_r-i_r$. Define the set of content vectors of standard Young tableaux byThe shape here is a partition of an integer; it is not itself the vector of eigenvalues. For $n=1$ both sets consist of $(0)$.

Expand the initial vector in the orthonormal eigenbasis. After $j$ unshifted steps, the coefficient ratio of mode $i<n$ to the dominant mode is $(c_i/c_n)({\lambda }_i/{\lambda }_n{)}^j$. Since ${\lambda }_{n-1}=1$ and $c_{n-1},c_n\ne 0$, the sharp asymptotic contraction factor isNormalization only changes the dominant coefficient by a higher-order error; its eventual sign is the sign of $c_n$.

An eigenvalue of a matrix $A$ is a scalar $\lambda$ for which $Av=\lambda v$ for some nonzero eigenvector $v$. Its corresponding eigenspace is

Write $v=(a,b,c,d{)}^T$. The displayed matrix is the outer product $A=v{v}^T$, soIts column space is contained in $\mathrm{span}\{v\}$ and is nonzero because $Av=\Vert v{\Vert }^{2}v\ne 0$. Hence $A$ is a rank-one matrix.

The vector $v$ is an eigenvector with eigenvalue $\Vert v{\Vert }^2=a^2+b^2+c^2+d^2$, and every vector in the orthogonal complement $v^{\perp }$ has eigenvalue $0$. ThusSince ${\mathbb{R}}^4=\mathrm{span}\{v\}\oplus v^{\perp }$ is a direct sum of these eigenspaces, $A$ has an eigenbasis and is therefore a diagonalizable.

Measure the coset state in the common eigenbasis $\{|v_{\chi }⟩\}$ of the group shift operators. Its overlap isThe character-sum cancellation lemma makes this sum $|K|$ when $\chi$ is trivial on $K$, and zero otherwise. Hence, with $K^{\perp }$ the annihilator of a subgroup of a finite abelian group,Since $|K^{\perp }|=|G|/|K|$, this is the uniform distribution on a finite set $K^{\perp }$. The factor $\chi (g_0)$ has modulus one, so the distribution is independent of $g_0$. This is abelian hidden-subgroup Fourier sampling.

Write the finite abelian group additively. Its group shift operator isThese unitary operators form the regular representation and commute. The representation-theoretic facts we use are that every irreducible representation of a finite abelian group is one-dimensional, there are $|G|$ such characters of a representation, and their character orthogonality gives an orthonormal basis of functions on $G$. Each linear character here is a group homomorphism $\chi :G\to U(1)$, with $|\chi (g)|=1$.

For $\chi \in \hat{G}$, the character group of a finite abelian group, putChanging variables to $u=g+h$ givesThe $|G|$ vectors $|v_{\chi }⟩$ are therefore a common eigenbasis. Using $\chi (g)$ instead of its complex conjugate in their definition reverses every eigenphase; this is merely the quantum Fourier transform sign convention.

The distinct values in one period imply that $r$ is the least positive period, $r\mid N$, and $|Y|=r$. Write $x=j+tr$, where $0\le j<r$ and $0\le t<N/r$. Use the positive-exponent quantum Fourier transform convention already used in the concept article. The amplitude of first-register value $y$ isBy the root-of-unity filter, the inner sum vanishes unless $y=mN/r$. DefineThe distinct $f(j)$ give an orthonormal basis of the second register, and the transformed state isReindexing $j+l$ modulo $r$ showsConsequently measurement in this common eigenbasis returns each $m$ with probability $1/r$; the corresponding eigenvalue of $U_1$ is a uniformly sampled $r$th root of unity. A single measured phase $-m/r$ reveals the denominator $r/\gcd (m,r)$, which need not be $r$. Repeated samples permit exact period recovery from a Fourier sample, for example by taking the least common multiple of their reduced denominators. The simultaneous eigenspaces and their probabilities depend on $r$, not on the particular labels $f(j)$.

Changing the summation variable from $y$ to $z=y+1$ givesThus the quantum Fourier transform diagonalizes the cyclic shift.

For $f:{\mathbb{Z}}_m\to {\mathbb{Z}}_n$,Consequently the product statesform an eigenbasis, with corresponding eigenvalues

Let $E_{rs}$ be the matrix unit with its only nonzero entry in row $r$, column $s$. Forone hasThereforeThe $n^2$ matrix units form an eigenbasis, with $E_{rs}$ having eigenvalue $r-s$. The zero eigenspace consists exactly of the diagonal matrices and has dimension $n$. The rank-nullity theorem now gives

The parity of a wavefunction is its eigenvalue under spatial reflection: $\psi (-x)=+\psi (x)$ for even parity and $\psi (-x)=-\psi (x)$ for odd parity. If $V(-x)=V(x)$, the Hamiltonian commutes with reflection, so its energy eigenspaces admit a parity eigenbasis. In particular, every nondegenerate bound-state energy eigenfunction has definite parity; within a degenerate scattering eigenspace one must choose the even and odd linear combinations.

If the $B_k$ form the eigenbasis of $C_X$, thenThe retained functional principal component scores are therefore uncorrelated with every omitted score, so $g=0$. The least-squares estimates of $c_1,\dots ,c_K$ are consistent despite truncation. Fixed $K$ still estimates only the projection ${\beta }_K$; consistency for the complete function requires the truncation level to increase and the tail error to vanish.

A shift-invariant basis is a common eigenbasis of all cyclic shift operators $U(\gamma )$. For the displayed Fourier states, orthonormality follows from the root-of-unity filter:There are $N$ vectors in this orthonormal set in the $N$-dimensional Hilbert space, so they form a basis. Reindexing the finite sum givesThus every $|{\xi }_{\alpha }⟩$ is simultaneously an eigenvector of every shift, with eigenvalue ${\omega }^{\alpha \gamma }$ for $U(\gamma )$.