应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】DL_POLY的config文件如何生成??
996/10返回列表首页上一页456789下一页末页
查看: 5089  |  回复: 98
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 模拟EPI ，点击这里进行查看

zyj8119

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by tpp001 at 2010-03-22 22:18:53:





program dlpolyinput
implicit none
integer i, numcoordinate, pbc, numatom
real, PARAMETER :: K=0
character *20 atomname, systemname
real x(200),y(200),z(200), lbox
open(20,file='CONF ...

谢谢啊，兄弟，学习了，我也编了程序，就是运行不出来。
赞一下(1人) 回复此楼
好好学习，天天向上。
51楼2010-03-22 22:20:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zyj8119 at 2010-03-22 22:20:18:

谢谢啊，兄弟，学习了，我也编了程序，就是运行不出来。

此贴很好，其实做这种程序类型软件，需要自己编程能力要好一点的，大家多多交流。。。
赞一下(1人) 回复此楼
好好学习，天天向上。
52楼2010-03-22 22:35:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
引用回帖:
Originally posted by tpp001 at 2010-03-22 22:20:00:
input-turn-config  
0         1       200
28.850
C1        -7.554        -6.98        -10.036
C2        -8.917        -7.191        -10.626
C2        -8.954        -8.558        -11.308
C2        -10.315        -8.796        -11.93
S14        -7.47        -5.332        -9.259
C1        1.742        -9.347         ...

这列数据，应该是使用EXCEL分列形成的吧，在MS里面导出的PDB文件不是这样的。
赞一下 回复此楼
好好学习，天天向上。
53楼2010-03-22 22:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by zyj8119 at 2010-03-22 22:40:45:

这列数据，应该是使用EXCEL分列形成的吧，在MS里面导出的PDB文件不是这样的。

恩,正确
回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
54楼2010-03-23 10:25:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by 老虎大王 at 2010-03-22 12:08:19:


不知道，应该是有的吧。你去他的网站上看。

不好意思,我这里还是有不少问题.....
现在我改用了3.09版本的,进行了一个npt的模拟,我的 体系是40个丁硫醇....  
我的体系运行出现下面错误...
------------------------------
DL_POLY_3 terminated due to error   445

error - r_14 > rcut in dihedrals_forces
--------------------------------
我把cutoff 加大(一直到盒子长的一半)还是不行....
下面是我 的control文件
integrator leapfrog verlet
temperature      300.00
pressure         0.001  
ensemble npt ber 0.2   0.2

steps               10000
multiple step         1
scale                10
print                50
stack               100
stats                10
traj             1000 100 0

timestep         0.0010
cutoff           9.5000
rvdw cutoff      9.5000
no electrostatics

job time              18000.00
close time            300.00
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
55楼2010-03-23 13:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3):谢谢 2010-03-23 15:19
出现这个错误，是说你的r_14太大了，超过了总的短程力截断半径rcut。你要考虑是否可以把dihedral_force 力场里面的 1-4 distance 减小，如果无法减小，就增大rcut。如果Rcut也无法再大了，已经到了盒子的一半，那么我建议你增大盒子，也就是多算一些分子。当然这一切都是在你没有其它错误的前提之下。其实你首先要仔细检查一下你的势参数的设置是否正确。
赞一下(2人) 回复此楼
56楼2010-03-23 14:10:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by 老虎大王 at 2010-03-23 14:10:23:
出现这个错误，是说你的r_14太大了，超过了总的短程力截断半径rcut。你要考虑是否可以把dihedral_force 力场里面的 1-4 distance 减小，如果无法减小，就增大rcut。如果Rcut也无法再大了，已经到了盒子的一半，那 ...

您好,为了找到我体系运行出错的原因,我用和以前一样的方法建立CONFIG,跑了一个包含100个水分子的纯水体系.在300k 1atm下.生成了HISTORY文件,我用WMD查看.
初始的100个水分子为下面形态:发现盒子边缘有些水分子的化学键断开了

运行结束后:发现有些水分子的化学键变的很长


我觉得可能是我的初始构型:就是建立CONFIG的方法不对. ....??
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
57楼2010-03-25 15:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+1):谢谢 2010-03-25 20:11
引用回帖:
Originally posted by tpp001 at 2010-03-25 15:55:14:

您好,为了找到我体系运行出错的原因,我用和以前一样的方法建立CONFIG,跑了一个包含100个水分子的纯水体系.在300k 1atm下.生成了HISTORY文件,我用WMD查看.
初始的100个水分子为下面形态:发现盒子边缘有些水分子 ...

方法应该没有什么问题，有时候，一个小小的格式错误都会导致问题
赞一下(2人) 回复此楼
好好学习，天天向上。
58楼2010-03-25 19:29:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by zyj8119 at 2010-03-25 19:29:06:

方法应该没有什么问题，有时候，一个小小的格式错误都会导致问题

不知道你用ms的pdb转化过来的config文件运行以后出现过这样的情况吗?望回复thanks
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
59楼2010-03-25 20:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+1):谢谢 2010-03-27 22:17
你第二张图我没有看懂。你计算键长的时候有没有用到周期性边界条件下的邻像变换？
赞一下(2人) 回复此楼
60楼2010-03-25 21:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tpp001 的主题更新
996/10返回列表首页上一页456789下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +15 低垂的野花 2026-07-31 20/1000 2026-08-01 14:10 by gltch
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么？ 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭