应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求，vasp计算自旋密度作图方法，急！
41/1返回列表
查看: 3049  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiehui0608

金虫 (小有名气)

[求助] 求，vasp计算自旋密度作图方法，急！ 已有1人参与

本人是vasp初学者，最近在算一个带有自旋的体系，想做一下它的自旋密度图，在网上看了很多，看到有人说用CHGCAR 文件通过VESTA软件可以做成，但是没有具体的步骤方法，还是不会做，有哪位高手指点一下，急需要！谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-03 08:44:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Athena0218

新虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Athena0218

新虫 (正式写手)

【答案】应助回帖

内容已删除
赞一下(1人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : chgsplit.pl
    • 2018-06-12 01:55:10, 1.6 K
3楼2018-06-12 01:55:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heyong0104

金虫 (正式写手)

硕士研究生
★ ★
xiehui0608(漫天飘雪代发): 金币+2, 谢谢交流 2018-06-13 15:08:55
用王老师的vaspkit可以处理得到。

发自小木虫IOS客户端
赞一下(1人) 回复此楼
待人要好，做事要专心，少说话，多做事。
4楼2018-06-13 15:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiehui0608 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂，是学校就去 +6 sllhht 2026-03-19 7/350 2026-03-19 19:50 by 制度的
[考研] 288求调剂 +14 于海海海海 2026-03-19 14/700 2026-03-19 19:40 by maocaozhuxi
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 一志愿南昌大学，327分，材料与化工085600 +3 Ncdx123456 2026-03-19 3/150 2026-03-19 13:18 by houyaoxu
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭