Dear Prof. Gale, I hope you are doing well.
I was testing GCMC simulations in GULP to estimate guest uptake in a porous material. Out of curiosity, I also tried the pGFN-FF model and was able to set up an initial input. However, when using the conv option, the MC moves are also applied to the framework atoms.
To keep the framework rigid, I froze all framework atoms using the three flags after their coordinates, but in that case the guests are often inserted too close to the material.
I have tested He and H2 as probe particles and several chemical potentials, trying to insert only a few guests, but the accepted configurations still show very short contacts. I am unsure whether this indicates a limitation of the model, a problem in how I am fixing the framework atoms, or something else.
So, my questions are: 1) is pGFN-FF intended to be used for GCMC simulations in GULP, and 2) is freezing the framework atoms with the coordinate flags the correct approach for rigid-framework GCMC simulations?
I have carefully reviewed the MC options in the manual, but I may be missing something here.
Best regards