  We introduce a lattice Boltzmann for simulating an immiscible binary fluid
mixture. Our collision rules are derived from a macroscopic thermodynamic
description of the fluid in a way motivated by the Cahn-Hilliard approach to
non-equilibrium dynamics. This ensures that a thermodynamically consistent
state is reached in equilibrium. The non-equilibrium dynamics is investigated
numerically and found to agree with simple analytic predictions in both the
one-phase and the two-phase region of the phase diagram.
