24小时热门版块排行榜    

查看: 5597  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:03 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:03 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-16 04:23 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 03:02 by DgNGHc3h5tPl
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 02:02 by DgNGHc3h5tPl
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
信息提示
请填处理意见