Correct way to force oxidation state?

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
cameronu
Newbie
Newbie
Posts: 3
Joined: Mon Mar 12, 2018 4:20 pm

Correct way to force oxidation state?

#1 Post by cameronu » Fri Aug 16, 2019 3:24 pm

Hi,

I am wishing to force the oxidation state of Sn in MASnI3 to be 4+ rather than the equilibrium 2+.

I have set NELECT = 48 rather than the expected 50 and explicitly set the Sn 5p occupancy to 0 in the POTCAR.

Is this reasonable or should I be setting the MAGMOM tag too?

Thanks in advance,

Cameron

cameronu
Newbie
Newbie
Posts: 3
Joined: Mon Mar 12, 2018 4:20 pm

Re: Correct way to force oxidation state?

#2 Post by cameronu » Wed Aug 21, 2019 4:41 pm

I have tested the charge per atom using bader charge analysis and noted that the charge for Sn has dropped from 13.146978 in the case of 2+ to 12.882354 in the case of 4+, so I am concerned this method did not work. If anybody knows of a good way to force the two electrons to be 'removed' from Sn please let me know.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Re: Correct way to force oxidation state?

#3 Post by admin » Thu Nov 07, 2019 8:28 am

The oxidation state of an atom is defined by surrounding atoms
and bonding to them. The change of the oxidation state can be achieved
with the change of the composition of the surrounding atomic sphere.
The variation of NELECT changes the number of electrons in the system
and does not lead to the local change on a specific atom.

Post Reply