24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 354  |  回复: 0

xzz1007

金虫 (小有名气)

[求助] 求助:关于内坐标优化结构时,梯度转化的问题

目前,我所研究的内容是:使用内坐标进行结构的优化。
已知:笛卡尔坐标系下的坐标(x,y,z)以及力(Fx,Fy,Fz)
如果使用内坐标优化的话,首先需要将笛卡尔坐标系下的坐标(x,y,z)以及力(Fx,Fy,Fz)
转化成内坐标的形式, 目前,坐标已经转化完毕了,但是相对应的力,不知道应该如何转化??
下面以甲烷为例,进一步说一下我的问题所在:

已知:CH4的笛卡尔坐标以及力

Cartesian coordinates:
---------------------------------------------------------------------
Center         Atomic                Coordinates (Angstroms)
Number      Number                   X          Y              Z
---------------------------------------------------------------------
1                  6                 0.000000   0.000000   0.000000
2                  1                 0.625276   0.625276   0.625276
3                  1                 -0.625276   -0.625276   0.625276
4                  1                  -0.625276   0.625276   -0.625276
5                  1                   0.625276   -0.625276   -0.625276
---------------------------------------------------------------------

Force::

------------------------------------------------------------------------------------------
Center         Atomic                                 Forces (Hartrees/Bohr)
Number       Number                   X                      Y                    Z
-------------------------------------------------------------------------------------------
1                   6              0.000000000   0.000000000    0.000000000
2                   1              0.000000000   0.000000000    0.000000425
3                   1               0.000000401   0.000000000    -0.000000142
4                   1              -0.000000200   -0.000000347   -0.000000142
5                   1              -0.000000200   0.000000347    -0.000000142
-------------------------------------------------------------------

我使用tinker中的一个子程序可以将笛卡尔坐标转化成内坐标
internal coordinates:

-------------------------------------------------------------------
            bond              angle                 tors
-------------------------------------------------------------------
1
2   1   1.08300
3   1   1.08300   2   109.4712
4   1   1.08300   2   109.4712   3   120.0000
5   1   1.08300   2   109.4712   3   -120.0000

-------------------------------------------------------------------
下面,我就不知道如何将力,也就是负的梯度,也传化成内坐标的形式呢??

[ 来自小组 减肥吧 ]
回复此楼
阳光总在风雨后,再坚持坚持吧.....
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 xzz1007 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 287求调剂 +7 land xuxu 2026-03-26 7/350 2026-03-27 00:44 by wxiongid
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 311求调剂 +5 lin0039 2026-03-26 5/250 2026-03-26 22:43 by 学员8dgXkO
[考研] 281求调剂 +3 亚克西good 2026-03-26 5/250 2026-03-26 19:48 by 不吃魚的貓
[考研] 279 分 求调剂 +3 睡个好觉_16 2026-03-24 3/150 2026-03-26 15:43 by zzll406
[考研] 296求调剂 +5 www_q 2026-03-20 5/250 2026-03-26 12:56 by 3Strings
[考研] 321求调剂 +3 璞玉~~ 2026-03-25 3/150 2026-03-25 19:07 by Zhanglab-TJU
[考研] 26考研-291分-厦门大学(085601)-柔性电子学院材料工程专业求调剂 +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 材料学硕,求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
信息提示
请填处理意见