Local Projector Research Articles

Embedding Quantum Many-Body Scars into Decoherence-Free Subspaces.

Quantum many-body scars are nonthermal excited eigenstates of nonintegrable Hamiltonians, which could support coherent revival dynamics from special initial states when scars form an equally spaced tower in the energy spectrum. For open quantum systems, engineering many-body scarred dynamics by a controlled coupling to the environment remains largely unexplored. Here, we provide a general framework to exactly embed quantum many-body scars into the decoherence-free subspaces of Lindblad master equations. The dissipative scarred dynamics manifest persistent periodic oscillations for generic initial states, and can be practically utilized to prepare scar states with potential quantum metrology applications. We construct the Liouvillian dissipators with the local projectors that annihilate the whole scar towers, and utilize the Hamiltonian part to rotate the undesired states out of the null space of dissipators. We demonstrate our protocol through several typical models hosting many-body scar towers and propose an experimental scheme to observe the dissipative scarred dynamics based on digital quantum simulations and resetting ancilla qubits.

Improving high-order VEM stability on badly-shaped elements

For the 2D and 3D Virtual Element Methods, a new approach to improve the conditioning of local and global matrices in the presence of badly-shaped polytopes is proposed. This new method defines the local projectors and the local degrees of freedom with respect to a set of scaled monomials recomputed on more well-shaped polytopes. This new approach is less computationally demanding than using the orthonormal polynomial basis. The effectiveness of our procedure is tested on different numerical examples characterized by challenging geometries of increasing complexity.

Fractionalization paves the way to local projector embeddings of quantum many-body scars

Fractionalization paves the way to local projector embeddings of quantum many-body scars

Analysis of divergence free conforming virtual elements for the Brinkman problem

In this paper, we develop stability analysis, including inverse inequality, [Formula: see text] norm equivalence and interpolation error estimates, for divergence free conforming virtual elements in arbitrary dimension. A local energy projector based on the local Stokes problem is suggested, which commutes with the divergence operator. After defining a discrete bilinear form and a stabilization involving only the boundary degrees of freedom (DoF) and parts of the interior DoF, a new divergence free conforming virtual element method is advanced for the Brinkman problem, which can be reduced to a simpler method due to the divergence free discrete velocity. An optimal convergence rate is derived for the discrete method. Furthermore, we achieve a uniform half convergence rate of the discrete method in consideration of the boundary layer phenomenon. Finally, some numerical results are provided to validate the convergence of the discrete method.

An unrefinement algorithm for planar THB-spline parameterizations

Unrefinement is a tool that allows to perform faster numerical simulations by controlling the level of precision in the specified area. We introduce an algorithm that creates a coarser geometry from an initial regular geometry, which is represented with respect to THB-splines, so that the coarser geometry approximates the initial geometry with a given level of the precision and is regular. The algorithm is based on the unrefinement of cells of the parameter domain, on which the geometry is defined, thus obtaining the new subdivision of the domain (the new subdomain hierarchy). The algorithm uses local linear projectors and optimization techniques in order to ensure the regularity and the approximation properties.

Accurate Local Homography-Based Projector Calibration Method

Accurate Local Homography-Based Projector Calibration Method

Compactly supported Wannier functions and strictly local projectors

Wannier functions that are maximally localized help in understanding many properties of crystalline materials. In the absence of topological obstructions, they are at least exponentially localized. In some cases such as flat-band Hamiltonians, it is possible to construct Wannier functions that are even more localized, so that they are compactly supported thus having zero support outside their corresponding locations. Under what general conditions is it possible to construct compactly supported Wannier functions? We answer this question in this paper. Specifically, we show that in 1d non-interacting tight-binding models, strict locality of the projection operator is a necessary and sufficient condition for a subspace to be spanned by a compactly supported orthogonal basis, independent of lattice translation symmetry. For any strictly local projector, we provide a procedure for obtaining such a basis. For higher dimensional systems, we discuss some additional conditions under which an occupied subspace is spanned by a compactly supported orthogonal basis, and show that the corresponding projectors are topologically trivial in many cases. We also show that a projector in arbitrary dimensions is strictly local if and only if for any chosen axis, its image is spanned by hybrid Wannier functions that are compactly supported along that axis.

Piecewise Divergence-Free Nonconforming Virtual Elements for Stokes Problem in Any Dimensions

Piecewise divergence-free nonconforming virtual elements are designed for Stokes problem in any dimensions. After introducing a local energy projector based on the Stokes problem and the stabilization, a divergence-free nonconforming virtual element method is proposed for Stokes problem. A detailed and rigorous error analysis is presented for the discrete method. An important property in the analysis is that the local energy projector commutes with the divergence operator. With the help of a divergence-free interpolation operator onto a generalized Raviart--Thomas element space, a pressure-robust nonconforming virtual element method is developed by simply modifying the right-hand side of the previous discretization. A reduced virtual element method is also discussed. Numerical results are provided to verify the theoretical convergence.

On D-invariant points and local Taylor interpolation on algebraic hypersurfaces in RN

We give the definition of D-invariant points on an irreducible algebraic hypersurface V in RN. We show that every regular point on irreducible quadratic hypersurface in RN is D-invariant. We prove that the local Taylor interpolation projector at a regular point a 2 V is an ideal projector if and only if a is D-invariant.

Local (T)HB-spline projectors via restricted hierarchical spline fitting

This paper is devoted to techniques for adaptive spline projection via quasi-interpolation, enabling the efficient approximation of given functions. We employ local least-squares fitting in restricted hierarchical spline spaces to establish novel projection operators for hierarchical splines of degree p. This leads to efficient spline projectors that require O(pd) floating point operations and O(1) evaluations of the given function per degree of freedom, while providing essentially the same accuracy as global approximation. Our spline projectors are based on a unifying framework for quasi-interpolation in hierarchical spline spaces. We present a detailed comparison with the scheme of Speleers and Manni (2016).

Bridging the hybrid high-order and virtual element methods

Abstract We present a unifying viewpoint on hybrid high-order and virtual element methods on general polytopal meshes in dimension $2$ or $3$, in terms of both formulation and analysis. We focus on a model Poisson problem. To build our bridge (i) we transcribe the (conforming) virtual element method into the hybrid high-order framework and (ii) we prove $H^m$ approximation properties for the local polynomial projector in terms of which the local virtual element discrete bilinear form is defined. This allows us to perform a unified analysis of virtual element/hybrid high-order methods, that differs from standard virtual element analyses by the fact that the approximation properties of the underlying virtual space are not explicitly used. As a complement to our unified analysis we also study interpolation in local virtual spaces, shedding light on the differences between the conforming and nonconforming cases.

Density functional theory plus dynamical mean-field theory with natural atomic orbital projectors

We introduce natural atomic orbitals as the local projector to define the correlated subspace for DFT + DMFT (density functional theory plus dynamical mean-field theory) calculation. The natural atomic orbitals are found to be stably constructed against the number and the radius of basis orbitals. It can also be self-consistently updated inside the DFT+DMFT loop. The spatial localization, electron occupation and the degree of correlation are investigated and compared with other conventional techniques. As a `natural' choice to describe the electron numbers, adopting natural atomic orbitals has advantage in terms of electron number counting. We further explore the reduction of computation cost by separating correlated orbitals into two subgroups based on the orbital occupancy. Our new recipe can serve as a useful choice for DFT+DMFT and related methods.

A mixed virtual element method for a pseudostress-based formulation of linear elasticity

In this paper we introduce and analyze a mixed virtual element method (mixed-VEM) for a pseudostress-displacement formulation of the linear elasticity problem with non-homogeneous Dirichlet boundary conditions. We follow a previous work by some of the authors, and employ a mixed formulation that does not require symmetric tensor spaces in the finite element discretization. More precisely, the main unknowns here are given by the pseudostress and the displacement, whereas other physical quantities such as the stress, the strain tensor of small deformations, and the rotation, are computed through simple postprocessing formulae in terms of the pseudostress variable. We first recall the corresponding variational formulation, and then summarize the main mixed-VEM ingredients that are required for our discrete analysis. In particular, we utilize a well-known local projector onto a suitable polynomial subspace to define a calculable version of our discrete bilinear form, whose continuous version requires information of the variables on the interior of each element. Next, we show that the global discrete bilinear form satisfies the hypotheses required by the Babuška–Brezzi theory. In this way, we conclude the well-posedness of our mixed-VEM scheme and derive the associated a priori error estimates for the virtual solutions as well as for the fully computable projections of them. Furthermore, we also introduce a second element-by-element postprocessing formula for the pseudostress, which yields an optimally convergent approximation of this unknown with respect to the broken H(div)-norm. In addition, this postprocessing formula can also be applied to the postprocessed stress tensor. Finally, several numerical results illustrating the good performance of the method and confirming the theoretical rates of convergence are presented.

Linear-Time Algorithm for Quantum 2SAT

A canonical result about satisfiability theory is that the 2-SAT problem can be solved in linear time, despite the NP-hardness of the 3-SAT problem. In the quantum 2-SAT problem, we are given a family of 2-qubit projectors $\Pi_{ij}$ on a system of $n$ qubits, and the task is to decide whether the Hamiltonian $H=\sum \Pi_{ij}$ has a 0-eigenvalue, or it is larger than $1/n^\alpha$ for some $\alpha=O(1)$. The problem is not only a natural extension of the classical 2-SAT problem to the quantum case, but is also equivalent to the problem of finding the ground state of 2-local frustration-free Hamiltonians of spin $\frac{1}{2}$, a well-studied model believed to capture certain key properties in modern condensed matter physics. While Bravyi has shown that the quantum 2-SAT problem has a classical polynomial-time algorithm, the running time of his algorithm is $O(n^4)$. In this paper we give a classical algorithm with linear running time in the number of local projectors, therefore achieving the best possible complexity.

A mixed virtual element method for the Brinkman problem

In this paper, we introduce and analyze a mixed virtual element method (mixed-VEM) for the two-dimensional Brinkman model of porous media flow with non-homogeneous Dirichlet boundary conditions. More precisely, we employ a dual-mixed formulation in which the only unknown is given by the pseudostress, whereas the velocity and pressure are computed via postprocessing formulae. We first recall the corresponding variational formulation, and then summarize the main mixed-VEM ingredients that are required for our discrete analysis. In particular, in order to define a calculable discrete bilinear form, whose continuous version involves deviatoric tensors, we propose two well-known alternatives for the local projector onto a suitable polynomial subspace, which allows the explicit integration of these terms. Next, we show that the global discrete bilinear form satisfies the hypotheses required by the Lax–Milgram lemma. In this way, we conclude the well-posedness of our mixed-VEM scheme and derive the associated a priori error estimates for the virtual solution as well as for the fully computable projection of it. Furthermore, we also introduce a second element-by-element postprocessing formula for the pseudostress, which yields an optimally convergent approximation of this unknown with respect to the broken ℍ(div)-norm. Finally, several numerical results illustrating the good performance of the method and confirming the theoretical rates of convergence are presented.

Simple proof of the detectability lemma and spectral gap amplification

The detectability lemma is a useful tool for probing the structure of gapped ground states of frustration-free Hamiltonians of lattice spin models. The lemma provides an estimate on the error incurred by approximating the ground space projector with a product of local projectors. We provide a new, simpler proof for the detectability lemma, which applies to an arbitrary ordering of the local projectors, and show that it is tight up to a constant factor. As an application we show how the lemma can be combined with a strong converse by Gao to obtain local spectral gap amplification: we show that by coarse-graining a local frustration-free Hamiltonian with a spectral gap $\gamma>0$ to a length scale $O(\gamma^{-1/2})$, one gets an Hamiltonian with an $\Omega(1)$ spectral gap.

Quantum 3-SAT Is QMA$_1$-Complete

Quantum satisfiability is a constraint satisfaction problem that generalizes classical boolean satisfiability. In the quantum $k$-SAT problem, each constraint is specified by a $k$-local projector and is satisfied by any state in its nullspace. Bravyi showed that quantum 2-SAT can be solved efficiently on a classical computer and that quantum $k$-SAT with $k\geq 4$ is QMA$_1$-complete [S. Bravyi, Efficient Algorithm for a Quantum Analogue of 2-SAT, eprint arXiv:quant-ph/0602108, 2006]. Quantum 3-SAT was known to be contained in QMA$_1$ [Bravyi, 2006], but its computational hardness was unknown until now. We prove that quantum 3-SAT is QMA$_1$-hard, and therefore complete for this complexity class.

A constructive quantum Lovasz local lemma for commuting projectors

The Quantum Satisfiability problem generalizes the Boolean satisfiability problem to the quantum setting by replacing classical clauses with local projectors. The Quantum Lov\'asz Local Lemma gives a sufficient condition for a Quantum Satisfiability problem to be satisfiable~\cite{ambainis2012quantum}, by generalizing the classical Lov\'asz Local Lemma. The next natural question that arises is: can a satisfying quantum state be \emph{efficiently} found, when these conditions hold? In this work we present such an algorithm, with the additional requirement that all the projectors commute. The proof follows the information theoretic proof given by Moser's breakthrough result in the classical setting~\cite{moser2009constructiveSTOC}. Similar results were independently published in~\cite{cubitt2011constructive, cubitt2013quantum}.

A constructive quantum Lovász local lemma for commuting projectors

The Quantum Satisfiability problem generalizes the Boolean satisfiability problem to the quantum setting by replacing classical clauses with local projectors. The Quantum Lovasz Local Lemma gives a...

Effortless quasi-interpolation in hierarchical spaces

We present a general and simple procedure to construct quasi-interpolants in hierarchical spaces. Such spaces are composed of a hierarchy of nested spaces and provide a flexible framework for local refinement. The proposed hierarchical quasi-interpolants are described in terms of the so-called truncated hierarchical basis. Assuming a quasi-interpolant is selected for each space associated with a particular level in the hierarchy, the hierarchical quasi-interpolants are obtained without any additional manipulation. The main properties (like polynomial reproduction) of the quasi-interpolants selected at each level are locally preserved in the hierarchical construction. We show how to construct hierarchical local projectors, and the local approximation order of the underling hierarchical space is also investigated. The presentation is detailed for the truncated hierarchical B-spline basis, and we discuss its extension to a more general framework.