2003/12/01 by Siu A. Chin, Chin, Siu A.
Mathematics · Physics and Astronomy · #Computational Physics (physics.comp-ph) #FOS: Physical sciences #Mathematical Physics (math-ph) #math-ph #math.MP #physics.comp-ph
paper · pdf · doi:10.48550/arxiv.physics/0312005
10 pages, no figure. Completely rewritten and expanded
arxiv created 2004/08/03 · arxiv updated 2009/12/01
The solution of many physical evolution equations can be expressed as an exponential of two or more operators acting on initial data. Accurate solutions can be systematically derived by decomposing the exponential in a product form. For time-reversible equations, such as the Hamilton or the Schrödinger equation, it is immaterial whether or not the decomposition coefficients are positive. In fact, most symplectic algorithms for solving classical dynamics contain some negative coefficients. For time-irreversible systems, such as the Fokker-Planck equation or the quantum statistical propagator, only positive-coefficient decompositions, which respect the time-irreversibility of the diffusion kernel, can yield practical algorithms. These positive time steps only, forward decompositions, are a highly effective class of factorization algorithms. This work introduce a framework for understanding the structure of these algorithms. By a suitable representation of the factorization coefficients, we show that specific error terms and order conditions can be solved \it analytically. Using this framework, we can go beyond the Sheng-Suzuki theorem and derive a lower bound for the error coefficient eVTV. By generalizing the framework perturbatively, we can further prove that it is not possible to have a sixth order forward algorithm by including only the commutator [VTV]≡[V,[T,V]]. The pattern of these higher order forward algorithms is that in going from the (2n)\rm th to the (2n+2)\rm th order, one must include a new commutator [VT2n-1V] in the decomposition process.