Energetic pulses in exciton-phonon molecular chains and conservative numerical methods for quasilinear Hamiltonian systems

Abstract

The phenomenon of coherent energetic pulse propagation in exciton-phonon molecular chains such as α-helix protein is studied using an ODE system model of Davydov-Scott type, both with numerical studies using a new unconditionally stable fourth-order accurate energy-momentum conserving time discretization and with analytical explanation of the main numerical observations. Impulsive initial data associated with initial excitation of a single amide-I vibration by the energy released by ATP hydrolysis are used as well as the best current estimates of physical parameter values. In contrast to previous studies based on a proposed long-wave approximation by the nonlinear Schrödinger (NLS) equation and focusing on initial data resembling the soliton solutions of that equation, the results here instead lead to approximation by the third derivative nonlinear Schrödinger equation, giving a far better fit to observed behavior. A good part of the behavior is indeed explained well by the linear part of that equation, the Airy PDE, while other significant features do not fit any PDE approximation but are instead explained well by a linearized analysis of the ODE system. A convenient method is described for construction of the highly stable, accurate conservative time discretizations used, with proof of its desirable properties for a large class of Hamiltonian systems, including a variety of molecular models.