generate a ML-FF for CO2 molecules

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
wuyang_lin
Newbie
Newbie
Posts: 2
Joined: Wed Mar 29, 2023 4:29 pm

generate a ML-FF for CO2 molecules

#1 Post by wuyang_lin » Wed Mar 29, 2023 4:33 pm

Dear VASP developers, 

 

We are trying to generate a new machine learning force field for CO2 adsorption on copper surfaces. We have finished the training for bulk Cu and several surfaces and the results compare well with ab initio MD and DFT optimizations. However, once we take CO2 molecules into account, problems arise: The ML-FF results show that an isolated CO2 molecules can not keep linear and the angle within the molecule is always between 165° and 180°. Besides, we tried to add a constraint condition in ICONST file to fix the angle to 180°, which will cause calculation fail at once and RATTLE algorithm convergence isn’t achieved. Last, if we use a DFT-optimized configuration with CO2 adsorbed on copper surface to train the force field, the results are also inaccurate and will make the force field less accurate for pure copper. 

 

We have also tried to simulate liquid CO2 but we again observe bending of molecules. 

 

Could you please help us find out what could be the reasons for getting bent CO2?

Many thanks,
Wuyang

alexey.tal
Global Moderator
Global Moderator
Posts: 228
Joined: Mon Sep 13, 2021 12:45 pm

Re: generate a ML-FF for CO2 molecules

#2 Post by alexey.tal » Thu Mar 30, 2023 8:05 am

Dear Wuyang,

Could you please provide all the relevant input and output files for your calculation according to the forum guidelines.

wuyang_lin
Newbie
Newbie
Posts: 2
Joined: Wed Mar 29, 2023 4:29 pm

Re: generate a ML-FF for CO2 molecules

#3 Post by wuyang_lin » Thu Mar 30, 2023 3:23 pm

Dear Alexey,

Thank you very much for you reply. Please find attached the relevant files.

Wuyang
You do not have the required permissions to view the files attached to this post.

alexey.tal
Global Moderator
Global Moderator
Posts: 228
Joined: Mon Sep 13, 2021 12:45 pm

Re: generate a ML-FF for CO2 molecules

#4 Post by alexey.tal » Mon Apr 03, 2023 9:59 am

In the provided files, the temperature grows from 300 K up to 1800 K, which indicates that the thermostat was not set correctly.
Indeed, in the INCAR files, the parameters for a Langevin thermostat are provided, however, MDALGO = 2 is set, which corresponds to the Nose-Hoover thermostat.
If you would like to perform an NVT simulation with the Nose-Hoover thermostat, you should set the SMASS tag.
Please see the wiki MDALGO and SMASS

Post Reply