24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 413  |  回复: 2
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 zxling 的 4 个金币

zxling

木虫 (小有名气)

[交流] 【求助】计算时出现“”段错误“是什么原因?

在用GW计算diamond时,在终端出现“段错误”是什么原因?
输入文件如下:
# Crystalline diamond
# Calculation of the GW correction to the direct band gap in Gamma
# Dataset 1: ground state calculation
# Dataset 2: calculation of the kss file
# Dataset 3: calculation of the screening (epsilon^-1 matrix for W)
# Dataset 4: calculation of the Self-Energy matrix elements (GW corrections)

ndtset      4

kptopt   1           
ngkpt    4 4 4      

# Dataset1: usual self-consistent ground-state calculation
# Definition of the k-point grid
nkpt1    10
nshiftk1  4
shiftk1  0.5 0.5 0.5  # This grid is the most economical
         0.5 0.0 0.0
         0.0 0.5 0.0
         0.0 0.0 0.5
prtden1  1         # Print out density


# Dataset2: calculation of kss file
# Definition of k-points
nkpt2    19             # A set of 19 k-points containing Gamma
nshiftk2  4
shiftk2  0.0 0.0 0.0  # This grid contains the Gamma point
         0.0 0.5 0.5
         0.5 0.0 0.5
         0.5 0.5 0.0
istwfk2  19*1                    # Option needed for Gamma
iscf2    -2             # Non self-consistent calculation
getden2  -1             # Read previous density file
nband2   9
nbandkss2 100        # Number of bands to store in KSS file

# Dataset3: Calculation of the screening (epsilon^-1 matrix)
optdriver3  3        # Screening calculation
getkss3     -1       # Obtain KSS file from previous dataset
nband3      25       # Bands to be used in the screening calculation
ecutwfn3    3.6      # Planewaves to be used to represent the wavefunctions
ecuteps3    6.0      # Dimension of the screening matrix
ppmfrq3    16.7 eV  # Imaginary frequency where to calculate the screening

# Dataset4: Calculation of the Self-Energy matrix elements (GW corrections)
optdriver4  4        # Self-Energy calculation
getkss4     -2       # Obtain KSS file from dataset 1
getscr4     -1       # Obtain SCR file from previous dataset
nband4      100      # Bands to be used in the Self-Energy calculation
ecutwfn4    5.0      # Planewaves to be used to represent the wavefunctions
ecutsigx4    6.0      # Dimension of the G sum in Sigma_x
                     # (the dimension in Sigma_c is controlled by npweps)
nkptgw4      1                # number of k-point where to calculate the GW correction
kptgw4                       # k-points
  0.000    0.000    0.000    # (Gamma)
bdgw4       4  5             # calculate GW corrections for bands from 4 to 5


# Definition of the unit cell: fcc
acell  3*6.722   # This is equivalent to  6.722
rprim  0.0  0.5  0.5   # FCC primitive vectors (to be scaled by acell)
       0.5  0.0  0.5   
       0.5  0.5  0.0

# Definition of the atom types
ntypat  1         # There is only one type of atom
znucl 6          # The keyword "znucl" refers to the atomic number of the
                  # possible type(s) of atom. The pseudopotential(s)
                  # mentioned in the "files" file must correspond
                  # to the type(s) of atom. Here, the only type is Silicon.
                        
# Definition of the atoms
natom 2           # There are two atoms
typat  1 1        # They both are of type 1, that is, Silicon.
xred              # Reduced coordinate of atoms
      0.0  0.0  0.0
      0.25 0.25 0.25

# Definition of the planewave basis set (at convergence 16 Rydberg 8 Hartree)
ecut 8.0          # Maximal kinetic energy cut-off, in Hartree

# Use only symmorphic operations
symmorphi 0

# Definition of the SCF procedure
nstep   10        # Maximal number of SCF cycles
diemac  12.0      # Although this is not mandatory, it is worth to
                  # precondition the SCF cycle. The model dielectric
                  # function used as the standard preconditioner
                  # is described in the "dielng" input variable section.
                  # Here, we follow the prescription for bulk silicon.
tolwfr  1.0d-10   

# This line added when defaults were changed (v5.3) to keep the previous, old behaviour
  iscf 5
当将晶格常数改变到7.3以上时,则在终端不会出现“段错误”,奇怪了,请问高手这个是什么原因?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mailoliver

金虫 (小有名气)

zxling(金币+5): 2010-03-30 09:05
zxling(金币+1): 2010-04-07 09:25
“段错误”与机器,程序和编译器都可能有关系。
如果,只是想用一下abinit就不要找原因了。你换一个编译器,或换一台机器试试。
2楼2010-03-25 07:45:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxling

木虫 (小有名气)

引用回帖:
Originally posted by mailoliver at 2010-03-25 07:45:23:
“段错误”与机器,程序和编译器都可能有关系。
如果,只是想用一下abinit就不要找原因了。你换一个编译器,或换一台机器试试。

谢谢!
3楼2010-03-30 09:05:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxling 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 372分,材料与化工,一志愿湖南大学,求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 土木304求调剂 +4 兔突突突, 2026-03-31 4/200 2026-04-04 13:34 by 1753564080
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 265求调剂 +11 yelck 2026-04-01 12/600 2026-04-01 19:12 by 549790059
[考研] 08工科275求调剂,可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
信息提示
请填处理意见