Для термодинамики -
http://webbook.nist.gov/chemistry/
Как расчитывать термодинамику - ну посмотри в википедии что-ли.
Для расчёта равновесных концентраций проще всего применять условную степень превращения, например:
2A + 3G = 4C + D
степень превращения x, x - от нуля (равновесие сдвинуто полностью влево) до 1 (равновесие сдвинуто полностью вправо)
Например в случае когда в начальный момент у нас A и G в стехиометрических соотношениях (допустим 0.2 и 0.3 M), имеем для концентрации каждого компонента при любой степени превращения
[A] = 0.2*(1-x)
[G] = 0.3*(1-x)
[C] = 0.4*x
[D] = 0.1*x
Если есть какие-то избытки или начальные концентрации продукта, надо учесть в качестве дополнительных слагаемых. Дальше дело техники. В равновесии
[C]^4*[D]/([A]^2*[G]^3) = K
Подставляем, получаем уравнение относительно одной переменной x, решаем, получаем равновесную степень превращения, расчитываем равновесные концентрации. Всё.
А в чём проблема с гетерофазной реакцией? Для термодинамических расчётов всё равно какая фаза, главное чтоб величины были для правильной фазы.