Adding user defined external potential

Documentation, Web site and code modifications

Moderators: baguetl, routerov

Locked
jbniser
Posts: 2
Joined: Wed Apr 25, 2018 1:28 pm

Adding user defined external potential

Post by jbniser » Wed Apr 25, 2018 2:36 pm

Not sure if it is worth a new topic.
I wish to introduce external potential of specific profiles to create inhomogeneity in charge density.
In quantum espresso I can do this easily by modifying the add_efield subroutine.
Would like to do the same in Abinit to be able to use some post-processing modules available with Abinit.
I could see that modifying vlocal can be a possibility but not sure at what level.
Kindly suggest possibilities.

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Adding user defined external potential

Post by ebousquet » Wed Apr 25, 2018 4:24 pm

Dear jbniser,
You could indeed add an additional potential potential in 67_common/setvtr.F90 and rhotov.F90. If you want to compare energy you'll also have to correct the energy accordingly in energy.F90.
All the best,
Eric

jbniser
Posts: 2
Joined: Wed Apr 25, 2018 1:28 pm

Re: Adding user defined external potential

Post by jbniser » Fri Apr 27, 2018 2:13 pm

Thanks Eric.

Looks like it is working.
I just added an external potential folded in the same way as vxc or vHartree, to vtrial wherever it was being calculated
in the three subroutines you had pointed out.

Locked