注册 | 登录
|||
空间分布函数SDF的计算及三维图示
2014–05–13 17:52:00 硕士生李耀对此文有所贡献
基本概念
分子动力学模拟中常以径向分布函数(radical distribution function, RDF)来描述粒子周围环境的分布特性, 并表征其短程有序性. 分布函数的概念可推广到三维$,$ 三维空间的分布函数被称为空间分布函数(spatial distribution function, SDF), 可定义为与粒子的局部数密度与平均数密度之比 Ω(x,y,z)=ρ(x,y,z)ρˉ, 数学实质是三维的某种密度分布.
一般来讲, 三维密度函数的图示方法有好几种, 其中主要的几种在量子化学的轨道图示中都有使用过, 如点云图, 径向分布图, 角度分布图, 等值面图, 切面图. 此外还有几种图示方法是计算流体力学中常用的, 在化学中使用较少.
计算方法
SDF计算过程与具体体系相关, 不易统一, Gromacs提供了一个g_spatial的工具, 但使用不是很方便. 很多人使用gOpenMol, 如 Distrubutions of Counterions around DNA. Molecular Dynamics Simulations Results. 我没有使用过gOpenMol, 无法评论.
下面是SDF的计算方法
中心分子: 确定要分析的目标分子
PBC平移: 将每一个目标分子平移至原点作为中心, 根据周期性边界条件(periodic boundary condition, PBC)对其他分子进行相应平移
最佳叠合/旋转一致: 对PBC平移后的构型进行最佳叠合或旋转, 使中心分子取向一致
统计密度: 利用中心分子取向一致后的构型统计其周围粒子的密度分布
几点说明
PBC平移的方法可参考我的博文 取整截断函数及其在PBC中的使用.
若每帧轨迹中含有N个目标分子, 则每帧轨迹最多可得到N个构型. 使用每个目标分子作为中心分子是为了充分利用每帧的数据, 如果你的轨迹文件时间很长, 每帧只使用一个中心分子也可以
最佳叠合/旋转一致过程我以前也简略 说过. 如果中心分子在MD过程中构型变化不大, 或其关键部分是刚性的, 变化不大, 那么只要选择中心分子上不共线的三个原子, 第一个原子作原点, 第二个原子处于x轴正向, 第三个原子处于xy平面, z轴垂直于三个原子所在的平面, 将所有PBC平移后的构型都旋转至标准取向即可. 严格来说, 应该对中心分子做最佳叠合, 使所有中心分子相对于参考中心分子的位置的偏差(RMSD)最小. 这不是一个简单的工作, 但可利用四元数方法解决.
多个分子的叠合显示可利用Jmol或VMD. VMD支持的多帧XYZ文件, 原子类型及数目必须固定, 有时使用不便; Jmol支持的XYZ文件更宽泛, 但无法操控太大的文件.
VMD处理大文件的能力很赞, 一百万个原子轻松, 再多点稍困难. 我最大操作四百万个原子, 比较卡了.
统计密度时, 一般是使用直角坐标, 将整个空间分为若干均匀的小格子, 然后统计每个小格子中出现的原子个数, 再除以小格子的体积, 平均, 得到原子的数密度, 将此密度与平均密度相比, 得到SDF.
统计密度时, 可只统计中心分子附近某一区域内的(可由RDF获知), 因离中心分子很远的地方SDF基本趋于均匀, 意义不大
一般需要对多帧轨迹进行计算统计, 直至算出的SDF数值收敛.
一般将SDF保存为cube文件. cube文件格式始于Gaussian, 但现在已基本成为标准格式了. 具体格式可参考Gaussian手册或卢天的博文 Gaussian型cube文件简介及读、写方法和简单应用.
cube文件的可视化程序很多, 大凡好点的分子可视化程序都支持, 常用的如GaussView, VMD, Jmol, ChemCraft, Chemirea, gOpenmol. 如果要同时显示多个等值面, GaussView可能不是一个好的选择, VMD可能最简单, 具体方法可参考 VMD 画cube文件等密度面图(轨道或电子密度等)
代码
CntAtm
RotMol
SDF
具体示例
以一帧含5000个水分子的构型为例, 计算水分子周围O原子的SDF.
目标分子, 水分子, 此帧5000个
以每个水分子为中心, 其余水分子PBC平移至[-L/2, L/2]范围内, 共得5000个构型. 其中一个如下
将中心水分子进行叠加显示, 可见最佳叠合/旋转一致前, 由于取向随机, 中心O原子周围H原子基本成球状分布,
中心水分子周围其他水分子分布也呈均匀分布
将中心水分子旋转一致后, 中心O原子周围H原子分布与单个水分子类似, 但由于水分子构型在MD过程中有变化, 不会完全重合
旋转一致后, 中心水分子周围其他水分子的分布不再均匀, 显示出各向异性, 这种各向异性从叠合后的构型图也可看出,
O原子分布
H原子分布
水分子分布
计算cube文件, 利用VMD显示等值面, 两个不同值的等值面(蓝色为5, 红色为10)如下, 水分子周围H氢键的四面体结构很明显
Archiver| 手机版| 科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-6-16 10:55
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社
深圳SEO优化公司潜江如何制作网站多少钱沈阳至尊标王荆门外贸网站设计推荐辽阳设计公司网站价格大连seo优化多少钱海东seo排名报价德州企业网站改版公司唐山至尊标王哪家好揭阳seo网站推广哪家好张北关键词按天计费推荐包头网站推广工具公司株洲品牌网站设计公司宁德网站定制公司那曲建网站哪家好晋城SEO按效果付费哪家好临夏网站优化按天扣费哪家好湛江网络推广推荐桂林如何制作网站哪家好威海百度网站优化排名哪家好渭南品牌网站设计保定企业网站制作公司思茅网站建设价格曲靖企业网站建设报价安阳网络营销报价莆田百搜标王公司自贡百度竞价多少钱拉萨模板制作推荐邵阳阿里店铺托管多少钱曲靖设计公司网站多少钱襄阳网站优化软件价格歼20紧急升空逼退外机英媒称团队夜以继日筹划王妃复出草木蔓发 春山在望成都发生巨响 当地回应60岁老人炒菠菜未焯水致肾病恶化男子涉嫌走私被判11年却一天牢没坐劳斯莱斯右转逼停直行车网传落水者说“没让你救”系谣言广东通报13岁男孩性侵女童不予立案贵州小伙回应在美国卖三蹦子火了淀粉肠小王子日销售额涨超10倍有个姐真把千机伞做出来了近3万元金手镯仅含足金十克呼北高速交通事故已致14人死亡杨洋拄拐现身医院国产伟哥去年销售近13亿男子给前妻转账 现任妻子起诉要回新基金只募集到26元还是员工自购男孩疑遭霸凌 家长讨说法被踢出群充个话费竟沦为间接洗钱工具新的一天从800个哈欠开始单亲妈妈陷入热恋 14岁儿子报警#春分立蛋大挑战#中国投资客涌入日本东京买房两大学生合买彩票中奖一人不认账新加坡主帅:唯一目标击败中国队月嫂回应掌掴婴儿是在赶虫子19岁小伙救下5人后溺亡 多方发声清明节放假3天调休1天张家界的山上“长”满了韩国人?开封王婆为何火了主播靠辱骂母亲走红被批捕封号代拍被何赛飞拿着魔杖追着打阿根廷将发行1万与2万面值的纸币库克现身上海为江西彩礼“减负”的“试婚人”因自嘲式简历走红的教授更新简介殡仪馆花卉高于市场价3倍还重复用网友称在豆瓣酱里吃出老鼠头315晚会后胖东来又人满为患了网友建议重庆地铁不准乘客携带菜筐特朗普谈“凯特王妃P图照”罗斯否认插足凯特王妃婚姻青海通报栏杆断裂小学生跌落住进ICU恒大被罚41.75亿到底怎么缴湖南一县政协主席疑涉刑案被控制茶百道就改标签日期致歉王树国3次鞠躬告别西交大师生张立群任西安交通大学校长杨倩无缘巴黎奥运