numat / RASPA2

NO LONGER UPDATED. Use the official repository.
https://github.com/iraspa/RASPA2
Other
77 stars 72 forks source link

Fixed bug in GCMC heat of adsorption calculation #30

Closed EfremBraun closed 6 years ago

EfremBraun commented 7 years ago

Using GCMC, RASPA calculates the heat of adsorption using formula 8 of doi.org/10.1021/ct700342k. In line 3723 of statistics.c, it can be seen that the system temperature is multiplied by the constant ENERGY_TO_KELVIN when it is already in units of Kelvin. This gives a heat of adsorption that is too high by 0.2*temperature, which is 0.5 kJ/mol at room temperature.