Printing the eignenvectors and values after a run has finished

CaesarSalad

Hello,

Is there a way to print the eigenvectors and eigenvalues of a system after a simulation has finished running?

Chris

eascrizzi

Hi Chris!

You can perform a NEWK properties calculation (see page 346 of CRYSTAL User's Manual) starting from your converged wavefunction and select the printing options for eigenvalues and eigenvectors.

Hope this answers your question

Eleonora

CaesarSalad

Hi Eleonora,

Thanks for your reply. So the calculation can simply be ran as:

NEWK
6 6
1 66 67

? Or do I have to run two separate calculations, one with "1 66" and one with "1 67" on the second line?

Chris

GiacomoAmbrogio

Hi Chris,

If you want to switch on printing options 66 and 67, your input should be:

NEWK
6 6
1 2
66 xxx
67 yyy

where xxx and yyy are the numbers "N" specified in the manual (see page 443, column input)

Let me know if you manage to get the correct printings

CaesarSalad

Hi Giacomo,

Thanks for the input. It seems that this would only print the eigenvalues. The eigenvectors won't be printed and I get this error at the end of the output file:

FERMI ENERGY AND DENSITY MATRIX CALCULATION ON COMPUTED EIGENVECTORS
DENSITY MATRIX AT SCF CYCLE ( 3+1)
SPIN LOCKING: NO ENERGY GAP COMPUTED
CORE DENSITY MATRIX CALCULATION
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTT NEWK TELAPSE 16.38 TCPU 14.78
ERROR **** properties **** END OF DATA IN INPUT DECK

My input is:

NEWK
6 6
1 2
66 -505
67 -505

..ran via crystal_properties binaries. Could it be because 67 works only by printing at the end of the SCF cycles?

Chris

GiacomoAmbrogio

Hi Chris,

The error comes from the missing END line that should close the NEWK block.

As for print option 67: it should work in properties as well as in crystal. That said, these print flags are not very well tested, and some of them may have broken in newer versions. I took a look in the code, and I can confirm that the -505 trick will not work here. Instead, you should set the value to 999 to get the maximum possible output (ie all the virtual states, or 499 for just the occupied states).

One more note: in the NEWK section this option was never extended to parallel execution, so you’ll need to run it on a single CPU core if you want it to print.

Hope this helps.