24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2097  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料考研求调剂总分280 +12 mkjlz1 2026-04-06 14/700 2026-04-06 22:48 by mkjlz1
[考研] 材料工程调剂 +3 小刘同学吖吖 2026-04-06 4/200 2026-04-06 20:56 by lbsjt
[考研] 一志愿国科大信工所,英二数二408总分293分求调剂 +3 ilcyuan 2026-04-02 4/200 2026-04-06 16:35 by likeihood
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 材料334求调剂 +19 Eecho# 2026-04-03 19/950 2026-04-06 08:37 by 小小树2024
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 电子信息调剂交叉学科有推荐吗 +6 jhtfeybgj 2026-04-01 9/450 2026-04-05 11:13 by 猪会飞
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 求调剂:085600材料与化工,考材科基,总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4 爱读书的小鳄鱼 2026-04-02 4/200 2026-04-02 18:35 by 求调剂zz
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 085600,321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
信息提示
请填处理意见