- Posts: 3
- Thank you received: 1
errors in the MOLDEN file generated by MRCC
- wzou
- Topic Author
- Offline
- New Member
Less
More
10 years 2 days ago #65
by wzou
errors in the MOLDEN file generated by MRCC was created by wzou
Dear MRCC developers,
There are two errors in the MOLDEN file generated by MRCC. Here is the input I used:
1, If the occupation number is specified (occ=... in the above input), the occupation numbers of the occupied orbitals in the MOLDEN file are all zero.
2, the ordering of orbital coefficients of spherical d, f, and g functions is wrong (the ordering of Cartesian functions are correct).
Here are some coefficients of MO #11-15 taken from the above calculation where the d functions have large contributions, and compared with the results of Gaussian.
Best regards,
Wenli
There are two errors in the MOLDEN file generated by MRCC. Here is the input I used:
Code:
# test
basis=cc-pVDZ
calc=scf
occ=3,0,0,0
mem=1GB
geom
B
H 1 r1
r1=1.2
1, If the occupation number is specified (occ=... in the above input), the occupation numbers of the occupied orbitals in the MOLDEN file are all zero.
2, the ordering of orbital coefficients of spherical d, f, and g functions is wrong (the ordering of Cartesian functions are correct).
Here are some coefficients of MO #11-15 taken from the above calculation where the d functions have large contributions, and compared with the results of Gaussian.
Code:
11 12 13 14 15
(sg) (pi) (pi) (dlta) (dlta)
energy 0.80906 0.87116 0.87116 0.90397 0.90397
(MOLDEN by MRCC)
10 d 0 0.00000 0.00000 0.00000 0.00000 1.00000
11 d+1 0.00000 0.00000 0.93042 0.00000 0.00000
12 d-1 -0.00000 0.00000 0.00000 1.00000 0.00000
13 d+2 0.00000 0.93042 0.00000 0.00000 0.00000
14 d-2 -0.61900 0.00000 0.00000 0.00000 0.00000
(by Gaussian)
10 d 0 0.61900 0.00000 0.00000 0.00000 0.00000
11 d+1 0.00000 0.93042 0.00000 0.00000 0.00000
12 d-1 0.00000 0.00000 0.93042 0.00000 0.00000
13 d+2 0.00000 0.00000 0.00000 1.00000 0.00000
14 d-2 0.00000 0.00000 0.00000 0.00000 1.00000
Best regards,
Wenli
Please Log in or Create an account to join the conversation.
- wzou
- Topic Author
- Offline
- New Member
Less
More
- Posts: 3
- Thank you received: 1
10 years 1 day ago #70
by wzou
Replied by wzou on topic errors in the MOLDEN file generated by MRCC
Dear developers,
The second error about spherical functions can be easily fixed. In the source code file integ.f, go to the line 5112, replace
by
and recompile the program.
Best regards,
Wenli
The second error about spherical functions can be easily fixed. In the source code file integ.f, go to the line 5112, replace
Code:
if(iang.ge.1) then
Code:
if(iang.eq.1) then
Best regards,
Wenli
The following user(s) said Thank You: kallay
Please Log in or Create an account to join the conversation.
- kallay
- Offline
- Administrator
- Mihaly Kallay
10 years 1 day ago #71
by kallay
Best regards,
Mihaly Kallay
Replied by kallay on topic errors in the MOLDEN file generated by MRCC
Dear Wenli,
Thank you very much for reporting this. We apologize for the errors, which have been fixed. A new program release will be available soon. For the meantime, you can simply fix the bugs as follows.
1) In file scf.f add the following lines to the end of subroutine occparser:
if(scftype.eq.1) then
do ii=1,nir
orbperir(ii+nir)=orbperir(ii)
enddo
endif
2) In file integ.f, in subroutine moldeninp replace line
if(iang.ge.1) then
with
if(iang.eq.1) then
Thank you very much for reporting this. We apologize for the errors, which have been fixed. A new program release will be available soon. For the meantime, you can simply fix the bugs as follows.
1) In file scf.f add the following lines to the end of subroutine occparser:
if(scftype.eq.1) then
do ii=1,nir
orbperir(ii+nir)=orbperir(ii)
enddo
endif
2) In file integ.f, in subroutine moldeninp replace line
if(iang.ge.1) then
with
if(iang.eq.1) then
Best regards,
Mihaly Kallay
Please Log in or Create an account to join the conversation.
Time to create page: 0.037 seconds