24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4881  |  回复: 16

15216105346

金虫 (著名写手)

[求助] 求计算键角、配位数的程序已有3人参与

本人做的工作是通过分子动力学模拟计算冰晶石Na3AlF6熔体(非晶体)中的O-F-O的键角,Al原子周围F的配位数。
目前我已得到平衡构像中所有原子的三维坐标,
现求助通过原子坐标数据来计算键角、配位数的程序!!!!
哪位大神有计算A-B-A型键角、配位数的程序,请发我一份!
我的邮箱:15216105346@163.com
不甚感激!
回复此楼

» 收录本帖的淘帖专辑推荐

动力学模拟

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
15216105346: 金币+40, ★★★很有帮助 2015-01-24 17:09:45
引用回帖:
4楼: Originally posted by 15216105346 at 2015-01-23 12:55:19
我的原子分布是没有周期性的,共1320个原子。
我想通过循环找到任意一个离Al原子最近的两个F原子坐标,这样计算这三个原子之间的距离a,b,c,再通过余弦定理求得夹角,
统计得到夹角的分布频率。
苦于不会编写 ...

建一循环,算出任意两个原子之间的距离,如果大于第一近邻的截断就不要,如果小于就保留为最近邻列表。
在每个原子的最近邻列表中,算出中心原子 与其它原子 的空间夹角,(数学问题,已知三角形三边长求夹角),最后算统计

不会编程会很麻烦,因为给你代码也需要你自己根据自己的数据格式进行调整。matlab会用不?

配位数可以通过lammps中的 compute  coord/atom来得到

» 本帖已获得的红花(最新10朵)

5楼2015-01-24 14:08:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbf3ng

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
配位数作个图就知道了,测量键角,pymol就可以,甚至自己可以算,三个点坐标算一个角的纯数学问题。

[ 发自小木虫客户端 ]
2楼2015-01-23 07:08:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congtaiupc

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
配位数可以通过分析RDF,然后再在Origin中作图,积分就是配位数。方法不只局限于一种~
用青春书写青春
3楼2015-01-23 09:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)

楼主您好,我也在学习熔体的分子动力学模拟,可以告诉我初始构象该怎么建吗?
6楼2015-09-14 15:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

15216105346

金虫 (著名写手)

引用回帖:
2楼: Originally posted by wbf3ng at 2015-01-23 07:08:21
配位数作个图就知道了,测量键角,pymol就可以,甚至自己可以算,三个点坐标算一个角的纯数学问题。

我的原子分布是没有周期性的,共1320个原子。
我想通过循环找到任意一个离Al原子最近的两个F原子坐标,这样计算这三个原子之间的距离a,b,c,再通过余弦定理求得夹角,
统计得到夹角的分布频率。
苦于不会编写程序。

配位数是怎么作图的啊!请详细介绍一下,大神!
4楼2015-01-23 12:55:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)

引用回帖:
6楼: Originally posted by 我是砖头 at 2015-09-14 15:41:13
楼主您好,我也在学习熔体的分子动力学模拟,可以告诉我初始构象该怎么建吗?

晶体cif文件导入,升温3000k融化,在降温到目标温度

[ 发自小木虫客户端 ]
7楼2015-09-14 22:55:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)

引用回帖:
7楼: Originally posted by 15216105346 at 2015-09-14 22:55:44
晶体cif文件导入,升温3000k融化,在降温到目标温度
...

谢谢
8楼2015-09-24 08:56:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjj318

金虫 (小有名气)

引用回帖:
7楼: Originally posted by 15216105346 at 2015-09-14 22:55:44
晶体cif文件导入,升温3000k融化,在降温到目标温度
...

cif文件是哪里来的,直接输入原子数目,势函数构建不可以吗
9楼2015-10-11 16:01:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elmer123

禁虫 (初入文坛)

本帖内容被屏蔽

10楼2016-04-25 15:50:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 15216105346 的主题更新
信息提示
请填处理意见