| 查看: 18413 | 回复: 82 | |||||||||||
| 本帖产生 2 个 模拟EPI ,点击这里进行查看 | |||||||||||
[交流]
【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
|
|||||||||||
资源收集 | 第一性原理 | 课件资源j | 分子动力学 |
化学教材 | 模拟 | 学习 | 乱七八糟 |
杂七杂八 |
» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
到新单位后,换了新的研究方向,没有团队,持续积累2区以上论文,能申请到面上吗
已经有8人回复
申请2026年博士
已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
2025冷门绝学什么时候出结果
已经有7人回复
请问有评职称,把科研教学业绩算分排序的高校吗
已经有6人回复
Bioresource Technology期刊,第一次返修的时候被退回好几次了
已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
上周分子模拟板块活动集锦PDF,下面提供一些候选人供大家评选“十大星虫”,欢迎投票
已经有16人回复
Autodock对接构象有问题,请教
已经有7人回复
求助:vasp结果文件中的mag值为负值
已经有4人回复
中国矿业大学分子模拟【个人认为很经典额】【无重复】(更新到第十四章)
已经有141人回复
MS分子动力学那个模块能处理断键和成键啊?
已经有9人回复
castep最常出错原因和解决办法(已搜索,无重复,谢谢评价哦)
已经有7人回复
Ubuntu 编译Lammps 单机串行错误,帮帮忙
已经有5人回复
【求助】谁会做这个图啊?
已经有6人回复
【求助】介孔炭吸附脱硫
已经有7人回复
【求助】加水后 碳纳米管不见了(已解决)
已经有4人回复
【求助】是不是monte carlo 模拟只适用于高分子方向的?
已经有6人回复
【求助】MS计算失败
已经有5人回复
» 抢金币啦!回帖就可以得到:
Ei期刊青年编委招募(工程设计方向)
+1/391
鲁东大学硕士研究生招生
+1/134
湘潭大学化学学院陈华杰教授课题组招收有机/高分子方向的博士研究生
+1/81
中国地质大学(武汉)—国家级青年人才杨明教授组-招收博士-新能源材料化学及催化材料
+1/77
坐标山东东营,诚征女友
+1/59
浙江农林大学森林食物资源挖掘与利用全国重点实验室2026年博士生招生
+1/44
校长团队招博士生和博士后
+1/34
国家青年人才叶立群教授课题组招收2026级博士研究生
+1/32
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/31
今年的函评专家费收到了吗
+1/23
限广州,征女友
+2/10
博士/硕士招生
+1/9
招若干有分子生物,细胞培养,动物实验背景的人员(中山大学)
+1/6
招若干有机合成人员 (中山大学)
+1/6
香港城市大学 Microsystems and Semiconductor Technology Lab 博士生招生启事
+1/5
中山大学院士团队王来源教授课题组招聘博士后
+2/4
CSC 因斯布鲁克大学计算机系 联培/攻博
+1/3
有没有一款可以听文献的APP
+1/2
招收2026年秋季入学博士生1名(北京科技大学 力学超材料/机器学习/增材制造相关方向)
+1/2
中国科学技术大学 精准智能化学重点实验室 武建昌课题组招聘博士,博士后
+1/1
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
5楼2011-04-25 09:49:04
6楼2011-04-25 09:54:03
7楼2011-04-25 10:02:52
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+2): 提问奖赏 2011-04-25 12:19:30
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+2): 提问奖赏 2011-04-25 12:19:30
|
这是复合物配体正确的坐标 HETATM 2623 C4 FER A 169 9.641 10.138 14.603 1.00 45.97 C HETATM 2624 C5 FER A 169 8.912 11.273 14.280 1.00 44.51 C HETATM 2625 C6 FER A 169 8.154 11.294 13.124 1.00 43.93 C HETATM 2626 C3 FER A 169 9.632 9.050 13.778 1.00 44.28 C HETATM 2627 C2 FER A 169 8.891 9.061 12.564 1.00 45.02 C HETATM 2628 C1 FER A 169 8.120 10.211 12.261 1.00 44.03 C HETATM 2629 C10 FER A 169 10.437 7.399 15.544 1.00 44.04 C HETATM 2630 C7 FER A 169 7.314 10.472 11.052 1.00 44.58 C HETATM 2631 C8 FER A 169 7.861 11.594 10.595 1.00 41.32 C HETATM 2632 C9 FER A 169 7.751 12.575 9.506 1.00 32.64 C HETATM 2633 O1 FER A 169 7.071 12.505 8.551 1.00 36.22 O HETATM 2634 O2 FER A 169 8.448 13.554 9.605 1.00 30.47 O HETATM 2635 O4 FER A 169 10.381 10.070 15.758 1.00 47.92 O HETATM 2636 O3 FER A 169 10.410 7.999 14.226 1.00 48.42 O 这是对接后配体坐标 HETATM 1 C 0 3.265 15.861 8.883 0.00 0.00 0.099 A HETATM 2 C 0 2.378 15.883 7.803 0.00 0.00 0.040 A HETATM 3 C 0 2.690 15.245 6.605 0.00 0.00 0.014 A HETATM 4 C 0 4.481 15.174 8.743 0.00 0.00 0.102 A HETATM 5 C 0 4.782 14.541 7.540 0.00 0.00 0.051 A HETATM 6 C 0 3.906 14.555 6.435 0.00 0.00 -0.033 A HETATM 7 O 0 2.997 16.490 10.083 0.00 0.00 -0.358 OA HETATM 8 H 0 3.803 16.440 10.620 0.00 0.00 0.218 HD HETATM 9 O 0 5.415 15.160 9.766 0.00 0.00 -0.353 OA HETATM 10 C 0 6.402 14.141 9.652 0.00 0.00 0.210 C HETATM 11 C 0 4.287 13.865 5.199 0.00 0.00 0.018 C HETATM 12 C 0 4.993 14.348 4.162 0.00 0.00 0.085 C HETATM 13 C 0 5.326 13.517 2.905 0.00 0.00 0.198 C HETATM 14 O 0 4.407 13.416 2.053 0.00 0.00 -0.645 OA HETATM 15 O 0 6.492 13.045 2.906 0.00 0.00 -0.645 OA 明显差别较大,我看了构象,对接后配体平躺着。这该怎么处理?请指教 |
8楼2011-04-25 10:29:26
9楼2011-04-25 10:37:32
10楼2011-04-25 10:56:52
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
12楼2011-04-25 11:02:47
13楼2011-04-25 11:04:21
14楼2011-04-25 11:18:55
15楼2011-04-25 11:52:19
16楼2011-04-25 12:14:02
★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:16:15
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:16:15
|
原帖位置http://muchong.com/bbs/viewthread.php?tid=3078961 不知道大家是否有做过像Gay-Berne potential那种跟分子方向和距离有关的site-site potential model 的模拟。对于这种非中心对称的两体作用势一般是怎么做截断的呢?又是如何做截断修正的呢? Lennard-Jones球对称势截断以后通过从截断半径到无穷远积分就可以了,可是如果分子不是球形的,比如是椭球状或者棒状,一般要怎么处理呢?? —————————————————— 目前我的处理方法和文献上看到的差不多,就是在4sigma长度处截断,但是没有做尾部修正,但是还是觉得这样做不是太准确,因为周期边界条件来实现体系无穷粒子的时候这个尾部势还是理论上不应忽略的。。 |
17楼2011-04-25 12:18:18
18楼2011-04-25 12:20:48
★ ★ ★
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:16:44
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:16:44
|
再在问一个。。NPT-MC模拟里面体积移动和位置移动的比例Frenkel和Smit在understanding Molecules Simulation 那本书里建议是1:N (N是格子中的分子数),他们说是体积移动的运算量通常是位置移动的N倍(非标度计算体积运动前后的总能量),因而为保证马尔柯夫链的对称性所以采取这个比例。 然而从我自己的测试看,即是比例是1:1对结果也没太明显的影响,所以我想请问这个比例是必然只能如此的吗,马尔柯夫链的对称性又是什么意思,印象中只听说过细致平衡的对称性。。 |
19楼2011-04-25 12:26:13
20楼2011-04-25 12:43:37
21楼2011-04-25 12:45:48
★ ★ ★
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:18:28
御剑江湖(金币+3): 提问奖赏 2011-04-25 15:18:28
22楼2011-04-25 12:48:05
23楼2011-04-25 12:51:31
24楼2011-04-25 13:38:04
25楼2011-04-25 14:07:18
26楼2011-04-25 14:14:08
27楼2011-04-25 14:25:31
28楼2011-04-25 14:47:13
29楼2011-04-25 14:48:19
|
分子模拟版块从今天起举行一个星期的“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”活动,在各位空闲之余期待各位去活动帖从大伙解答问题,除了解答后有金币重奖,活动结束后,我们将做出此活动专辑,形成一份本活动PDF文件,供虫子们下载,活动将由分子模拟版主委员会及群众投票评出”十大星虫“。 详见:http://muchong.com/bbs/viewthread.php?tid=3119249&fpage=1 |
30楼2011-04-25 15:38:09
31楼2011-04-25 15:39:15
32楼2011-04-25 16:04:47
33楼2011-04-25 16:32:26
34楼2011-04-25 17:36:41
35楼2011-04-25 17:43:37
36楼2011-04-25 18:08:09
37楼2011-04-25 18:15:15
★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 分享资源奖 2011-04-25 21:55:08
小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+3): 分享资源奖 2011-04-25 21:55:08
|
刚从网上搜到的,大家看看 CPMD:是一个Car-Parrinello 分子动力学程序包 很长一段时间以来,一直有人在问BOMD与CPMD究竟有什么不同?这里我就简单说一些二者的区别。 实际上,我想大部分人理解和接触的第一原理分子动力学方法以CPMD居多。CPMD,就是Car和Parrinello两个人作出的基于密度泛函的分子动力学方法,其特点是在引入电子虚拟质量,将电子运动耦合到了运动方程中,每一步分子动力学计算后,对电子结构的计算就不再需要自洽场迭代这一过程,因此可以大大节省计算资源,在计算机还不是那么好的80,90年代是弥足珍贵的。CPMD的建立开创了第一原理分子动力学方法的新时代,使其真正开始了实用化进程。 而BOMD,顾名思义,就是Born-Oppenheimer分子动力学。Born-Oppenheimer近似,也就是绝热近似,指的是将电子和离子的求解分离开来,只处理离子的动力学部分,而认为电子可以快速跟上电子的运动。其特点是每一步分子动力学计算之后都需要对电子结构进行自洽场迭代,使电子达到基态。正式由于这个自洽场迭代过程需要的计算量巨大,致使其一直没有达到广泛运用,直到90年代中后期,计算机技术的发展,才使BOMD开始逐步被人们重视。 CPMD和BOMD各有优劣。CPMD虽然计算速度更快,不需要进行自洽场迭代计算,计算量小,但是由于计算中并没有使体系真正达到基态,而只是尽量靠近基态,因此其准确性对电子的虚拟质量这个参数的选取依赖程度很大,一旦计算参数不对,得到的体系很可能远远偏离真实的势能面,得不到正确的动力学轨迹。同时,为了保证其尽量靠近基态,分子动力学时间步长一般选得较小。而BOMD虽然计算量大,但是由于每一步都保证系统达到基态,因此其分子动力学步长可以取得较大,一般1fs到5fs都有可能。综合来说,如果能结合二者的优势,第一原理分子动力学计算效率将大大提高,这也成为近年来其发展的重要方向。 官网:http://www.cpmd.org/ 简介: What is CPMD ? The CPMD code is a plane wave/pseudopotential implementation of Density Functional Theory, particularly designed for ab-initio molecular dynamics. Its first version was developed by Jurg Hutter at IBM Zurich Research Laboratory starting from the original Car-Parrinello codes. During the years many people from diverse organizations contributed to the development of the code and of its pseudopotential library: Michele Parrinello, Jurg Hutter, D. Marx, P. Focher, M. Tuckerman, W. Andreoni, A. Curioni, E. Fois, U. Roetlisberger, P. Giannozzi, T. Deutsch, A. Alavi, D. Sebastiani, A. Laio, J. VandeVondele, A. Seitsonen, S. Billeter and others. The current version, 3.13, is copyrighted jointly by IBM Corp and by Max Planck Institute, Stuttgart, and is distributed free of charge to non-profit organizations ( see download ). Profit organizations interested at the code should contact us . CPMD runs on many different computer architectures and it is well parallelized (MPI and Mixed MPI/SMP). CPMD main characteristics works with norm conserving or ultrasoft pseudopotentials LDA, LSD and the most popular gradient correction schemes; free energy density functional implementation isolated systems and system with periodic boundary conditions; k-points molecular and crystal symmetry wavefunction optimization: direct minimization and diagonalization geometry optimization: local optimization and simulated annealing molecular dynamics: constant energy, constant temperature and constant pressure path integral MD response functions excited states many electronic properties time-dependent DFT (excitations, molecular dynamics in excited states) coarse-grained non-Markovian metadynamics CPMD 是一个Car-Parrinello 分子动力学程序包。它也基于平面波赝势框架,支持模守恒和超软两种赝势。CPMD 软件中包括自由能密度泛函、路径积分分子动力学和meta-dynamics 等方法。它可以计算响应函数、激发态和其他许多电子性质。 |
38楼2011-04-25 18:25:09
39楼2011-04-25 19:36:46
40楼2011-04-25 21:11:44
41楼2011-04-25 21:27:58
42楼2011-04-25 21:36:55
43楼2011-04-25 21:38:47
44楼2011-04-25 21:50:27
★ ★ ★ ★ ★
zh1987hs(金币+5, 模拟EPI+1): 谢谢 2011-04-26 14:26:42
zh1987hs(金币+5, 模拟EPI+1): 谢谢 2011-04-26 14:26:42
|
要单独制作一个UserControl,这个控件里面是由多个Panel组合而成的。 在加载这个控件的窗口要有一个专门的控制组件,该组件记录边栏的状态,如隐藏、左停靠、右停靠、下停靠等。根据控件自己的点击事件来设置控件的Dock属性。 控件内部的分割,可用多个按钮加多个Panel进行设计,内部要根据VS的设计接口编辑Designer来实现所见即所得的编辑。 如果嫌麻烦,也可以直接使用DevExpress控件库,该库里面包含你所展示的窗口的所有功能,而且还支持换肤。可以省去很多不必要的时间。 |
45楼2011-04-25 22:01:08
46楼2011-04-25 23:19:36
47楼2011-04-26 09:22:38
48楼2011-04-26 09:48:18
49楼2011-04-26 10:35:13
50楼2011-04-26 11:23:19













回复此楼





