求助QE一般能够计算的原子总数

相信永远的雪

请问QE是否能够计算600多个原子的结构呢？我使用400个核，发现1小时1个电子步都没有动。

VASP可以很轻松计算这个case。

QE输入文件如下：

&CONTROL
calculation = 'relax'
nstep  = 1000 ! ion relax step
title = 'pwscf'
restart_mode = 'from_scratch'
outdir = './'
prefix = 'pwscf'
wf_collect = .false.
pseudo_dir = './'
pseudo_dir   = 'pseudo/ONCV',
etot_conv_thr = 1.0D-6
forc_conv_thr = 1.0D-6
tstress = .true.
tprnfor = .true.
/
&SYSTEM
ibrav = 0
nat=         656
ntyp=           6
ecutwfc = 50
ecutrho = 400
occupations = 'smearing'
degauss = 0.001
smearing = 'gaussian'
nosym =.TRUE. !no symmetry
input_dft='PBE'
vdw_corr='DFT-D3'
/
&ELECTRONS
electron_maxstep = 250
conv_thr = 1.0D-10
diagonalization = 'david' ! 'cg'
mixing_mode   = 'plain', ! default: 'plain'
/
  &IONS
ion_dynamics = 'bfgs'
!pot_extrapolation = 'second_order'
!wfc_extrapolation = 'second_order'
/
&CELL
cell_dynamics = 'bfgs'
!press_conv_thr = 0.1
/
ATOMIC_SPECIES
C  12.01   C_ONCV_PBE-1.2.upf
N  14.01   N_ONCV_PBE-1.2.upf
H  1.008   H_ONCV_PBE-1.2.upf
Pb 207.2   Pb_ONCV_PBE-1.2.upf
I  126.9   I_ONCV_PBE-1.2.upf
Cs 132.9   Cs_ONCV_PBE-1.2.upf

K_POINTS crystal
1
0 0 0 1 #band edge at gamma point

CELL_PARAMETERS (angstrom)
  24.7917575836    0   0
  0  24.7502574874    0
  0   0   25.1272659217

ATOMIC_POSITIONS (crystal)
  C    0.1250191180    0.1250365080    0.3901285530
  C    0.1250191180    0.1250365080    0.8901285530

...

请大家帮忙指出问题/意见，谢谢！

876449830

QE好像不太适合做优化。

丁越

400个核可真是太豪横了。把模型减小，这么多原子算除了cp2k能轻松肝动外，其他的程序就挺费劲了。

ChemG

QE优化确实比较慢，而且还没用对称性肯定更慢了，这么大的体系用黑科技CP2K多好。

leeru

话说600原子的话，100核的小集群也不是很快

leeru

我是用的cp2k，自己组的集群就100核

相信永远的雪

丁越 发表于 2022-7-19 09:22
400个核可真是太豪横了。把模型减小，这么多原子算除了cp2k能轻松肝动外，其他的程序就挺费劲了。

VASP还是比较快的。CP2K也可以。但是我需要用到一个特殊的方法（electron occupation setting），这个方法h好像只能在QE上用。

gog

相信永远的雪 发表于 2022-7-23 07:50
VASP还是比较快的。CP2K也可以。但是我需要用到一个特殊的方法（electron occupation setting），这个方 ...

没法建议。体系这么大，还不用对称性。

相信永远的雪

ChemG 发表于 2022-7-19 09:26
QE优化确实比较慢，而且还没用对称性肯定更慢了，这么大的体系用黑科技CP2K多好。

请问可以用VASP优化结构，然后把优化后的结构用QE做静态计算吗？
因为我需要设置electron occupation，VASP没有这个功能，QE有。

wth1219

亲测512个原子的体系、80 Ry的CUTOFF，72核的单点计算还是算得动的。

七尺贱

400核？！！！这么多核并行

卡开发发

相信永远的雪 发表于 2022-7-23 11:39
请问可以用VASP优化结构，然后把优化后的结构用QE做静态计算吗？
因为我需要设置electron occupation，V ...

你要electron occupation的什么选项？如果是展宽设置其实其他程序大体上都差不多，如果要约束占据数，VASP有FERWE和FERDO的字段。

wth1219

七尺贱 发表于 2022-7-23 13:03
400核？！！！这么多核并行

非常怀疑就是因为并行的核太多了才跑不动的……

七尺贱

wth1219 发表于 2022-7-23 13:25
非常怀疑就是因为并行的核太多了才跑不动的……

好奢侈，没有这么多计算资源可以用，最多也就72核并行

相信永远的雪

卡开发发 发表于 2022-7-23 13:03
你要electron occupation的什么选项？如果是展宽设置其实其他程序大体上都差不多，如果要约束占据数，VAS ...

谢谢！