24小时热门版块排行榜

返回列表

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

zhy0512

新虫 (初入文坛)

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

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

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

发自小木虫Android客户端

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

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

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

jerkwin

专家顾问 (正式写手)

专家经验: +14
应助: 454 (硕士)
金币: 20699.1
帖子: 813
在线: 2648.3小时
虫号: 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