News:

SMF - Just Installed!

Main Menu
Menu

Show posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.

Show posts Menu

Topics - Ambroise

#1
Dear RASPA community,

I am simulating water adsorption in a rigid hydrophobic Zeolite. I am puzzling to grasp the relationship that links the external pressure "ExternalPressure" and the fugacity coefficent "FugacityCoefficient" in RASPA input script with the chemical potential calculated by RASPA to perform a GCMC simulation.

Instead of using the peng-Robinson state equation to calculate the fugacity, I input indeed myself the fugacity properly calculated before running a simulation for my water model. However, by using the widom insertion method to calculate the chemical potential of a box with only water (for tests), there is a mismatch with the expect chemical potential corresponding to the input fugacity coefficient.

For information, I used the relationship from statistical physics to compute the fugacity coefficient phi for input in RASPA:

phi = fugacity/P_ext = exp((mu_real(P_ext,T)-mu_perfect_gas(P_ext,T))/RT)

Is this relation used in RASPA or am I wrong?

Any help would be greatly appreciated!  :)