计算化学公社

 找回密码 Forget password
 注册 Register
Views: 685|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[吐槽] 越来越受不了高斯的cube文件了!!!我有更好方案

[复制链接 Copy URL]

592

帖子

12

威望

4799

eV
积分
5631

Level 6 (一方通行)

密度泛函·小卒

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 yjcmwgk 于 2024-12-28 16:37 编辑

用高斯十几年了,今天突然特别想一脚踹烂这个cube文件格式!
我希望的cube文件应该这样:
  1. U { A } U_END           !U就是Unit。这一行如果是A,表示本文件中坐标是埃为单位,如果是B,表示本文件中坐标是波尔为单位

  2. T { LOL } T_END         !T就是type,这个cube文件是干啥的。这里写的是LOL。当然还可以是EP(静电势)、D(电子密度)、IRI、NICS、S(也就是Spin)、……等等等等

  3. M { 3                   !M就是分子的意思,3表示分子里有三个原子
  4. x1  y1  z1  原子序数
  5. x2  y2  z2  原子序数
  6. x3  y3  z3  原子序数     !这里每一行都是四个数值,前三个是原子核所在的坐标,是双精度实数,最后一个是原子序数,是整数
  7. } M_END

  8. C {
  9.     a0  b0   c0
  10.     na  va1  va2  va3  
  11.     nb  vb1  vb2  vb3
  12.     nc  vc1  vc2  vc3   !这里的C就是cube的意思,a0 b0 c0是起始采样点。nx  ny  nz是三个整数,表示各个边上有多少个采样点。每一行有一个v表达的方向矢量
  13. } C_END

  14. D {
  15.     x0  y0  z0  num0
  16.     x1  y1  z1  num1
  17.     x2  y2  z2  num2
  18.     x3  y3  z3  num3
  19.     x4  y4  z4  num4
  20.     x5  y5  z5  num5
  21.     ......
  22.     一直往下写就是了,一个采样点对应一个数字,写吧    !na*nb*nc有多大,这里就有多少行!
  23. } D_END                 !这里的D就是data的意思

  24. N { 2                   !这里的N就是note也就是注释的意思,2就是有2行注释,你有几行注释就写几,注释内容你想写啥就写啥,不是给机器看的,是给人看的。
  25.    去你的Gaussian cube文件吧!
  26.    我这样更清晰得多!
  27. } N_END

复制代码


评分 Rate

参与人数
Participants 2
eV +9 收起 理由
Reason
Serious + 4 你太可爱
丁越 + 5 哈哈哈哈哈

查看全部评分 View all ratings

一出生响亮登场,十几岁快乐成长,
二十岁天天向上,三十岁基本定向,
四十岁拼命打创,五十岁回首一望,
六十岁告老还乡,七十岁搓搓麻将,
八十岁躺在床上,九十岁挂在墙上,
人生一世,匆匆忙忙,生得嘹亮,走得凄凉!
生活就像五味瓶,酸甜苦辣难消停!
该吃吃,该喝喝,遇事别去心里搁,
想哭哭,想笑笑,烦恼就往云外抛,
记住甜,忘掉苦,亲人朋友好相处,
心情好,最重要,自我麻痹乐逍遥!

326

帖子

3

威望

1752

eV
积分
2138

Level 5 (御坂)

2#
发表于 Post on 2025-3-15 23:09:41 | 只看该作者 Only view this author
好有趣的提议。确实有必要让文件格式这么标准化一下了。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 04:00 , Processed in 0.166056 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list