  In this paper, we generalize the normalized gradient flow method which was
first applied to computing the least energy ground state to compute the least
action ground state. A continuous normalized gradient flow (CNGF) will be
presented and the action diminishing property will be proved to provide a
mathematical justification of the gradient flow with discrete normalization
(GFDN). Then we use backward-forward Euler method to further discretize the
GFDN in time which leads to the GFDN-BF scheme. It is shown that the GFDN-BF
scheme preserves the positivity and diminishes the action unconditionally. We
compare it with other three schemes which are modified from corresponding ones
designed for the least energy ground state and the numerical results show that
the GFDN-BF scheme performs much better than the others in accuracy, efficiency
and robustness for large time steps. Extensive numerical results of least
action ground states for several types of potentials are provided. We also use
our numerical results to verify some existing results and lead to some
conjectures.
