A lattice Boltzmann method for miscible gases is presented. In this model, the standard lattice Boltzmann method is employed for each species composing the mixture. Diffusion interaction among species is taken into account by means of a force derived from kinetic theory of gases. Transport coefficients expressions are recovered from the kinetic theory. Species with dissimilar molar masses are simulated by also introducing a force. Finally, mixing dynamics is recovered as shown in different applications: an equimolar counterdiffusion case, Loschmidt’s tube experiment, and an opposed jets flow simulation. Since collision is not altered, the present method can easily be introduced in any other lattice Boltzmann algorithms.