From NWChem
Viewed 2548 times, With a total of 6 Posts

Clicked A Few Times
Threads 12
Posts 18


6:47:53 AM PDT  Mon, Jul 15th 2013 

Hi,
By reading the previous posts, I know that the Molecular orbitals after hf or dft calculation are orthonormal, but eigenvectors probably not. The following relation can be verified:
C.transpose * S * C = I
(C is the eigenvector, and S is the overlap matrix between atomic orbitals.)
I code up some scripts to extract the matrix from NWChem, and try to verify this expression.
1. I turn on the command "print "overlap"" in input NWChem script, then write up a script to extract the overlap matrix from the output.
2. I use the utility mov2asc to extract the eigenvectors from the .movecs file.
3. Use the C and S obtained in previous two steps to verify the above relation.
Following the above steps, I failed to get the identity matrix by evaluating C.transpose * S * C.
What is the possible problem of my procedure?
Does anyone can give some suggestions about this?




Edoapra Forum:Admin, Forum:Mod, bureaucrat, sysop


Forum Vet
Threads 9
Posts 1541


4:43:54 PM PDT  Thu, Jul 18th 2013 

Quote:Alexx Jul 15th 6:47 amHi,
1. I turn on the command "print "overlap"" in input NWChem script, then write up a script to extract the overlap matrix from the output.
The correct print instruction for dumping the overlap matrix is
print "ao overlap"
"print overlap" will print the overlap integrals, instead.
Cheers, Edo

Edited On 4:44:41 PM PDT  Thu, Jul 18th 2013 by Edoapra




Clicked A Few Times
Threads 12
Posts 18


7:11:01 AM PDT  Mon, Aug 12th 2013 

Hi, Edo
I can verify the expression for the single atom now. But I still can't get it right for the molecules like H2O.
For water molecule, the diagonal elements of C.transpose*S*C are all ones, but for some offdiagonal elements, the value is not zero, some of which are even around 0.0001.
For heavy molecules like uranium compounds, some offdiagonal elements are even bigger, around 0.1
Can you give some hints about this?
Quote:Edoapra Jul 18th 5:43 pmQuote:Alexx Jul 15th 6:47 amHi,
1. I turn on the command "print "overlap"" in input NWChem script, then write up a script to extract the overlap matrix from the output.
The correct print instruction for dumping the overlap matrix is
print "ao overlap"
"print overlap" will print the overlap integrals, instead.
Cheers, Edo




Edoapra Forum:Admin, Forum:Mod, bureaucrat, sysop


Forum Vet
Threads 9
Posts 1541


9:58:26 AM PDT  Mon, Aug 12th 2013 

Increase number of digits in Overlap Matrix printout

Alexx
Have you tried to increase the number of digits in the Overlap Matrix printout
as suggested in
http://www.nwchemsw.org/index.php/Special:AWCforum/st/id910/more_significant_digits.html




Clicked A Few Times
Threads 12
Posts 18


10:29:33 AM PDT  Mon, Aug 12th 2013 

Hi,
I find that "print overlap" prints overlap matrix elements, which does the same thing as "print ao overlap" does, and "print overlap" prints overlap elements with more precision, which is 8 digits.
Do I need even more digits to do that?
thanks,




Edoapra Forum:Admin, Forum:Mod, bureaucrat, sysop


Forum Vet
Threads 9
Posts 1541


12:31:15 PM PDT  Mon, Aug 12th 2013 

How many digits do you have for the eigenvectors?
Cheers, Edo




Clicked A Few Times
Threads 12
Posts 18


12:40:32 PM PDT  Mon, Aug 12th 2013 

Hi, Edo !
The number of digits for the eigenvectors is 15.
Below is a small section of the output:
0.182031731624931E04 0.103280786104367E04 0.315402471342734E04
0.381586117380476E04 0.789728509277435E04 0.872883749498206E04
0.462811845480190E04 0.223511333780241E03 0.696050502583492E03
0.110518876210707E03 0.209263633645397E02 0.141537941674604E03
0.636958647241003E02 0.655655688213226E02 0.109112070692797E02
thanks,
Quote:Edoapra Aug 12th 1:31 pmHow many digits do you have for the eigenvectors?
Cheers, Edo



AWC's:
2.5.10 MediaWiki  Stand Alone Forum Extension
Forum theme style by: AWC