-
Robustness of near-thermal dynamics on digital quantum computers
Authors:
Eli Chertkov,
Yi-Hsiang Chen,
Michael Lubasch,
David Hayes,
Michael Foss-Feig
Abstract:
Understanding the impact of gate errors on quantum circuits is crucial to determining the potential applications of quantum computers, especially in the absence of large-scale error-corrected hardware. We put forward analytical arguments, corroborated by extensive numerical and experimental evidence, that Trotterized quantum circuits simulating the time-evolution of systems near thermal equilibriu…
▽ More
Understanding the impact of gate errors on quantum circuits is crucial to determining the potential applications of quantum computers, especially in the absence of large-scale error-corrected hardware. We put forward analytical arguments, corroborated by extensive numerical and experimental evidence, that Trotterized quantum circuits simulating the time-evolution of systems near thermal equilibrium are substantially more robust to both quantum gate errors and Trotter (discretization) errors than is widely assumed. In Quantinuum's trapped-ion computers, the weakly entangling gates that appear in Trotterized circuits can be implemented natively, and their error rate is smaller when they generate less entanglement; from benchmarking, we know that the error for a gate $\exp[-i (Z\otimes Z) τ]$ decreases roughly linearly with $τ$, up to a small offset at $τ= 0$. We provide extensive evidence that this scaling, together with the robustness of near-thermal dynamics to both gate and discretization errors, facilitates substantial improvements in the achievable accuracy of Trotterized dynamics on near-term quantum computers. We make heavy use of a new theoretical tool -- a statistical ensemble of random product states that approximates a thermal state, which can be efficiently prepared with low noise on quantum computers. We outline how the random product state ensemble can be used to predict, optimize, and design Hamiltonian simulation experiments on near-thermal quantum systems.
△ Less
Submitted 14 October, 2024;
originally announced October 2024.
-
Nonlinear dynamics as a ground-state solution on quantum computers
Authors:
Albert J. Pool,
Alejandro D. Somoza,
Conor Mc Keever,
Michael Lubasch,
Birger Horstmann
Abstract:
For the solution of time-dependent nonlinear differential equations, we present variational quantum algorithms (VQAs) that encode both space and time in qubit registers. The spacetime encoding enables us to obtain the entire time evolution from a single ground-state computation. We describe a general procedure to construct efficient quantum circuits for the cost function evaluation required by VQA…
▽ More
For the solution of time-dependent nonlinear differential equations, we present variational quantum algorithms (VQAs) that encode both space and time in qubit registers. The spacetime encoding enables us to obtain the entire time evolution from a single ground-state computation. We describe a general procedure to construct efficient quantum circuits for the cost function evaluation required by VQAs. To mitigate the barren plateau problem during the optimization, we propose an adaptive multigrid strategy. The approach is illustrated for the nonlinear Burgers equation. We classically optimize quantum circuits to represent the desired ground-state solutions, run them on IBM Q System One and Quantinuum System Model H1, and demonstrate that current quantum computers are capable of accurately reproducing the exact results.
△ Less
Submitted 27 September, 2024; v1 submitted 25 March, 2024;
originally announced March 2024.
-
Noisy intermediate-scale quantum simulation of the one-dimensional wave equation
Authors:
Lewis Wright,
Conor Mc Keever,
Jeremy T. First,
Rory Johnston,
Jeremy Tillay,
Skylar Chaney,
Matthias Rosenkranz,
Michael Lubasch
Abstract:
We design and implement quantum circuits for the simulation of the one-dimensional wave equation on the Quantinuum H1-1 quantum computer. The circuit depth of our approach scales as $O(n^{2})$ for $n$ qubits representing the solution on $2^n$ grid points, and leads to infidelities of $O(2^{-4n} t^{2})$ for simulation time $t$ assuming smooth initial conditions. By varying the qubit count we study…
▽ More
We design and implement quantum circuits for the simulation of the one-dimensional wave equation on the Quantinuum H1-1 quantum computer. The circuit depth of our approach scales as $O(n^{2})$ for $n$ qubits representing the solution on $2^n$ grid points, and leads to infidelities of $O(2^{-4n} t^{2})$ for simulation time $t$ assuming smooth initial conditions. By varying the qubit count we study the interplay between the algorithmic and physical gate errors to identify the optimal working point of minimum total error. Our approach to simulating the wave equation can readily be adapted to other quantum processors and serve as an application-oriented benchmark.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
Towards adiabatic quantum computing using compressed quantum circuits
Authors:
Conor Mc Keever,
Michael Lubasch
Abstract:
We describe tensor network algorithms to optimize quantum circuits for adiabatic quantum computing. To suppress diabatic transitions, we include counterdiabatic driving in the optimization and utilize variational matrix product operators to represent adiabatic gauge potentials. Traditionally, Trotter product formulas are used to turn adiabatic time evolution into quantum circuits and the addition…
▽ More
We describe tensor network algorithms to optimize quantum circuits for adiabatic quantum computing. To suppress diabatic transitions, we include counterdiabatic driving in the optimization and utilize variational matrix product operators to represent adiabatic gauge potentials. Traditionally, Trotter product formulas are used to turn adiabatic time evolution into quantum circuits and the addition of counterdiabatic driving increases the circuit depth per time step. Instead, we classically optimize a parameterized quantum circuit of fixed depth to simultaneously capture adiabatic evolution together with counterdiabatic driving over many time steps. The methods are applied to the ground state preparation of quantum Ising chains with transverse and longitudinal fields. We show that the classically optimized circuits can significantly outperform Trotter product formulas. Additionally, we discuss how the approach can be used for combinatorial optimization.
△ Less
Submitted 20 June, 2024; v1 submitted 9 November, 2023;
originally announced November 2023.
-
A Modular Engine for Quantum Monte Carlo Integration
Authors:
Ismail Yunus Akhalwaya,
Adam Connolly,
Roland Guichard,
Steven Herbert,
Cahit Kargi,
Alexandre Krajenbrink,
Michael Lubasch,
Conor Mc Keever,
Julien Sorci,
Michael Spranger,
Ifan Williams
Abstract:
We present the Quantum Monte Carlo Integration (QMCI) engine developed by Quantinuum. It is a quantum computational tool for evaluating multi-dimensional integrals that arise in various fields of science and engineering such as finance. This white paper presents a detailed description of the architecture of the QMCI engine, including a variety of distribution-loading methods, a novel quantum ampli…
▽ More
We present the Quantum Monte Carlo Integration (QMCI) engine developed by Quantinuum. It is a quantum computational tool for evaluating multi-dimensional integrals that arise in various fields of science and engineering such as finance. This white paper presents a detailed description of the architecture of the QMCI engine, including a variety of distribution-loading methods, a novel quantum amplitude estimation method that improves the statistical robustness of QMCI calculations, and a library of statistical quantities that can be estimated. The QMCI engine is designed with modularity in mind, allowing for the continuous development of new quantum algorithms tailored in particular to financial applications. Additionally, the engine features a resource mode, which provides a precise resource quantification for the quantum circuits generated. The paper also includes extensive benchmarks that showcase the engine's performance, with a focus on the evaluation of various financial instruments.
△ Less
Submitted 11 August, 2023;
originally announced August 2023.
-
Realization of quantum signal processing on a noisy quantum computer
Authors:
Yuta Kikuchi,
Conor Mc Keever,
Luuk Coopmans,
Michael Lubasch,
Marcello Benedetti
Abstract:
Quantum signal processing (QSP) is a powerful toolbox for the design of quantum algorithms and can lead to asymptotically optimal computational costs. Its realization on noisy quantum computers without fault tolerance, however, is challenging because it requires a deep quantum circuit in general. We propose a strategy to run an entire QSP protocol on noisy quantum hardware by carefully reducing ov…
▽ More
Quantum signal processing (QSP) is a powerful toolbox for the design of quantum algorithms and can lead to asymptotically optimal computational costs. Its realization on noisy quantum computers without fault tolerance, however, is challenging because it requires a deep quantum circuit in general. We propose a strategy to run an entire QSP protocol on noisy quantum hardware by carefully reducing overhead costs at each step. To illustrate the approach, we consider the application of Hamiltonian simulation for which QSP implements a polynomial approximation of the time evolution operator. We test the protocol by running the algorithm on the Quantinuum H1-1 trapped-ion quantum computer powered by Honeywell. In particular, we compute the time dependence of bipartite entanglement entropies for Ising spin chains and find good agreements with exact numerical simulations. To make the best use of the device, we determine optimal experimental parameters by using a simplified error model for the hardware and numerically studying the trade-off between Hamiltonian simulation time, polynomial degree, and total accuracy. Our results are the first step in the experimental realization of QSP-based quantum algorithms.
△ Less
Submitted 27 September, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Barren plateaus in quantum tensor network optimization
Authors:
Enrique Cervero Martín,
Kirill Plekhanov,
Michael Lubasch
Abstract:
We analyze the barren plateau phenomenon in the variational optimization of quantum circuits inspired by matrix product states (qMPS), tree tensor networks (qTTN), and the multiscale entanglement renormalization ansatz (qMERA). We consider as the cost function the expectation value of a Hamiltonian that is a sum of local terms. For randomly chosen variational parameters we show that the variance o…
▽ More
We analyze the barren plateau phenomenon in the variational optimization of quantum circuits inspired by matrix product states (qMPS), tree tensor networks (qTTN), and the multiscale entanglement renormalization ansatz (qMERA). We consider as the cost function the expectation value of a Hamiltonian that is a sum of local terms. For randomly chosen variational parameters we show that the variance of the cost function gradient decreases exponentially with the distance of a Hamiltonian term from the canonical centre in the quantum tensor network. Therefore, as a function of qubit count, for qMPS most gradient variances decrease exponentially and for qTTN as well as qMERA they decrease polynomially. We also show that the calculation of these gradients is exponentially more efficient on a classical computer than on a quantum computer.
△ Less
Submitted 2 April, 2023; v1 submitted 1 September, 2022;
originally announced September 2022.
-
Classically optimized Hamiltonian simulation
Authors:
Conor Mc Keever,
Michael Lubasch
Abstract:
Hamiltonian simulation is a promising application for quantum computers to achieve a quantum advantage. We present classical algorithms based on tensor network methods to optimize quantum circuits for this task. We show that, compared to Trotter product formulas, the classically optimized circuits can be orders of magnitude more accurate and significantly extend the total simulation time.
Hamiltonian simulation is a promising application for quantum computers to achieve a quantum advantage. We present classical algorithms based on tensor network methods to optimize quantum circuits for this task. We show that, compared to Trotter product formulas, the classically optimized circuits can be orders of magnitude more accurate and significantly extend the total simulation time.
△ Less
Submitted 2 June, 2023; v1 submitted 23 May, 2022;
originally announced May 2022.
-
Variational quantum amplitude estimation
Authors:
Kirill Plekhanov,
Matthias Rosenkranz,
Mattia Fiorentini,
Michael Lubasch
Abstract:
We propose to perform amplitude estimation with the help of constant-depth quantum circuits that variationally approximate states during amplitude amplification. In the context of Monte Carlo (MC) integration, we numerically show that shallow circuits can accurately approximate many amplitude amplification steps. We combine the variational approach with maximum likelihood amplitude estimation [Y.…
▽ More
We propose to perform amplitude estimation with the help of constant-depth quantum circuits that variationally approximate states during amplitude amplification. In the context of Monte Carlo (MC) integration, we numerically show that shallow circuits can accurately approximate many amplitude amplification steps. We combine the variational approach with maximum likelihood amplitude estimation [Y. Suzuki et al., Quantum Inf. Process. 19, 75 (2020)] in variational quantum amplitude estimation (VQAE). VQAE typically has larger computational requirements than classical MC sampling. To reduce the variational cost, we propose adaptive VQAE and numerically show in 6 to 12 qubit simulations that it can outperform classical MC sampling.
△ Less
Submitted 28 February, 2022; v1 submitted 8 September, 2021;
originally announced September 2021.
-
Filtering variational quantum algorithms for combinatorial optimization
Authors:
David Amaro,
Carlo Modica,
Matthias Rosenkranz,
Mattia Fiorentini,
Marcello Benedetti,
Michael Lubasch
Abstract:
Current gate-based quantum computers have the potential to provide a computational advantage if algorithms use quantum hardware efficiently. To make combinatorial optimization more efficient, we introduce the Filtering Variational Quantum Eigensolver (F-VQE) which utilizes filtering operators to achieve faster and more reliable convergence to the optimal solution. Additionally we explore the use o…
▽ More
Current gate-based quantum computers have the potential to provide a computational advantage if algorithms use quantum hardware efficiently. To make combinatorial optimization more efficient, we introduce the Filtering Variational Quantum Eigensolver (F-VQE) which utilizes filtering operators to achieve faster and more reliable convergence to the optimal solution. Additionally we explore the use of causal cones to reduce the number of qubits required on a quantum computer. Using random weighted MaxCut problems, we numerically analyze our methods and show that they perform better than the original VQE algorithm and the Quantum Approximate Optimization Algorithm (QAOA). We also demonstrate the experimental feasibility of our algorithms on a Honeywell trapped-ion quantum processor.
△ Less
Submitted 10 February, 2022; v1 submitted 18 June, 2021;
originally announced June 2021.
-
A Quantum Inspired Approach to Exploit Turbulence Structures
Authors:
Nikita Gourianov,
Michael Lubasch,
Sergey Dolgov,
Quincy Y. van den Berg,
Hessam Babaee,
Peyman Givi,
Martin Kiffner,
Dieter Jaksch
Abstract:
Understanding turbulence is the key to our comprehension of many natural and technological flow processes. At the heart of this phenomenon lies its intricate multi-scale nature, describing the coupling between different-sized eddies in space and time. Here we introduce a new paradigm for analyzing the structure of turbulent flows by quantifying correlations between different length scales using me…
▽ More
Understanding turbulence is the key to our comprehension of many natural and technological flow processes. At the heart of this phenomenon lies its intricate multi-scale nature, describing the coupling between different-sized eddies in space and time. Here we introduce a new paradigm for analyzing the structure of turbulent flows by quantifying correlations between different length scales using methods inspired from quantum many-body physics. We present results for interscale correlations of two paradigmatic flow examples, and use these insights along with tensor network theory to design a structure-resolving algorithm for simulating turbulent flows. With this algorithm, we find that the incompressible Navier-Stokes equations can be accurately solved within a computational space reduced by over an order of magnitude compared to direct numerical simulation. Our quantum-inspired approach provides a pathway towards conducting computational fluid dynamics on quantum computers.
△ Less
Submitted 4 July, 2022; v1 submitted 10 June, 2021;
originally announced June 2021.
-
Variational inference with a quantum computer
Authors:
Marcello Benedetti,
Brian Coyle,
Mattia Fiorentini,
Michael Lubasch,
Matthias Rosenkranz
Abstract:
Inference is the task of drawing conclusions about unobserved variables given observations of related variables. Applications range from identifying diseases from symptoms to classifying economic regimes from price movements. Unfortunately, performing exact inference is intractable in general. One alternative is variational inference, where a candidate probability distribution is optimized to appr…
▽ More
Inference is the task of drawing conclusions about unobserved variables given observations of related variables. Applications range from identifying diseases from symptoms to classifying economic regimes from price movements. Unfortunately, performing exact inference is intractable in general. One alternative is variational inference, where a candidate probability distribution is optimized to approximate the posterior distribution over unobserved variables. For good approximations, a flexible and highly expressive candidate distribution is desirable. In this work, we use quantum Born machines as variational distributions over discrete variables. We apply the framework of operator variational inference to achieve this goal. In particular, we adopt two specific realizations: one with an adversarial objective and one based on the kernelized Stein discrepancy. We demonstrate the approach numerically using examples of Bayesian networks, and implement an experiment on an IBM quantum computer. Our techniques enable efficient variational inference with distributions beyond those that are efficiently representable on a classical computer.
△ Less
Submitted 4 November, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
Hardware-efficient variational quantum algorithms for time evolution
Authors:
Marcello Benedetti,
Mattia Fiorentini,
Michael Lubasch
Abstract:
Parameterized quantum circuits are a promising technology for achieving a quantum advantage. An important application is the variational simulation of time evolution of quantum systems. To make the most of quantum hardware, variational algorithms need to be as hardware-efficient as possible. Here we present alternatives to the time-dependent variational principle that are hardware-efficient and do…
▽ More
Parameterized quantum circuits are a promising technology for achieving a quantum advantage. An important application is the variational simulation of time evolution of quantum systems. To make the most of quantum hardware, variational algorithms need to be as hardware-efficient as possible. Here we present alternatives to the time-dependent variational principle that are hardware-efficient and do not require matrix inversion. In relation to imaginary time evolution, our approach significantly reduces the hardware requirements. With regards to real time evolution, where high precision can be important, we present algorithms of systematically increasing accuracy and hardware requirements. We numerically analyze the performance of our algorithms using quantum Hamiltonians with local interactions.
△ Less
Submitted 27 July, 2021; v1 submitted 25 September, 2020;
originally announced September 2020.
-
Parallel time-dependent variational principle algorithm for matrix product states
Authors:
Paul Secular,
Nikita Gourianov,
Michael Lubasch,
Sergey Dolgov,
Stephen R. Clark,
Dieter Jaksch
Abstract:
Combining the time-dependent variational principle (TDVP) algorithm with the parallelization scheme introduced by Stoudenmire and White for the density matrix renormalization group (DMRG), we present the first parallel matrix product state (MPS) algorithm capable of time evolving one-dimensional (1D) quantum lattice systems with long-range interactions. We benchmark the accuracy and performance of…
▽ More
Combining the time-dependent variational principle (TDVP) algorithm with the parallelization scheme introduced by Stoudenmire and White for the density matrix renormalization group (DMRG), we present the first parallel matrix product state (MPS) algorithm capable of time evolving one-dimensional (1D) quantum lattice systems with long-range interactions. We benchmark the accuracy and performance of the algorithm by simulating quenches in the long-range Ising and XY models. We show that our code scales well up to 32 processes, with parallel efficiencies as high as 86%. Finally, we calculate the dynamical correlation function of a 201-site Heisenberg XXX spin chain with $1/r^2$ interactions, which is challenging to compute sequentially. These results pave the way for the application of tensor networks to increasingly complex many-body systems.
△ Less
Submitted 12 June, 2020; v1 submitted 12 December, 2019;
originally announced December 2019.
-
Variational quantum algorithms for nonlinear problems
Authors:
Michael Lubasch,
Jaewoo Joo,
Pierre Moinier,
Martin Kiffner,
Dieter Jaksch
Abstract:
We show that nonlinear problems including nonlinear partial differential equations can be efficiently solved by variational quantum computing. We achieve this by utilizing multiple copies of variational quantum states to treat nonlinearities efficiently and by introducing tensor networks as a programming paradigm. The key concepts of the algorithm are demonstrated for the nonlinear Schrödinger equ…
▽ More
We show that nonlinear problems including nonlinear partial differential equations can be efficiently solved by variational quantum computing. We achieve this by utilizing multiple copies of variational quantum states to treat nonlinearities efficiently and by introducing tensor networks as a programming paradigm. The key concepts of the algorithm are demonstrated for the nonlinear Schrödinger equation as a canonical example. We numerically show that the variational quantum ansatz can be exponentially more efficient than matrix product states and present experimental proof-of-principle results obtained on an IBM Q device.
△ Less
Submitted 17 December, 2019; v1 submitted 21 July, 2019;
originally announced July 2019.
-
Bosonic Fractional Quantum Hall States on a Finite Cylinder
Authors:
Paolo Rosson,
Michael Lubasch,
Martin Kiffner,
Dieter Jaksch
Abstract:
We investigate the ground state properties of a bosonic Harper-Hofstadter model with local interactions on a finite cylindrical lattice with filling fraction $ν=1/2$. We find that our system supports topologically ordered states by calculating the topological entanglement entropy, and its value is in good agreement with the theoretical value for the $ν=1/2$ Laughlin state. By exploring the behavio…
▽ More
We investigate the ground state properties of a bosonic Harper-Hofstadter model with local interactions on a finite cylindrical lattice with filling fraction $ν=1/2$. We find that our system supports topologically ordered states by calculating the topological entanglement entropy, and its value is in good agreement with the theoretical value for the $ν=1/2$ Laughlin state. By exploring the behaviour of the density profiles, edge currents and single-particle correlation functions, we find that the ground state on the cylinder shows all signatures of a fractional quantum Hall state even for large values of the magnetic flux density. Furthermore, we determine the dependence of the correlation functions and edge currents on the interaction strength. We find that depending on the magnetic flux density, the transition towards Laughlin-like behaviour can be either smooth or happens abruptly for some critical interaction strength.
△ Less
Submitted 11 March, 2019; v1 submitted 11 January, 2019;
originally announced January 2019.
-
Multigrid Renormalization
Authors:
Michael Lubasch,
Pierre Moinier,
Dieter Jaksch
Abstract:
We combine the multigrid (MG) method with state-of-the-art concepts from the variational formulation of the numerical renormalization group. The resulting MG renormalization (MGR) method is a natural generalization of the MG method for solving partial differential equations. When the solution on a grid of $N$ points is sought, our MGR method has a computational cost scaling as…
▽ More
We combine the multigrid (MG) method with state-of-the-art concepts from the variational formulation of the numerical renormalization group. The resulting MG renormalization (MGR) method is a natural generalization of the MG method for solving partial differential equations. When the solution on a grid of $N$ points is sought, our MGR method has a computational cost scaling as $\mathcal{O}(\log(N))$, as opposed to $\mathcal{O}(N)$ for the best standard MG method. Therefore MGR can exponentially speed up standard MG computations. To illustrate our method, we develop a novel algorithm for the ground state computation of the nonlinear Schrödinger equation. Our algorithm acts variationally on tensor products and updates the tensors one after another by solving a local nonlinear optimization problem. We compare several different methods for the nonlinear tensor update and find that the Newton method is the most efficient as well as precise. The combination of MGR with our nonlinear ground state algorithm produces accurate results for the nonlinear Schrödinger equation on $N = 10^{18}$ grid points in three spatial dimensions.
△ Less
Submitted 15 July, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Tensor network states in time-bin quantum optics
Authors:
Michael Lubasch,
Antonio A. Valido,
Jelmer J. Renema,
W. Steven Kolthammer,
Dieter Jaksch,
Myungshik S. Kim,
Ian Walmsley,
Raúl García-Patrón
Abstract:
The current shift in the quantum optics community towards large-size experiments -- with many modes and photons -- necessitates new classical simulation techniques that go beyond the usual phase space formulation of quantum mechanics. To address this pressing demand we formulate linear quantum optics in the language of tensor network states. As a toy model, we extensively analyze the quantum and c…
▽ More
The current shift in the quantum optics community towards large-size experiments -- with many modes and photons -- necessitates new classical simulation techniques that go beyond the usual phase space formulation of quantum mechanics. To address this pressing demand we formulate linear quantum optics in the language of tensor network states. As a toy model, we extensively analyze the quantum and classical correlations of time-bin interference in a single fiber loop. We then generalize our results to more complex time-bin quantum setups and identify different classes of architectures for high-complexity and low-overhead boson sampling experiments.
△ Less
Submitted 12 June, 2018; v1 submitted 28 December, 2017;
originally announced December 2017.
-
Systematic construction of density functionals based on matrix product state computations
Authors:
Michael Lubasch,
Johanna I. Fuks,
Heiko Appel,
Angel Rubio,
J. Ignacio Cirac,
Mari-Carmen Bañuls
Abstract:
We propose a systematic procedure for the approximation of density functionals in density functional theory that consists of two parts. First, for the efficient approximation of a general density functional, we introduce an efficient ansatz whose non-locality can be increased systematically. Second, we present a fitting strategy that is based on systematically increasing a reasonably chosen set of…
▽ More
We propose a systematic procedure for the approximation of density functionals in density functional theory that consists of two parts. First, for the efficient approximation of a general density functional, we introduce an efficient ansatz whose non-locality can be increased systematically. Second, we present a fitting strategy that is based on systematically increasing a reasonably chosen set of training densities. We investigate our procedure in the context of strongly correlated fermions on a one-dimensional lattice in which we compute accurate training densities with the help of matrix product states. Focusing on the exchange-correlation energy, we demonstrate how an efficient approximation can be found that includes and systematically improves beyond the local density approximation. Importantly, this systematic improvement is shown for target densities that are quite different from the training densities.
△ Less
Submitted 25 August, 2016; v1 submitted 21 March, 2016;
originally announced March 2016.
-
Algorithms for finite Projected Entangled Pair States
Authors:
Michael Lubasch,
J. Ignacio Cirac,
Mari-Carmen Bañuls
Abstract:
Projected Entangled Pair States (PEPS) are a promising ansatz for the study of strongly correlated quantum many-body systems in two dimensions. But due to their high computational cost, developing and improving PEPS algorithms is necessary to make the ansatz widely usable in practice. Here we analyze several algorithmic aspects of the method. On the one hand, we quantify the connection between the…
▽ More
Projected Entangled Pair States (PEPS) are a promising ansatz for the study of strongly correlated quantum many-body systems in two dimensions. But due to their high computational cost, developing and improving PEPS algorithms is necessary to make the ansatz widely usable in practice. Here we analyze several algorithmic aspects of the method. On the one hand, we quantify the connection between the correlation length of the PEPS and the accuracy of its approximate contraction, and discuss how purifications can be used in the latter. On the other, we present algorithmic improvements for the update of the tensor that introduce drastic gains in the numerical conditioning and the efficiency of the algorithms. Finally, the state-of-the-art general PEPS code is benchmarked with the Heisenberg and quantum Ising models on lattices of up to $21 \times 21$ sites.
△ Less
Submitted 4 September, 2014; v1 submitted 13 May, 2014;
originally announced May 2014.
-
Unifying Projected Entangled Pair States contractions
Authors:
Michael Lubasch,
J. Ignacio Cirac,
Mari-Carmen Bañuls
Abstract:
The approximate contraction of a Projected Entangled Pair States (PEPS) tensor network is a fundamental ingredient of any PEPS algorithm, required for the optimization of the tensors in ground state search or time evolution, as well as for the evaluation of expectation values. An exact contraction is in general impossible, and the choice of the approximating procedure determines the efficiency and…
▽ More
The approximate contraction of a Projected Entangled Pair States (PEPS) tensor network is a fundamental ingredient of any PEPS algorithm, required for the optimization of the tensors in ground state search or time evolution, as well as for the evaluation of expectation values. An exact contraction is in general impossible, and the choice of the approximating procedure determines the efficiency and accuracy of the algorithm. We analyze different previous proposals for this approximation, and show that they can be understood via the form of their environment, i.e. the operator that results from contracting part of the network. This provides physical insight into the limitation of various approaches, and allows us to introduce a new strategy, based on the idea of clusters, that unifies previous methods. The resulting contraction algorithm interpolates naturally between the cheapest and most imprecise and the most costly and most precise method. We benchmark the different algorithms with finite PEPS, and show how the cluster strategy can be used for both the tensor optimization and the calculation of expectation values. Additionally, we discuss its applicability to the parallelization of PEPS and to infinite systems (iPEPS).
△ Less
Submitted 7 April, 2014; v1 submitted 26 November, 2013;
originally announced November 2013.
-
Generalized Householder Transformations for the Complex Symmetric Eigenvalue Problem
Authors:
J. H. Noble,
M. Lubasch,
U. D. Jentschura
Abstract:
We present an intuitive and scalable algorithm for the diagonalization of complex symmetric matrices, which arise from the projection of pseudo--Hermitian and complex scaled Hamiltonians onto a suitable basis set of "trial" states. The algorithm diagonalizes complex and symmetric (non--Hermitian) matrices and is easily implemented in modern computer languages. It is based on generalized Householde…
▽ More
We present an intuitive and scalable algorithm for the diagonalization of complex symmetric matrices, which arise from the projection of pseudo--Hermitian and complex scaled Hamiltonians onto a suitable basis set of "trial" states. The algorithm diagonalizes complex and symmetric (non--Hermitian) matrices and is easily implemented in modern computer languages. It is based on generalized Householder transformations and relies on iterative similarity transformations T -> T' = Q^T T Q, where Q is a complex and orthogonal, but not unitary, matrix, i.e, Q^T equals Q^(-1) but Q^+ is different from Q^(-1). We present numerical reference data to support the scalability of the algorithm. We construct the generalized Householder transformations from the notion that the conserved scalar product of eigenstates Psi_n and Psi_m of a pseudo-Hermitian quantum mechanical Hamiltonian can be reformulated in terms of the generalized indefinite inner product [integral of the product Psi_n(x,t) Psi_m(x,t) over dx], where the integrand is locally defined, and complex conjugation is avoided. A few example calculations are described which illustrate the physical origin of the ideas used in the construction of the algorithm.
△ Less
Submitted 9 September, 2013; v1 submitted 24 January, 2013;
originally announced January 2013.
-
Adiabatic Preparation of a Heisenberg Antiferromagnet Using an Optical Superlattice
Authors:
Michael Lubasch,
Valentin Murg,
Ulrich Schneider,
J. Ignacio Cirac,
Mari-Carmen Bañuls
Abstract:
We analyze the possibility to prepare a Heisenberg antiferromagnet with cold fermions in optical lattices, starting from a band insulator and adiabatically changing the lattice potential. The numerical simulation of the dynamics in 1D allows us to identify the conditions for success, and to study the influence that the presence of holes in the initial state may have on the protocol. We also extend…
▽ More
We analyze the possibility to prepare a Heisenberg antiferromagnet with cold fermions in optical lattices, starting from a band insulator and adiabatically changing the lattice potential. The numerical simulation of the dynamics in 1D allows us to identify the conditions for success, and to study the influence that the presence of holes in the initial state may have on the protocol. We also extend our results to two-dimensional systems.
△ Less
Submitted 6 November, 2011; v1 submitted 8 June, 2011;
originally announced June 2011.
-
Dynamical enhancement of spatial entanglement in massive particles
Authors:
Michael Lubasch,
Florian Mintert,
Sandro Wimberger
Abstract:
We discuss dynamical enhancement of entanglement in a driven Bose-Hubbard model and find an enhancement of two orders of magnitude which is robust against fluctuations in experimental parameters.
We discuss dynamical enhancement of entanglement in a driven Bose-Hubbard model and find an enhancement of two orders of magnitude which is robust against fluctuations in experimental parameters.
△ Less
Submitted 21 September, 2010;
originally announced September 2010.
-
Detection of avoided crossings by fidelity
Authors:
Patrick Plötz,
Michael Lubasch,
Sandro Wimberger
Abstract:
The fidelity, defined as overlap of eigenstates of two slightly different Hamiltonians, is proposed as an efficient detector of avoided c rossings in the energy spectrum. This new application of fidelity is motivated for model systems, and its value for analyzing complex qua ntum spectra is underlined by applying it to a random matrix model and a tilted Bose-Hubbard system.
The fidelity, defined as overlap of eigenstates of two slightly different Hamiltonians, is proposed as an efficient detector of avoided c rossings in the energy spectrum. This new application of fidelity is motivated for model systems, and its value for analyzing complex qua ntum spectra is underlined by applying it to a random matrix model and a tilted Bose-Hubbard system.
△ Less
Submitted 2 February, 2011; v1 submitted 23 September, 2009;
originally announced September 2009.
-
Structure, Time Propagation and Dissipative Terms for Resonances
Authors:
U. D. Jentschura,
A. Surzhykov,
M. Lubasch,
J. Zinn-Justin
Abstract:
For odd anharmonic oscillators, it is well known that complex scaling can be used to determine resonance energy eigenvalues and the corresponding eigenvectors in complex rotated space. We briefly review and discuss various methods for the numerical determination of such eigenvalues, and also discuss the connection to the case of purely imaginary coupling, which is PT-symmetric. Moreover, we show…
▽ More
For odd anharmonic oscillators, it is well known that complex scaling can be used to determine resonance energy eigenvalues and the corresponding eigenvectors in complex rotated space. We briefly review and discuss various methods for the numerical determination of such eigenvalues, and also discuss the connection to the case of purely imaginary coupling, which is PT-symmetric. Moreover, we show that a suitable generalization of the complex scaling method leads to an algorithm for the time propagation of wave packets in potentials which give rise to unstable resonances. This leads to a certain unification of the structure and the dynamics. Our time propagation results agree with known quantum dynamics solvers and allow for a natural incorporation of structural perturbations (e.g., due to dissipative processes) into the quantum dynamics.
△ Less
Submitted 15 February, 2008; v1 submitted 7 November, 2007;
originally announced November 2007.
-
Quantum Dot Potentials: Symanzik Scaling, Resurgent Expansions and Quantum Dynamics
Authors:
A. Surzhykov,
M. Lubasch,
J. Zinn-Justin,
U. D. Jentschura
Abstract:
This article is concerned with a special class of the ``double-well-like'' potentials that occur naturally in the analysis of finite quantum systems. Special attention is paid, in particular, to the so-called Fokker-Planck potential, which has a particular property: the perturbation series for the ground-state energy vanishes to all orders in the coupling parameter, but the actual ground-state e…
▽ More
This article is concerned with a special class of the ``double-well-like'' potentials that occur naturally in the analysis of finite quantum systems. Special attention is paid, in particular, to the so-called Fokker-Planck potential, which has a particular property: the perturbation series for the ground-state energy vanishes to all orders in the coupling parameter, but the actual ground-state energy is positive and dominated by instanton configurations of the form exp(-a/g), where a is the instanton action. The instanton effects are most naturally taken into account within the modified Bohr-Sommerfeld quantization conditions whose expansion leads to the generalized perturbative expansions (so-called resurgent expansions) for the energy values of the Fokker-Planck potential. Until now, these resurgent expansions have been mainly applied for small values of coupling parameter g, while much less attention has been paid to the strong-coupling regime. In this contribution, we compare the energy values, obtained by directly resumming generalized Bohr-Sommerfeld quantization conditions, to the strong-coupling expansion, for which we determine the first few expansion coefficients in powers of g^(-2/3). Detailed calculations are performed for a wide range of coupling parameters g and indicate a considerable overlap between the regions of validity of the weak-coupling resurgent series and of the strong-coupling expansion. Apart from the analysis of the energy spectrum of the Fokker-Planck Hamiltonian, we also briefly discuss the computation of its eigenfunctions. These eigenfunctions may be utilized for the numerical integration of the (single-particle) time-dependent Schroedinger equation and, hence, for studying the dynamical evolution of the wavepackets in the double-well-like potentials.
△ Less
Submitted 27 November, 2006; v1 submitted 1 September, 2006;
originally announced September 2006.