How to describe the effect of electric field on interlayer distcance

[yuc111]

Hi, all!

I would like to use VASP to calculate the interlayer distance with the effect of external electric field in bilayer BN. in this regard, I need to relax the lattice by using EFIELD or EFIELD_PEAD block in INCAR, but I get the different trend compared with experiments. So I am wondering is there any correct approaches, which can correctly describe the enfluence of the real electric field force? for example AIMD calculation under the elctric field?

Thank you for your time!

Regards,
Aster

[alexey.tal]

Dear Aster,

I moved your post to the "From users for users" section of our forum, since it is more of a question for the user community. Perhaps someone familiar with this problem can share their experience.
I don't have expertise in this topic, but it seems to me that it should be a relaxation calculation in the presence of an electric field, so I don't see why running AIMD would be a better way to relax the structure under the electric field. On the other hand, if the temperature effects play a major role for this system, AIMD would be a way to include the atomic vibrations.
If you provide more information, i.e., input and output files I can check if I see something that appears wrong in the way you set up your calculation.
Here are some things that one should be generally aware of for this type of calculations:

• Considering that it is a bilayer BN, one would need to include Van der Waals corrections to reproduce the experimental interlayer distance.
• One would likely need to include the dipole corrections for such a calculation.
• It is necessary to include a substantial amount of vacuum to avoid spurious interaction between periodic images.

Best wishes,
Alexey