关于MDAnalysis中Density对象产出内容及bins含义的技术问询

一、Density对象核心产出内容 #

• 核心为三维密度场数据：以numpy数组形式存储，对应空间中每个网格点的原子数密度或质量密度（由初始化DensityAnalysis时的density_type参数决定，默认是原子数密度）
• 附带关键元数据：
  • origin：网格的起始坐标（三维数组）
  • spacing：每个网格单元（bin）在三个轴向上的边长（三维数组）
  • extent：密度场覆盖的空间边界范围（六元组，格式为[x_min, x_max, y_min, y_max, z_min, z_max]）

二、(113, 113, 113) bins的含义 #

• 这是三维空间被划分的网格单元数量：x、y、z三个轴分别被分成113个连续的网格单元，总网格数为113×113×113
• 生成逻辑：DensityAnalysis会根据分析原子组的空间分布范围，结合指定的bin边长（或自动计算的默认值）来确定每个轴的bin数，计算公式近似为：

  单轴bin数 = 向上取整(该轴方向原子分布的总长度 / bin边长)
  

• 实例说明：若某轴方向原子的分布范围是22.6nm，默认bin边长为0.2nm，那么22.6 / 0.2 = 113，就会得到该轴113个bin的结果

三、Density对象的实用解读方法 #

• 提取数值数据：通过density.grid直接获取三维密度数组，可用于计算密度峰值、区域平均密度等量化分析
• 可视化展示：配合matplotlib的contourf（切面密度图）或plot_surface，或者Mayavi的三维渲染工具，直观呈现密度分布
• 单位转换：使用density.convert_units(target_unit)方法切换密度单位，比如从atoms/nm³转换为g/cm³

内容的提问来源于stack exchange，提问作者Jovey Osagie