 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 3054 | 回復(fù): 3 | |||
xiehui0608金蟲 (小有名氣)
|
[求助]
求，vasp計(jì)算自旋密度作圖方法，急！ 已有1人參與
|
| 本人是vasp初學(xué)者，最近在算一個(gè)帶有自旋的體系，想做一下它的自旋密度圖，在網(wǎng)上看了很多，看到有人說用CHGCAR 文件通過VESTA軟件可以做成，但是沒有具體的步驟方法，還是不會(huì)做，有哪位高手指點(diǎn)一下，急需要！謝謝！ |
回復(fù)此樓» 猜你喜歡
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
-
一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！
已經(jīng)有3人回復(fù)
-
310求調(diào)劑
已經(jīng)有3人回復(fù)
-
機(jī)械專碩299求調(diào)劑至材料
已經(jīng)有4人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有12人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
301求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助：VASP 總態(tài)密度可以由兩個(gè)自旋DOS絕對(duì)值相加不？
已經(jīng)有10人回復(fù)
- vasp計(jì)算態(tài)密度結(jié)果出問題
已經(jīng)有11人回復(fù)
- 求助vasp中的自旋軌道耦合
已經(jīng)有3人回復(fù)
- vasp如何計(jì)算差分電荷密度
已經(jīng)有7人回復(fù)
- vasp自洽計(jì)算K點(diǎn)不匹配問題，急�。�
已經(jīng)有15人回復(fù)
- VASP計(jì)算態(tài)密度出現(xiàn)問題
已經(jīng)有12人回復(fù)
- 【求助】VASP計(jì)算體系性能沒有f能帶電子
已經(jīng)有24人回復(fù)
- 【求助】求VASP計(jì)算ZnO能帶和態(tài)密度的詳細(xì)輸入設(shè)置
已經(jīng)有15人回復(fù)
- 【求助】vasp 離子馳豫與電子自旋的關(guān)系
已經(jīng)有11人回復(fù)
- 【求助】VASP計(jì)算電子態(tài)密度
已經(jīng)有34人回復(fù)
- 【求助】怎樣用高斯計(jì)算自旋密度
已經(jīng)有15人回復(fù)
Athena0218
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 296.1
- 散金: 710
- 紅花: 8
- 帖子: 438
- 在線: 65.4小時(shí)
- 蟲號(hào): 5855744
- 注冊(cè): 2017-03-06
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
xiehui0608(漫天飄雪代發(fā)): 金幣+5, 謝謝交流 2018-06-13 15:08:45
xiehui0608(漫天飄雪代發(fā)): 金幣+5, 謝謝交流 2018-06-13 15:08:45
|
在網(wǎng)上下載chgsplit.pl腳本，將CHGCAR生成CHGCAR_mag以及CHGCAR_tot.直接用VESTA打開CHGCAR_mag，這個(gè)是Alpha-Beta。調(diào)節(jié)VESTA工具欄objects的property里的isosurface。調(diào)節(jié)這個(gè)值顯示的更清楚。懶得傳附件了，chgsplit.pl腳本直接復(fù)制過來(lái). #!/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
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 296.1
- 散金: 710
- 紅花: 8
- 帖子: 438
- 在線: 65.4小時(shí)
- 蟲號(hào): 5855744
- 注冊(cè): 2017-03-06
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
| 內(nèi)容已刪除 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : chgsplit.pl
2018-06-12 01:55:10, 1.6 K
3樓2018-06-12 01:55:58
heyong0104
金蟲 (正式寫手)
碩士研究生
- 應(yīng)助: 1 (幼兒園)
- 金幣: 205.6
- 散金: 581
- 紅花: 2
- 帖子: 473
- 在線: 72小時(shí)
- 蟲號(hào): 4468379
- 注冊(cè): 2016-03-05
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2018-06-13 15:06:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiehui0608 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|