24小时热门版块排行榜    

查看: 1978  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

03010328

银虫 (正式写手)

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

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

fhh2626

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
月只蓝: 金币+1, 应助指数+1, 感谢指导! 2013-11-22 15:32:13
03010328: 金币+3, 有帮助, 谢谢! 2013-11-25 09:19:39
而且为什么要自己编程呢,你要编出个可以运行的MD程序没个一年半载的根本不大可能

你说的考虑DNA链的角张力,用NAMD中的user defined force就可以实现啊
14楼2013-09-24 17:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

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的回帖
信息提示
请填处理意见