The numerical simulations of bubble motion in dilute surfactant solutions reported previously by two of the authors contained a serious numerical inaccuracy. In agreement with experiments, single bubbles released from rest were predicted to reach a maximum speed before slowing to a terminal speed. However, subsequent experiments demonstrated that, in the simulations, the bubbles reached their terminal speed too quickly. The source of the discrepancy is an inaccuracy associated with the numerical algorithm used to solve the surfactant transport equation on the bubble surface. After correcting the problem, the simulations agree much better with the experiments. (C) 2004 Elsevier Inc. All rights reserved.