No abstract

Abstract. We show a novel systematic way to construct conservative finite difference schemes for quasilinear first-order system of ordinary differential equations with conserved quantities. In particular, this includes both autonomous and non-autonomous dynamical systems with conserved quantities of arbitrary forms, such as time-dependent conserved quantities. Sufficient conditions to construct conservative schemes of arbitrary order are derived using the multiplier method. General formulas for first-order conservative schemes are constructed using divided difference calculus. New conservative schemes are found for various dynamical systems such as Euler's equation of rigid body rotation, Lotka-Volterra systems, the planar restricted three-body problem and the damped harmonic oscillator.

We present the multiplier method of constructing conservative finite difference schemes for ordinary and partial differential equations. Given a system of differential equations possessing conservation laws, our approach is based on discretizing conservation law multipliers and their associated density and flux functions. We show that the proposed discretization is consistent for any order of accuracy when the discrete multiplier has a multiplicative inverse. Moreover, we show that by construction, discrete densities can be exactly conserved. In particular, the multiplier method does not require the system to possess a Hamiltonian or variational structure. Examples, including dissipative problems, are given to illustrate the method. In the case when the inverse of the discrete multiplier becomes singular, consistency of the method is also established for scalar ODEs provided the discrete multiplier and density are zero-compatible.

We show the arbitrarily long-term stability of conservative methods for autonomous ODEs. Given a system of autonomous ODEs with conserved quantities, if the preimage of the conserved quantities possesses a bounded locally finite neighborhood, then the global error of any conservative method with the uniformly bounded displacement property is bounded for all time, when the uniform time step is taken sufficiently small. On finite precision machines, the global error still remains bounded and independent of time until some arbitrarily large time determined by machine precision and tolerance. The main result is proved using elementary topological properties for discretized conserved quantities which are equicontinuous. In particular, long-term stability is also shown using an averaging identity when the discretized conserved quantities do not explicitly depend on time steps. Numerical results are presented to illustrate the long-term stability result.

We study the effect of an environment consisting of noninteracting two level systems on Landau-Zener transitions with an interest on the performance of an adiabatic quantum computer. We show that if the environment is initially at zero temperature, it does not affect the transition probability. An excited environment, however, will always increase the probability of making a transition out of the ground state. For the case of equal intermediate gaps, we find an analytical upper bound for the transition probability in the limit of large number of environmental spins. We show that such an environment will only suppress the probability of success for adiabatic quantum computation by at most a factor close to 1/2.