Title: Explicit block encodings of boundary value problems for many-body elliptic operators
Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.  more » « less
Award ID(s):
2120757
PAR ID:
10678970
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Quantum
Date Published:
Journal Name:
Quantum
Volume:
9
ISSN:
2521-327X
Page Range / eLocation ID:
1764
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Some linear integro-differential operators have old and classical representations as the Dirichlet-to-Neumann operators for linear elliptic equations, such as the 1/2-Laplacian or the generator of the boundary process of a reflected diffusion. In this work, we make some extensions of this theory to the case of a nonlinear Dirichlet-to-Neumann mapping that is constructed using a solution to a fully nonlinear elliptic equation in a given domain, mapping Dirichlet data to its normal derivative of the resulting solution. Here we begin the process of giving detailed information about the Lévy measures that will result from the integro-differential representation of the Dirichlet-to-Neumann mapping. We provide new results about both linear and nonlinear Dirichlet-to-Neumann mappings. Information about the Lévy measures is important if one hopes to use recent advancements of the integro-differential theory to study problems involving Dirichlet-to-Neumann mappings. 
    more » « less
  2. Grid-free Monte Carlo methods based on thewalk on spheres (WoS)algorithm solve fundamental partial differential equations (PDEs) like the Poisson equation without discretizing the problem domain or approximating functions in a finite basis. Such methods hence avoid aliasing in the solution, and evade the many challenges of mesh generation. Yet for problems with complex geometry, practical grid-free methods have been largely limited to basic Dirichlet boundary conditions. We introduce thewalk on stars (WoSt)algorithm, which solves linear elliptic PDEs with arbitrary mixed Neumann and Dirichlet boundary conditions. The key insight is that one can efficiently simulate reflecting Brownian motion (which models Neumann conditions) by replacing the balls used by WoS withstar-shapeddomains. We identify such domains via the closest point on the visibility silhouette, by simply augmenting a standard bounding volume hierarchy with normal information. Overall, WoSt is an easy modification of WoS, and retains the many attractive features of grid-free Monte Carlo methods such as progressive and view-dependent evaluation, trivial parallelization, and sublinear scaling to increasing geometric detail. 
    more » « less
  3. Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to typical eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols that enable the use of quantum computers to study this question. We present protocols both for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we present both a system-size-linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we formulate an exact state preparation approach based on matrix product states that yields quasipolynomial-depth circuits, as well as a variational approach with a polynomial-depth ansatz circuit. We also provide proof of principle state-preparation demonstrations on superconducting quantum hardware. 
    more » « less
  4. A bstract There are two parts to this work: first, we study the error correction properties of the real-space renormalization group (RG). The long-distance operators are the (approximately) correctable operators encoded in the physical algebra of short-distance operators. This is closely related to modeling the holographic map as a quantum error correction code. As opposed to holography, the real-space RG of a many-body quantum system does not have the complementary recovery property. We discuss the role of large N and a large gap in the spectrum of operators in the emergence of complementary recovery. Second, we study the operator algebra exact quantum error correction for any von Neumann algebra. We show that similar to the finite dimensional case, for any error map in between von Neumann algebras the Petz dual of the error map is a recovery map if the inclusion of the correctable subalgebra of operators has finite index. 
    more » « less
  5. This work contributes to nonlocal vector calculus as an indispensable mathematical tool for the study of nonlocal models that arises in a variety of applications. We define the nonlocal half-ball gradient, divergence and curl operators with general kernel functions (integrable or fractional type with finite or infinite supports) and study the associated nonlocal vector identities. We study the nonlocal function space on bounded domains associated with zero Dirichlet boundary conditions and the half-ball gradient operator and show it is a separable Hilbert space with smooth functions dense in it. A major result is the nonlocal Poincaré inequality, based on which a few applications are discussed, and these include applications to nonlocal convection–diffusion, nonlocal correspondence model of linear elasticity and nonlocal Helmholtz decomposition on bounded domains. 
    more » « less