NaF Berry Phase calculation

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Locked
Message
Author
tilak
Newbie
Newbie
Posts: 25
Joined: Thu Sep 09, 2010 1:40 pm

NaF Berry Phase calculation

#1 Post by tilak » Fri Jan 17, 2014 5:14 pm

Dear VASP users and experts,

In "Berry Phase Type calculation" on NaF, if I use IGPAR=1, I have got:
For undistorted case:
-----------------------
e<r>_ev=( 0.00000 0.00000 0.00000 ) e*Angst
e<r>_bp=( 0.00000 0.00000 0.00000 ) e*Angst

Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.00000 ) e*Angst
ionic dipole moment: p[ion]=( 2.25510 2.25510 2.25510 ) e*Angst
----------------------------------------

For distorted case:
----------------
e<r>_ev=( 0.00000 0.00000 0.00004 ) e*Angst
e<r>_bp=( 0.00000 0.18033 0.18033 ) e*Angst

Total electronic dipole moment: p[elc]=( 0.00000 0.18033 0.18037 ) e*Angst
ionic dipole moment: p[ion]=( 2.25510 2.25510 1.93939 ) e*Angst
--------------------------------------
From these above two set of data, I can have "Born Effective Charge (BEC)": 0.9961864219 elect., consistent with manual.


BUT, If use IGPAR=3, I have got:

For distorted case, (undistorted case as same before):
------------------
e<r>_ev=( 0.00000 0.00000 0.00004 ) e*Angst
e<r>_bp=( -0.18033 -0.18033 0.00000 ) e*Angst

Total electronic dipole moment: p[elc]=( -0.18033 -0.18032 0.00004 ) e*Angst
ionic dipole moment: p[ion]=( 2.25510 2.25510 1.93939 ) e*Angst
-----------------------------------

If you look into the Total Electronic dipole moment, it has BP contribution, with same mag. but different sign, and If you use it for BEC calculation, it will be: -14.9967481117 elect., not same than before! How do I explain?

Thank you very much for your comments,
Tilak
Last edited by tilak on Fri Jan 17, 2014 5:14 pm, edited 1 time in total.

support_vasp
Global Moderator
Global Moderator
Posts: 1817
Joined: Mon Nov 18, 2019 11:00 am

Re: NaF Berry Phase calculation

#2 Post by support_vasp » Thu Sep 12, 2024 9:03 am

Hi,

We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.

Best wishes,

VASP


Locked