计算化学公社

标题: SDF/空间分布函数 [打印本页]

作者
Author: tjuchan 时间: 2016-5-15 22:04
标题: SDF/空间分布函数
采用MC跑出的结构，想分析分子周围某个原子的空间分布函数。由于结构文件（*.pdb）的格式问题，于是想自己编写计算SDF的程序。1、读取每个点坐标，固定不共线三点作为参考点，将它的坐标移动至原点，x轴和xy平面上。作为中心点；如图1（DMF分子，因为有二面角项，所以不可能所有分子都完美重合）
2、其它点随之变化坐标，采用周期边界将所有点转移到盒子内；如图2
问题是如何计算空间原子的密度呢？
已知空间各点的坐标，如何转化为cube文件呢？

作者
Author: sobereva 时间: 2016-5-16 11:23
就把空间划分为大小相同的格子，统计每个格子内在整个轨迹中出现的粒子数目，然后除以格子体积和帧数，就行了。cube文件的格式参照
Gaussian型cube文件简介及读、写方法和简单应用
http://sobereva.com/125

欢迎光临计算化学公社 (http://bbs.keinsci.com/)