Dear all,
- You mentioned that " Axel…the latest version of lammps I am using…that might be the reason for error happened while patch applied… "
Yongbiao answer: So, what is the solution. Should I use a old stable version? Or something else I can do?
- You mentioned that "
“nvt+gcmc results should agrees with the result of npt and nvt”.
I have no idea what this statement means.
- How do you check whether a gcmc simulation has converged ?
- Do you check the chemical potential ?
- What do you want out of gcmc simulation ?
These are some of the questions that you should start asking !."
Yongbiao answer:
By saying “nvt+gcmc results should agrees with the result of npt and nvt”, I mean,
If you perform nvt or npt simulation for a system, say, CO2 gas, you can obtain a density under a temperature T and a pressure P.
Then you turn to nvt+gcmc, the same density should be obtained under the same T and P. One difference is that in nvt+gcmc the number of molecules is fluctuated.
-
How do you check whether a gcmc simulation has converged ?
Answer: I check by compare the density, the pressure, the temperature to that of nvt simulation or npt simulation. For monoatomic systems, this check agrees very well. But in the systems of water or CO2, the temperature and pressure goes to very very high, the density goes to 0 or very very high. This can not be correct. -
Do you check the chemical potential ?
Answer: Yes, I checked. I also tried different values of potentials, just for tests. None of them works.
- What do you want out of gcmc simulation ?
Answer: I what to study some adsorption issues. But the first step must be making sure that gcmc works well. So I choose the simple systems such as Argon, CO2 and water to test the code. Because for these systems the results can be found in the liaterature. If I can not reproduce these results, it means my adsorption results can not be correct.
I have been learniing this method for long time. At the beginning I tested the SPCE water system, it does not work. So I thought maybe I should start from a evern simpler system, then I turned to the argon system. I found that nvt+gcmc works very well for argon systems. I thought I have known how to run gcmc simulation with lammps. Then I turned to CO2 for further tests. But it does not work. I spent much time on it but could not find why. I do not think the theory is difficult. Anyway, an example will help a lot. That is the reason why I am asking for an example. I think an example will tell me what my problem is.
I will be very grateful to your help!
Best,
Yongbiao