QuantumComputinglib learn · inspect · formalize
Checked on this commit 4,524 public declarations commit 4e2c5b2c0be5 Build record

Quantum Scientific Computing / Polynomial Signals

Polynomial Signals

Exercise PLAN-504

advanced · proof · test_public

primary source nl only. Additive answer successor with distinct content-bound review; historical rejected answer retained only for audit. No new Lean verification is claimed for this NL view.

Assumptions and conventions

  • A may be non-normal and non-diagonalizable; only its singular value decomposition is used
  • sigma_i are the singular values; v_i right singular vectors, w_i left singular vectors
  • Chebyshev first-kind polynomials T_m

Question

Let A be a general (not necessarily normal) n-qubit operator with singular value decomposition sending each right singular vector v_i to sigma_i times the left singular vector w_i. A block-encoding unitary U realizes A in its flag block, and the qubitization construction attaches to each index i a plane spanned by the flagged right vector and a flag-orthogonal companion, mapped by U to a plane spanned by the flagged LEFT vector and its companion. Prove that the EVEN qubitization output assembles into the operator sum over i of v_i times T_{2k}(sigma_i) times the adjoint of v_i, whereas the ODD output assembles into sum over i of w_i times T_{2k+1}(sigma_i) times the adjoint of v_i. Explain why the parity of the Chebyshev order decides whether the left factor is a right singular vector or a left singular vector.

Hint: mathematical route

Even paired qubitization has clean block V T_{2k}(Sigma)V†; an extra directed UR gives W T_{2k+1}(Sigma)V†. The guaranteed maps are U:H_i→H'_i and U†:H'_i→H_i; generally nonzero flag-orthogonal terms disappear only under the stated clean projection.

  1. From the exact clean block and an SVD, define directed companion matrices for U:H_i→H'_i and U†:H'_i→H_i, without a bidirectional-U or distinct-plane premise; singular-value1 is handled directly when leakage vanishes.
  2. R is diag(1,-1) in both companion bases, so E=U†RUR acts as M(sigma_i)^2 on the input plane; its kth power has clean coefficient T_{2k}(sigma_i) and generally nonzero orthogonal leakage.
  3. Projecting out leakage and summing over the right singular basis yields the even operator sum_i |v_i>T_{2k}(sigma_i)<v_i|.
  4. The extra directed UR has between-basis matrix M(sigma_i), so UR E^k has clean coefficient T_{2k+1}(sigma_i) on the left singular vector, with orthogonal leakage retained before projection.
  5. Projecting and summing gives the odd operator sum_i |w_i>T_{2k+1}(sigma_i)<v_i|; the zero-singular-value odd term vanishes, and the singular-value1 case needs no normalized leakage companion.
  6. Parity distinguishes return to the right plane from a final directed map to the left plane; it is not a universal-U swap rule or a deterministic clean-state statement, and no independent gate/error guarantee is inferred.
Reference solution

Interpret the requested output operators as the flag-to-flag blocks of the qubitization unitaries, as in the given block-encoding model, not as their entire output states. Let J map a system vector into the clean flag state, let P=JJ† and R=2P-I, and use the question's exact clean block J†UJ=A. Hence A is a contraction: ||Ax||≤||UJx||=||x||, so each singular value sigma_i lies in[0,1]. No normality, diagonalizability or Hermiticity of A is assumed. If one instead had only a scaled or approximate encoding of a separate target, this exact premise would have to be replaced by its actual normalized block and an error analysis; the present identities do not automatically transform the unscaled target.

Choose orthonormal right/left singular bases satisfying Av_i=sigma_i w_i and A†w_i=sigma_i v_i. Put p_i=Jv_i and p'_i=Jw_i. For sigma_i<1 let s_i=sqrt(1-sigma_i^2), q_i=(U†p'_i-sigma_i p_i)/s_i, and q'_i=(Up_i-sigma_i p'_i)/s_i. Exact compression and unitarity imply that q_i,q'_i are unit vectors orthogonal to the entire clean subspace. The directed input/output planes are H_i=span{p_i,q_i} and H'_i=span{p'_i,q'_i}. Directly, Up_i=sigma_i p'_i+s_i q'_i, Uq_i=s_i p'_i-sigma_i q'_i, U†p'_i=sigma_i p_i+s_i q_i, U†q'_i=s_i p_i-sigma_i q_i. Thus U:H_i→H'_i and U†:H'_i→H_i have matrix S_i=[[sigma_i,s_i],[s_i,-sigma_i]] between the respective ordered bases. These are directed maps: U† is the guaranteed reverse, not a second application of U. The geometric planes may coincide or differ; neither possibility changes the argument.

In both bases R has matrix Z=diag(1,-1). Set E=U†RUR and M_i=S_i Z=[[sigma_i,-s_i],[s_i,sigma_i]]. Then E on H_i has matrix M_i^2, while UR maps H_i to H'_i with matrix M_i. Matrix powers of this scalar rotation yield the full state actions E^k p_i=T_{2k}(sigma_i)p_i+s_i U_{2k-1}(sigma_i)q_i for k≥1, UR E^k p_i=T_{2k+1}(sigma_i)p'_i+s_i U_{2k}(sigma_i)q'_i for k≥0. Here U_l is the second-kind polynomial, not the bare encoding unitary U, and E^0 p_i=p_i is treated directly. These formulas retain the generally nonzero flag-orthogonal components. In the odd expression the R in the extra UR matters on the leakage; it must not simply be omitted as an unqualified extra-U step.

Applying the clean projection kills q_i and q'_i. For every right basis vector v_i we therefore have J†E^kJ v_i=T_{2k}(sigma_i)v_i, J†UR E^kJ v_i=T_{2k+1}(sigma_i)w_i. Since the right vectors form an orthonormal basis, summing these actions gives exactly J†E^kJ=sum_i |v_i>T_{2k}(sigma_i)<v_i|, J†UR E^kJ=sum_i |w_i>T_{2k+1}(sigma_i)<v_i|. The even paired iterate returns to the right plane and therefore uses v_i on both outer sides; the odd extra directed UR ends in the left plane and therefore uses w_i on the left. Parity refers to the alternating U/U† construction and this final directed map, not a universal rule that each U swaps two planes in both directions.

At sigma_i=0 the companion definitions remain valid because s_i=1; the odd clean coefficient is T_{2k+1}(0)=0, whereas the even one is T_{2k}(0)=(-1)^k. At sigma_i=1 the leakage norms vanish, so use the direct relations Up_i=p'_i and U†p'_i=p_i instead of dividing by zero. The input clean line is fixed by E, while UR sends it to the output clean line: the output line need not be the same physical line as the input. This gives the required even and odd coefficients1. It does not assert that an arbitrary general U has a single invariant clean eigenvector. Degenerate singular-vector choices do not change the resulting operators: T_{2k}(x)=h_k(x^2) yields the even operator h_k(A†A), while T_{2k+1}(x)=x g_k(x^2) yields the odd operator A g_k(A†A). In particular the odd expression does not depend on arbitrary choices of zero-singular-value left vectors.

Two concrete checks explain why the corrected distinctions matter. In the source's Eq10.79 example, write the ambient ancilla-major basis as e_0,e_1,e_2,e_3. Its U sends e_0→e_2→e_1→e_0 and fixes e_3. For the nilpotent A and its sigma=0 right/left vectors, H=span{e_0,e_2} and H'=span{e_1,e_2}; yet U(H')=span{e_0,e_1} is not H. This disproves a bidirectional-U swap claim, while U† does return H' to H. Separately let A=(1/2)I_2 and use U=[[tI_2,sI_2],[sI_2,-tI_2]], t=1/2, s=sqrt(3)/2. For a unit v, E Jv=(-1/2)|0>v+(sqrt(3)/2)|1>v. Its clean block coefficient is -1/2, but the unit-norm output is not the unnormalized clean vector alone.

For a unit input psi, if F denotes one of the displayed clean operators, projection succeeds with probability ||F psi||^2; conditional on a nonzero success the normalized system vector is F psi/||F psi||. A zero branch cannot be normalized. The clean block identities alone are not deterministic state-preparation claims or gate/query/finite-bit resource bounds. The source's general-SVD derivation in Eqs10.61–10.75 establishes these formulas directly; its printed Proposition10.7 has a narrower Hermitian label, which is not used as an extra premise or as a substitute for this derivation.

Assessment rubric
  1. Uses E=U†RUR and shows the even clean block is sum_i |v_i>T_{2k}(sigma_i)<v_i|, with exact flag-to-flag meaning rather than a whole-state equality. (1 weight)
  2. Tracks the extra directed UR, not an omitted R or a bidirectional-U swap, and obtains the odd clean block sum_i |w_i>T_{2k+1}(sigma_i)<v_i|. (1 weight)
  3. Derives the scalar amplitudes from the directed singular-vector companion matrices, retains flag-orthogonal leakage before projection, and treats k=0 and sigma_i=0/1 without zero-norm division. (1 weight)
  4. Explains parity by alternating U:H_i→H'_i and U†:H'_i→H_i with a final UR for odd order; distinguishes right-even from balanced-odd assembly without requiring distinct geometric planes or making deterministic full-output, scaling/error or resource claims. (1 weight)

Graph-grounded prerequisites

These are support/source anchors, not an automatically extracted theorem-dependency proof.

  • Block Encoding (frozen curriculum anchor; see graph registry)
  • Matrix analysis (frozen curriculum anchor; see graph registry)
  • Operator and semigroup theory (frozen curriculum anchor; see graph registry)

Download the shared core registry

Lean view and evidence boundary

No aligned Lean view has been admitted for this native core in this textbook checkpoint. This is not a claim that the underlying mathematics lacks a formal proof. Contributions must preserve this exact target contract and split.

Read the Lean track and existing diagnostic warmups

Sources and provenance

Contract: 8ad3d619d213308a8f49f558407d3deeb56a1335837f03a1bc637c6af949ab26

Item digest: 92868fd70b2d449319fb8c5e9c10064c9f85cad43b5a9e6c45372d19207914ca

Graph snapshot: qcb-snapshot-ffa644f0da3031a1ce754b9a

Authored support steps: 6. This is not measured graph distance or model difficulty.

Answer revision lineage

{ "predecessor_item_digest": "c381cbd3f7121acaa8f6a1c1901efc4466ddd3eb9446022ad1a5cef85a9c92ea", "predecessor_key": "qcb-core-plan-504-svt-parity-left-right-vectors", "repair_record": "shared/repairs/plan-504-projected-svt-direction-repair.json", "repair_record_digest": "a3c2b4ca1fdb90c64ae228ce9af51f3f4886b61dc9cadf8b9436f12d1575f52f", "review_record": "reports/textbook-v1/repairs/plan-504-projected-svt-direction-repair-review.json", "review_digest": "072f8c4d107365253414585c695b6d94da57844a70baf1ca6207b90d78f565dd" }