This thread is locked.Only browsing is available.
Top Page > Browsing
NEB restart
Date: 2022/11/10 00:34
Name: Mehdi Vejdanihemmat   <m.vejdanih@gmail.com>

Dear All,

when I want to restart a stopped NEB run, I face below error like this person: https://www.openmx-square.org/forum/patio.cgi?mode=view&no=1887


"Format error for NEB1.Atoms.SpeciesAndCoordinates"

As I know there is no error with the format mentioned but I don't know why that happens. thank you in advance if you help me.

regards,
Mehdi
メンテ
Page: [1]

Re: NEB restart ( No.1 )
Date: 2022/11/11 11:29
Name: T. Ozaki

Hi,

Can I assume that you restarted with the dat# file generated with the first calculation?
Since I was not able to reproduce the issue on my trial, could you share an input file
for a minimal case which reproduces the issue?

Thank you for your cooperation in advance.

Regards,

TO
メンテ
Re: NEB restart ( No.2 )
Date: 2022/11/12 22:24
Name: Mehdi Vejdanihemmat  <m.vejdanih@gmail.com>

Dear Prof. Ozaki,

Thank you for your reply. I removed Hubbard related keywords and "off" option of each atomic coordinate (I used DFT instead of DFT+U) in dat#, then run that. The problem was solved (I mean the restart file (dat#) was run without the mentioned error). I think the related source code must be corrected and consider DFT+U method.

Best regards,
Mehdi Vejdanihemmat
メンテ
Re: NEB restart ( No.3 )
Date: 2022/12/01 17:50
Name: T. Ozaki

Hi,

I could reproduce your issue in case of the plus U method.
This is a bug in fact.
At this moment, could you eliminate the "on" or "off" flat in the NEB1.Atoms.SpeciesAndCoordinates
such as


<NEB1.Atoms.SpeciesAndCoordinates
  1    C  -0.66901031202659    0.00000000173238  -2.04470751948644    2.0    2.0  off
  2    C    0.66883867997214  -0.00000000248980  -2.04470705958138    2.0    2.0  off
  3    H    1.24206158038740  -0.92699468281945  -2.08949596495434    0.5    0.5  off
  4    H    1.24206156687189    0.92699467625303  -2.08949596488917    0.5    0.5  off
  5    H  -1.24218856731416  -0.92701843063594  -2.08950378204356    0.5    0.5  off
  6    H  -1.24218856604257    0.92701843424971  -2.08950377898137    0.5    0.5  off
  7    C  -0.66901032130227    0.00000000426894    2.04470750393403    2.0    2.0  off
  8    C    0.66883865524526    0.00000000170855    2.04470706094402    2.0    2.0  off
  9    H    1.24206157087749  -0.92699464416865    2.08949596852624    0.5    0.5  off
  10    H    1.24206157402616    0.92699464522607    2.08949596828241    0.5    0.5  off
  11    H  -1.24218855804748  -0.92701841093133    2.08950377953822    0.5    0.5  off
  12    H  -1.24218857663616    0.92701840471968    2.08950378821700    0.5    0.5  off
NEB1.Atoms.SpeciesAndCoordinates>

into

<NEB1.Atoms.SpeciesAndCoordinates
  1    C  -0.66901031202659    0.00000000173238  -2.04470751948644    2.0    2.0
  2    C    0.66883867997214  -0.00000000248980  -2.04470705958138    2.0    2.0
  3    H    1.24206158038740  -0.92699468281945  -2.08949596495434    0.5    0.5
  4    H    1.24206156687189    0.92699467625303  -2.08949596488917    0.5    0.5
  5    H  -1.24218856731416  -0.92701843063594  -2.08950378204356    0.5    0.5
  6    H  -1.24218856604257    0.92701843424971  -2.08950377898137    0.5    0.5
  7    C  -0.66901032130227    0.00000000426894    2.04470750393403    2.0    2.0
  8    C    0.66883865524526    0.00000000170855    2.04470706094402    2.0    2.0
  9    H    1.24206157087749  -0.92699464416865    2.08949596852624    0.5    0.5
  10    H    1.24206157402616    0.92699464522607    2.08949596828241    0.5    0.5
  11    H  -1.24218855804748  -0.92701841093133    2.08950377953822    0.5    0.5
  12    H  -1.24218857663616    0.92701840471968    2.08950378821700    0.5    0.5
NEB1.Atoms.SpeciesAndCoordinates>


The eliminatio of "on" or "off" should be done for all the keywords:
NEB1.Atoms.SpeciesAndCoordinates, NEB2.Atoms.SpeciesAndCoordinates, ...,
Then, you can restart with the dat# file normally even for the plus U calculation.
This will be fixed in the next release.

Thank you very for your report and coorporation in advance.

Regards,

TO

メンテ

Page: [1]