A lattice Boltzmann (LB) method is presented for solving the energy conservation equation in two phases when the phase change effects are included in the model. This approach employs multiple distribution functions, one for a pseudotemperature scalar variable and the rest for the various species. A nonideal equation of state (EOS) is introduced by using a pseudopotential LB model. The evolution equation for the pseudotemperature variable is constructed in such a manner that in the continuum limit one recovers the well known macroscopic energy conservation equation for the mixtures. Heats of reaction, the enthalpy change associated with the phase change, and the diffusive transport of enthalpy are all taken into account; but the dependence of enthalpy on pressure, which is usually a small effect in most nonisothermal flows encountered in chemical reaction systems, is ignored. The energy equation is coupled to the LB equations for species transport and pseudopotential interaction forces through the EOS by using the filtered local pseudotemperature field. The proposed scheme is validated against simple test problems for which analytical solutions can readily be obtained.