| 查看: 19345 | 回复: 82 | ||||||||||||
| 本帖产生 2 个 模拟EPI ，点击这里进行查看 | ||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||||||||
[交流]
【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
|
||||||||||||
 为活跃本版气氛，加强本版专业实力， 分子模拟板块将于2011年4月25日---2011年4月29日举行“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”活动 活动主题：“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”活动 活动内容： 活动涉及内容包括经典/第一性原理分子动力学，Monte-Carlo，分子力学等方法的基本原理，相关软件或代码（如Discover-Studio，sybyl，autodock，Gromacs，NAMD，amber，CHARMM，LAMMPS，materials studio，DL_POLY,XMD,Materials-Explorer，CPMD 及其他数据可视化软件VMD,atomeye,xmovie,pymol等）的使用模拟出现的问题及各种资源共享，以及它们在纳米材料，表面/界面，块体材料和溶液，高分子，生物体系中的实际应用。 欢迎熟悉其中一个或多个方面的虫子参与本次活动。 活动形式：采用现场提问，现场解答的形式，活动结束后，将做出此活动专辑，形成一份PDF文件，供虫子们下载，活动结束后将由分子模拟版主委员会及群众投票评出”十大星虫“并给予适当金币奖赏。 奖励办法：500金币做为此活动奖赏，无论提问或是回答都得奖赏，只要回帖就领金币，提问奖1-3个金币，回答奖5-10个金币，对于经典回帖，经典解答，视情况分子模拟板块将再加奖赏金币，活动将由分子模拟版主委员会及群众投票评出”十大星虫“。 [ Last edited by ghcacj on 2011-4-25 at 09:25 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 | 第一性原理 | 课件资源j | 分子动力学 |
| 化学教材 | 模拟 | 学习 | 乱七八糟 |
| 杂七杂八 |
» 猜你喜欢
291求调剂
已经有15人回复
-
材料学硕333求调剂
已经有17人回复
-
070300化学求调剂
已经有5人回复
-
070300一志愿211，312分求调剂院校
已经有6人回复
-
物理学调剂
已经有3人回复
-
08工科，295，接受跨专业调剂
已经有5人回复
-
材料科学与工程调剂
已经有12人回复
-
本2一志愿C9-333分，材料科学与工程，求调剂
已经有5人回复
-
一志愿西交288化工专硕求调剂
已经有4人回复
-
308求调剂
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 上周分子模拟板块活动集锦PDF，下面提供一些候选人供大家评选“十大星虫”，欢迎投票
已经有16人回复
- Autodock对接构象有问题，请教
已经有7人回复
- 求助：vasp结果文件中的mag值为负值
已经有4人回复
- 中国矿业大学分子模拟【个人认为很经典额】【无重复】（更新到第十四章）
已经有141人回复
- MS分子动力学那个模块能处理断键和成键啊？
已经有9人回复
- castep最常出错原因和解决办法（已搜索，无重复，谢谢评价哦）
已经有7人回复
- Ubuntu 编译Lammps 单机串行错误，帮帮忙
已经有5人回复
- 【求助】谁会做这个图啊？
已经有6人回复
- 【求助】介孔炭吸附脱硫
已经有7人回复
- 【求助】加水后 碳纳米管不见了（已解决）
已经有4人回复
- 【求助】是不是monte carlo 模拟只适用于高分子方向的？
已经有6人回复
- 【求助】MS计算失败
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
五邑大学（江门-利物浦公共卫生研究所）招收药学、生物医药专业调剂生
+1/129
-
2026江西科技师范大学申亮组招硕士研究生
+1/86
-
南方医科大学中药学院 调剂硕士生一名 中药学专业
+1/81
-
坐标广州，征女友
+2/40
-
中山大学潘越教授现新增1个分子医学学术学位硕士名额
+1/39
-
球差相关知识
+1/37
-
国家工程物理交叉科学研究中心招收凝聚态物理及相关专业硕士研究生一名
+1/34
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收博士生
+1/17
-
山西师范大学2026年化学硕士点招收调剂生
+1/14
-
中科院和北京工商大学招收2026博士（化学背景）
+1/11
-
武汉工程大学化学与环境工程学院 2026 年硕士研究生调剂公告
+1/9
-
墨子实验室理论计算和模拟研究组诚聘海内外优秀人才
+1/5
-
墨尔本大学伯明翰大学联合全奖国际博士生
+1/4
-
齐齐哈尔大学化学与化学工程学院招收专硕调剂
+1/4
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/4
-
南京林业大学化学工程学院柏惺峰课题组招收硕士研究生1名
+1/3
-
美国田纳西理工 | 人工智能/机器人/自动驾驶方向（全奖）
+1/3
-
2026年北京石油化工学院环境学科-水污染控制工程方向课题组欢迎您
+1/1
-
三峡大学材化学院磷硅精细化工课题组招收化学类、材料类、材料与化工类硕士
+1/1
-
调剂来这边
+1/1
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zh1987hs(金币+10): 谢谢 2011-04-26 14:28:01
zh1987hs(金币+10): 谢谢 2011-04-26 14:28:01
|
例一 units metal # 单位为lammps 中的metel 类型 boundary p p p # 周期性边界条件 atom_style atomic # 原子模式 lattice fcc 3.61 # Cu 的晶格常数3.61 region box block 0 4 0 4 0 4 # x,y,z 各方向上的晶胞重复单元数，也即区域大小 create_box 1 box # 将上述区域指定为模拟的盒子 create_atoms 1 box # 将原子按晶格填满盒子 pair_style eam # 选取 Cu 的EAM 势作为模型 pair_coeff * * Cu_u3.eam # EAM 势文件名称 run 0 # 运行0 步，仅为启动lammps 的热力学数据计算 variable E equal pe # 定义变量 E 为系统总势能 variable N equal atoms # 定义变量 N 为系统总原子数 print "the number of atoms & system energy now are $N $E" # 打印信息 create_atoms 1 single 2.45 2.05 2.05 # 在该位置插入一个原子 min_style sd # 能量最小化模式，sd minimize 1.0e-12 1.0e-12 1000 1000 # 能量最小化参数，指数越大最小化程度越深 print "interstitial introduced, minimized: $N atoms, energy is $E" fix 1 all nvt 100 100 100 drag 0.2 # nvt 系综，原子数、体积和温度保持不变;T=100K timestep 0.005 # 步长 0.005fs run 1000 # 运行 1000 步 print "nvt performed, temperature up: $N atoms, total energy is $E" fix 1 all nvt 100 0.0001 100 drag 0.2 # nvt 系综，温度由100K 到0.0001K run 1000 # 运行 1000 步 print "nvt performed, temperature down: $N atoms, total energy is $E" compute 3 all pe/atom # 计算每个原子的势能 compute 4 all ke/atom # 计算每个原子的动能 compute 5 all coord/atom 3.0 # 计算每个原子的近邻原子数 dump 1 all custom 1 dump.atom id xs ys zs c_3 c_4 c_5 # 将信息写入dump.atom min_style sd minimize 1.0e-12 1.0e-12 10000 10000 # 再次能量最小化 print "the final state: $N atoms, total energy is $E" # 打印信息 ~ 例二 boundary p s s # 边界条件，拉伸方向是周期性，其余是自由边界；如果是薄膜拉伸则是两个周期性，块体则是三个周期性 units metal #单位制定义为metal atom_style atomic #原子类型自动 neighbor 2.0 bin #截断半径相关的东西 neigh_modifydelay 1 check yes # 邻近原子列表更新速度 #create geometry lattice fcc 3.61 # 定义晶胞为fcc，晶格常数3.61A region box block 0 30 0 3 0 3 #定义一个长方体区域叫box，长30，宽和高是3 create_box 1 box #创建了这样一个box create_atoms 1 box #在box里创建了一种原子 mass 1 63.546 #定义这种原子的质量是63.546 # potentials pair_style eam # 定义势函数是EAM pair_coef * * Cu_u3.eam #势所需要的参数在此文件里 #define groups region 1 block INF 1 INF INF INF INF # 定义了一个叫1的区域 group left region 1 #定义此区域里的原子叫left region 2 block 29 INF INF INF INF INF #定义了一个叫2的区域 group right region 2 #定义此区域里的原子叫right group boundary union left right #定义left+ right = boundary group mobile subtract all left # 定义mobile= all - left # initialvelocities velocity left set 0.0 0.0 0.0 # 设置原子初速度为0 computep all pressure thermo_temp # 计算应力，计算结果记为p variable pressx equal c_p[1] #定义变量pressx=c_p[1],c_p[1]的意思是p里第一个值 variable pressy equal c_p[2] #定义变量pressy=c_p[2],c_p[2]的意思是p里第二个值 variable pressz equal c_p[3] #定义变量pressz=c_p[3],c_p[3]的意思是p里第三个值 thermo_style custom step temp etotal press v_pressx v_pressy v_pressz vol dump 1 all atom 1000 dump.tensile # 输出结果到dump.tensile minimize 1.0e-6 1.0e-6 1000 1000 # 驰豫 fix 1 left setforce 0.0 NULL NULL #固定左边的原子 fix 2 all npt 1.0 1.0 1.0 aniso 0.0 0.0 NULL NULL NULL NULL 10.0 drag 1.0 thermo 1000 # 每1000步输出一次结果 timestep 0.002 #时间步长0.002ps run 10000 #运行10000步 unfix2 fix 3 all nve fix 4 all temp/rescale 100 1.0 1.0 5 1.0 fix 5 all deform 1000 x erate 0.002 units box #设置沿x方向拉伸，应变率为0.002(1/ps) compute s all stress/atom #计算每个原子上的应力，计算结果记为s compute strea all reduce sum c_s[1] #将s里第一个值求和，计算结果记为strea variable stressa equal c_strea/vol # 定义变量stressa= strea/vol, vol是体积 compute streb all reduce sum c_s[2] variable stressb equalc_streb/vol compute strec all reduce sum c_s[3] variable stressc equal c_strec/vol thermo_style custom step tempetotal press v_stressa v_stressb v_stressc vol run 300000 |
11楼2011-04-25 11:02:36
2楼2011-04-25 09:27:16
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+3): 只要提问就奖 2011-04-25 09:39:08
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+3): 只要提问就奖 2011-04-25 09:39:08
|
顶lz 我有个问题困扰了我半个月了。 非晶态淬火的相关问题。 比如非晶碳。 文献上说“液相--快淬”制非晶。 我用discover操作，不能改变键的配位，不能断键、成键。 所以用Dmol3来做，可是里面不知道怎么淬火！具体怎么操作。 Dmol3里的dynamics只有一个温度选项。而且计算非常慢。怎么淬火？ （1）Dmol3要成键，段键，改变配位，怎样设置参数？ （2）Dmol3怎样淬火？是否先用discover淬火，然后用Dmol3成键？ |
3楼2011-04-25 09:37:14
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 分享资源重奖，欢迎再来 2011-04-25 09:46:01
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 分享资源重奖，欢迎再来 2011-04-25 09:46:01
|
原帖地址http://muchong.com/bbs/viewthread.php?tid=3119236&pid=118825&page=1#pid118825 中国矿业大学分子模拟【个人认为很经典额】课件 对照ms的discover模块，很容易上手， 各种原理，各种解释。 新手必备，老鸟提高必备！ 能让人形成系统的知识体系！ 废话不多说，上连接，总共11个文件 下载链接无需注册，直接下载，方便虫子们！ 1、http://ishare.iask.sina.com.cn/f/14872882.html 2、http://ishare.iask.sina.com.cn/f/14872883.html 3、http://ishare.iask.sina.com.cn/f/14872884.html 4、http://ishare.iask.sina.com.cn/f/14872888.html 5、http://ishare.iask.sina.com.cn/f/14872889.html 6、http://ishare.iask.sina.com.cn/f/14872998.html 7、http://ishare.iask.sina.com.cn/f/14873001.html 8、http://ishare.iask.sina.com.cn/f/14873002.html 9、http://ishare.iask.sina.com.cn/f/14872999.html 10、http://ishare.iask.sina.com.cn/f/14873000.html 11、http://ishare.iask.sina.com.cn/f/14873019.html 附送MS castep 计算化学高级教程 、中国石油大学+创腾培训.pdf http://ishare.iask.sina.com.cn/f/14873015.html http://ishare.iask.sina.com.cn/f/14873018.html |
4楼2011-04-25 09:41:37
