songyu 发表于 2024-7-4 11:34

谢谢

乐平 发表于 2024-7-4 12:26
这个脚本在 VTST 官网上就能下载到

https://theory.cm.utexas.edu/vtsttools/download.html

好的，谢谢

单单 发表于 2024-7-1 14:19
请问dist.pl脚本可以发给我一个吗？谢谢

这个脚本在 VTST 官网上就能下载到

https://theory.cm.utexas.edu/vtsttools/download.html

vtstscripts.tgz   解压之后里面全是各种有用的脚本

单单 发表于 2024-7-1 14:19
请问dist.pl脚本可以发给我一个吗？谢谢

1. #!/usr/bin/env perl
   
2. #;-*- Perl -*-
   
3. 
   
4. use FindBin qw($Bin);
   
5. use lib "$Bin";
   
6. use Vasp;
   
7. 
   
8. # This program prints the difference in each coordinate for two POSCAR files.
   
9. 
   
10. # variables needed by the script
    
11. 
    
12. @args=@ARGV;
    
13. @args==2 || die "usage: dist.pl <POSCAR 1 filename> <POSCAR 2 filename> ";
    
14. 
    
15. $poscarfile1 = $args[0];
    
16. $poscarfile2 = $args[1];
    
17. 
    
18. ($coordinates1,$basis,$lattice,$num_atoms,$total_atoms,$selectiveflag,$selective)
    
19.   = read_poscar($poscarfile1);
    
20. 
    
21. ($coordinates2,$basis,$lattice,$num_atoms,$total_atoms,$selectiveflag,$selective)
    
22.   = read_poscar($poscarfile2);
    
23. 
    
24. $difference = pbc_difference($coordinates1,$coordinates2,$total_atoms);
    
25. $cartesian = dirkar($difference,$basis,$lattice,$total_atoms);
    
26. $distance = magnitude($cartesian,$total_atoms);
    
27. 
    
28. print $distance,"\n";
    

复制代码

请问dist.pl脚本可以发给我一个吗？谢谢

wypkdhd 发表于 2023-12-20 10:39
其实你还有两个参数有问题，既然是氧化物了。
If you have no a priori knowledge of your system, for in ...

非常感谢您的解答，参数这一部分我再仔细看看

其实你还有两个参数有问题，既然是氧化物了。
If you have no a priori knowledge of your system, for instance, if you do not know whether your system is an insulator, semiconductor or metal then always use Gaussian smearing ISMEAR=0 in combination with a small SIGMA=0.03-0.05.
插点数为距离/0.8四舍五入。
iopt对应算法，具体去vtst官网自己看，iopt=7没iopt=3安全。

wypkdhd 发表于 2023-12-20 08:33
iopt=3 插点数为5

是把这两个参数改成这样嘛？我尝试一下，非常感谢！顺便方便请教一下这样修改的理由嘛？

iopt=3 插点数为5