24小时热门版块排行榜    

查看: 2057  |  回复: 17

03010328

银虫 (正式写手)

[求助] MD进行分子模拟关于计算键角张力的问题 已有1人参与

我打算用分子动力学模拟基于DNA连接的纳米粒子的自组装,但是在构建初始化模型计算力场的时候,由于要考虑到DNA分子的刚性,所以在DNA链上必须要引入键角张力,其势能函数与弹簧振子的弹性势能一样,不过将距离换成键角,而通过粗粒化模型我们只知道相邻三个分子的三维空间坐标,现在问题来了,因为力是势能函数对位移的一阶负倒数,虽然可以通过余弦定理将键角张力势能函数中的键角转换成位移,可是如何将每个原子的受力求解出来,请神人赐教。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

梦のしずく
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
月只蓝: 金币+2, 感谢指导! 2013-11-22 15:31:27
与计算其他的力没有什么两样,对位移求导数
2楼2013-09-23 09:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

03010328

银虫 (正式写手)

引用回帖:
2楼: Originally posted by sciencejoy at 2013-09-23 09:56:55
与计算其他的力没有什么两样,对位移求导数

共价键中只有一个键,相应也只有一个位移,很容导出力;但是键角最终用位移表示出来后会含有三个位移,此时要求力的话,应该对哪一个位移求导呢?
梦のしずく
3楼2013-09-23 14:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fhh2626

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流! 2013-11-22 15:31:34
你自己编程计算?应该是将直角坐标转化为欧拉坐标之后对角度求导
4楼2013-09-23 15:35:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

03010328

银虫 (正式写手)

引用回帖:
4楼: Originally posted by fhh2626 at 2013-09-23 15:35:14
你自己编程计算?应该是将直角坐标转化为欧拉坐标之后对角度求导

因为没有源程序,老板让自己尝试着写,比如共价键、范德华力和LJ作用力对位移求导直可以导出来,但是键角张力是两个共价键的夹角的函数,即使用余弦定理通过三段位移表示出了键角,从而使键角张力变成了位移的函数,可是问题很明显啊,这时函数中出现了三个位移,我不知道对哪一个求导了?还望不吝赐教。
梦のしずく
5楼2013-09-23 17:05:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

03010328

银虫 (正式写手)

引用回帖:
4楼: Originally posted by fhh2626 at 2013-09-23 15:35:14
你自己编程计算?应该是将直角坐标转化为欧拉坐标之后对角度求导

就是这个问题整得头疼。
MD进行分子模拟关于计算键角张力的问题
1-3.jpg

梦のしずく
6楼2013-09-23 17:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

【答案】应助回帖


月只蓝: 金币+1, 感谢指导! 2013-11-22 15:31:46
引用回帖:
6楼: Originally posted by 03010328 at 2013-09-23 17:12:45
就是这个问题整得头疼。

1-3.jpg
...

求谁受力对谁求导
7楼2013-09-23 20:47:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

03010328

银虫 (正式写手)

引用回帖:
7楼: Originally posted by sciencejoy at 2013-09-23 20:47:02
求谁受力对谁求导...

我的第一想法是这样做,不过计算量会很大了。
MD进行分子模拟关于计算键角张力的问题-1
1-4.jpg

梦のしずく
8楼2013-09-23 21:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

引用回帖:
8楼: Originally posted by 03010328 at 2013-09-23 21:06:44
我的第一想法是这样做,不过计算量会很大了。

1-4.jpg
...

当然不可以
9楼2013-09-23 22:37:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

03010328

银虫 (正式写手)

引用回帖:
9楼: Originally posted by sciencejoy at 2013-09-23 22:37:59
当然不可以...

话说一半,那你觉得应该怎么求导,想不到了,网上也搜不到,多多指点。
梦のしずく
10楼2013-09-24 12:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 03010328 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +12 不知道叫什么! 2026-03-15 14/700 2026-03-16 16:45 by 我的船我的海
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 3/150 2026-03-16 16:44 by 我的船我的海
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
请填处理意见