Physics – Computational Physics
Scientific paper
2003-12-01
Physics
Computational Physics
10 pages, no figure. Completely rewritten and expanded
Scientific paper
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\"odinger 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 $e_{VTV}$. 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]\equiv[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 $[VT^{2n-1}V]$ in the decomposition process.
No associations
LandOfFree
The Structure of Positive Decompositions of Exponential Operators does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with The Structure of Positive Decompositions of Exponential Operators, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and The Structure of Positive Decompositions of Exponential Operators will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-688210