Section 1
矩阵的行列数:
行数=列数,其中行数等于:各原子s轨道个数+p轨道个数*3之和。
如 CH4 C H 4 ,C有1s,2s,2p轨道,H有1s轨道,则矩阵的行列数为1*(1+1+1*3)+4*1=9
任务目标1:
根据自己要做的分子,先计算矩阵的行列数n。将1到 n2 n 2 顺次填入这个矩阵,并显示到控制台上。
Section 2
输入文件格式::
nBasis
ng am
alpha[1],c[1]
.
.
.
alpha[ng],c[ng]
ng,am
alpha[1],c[1]
.
.
.
alpha[ng],c[ng]
...共有nBasis组
每组是一个shell,依次递增
ng是the number of Gaussian Function(高斯函数的个数)am是角动量
每个shell类型应包含:
type shell_type !shell has the same alpha,coeff,iatom,am
integer :: am !quantum number l
integer :: ng !number of gauss function
integer :: iatom !number of atom
integer :: nicart !( ( am_read + 1 ) * ( am_read + 2 ) ) / 2, number of orbital
Integer basis1,basis2
real*8,allocatable :: alpha(:)
real*8,allocatable :: coeff(:)
end type