计算化学公社

 找回密码 Forget password
 注册 Register

求助ADMP相关问题

查看数: 10749 | 评论数: 17 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-11-24 15:26

正文摘要:

本帖最后由 niubaotian85 于 2017-11-24 15:32 编辑 这是我在一篇文献上看到的一个分析手段,计算时在高斯中输入#p MP2/6-311++G** ADMP(maxpoint=1000,stepsize=2000,fullscf) temperature=100.0 pressure=1 ...

回复 Reply

ACS_Nano 发表于 Post on 2025-6-4 22:18:37
sobereva 发表于 2025-6-4 17:58
你不设置显示键长当然不会默认显示

参考北京科音高级量子化学培训班(http://www.keinsci.com/KAQC) ...

知道了,谢谢sob老师。
sobereva 发表于 Post on 2025-6-4 17:58:57
ACS_Nano 发表于 2025-6-4 17:52
您好,麻烦问下你最后是怎么解决查看resutls -> Trajectory -> Plots的,我的情况和你一样,GV版本是5.0. ...

你不设置显示键长当然不会默认显示

参考北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)极其详细讲AIMD模拟部分的幻灯片:


gview 5已经很老了,没有这个功能。如果没有gview 6,用下文的方法转成xyz轨迹,然后在VMD里得到键长的变化

GauMD2xyz:将Gaussian的从头算动力学轨迹转换为xyz轨迹的程序
http://sobereva.com/126

ACS_Nano 发表于 Post on 2025-6-4 17:52:56
niubaotian85 发表于 2017-11-24 22:57
读取键长后得到的曲线很不一样。文献中是振动的曲线,我的是一条上升的近乎直线,如果可以,请帮我看一看

您好,麻烦问下你最后是怎么解决查看resutls -> Trajectory -> Plots的,我的情况和你一样,GV版本是5.0.9,打开ADMP的log文件,但是只有resutls -> Trajectory,一打开就是很多能量啥的,没有键长这个。非常感谢。
niubaotian85 发表于 Post on 2017-11-25 21:01:34
yflchx 发表于 2017-11-25 20:44
好像他们做的是恒温模拟,要用到Nuclear Kinetic Energy Thermostat Option,你写temperature的话根本没 ...

十分感谢,谢谢谢谢
yflchx 发表于 Post on 2017-11-25 20:44:52
本帖最后由 yflchx 于 2017-11-25 20:50 编辑
niubaotian85 发表于 2017-11-25 19:47
是这样的,文献里没有输入,我不知道输入对不对,MP2方法和对应的基组是没有问题的,他做了100k和300k的 ...

好像他们做的是恒温模拟,要用到Nuclear Kinetic Energy Thermostat Option,你写temperature的话根本没用。

去看G09的IOP手册:IOP(1/80-82,89,91,92),在IOP手册里面搜索thermostat,以防我给你写的不对。

我只能帮到这里了,八九年前我用过ADMP,现在没有太多空闲时间去研究这个。

实在不行,一个Email过去,请求作者帮帮忙。
yflchx 发表于 Post on 2017-11-25 18:27:43
niubaotian85 发表于 2017-11-25 16:26
因为我没做过ADMP,可能问的不专业,希望见谅。就是,作者在文献中讨论动力学模拟的时候,看不同温度下是否 ...

这个不好做判断,我没看过原文。不知道作者是如何讨论的。

假如输入来自文献的话,MP2、temperature、pressure我感觉很奇怪。
niubaotian85 发表于 Post on 2017-11-25 16:26:03
本帖最后由 niubaotian85 于 2017-11-25 16:42 编辑
yflchx 发表于 2017-11-25 09:18
人家描绘的是两个H之间的键长,你描绘的是Be(5)和H(6)的,两码事。
因为我没做过ADMP,可能问的不专业,希望见谅。就是,作者在文献中讨论动力学模拟的时候,看不同温度下是否被吸附,那是不是应该考虑Be与H2间的距离呢?考虑H2间原子的距离能说明它们被吸附上了没有么?文献中说了一句the trajectories of H2 molecules from Be atom in C2H2Be(2H2) at 100K.也让我感觉像是说与Be原子的距离。而且我发现Be原子和H2的距离是在1.5埃左右。所以以上的疑问,外加我做出的结果,如果是H2间H-H的距离又和文献中不一样,所以还是想问问你,是不是我哪里有问题。非常感谢,真是太麻烦你了。
niubaotian85 发表于 Post on 2017-11-25 16:15:06
本帖最后由 niubaotian85 于 2017-11-25 16:23 编辑
yflchx 发表于 2017-11-25 09:18
人家描绘的是两个H之间的键长,你描绘的是Be(5)和H(6)的,两码事。

哦,是这样,但是尝试后发现,我做的是在0.7-0.82振动,但是文献是在1.5-1.8埃间振动,这不合理啊,一样的方法,却答案不一样,是为什么呢?
yflchx 发表于 Post on 2017-11-25 09:18:03
niubaotian85 发表于 2017-11-24 22:57
读取键长后得到的曲线很不一样。文献中是振动的曲线,我的是一条上升的近乎直线,如果可以,请帮我看一看

人家描绘的是两个H之间的键长,你描绘的是Be(5)和H(6)的,两码事。
niubaotian85 发表于 Post on 2017-11-24 22:38:03
yflchx 发表于 2017-11-24 21:41
为什么不试试读取.log或.out呢。

我试过了,好着呢。最后一个是close,倒数第二个是:Plot Molecular  ...

太感谢了
yflchx 发表于 Post on 2017-11-24 21:41:05
niubaotian85 发表于 2017-11-24 21:32
真的非常感谢,麻烦你了,我用的是GV5.0.9,我用的fchk文件,读取时有Trajectory这一项,但是鼠标放在上 ...

为什么不试试读取.log或.out呢。

我试过了,好着呢。最后一个是close,倒数第二个是:Plot Molecular Property...点开就是定义描绘键长。
niubaotian85 发表于 Post on 2017-11-24 21:32:39
yflchx 发表于 2017-11-24 20:39
你用的GV版本比较低,用GV5.0.8可以。

真的非常感谢,麻烦你了,我用的是GV5.0.9,我用的fchk文件,读取时有Trajectory这一项,但是鼠标放在上面没有其他可以选择的了,只能直接点开,所以我才怀疑是不是我命令行错了。请问该怎么操作。
yflchx 发表于 Post on 2017-11-24 20:39:25
niubaotian85 发表于 2017-11-24 20:31
谢谢你,但是我在读取的时候是fchk,在result里可以选择Trajectory,然后没有可以选择的,直接点开,是十 ...

你用的GV版本比较低,用GV5.0.8可以。

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

GMT+8, 2025-8-12 22:54 , Processed in 0.901186 second(s), 25 queries , Gzip On.

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