| 查看: 2810 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
VASP2BOLTZTRAP求助
|
||
|
应用 vasp2boltztrap出现错误 /usr/local/lib/python2.7/dist-packages/ase/lattice/spacegroup.py:5: UserWarning: Moved to ase.spacegroup warnings.warn('Moved to ase.spacegroup') Number of electrons not found. Please set the number manually in hte.intrans. get_vasp_bandstructure(): could not find OUTCAR file <./OUTCAR> Traceback (most recent call last): File "si.py", line 39, in <module> vasp2boltz.write_bandstructure_boltztrap(bs) File "/home/weiyidan/computing/电导率/vasp2boltz.py", line 100, in write_bandstructure_boltztrap if bandstructure['nspin']!=1: TypeError: 'NoneType' object has no attribute '__getitem__' poscar如下 NONDEFECT 1.00000000000000 7.6802000999000004 0.0000000000000000 0.0000000000000000 0.0000000000000000 10.8613996505999992 0.0000000000000000 0.0000000000000000 0.0000000000000000 7.6802000999000004 Si 32 Direct 0.0000000000000000 0.1250000000000000 0.5000000000000000 0.7500000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.7500000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.1250000000000000 0.0000000000000000 0.2500000000000000 0.3750000000000000 0.7500000000000000 0.5000000000000000 0.1250000000000000 0.0000000000000000 0.2500000000000000 0.3750000000000000 0.2500000000000000 0.2500000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.2500000000000000 0.2500000000000000 0.5000000000000000 0.2500000000000000 0.7500000000000000 0.7500000000000000 0.0000000000000000 0.5000000000000000 0.7500000000000000 0.3750000000000000 0.7500000000000000 0.5000000000000000 0.1250000000000000 0.5000000000000000 0.7500000000000000 0.3750000000000000 0.2500000000000000 0.0000000000000000 0.6250000000000000 0.5000000000000000 0.7500000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.7500000000000000 0.7500000000000000 0.0000000000000000 0.7500000000000000 0.2500000000000000 0.2500000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.6250000000000000 0.0000000000000000 0.2500000000000000 0.8750000000000000 0.7500000000000000 0.5000000000000000 0.6250000000000000 0.0000000000000000 0.2500000000000000 0.8750000000000000 0.2500000000000000 0.2500000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.7500000000000000 0.2500000000000000 0.5000000000000000 0.7500000000000000 0.7500000000000000 0.7500000000000000 0.5000000000000000 0.5000000000000000 0.7500000000000000 0.8750000000000000 0.7500000000000000 0.5000000000000000 0.6250000000000000 0.5000000000000000 0.7500000000000000 0.8750000000000000 0.2500000000000000 |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有6人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有4人回复
-
今年春晚有几个节目很不错，点赞！
已经有6人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Vasp+Boltztrap计算自旋极化材料求助
已经有7人回复
- 求助boltztrap里的vasp2boltz.py怎么使用？
已经有12人回复
5楼2016-12-27 13:09:32
|
我的 py # Example of how to use vasp2boltz.py # Create a .py file with a content similar to the following # It is necessary to have ASE installed: # https://wiki.fysik.dtu.dk/ase/ from ase import io from ase.lattice.spacegroup import Spacegroup import vasp2boltz ao = io.read('POSCAR') # If you do not have spglib installed, you need to add the spacegroup no. # e.g. sg = Spacegroup(216) for the 'F -4 3 m' space group # If you use the conventional cell rather than the primitive unit cell in your # POSCAR file for e.g. I or F spacegroups, you need to add this information. # Example 1: sg = Spacegroup(10) ao.info = {'spacegroup': sg} # Example 2. Here the unit cell in the POSCAR file is cubic (conventional): #sg = Spacegroup(10) #ao.info = {'spacegroup': sg, 'unit_cell': 'conventional'} # Read number of electrons from OUTCAR. foundnelect = False try: for line in open('OUTCAR', 'r'): if 'NELECT' in line: nelect = float(line.split()[2]) foundnelect = True except: pass if not foundnelect: print 'Number of electrons not found. Please set the number manually in hte.intrans.' nelect = 1.0 # The remaining part takes care of writing the files required by BoltzTraP. bs = vasp2boltz.get_vasp_bandstructure() vasp2boltz.write_bandstructure_boltztrap(bs) vasp2boltz.write_structure_boltztrap(ao) vasp2boltz.write_intrans_boltztrap(n_electrons = nelect) |
2楼2016-12-14 18:24:17
3楼2016-12-15 13:43:44
4楼2016-12-27 13:08:16