In this paper the thermodiffusive behavior of an equimolar binary mixture subject to repulsive potentials of the form (/r) n is investigated by using nonequilibrium molecular dynamics (NEMD) and the thermodiffusion (Soret) coefficient, ST, is computed in a wide range of temperatures. With the aim to contribute to the study of the dependence of the Soret coefficient on the interaction potential, the exponent n of the potential is varied from 1 to 12, that is from a pseudocoulombian to a pseudohard-sphere interaction. The steady state equation is integrated for the composition function under reasonable assumptions and it is shown that in some cases the request for it to be linear cannot be satisfied. For this reason nonlinear functions are used to fit the NEMD composition data. The simulations indicate a negligible dependence of ST on the composition (in the composition range here considered) while the dependence on the temperature is more marked. The computed values of ST as a function of the temperature are fitted with analytical functions. It is found that with n3 (medium and short range interaction) the model system behaves like a dilute gas mixture with the Soret coefficient varying with the temperature almost like 1/T. In the case of n=1 (long range interaction), ST has a more complex dependence on T: in particular it shows a change of sign. The analytical fitting functions, ST (T), are used in the integrated steady state equation thus obtaining the steady state composition profile and its comparison with the NEMD results indicates the grounding of the approach here proposed. © 2008 American Institute of Physics.