自编程序

SGROUP
功能
 
用于判断或指定周期体系的空间群,包括一维棒群,二维层群,三维空间群。对于三维体系,产生的数据文件可以用通过Matlab显示Wigner-Seitz晶胞和(约化)布里渊区。
说明见使用手册。
下载
 
可执行程序(Windows),使用手册,和示例
(版本号:3.0.06)
运行方法双击work目录下的run.bat。细节见使用手册。
示例 
 例1:一氧化碳吸附于MgO表面。使用惯用晶胞(ifprim = .f.)。
注意:FORTRAN语言要求“$”之前一定要有一个以上的空格,否则无法读取输入数据。
# CO on MgO (001) two layer slab: one-side adsorption

 $lattice
 centring = 'P'
    idimn =   2

       aa =   2.97692
       bb =   2.97692
    alpha =  90.0
     beta =  90.0
    gamma =  90.0

   lprint =  -1
   ifprim = .f.
 $end
O   0.00   0.00   4.5625
C   0.00   0.00   3.4125
Mg  0.00   0.00   1.0525
O   0.50   0.50   1.0525
Mg  0.00   0.00  -1.0525
O   0.50   0.50  -1.0525
 

计算结果:test1.txt
这是55号层群,对应99号空间群。

例2:指定空间群(ITNO = 225;面心立方),最简输出,并用Matlab绘制布里渊区和约化布里渊区(使用程序产生的数据文件WS3D.M和IBZ3D.M)。

# cubic fcc lattice (ITNO=225)
 $lattice
    IDIMN = 3
     ITNO = 225

       aa = 2.0

    iws3d = 0,0,6,1,0
    ibz3d = 3,10

 $end
X 0 0 0
 

布里渊区和约化布里渊区如图:

升级 
 (暂无)

 

 

Level Fitting (LF)
功能LEVEL产生的数据文件计算双原子分子的光谱常数。已对LEVEL 7.4 - 8.0进行了测试。
说明




 
1. 程序需要LEVEL的输出文件fort.7和fort.8。拟合转动常数还需要fort.9。为此,需要在LEVEL输入文件定义LPPOT=-1,LXPCT=-1,LCDC=2。
2. 可以拟合的光谱常数有:Re,Te,We,WeXe,De,和Be。
3. 由于较高的振动能级非谐性较强,因此振动能级最多只读取前11个(Vmax=10)。
4. 程序产生的FORT.DAT文件可直接用于Origin 6.x/7.x绘图。
   方法:
   打开Origin。
   从菜单选择File/Import/Single ASCII,找到FORT.DAT文件,打开并绘图。
注意事项
 
LEVEL和Level Fitting能够正常结束,既不表明计算结果正确,也不意味着输入的参数合适。用户应仔细阅读LEVEL的使用手册。建议多调试RMIN和RMAX这两个参数。
下载源代码
 可执行程序(Windows)和示例
运行方法(Win) 
 1. 把编译好的可执行程序“level77.exe”放在bin目录下,若是其它文件名或其它目录,请修改run77.bat。
2. 编写
LEVEL的输入文件。注意:本例提供的输入文件case1.5是LEVEL 7.7格式的,不能用于低版本LEVEL
3. 用鼠标双击run77.bat,进行计算。
4. 拟合的光谱数据位于“RESULTS.OUT”文件中。本例的输出为:
  =================FITTING==================
  === SUMMARY OF SPECTROSCOPIC CONSTANTS ===
  ==========================================

THE BOND LENGTH     Re   =         2.1429 ANGS.
THE ENERGY          Te   =      9907.2439 CM^-1
THE DISSOCIATION E. De   =                N.A.
THE FREQUENCY       We   =       394.7918 CM^-1
                    WeXe =         1.0962 CM^-1
MAX NO. OF LEVELS   Vmax =        10
ROT-FREQUENCY       Be   =         0.1684 CM^-1


*** NOTE ***
 ONLY THE FIRST 11 VIBRATIONAL LEVELS (Vmax=10)
 CAN BE FITTED AT MOST.
升级 
 2016.04.30解决了在Linux系统中无法运行的问题。
2016.12.03:拟合ωe,ωeXe,ωeYe,ωeZe。改善了非谐性较强的势能曲线的拟合。

 

 

Molden2AIM
主页http://zorkzou.github.io/  https://zorkzou.github.io/
功能把Molden文件转化为WFN格式,用于波函分析。从3.0版开始,可以输出.47文件格式,用于NBO分析。从4.0版开始,可以输出.WFX文件格式。

 

 

NBO2Molden
主页http://zorkzou.github.io/  https://zorkzou.github.io/
功能把NBO绘图文件转化为Molden格式。

 

 

Molbo
主页http://zorkzou.github.io/  https://zorkzou.github.io/
功能把Molpro输出文件转化为NBO程序.47输入文件格式,进行NBO分析。同时还计算Mayer键级。

 

 

UniMoVib
主页http://zorkzou.github.io/  https://zorkzou.github.io/
功能利用其他程序计算的振动频率或Hessian等数据,进行最高点群的对称性分析,振动频率计算,和热化学计算,等。具有30多个程序的接口。