Quantum methods encode two phase materials for faster simulations
Quantum Block Encodings for Periodic Two-Phase Finite Element Operators: 2D Poisson and 2D Elasticity
Computational Engineering, Finance, and Science
Summary
Simulating how materials like metals or composites behave under forces needs solving complex math problems. The authors developed a way to translate important material properties into quantum computing language exactly, which could make certain material simulations more efficient. They designed special quantum routines that work well even for materials made of two different parts with repeating patterns. Their method ensures the complexity stays manageable no matter how finely the material is divided for simulation.
What this means in practice
- •For quantum algorithm developers: Implement exact quantum operators for simulating complex material behaviors with two-phase periodic structures using finite element models.
- •For quantum simulation engineers: Optimize quantum resources for elastic and Poisson problem solvers in materials with repeating microstructures, ensuring predictable computational costs.
Authors
Krishnan Suresh
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 $ν$ and every mesh resolution, and the subnormalization is the closed form $α= E(33+ν)/\bigl[6(1-ν^{2})\bigr]$, exceeding $\|\mathbf{K}\|_{2}$ by the resolution-independent factor $(33+ν)/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 $ν= 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 $α= \|\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