  A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815,
(1993)] lattice Boltzmann model for simulating fluids with multiple components
and interparticle forces is described in detail. Macroscopic equations
governing the motion of each component are derived by using Chapman-Enskog
method. The mutual diffusivity in a binary mixture is calculated analytically
and confirmed by numerical simulation. The diffusivity is generally a function
of the concentrations of the two components but independent of the fluid
velocity so that the diffusion is Galilean invariant. The analytically
calculated shear kinematic viscosity of this model is also confirmed
numerically.
