iRASPA Community

RASPA => Output files => Topic started by: Daqiang on April 07, 2018, 03:08:37 PM

Title: ENERGY DRIFT
Post by: Daqiang on April 07, 2018, 03:08:37 PM
Hi Dr Dubbeldam:
   I used the RASPA2.0.3 to simulate the isotherm of acetylene. However, when the pressure is higher than about 1500 Pa, I got a warning in the output file.
   
Simulation finished,  1 warnings
WARNING: ENERGY DRIFT (INTERNAL CONSISTENCY ERROR IN THE CODE), THE SIMULATION RESULTS ARE WRONG!!

How can I resolve the warning? Adding the NumberOfCycles?

Thanks,

Daqiang
Title: Re: ENERGY DRIFT
Post by: David Dubbeldam on April 07, 2018, 04:28:20 PM
For any error, always retry with the latest version (which is RASPA 2.0.30).
Title: Re: ENERGY DRIFT
Post by: Daqiang on April 10, 2018, 08:50:56 AM
Thank David. Using the latest version, I have resolve the issue.
Title: Re: ENERGY DRIFT
Post by: Somananda Sanyal on March 04, 2019, 01:49:48 PM
Dear Dr. Dubbeldam,

    I am getting the same error with the latest version of RASPA, i.e. 2.0.35.
How to resolve the warning and if I may know too why is it occurring?

Thanks.
Somananda
Title: Re: ENERGY DRIFT
Post by: David Dubbeldam on March 06, 2019, 02:29:19 PM
This can possibly occur when not all potentials are well-defined. An example is when the hydrogens do NOT have a Lennard-Jones potential, but the connecting atom has a VDW radius that encompasses the hydrogen (or it should). But if the VDW radius is too small, then at high density the atoms will start to overlap leading to numerical problems.
Title: Re: ENERGY DRIFT
Post by: kokojiang on February 15, 2023, 08:10:07 AM
Dear Dr. Dubbeldam,

    I am getting the same error too. Like what you explained, the error happened as the number of adsorbates increased. When the number of adsorbates is low, it works well. But when the number of adsorbates is large, like 2500, the error appeared. In the forcefield I use, hydrogens did have no lennard-jones potential. How can I solve this warning? Should I change the forcefield or just run a larger simulation cycle?

Thank you!.
Jiang
Title: Re: ENERGY DRIFT
Post by: David Dubbeldam on June 13, 2023, 02:00:25 PM
I would modify the force field. One of the options is to put a hard-sphere potential around the hydrogen to shield it from overlap.