计算化学公社

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

[GROMACS] 大家好,我在用gromacs进行模拟时,出现如下问题,本人小白,求高人指点

[复制链接 Copy URL]

13

帖子

0

威望

88

eV
积分
101

Level 2 能力者

本帖最后由 zouguo 于 2024-5-2 14:30 编辑

1. 主要目的:用gromacs对小分子金属配合物进行分子动力学研究,查看其聚集态堆积模式(盒子里只有这一种金属配合物的存在,也不添加溶剂分子)
2. 首先用高斯构建分子模型,在进行opt freq计算得分子基态最优构型,得到fchk文件,用gv6打开,再另存为mol2文件或pdb文件
3. 用sob老师的sobtop将mol2文件转化成跑gromacs要用的.top.itp.gro文件(打开sobtop,将mol2拖进去,分别输入1、1、0、4回车,在填入路径,回车,即可得到想要的文件s,不得不说Sob,YYDS呀),准备用以下代码,来完成模拟:
gmx editconf -f 3pt.gro -o 1_newbox.gro -c -d 1.0 -bt cubic (没问题)
gmx grompp -f em.mdp -c 1_newbox.gro -p1.top -o em.tpr -maxwarn 1 (出错了)
gmx mdrun -v -deffnm em
gmx grompp -f nvt.mdp -c em.gro -ptopol.top -r 3pt.gro-o nvt.tpr -maxwarn 2
gmx mdrun -deffnm nvt
gmx grompp -f npt.mdp -c nvt.gro -t nvt.cpt-p topol.top -r 3pt.gro-o npt.tpr -maxwarn 2
gmx mdrun -deffnm npt
gmx editconf -f npt.gro -o npt.pdb

4. 进入到上述文件所在文件夹,运行win10下Ubuntu子系统,给分子添加盒子:gmx editconf -f 1.gro -o 1_newbox.gro -c -d 1.0 -bt cubic
   这一步倒没什么问题,生成了1_newbox.gro文件,且能够用VMD打开,输入Pbc box,回车后能看到盒子


5. 将1_newbox.gro按照能量最小化参数em.mdp产生em.tpr文件,这一步出错,且产生了一个名为mdout.mdp的文件,详情如下:
zoug@DESKTOP-EOKETI3:/mnt/g/Desktop/mdp7$ gmx grompp -f em.mdp -c 1_newbox.gro -p 1.top -o em.tpr -maxwarn 1
              :-) GROMACS - gmx grompp, 2021.4-Ubuntu-2021.4-2 (-:

                            GROMACS is written by:
     Andrey Alekseenko              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
    Gilles Gouaillardet             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 2021.4-Ubuntu-2021.4-2
Executable:   /usr/bin/gmx
Data prefix:  /usr
Working dir:  /mnt/g/Desktop/mdp7
Command line:
  gmx grompp -f em.mdp -c 1_newbox.gro -p 1.top -o em.tpr -maxwarn 1

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 2096099290

-------------------------------------------------------
Program:     gmx grompp, version 2021.4-Ubuntu-2021.4-2
Source file: src/gromacs/gmxpreprocess/gmxcpp.cpp (line 298)

Fatal error:
Topology include file "G:\Desktop\mdp7\1.itp" not found

For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
-------------------------------------------------------
所有相关文件已上传,麻烦高人指点迷津,谢谢!
祝好
gromacs小白
20240502

1.gro

6.97 KB, 下载次数 Times of downloads: 0

1.itp

133.38 KB, 下载次数 Times of downloads: 0

1.mol2

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

1.top

343 Bytes, 下载次数 Times of downloads: 4

1_newbox.gro

6.76 KB, 下载次数 Times of downloads: 0

em.mdp

1.05 KB, 下载次数 Times of downloads: 4

网上copy后修改的

mdout.mdp

10 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

2#
发表于 Post on 2024-5-2 23:22:14 | 只看该作者 Only view this author
其它文字我没看,红字提示明显说明路径没弄对。这种基本问题凭提示解决
北京科音自然科学研究中心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

13

帖子

0

威望

88

eV
积分
101

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-5-6 09:39:49 | 只看该作者 Only view this author
sobereva 发表于 2024-5-2 23:22
其它文字我没看,红字提示明显说明路径没弄对。这种基本问题凭提示解决

谢谢sob老师!

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 23:58 , Processed in 0.165988 second(s), 23 queries , Gzip On.

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