| 查看: 3245 | 回复: 3 | |||
xiehui0608金虫 (小有名气)
|
[求助]
求，vasp计算自旋密度作图方法，急！ 已有1人参与
|
| 本人是vasp初学者，最近在算一个带有自旋的体系，想做一下它的自旋密度图，在网上看了很多，看到有人说用CHGCAR 文件通过VESTA软件可以做成，但是没有具体的步骤方法，还是不会做，有哪位高手指点一下，急需要！谢谢！ |
回复此楼» 猜你喜欢
静等基金结果
已经有8人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有9人回复
-
好奇怪的filecode
已经有5人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：VASP 总态密度可以由两个自旋DOS绝对值相加不？
已经有10人回复
- vasp计算态密度结果出问题
已经有11人回复
- 求助vasp中的自旋轨道耦合
已经有3人回复
- vasp如何计算差分电荷密度
已经有7人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- VASP计算态密度出现问题
已经有12人回复
- 【求助】VASP计算体系性能没有f能带电子
已经有24人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【求助】vasp 离子驰豫与电子自旋的关系
已经有11人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】怎样用高斯计算自旋密度
已经有15人回复
Athena0218
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 296.1
- 散金: 710
- 红花: 8
- 帖子: 438
- 在线: 65.4小时
- 虫号: 5855744
- 注册: 2017-03-06
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★
xiehui0608(漫天飘雪代发): 金币+5, 谢谢交流 2018-06-13 15:08:45
xiehui0608(漫天飘雪代发): 金币+5, 谢谢交流 2018-06-13 15:08:45
|
在网上下载chgsplit.pl脚本，将CHGCAR生成CHGCAR_mag以及CHGCAR_tot.直接用VESTA打开CHGCAR_mag，这个是Alpha-Beta。调节VESTA工具栏objects的property里的isosurface。调节这个值显示的更清楚。懒得传附件了，chgsplit.pl脚本直接复制过来. #!/usr/bin/env perl #;-*- Perl -*- @ARGV>=1 || die "usage: chgsplit.pl <CHGCAR> \n"; open (IN1,$ARGV[0]) || die ("Can't open file $!"  ;open (OUT,">CHGCAR_tot" ;# read the header of the CHGCAR for ($i=0; $i<6; $i++) { $line1 = <IN1>; $header1 .= $line1; } $atoms1 = <IN1>; $header1 .= $atoms1; @atoms1 = split(/\s+/,$atoms1); $sum1 += $_ for @atoms1; print "Atoms in file: ".$sum1."\n"; for ($i=0; $i<$sum1+2; $i++) { $header1 .= <IN1>; } # start reading the total charge density (spin up plus spin down) part $points1 = <IN1>; $header1 .= $points1; @points1 = split(/\s+/,$points1); $psum1 = 1; for ($i=1; $i<@points1; $i++) { $psum1 *= $points1[$i]; } print "Points in total charge density: ".$psum1."\n"; print OUT $header1; for ($i=0; $i<$psum1/5; $i++) { $line1 = <IN1>; @line1 = split(/\s+/,$line1); printf OUT "%18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; } close(OUT); # start reading the magnetization density (spin up minus spin down) part open (OUT,">CHGCAR_mag" ;$line1 = <IN1>; while($line1 != $points1){ $line1 = <IN1>; } $points2 = $line1; @points2 = split(/\s+/,$points2); $psum2 = 1; for ($i=1; $i<@points2; $i++) { $psum2 *= $points2[$i]; } print "Points in magnetization density: ".$psum2."\n"; if ($psum1 != $psum2) {die ("Number of points not same in two parts!" ;}print OUT $header1; for ($i=0; $i<$psum1/5; $i++) { $line1 = <IN1>; @line1 = split(/\s+/,$line1); printf OUT "%18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; } close(OUT); close(IN1); . |
2楼2018-06-12 01:53:43
Athena0218
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 296.1
- 散金: 710
- 红花: 8
- 帖子: 438
- 在线: 65.4小时
- 虫号: 5855744
- 注册: 2017-03-06
- 专业: 理论和计算化学
【答案】应助回帖
| 内容已删除 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : chgsplit.pl
2018-06-12 01:55:10, 1.6 K
3楼2018-06-12 01:55:58
heyong0104
金虫 (正式写手)
硕士研究生
- 应助: 1 (幼儿园)
- 金币: 205.6
- 散金: 581
- 红花: 2
- 帖子: 473
- 在线: 72小时
- 虫号: 4468379
- 注册: 2016-03-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2018-06-13 15:06:19
100