Split-operator technique for propagating phase space functions: Exploring chaotic, dissipative and relativistic dynamics
Introduction
The split-operator technique (SOT) has been considered one of the most plausible alternatives to perform the time-evolution of localized initial states in quantum systems. It has been successfully applied for the solution of the Schrödinger equation in coordinates as well as phase space representations [1], [2], [3], [4]. It has become an ideal for propagating initial states in many branches of physics, for instance, applications involving wave-packet dynamics in molecular systems [5], [6], [7] and frequent use of quantum molecular dynamics [8], [9]. In the recent past, Finite difference has adopted the SOT scheme to simulate propagation phenomena in quantum mechanics and quantum field theory [10], [11]. Moreover, the versatility of this hybrid model becomes an advantage of using both modeling schemes; Noticeably, to simulate the scattering of sound in the time-dependent domain [12], and the propagation of waves through optical fibers [13]. Recently, the SOT has exploited the power of high performance computing and it has been implemented in modern hardware architectures for solving time-dependent Schrödinger equations and time-dependent Dirac equations [14] as well as the Klein–Gordon equation [15]. In contrast to quantum mechanics, little is known about the benefits of this approach in the context of classical mechanics, in fact, there is only one work in the literature where this issue has been addressed in detail. More precisely, the pioneering work of Dattoli et al. [16] has developed a numerical method based on SOT for propagating the Liouville equation, and it was successfully applied for the case of a one-dimensional Hamiltonian system. Nevertheless, the numerical approach has not been widely noticed in many significant situations of classical dynamics. Possibly, it is due to their work having a limitation of implementation, the initial state should be known in an exact form as an inevitable step to perform the propagation in phase space. Furthermore, the traditional propagation methods in classical mechanics continue being the favorites to propagate initial states, since they involve mainly three stages, namely: (i) the initial state is frequently sampled on phase space accordingly with a particular distribution function, (ii) the corresponding classical trajectories are propagated at final time using a numerical integrator for solving the classical equations of motion; thus, those propagated trajectories contribute to the final state.(iii) In general, the propagated state should be smoothed in a grid, since all classical trajectories contribute in a different way to each phase space point. This methodology is commonly used in classical molecular dynamics [17], [18], semiclassical propagation of wave-packets and Wigner functions [19], [20], [21], as well as studies about quantum–classical correspondence where a comparison with classical calculations is necessary [22], [23]. However, it is worth mentioning that some complications can arise in the traditional propagation methods, for instance, the sampling method cannot always cover completely the initial state, and the convergence of methods as well as the propagated state will be a compromise. The numerical integration generally is time consuming when non-linearities appear in the classical equations of motion, and the propagated state requires a lot of classical trajectories in order to obtain a good resolution in phase space.
The present paper is concerned with classical phase space propagation based on the split-operator technique, and one of the motivations of this work is to fill the knowledge gap on this isolated topic. In order to achieve this, we will adopt the SOT to different scenarios of classical (relativistic) mechanics involving; regular, chaotic motion, and dissipative dynamics, and subsequently, propagation in high-dimensional phase spaces. We propose a numerical method for propagating phase space distribution functions with several advantages in contrast to the traditional methods. More precisely, (i) our approach is directly based on computations on a grid of size representing the whole phase space of classical systems, thus the resolution in phase space is defined at the beginning of the calculation, and it can be managed with the actual computational resources. (ii) The initial state should be storage on grid, but it can be computed in exact form or sampled as is done by the traditional methods, and finally (iii) the numerical propagation is accomplished by the fast Fourier transform where the computation method takes the order of operations that will reduce significantly the computational cost of the numerical simulations.
This paper is organized as follows: Section 2 explains the theoretical background of propagation phase space distribution functions according to the Liouville equation, and the numerical solution based on fast Fourier transform. In Section 3 we describe the method for propagating Hamiltonian systems with -kick interaction. In Sections 4 Propagation of non-Hermitian quantum dynamics in the classical limit, 5 Propagation of a relativistic particle is considered for the time evolution of phase space distribution functions for dissipative dynamics associated with the classical limit of non-Hermitian quantum systems, and relativistic systems, respectively. The Section 6 deal with the phase space propagation for an interacting one-dimensional classical systems. Finally, we conclude in Section 7.
Section snippets
Liouville equation and the numerical schema
The time-evolution of phase space distribution function for a one-dimensional system is governed by the Liouville equation where is the phase space distribution function and is the Liouville operator. The latter is given by where is the classical Hamiltonian function. The solution to Eq. (1) can formally be written as , where is the phase space distribution function at initial time, and the evolution
Propagation of -kicked systems
In the literature, impulsively driven systems are often referred to -kicked systems, due to the fact that the duration of each impulse is infinitely small compared to the period of the unperturbed system. These systems have become an essential tool for studying many phenomena from both theoretical and experimental perspectives, for instance, one of the primary approaches is quantum chaos in cold atoms [24], [25], which ensures that it can provide valuable insights into non-linear-Hamiltonian
Propagation of non-Hermitian quantum dynamics in the classical limit
For decades considerable attention has been paid to describe the non-Hermitian Hamiltonian systems in the context of quantum mechanics. In particular, existing theoretical approaches provide an instructive description of open quantum systems, specifically these are subjected to systems involving decay, scattering and transport phenomena [30], [31], [32]. However, little is known about the quantum–classical correspondence of non-Hermitian Hamiltonian systems, and recent theoretical studies have
Propagation of a relativistic particle
Recently, the problem of phase space evolution for a relativistic system has been addressed in the framework of Hamiltonian mechanics [40], more precisely, for the relativistic harmonic oscillator. In this case, the propagation is achieved by considering the method proposed in Ref. [16] and it has successfully applied for an initially Gaussian distribution. Nonetheless, it is in general not amenable for non analytical initial states. In order to illustrate our numerical method, we propagate
Propagation of an interacting Hamiltonian system
In the recent past, many theoretical efforts have been carried out to understand one of the most intriguing quantum phenomena, entanglement, and many authors have explored different measures in order to characterize it. The most common measures for interacting (biparticle) systems are the entropies, which gives the measure of non-separability of individual subsystems. In spirit, to understand the classical counterpart of this phenomenon very few studies have been carried out in the past [41],
Conclusions
We have adapted the SOT for propagating phase space distribution functions in scenarios previously not considered in the literature. The numerical approach is rather general and allows to propagate almost any localized initial state, and it is not necessary to know the state in exact form as is the case of Ref. [16]. We have considered the time evolution of initial states in the framework of classical and relativistic mechanics, in particular we have shown the performance of the method in
Acknowledgments
The first author thanks the Grupo de Materiales Magnéticos & Nanoestructuras of Universidad del Quindío for their generous provision of computer facilities, and he acknowledges the financial support from Universidad del Quindío.
References (44)
- et al.
J. Comput. Phys.
(1982) - et al.
Physica D
(1998) - et al.
Comput. Phys. Commun.
(2011) - et al.
J. Comput. Phys.
(2009) Chem. Phys.
(2002)- et al.
Commun. Nonlinear Sci. Numer. Simul.
(2003) - et al.
Comput. Phys. Commun.
(2009) - et al.
Physica A
(2004) - et al.
Phys. Rev. Lett.
(1991) - et al.
J. Chem. Phys.
(1982)
J. Chem. Phys.
Comput. Phys. Commun.
J. Chem. Phys.
J. Chem. Phys.
Annu. Rev. Phys. Chem.
La Rivista del Nuovo Cimento
Z. Kristallogr.
J. Lightwave Technol.
Phys. Rev. E
J. Chem. Phys.
J. Phys. Chem. A
Cited by (7)
Phase-space studies of backscattering diffraction of defective Schrödinger cat states
2021, Scientific ReportsPhase-space approach to time evolution of quantum states in confined systems. The spectral split-operator method
2020, Advances in Intelligent Systems and ComputingThe Effect of Elastic and Inelastic Scattering on Electronic Transport in Open Systems
2019, International Journal of Applied Mathematics and Computer ScienceThe Phase-Space Approach to time Evolution of Quantum States in Confined Systems: The Spectral Split-Operator Method
2019, International Journal of Applied Mathematics and Computer ScienceDissipative transport of thermalized electrons through a nanodevice
2017, Physical Review B