How to deal with fractional occupation ?

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
mazuju
Posts: 15
Joined: Wed Jul 14, 2010 2:05 am

How to deal with fractional occupation ?

Post by mazuju » Tue Nov 16, 2010 3:08 am

Hi,

In compound ABC, A atoms and B atoms are at the same position, and half occupied.
How to write the position information in input file?

Best,
Ma ZJ

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: How to deal with fractional occupation ?

Post by jzwanzig » Tue Nov 16, 2010 11:36 am

The simple but very approximate way is to use the Virtual Crystal Approximation: see input variable mixalch and related variables. I think this approach is only allowed in abinit with norm-conserving pseudopotentials, not with PAW. The more "realistic" approach but more time-consuming is to make a supercell of your structure and simply populate sites by hand with the different atoms in the right ratio.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

mazuju
Posts: 15
Joined: Wed Jul 14, 2010 2:05 am

Re: How to deal with fractional occupation ?

Post by mazuju » Tue Nov 16, 2010 1:40 pm

jzwanzig wrote:The simple but very approximate way is to use the Virtual Crystal Approximation: see input variable mixalch and related variables. I think this approach is only allowed in abinit with norm-conserving pseudopotentials, not with PAW. The more "realistic" approach but more time-consuming is to make a supercell of your structure and simply populate sites by hand with the different atoms in the right ratio.


Thank you for your advice.

Locked