计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: yihanxu
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 如何解决Domain decomposition does not work with nstlist=0呢?

[复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

16#
发表于 Post on 2019-3-31 12:30:37 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-31 08:22
谢谢老师。如果是gmx2018.4且设置如下,buffer距离是0吗?
pbc = xyz
rlist = 1.2

一看手册就知道,你手动设这个没意义,程序是自动确定的
With the cutoff-scheme=Verlet cutoff-scheme, this is by default set by the verlet-buffer-tolerance option and the value of rlist is ignored.

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

17#
 楼主 Author| 发表于 Post on 2019-4-15 10:54:22 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-4-15 21:14 编辑

@sobereva 老师好,我想研究真空中分子间相互作用,参考了Actos四聚体分子自组装的MD模拟。体系包含43个S4(2-)阴离子、86个Li+阳离子、和1个聚合物链。但是提交md任务后会出现Fatal error: Domain decomposition does not work with nstlist=0,请问该怎么解决呢?在其他条件都不变的情况下,我又分别试了只放1个聚合物链、1个聚合物链+1组离子(1个S4+2个Li)、1个聚合物链+5组离子,都没有报错,但是1个聚合物链+10组离子就会发生相同的错误了。请老师指点。谢谢。

附md.mdp
define =
integrator = md
dt         = 0.002   ; ps
nsteps     = 3000000
comm-grps  = system
comm-mode  = Angular
energygrps =
;
nstxout = 0
nstvout = 0
nstfout = 0
nstlog  = 500
nstenergy = 1000
nstxout-compressed = 1000
compressed-x-grps  = system
;
pbc = no
nstlist = 0
rlist = 0
cutoff-scheme = group
coulombtype   = cut-off
rcoulomb      = 0
vdwtype       = cut-off
rvdw          = 0
;
Tcoupl  = V-rescale
tau_t   = 0.2
tc_grps = system
ref_t   = 298.15
;
gen_vel  = no
gen_temp = 298.15
gen_seed = -1
;
freezegrps  =
freezedim   =
constraints = hbonds





onePANI_multiS.png (478.48 KB, 下载次数 Times of downloads: 33)

初始结构

初始结构

onePANI_multiS - Copy.gro

17.73 KB, 下载次数 Times of downloads: 2

onePANI_5S.gro

7.49 KB, 下载次数 Times of downloads: 2

md.log

14.7 KB, 下载次数 Times of downloads: 2

live each day as if it was your last

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

18#
 楼主 Author| 发表于 Post on 2019-4-16 11:22:34 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-4-15 21:23 编辑

@sobereva 老师好,因为上一楼的error我没能解决,所以又参考了纳米水球MD模拟的设置来研究真空中分子间相互作用,但还是没成功,可以请老师帮忙看看可以怎么办吗?谢谢。

onePANI_multiS.gro

17.32 KB, 下载次数 Times of downloads: 2

md.log

14.97 KB, 下载次数 Times of downloads: 1

md.mdp

743 Bytes, 下载次数 Times of downloads: 2

md.sbatch

470 Bytes, 下载次数 Times of downloads: 1

live each day as if it was your last

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

19#
 楼主 Author| 发表于 Post on 2019-4-17 09:59:09 | 只看该作者 Only view this author
mdrun用了-nt 1,可以计算了
live each day as if it was your last

326

帖子

0

威望

5312

eV
积分
5638

Level 6 (一方通行)

20#
发表于 Post on 2020-4-18 20:48:57 | 只看该作者 Only view this author
yihanxu 发表于 2019-4-17 09:59
mdrun用了-nt 1,可以计算了

lz您好,我们碰到了类似的问题,但不能一直一个核跑吧,您后来找到更合适的解决办法了吗

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

21#
 楼主 Author| 发表于 Post on 2020-4-19 04:49:04 | 只看该作者 Only view this author
mol 发表于 2020-4-18 06:48
lz您好,我们碰到了类似的问题,但不能一直一个核跑吧,您后来找到更合适的解决办法了吗

很可惜 没找到其他好办法 有的情况试其他的核数 说不准几个核可以
live each day as if it was your last

191

帖子

0

威望

1098

eV
积分
1289

Level 4 (黑子)

22#
发表于 Post on 2020-12-6 21:30:29 | 只看该作者 Only view this author
sobereva 发表于 2019-3-23 03:50
你的这个链本身是有限长度的,还是无限长度的跨越盒子那种?如果是前者,em.gro里就不应该有盒子信息,而且 ...

请问老师链本身是有限长度的,em.gro里就不应该有盒子信息——是删掉gro文件最后一行解决吗?删掉后还是报Domain decomposition does not work with nstlist=0,请问如何解决呢,谢谢
科音学员

191

帖子

0

威望

1098

eV
积分
1289

Level 4 (黑子)

23#
发表于 Post on 2020-12-6 21:30:56 | 只看该作者 Only view this author
请问同学你最后是怎么解决的呢,谢谢
科音学员

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

24#
 楼主 Author| 发表于 Post on 2020-12-7 04:58:38 | 只看该作者 Only view this author
一条君 发表于 2020-12-6 07:30
请问同学你最后是怎么解决的呢,谢谢

mdrun用了-nt 1,可以计算了
live each day as if it was your last

53

帖子

0

威望

787

eV
积分
840

Level 4 (黑子)

25#
发表于 Post on 2021-4-13 22:06:07 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-23 09:04
是有限长度的。谢谢老师,我再仔细检查检查。构建邻居列表时用的group方式,这方式基于电荷组,但是GAFF ...

请问你当时是怎么解决这个问题的?我现在也遇到了,谷歌目前也没找到解决方法,是不是和并行有关系?

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

26#
 楼主 Author| 发表于 Post on 2021-4-14 10:39:01 | 只看该作者 Only view this author
王采岑 发表于 2021-4-13 08:06
请问你当时是怎么解决这个问题的?我现在也遇到了,谷歌目前也没找到解决方法,是不是和并行有关系?

mdrun用了-nt 1,可以计算了
live each day as if it was your last

53

帖子

0

威望

787

eV
积分
840

Level 4 (黑子)

27#
发表于 Post on 2021-4-14 10:44:38 | 只看该作者 Only view this author
yihanxu 发表于 2021-4-14 10:39
mdrun用了-nt 1,可以计算了

我用-nt 1 会出现下面的错误。
Fatal error:
The total number of threads requested (1) is not divisible by the number of
thread-MPI ranks requested (6)
我刚发了一个帖子问,在谷歌上搜到的是要用少于8个线程,可以不调用Domain decomposition,这个应该是MPI的问题,我还没弄清楚,对inux的知识少了,等社长或其他大佬回复,再来

384

帖子

0

威望

3364

eV
积分
3748

Level 5 (御坂)

28#
 楼主 Author| 发表于 Post on 2021-4-14 14:14:54 | 只看该作者 Only view this author
王采岑 发表于 2021-4-13 20:44
我用-nt 1 会出现下面的错误。
Fatal error:
The total number of threads requested (1) is not divis ...

嗯嗯 要是他们回复了请艾特我 谢谢!你可以先试试其他数目
live each day as if it was your last

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 22:21 , Processed in 0.622618 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list