We present a simplified analog quantum simulation protocol for preparing quantum states that embe... more We present a simplified analog quantum simulation protocol for preparing quantum states that embed solutions of parabolic partial differential equations, including the heat, Black-Scholes and Fokker-Planck equations. The key idea is to approximate the heat equations by a system of hyperbolic heat equations that involve only first-order differential operators. This scheme requires relatively simple interaction terms in the Hamiltonian, which are the electric and magnetic dipole moment-like interaction terms that would be present in a Jaynes-Cummings-like model. For a ddimensional problem, we show that it is much more appropriate to use a single d-level quantum system-a qudit-instead of its qubit counterpart, and d + 1 qumodes. The total resource cost is efficient in d and precision error, and has potential for realisability for instance in cavity and circuit QED systems.
Matrix geometric means between two positive definite matrices can be defined equivalently from di... more Matrix geometric means between two positive definite matrices can be defined equivalently from distinct perspectives -- as solutions to certain nonlinear systems of equations, as points along geodesics in Riemannian geometry, and as solutions to certain optimisation problems. This diversity already suggests the potential for varied applications, as well as acting as a bridge between different domains. Here we devise new quantum subroutines to efficiently prepare quantum unitary operators that embed the standard matrix geometric mean and its generalisations called the weighted matrix geometric mean. This enables the construction of solutions to the algebraic Riccati equation, which is an important class of nonlinear systems of equations that appears in machine learning, optimal control, estimation, and filtering. Using these subroutines, we present a new class of quantum learning algorithms called quantum geometric mean metric learning. This has applications in efficiently finding the best distance measure and solving classification problems in the weakly supervised limit and for anomaly detection, for both classical and quantum problems. We also show how our method can be generalised to a particular p^th-order system of nonlinear equations. These quantum subroutines for matrix geometric means are also useful in other areas of quantum information. For example, we show how to use them in the estimation of geometric Renyi relative entropies and the Uhlmann fidelity by means of the Fuchs--Caves observable. In particular, our quantum algorithms for estimating the Uhlmann and Matsumoto fidelities have optimal dependence on the precision. Finally, we provide a BQP-complete problem based on matrix geometric means that can be solved by our subroutines, thus characterising their computational capability.
Quantum simulators were originally proposed for simulating one partial differential equation (PDE... more Quantum simulators were originally proposed for simulating one partial differential equation (PDE) in particular-Schrödinger's equation. Can quantum simulators also efficiently simulate other PDEs? While most computational methods for PDEs-both classical and quantum-are digital (PDEs must be discretised first), PDEs have continuous degrees of freedom. This suggests that an analog representation can be more natural. While digital quantum degrees of freedom are usually described by qubits, the analog or continuous quantum degrees of freedom can be captured by qumodes. Based on a method called Schrödingerisation, we show how to directly map D-dimensional linear PDEs onto a (D + 1)-qumode quantum system where analog or continuous-variable Hamiltonian simulation on D + 1 qumodes can be used. This very simple methodology does not require one to discretise PDEs first, and it is not only applicable to linear PDEs but also to some nonlinear PDEs and systems of nonlinear ODEs. We show some examples using this method, including the Liouville equation, heat equation, Fokker-Planck equation, Black-Scholes equations, wave equation and Maxwell's equations. We also devise new protocols for linear PDEs with random coefficients, important in uncertainty quantification, where it is clear how the analog or continuous-variable framework is most natural. This also raises the possibility that some PDEs may be simulated directly on analog quantum systems by using Hamiltonians natural for those quantum systems.
Quantum simulation is known to be capable of simulating certain dynamical systems in continuous t... more Quantum simulation is known to be capable of simulating certain dynamical systems in continuous time-Schrödinger's equations being the most direct and well-known-more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger's equations via a method called Schrödingerisation [1]. Building on the observation that iterative methods in linear algebra, and more generally discrete linear dynamical systems, can be viewed as discrete time approximations of dynamical systems which evolve continuously in time, we can apply the Schrödingerisation technique. Thus quantum simulation can be directly applied to the continuous-time limits of some of the simplest iterative methods. This applies to general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we introduce the quantum Jacobi and quantum power methods for solving the quantum linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix respectively. The proposed quantum simulation can be performed on either discrete-variable quantum systems or on hybrid continuousvariable and discrete-variable quantum systems. This framework provides an interesting alternative method to solve linear algebra problems using quantum simulation.
Quantum simulation is known to be capable of simulating certain dynamical systems in continuous t... more Quantum simulation is known to be capable of simulating certain dynamical systems in continuous time-Schrödinger's equations being the most direct and well-known-more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger's equations via a method called Schrödingerisation [1]. Building on the observation that iterative methods in linear algebra, and more generally discrete linear dynamical systems, can be viewed as discrete time approximations of dynamical systems which evolve continuously in time, we can apply the Schrödingerisation technique. Thus quantum simulation can be directly applied to the continuous-time limits of some of the simplest iterative methods. This applies to general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we introduce the quantum Jacobi and quantum power methods for solving the quantum linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix respectively. The proposed quantum simulation can be performed on either discrete-variable quantum systems or on hybrid continuousvariable and discrete-variable quantum systems. This framework provides an interesting alternative method to solve linear algebra problems using quantum simulation.
This paper explores the feasibility of quantum simulation for partial differential equations (PDE... more This paper explores the feasibility of quantum simulation for partial differential equations (PDEs) with physical boundary or interface conditions. Semi-discretisation of such problems does not necessarily yield Hamiltonian dynamics and even alters the Hamiltonian structure of the dynamics when boundary and interface conditions are included. This seemingly intractable issue can be resolved by using a recently introduced Schrödingerisation method [JLY22a,JLY22b]-it converts any linear PDEs and ODEs with non-Hermitian dynamics to a system of Schrödinger equations, via the so-called warped phase transformation that maps the equation into one higher dimension. We implement this method for several typical problems, including the linear convection equation with inflow boundary conditions and the heat equation with Dirichlet and Neumann boundary conditions. For interface problems we study the (parabolic) Stefan problem, linear convection, and linear Liouville equations with discontinuous and even measure-valued coefficients. We perform numerical experiments to demonstrate the validity of this approach, which helps to bridge the gap between available quantum algorithms and computational models for classical and quantum dynamics with boundary and interface conditions.
We construct quantum algorithms to compute the solution and/or physical observables of nonlinear ... more We construct quantum algorithms to compute the solution and/or physical observables of nonlinear ordinary differential equations (ODEs) and nonlinear Hamilton-Jacobi equations (HJE) via linear representations or exact mappings between nonlinear ODEs/HJE and linear partial differential equations (the Liouville equation and the Koopman-von Neumann equation). The connection between the linear representations and the original nonlinear system is established through the Dirac delta function or the level set mechanism. We compare the quantum linear systems algorithms based methods and the quantum simulation methods arising from different numerical approximations, including the finite difference discretisations and the Fourier spectral discretisations for the two different linear representations, with the result showing that the quantum simulation methods usually give the best performance in time complexity. We also propose the Schrödinger framework to solve the Liouville equation for the HJE, since it can be recast as the semiclassical limit of the Wigner transform of the Schrödinger equation. Comparsion between the Schrödinger and the Liouville framework will also be made.
Solving the time-dependent Schrödinger equation is an important application area for quantum algo... more Solving the time-dependent Schrödinger equation is an important application area for quantum algorithms. We consider Schrödinger's equation in the semi-classical regime. Here the solutions exhibit strong multiple-scale behavior due to a small parameter ℏ, in the sense that the dynamics of the quantum states and the induced observables can occur on different spatial and temporal scales. Such a Schrödinger equation finds many applications, including in Born-Oppenheimer molecular dynamics and Ehrenfest dynamics. This paper considers quantum analogues of pseudo-spectral (PS) methods on classical computers. Estimates on the gate counts in terms of ℏ and the precision ε are obtained. It is found that the number of required qubits, m, scales only logarithmically with respect to ℏ. When the solution has bounded derivatives up to order ℓ, the symmetric Trotting method has gate complexity O((εℏ)−12polylog(ε−32ℓℏ−1−12ℓ)), provided that the diagonal unitary operators in the pseudo-spectral ...
We present a simplified analog quantum simulation protocol for preparing quantum states that embe... more We present a simplified analog quantum simulation protocol for preparing quantum states that embed solutions of parabolic partial differential equations, including the heat, Black-Scholes and Fokker-Planck equations. The key idea is to approximate the heat equations by a system of hyperbolic heat equations that involve only first-order differential operators. This scheme requires relatively simple interaction terms in the Hamiltonian, which are the electric and magnetic dipole moment-like interaction terms that would be present in a Jaynes-Cummings-like model. For a ddimensional problem, we show that it is much more appropriate to use a single d-level quantum system-a qudit-instead of its qubit counterpart, and d + 1 qumodes. The total resource cost is efficient in d and precision error, and has potential for realisability for instance in cavity and circuit QED systems.
Matrix geometric means between two positive definite matrices can be defined equivalently from di... more Matrix geometric means between two positive definite matrices can be defined equivalently from distinct perspectives -- as solutions to certain nonlinear systems of equations, as points along geodesics in Riemannian geometry, and as solutions to certain optimisation problems. This diversity already suggests the potential for varied applications, as well as acting as a bridge between different domains. Here we devise new quantum subroutines to efficiently prepare quantum unitary operators that embed the standard matrix geometric mean and its generalisations called the weighted matrix geometric mean. This enables the construction of solutions to the algebraic Riccati equation, which is an important class of nonlinear systems of equations that appears in machine learning, optimal control, estimation, and filtering. Using these subroutines, we present a new class of quantum learning algorithms called quantum geometric mean metric learning. This has applications in efficiently finding the best distance measure and solving classification problems in the weakly supervised limit and for anomaly detection, for both classical and quantum problems. We also show how our method can be generalised to a particular p^th-order system of nonlinear equations. These quantum subroutines for matrix geometric means are also useful in other areas of quantum information. For example, we show how to use them in the estimation of geometric Renyi relative entropies and the Uhlmann fidelity by means of the Fuchs--Caves observable. In particular, our quantum algorithms for estimating the Uhlmann and Matsumoto fidelities have optimal dependence on the precision. Finally, we provide a BQP-complete problem based on matrix geometric means that can be solved by our subroutines, thus characterising their computational capability.
Quantum simulators were originally proposed for simulating one partial differential equation (PDE... more Quantum simulators were originally proposed for simulating one partial differential equation (PDE) in particular-Schrödinger's equation. Can quantum simulators also efficiently simulate other PDEs? While most computational methods for PDEs-both classical and quantum-are digital (PDEs must be discretised first), PDEs have continuous degrees of freedom. This suggests that an analog representation can be more natural. While digital quantum degrees of freedom are usually described by qubits, the analog or continuous quantum degrees of freedom can be captured by qumodes. Based on a method called Schrödingerisation, we show how to directly map D-dimensional linear PDEs onto a (D + 1)-qumode quantum system where analog or continuous-variable Hamiltonian simulation on D + 1 qumodes can be used. This very simple methodology does not require one to discretise PDEs first, and it is not only applicable to linear PDEs but also to some nonlinear PDEs and systems of nonlinear ODEs. We show some examples using this method, including the Liouville equation, heat equation, Fokker-Planck equation, Black-Scholes equations, wave equation and Maxwell's equations. We also devise new protocols for linear PDEs with random coefficients, important in uncertainty quantification, where it is clear how the analog or continuous-variable framework is most natural. This also raises the possibility that some PDEs may be simulated directly on analog quantum systems by using Hamiltonians natural for those quantum systems.
Quantum simulation is known to be capable of simulating certain dynamical systems in continuous t... more Quantum simulation is known to be capable of simulating certain dynamical systems in continuous time-Schrödinger's equations being the most direct and well-known-more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger's equations via a method called Schrödingerisation [1]. Building on the observation that iterative methods in linear algebra, and more generally discrete linear dynamical systems, can be viewed as discrete time approximations of dynamical systems which evolve continuously in time, we can apply the Schrödingerisation technique. Thus quantum simulation can be directly applied to the continuous-time limits of some of the simplest iterative methods. This applies to general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we introduce the quantum Jacobi and quantum power methods for solving the quantum linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix respectively. The proposed quantum simulation can be performed on either discrete-variable quantum systems or on hybrid continuousvariable and discrete-variable quantum systems. This framework provides an interesting alternative method to solve linear algebra problems using quantum simulation.
Quantum simulation is known to be capable of simulating certain dynamical systems in continuous t... more Quantum simulation is known to be capable of simulating certain dynamical systems in continuous time-Schrödinger's equations being the most direct and well-known-more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger's equations via a method called Schrödingerisation [1]. Building on the observation that iterative methods in linear algebra, and more generally discrete linear dynamical systems, can be viewed as discrete time approximations of dynamical systems which evolve continuously in time, we can apply the Schrödingerisation technique. Thus quantum simulation can be directly applied to the continuous-time limits of some of the simplest iterative methods. This applies to general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we introduce the quantum Jacobi and quantum power methods for solving the quantum linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix respectively. The proposed quantum simulation can be performed on either discrete-variable quantum systems or on hybrid continuousvariable and discrete-variable quantum systems. This framework provides an interesting alternative method to solve linear algebra problems using quantum simulation.
This paper explores the feasibility of quantum simulation for partial differential equations (PDE... more This paper explores the feasibility of quantum simulation for partial differential equations (PDEs) with physical boundary or interface conditions. Semi-discretisation of such problems does not necessarily yield Hamiltonian dynamics and even alters the Hamiltonian structure of the dynamics when boundary and interface conditions are included. This seemingly intractable issue can be resolved by using a recently introduced Schrödingerisation method [JLY22a,JLY22b]-it converts any linear PDEs and ODEs with non-Hermitian dynamics to a system of Schrödinger equations, via the so-called warped phase transformation that maps the equation into one higher dimension. We implement this method for several typical problems, including the linear convection equation with inflow boundary conditions and the heat equation with Dirichlet and Neumann boundary conditions. For interface problems we study the (parabolic) Stefan problem, linear convection, and linear Liouville equations with discontinuous and even measure-valued coefficients. We perform numerical experiments to demonstrate the validity of this approach, which helps to bridge the gap between available quantum algorithms and computational models for classical and quantum dynamics with boundary and interface conditions.
We construct quantum algorithms to compute the solution and/or physical observables of nonlinear ... more We construct quantum algorithms to compute the solution and/or physical observables of nonlinear ordinary differential equations (ODEs) and nonlinear Hamilton-Jacobi equations (HJE) via linear representations or exact mappings between nonlinear ODEs/HJE and linear partial differential equations (the Liouville equation and the Koopman-von Neumann equation). The connection between the linear representations and the original nonlinear system is established through the Dirac delta function or the level set mechanism. We compare the quantum linear systems algorithms based methods and the quantum simulation methods arising from different numerical approximations, including the finite difference discretisations and the Fourier spectral discretisations for the two different linear representations, with the result showing that the quantum simulation methods usually give the best performance in time complexity. We also propose the Schrödinger framework to solve the Liouville equation for the HJE, since it can be recast as the semiclassical limit of the Wigner transform of the Schrödinger equation. Comparsion between the Schrödinger and the Liouville framework will also be made.
Solving the time-dependent Schrödinger equation is an important application area for quantum algo... more Solving the time-dependent Schrödinger equation is an important application area for quantum algorithms. We consider Schrödinger's equation in the semi-classical regime. Here the solutions exhibit strong multiple-scale behavior due to a small parameter ℏ, in the sense that the dynamics of the quantum states and the induced observables can occur on different spatial and temporal scales. Such a Schrödinger equation finds many applications, including in Born-Oppenheimer molecular dynamics and Ehrenfest dynamics. This paper considers quantum analogues of pseudo-spectral (PS) methods on classical computers. Estimates on the gate counts in terms of ℏ and the precision ε are obtained. It is found that the number of required qubits, m, scales only logarithmically with respect to ℏ. When the solution has bounded derivatives up to order ℓ, the symmetric Trotting method has gate complexity O((εℏ)−12polylog(ε−32ℓℏ−1−12ℓ)), provided that the diagonal unitary operators in the pseudo-spectral ...
This thesis is an exploration of the power of photonic resources, as viewed from several differen... more This thesis is an exploration of the power of photonic resources, as viewed from several different but related perspectives. They range from quantum computation, precision parameter estimation to the thermodynamics of relativistic quantum systems, as applied to cosmology in particular. The use of photonic states allows us to address several important questions about the resources required in quantum mechanical processes.
In chapter 1, we propose a new quantum computational model, called the `power of one qumode', that relies mainly on a single-mode photonic squeezed state. In particular, we show the amount of squeezing can quantitatively relate the resource requirements of factoring to the problem of finding the trace of large unitary matrices, a result with consequences for understanding how powerful quantum computation really is. Furthermore, we can connect squeezing to other known resources like precision, energy, qudit dimensionality and qubit number, which is a useful stepping stone to finding the resources that enable quantum computation.
In chapter 2, we exploit the quantum mechanical properties of photonic states for use in precision parameter estimation of general linear optical processes, which is useful for a diverse number of applications, from characterising an unknown process in a photonic quantum computer to biological imaging. We introduce a formalism that quantifies this improvement in precision. We also provide conditions under which one can easily check for photonic states that are optimal to use in this context, which is a potentially important result for future experimental efforts.
In chapter 3, we explore the connection between two-mode squeezed states, commonly used in quantum optics, and relativistic quantum processes, in particular in cosmology. Using this connection, we apply recently developed tools from the thermodynamics of quantum systems perturbed far from equilibrium to address an old question of entropy production in cosmology from a surprising new angle.
Uploads
Papers by Nana Liu
This diversity already suggests the potential for varied applications, as well as acting as a bridge between different domains. Here we devise new quantum subroutines to efficiently prepare quantum unitary operators that embed the standard matrix geometric mean and its generalisations called the weighted matrix geometric mean. This enables the construction of solutions to the algebraic Riccati equation, which is an important class of nonlinear systems of equations that appears in machine learning, optimal control, estimation, and filtering.
Using these subroutines,
we present a new class of quantum learning algorithms called quantum geometric mean metric learning. This has applications in efficiently finding the best distance measure and solving classification problems in the weakly supervised limit and for anomaly detection, for both classical and quantum problems. We also show how our method can be generalised to a particular p^th-order system of nonlinear equations.
These quantum subroutines for matrix geometric means are also useful in other areas of quantum information. For example, we show how to use them in the estimation of geometric Renyi relative entropies and the Uhlmann fidelity by means of the Fuchs--Caves observable.
In particular, our quantum algorithms for estimating the Uhlmann and Matsumoto fidelities have optimal dependence on the precision.
Finally, we provide a BQP-complete problem based on matrix geometric means that can be solved by our subroutines, thus characterising their computational capability.
This diversity already suggests the potential for varied applications, as well as acting as a bridge between different domains. Here we devise new quantum subroutines to efficiently prepare quantum unitary operators that embed the standard matrix geometric mean and its generalisations called the weighted matrix geometric mean. This enables the construction of solutions to the algebraic Riccati equation, which is an important class of nonlinear systems of equations that appears in machine learning, optimal control, estimation, and filtering.
Using these subroutines,
we present a new class of quantum learning algorithms called quantum geometric mean metric learning. This has applications in efficiently finding the best distance measure and solving classification problems in the weakly supervised limit and for anomaly detection, for both classical and quantum problems. We also show how our method can be generalised to a particular p^th-order system of nonlinear equations.
These quantum subroutines for matrix geometric means are also useful in other areas of quantum information. For example, we show how to use them in the estimation of geometric Renyi relative entropies and the Uhlmann fidelity by means of the Fuchs--Caves observable.
In particular, our quantum algorithms for estimating the Uhlmann and Matsumoto fidelities have optimal dependence on the precision.
Finally, we provide a BQP-complete problem based on matrix geometric means that can be solved by our subroutines, thus characterising their computational capability.
In chapter 1, we propose a new quantum computational model, called the `power of one qumode', that relies mainly on a single-mode photonic squeezed state. In particular, we show the amount of squeezing can quantitatively relate the resource requirements of factoring to the problem of finding the trace of large unitary matrices, a result with consequences for understanding how powerful quantum computation really is. Furthermore, we can connect squeezing to other known resources like precision, energy, qudit dimensionality and qubit number, which is a useful stepping stone to finding the resources that enable quantum computation.
In chapter 2, we exploit the quantum mechanical properties of photonic states for use in precision parameter estimation of general linear optical processes, which is useful for a diverse number of applications, from characterising an unknown process in a photonic quantum computer to biological imaging. We introduce a formalism that quantifies this improvement in precision. We also provide conditions under which one can easily check for photonic states that are optimal to use in this context, which is a potentially important result for future experimental efforts.
In chapter 3, we explore the connection between two-mode squeezed states, commonly used in quantum optics, and relativistic quantum processes, in particular in cosmology. Using this connection, we apply recently developed tools from the thermodynamics of quantum systems perturbed far from equilibrium to address an old question of entropy production in cosmology from a surprising new angle.