Adding user defined external potential

Moderators: jbeuken, jzwanzig, gonze, mcote

Adding user defined external potential

Postby 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.
jbniser
 
Posts: 2
Joined: Wed Apr 25, 2018 1:28 pm

Re: Adding user defined external potential

Postby 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
ebousquet
 
Posts: 158
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Adding user defined external potential

Postby 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.
jbniser
 
Posts: 2
Joined: Wed Apr 25, 2018 1:28 pm


Return to Proposed code modifications

Who is online

Users browsing this forum: No registered users and 1 guest

cron