24小时热门版块排行榜    

查看: 2460  |  回复: 16

cj4566

木虫 (正式写手)

[求助] 氧化锌吸附染料的pwscf计算

我在对氧化锌吸附体系(8*8*6超晶格吸附染料)进行relax计算的时候,输出文件卡在开头一直没变,大家看看是什么问题啊,谢谢了。
输出文件结尾:不知道charge为什么会是负的。
     number of k points=     1
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   2.0000000

     Dense  grid: 17475363 G-vectors     FFT dimensions: ( 250, 250, 640)

     Smooth grid:  4421473 G-vectors     FFT dimensions: ( 160, 160, 400)

     Largest allocated arrays     est. size (Mb)     dimensions
        Kohn-Sham Wavefunctions       706.11 Mb     (   8648, 5351)
        NL pseudopotentials          2054.32 Mb     (   8648,15568)
        Each V/rho on FFT grid          9.54 Mb     ( 625000)
        Each G-vector array             2.08 Mb     ( 273054)
        G-vector shells                 0.96 Mb     ( 125618)
     Largest temporary arrays     est. size (Mb)     dimensions
        Auxiliary wavefunctions      2824.43 Mb     (   8648,21404)
        Each subspace H/S matrix      279.52 Mb     (4280,4280)
        Each matrix   1271.12 Mb     (  15568, 5351)
        Arrays for rho mixing          76.29 Mb     ( 625000,   8)

     Initial potential from superposition of free atoms
     Check: negative starting charge=   -0.562288

     starting charge **********, renormalised to **********

     negative rho (up, down):  0.562E+00 0.000E+00
     Starting wfc are 6080 randomized atomic wfcs
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xueht987: 金币+1, 鼓励交流! 2012-07-14 19:22:55
是不是内存不够了呢?看样子你的任务没有几十G内存是不够的
2楼2012-07-13 09:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
xueht987: 金币+1, 鼓励交流! 2012-07-14 19:23:07
多少个原子?用了多少个核?几百个原子的话,没有100个cpu核就几乎算不动的。
3楼2012-07-13 10:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xueht987: 金币+1, 鼓励交流! 2012-07-14 19:23:24
不是CPU的问题,是内存远远不够。一般cluster都是一个CPU对应一定的内存,但是每个机器又有内存的上限。这相当于每个core需要接近8G的内存。你最好是试试256G内存的机器。同时也可以减小k点密度,减小ecutwfc,还有就是缩小cell的尺寸等等。
frank_zhan,zy2zhan@gmail.com
4楼2012-07-13 10:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)

最好能把你的.in文件也贴出来,这样大家才能知道确切的问题。
frank_zhan,zy2zhan@gmail.com
5楼2012-07-13 10:20:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

大家好热心啊,下面是输入文件:
&CONTROL
                 calculation  = 'relax',
                 restart_mode = 'from_scratch',
                 prefix       = "ZnO-886-218",
                 pseudo_dir   = "/pkg/suse11/quantum/5.0/pseudo",
                 outdir       = "/home/n7388683/QE/tmp",                 
/
&SYSTEM
                       ibrav = 4,
                   celldm(1) = 49.23871794442488,
                   celldm(2) = 1.201796131409272,
                   celldm(3) = 2.53237849247774,
                         nat = 1283,
                        ntyp = 6,
                     ecutwfc = 25,
                     ecutrho = 250,
/
&ELECTRONS
                 conv_thr    = 1.D-6,
                 mixing_beta = 0.6,
/
&IONS
                ion_dynamics = 'bfgs',           
/
ATOMIC_SPECIES
   Zn   65.3799972534  Zn.pbe-van.UPF
    O   15.9989995956  O.pbe-van_ak.UPF
    C   12.0109996796  C.pbe-van_bm.UPF
    H    1.0080000162  H.pbe-van_ak.UPF
    N   14.0069999695  N.pbe-van_ak.UPF
    S   32.0600013733  S.pbe-van_bm.UPF
ATOMIC_POSITIONS crystal
。。。
K_POINTS automatic
  1 1 1   0 0 0

一共1083个原子,太多了这里就省了,用的是64个cpu,480gb内存,下面是脚本:
#!/bin/bash -l

#PBS -N MPI
#PBS -l walltime=100:00:00
#PBS -l select=4:ncpus=16:mem=120gb

module load quantum/5.0.MPT
cd $PBS_O_WORKDIR

mpiexec -n 64 pw.x < ZnO-886-218.rx.in > ZnO-886-218.rx.out

是不是因为体系太大啊,谢谢了!!!
6楼2012-07-13 10:56:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
6楼: Originally posted by cj4566 at 2012-07-13 10:56:15
大家好热心啊,下面是输入文件:
&CONTROL
                 calculation  = 'relax',
                 restart_mode = 'from_scratch',
                 prefix       = "ZnO-886-218",
    ...

是1283个原子
7楼2012-07-13 10:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
3楼: Originally posted by goldenfisher at 2012-07-13 10:01:16
多少个原子?用了多少个核?几百个原子的话,没有100个cpu核就几乎算不动的。

按您这么说,我肯定得需要几百个cpu了
8楼2012-07-13 10:58:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


xueht987: 金币+1, 鼓励交流! 2012-07-14 19:23:57
引用回帖:
7楼: Originally posted by cj4566 at 2012-07-13 10:57:54
是1283个原子...

1k多个原子?多到逆天了……一般DFT最多也就见过一两百原子的,并且就已经要求是非常强悍的机器了……这1k多个……64个核心算肯定不够……
思想重于技巧,内涵重于表象
9楼2012-07-13 12:16:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
9楼: Originally posted by souledge at 2012-07-13 12:16:03
1k多个原子?多到逆天了……一般DFT最多也就见过一两百原子的,并且就已经要求是非常强悍的机器了……这1k多个……64个核心算肯定不够……...

嘿嘿  都逆天了阿,看来软件开发的空间还很大嘛。我看到下面文章中的体系也有上千原子啊:
http://pubs.rsc.org/en/content/articlelanding/2012/jm/c1jm13605d
10楼2012-07-13 12:23:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见