计算化学公社

标题: 越来越受不了高斯的cube文件了！！！我有更好方案 [打印本页]

---

作者
Author: yjcmwgk    时间: 2024-12-28 15:20
标题: 越来越受不了高斯的cube文件了！！！我有更好方案


用高斯十几年了，今天突然特别想一脚踹烂这个cube文件格式！
我希望的cube文件应该这样：

1. U { A } U_END           !U就是Unit。这一行如果是A，表示本文件中坐标是埃为单位，如果是B，表示本文件中坐标是波尔为单位
   
2. 
   
3. T { LOL } T_END         !T就是type，这个cube文件是干啥的。这里写的是LOL。当然还可以是EP（静电势）、D（电子密度）、IRI、NICS、S（也就是Spin）、……等等等等
   
4. 
   
5. M { 3                   !M就是分子的意思，3表示分子里有三个原子
   
6. x1  y1  z1  原子序数
   
7. x2  y2  z2  原子序数
   
8. x3  y3  z3  原子序数     !这里每一行都是四个数值，前三个是原子核所在的坐标，是双精度实数，最后一个是原子序数，是整数
   
9. } M_END
   
10. 
    
11. C {
    
12.     a0  b0   c0
    
13.     na  va1  va2  va3  
    
14.     nb  vb1  vb2  vb3
    
15.     nc  vc1  vc2  vc3   !这里的C就是cube的意思，a0 b0 c0是起始采样点。nx  ny  nz是三个整数，表示各个边上有多少个采样点。每一行有一个v表达的方向矢量
    
16. } C_END
    
17. 
    
18. D {
    
19.     x0  y0  z0  num0
    
20.     x1  y1  z1  num1
    
21.     x2  y2  z2  num2
    
22.     x3  y3  z3  num3
    
23.     x4  y4  z4  num4
    
24.     x5  y5  z5  num5
    
25.     ......
    
26.     一直往下写就是了，一个采样点对应一个数字，写吧    ！na*nb*nc有多大，这里就有多少行！
    
27. } D_END                 !这里的D就是data的意思
    
28. 
    
29. N { 2                   !这里的N就是note也就是注释的意思，2就是有2行注释，你有几行注释就写几，注释内容你想写啥就写啥，不是给机器看的，是给人看的。
    
30.    去你的Gaussian cube文件吧！
    
31.    我这样更清晰得多！
    
32. } N_END
    
33. 
    

复制代码

---

作者
Author: Stardust0831    时间: 2025-3-15 23:09
好有趣的提议。确实有必要让文件格式这么标准化一下了。

---

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

Powered by Discuz! X3.3