Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE

Simulation of quantum dynamics by difference schemes

A. P. Makkaveev, D. I. Pomozov

Published Jun 1, 2005• 0 citations

Abstract

We represent program for the simulation of evolution one and the two-particle quantum systems on the classical computer. Program is the part of the work of our department on the program packet of simulation. Simulation is produced by the numerical solution of the Schroedinger equation with the use of finite difference schemes.

DOI: 10.1117/12.620524

Related papers

Quantum finite difference solvers for physics simulation

Anthony Chagneau, Laëticia Nathoo, Jérémy Alloul +1 more

Abstract Physics systems are becoming increasingly complex and require more and more computing time. Quantum computing, which has shown its efficiency on some problems, such as the factorisation of a number with Shor's algorithm, may be the solution to reduce these computation times. Here, the authors propose two quantum numerical schemes for the simulation of physics phenomena, based on the finite difference method. The aim is to see if quantum versions of standard numerical schemes offer an advantage over their classical counterparts, either in accuracy, stability or computation time. First, the authors will present the different phenomena studied as well as the classical solution methods chosen. The authors will then describe the implementation of the quantum numerical schemes and present some results obtained on the different physics phenomena beforehand and then compare both approaches, classical and quantum.

IET Quantum CommunicationMar 3, 20234 citations

Numerical Quantum Dynamics

Emil Kieri

We consider computational methods for simulating the dynamics of molecular systems governed by the time-dependent Schrödinger equation. Solving the Schrödinger equation numerically poses a challenge due to its often highly oscillatory solutions, and to the exponential growth of work and memory with the number of particles in the system. Two different classes of problems are studied: the dynamics of the nuclei in a molecule and the dynamics of an electron in orbit around a nucleus. For the first class of problems we present new computational methods which exploit the relation between quantum and classical dynamics in order to make the computations more efficient. For the second class of problems, the lack of regularity in the solution poses a computational challenge. Using knowledge of the non-smooth features of the solution we construct a new method with two orders higher accuracy than what is achieved by direct application of a difference stencil.

Jan 1, 20130 citations

Controlled Quantum Packets

Salvatore DeMartino, Silvio DeSiena

We look at time evolution of a physical system from the point of view of dynamical control theory. Normally we solve motion equation with a given external potential and we obtain time evolution. Standard examples are the trajectories in classical mechanics or the wave functions in Quantum Mechanics. In the control theory, we have the configurational variables of a physical system, we choose a velocity field and with a suited strategy we force the physical system to have a well defined evolution. The evolution of the system is the 'premium' that the controller receives if he has adopted the right strategy. The strategy is given by well suited laboratory devices. The control mechanisms are in many cases non linear; it is necessary, namely, a feedback mechanism to retain in time the selected evolution. Our aim is to introduce a scheme to obtain Quantum wave packets by control theory. The program is to choose the characteristics of a packet, that is, the equation of evolution for its centre and a controlled dispersion, and to give a building scheme from some initial state (for example a solution of stationary Schroedinger equation). It seems natural in this view to use stochastic approach to Quantum Mechanics, that is, Stochastic Mechanics [S.M.]. It is a quantization scheme different from ordinary ones only formally. This approach introduces in quantum theory the whole mathematical apparatus of stochastic control theory. Stochastic Mechanics, in our view, is more intuitive when we want to study all the classical-like problems. We apply our scheme to build two classes of quantum packets both derived generalizing some properties of coherent states.

Jan 1, 19960 citations

Computational Methods for Simulating Quantum Computers

Hans A. De Raedt, Kristel F. L. Michielsen

This review gives a survey of numerical algorithms and software to simulate quantum computers. It covers the basic concepts of quantum computation and quantum algorithms and includes a few examples that illustrate the use of simulation software for ideal and physical models of quantum computers.

arXiv:quant-ph/0406210Jun 27, 200436 citations

Adaptive kernel methods to simulate quantum phase space flow

Arnaldo Donoso, Laboratorio de F ́ isica

A technique for simulating quantum dynamics in phase space i discussed. It makes use of ensembles of classical trajectories to approximate the distribution functions and their derivatives by implementing Adaptive Kernel Density Estimation. It is found to improve the accuracy and stability of the simulations compared to more conventional particle methods. Formulation of the method in higher dimensions is straightforward. Complete analytical description of physical systems exist only for the most simple ideal cases. The vast majority of physical problems require the solution of approximate equations most fre- quently by using numerical calculations. To this end, much effort has been devoted to the devel- opment of suitable techniques to be used in computer simulations. The increasing computational power makes possible the description of more complicated, that is, less ideal systems allowing for more precise predictions of their properties as well as for a better understanding of the underlying phenomena. Among the wealth of numerical approaches to solve dynamical equations by means of discretization, particle methods have gained attention in the last few decades as a conceptu- ally appealing alternative. Conventional numerical solvers rely on a regular discretization of the independent variables to facilitate the evaluation of derivatives to a given order of accuracy. As a consequence, a mesh needs to be defined as the base of the simulations, parts of which remain empty at any given simulation time even with the most sophisticated mesh schemes. The cen- tral idea of particle methods consists of representing functions by sets of particles conveniently located and weighted such that the function is appropriately sampled by the particles. Evolution of the functions is then accomplished by an equivalent evolution of the particles. This meshless or Lagrangian type of discretization has some intrinsic computational efficiency advantages because there is no wasted space. The advantage, however, comes at the expense of having to evaluate derivatives of functions on unstructured arrays of sampling points. The issues about convergence and stability of the resulting algorithms have been the subject of active and fruitfull research. Currently there are some well established techniques to solve dynamical equations with particle methods. Our interest is to provide an alternative ingredient to further improve the performance of numerical solutions of dynamical problems of physical relevance. In the present work we discuss, in the context of particle methods, the implementation of a technique for representing functions and derivatives with adaptive parameters. Adaptiveness require the parameters to be automatically ad- justed to optimal values according to the location of the particles. As shown below, the proposed scheme is conceptually simple and has been tested in a number of different situations. In some simple test models the present work focuses on quantum dynamics, previously treated with more

Jan 1, 20060 citations

Alternative descriptions of wave and particle aspects of the harmonic oscillator

Dieter Schuch

The dynamical properties of the wave and particle aspects of the harmonic oscillator can be studied with the help of the time-dependent Schroedinger equation (SE). Especially the time-dependence of maximum and width of Gaussian wave packet solutions allow to show the evolution and connections of those two complementary aspects. The investigation of the relations between the equations describing wave and particle aspects leads to an alternative description of the considered systems. This can be achieved by means of a Newtonian equation for a complex variable in connection with a conservation law for a nonclassical angular momentum-type quantity. With the help of this complex variable, it is also possible to develop a Hamiltonian formalism for the wave aspect contained in the SE, which allows to describe the dynamics of the position and momentum uncertainties. In this case the Hamiltonian function is equivalent to the difference between the mean value of the Hamiltonian operator and the classical Hamiltonian function.

Jan 1, 19930 citations