A numerical method was proposed to propagate the quantum system with a time-dependent Hamiltonian. The propagator is composed of the time-independent-Hamiltonian-type multisplit forms, which becomes possible by realizing that the evolution operator has a close link with the time-shift operator. The present algorithm can achieve a high order of accuracy through the recurrence of multisplit forms with predetermined coefficients, and it has the advantage that only the Hamiltonian itself at different time instants is involved, while other currently available schemes need to perform either a time integration or a time differential of the Hamiltonian. This essential difference has made our proposal highly competitive. A one-dimensional model system is considered to verify the accuracy and efficiency of the numerical scheme.
A numerical method was proposed to propagate the quantum system with a time-dependent Hamiltonian. The propagator is composed of the time-independent-Hamiltonian-type multisplit forms, which becomes possible by realizing that the evolution operator has a close link with the time-shift operator. The present algorithm can achieve a high order of accuracy through the recurrence of multisplit forms with predetermined coefficients, and it has the advantage that only the Hamiltonian itself at different time instants is involved, while other currently available schemes need to perform either a time integration or a time differential of the Hamiltonian. This essential difference has made our proposal highly competitive. A one-dimensional model system is considered to verify the accuracy and efficiency of the numerical scheme.
We propose a new method for simulating certain type of time-dependent Hamiltonian H(t)=∑i=1mγi(t)Hi where γi(t) (and its higher order derivatives) is bounded, computable function of time t, and each Hi is time-independent, and could be efficiently simulated. Our quantum algorithms are based on high-order Runge-Kutta method and forward Euler method, where the time interval is divided into subintervals. Then in an iterative manner, the evolution operator at given time step is built upon the evolution operator at previous time step, utilizing algorithmic operations from the recently introduced quantum singular value transformation framework.
Effective modeling and numerical spectral-based propagation schemes are proposed for addressing the challenges in time-dependent quantum simulations of systems ranging from atoms, molecules, and nanostructures to emerging nanoelectronic devices. While time-dependent Hamiltonian problems can be formally solved by propagating the solutions along tiny simulation time steps, a direct numerical treatment is often considered too computationally demanding. In this paper, however, we propose to go beyond these limitations by introducing high-performance numerical propagation schemes to compute the solution of the time-ordered evolution operator. In addition to the direct Hamiltonian diagonalizations that can be efficiently performed using the new eigenvalue solver FEAST, we have designed a Gaussian propagation scheme and a basis-transformed propagation scheme (BTPS) which allow to reduce considerably the simulation times needed by time intervals. It is outlined that BTPS offers the best computational efficiency allowing new perspectives in time-dependent simulations. Finally, these numerical schemes are applied to study the ac response of a (5,5) carbon nanotube within a three-dimensional real-space mesh framework.
The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.
We describe an efficient path integral scheme for calculating the propagator of an arbitrary quantum system, as well as that of a stochastic system in special cases where the Fokker–Planck equation obeys strict detailed balance. The basic idea is to split the respective Hamiltonian into two exactly solvable parts and then to employ a symmetric decomposition of the time evolution operator, which is exact up to a high order in the time step. The resulting single step propagator allows rather large time steps in a path integral and leads to convergence with fewer time slices. Because it involves no system-specific reference system, the algorithm is amenable to all known numerical schemes available for evaluating quantum path integrals. In this way one obtains a highly accurate method, which is simultaneously fast, stable, and computationally simple. Numerical applications to the real time quantum dynamics in a double well and to the stochastic dynamics of a bistable system coupled to a harmonic mode show our method to be superior over the approach developed by the Makri group in their quasiadiabatic propagator representation, to say nothing about the propagation scheme based on the standard Trotter splitting.
This talk is about a method for solving the numerical initial value problem of the time-dependent Schrodinger equation with a Hamiltonian which may depend on time. It is only for time independent Hamiltonians that I was able to describe the accuracy and stability properties in the form of proven theorems. All practical experience shows, that the time-dependent case behaves equally well. In this method, the Hamiltonian comes in only as a hermitian linear operator and not as a mathematical expression that can be reordered or otherwise transformed in order to speed up the necessary calculations. In all cases where it is a feasible computational task to apply the Hamilton operator to a generic state vector, my suggestion is to try this method as a workhorse. It can be considered a ‘Simpson’s rule’ for quantum dynamics. There is no way for me to explain the matter from a technical point of view in 10 minutes, so I’ll try to convey the basic idea of the method. This idea is to treat quantum dynamics in closer analogy to classical dynamics as one usually does. One may argue, that there is a fundamental difference between classical and quantum dynamics since the Schrodinger equation is first order with respect to the time derivative, whereas classical mechanics is inherently second order, since forces determine accelerations and not velocities. This is not so clear a distinction as one might think, since also classical mechanics can be given the first-order form of canonical equations. I propose the opposite of such an approach, namely to give quantum dynamics a second-order form by simply differentiating the time-dependent Schrodinger equation with respect to time. To be sure, this is only a heuristic device for defining a computationally feasible approximation and not a proposal to change quantum mechanics. For the resulting equation for the second time derivative of the wave function, I then carry over and use the integration methods from the field of granular matter dynamics which I have been using in an industrial project concerning simulation of the toning process in electro-photographic copiers/printers over some years [4]. Such methods come in
In this study, a variety of methods are tested and compared for the numerical solution of the Schrödinger equation for few-body systems with explicitely time-dependent Hamiltonians, with the aim to find the optimal one. The configuration interaction method, generally applied to find stationary eigenstates accurately and without approximations to the wavefunction's structure, may also be used for the time-evolution, which results in a large linear system of ordinary differential equations. The large basis sizes typically present when the configuration interaction method is used calls for efficient methods for the time evolution. Apart from efficiency, adaptivity (in the time domain) is the other main focus in this study, such that the time step is adjusted automatically given some requested accuracy. A method is suggested here, based on an exponential integrator approach, combined with different ways to implement the adaptivity, which was found to be faster than a broad variety of other methods that were also considered.