计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1906|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求助:生成em.tpr文件时出现严重警告,导致错误

[复制链接 Copy URL]

23

帖子

0

威望

49

eV
积分
72

Level 2 能力者

运行如下命令时:gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr
出现以下严重警告,最终导致报错:
WARNING 1 [file topol.top, line 11187]:  You are using Ewald electrostatics in a system with net charge. This can
  lead to severe artifacts, such as ions moving into regions with low
  dielectric, due to the uniform background charge. We suggest to
  neutralize your system with counter ions, possibly in combination with a
  physiological salt concentration.
报错内容:
Fatal error:
Too many warnings (1).
If you are sure all warnings are harmless, use the -maxwarn option.

请问我可以直接利用-maxwarn命令跳过这个警告吗?
结尾附上运行完整内容:
[root@192 CXCL8_5d]# gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr
                      :-) GROMACS - gmx grompp, 2020.6 (-:

                            GROMACS is written by:
     Emile Apol      Rossen Apostolov      Paul Bauer     Herman J.C. Berendsen
    Par Bjelkmar      Christian Blau   Viacheslav Bolnykh     Kevin Boyd   
Aldert van Buuren   Rudi van Drunen     Anton Feenstra       Alan Gray     
  Gerrit Groenhof     Anca Hamuraru    Vincent Hindriksen  M. Eric Irrgang  
  Aleksei Iupinov   Christoph Junghans     Joe Jordan     Dimitrios Karkoulis
    Peter Kasson        Jiri Kraus      Carsten Kutzner      Per Larsson   
  Justin A. Lemkul    Viveca Lindahl    Magnus Lundborg     Erik Marklund   
    Pascal Merz     Pieter Meulenhoff    Teemu Murtola       Szilard Pall   
    Sander Pronk      Roland Schulz      Michael Shirts    Alexey Shvetsov  
   Alfons Sijbers     Peter Tieleman      Jon Vincent      Teemu Virolainen
Christian Wennberg    Maarten Wolf      Artem Zhmurov   
                           and the project leaders:
        Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel

Copyright (c) 1991-2000, University of Groningen, The Netherlands.
Copyright (c) 2001-2019, The GROMACS development team at
Uppsala University, Stockholm University and
the Royal Institute of Technology, Sweden.
check out http://www.gromacs.org for more information.

GROMACS is free software; you can redistribute it and/or modify it
under the terms of the GNU Lesser General Public License
as published by the Free Software Foundation; either version 2.1
of the License, or (at your option) any later version.

GROMACS:      gmx grompp, version 2020.6
Executable:   /opt/gromacs-2020.6/bin/gmx_mpi
Data prefix:  /opt/gromacs-2020.6
Working dir:  /home/moi/gromacs_MD/1/CXCL8_5d
Command line:
  gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr

Ignoring obsolete mdp entry 'title'
Ignoring obsolete mdp entry 'ns_type'

NOTE 1 [file em.mdp]:
  With Verlet lists the optimal nstlist is >= 10, with GPUs >= 20. Note
  that with the Verlet scheme, nstlist has no effect on the accuracy of
  your simulation.

Setting the LD random seed to -549787922
Generated 3003 of the 3003 non-bonded parameter combinations
Generating 1-4 interactions: fudge = 0.5
Generated 3003 of the 3003 1-4 parameter combinations
Excluding 3 bonded neighbours molecule type 'Protein_chain_A'
Excluding 3 bonded neighbours molecule type 'ligand'
Excluding 2 bonded neighbours molecule type 'SOL'
Excluding 1 bonded neighbours molecule type 'CL'

NOTE 2 [file topol.top, line 11187]:
  System has non-zero total charge: 0.001000
  Total charge should normally be an integer. See
  http://www.gromacs.org/Documentation/Floating_Point_Arithmetic
  for discussion on how close it should be to an integer.




WARNING 1 [file topol.top, line 11187]:
  You are using Ewald electrostatics in a system with net charge. This can
  lead to severe artifacts, such as ions moving into regions with low
  dielectric, due to the uniform background charge. We suggest to
  neutralize your system with counter ions, possibly in combination with a
  physiological salt concentration.


++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++
J. S. Hub, B. L. de Groot, H. Grubmueller, G. Groenhof
Quantifying Artifacts in Ewald Simulations of Inhomogeneous Systems with a Net
Charge
J. Chem. Theory Comput. 10 (2014) pp. 381-393
-------- -------- --- Thank You --- -------- --------

Analysing residue names:
There are:    70    Protein residues
There are:     1      Other residues
There are:  6990      Water residues
There are:     4        Ion residues
Analysing Protein...
Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...
Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...
Number of degrees of freedom in T-Coupling group rest is 45582.00
Calculating fourier grid dimensions for X Y Z
Using a fourier grid of 52x52x52, spacing 0.117 0.117 0.117
Estimate for the relative computational load of the PME mesh part: 0.18
This run will generate roughly 2 Mb of data

There were 2 notes

There was 1 warning

-------------------------------------------------------
Program:     gmx grompp, version 2020.6
Source file: src/gromacs/gmxpreprocess/grompp.cpp (line 2353)

Fatal error:
Too many warnings (1).
If you are sure all warnings are harmless, use the -maxwarn option.

For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
-------------------------------------------------------
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.





206

帖子

0

威望

1153

eV
积分
1359

Level 4 (黑子)

终身学习

2#
发表于 Post on 2024-1-26 14:40:04 | 只看该作者 Only view this author
系统的总电荷未中和,可以参考这里的回答:
http://bbs.keinsci.com/thread-19226-1-1.html
Open source enables open science.

23

帖子

0

威望

49

eV
积分
72

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-1-26 14:45:07 | 只看该作者 Only view this author
Seyilaxa 发表于 2024-1-26 14:40
系统的总电荷未中和,可以参考这里的回答:
http://bbs.keinsci.com/thread-19226-1-1.html

请问如果提示System has non-zero total charge: 0.001000
是不是说明系统总电荷量小呢?可不可以用maxwarn来忽略这个警告呢?

206

帖子

0

威望

1153

eV
积分
1359

Level 4 (黑子)

终身学习

4#
发表于 Post on 2024-1-26 15:58:30 | 只看该作者 Only view this author
MOI_001 发表于 2024-1-26 14:45
请问如果提示System has non-zero total charge: 0.001000
是不是说明系统总电荷量小呢?可不可以用maxw ...

可以检查一下具体是哪一部分的电荷出现问题,一般影响不大可以-maxwarn忽略掉,看看模拟结果是否有异常。实在担心的话在拓扑文件里把-0.001的电荷随便分配给一个原子就行。
Open source enables open science.

23

帖子

0

威望

49

eV
积分
72

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2024-1-26 22:05:41 | 只看该作者 Only view this author
Seyilaxa 发表于 2024-1-26 15:58
可以检查一下具体是哪一部分的电荷出现问题,一般影响不大可以-maxwarn忽略掉,看看模拟结果是否有异常。 ...

好的 感谢

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6#
发表于 Post on 2024-1-27 00:44:35 | 只看该作者 Only view this author
参考北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/workshop/KGMX_content.html)的ppt:
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 20:43 , Processed in 0.172356 second(s), 23 queries , Gzip On.

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