24小时热门版块排行榜    

查看: 2845  |  回复: 18

Johnnase

木虫 (正式写手)

引用回帖:
9楼: Originally posted by lkd-syl at 2014-06-19 20:48:15
你可以把输入文件贴出来看看。。。。...

我用的是Discover模块,选的是NPT系统,1290K,Verlocity Scaling控温、Andersen控压方法,一共运行了30Ps。
Inp文件如下:#BIOSYM btcl 3
#
#  Input File For Discover Generated By Materials Studio
#  Input Client Model Document: D:\MS mission\Alumium droplet Files\Documents\Al.xsd
#  Job: [RJY9Z] - Al Disco Dynamics
#
  autoEcho off
#
# Begin Forcefield Section
  begin forcefield = compass
#
# Nonbond section:
  forcefield nonbond \
      -separate_coulomb \
    vdw \
    summation_method = atom_based \
      cutoff = 9.50 \
      spline_width = 1.00 \
      buffer_width = 0.50 \
    coulomb \
    dielectric_value = 1.0000
#
# vdW long range energy correction
  vdWTailCorrection 9.50
# End Forcefield Section
#
# Dynamics Section:
  dynamics \
    time = 30000.00 \
    timestep = 1.00 \
    initial_temperature = 1290.00 \
    ensemble = NPT \
    temperature_control_method = velocity_scaling \
    temperature_window = 100.00 \
    temperature = 1290.00 \
    press_choice = pressure \
    pressure = 0.00010 \
    pressure_control_method = andersen_cp \
    cell_mass = 20.00 \
    deviation = 5000.00 \
    execute frequency = 5000.00 \
      command = {print history}

这中情况下得到的液体胞混乱度比较好,但是能量如之前回复您的那张图所示,是随着模拟时间增大后保持不变的。
请问可否告知您在对超晶胞做分子动力学时候的具体设置:主要是控温方法和控压方法,模拟时长,以及得到的能量随时间变化的曲线,我希望知道什么样的能量——时间变化曲线是正确的。
不胜感激!
Jugglers and singersr equire applause!
11楼2014-06-20 15:46:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lkd-syl

金虫 (著名写手)

快乐家族

引用回帖:
11楼: Originally posted by Johnnase at 2014-06-19 19:46:57
我用的是Discover模块,选的是NPT系统,1290K,Verlocity Scaling控温、Andersen控压方法,一共运行了30Ps。
Inp文件如下:#BIOSYM btcl 3
#
#  Input File For Discover Generated By Materials Studio
#  ...

你用那种控温方式是在看别人做的用的吗?  我见过用Brendsen的,另外那个能量升高后稳定在一个定值波动,这个没什么问题吧,因为你做的是一个升温动力学,能量升高也是理所当然的,,,
Nothingimpossible!!
12楼2014-06-21 11:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Johnnase

木虫 (正式写手)

引用回帖:
12楼: Originally posted by lkd-syl at 2014-06-21 11:30:11
你用那种控温方式是在看别人做的用的吗?  我见过用Brendsen的,另外那个能量升高后稳定在一个定值波动,这个没什么问题吧,因为你做的是一个升温动力学,能量升高也是理所当然的,,,...

首先感谢您的回复!
我各种控温的方式都尝试过了,包括nose、andersen和berendsen的,但是就混乱度而言我发现用verlocity scale的最好,其他几种控温方式混乱度比较差,就如同之前那张图上展示的。我看到过网上有个教程用的是nose控温方式,网址链接如下:http://www.docin.com/p-484673708.html
还有,您说我的是做的升温动力学?那么应当怎样设置为恒温动力学呢?我其实也是想做恒温动力学的。还请指导一下,不胜感激!
另外还有一点就是,选择NPT的话,那么压力应当控制为一个大气压吗?还是应当是别的数值?
Jugglers and singersr equire applause!
13楼2014-06-21 12:22:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lkd-syl

金虫 (著名写手)

快乐家族

引用回帖:
13楼: Originally posted by Johnnase at 2014-06-20 16:22:07
首先感谢您的回复!
我各种控温的方式都尝试过了,包括nose、andersen和berendsen的,但是就混乱度而言我发现用verlocity scale的最好,其他几种控温方式混乱度比较差,就如同之前那张图上展示的。我看到过 ...

我说的意思你可能没明白,我上边说的升温过程是从一个可能是室温(建立超胞的状态,此状态晶格规则排列,具体条件是什么我还没弄明白,但后边目标温度在1290K的话这个是不影响结果的,不管初始时300还是400K,,,)的构型,然后直接升到1290K,在此温度下进行的一个弛豫过程,而不是说从300K慢慢升到1290K的过程,1290K相比于之前的(比如300K)就是一个升温过程,此过程能量是升高的,关于压力的话好像也没有明确的界定,没有特殊要求的话应该还是在一个标准条件下进行的,,,
Nothingimpossible!!
14楼2014-06-21 12:40:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

聆听世界

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by lkd-syl at 2014-05-12 16:25:24
做倒是好几个模块都可以做的,我做的铁液凝固,用的GULP,貌似castep,dmol3都可以,关键要看你想做什么了,,,...

你好,考虑钢液与Al2O3粒子间润湿可以用GULP做吗,我是新手,望解答
15楼2014-11-02 15:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)

引用回帖:
5楼: Originally posted by Johnnase at 2014-05-12 20:49:58
我其实想做铝液液滴在石墨基体上的润湿模拟,想先建立铝液液滴的模型。...

你好!
很开心发现你跟我是同个研究方向,我是中南的,看到你的帖子,液态铝在石墨上润湿,想必师兄你是搞电解铝的吧!我也是,能不能把你的石墨模型结构cif文件和模拟所用的gin文件给我一份啊! 不甚感激
16楼2015-02-17 15:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Johnnase

木虫 (正式写手)

引用回帖:
16楼: Originally posted by 15216105346 at 2015-02-17 15:20:32
你好!
很开心发现你跟我是同个研究方向,我是中南的,看到你的帖子,液态铝在石墨上润湿,想必师兄你是搞电解铝的吧!我也是,能不能把你的石墨模型结构cif文件和模拟所用的gin文件给我一份啊! 不甚感激...

我们专业应该是同一个方向,不知道你具体是做什么的?
石墨的模型在MS晶体库里面就有,你只要用适当的方法就可以切成方形的就可以了。
Jugglers and singersr equire applause!
17楼2015-03-19 10:31:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蜉蝣生物

金虫 (小有名气)

引用回帖:
9楼: Originally posted by lkd-syl at 2014-06-19 20:48:15
你可以把输入文件贴出来看看。。。。...

师兄我想问一下建立超晶胞后需要对其进行GULP的几何优化么
18楼2017-10-31 22:40:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钱某人

新虫 (正式写手)

温度升高能量增大不是正常的吗?

发自小木虫Android客户端
19楼2017-11-01 01:13:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Johnnase 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 我感觉当老师好累啊。。 +10 fairy1122 2024-11-15 10/500 2024-11-15 23:13 by 鱼翔浅底1
[考博] 随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利 +6 Ricoch4t 2024-11-13 15/750 2024-11-15 19:24 by 青青之洲
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[论文投稿] 职称论文 +10 xiaoqi_66 2024-11-09 13/650 2024-11-15 16:45 by sugarmei
[精细化工] 同一个反应相同的反应条件,是不是反应结果相差不大? +5 青霉素 2024-11-11 5/250 2024-11-15 13:43 by zyqchem
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[教师之家] 北大教授何怀宏曾如此描述他的同行 +12 zju2000 2024-11-09 12/600 2024-11-15 12:13 by PJ女神
[硕博家园] 和导师撕破脸了 +16 圆满123456 2024-11-13 20/1000 2024-11-15 10:09 by 小杰1572
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+3 猪小耍 2024-11-13 12/600 2024-11-15 08:41 by 北京莱茵润色
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 二审审一年的佛系编辑 +10 thefoxrain 2024-11-09 15/750 2024-11-12 19:27 by lide966
[论文投稿] 论文接受后,编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
[硕博家园] 同步辐射 +4 小张要加油努力 2024-11-10 7/350 2024-11-11 13:47 by mpdfwxgui
[论文投稿] renewable energy 添加作者 15+3 Hebauwww12 2024-11-09 4/200 2024-11-11 08:42 by 北京莱茵润色
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见