Optimal control theory (OCT) applied to driving molecular systems by means of femtosecond pulses is now a mature area, but many of its intricacies are as yet unexplored. As a numerical tool, the many variations on the basic method differ not only in computer efficiency but in the type of solutions obtained. In this paper we survey this diversity, focusing on the use of multiphoton IR laser excitation to control either (1) the state selectivity or (2) the photodissociation in a ID Morse potential. We compare two distinct algorithms, the Krotov method and the gradient method. The former method generates large changes in the field at each iteration, white the latter does not. As a result, the Krotov method virtually always leads to pulses that are very different from the initial guess, while with the gradient method this is not always the case. We then analyze the effect of changing the final time, T, and find that it also can have a profound effect on the nature of the optimal solutions. Finally, we compare the solutions obtained using two different projectors to describe the bond-breaking process : a coordinate projector and a projector over scattering states. Again we observe that the optimal pulses and the dynamics they generate are markedly different in the two cases. This ambiguity in the definition of the optimal pulses may be viewed as a shortcoming of the approach, or alternatively it may be viewed as giving the method extra flexibility.