Page 1 of 1

BORN EFFECTIVE CHARGES

Posted: Tue Sep 29, 2020 9:44 am
by wwlphys
If I get a born effective charge:

BORN EFFECTIVE CHARGES (including local field effects)
------------------------------------------------------
ion 1
1 2.98061 0.03000 0.00000
2 -0.04000 2.98061 0.00000
3 0.00000 -0.00000 2.95148

when the displacement of the first ion is (0.01, 0.0, 0.0),is the dipole moment 0.0003 or -0.0004?

Re: BORN EFFECTIVE CHARGES

Posted: Tue Sep 29, 2020 2:28 pm
by henrique_miranda
You can find this information in the code: pead.F
Here is the point at which the Born effective charges are computed (when using PEAD):

Code: Select all

BORN_SCF(IDIR,:,:)=(TIFOR(:,:)-F_INIT(:,:))/EFIELD_PEAD(IDIR)
F_INIT(3,NIONS) and TIFOR(3,NIONS) contain the initial and final forces (after an electric field is applied).

Here is the line where they are printed:

Code: Select all

WRITE (IO%IU6,'(I5,3F12.5)') IDIR,BORN_SCF(IDIR,:,N)
so the first number on each line (1,2,3) corresponds to the (x,y,z) direction of the electric field.
To answer your question:
If you displace ion 1 along (1.0, 0.0, 0.0) you should get a dipole moment of (2.98061,-0.04000,0.0).

Re: BORN EFFECTIVE CHARGES

Posted: Wed Sep 30, 2020 11:35 am
by wwlphys
Dear Henrique,
Thank you!