Related Experiment Video
Updated: Feb 3, 2026

Excitonic Hamiltonians for Calculating Optical Absorption Spectra and Optoelectronic Properties of Molecular Aggregates and Solids
Published on: May 27, 2020
Stability of velocity-Verlet- and Liouville-operator-derived algorithms to integrate non-Hamiltonian systems
1The Institute for Solid State Physics, The University of Tokyo, Kashiwanoha 5-1-5, Kashiwa, Chiba 277-8581, Japan.
Abstract:
We investigate the difference between the velocity Verlet and the Liouville-operator-derived (LOD) algorithms by studying two non-Hamiltonian systems, one dissipative and the other conservative, for which the Jacobian of the transformation can be determined exactly. For the two systems, we demonstrate that (1) the velocity Verlet scheme fails to integrate the former system while the first- and second-order LOD schemes succeed and (2) some first-order LOD fails to integrate the latter system while the velocity Verlet and the other first- and second-order schemes succeed. We have shown that the LOD schemes are stable for the former system by determining the explicit forms of the shadow Hamiltonians which are exactly conserved by the schemes. We have shown that the Jacobian of the velocity Verlet scheme for the former system and that of the first-order LOD scheme for the latter system are always smaller than the exact values, and therefore, the schemes are unstable. The decomposition-order dependence of LOD schemes is also considered.
Related Concept Videos
Velocity and Position by Integral Method
Consider an example to calculate the velocity and position from the acceleration function. A motorboat is traveling at a constant velocity of 5.0 m/s when it starts to decelerate to arrive at the dock. Its acceleration is...
Second Derivatives and Laplace Operator
Consider a scalar function. The curl of its...
BIBO stability of continuous and discrete -time systems
To determine the BIBO stability, the convolution integral is utilized when a bounded continuous-time input is applied to a Linear Time-Invariant (LTI) system....
Nuclear Stability
To hold positively charged protons together...
RNA Stability
Trial and Error and Algorithm

