W. Forner, ACCURACY OF DAVYDOV VERTICAL-BAR-D-1] APPROXIMATION FOR SOLITON DYNAMICS IN PROTEINS, Physical review. B, Condensed matter, 53(10), 1996, pp. 6291-6317
For the Davydov Hamiltonian several special cases are known which can
be solved analytically. Starting from these cases we show that the ini
tial state for a simulation using Davydov's \D-1] approximation has to
be constructed from a given set of initial lattice displacements and
momenta in the form of a coherent state with its amplitudes independen
t of the lattices site, corresponding to Davydov's \D-2] approximation
. The site dependences in the \D-1] ansatz evolve from this initial st
ate exclusively via the equations of motion. Starting the \D-1] simula
tion from an ansatz with site-dependent coherent-state amplitudes lead
s to an evolution which is different from the analytical solutions for
the special cases. Thus also in applications of the \D-1] ansatz to p
olyacetylene \D-2]-type initial states always have to be used in contr
ast to our previous suggestion [W. Forner, J. Phys. Condens. Matter 6,
9105 (1994)]. Further we expand the known exact solutions in Taylor s
eries in time and compare expectation values in different orders with
the exact results. We find that for an approximation up to third order
in time (for the wave function) norm, and total energy, as well as di
splacements and momenta are reasonably correct for a time up to approx
imate to 0.12-0.14 ps, depending somewhat on the coupling strength for
the transportless case. We performed long-time simulations using the
\D-1] approximation where we computed expectation values of the releva
nt operators with the state ((H) over cap/J)\D-1] and the deviation \d
elta] from the exact solution over long times, namely 10 ns. We found
that in the very long-time scale the \D-1] ansatz is very close to an
exact solution. Further we report results from an investigation of the
very short-time behavior of the \D-1] state compared with that of an
expansion of the exact solution in powers of time t. Within a time of
roughly 0.10-0.15 ps the second- and third-order corrections turned ou
t to be not very important. This is due to the fact that our first-ord
er state contains already some terms of the expansion, summed up to in
finite order. We found good agreement of the results obtained with our
expansion and those from the corresponding \D-1] simulations within t
he time of about 0.10 ps. Altogether we have shown that the \D-1] stat
e, although of approximative nature, is very close to an exact solutio
n of the Davydov model on time scales from some fs up to ns.