Numerical Modelling of Some Quantum Systems (original) (raw)

An introduction to the problem of bridging quantum and classical dynamics

The European Physical Journal Special Topics, 2015

Simulating the exact quantum dynamics of realistic interacting systems is presently a task beyond reach but for the smallest of them, as the numerical cost for solving the time-dependent Schrödinger equation scales exponentially with the number of degrees of freedom. Mixed quantum-classical methods attempt to solve this problem by starting from a full quantum description of the system and subsequently partitioning the degrees of freedom in two subsets: the quantum subsystem and the bath. A classical limit is then taken for the bath while preserving, at least approximately, the quantum evolution of the subsystem. A key, as yet not fully resolved, theoretical question is how to do so by constructing a consistent description of the overall dynamics. An exhaustive review of this class of methods is beyond the scope of this paper and we shall limit ourselves to present, as an example, a specific approach, known as the LANDM-Map method. The method stems from an attempt at taking a rigorous limit for the classical degrees of freedom starting from a path integral formulation of the full quantum problem. The results that we discuss are not new, but our intent here is to present them as an introduction to the problem of mixed quantum classical dynamics. We shall also indicate a broad classification of the available approaches, their limitations, and some open questions in this field.

Modelling of quantum mechanical control systems

Mathematical Modelling

A program of studies of quantum mechanical control systems is initiated. Following the historical development of quantum mechanics, the quantum control model is obtained from a corresponding classical structure. Second order linear and bilinear control systems as well as first order linear control systems are investigated; it is further demonstrated that the analysis may be extended to some nonlinear control problems. The results derived for these systems form an interesting example of the general theory of quantization.

Quantum simulation of partial differential equations: Applications and detailed analysis

Physical Review A

We study a recently introduced simple method [S. Jin, N. Liu, and Y. Yu, Quantum simulation of partial differential equations via Schrödingerisation, arXiv:2212.13969] for solving general linear partial differential equations with quantum simulation. This method converts linear partial differential equations into a Hamiltonian system, using a simple transformation called the warped phase transformation. Here we provide a more-in-depth technical discussion and expand on this approach in a more detailed and pedagogical way. We apply this to examples of partial differential equations, including heat, convection, Fokker-Planck, linear Boltzmann, and Black-Scholes equations. This approach can also be extended to general linear partial differential equations, including the Vlasov-Fokker-Planck equation and the Liouville representation equation for nonlinear ordinary differential equations. Extension to higher-order time derivatives is also possible.

Quantum simulation of discrete linear dynamical systems and simple iterative methods in linear algebra via Schrödingerisation

ArXiv:2304.02865, 2023

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 in the semi-classical regime

Quantum

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 ...

Quantum Simulation for Quantum Dynamics with Artificial Boundary Conditions

arXiv (Cornell University), 2023

Quantum dynamics, typically expressed in the form of a time-dependent Schrödinger equation with a Hermitian Hamiltonian, is a natural application for quantum computing. However, when simulating quantum dynamics that involves the emission of electrons, it is necessary to use artificial boundary conditions (ABC) to confine the computation within a fixed domain. The introduction of ABCs alters the Hamiltonian structure of the dynamics, and existing quantum algorithms can not be directly applied since the evolution is no longer unitary. The current paper utilizes a recently introduced Schrödingerisation method [JLY22a, JLY22b] that converts non-Hermitian dynamics to a Schrödinger form, for the artificial boundary problems. We implement this method for three types of ABCs, including the complex absorbing potential technique, perfectly matched layer methods, and Dirichlet-to-Neumann approach. We analyze the query complexity of these algorithms, and perform numerical experiments to demonstrate the validity of this approach. This helps to bridge the gap between available quantum algorithms and computational models for quantum dynamics in unbounded domains.

Discrete time quantum mechanics

Computers & Mathematics with Applications, 1994

paper summarizes a research program that has been underway for a decade. The objective is to find a fast and accurate scheme for solving quantum problems which does not involve a Monte Carlo algorithm.