Electric field calculations

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
karthik_hj
Newbie
Newbie
Posts: 23
Joined: Wed Jul 20, 2022 2:24 pm

Electric field calculations

#1 Post by karthik_hj » Wed Jul 20, 2022 2:32 pm

I gave the electric field to an already relaxed structure and it shows 'electronic self consistency is not achieved'. I even tried NELM=300 but still the ground state energy has not been converged. If I put ISIF 2 tag along with electric field tag, calculations ends abruptly showing error. Is there anything else we can do. Would it be right if I take the CONTCAR after giving electric field and relax it again using ISIF tags?

Or is there another procedure to do the same

I put idipol ldipol and electric field tags in the INCAR . Should I put dipol tag and Efield_pead tag also

I want to find the effect of electric field on the electronic structure of a 2d metallic system

henrique_miranda
Global Moderator
Global Moderator
Posts: 414
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Electric field calculations

#2 Post by henrique_miranda » Wed Jul 20, 2022 10:05 pm

What is the error that you obtain when you say that the calculation ends abruptly?

Could you share the INCAR, POSCAR, KPOINTS, POTCAR and OUTCAR files?

karthik_hj
Newbie
Newbie
Posts: 23
Joined: Wed Jul 20, 2022 2:24 pm

Re: Electric field calculations

#3 Post by karthik_hj » Fri Jul 22, 2022 1:50 pm

The error is gone now when i put dipol tag also along with isif-2. Is it the correct procedure? What does an efield_pead tag do ? Could you provide a workflow for the electric field calculations for a structure?

henrique_miranda
Global Moderator
Global Moderator
Posts: 414
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Electric field calculations

#4 Post by henrique_miranda » Fri Jul 22, 2022 3:42 pm

The error is gone now when i put dipol tag also along with isif-2. Is it the correct procedure?
I don't know why setting ISIF would have any influence on the self-consistent computation with an electric field.
What does an efield_pead tag do ?
It applies an electric field.
wiki/index.php/EFIELD_PEAD
Could you provide a workflow for the electric field calculations for a structure
I don't understand what you mean

karthik_hj
Newbie
Newbie
Posts: 23
Joined: Wed Jul 20, 2022 2:24 pm

Re: Electric field calculations

#5 Post by karthik_hj » Wed Jul 27, 2022 5:51 am

I gave electric field using
LDIPOL = .TRUE.
EFIELD = 1.0
IDIPOL = 3
DIPOL = 0.5 0.5 0.5
Is it necessary to give EFIELD_PEAD tag? What does it do?

What i mean by workflow is the procedure to do electric field calculations. Please see whether which workflow is correct

Workflow 1
Give Electric field tags in INCAR and relax the structure with NSW=0
Give Electric field tags in INCAR and find dos and similarly bandstructure NSW=0


Workflow 2
Give Electric field tags in INCAR and relax the structure with NSW=0
Do not give Electric field tags in INCAR to find dos and bandstructure

henrique_miranda
Global Moderator
Global Moderator
Posts: 414
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Electric field calculations

#6 Post by henrique_miranda » Thu Jul 28, 2022 3:54 pm

The EFIELD tag refers to the electric field applied using the method described in https://doi.org/10.1103/PhysRevB.46.16067
This works for slabs or molecules as mentioned on this page:
wiki/index.php/Monopole_Dipole_and_Quad ... orrections

The EFIELD_PEAD tag refers to the electric field applied using the method described in https://doi.org/10.1103/PhysRevB.63.155107
This also works for bulk systems and its described on this page:
wiki/index.php/Berry_phases_and_finite_ ... ric_fields

You should not mix these two INCAR tags.

Workflow 1 is the correct one.

karthik_hj
Newbie
Newbie
Posts: 23
Joined: Wed Jul 20, 2022 2:24 pm

Re: Electric field calculations

#7 Post by karthik_hj » Mon Aug 01, 2022 5:03 am

Can these methods be only used for insulators/semiconductors? Because when I used EFIELD_PEAD for metals ,VASP exited.

henrique_miranda
Global Moderator
Global Moderator
Posts: 414
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Electric field calculations

#8 Post by henrique_miranda » Mon Aug 01, 2022 5:12 am

EFIELD_PEAD only works for systems with a gap:
https://www.vasp.at/wiki/index.php/EFIELD_PEAD

The maximum magnitude of the electric field you can apply is limited by the size of the gap of the system.

Post Reply