交互式晶体结构三维可视化与分析工具

资 源 简 介

该项目致力于构建一个功能全面的MATLAB应用程序，用于读取、解析并在三维空间中渲染复杂的晶体结构。核心功能模块包括：1. 多格式数据导入引擎，支持解析标准的晶体信息文件（CIF）以及计算材料学常用的POSCAR格式，能够自动提取晶格常数、空间群信息及原子占位坐标；2. 坐标变换与扩胞算法，实现从分数坐标到笛卡尔坐标的精确转换，支持基于周期性边界条件的超胞（Supercell）构建，以及根据对称性操作生成完整的单胞原子排布；3. 高级三维渲染引擎，利用MATLAB的patch、surface及光照模型（lighting models），以球棍模型（Ball-and-Stick）或空间填充模型（Space-filling）逼真展示原子和化学键，根据元素周期表自动匹配原子颜色与半径；4. 交互式分析工具，用户可在GUI界面中实时旋转、缩放视图，动态调整成键阈值以构建或断开化学键，测量选定原子间的键长与键角，并支持导出高分辨率的矢量图或位图用于学术发表。该工具旨在为材料科学、化学及物理研究人员提供一个轻量级且高度可定制的结构分析与绘图解决方案。

详 情 说 明

MATLAB交互式晶体结构三维可视化平台

项目简介

该项目是一个基于MATLAB构建的轻量级晶体结构分析与可视化应用程序。它利用MATLAB App Designer的图形组件（uifigure）和强大的三维绘图引擎，为用户提供了一个直观的界面来加载、查看和分析晶体结构。程序支持从标准结构文件导入数据，构建超胞，并在三维空间中以球棍模型渲染原子和化学键，同时提供了交互式的键长测量和高分辨率图片导出功能。

功能特性

• 多源数据支持：程序内置了文件导入逻辑，能够区分并处理CIF（晶体信息文件）和POSCAR（VASP结构文件）格式，同时也包含一个内置的钙钛矿（CaTiO3）示例数据用于快速演示。
• 交互式三维渲染：

* 使用光照模型（Gouraud shading）和材质设置（Dull material）创建具有立体感的原子显示。 * 支持鼠标操作进行三维场景的旋转、缩放和平移。

• 超胞构建（Supercell）：用户可以通过界面动态设置a、b、c三个方向的扩胞倍数（1x到10x），程序会自动计算扩胞后的原子坐标并刷新视图。
• 化学键分析：

* 提供动态滑块调整成键阈值（Bond Cutoff），实时断开 or 生成化学键。 * 性能保护机制：为防止卡顿，代码中包含逻辑判断，仅在原子总数小于500时自动计算并绘制化学键。

• 交互式测量：实现了点击原子事件响应。用户依次点击两个原子后，会在信息面板显示选中原子的元素类型、坐标以及它们之间的欧氏距离（键长）。
• 图形导出：集成exportgraphics功能，支持将当前三维视图导出为高分辨率的PNG、JPG位图或EPS矢量图。

系统要求

• MATLAB R2020a 或更高版本（推荐使用最新版本以获得最佳的uifigure和uiaxes性能）。
• 无需额外的工具箱，但部分高级文件解析可能依赖基础数学库。

使用方法

1. 启动程序：在MATLAB命令行窗口直接输入 main 并回车即可启动应用，无需额外参数。
2. 加载数据：

* 点击“加载示例(CaTiO3)”可直接查看内置结构。 * 点击“导入结构文件”可选择本地的CIF或POSCAR文件。

1. 视图操作：

* 在左侧绘图区按住鼠标左键旋转，右键平移，滚轮缩放。 * 右侧“视图控制”面板可调整键长阈值或隐藏化学键。

1. 构建超胞：在“超胞构建”面板输入各方向倍数，点击“应用扩胞”。
2. 测量距离：直接在3D视图中点击任意原子（原子会高亮或在日志中反馈），连续点击两个原子后，右侧信息区将显示距离。
3. 保存图片：点击“导出高分辨率图片”并选择保存路径。

代码实现与逻辑分析

该项目完全包含在 main.m 文件中，采用单一入口函数配合嵌套函数的架构。通过共享的 app 结构体在不同回调函数间传递数据（app.Data）和状态（app.State）。

1. 界面架构与初始化

程序摒弃了传统的GUIDE，采用现代的 uifigure 和 uigridlayout 进行布局。

• 布局逻辑：界面被划分为左右两栏。左侧放置 uiaxes 用于三维绘图，关闭了默认的坐标轴网格，并配置了 infinite 光源；右侧放置控制面板，包含文件、扩胞、视图、信息和导出五个功能区。
• 交互定义：显式定义了 uiaxes 的交互模式（旋转、缩放、平移），并禁用了默认的数据提示（DataTip），以便接管鼠标点击事件用于自定义的距离测量。

2. 数据处理与坐标变换

核心数据存储在 app.Data 中，包含晶格矩阵（Lattice）和原子信息表（Element, Frac）。

• 超胞算法：在 renderScene 函数中实现。通过三重循环遍历 0 到 SuperCell-1 的整数位移，将位移量加到原始的分数坐标上。
• 坐标转换：程序严格遵循晶体学定义，将扩胞后的分数坐标乘以晶格矩阵转换为笛卡尔坐标，公式逻辑为：Cartesian = Fractional * Lattice。
• 渲染数据缓存：为了支持后续的交互点击，计算出的所有原子笛卡尔坐标和元素类型会被缓存到 app.Gui.CurrentRenderData 中。

3. 三维渲染引擎

渲染逻辑集中在 renderScene 函数，采用“清空-重绘”策略。

• 原子绘制：利用 sphere 函数生成球体网格数据。代码实现了基于原子类型的批量处理：先获取唯一元素列表，针对每种元素获取对应的半径和颜色，然后利用 surface 函数在指定坐标处绘制球体。为了实现交互，每个原子的 surface 对象都绑定了 ButtonDownFcn 回调，并将索引存储在 Tag 属性中。
• 化学键绘制：

* 距离计算：使用 pdist 函数计算所有原子间的欧氏距离，通过 squareform 转换为矩阵。 * 阈值筛选：根据用户设定的 BondCutoff 筛选出符合成键条件的原子对。 * 几何生成：对于每一根键，代码通过计算两点间的向量差，结合 cylinder 函数生成的圆柱网格，通过旋转矩阵变换将标准圆柱对齐到两个原子之间。

4. 交互式测量逻辑

• 点击识别：当用户点击原子时，触发 onAtomClick。函数通过读取被点击对象的 Tag 属性获取原子在数据表中的索引。
• 状态管理：使用 app.State.SelectedAtoms 维护当前选中的原子队列。
• 距离计算：当选中原子数量达到2个时，利用 norm 函数计算两个笛卡尔坐标向量差的模（即距离），并将结果输出到文本区域。

5. 文件解析与异常处理

• 文件流向：onLoadFile 根据文件后缀名（.cif 或 .vasp/.POSCAR）分流到相应的解析函数（parseCIF 或 parsePOSCAR）。
• 容错机制：包含 try-catch 块，若解析过程出错或解析后原子数据为空，会弹出 uialert 警告窗口提示用户，防止程序崩溃。