| 查看: 724 | 回复: 7 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请问如何计算LJ团簇的稳定构型？
|
|||
|
各位前辈好： 小弟分子动力学期末作业的题目是用共轭梯度法计算LJ团簇的稳定构型，但是LJ势只与两两分子间的距离有关，共轭梯度法只能得到势能最低时两两分子间的距离，（1）如何得到构型（每个分子的坐标）呢？。（2）计算初值也不应该是随便设的吧……势能Local minimum比较多，（3）如何才能找到Global minimum呢？小弟初学，希望各位前辈不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有12人回复
-
359求调剂
已经有7人回复
-
070300化学专业279调剂
已经有10人回复
-
环境工程297分求调剂一志愿杭高院
已经有9人回复
-
一志愿武理材料工程302调剂环化或化工
已经有9人回复
-
309分085801求调剂
已经有5人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有5人回复
-
材料调剂
已经有4人回复
-
材料调剂
已经有3人回复
-
311（085601）求调剂
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
 |
2楼2009-12-20 22:50:41
3楼2009-12-21 01:53:44
4楼2009-12-21 08:43:02
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+6,VIP+0):呵呵 够有耐心的 12-21 10:43
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+6,VIP+0):呵呵 够有耐心的 12-21 10:43
|
你已经有了两两分子之间的距离,这是比原子位置更详尽的结构信息,又何必再去找他们的空间位置呢? 知道了空间位置, 又有什么用处. 如果你真的要还原, 我不知道有没有简明的算法, 但我这里倒有个笨方法: 1. 令一号原子位于(0,0,0) 2. 令二号原子位于(x2,0,0),x2是1号原子和2号原子之间的距离. 3. 令其它原子的位置在(xi,yi,zi), 这就有3(N-2)=3N-6个未知数.N是原子个数. 4. 利用你的已知条件,建立方程组, 如第i号原子: (xj-xi)^2+(yj-yi)^2+(zj-zi)^2=rji^2 (j=i+1,N, i 不等于j) 你这样可以建立(1+N*N)/2个方程. 由(1+N*N)/2>3N-6可知,N<3或N>4的情况下,你可以解出所有原子位置. 现在假如你有5个原子,按上述方法,你首先按上述1,2步,定出1号和2号原子坐标,这时候坐标系也就给定了. 第三步: 对第3号原子有: x3*x3+y3*Y3+z3*z3=r13*r13 ,也就是说,3号原子位于以一号原子为中心,以r13为半径的圆周上. 同理,3号原子也位于以2号原子为中心,以r23为中心的圆周上, 即(x3-x2)^2+(y3-y2)^2+(z3-z2)^2=r23*r23, 注意y2=z2=0. 那么, 3号原子一定位于这两个圆周的交点,画图可知,有两个解. 第四步: 对第4号原子重复上述过程, 它的位置也有两个. 现在即使加上3号和4号原子的距离信息,也不能确定它们各自的位置(因为前面已经证明, 只有N>4才可能定出所有原子位置).由方程个数也可看出,现在有5个方程(r13,r23,r14,r24,r34), 却有六个未知数(x3,y3,z3,x4,y4,z4) 第五步:现在考虑第五个原子,重复上述过程, 你可以多建好多个方程(r15,r25,r35,r45....) 却只多出三个未知数,你现在可以解出1至5号原子的位置了..... 然后,你再固定1,2号原子不变, 另选三个原子当作3,4,5号原子,再解出它们的位置...(当然你也可以一股脑解出所有原子的位置), 前提是你的各原子之间的两两距离是合理的,方程是有解的. 最后你按所有原子位置,再算一下两两距离,看跟你的两两距离一样否? 这个问题,应该也可以用最优化的方法去求解，我就不再说下去了。打字太烦了。 [ Last edited by 老虎大王 on 2009-12-21 at 11:43 ] |
5楼2009-12-21 10:20:18
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
6楼2009-12-21 10:34:23
7楼2009-12-22 07:35:48
漂泊四方
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1321.2
- 散金: 50
- 红花: 6
- 帖子: 198
- 在线: 115.2小时
- 虫号: 897145
- 注册: 2009-11-08
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 很好的资源 12-28 20:21
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 很好的资源 12-28 20:21
|
你的目的是找LJ团簇的全局优化结构吧？按照我的理解，共轭梯度法只是一种局域优化的算法？所以只用这个算法找不到Global minimum. 有很多全局搜索的算法可以提供给你使用。比如SA，GA，AIOA，basin-hopping 等。 由于LJ是一个最简单的团簇，简单暴力的办法可能更直接些。你可以完全随机生成新结构，然后用你的共轭梯度算法进行局域优化，而忽略LJ团簇结构的特性。当随机次数足够的时候，你可以认为你的你的结构能量是全局最优。实际上是没有办法验证你找到的结构时候全局最优的，最好的方法是和前人的结果做比较。 就LJ 团簇， 你可以参照下我师兄以前做的结果和方法。 http://staff.ustc.edu.cn/~clj/ 他主页提供代码 ^_^。 |
8楼2009-12-28 14:23:19
10