Block-encoding operators are one of the essential components in quantum algorithms based on Quantum Signal Processing. Their gate complexity largely determines the overall gate complexity of the full algorithm. Using variational methods, we compile single-ancilla block-encoding unitaries with near-optimal resource requirements for a large range of input matrices. We find that the number of variational parameters in the parameterized quantum circuit approaches the number of free parameters in the input matrices, depending on whether they are real, complex and/or hermitian. Additionally, symmetries present in the input matrix can be incorporated into the ansatz circuit, reducing the parameter count further and enhancing optimizability. While performing a variational compilation of block-encodings ceases to be computationally feasible for large system sizes, the constructed operators can be used as components of larger block-encodings via a linear combination of block-encodings.
Block-encoding operators are one of the essential components in quantum algorithms based on Quantum Signal Processing. Their gate complexity largely determines the overall gate complexity of the full algorithm. Using variational methods, we compile single-ancilla block-encoding unitaries with near-optimal resource requirements for a large range of input matrices. We find that the number of variational parameters in the parameterized quantum circuit approaches the number of free parameters in the input matrices, depending on whether they are real, complex and/or hermitian. Additionally, symmetries present in the input matrix can be incorporated into the ansatz circuit, reducing the parameter count further and enhancing optimizability. While performing a variational compilation of block-encodings ceases to be computationally feasible for large system sizes, the constructed operators can be used as components of larger block-encodings via a linear combination of block-encodings.
One of the challenges in quantum computing is the synthesis of unitary operators into quantum circuits with polylogarithmic gate complexity. Exact synthesis of generic unitaries requires an exponential number of gates in general. We propose a novel approximate quantum circuit synthesis technique by relaxing the unitary constraints and interchanging them for ancilla qubits via block encodings. This approach combines smaller block encodings, which are easier to synthesize, into quantum circuits for larger operators. Due to the use of block encodings, our technique is not limited to unitary operators and can be applied for the synthesis of arbitrary operators. We show that operators which can be approximated by a canonical polyadic expression with a polylogarithmic number of terms can be synthesized with polylogarithmic gate complexity with respect to the matrix dimension.
Implementing polynomial functions of Hermitian matrices on quantum hardware is a foundational task in quantum computing, critical for accurate Hamiltonian simulation, quantum linear system solving, high-fidelity state preparation, machine learning kernels, and other advanced quantum algorithms. Existing state-of-the-art techniques, including Qubitization, Quantum Singular Value Transformation (QSVT), and Quantum Signal Processing (QSP), rely heavily on block-encoding the Hermitian matrix. These methods are often constrained by the complexity of preparing the block-encoded state, the overhead associated with the required ancillary qubits, or the challenging problem of angle synthesis for the polynomial's phase factors, which limits the achievable circuit depth and overall efficiency. In this work, we propose a novel and resource-efficient approach to implement arbitrary polynomials of a Hermitian matrix by leveraging the Generalized Quantum Signal Processing (GQSP) framework. Our method circumvents the need for block-encoding and avoids the compounding post-selection overheads characteristic of LCU-based constructions, achieving a stable, degree-independent success probability. We derive closed-form expressions for symmetric polynomial expansions and demonstrate how linear combinations of GQSP circuits can realize the desired transformation. This approach reduces resource overhead and opens new pathways for quantum algorithm design for functions of Hermitian matrices, particularly in settings where the Hermitian operator arises naturally from symmetric combinations of unitaries.
In quantum computing, every operation must be unitary to ensure reversibility and preserve probability amplitudes, which is a fundamental requirement for any quantum process. However, many real-world computations in areas such as machine learning, signal processing, and dynamical modeling inherently involve non-unitary matrices. Since these cannot be directly implemented on quantum hardware, a practical workaround is required. One promising solution is block-encoding, a framework that embeds a non-unitary matrix A in C^{n x n} into a larger unitary matrix U_A in C^{2n x 2n}, enabling the simulation of non-unitary actions within quantum circuits through ancilla preparation and post-selection. This thesis presents a novel and efficient method for block-encoding based on Gram-Schmidt orthogonalization. The algorithm constructs a unitary U_A with A embedded as the top-left block and remaining components filled to ensure global unitarity. Compared to traditional Singular Value Decomposition (SVD)-based approaches, the proposed method achieves at least 20% improvement in execution time, without compromising output fidelity or numerical stability. A corresponding quantum circuit is developed to extract A|psi> (that is, the transformed version of the input quantum state |psi> after applying the linear operator A) using post-selection on ancilla qubits, with the success probability determined by
Quantum signal processing combined with quantum eigenvalue transformation has recently emerged as a unifying framework for several quantum algorithms. In its standard form, it consists of two separate routines: block encoding, which encodes a Hamiltonian in a larger unitary, and signal processing, which achieves an almost arbitrary polynomial transformation of such a Hamiltonian using rotation gates. The bottleneck of the entire operation is typically constituted by block encoding, and in recent years, several problem-specific techniques have been introduced to overcome this problem. Within this framework, we present a procedure to block-encode a Hamiltonian based on its matrix product operator (MPO) representation. More specifically, we encode every MPO tensor in a larger unitary of dimension D+2, where $D=\ensuremath{\lceil}{log}_{2}(\ensuremath{\chi})\ensuremath{\rceil}$ is the number of subsequently contracted qubits that scales logarithmically with the virtual bond dimension $\ensuremath{\chi}$. Given any system of size L, our method requires L+D ancillary qubits in total, while the number of one- and two-qubit gates decomposing the block-encoding circuit scales as $O(L{\ensuremath{\chi}}^{2})$. Moreover, block encodings generally require postselection. Our model exhibits a success probability decaying exponentially with L. Nevertheless, we present a way to mitigate this limitation for finite systems.
Quantum algorithms based on Quantum Signal Processing (QSP) offer the potential for speedups across a broad range of applications, with block encodings serving as the central input model. In this framework, non-unitary matrices are embedded into larger unitary operators, and the cost of constructing these encodings often dominates the overall gate complexity. In this work, we introduce Tag-and-Restore Encoding (TARE), a block-encoding method for linear combinations of Pauli strings. In this method coefficient magnitudes are absorbed into a unitary built from a set of mutually anti-commuting Pauli strings acting on the system register. These Pauli strings are then mapped to the target Pauli strings through appropriate transformations, yielding a block encoding of the target operator. The ancilla register size scales logarithmically with the number of Pauli strings and can be extended to larger registers providing a width/depth tradeoff. We evaluate TARE through numerical simulations of the transverse-field Ising model, the Jordan-Wigner image of a fermionic star Hamiltonian, and random Pauli-string operators. Compared with standard Linear Combination of Unitaries (LCU), TARE substantially reduces the T-gate count while improving circuit depth in several cases. These results suggest that TARE can provide resource-efficient block encodings for a wide range of relevant systems.
Quantum circuits naturally implement unitary operations on input quantum states. However, non-unitary operations can also be implemented through "block encodings", where additional ancilla qubits are introduced and later measured. While block encoding underlies a host of well-established theoretical applications, its implementation in practice has been prohibitively expensive for current quantum hardware. In this paper, we present practical and explicit block encoding circuits implementing matrix polynomial transformations of a target matrix. With standard approaches, block-encoding a degree-d matrix polynomial requires a circuit depth scaling as d times the depth for block-encoding the original matrix alone. By leveraging the recently introduced Fast One-Qubit Controlled Select Linear Combination of Unitaries (FOQCS-LCU) framework, we show that the additional circuit-depth overhead required for encoding matrix polynomials can be reduced to scale linearly in d with no dependence on system size or the cost of block encoding the original matrix. Moreover, we demonstrate that the FOQCS-LCU circuits and their associated matrix polynomial transformations can be controlled with negligible overhead, enabling efficient applications such as Hadamard tests. Finally, we provide explicit hardware-efficient circuits for representative non-unitary operators, namely quantum spin models, together with detailed non-asymptotic gate counts and circuit depths.