Request for Assistance with CINEB Calculations.
Posted: Sat Jun 10, 2023 2:20 pm
Dear Professor,
I want to use the CINEB method (VASP) to calculate the adsorption barrier of a thiophene molecule on a Cu (100) stepped surface. I have prepared the initial state where the molecule is located very far from the surface, approximately 9Å (gas phase), and the final state where the molecule is adsorbed in the most stable configuration (complex phase). Additionally, I have used the modemake.pl and nebmake.pl commands to generate intermediate states (5 images). The images were generated perfectly, with the distance between the surface and the molecule decreasing as we progress from the first to the last image.
However, I encountered a problem when I started the CINEB calculations to find the potential energy. The molecules dissociated spontaneously, even when I fixed the positions of the molecules. So I stopped the calculations. Or maybe the molecule's geometry will return to normal at the end of the relaxation?
I have attached a file to this message containing the scripts and the states related to these calculations.
I would be extremely grateful if you could review my question and provide me with a favorable response.
Thank you for your attention, and I look forward to your kind reply.
Yours sincerely,
I want to use the CINEB method (VASP) to calculate the adsorption barrier of a thiophene molecule on a Cu (100) stepped surface. I have prepared the initial state where the molecule is located very far from the surface, approximately 9Å (gas phase), and the final state where the molecule is adsorbed in the most stable configuration (complex phase). Additionally, I have used the modemake.pl and nebmake.pl commands to generate intermediate states (5 images). The images were generated perfectly, with the distance between the surface and the molecule decreasing as we progress from the first to the last image.
However, I encountered a problem when I started the CINEB calculations to find the potential energy. The molecules dissociated spontaneously, even when I fixed the positions of the molecules. So I stopped the calculations. Or maybe the molecule's geometry will return to normal at the end of the relaxation?
I have attached a file to this message containing the scripts and the states related to these calculations.
I would be extremely grateful if you could review my question and provide me with a favorable response.
Thank you for your attention, and I look forward to your kind reply.
Yours sincerely,