24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2094  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京2,材料与化工308求调剂 +13 熊二想上岸 2026-04-04 14/700 2026-04-06 19:42 by guanxin1001
[考研] 08600生物与医药-327 +8 18755400796 2026-04-05 8/400 2026-04-06 11:59 by lijunpoly
[考研] 生物学求调剂 +5 15064154688 2026-04-03 5/250 2026-04-06 11:56 by lijunpoly
[考研] 315求调剂 +14 欣喜777 2026-04-04 15/750 2026-04-06 06:25 by houyaoxu
[考研] 求调剂 +11 熊二想上岸 2026-04-04 11/550 2026-04-05 22:21 by 醉翁wl
[考研] 材料专硕322分 +10 哈哈哈吼吼吼哈 2026-04-04 10/500 2026-04-05 21:22 by 学员8dgXkO
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 求调剂求调剂 +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 288求调剂 一志愿哈工大 材料与化工 +39 洛神哥哥 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
信息提示
请填处理意见