Register to reply 
How to calculate pseudopotentials?by jd1828af
Tags: pseudopotentials 
Share this thread: 
#1
Oct1305, 08:04 AM

P: 1

Can anyone explain how pseudopotentials are calculated using density functional theory?



#2
Oct1905, 09:25 AM

P: 10

In the book of
Planewaves, pseudopotentials and the LAPWmethod, D. Singh, Kluwer Academic Publishing (1994), ISBN 0792394217 it gives some introduction. 


#3
Nov1405, 09:32 AM

P: 4,006

regards marlon 


#4
Nov1405, 09:35 AM

P: 4,006

How to calculate pseudopotentials?
as an addendum, i especially recommend the reference nr 5 (Goedecker, phys rev B, vol 42, 8858)in the above paper for a more introductory explanation
regards and enjoy marlon 


#5
Nov1605, 08:36 AM

P: 4,006

One of the first REAL ab initio methods (ie no fitting with experimental data is required what so ever)to calculate pseudopotentials is : phys rev , vol 116, 287 (1959)
marlon 


#6
Nov2805, 05:46 PM

P: 44

does anyone know of where I could find some calculated pseudopotentials for a solid. Ill go for just about any solid, although CO2 is what I would really need. Im just trying to see if my calculations are anywhere close to what they should be.



#7
Dec205, 01:28 PM

P: 4,006

I do not understand how you can check your calculations without having the pseudo potential file to model the [valence electron] <> [nucleus + core electrons]interactions ??? How exactly do you procede ? I mean, assuming you are doing Hartree Fock or DFT calculations : 1) what software are you using ? Usually, the pseudopotential files (like the ) can be found on the server of the software developer ; like in the case of SIESTA or ABINIT. 2) What basis set are you using (like doubble zeta gaussians) or do you procede with plane waves ? 3) In DFT case, what approximation do you use for the exchange/correlation functional ? LDA, GGA, Hybrid models ??? For metals, you must use the LYPcorrelation functional of Parr et al (1988) together with the PBE exchange functional of Perdew, Burke and Ernzerhof (1996) 4) Normally, here you would chose your pseudo potential file. If you take one from literature, be sure that it is "transferable". Also, remember that the pseudopotential file that you choose will be partially determined by your exchange/correlation functional. You need to do a lot of benchmarking to get the right one. The mentioned functionals can all be found in the previous references that i gave you. marlon 


#8
Dec1005, 09:05 AM

P: 4,006

Ok, this maybe a bit late but i forgot to refer you to the following site on Vanderbilt's pseudopotentials :
here regards marlon 


#9
Dec1705, 12:13 PM

P: 44

Yes, I am doing hartreefock calculations using abinit. I just need to show that the program is working correctly before I try using any of the other methods. I was getting very strange energies which I found was caused by the cell volume being a little too small.



#10
Feb106, 07:05 AM

P: 4,006

Just to be complete, let me again show you the general way to proceed after you have selected all necessary input data (pseudo potentials, exchange correlation functional, atomic lattice and atomic positions) 1) perform a convergence test with respect to the energy cutoff value 2) perform a convergence test with respect to the selected kmesh. The bigger the unitcell, the smaller the kpoint mesh (due to the inverse connection between Wigner Seitz unit vectors and (reciprocal) Brillouin unit vectors). 3) Select the mesh and E cutoff of the above two tests and perform an relaxation of the atomic positions. 4) perform a relaxation of the lattice with the optimized atomic positions as input for this test. 5) The outcome of test 4 will give you the best structure (both latticedimensions and atomic positions) and be sure that after test 4, the atomic forces are below 0.05 eV/angström !!! If you have done the above tests, your problems should be solved. If they are not, there is something wrong with the pseudopotential that you used. Perhaps it is not transferrable, have you checked ? Select another pseudopotential file from the Vanderbilt website. Those files are very reliable. Another option is an incorrect input for the lattice and atomic positions. To be sure of that, use the input from the "Web of Elements" website. Finally, if you can, try to compare your results with a SIESTA simulation (if you know how to work with SIESTA ofcourse). regards marlon 


#11
Mar308, 05:15 AM

P: 70

My system is rather large  performing many many separate calculations will certainly take a very long time. Basically, I cut corners by increasing kpts and energy cutoff value simultaneously until the quantities that I wanted to obtain had converged. 


Register to reply 
Related Discussions  
How to calculate.  Calculus & Beyond Homework  4  
How do I calculate this?  Introductory Physics Homework  1  
Calculate e  Calculus  10  
How to Calculate the AU  Astronomy & Astrophysics  5 