  This paper (published in 2013) introduces a refinement of the equations for
the pair excitation function used in our previous work with D. Margetis. The
new equations are Euler-Lagrange equations, and the solutions conserve energy
and the number of particles.
