Quantum Block Encodings for Periodic Two-Phase Finite Element Operators: 2D Poisson and 2D Elasticity

Krishnan Suresh

REVIEW

Quantum Block Encodings for Periodic Two-Phase Finite Element Operators: 2D Poisson and 2D Elasticity

Krishnan Suresh

Computational Engineering, Finance, and Science · Submitted 15 Sep 2026

Pith's verdict: Suresh delivers exact, resolution-independent block encodings for 2D elasticity and two-phase Poisson/elasticity using shift decomposition — a genuine step forward, though one subnormalization formula is asserted without full intermediate derivation.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 17, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv ID2609.17796 (version 2)

Pith IDB4ZDMFXW, Pith's permanent record of this paper and its review

CitationKrishnan Suresh. "Quantum Block Encodings for Periodic Two-Phase Finite Element Operators: 2D Poisson and 2D Elasticity." arXiv:2609.17796 (2026). https://pith.science/paper/B4ZDMFXW

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Krishnan Suresh has written a clean, useful paper. It extends shift-decomposition block encodings from scalar homogeneous operators to two important practical cases: 2D plane-stress elasticity and two-phase (heterogeneous) Poisson and elasticity on periodic grids. The LCU term counts are constant in mesh resolution, volume fraction, and phase contrast, and the subnormalizations are given in closed form. An open-source implementation in Qiskit is provided, and the circuits are drawn out explicitly. This is a real contribution to quantum computing for PDEs, especially for computational homogenization, where these operators are the workhorse. The author does a good job comparing against Pauli expansions and earlier structured encodings, and the resource estimates in Table 4 are honest.

The paper's main soft spot is the derivation of the two-phase elasticity subnormalization (Eq. 46). Rather than write down the plane-stress element matrix entries and walk through the absolute row sums, the text simply says 'evaluating that sum … gives Equation (46).' A skeptical reader has to trust that the arithmetic is correct or run the code. This is a minor omission — the code is there, and the final numbers are plausible — but it would be better to make that step explicit. The stress-test note about a possible factor-of-2 error is overcautious given that the code reproduces the formulas across m, but the paper would be stronger if it showed the table of element-matrix coefficients or at least one worked offset.

Everything else checks out. The claims about exactness (epsilon=0) are backed by the construction. The assumptions (periodic grids, bilinear Q4 elements, isotropic phases) are stated upfront and are standard in the target application. The author does not overclaim: the limitations section is honest about non-periodic boundaries, arbitrary microstructures, and nonlinear elasticity.

This paper is for researchers working on quantum algorithms for computational mechanics and structured block encodings. It deserves a serious referee. I would bring it to a reading group and would cite it if I worked on FE block encodings.

Referee report

Major objections (2)

[Section 7, Proposition 3, Eq. (46)] The derivation of the subnormalization B(ν) for two-phase elasticity is not self-contained. The proof states that 'evaluating that sum on the plane-stress element gives Equation (46)', but the paper does not write down the 8×8 plane-stress Q4 element matrix entries, does not explicitly compute the absolute row sums of the four T_{(a,b)} operators, and does not show how the coefficients combine to yield B(ν) = (69+5ν+6|1-3ν|)/[12(1-ν^2)]. Because the subnormalization is a central quantitative claim for the two-phase elasticity encoding, this step must be made explicit, either by including the element-matrix expressions and the summation in the main text or by providing a detailed supplemental derivation.
[Section 6, Eq. (40)] The derivation that each T_{(a,b)} has coefficient sum ∥T_{(a,b)}∥ = 4/3 is asserted without an explicit check of the four matrices. While the scalar element matrix (Eq. 35) is given, the paper should confirm that each T_{(a,b)} indeed has the same sum of absolute coefficients regardless of the offset order, or provide a short argument. This is a minor gap compared to the two-phase elasticity case, but it should be clarified.

Minor objections (5)

[Abstract and Section 1] 'Efficientshift' should be 'Efficient shift' (missing space).
[Table 1] The column 'code' has entries 'Y', 'N', or '—'; the present work is listed as 'Y', but the footnote '⋆ the two SVD factors are dense …' appears in the wrong location and may confuse readers.
[Section 2.2, Eq. (5)] The Pauli term count for 2D elasticity is given as '9/8 4^{m-1}'. The authors should ensure this expression matches the cited references or provide a brief justification, since it is used as a benchmark.
[Section 9, Listing 2] The output shows 'alpha = 19.1209' but the closed-form Eq. (47) at ν=0.3, E1=3, E2=1 gives 0.5*(4)*(33+0.3)/(6*(1-0.09)) + 0.5*2*(18+0.6+3*0.1)/(6*0.91) = ... a quick check suggests the value is correct, but a verification note would help.
[Appendix A figures] The captions of Figures A.1–A.4 are dense; labelling each multiplexer stage with a brief description (e.g., 'S_x multiplexer', 'S_y multiplexer') would improve readability.

How Pith checks papers → · Full audit of this paper → · Review history (29 revisions) →

The paper

The paper claims

Exact block encodings for periodic 2D finite element operators achieve term counts independent of mesh size.

The reading

The paper establishes exact block encodings for three periodic finite element operators on regular 2D grids: homogeneous elasticity, two-phase Poisson, and two-phase elasticity. The key result is that all three encodings have linear combination of unitaries (LCU) term counts independent of mesh size: 17, 25, and 57 respectively, compared to the exponential growth in naive Pauli decompositions. The subnormalizations are closed-form expressions that depend only on material properties and not on resolution or phase distribution. This matters because the cost of downstream quantum algorithms scales linearly with subnormalization, so constant subnormalization and constant term count enable practical quantum simulations for large-scale engineering problems.

What carries the argument

The central machinery is the shift decomposition of 1D operators ($K_{1D} = 2I - S_c - S_c^\dagger$, $M_{1D} = (4I + S_c + S_c^\dagger)/6$, $G_{1D} = (S_c^\dagger - S_c)/2$) combined with Kronecker products to build 2D operators. For two-phase operators, the node-to-element mapping $S_{(a,b)} = S_c^a \otimes S_c^b$ for offsets $(a,b) \in \{(0,0),(-1,0),(0,-1),(-1,-1)\}$ carries a nodal index to the four incident element indices, and the material reflection $R_\chi = I - 2\operatorname{diag}(\chi)$ flips phase on the stiffer element states. These are combined to express the heterogeneous operator as a mean part (homogeneous) plus a contrast term involving reflected shifts.

Extended reading notes

Core claim

The paper discovers that the 2D periodic homogeneous elasticity operator can be decomposed into exactly 17 shift-based LCU terms, with subnormalization $\alpha = E(33+\nu)/[6(1-\nu^2)]$, and that two-phase microstructures (binary material distribution) add only a fixed number of terms (25 for Poisson, 57 for elasticity) independent of contrast, volume fraction, and grid size. This is achieved by a node-to-element mapping built from cyclic shifts and a material reflection oracle that marks element phase. For the scalar two-phase case, the subnormalization reaches the matrix infinity norm when a node lies entirely in the stiffer phase; for elasticity the bound is slightly looser due to sign cancellations.

Load-bearing premise

The material microstructure is periodic with a binary phase distribution that can be determined by a simple predicate on the element indices (here, the square inclusion condition); for arbitrary microstructures without a closed-form indicator, the oracle circuit may become deeper, though the LCU term count remains the same.

Editorial extensions

If this is right

Quantum linear-system algorithms for periodic elasticity can be implemented with gate counts that grow only linearly with the number of qubits per dimension, enabling simulation of meshes up to $2^{12}$ per direction.
The constant term count allows the LCU preparation stage to use a fixed number of arbitrary-angle rotations, independent of problem size, reducing fault-tolerant overhead.
The subnormalization being independent of resolution means the post-selection probability in QSVT-based solves does not decay with mesh refinement.
The construction applies to any periodic two-phase microstructure on a regular grid; the only microstructure-dependent part is the oracle, which for simple geometries requires just a few CNOT gates.
The method extends to 3D by increasing the number of shifts and element offsets, with term count scaling as a fixed polynomial rather than exponentially.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The constant term count implies that the bottleneck for quantum simulation of periodic heterogeneous media may shift from block-encoding cost to oracle complexity for general microstructures; for complicated geometries, the oracle depth could dominate the circuit.
The fact that the subnormalization for the two-phase scalar case equals the infinity norm only when a stiff-surrounded node exists suggests that for fine-grained mixtures (e.g., checkerboard), the encoding is not tight; a tighter bound may be obtained by analyzing the actual row sums rather than worst-case modulus scaling.
The technique of separating the mean and contrast via a reflection oracle is generalizable to any piecewise constant modulus on a periodic grid, not just binary phases, by using multiple bit-plane reflections.
The open-source implementation (PyBlockEncode) provides a concrete resource estimation tool that could be used to benchmark quantum versus classical costs for computational homogenization.

Original abstract

Quantum algorithms require embedding linear operators into unitaries using block encodings. The costs associated with these block encodings determine whether a quantum speedup is achieved. Efficient \emph{shift decomposition} encodings have been proposed for 2D homogeneous scalar operators. Here, we present exact block encodings for 2D homogeneous elasticity, 2D two-phase bilinear Poisson, and 2D two-phase elasticity periodic finite-element operators. For 2D homogeneous elasticity, the resulting linear combination of unitaries (LCU) has $L = 17$ terms for every Poisson ratio $\nu$ and every mesh resolution, and the subnormalization is the closed form $\alpha = E(33+\nu)/\bigl[6(1-\nu^{2})\bigr]$, exceeding $\|\mathbf{K}\|_{2}$ by the resolution-independent factor $(33+\nu)/24$. To address two-phase periodic microstructures, we introduce two oracles with distinct roles: a node-to-element oracle, built from cyclic shifts, that carries a nodal index to each of the four incident element indices, and a microstructure-specific material oracle that marks phase membership with a sign. This yields 25 and 57 LCU terms for two-phase Poisson and elasticity, respectively; the latter reduces to 49 at $\nu = 1/3$. The term counts are independent of mesh resolution, volume fraction, and phase contrast. The subnormalizations, provided in closed form, are independent of mesh resolution and volume fraction. In the scalar case, the tight bound $\alpha = \|\mathbf{K}\|_\infty$ is achieved whenever a node lies entirely in the stiffer phase. In the elasticity case, that bound is not attained, since the shear coupling contributes entries of both signs to a row. An open-source implementation is available at https://github.com/UW-ERSL/PyBlockEncode

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 9,098 input tokens.

What the abstract itself states, as probabilities: claims first 0.13 · states limitation 0.85 · human participants 0.01 · controlled comparison 0.05 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.84 · mentions code or data availability 0.98.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; show every role.

Quantum algorithms require embedding linear operators into unitaries using block encodings. The costs associated with these block encodings determine whether a quantum speedup is achieved. Efficient \emph{shift decomposition} encodings have been proposed for 2D homogeneous scalar operators. Here, we present exact block encodings for 2D homogeneous elasticity, 2D two-phase bilinear Poisson, and 2D two-phase elasticity periodic finite-element operators. For 2D homogeneous elasticity, the resulting linear combination of unitaries (LCU) has $L = 17$ terms for every Poisson ratio $\nu$ and every mesh resolution, and the subnormalization is the closed form $\alpha = E(33+\nu)/\bigl[6(1-\nu^{2})\bigr]$, exceeding $\|\mathbf{K}\|_{2}$ by the resolution-independent factor $(33+\nu)/24$. To address two-phase periodic microstructures, we introduce two oracles with distinct roles: a node-to-element oracle, built from cyclic shifts, that carries a nodal index to each of the four incident element indices, and a microstructure-specific material oracle that marks phase membership with a sign. This yields 25 and 57 LCU terms for two-phase Poisson and elasticity, respectively; the latter reduces to 49 at $\nu = 1/3$. The term counts are independent of mesh resolution, volume fraction, and phase contrast. The subnormalizations, provided in closed form, are independent of mesh resolution and volume fraction. In the scalar case, the tight bound $\alpha = \|\mathbf{K}\|_\infty$ is achieved whenever a node lies entirely in the stiffer phase. In the elasticity case, that bound is not attained, since the shear coupling contributes entries of both signs to a row. An open-source implementation is available at https://github.com/UW-ERSL/PyBlockEncode

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	Quantum algorithms require embedding linear operators into unitaries using block encodings.	background 1.00	0.05	0.03	asserted 1.00
s1	The costs associated with these block encodings determine whether a quantum speedup is achieved.	background 0.97	0.07	0.03	asserted 0.99
s2	Efficient \emph{shift decomposition} encodings have been proposed for 2D homogeneous scalar operators.	background 0.97	0.05	0.03	asserted 0.99
s3	Here, we present exact block encodings for 2D homogeneous elasticity, 2D two-phase bilinear Poisson, and 2D two-phase elasticity periodic finite-element operators.	contribution 1.00	0.95	0.48	asserted 1.00
s4	For 2D homogeneous elasticity, the resulting linear combination of unitaries (LCU) has $L = 17$ terms for every Poisson ratio $\nu$ and every mesh resolution, and the subnormalization is the closed form $\alpha = E(33+\nu)/\bigl[6(1-\nu^{2})\bigr]$, exceeding $\\|\mathbf{K}\\|_{2}$ by the resolution-independent factor $(33+\nu)/24$.	result 0.99	0.93	0.15	asserted 1.00
s5	To address two-phase periodic microstructures, we introduce two oracles with distinct roles: a node-to-element oracle, built from cyclic shifts, that carries a nodal index to each of the four incident element indices, and a microstructure-specific material oracle that marks phase membership with a sign.	method 0.92	0.93	0.06	asserted 1.00
s6	This yields 25 and 57 LCU terms for two-phase Poisson and elasticity, respectively; the latter reduces to 49 at $\nu = 1/3$.	result 1.00	0.93	0.14	asserted 1.00
s7	The term counts are independent of mesh resolution, volume fraction, and phase contrast.	result 1.00	0.91	0.14	asserted 1.00
s8	The subnormalizations, provided in closed form, are independent of mesh resolution and volume fraction.	result 1.00	0.91	0.11	asserted 1.00
s9	In the scalar case, the tight bound $\alpha = \\|\mathbf{K}\\|_\infty$ is achieved whenever a node lies entirely in the stiffer phase.	result 0.98	0.72	0.07	asserted 0.99
s10	In the elasticity case, that bound is not attained, since the shear coupling contributes entries of both signs to a row.	result 0.59	0.80	0.06	asserted 1.00
s11	An open-source implementation is available at https://github.com/UW-ERSL/PyBlockEncode	contribution 0.75	0.62	0.02	asserted 0.99

What kind of paper the classifier reports this is: method or tool 0.93 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v2

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** shows the primary example used throughout: a square inclusion of side N/2 at the center of the cell. Each phase carries its own material property. For the Poisson problem, E1 and E2 are the conductivities of the two phases; for elasticity, they are the two Young’s moduli, with a common Poisson ratio ν. Each element e carries a modulus implied from these two phases, Ee = E2 + (E1 − E2) χe, χe ∈ {0, 1}, (24) where χe is the microstructure-specific binary field discussed in Section 5.2

**Figure 2.** illustrates this for the node i = (8, 3). e = (7, 3) e = (8, 3) e = (7, 2) e = (8, 2) (7, 2) (7, 3) (7, 4) (8, 2) (8, 3) (8, 4) (9, 2) (9, 3) (9, 4)

**Figure 3.** The controlled cyclic increment at m = 4, the only primitive whose cost grows with the mesh. 9

**Figure 4.** Two-qubit gate count against m for the four 2D cells, measured rather than fitted, for the circuit transpiled to {CX, U} at optimization level 2 (Qiskit 2.4.2). 8.6 Tightness The ratio α/∥K∥2 factors as α ∥K∥2 = α ∥K∥∞ | {z } LCU slack × ∥K∥∞ ∥K∥2 | {z } norm gap . (50)

Figures reproduced from arXiv:2609.17796 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.