VASP中怎样才能做出某个单原子的自旋密度图？

小兔

下图是我在文献中所看到的自旋密度图，文献作者只把特定的Ti上的自旋轨道图给表示出来了，而我是用vaspkit的32功能生成SPIN.vasp文件，在VESTA中打开，确实正个体系的都有，如图2，想问问大家真个应该怎么画出特定原子的自旋密度图？

wzkchem5

文献没有画单个原子的自旋密度图，只不过文献的自旋密度刚好有一部分局域在一个原子上，而你的自旋密度恰好不是这样的而已

Shana

Band decomposed charge densities.https://www.vasp.at/wiki/index.p ... ed_charge_densities
他画的局域轨道的电荷密度，一般是某一个轨道或者某个能量区间的电荷密度，你找一下电子或者空穴局域化的轨道，然后输出这一个轨道的电荷密度就可以了。如果你的轨道不局域的话，是不行的。

小兔

Shana 发表于 2022-5-15 11:55
Band decomposed charge densities.https://www.vasp.at/wiki/index.php/Band_decomposed_charge_densities ...

好的，谢谢老师，我再试一下

小兔

wzkchem5 发表于 2022-5-14 18:50
文献没有画单个原子的自旋密度图，只不过文献的自旋密度刚好有一部分局域在一个原子上，而你的自旋密度恰好 ...

嗯嗯，文献的确是这样的，但是老师，到底怎样才可以把自旋密度局域到一个原子上呢？我尝试文中说的，把特定原子的MGMGOM调大，然后增加了一个电子，但出来的图还是这样。下面是我的INCAR

小兔

小兔 发表于 2022-5-16 22:26
嗯嗯，文献的确是这样的，但是老师，到底怎样才可以把自旋密度局域到一个原子上呢？我尝试文中说的，把特 ...

C:\Users\Administrator\Desktop

小兔

小兔 发表于 2022-5-16 22:27

SYSTEM = ana
ISTART = 0
ICHARG = 2
ISPIN = 2
ISMEAR = 0
SIGMA = 0.5
GGA = PE
IVDW = 11
ENCUT = 400
EDIFF = 1E-4
NELM = 500
IBRION = 1
POTIM = 0.2
EDIFFG = 0.05
NSW = 1000
LREAL = .FALSE.
#PREC = Accurate
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 2 2 2 2
LDAUU = 0.0 4.0 2.0 2.0 2.0
LDAUJ = 0.0 0.0
LDAUPRINT = 0
LMAXMIX = 4
ALGO = FAST
EMIN = -20
EMAX = 20
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001
NELECT = 780
MAMGOM = 96*2 1*2 1*0 1*1 30*0

小兔

wzkchem5 发表于 2022-5-14 18:50
文献没有画单个原子的自旋密度图，只不过文献的自旋密度刚好有一部分局域在一个原子上，而你的自旋密度恰好 ...

希望老师指教

wzkchem5

小兔 发表于 2022-5-16 15:26
嗯嗯，文献的确是这样的，但是老师，到底怎样才可以把自旋密度局域到一个原子上呢？我尝试文中说的，把特 ...

你的体系和文献严格一致吗？
自旋密度是否局域到一个原子，是由体系本身的性质决定的，而不是你想局域就能局域的。如果这个体系的自旋本来就该是离域的，那你不管怎么算也不可能算出局域的自旋密度，而且也没必要非得算出局域的自旋密度。同理，文献算出来是局域的而你算出来是离域的，也是完全正常的，除非你的体系和计算方法都和实验严格一致，那样你和文献算得不一样才成为一个问题。

小兔

wzkchem5 发表于 2022-5-16 23:34
你的体系和文献严格一致吗？
自旋密度是否局域到一个原子，是由体系本身的性质决定的，而不是你想局域就 ...

是的是的老师，文献中计算的是TiO2上负载Pt，而我是想通过这种方法来计算WO3上负载Cu时，在光催化条件下，电子从体相到表面。在由表面到负载金属上的能量。

chen0201

vaspkit没有32功能啊？只有62，请问是62功能吗？另外请问是需要一个INPUT.in文件吗？

阿司匹林泡腾片

chen0201 发表于 2023-5-4 16:41
vaspkit没有32功能啊？只有62，请问是62功能吗？另外请问是需要一个INPUT.in文件吗？

vaspkit-312