24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhy0512

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 355.6
帖子: 12
在线: 11.2小时
虫号: 4029748

[交流] 对相关势函数求原子受力

已知各原子的三维坐标，对二体势（例如Morse势）求原子受力，可直接用势函数对位移求导。但三体势（例如键角能）如何求个原子受力？
查了好久文献资料都没找到相关内容，还希望各位前辈指教。
若有相关的文献和书籍，还望推荐。万分感谢。

发自小木虫Android客户端

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有6人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有12人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有16人回复
初秋的晨风已经有3人回复
面上没中，邀请各位路过的虫友分析一下分数已经有11人回复
Ei源刊怎么投已经有4人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复

1楼 2015-11-23 11:57:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

专家经验: +14
应助: 454 (硕士)
金币: 20659.1
帖子: 813
在线: 2648.7小时
虫号: 1023452

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

25楼: Originally posted by zhy0512 at 2015-11-24 13:33:14
谢谢。还想请教一下。
公式推导是将cos用余弦公式展开，然后再用链式求导吗？

...

仔细看看下面的注释吧
!  i<---j--->k Tht = Ang(i j k)
!  E = E(Tht), Rji*Rjk*cosTht = Dot(rji, rjk) = Xi*Xk+Yi*Yk+Zi*Zk
!  Rjk*gradi[Rji*cosTht] = gradi[Xi*Xk+Yi*Yk+Zi*Zk]
!  cosTht*gradi[Rji]-Rji*sinTht*gradi[Tht] = ^rjk
! -Rji*sinTht*gradi[Tht] = ^rjk - cosTht ^rji
!  gradi[Tht] = -(^rjk - cosTht ^rji)/(Rji*sinTht)

!  Fi = -gradi E(Tht) = -{ dE/dXi, dE/dYi, dE/dZi }
!    = -dE/dTht gradi[Tht]
!    = dE/dTht (^rjk - cosTht ^rji)/(Rji*sinTht)
!  Fk = dE/dTht (^rji - cosTht ^rjk)/(Rjk*sinTht)
!  Fj = -Fi-Fk