Restrictions of Quadratic Forms to Lagrangian Planes,Quadratic Matrix Equations,and Gyroscopic Stabilization

We discuss the symplectic geometry of linear Hamiltonian systems with nondegenerate Hamiltonians. These systems can be reduced to linear second-order differential equations characteristic of linear oscillation theory. This reduction is related to the problem on the signatures of restrictions of quadratic forms to Lagrangian planes. We study vortex symplectic planes invariant with respect to linear Hamiltonian systems. These planes are determined by the solutions of quadratic matrix equations of a special form. New conditions for gyroscopic stabilization are found.     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Dirac constraint analysis and symplectic structure of anti-self-dual Yang-Mills equations

We present the explicit form of the symplectic structure of anti-self-dual Yang-Mills (ASDYM) equations in Yang’s J- and K-gauges in order to establish the bi-Hamiltonian structure of this completely integrable system. Dirac’s theory of constraints is applied to the degenerate Lagrangians that yield the ASDYM equations. The constraints are second class as in the case of all completely integrable systems which stands in sharp contrast to the situation in full Yang-Mills theory. We construct the Dirac brackets and the symplectic 2-forms for both J- and K-gauges. The covariant symplectic structure of ASDYM equations is obtained using the Witten-Zuckerman formalism. We show that the appropriate component of the Witten-Zuckerman closed and conserved 2-form vector density reduces to the symplectic 2-form obtained from Dirac’s theory. Finally, we present the Bäcklund transformation between the J- and K-gauges in order to apply Magri’s theorem to the respective two Hamiltonian structures.     

On Covariant Phase Space and the Variational Bicomplex

The notion of a phase space in classical mechanics is well known. The extension of this concept to field theory however, is a challenging endeavor, and over the years numerous proposals for such a generalization have appeared in the literature. In this paper We review a Hamiltonian formulation of Lagrangian field theory based on an extension to infinite dimensions of J.-M. Souriau's symplectic approach to mechanics. Following G. Zuckerman, we state our results in terms of the modern geometric theory of differential equations and the variational bicomplex. As an elementary example, we construct a phase space for the Monge–Ampere equation.     

Canonical measures on the moduli spaces of compact Riemann surfaces

We study some explicit relations between the canonical line bundle and the Hodge bundle over moduli spaces for low genus. This leads to a natural measure on the moduli space of every genus which is related to the Siegel symplectic metric on Siegel upper half-space as well as to the Hodge metric on the Hodge bundle.     

Application of symplectic algorithms to QCT calculation:H+H_2 system

The basic theory of symplectic algorithm was introduced. A comparison between Runge-Kutta method and symplectic integration method was preformed in the simulation of the long time behavior of H + H2 system on BKMP potential energy surface. Our results reveal a dis-sipative behavior in the integral of ordinary differential equation by the fourth order Runge-Kutta method, which causes incorrect simulation results in QCT calculations. However, when the symplectic integration method is applied, the dissipative behavior is not found in the same system. When the initial state is the same, the energy deviation of fourth order symplectic integral method is almost one percent of that of fourth order Runge-Kutta method in a 60000-step simulation, and that of sixth order symplectic integral method is much less. These results show that the symplectic integral methods are always the better choice in the integral calculation of the long time behavior in maintaining energy conservation.     

Exponentially fitted symplectic methods for the numerical integration of the Schrödinger equation

The computation of the energy eigenvalues of the one-dimensional time-independent Schrödinger equation is considered. Exponentially fitted and trigonometrically fitted symplectic integrators are obtained, by modification of the first and second order Yoshida symplectic methods. Numerical results are obtained for the one-dimensional harmonic oscillator and Morse potential.AMS subject classification: 65L15Funding by research project 71239 of Prefecture of Western Macedonia and the E.U. is gratefully acknowledged.