Volume + Position optimization (fixed cell shape)

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

Locked
Message
Author
ralf_tonner1
Newbie
Newbie
Posts: 17
Joined: Mon Dec 21, 2020 2:19 pm

Volume + Position optimization (fixed cell shape)

#1 Post by ralf_tonner1 » Thu Nov 03, 2022 6:18 pm

Dear admin,
In the vasp wiki, I couldn't find ISIF tag which can be used to perform both position and volume optimization while keeping the cell shape fixed. After going through the source code, I updated 2 lines in 'src/main.F' file as following

Code: Select all

src/main.F
...
[line 3996 (VASP-6.3.2)]
IF (DYN%ISIF<5) FACT=10*DYN%POTIM*EVTOJ/AMTOKG * 1E-10_q -->  IF ((DYN%ISIF<5).OR.(DYN%ISIF==20)) FACT=10*DYN%POTIM*EVTOJ/AMTOKG * 1E-10_q
...
[line 4026 (VASP-6.3.2)]
IF (DYN%ISIF==7) THEN --> IF ((DYN%ISIF==7).OR.(DYN%ISIF==20)) THEN
...
And now I can use ISIF=20 for position-volume optimization while keeping the cell shape fixed.
This worked for me, at least in simple geometry optimization (not molecular dynamics). I was wondering if this can be the right way to go.

Thanks

jonathan_lahnsteiner2
Global Moderator
Global Moderator
Posts: 149
Joined: Fri Jul 01, 2022 2:17 pm

Re: Volume + Position optimization (fixed cell shape)

#2 Post by jonathan_lahnsteiner2 » Thu Nov 03, 2022 8:04 pm

Dear ralf_tonner1,

Thank you for the information.
We are aware that there is no ISIF flag with the functionality you mention.
If a user desires to update positions and box volume ( without altering it's shape) it is recommended to use the ICONST file.
As described on this page
wiki/index.php/ICONST
under the header Restrictions on the volume and/or shape of the simulation cell
you can use the ICONST file to get the functionality you are looking for.

All the best Jonathan

ralf_tonner1
Newbie
Newbie
Posts: 17
Joined: Mon Dec 21, 2020 2:19 pm

Re: Volume + Position optimization (fixed cell shape)

#3 Post by ralf_tonner1 » Fri Nov 04, 2022 1:30 pm

Dear Jonathan,
Thank you very much. But if I get it right, the ICONST file is read only for molecular dynamics (IBRION=0). For simple geometry relaxation (IBRION=1 or 2) this will not work. or Am I missing something?

Thanks

jonathan_lahnsteiner2
Global Moderator
Global Moderator
Posts: 149
Joined: Fri Jul 01, 2022 2:17 pm

Re: Volume + Position optimization (fixed cell shape)

#4 Post by jonathan_lahnsteiner2 » Mon Nov 07, 2022 7:29 am

Dear ralf_tonner1,

No you are completely right. I missed that you are looking for way
doing the optimization with BRION=1,2. But the feature you are looking for will be added in one
of the coming vasp releases.

Thank you for the information.
All the Best Jonathan

matthewkuner
Newbie
Newbie
Posts: 12
Joined: Wed Oct 12, 2022 8:17 pm

Re: Volume + Position optimization (fixed cell shape)

#5 Post by matthewkuner » Fri Jan 27, 2023 1:00 am

As a follow up: Adding in a setting for shape-controlled (with volume and atomic position DoF) would be hugely useful in my DFT work. I look forward to seeing this in the future!

matthewkuner
Newbie
Newbie
Posts: 12
Joined: Wed Oct 12, 2022 8:17 pm

Re: Volume + Position optimization (fixed cell shape)

#6 Post by matthewkuner » Wed Mar 08, 2023 6:59 pm

@ jonathan_lahnsteiner2, in addition to the upcoming release for shape-controlled relaxations that you mentioned above, are there any plans to introduce a relaxation setting that would fix only the angles between the lattice vectors while running DFT? I think this would be greatly useful for studying metastable structures

jonathan_lahnsteiner2
Global Moderator
Global Moderator
Posts: 149
Joined: Fri Jul 01, 2022 2:17 pm

Re: Volume + Position optimization (fixed cell shape)

#7 Post by jonathan_lahnsteiner2 » Thu Mar 09, 2023 2:05 pm

Dear matthewkuner,

Yes I agree, this would be a very useful feature. The implementations of the
cell optimization was already done. Currently we are still testing the implementation.
In future versions it will be possible to keep the angles between the lattice vectors constant.
For the usage of the new feature, the vasp wiki will be updated as soon as the feature is available.

All the best Jonathan

ralf_tonner1
Newbie
Newbie
Posts: 17
Joined: Mon Dec 21, 2020 2:19 pm

Re: Volume + Position optimization (fixed cell shape)

#8 Post by ralf_tonner1 » Tue Aug 08, 2023 6:19 pm

Dear Admin,
I wanted to inquire about the implementation of geometry optimization with volume and position optimization using the fixed cell shape. I see in the source code that this feature has already been implemented in VASP 6.4.2 with ISIF=8 option. This information is not updated in the VASP wiki yet. Could you kindly confirm or provide any insights on whether this option is available for use?

Thanks

jonathan_lahnsteiner2
Global Moderator
Global Moderator
Posts: 149
Joined: Fri Jul 01, 2022 2:17 pm

Re: Volume + Position optimization (fixed cell shape)

#9 Post by jonathan_lahnsteiner2 » Thu Aug 10, 2023 9:44 am

Dear ralf_tonner1,

Thank you for the information. As you state the feature is implemented and can be used.
The vasp wiki was updated accordingly.

All the best Jonathan

ralf_tonner1
Newbie
Newbie
Posts: 17
Joined: Mon Dec 21, 2020 2:19 pm

Re: Volume + Position optimization (fixed cell shape)

#10 Post by ralf_tonner1 » Thu Aug 17, 2023 7:57 am

To admin:
Nice. Quick comment.
1. I believe the ISIF=8 feature is only valid for VASP>=6.4.2. If so, probably this information can be added to the wiki. I think, using vasp 5.4.4 and 6.3.x are still quite common.
2. Minor: "ISIF=0|1|2|3|4|5|6|7|8" <= At the top of the ISIF wiki page.

Thanks

jonathan_lahnsteiner2
Global Moderator
Global Moderator
Posts: 149
Joined: Fri Jul 01, 2022 2:17 pm

Re: Volume + Position optimization (fixed cell shape)

#11 Post by jonathan_lahnsteiner2 » Mon Aug 21, 2023 1:25 pm

Dear ralf_tonner,

Thank you for the reminder. The vasp wiki is updated accordingly.

Best Jonathan

Locked