24小时热门版块排行榜    

查看: 5518  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 我不理解! +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见