We introduce the “tensor-train split-operator Fourier transform” (TT-SOFT) method for simulations of multidimensional nonadiabatic quantum dynamics. TT-SOFT is essentially the grid-based SOFT method implemented in dynamically adaptive tensor-train representations. In the same spirit of all matrix product states, the tensor-train format enables the representation, propagation, and computation of observables of multidimensional wave functions in terms of the grid-based wavepacket tensor components, bypassing the need of actually computing the wave function in its full-rank tensor product grid space. We demonstrate the accuracy and efficiency of the TT-SOFT method as applied to propagation of 24-dimensional wave packets, describing the S 1 / S 2 interconversion dynamics of pyrazine after UV photoexcitation to the S 2 state. Our results show that the TT-SOFT method is a powerful computational approach for simulations of quantum dynamics of polyatomic systems since it avoids the exponential scaling problem of full-rank grid-based representations.
Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.
Abstract In the present paper, we propose and analyse a class of tensor methods for the efficient numerical computation of the dynamics and spectrum of high-dimensional Hamiltonians. We focus on the complex-time evolution problems. We apply the quantized-TT (QTT) matrix product states type tensor approximation that allows to represent N-d tensors generated by the grid representation of d-dimensional functions and operators with log-volume complexity, O(d log N), where N is the univariate discretization parameter in space. Making use of the truncated Cayley transform method allows us to recursively separate the time and space variables and then introduce the efficient QTT representation of both the temporal and the spatial parts of the solution to the high-dimensional evolution equation. We prove the exponential convergence of the m-term time-space separation scheme and describe the efficient tensor-structured preconditioners for the arising system with multidimensional Hamiltonians. For the class of "analytic" and low QTT-rank input data, our method allows to compute the solution at a fixed point in time t=T>0 with an asymptotic complexity of order O(d log N ln^q (1/ε)), where ε>0 is the error bound and q is a fixed small number. The time-and-space separation method via the QTT-Cayley-transform enables us to construct a global m-term separable (x,t)-representation of the solution on a very fine time-space grid with complexity of order O(dm^4 log N_t log N), where N_t is the number of sampling points in time. The latter allows efficient energy spectrum calculations by FFT (or QTT-FFT) of the autocorrelation function computed on a sufficiently long time interval [0,T]. Moreover, we show that the spectrum of the Hamiltonian can also be represented by the poles of the t-Laplace transform of a solution. In particular, the approach can be an option to compute the dynamics and the spectrum in the time-dependent molecular Schrödinger equation.
In the present paper, we propose and analyse a class of tensor methods for the efficient numerical computation of dynamics and spectrum of high-dimensional Hamiltonians. The main focus is on the complex-time evolution problems. We apply the recent quantized-TT (QTT) matrix product states tensor approximation that allows to represent N -d tensors generated by grid representation of d-dimensional functions and operators with log-volume complexity, O(d logN), where N is the univariate discretization parameter in space. Making use of the truncated Cayley transform method allows us to recursively separate the time and space variables and then introduce the efficient QTT representation of both the temporal and spatial parts of solution to the high-dimensional evolution equation. We prove the exponential convergence of the m-term time-space separation scheme and describe the efficient tensor-structured preconditioners for the arising system with multidimensional Hamiltonians. For a class of “analytic” and low QTT-rank input data, our method allows to compute the solution at a fix point in time t = T > 0 with the asymptotic complexity of order O(d logN ln 1e ), where e > 0 is the error bound and q a fixed small number. The time-and-space separation method via the QTT-Cayley-transform enables us to construct the global m-term separable (x, t)-representation of a solution on very fine time-space grid with complexity of order O(dm4 logNt logN), where Nt is the number of sampling points in time. The latter allows the efficient energy spectrum calculations by FFT (or QTTFFT) of autocorrelation function computed on sufficiently long time interval [0, T ]. Moreover, we show that the spectrum of Hamiltonian can also be represented by the poles of the t-Laplace transform of a solution. In particular, the approach can be an option to compute dynamics and spectrum in the time-dependent molecular Schrodinger equation. AMS Subject Classification: 65F30, 65F50, 65N35, 65F10
We develop a dynamical tensor-train method for the Wigner--Fokker--Planck equation arising in continuous-variable open quantum systems. The Wigner function is evolved directly in tensor-train form using a projector-splitting integrator, while the nonlocal Wigner pseudo-differential operator is treated through low-rank representations of the potential difference and matrix-free exponential actions. The resulting method avoids full tensor-product computations in high-dimensional phase space. Numerical experiments in one, two, and three spatial dimensions demonstrate the accuracy and scalability of the method, including simulations in a six-dimensional phase space.
Quantification of stochastic or quantum systems by a joint probability density or wave function is a notoriously difficult computational problem, since the solution depends on all possible states (or realizations) of the system. Due to this combinatorial flavor, even a system containing as few as ten particles may yield as many as
1010 discretized states. None of even modern supercomputers are capable to cope with this curse of dimensionality straightforwardly, when the amount of quantum particles, for example, grows up to more or less interesting order of hundreds. A traditional approach for a long time was to avoid models formulated in terms of probabilistic functions, and simulate particular system realizations in a randomized process. Since different times in different communities, data-sparse methods came into play. Generally, they aim to define all data points indirectly, by a map from a low amount of representers, and recast all operations (e.g. linear system solution) from the initial data to the effective parameters. The most advanced techniques can be applied (at least, tried) to any given array, and do not rely explicitly on its origin. The current work contributes further progress to this area in the particular direction: tensor product methods for separation of variables. The separation of variables has a long history, and is based on the following elementary concept: a function of many variables may be expanded as a product of univariate functions. On the discrete level, a function is encoded by an array of its values, or a tensor. Therefore, instead of a huge initial array, the separation of variables allows to work with univariate factors with much less efforts. The dissertation contains a short overview of existing tensor representations: canonical PARAFAC, Hierarchical Tucker, Tensor Train (TT) formats, as well as the artificial tensorisation, resulting in the Quantized Tensor Train (QTT) approximation method. The contribution of the dissertation consists in both theoretical constructions and practical numerical algorithms for high-dimensional models, illustrated on the examples of the Fokker-Planck and the chemical master equations. Both arise from stochastic dynamical processes in multiconfigurational systems, and govern the evolution of the probability function in time. A special focus is put on time propagation schemes and their properties related to tensor product methods. We show that these applications yield large-scale systems of linear equations, and prove analytical separable representations of the involved functions and operators. We propose a new combined tensor format (QTT-Tucker), which descends from the TT format (hence TT algorithms may be generalized smoothly), but provides complexity reduction by an order of magnitude. We develop a robust iterative solution algorithm, constituting most advantageous properties of the classical iterative methods from numerical analysis and alternating density matrix renormalization group (DMRG) techniques from quantum physics. Numerical experiments confirm that the new method is preferable to DMRG algorithms. It is as fast as the simplest alternating schemes, but as reliable and accurate as the Krylov methods in linear algebra.
Methods for efficient simulations of multidimensional quantum dynamics are essential for theoretical studies of chemical systems where quantum effects are important, such as those involving rearrangements of protons or electronic configurations. Here, we introduce the functional tensor-train Chebyshev (FTTC) method for rigorous nuclear quantum dynamics simulations. FTTC is essentially the Chebyshev propagation scheme applied to the initial state represented in a continuous analogue tensor-train format. We demonstrate the capabilities of FTTC as applied to simulations of proton quantum dynamics in a 50-dimensional model of hydrogen-bonded DNA base pairs.