## LCALCEPS, IVDW, TOTAL ENERGY, and BORN EFFECTIVE CHARGES

Problems running VASP, runtime errors, technical questions.
Please check whether the answer to your question is given in
the VASP online manual or has been discussed in this forum
previously!

Moderators: Moderator, Global Moderator

### LCALCEPS, IVDW, TOTAL ENERGY, and BORN EFFECTIVE CHARGES

Hello,

I am using VASP to calculate the BORN EFFECTIVE CHARGES, with:

LCALCEPS = .TRUE.
EFIELD_PEAD = 0.0075 0.0075 0.0075

When I combined it with vdW corrections, either with TS or DFT-D2/D3, the final total energy writen in the OUTCAR or OSZICAR is:
TOTAL_ENERGY = TOTAL_ENERGY(SCF) + 2 * E_vdW
istead of just:
TOTAL_ENERGY = TOTAL_ENERGY(SCF) + 1 * E_vdW

It adds the vdW energy twice!
my questions:

1) Is this a bug?
2) If yes, what about the born charges. Are they ok? or they affected from this bug too?

Thanks from advanced!
Best regards,
Ido
ido
Newbie

Posts: 9
Joined: Mon Oct 29, 2012 8:44 am

### Re: LCALCEPS, IVDW, TOTAL ENERGY, and BORN EFFECTIVE CHARGES

1) Will be checked.
2) Born effective charges (BEC) depend only on the electronic structure.
Because vdW corrections do not change the electronic structure
calculated BEC are valid (the combination of of LCALCEPS=.TRUE.
and vdW corrections is not reasonable and will be blocked).
admin
Administrator

Posts: 2907
Joined: Tue Aug 03, 2004 8:18 am

### Re: LCALCEPS, IVDW, TOTAL ENERGY, and BORN EFFECTIVE CHARGES

Many Thanks, admin!!!

Regarding point 2), when vdW corrections will be blocked in BEC calculation, notice that it calculates forces differences F(E)-F(0) with respect to electric field E. So if the structure is vdW oprtimized than the vdW forces of the optimized structure should be included in the total forces F(0) (which are usually very small, below the threshold), as far as I understood it, since F(0) are the corresponding forces of the optimized structure with the vdW corrections.

Best regards,
Ido
ido
Newbie

Posts: 9
Joined: Mon Oct 29, 2012 8:44 am

Return to Bugreports

### Who is online

Users browsing this forum: No registered users and 1 guest