Re: Density matrices for NC DFT ( No.1 ) 
 Date: 2017/06/17 22:26
 Name: Daniil
 I browsed the code and found this comment in several files:
CDM[0] Re alpha alpha density matrix CDM[1] Re beta beta density matrix CDM[2] Re alpha beta density matrix CDM[3] Im alpha beta density matrix iDM[0][0] Im alpha alpha density matrix iDM[0][1] Im beta beta density matrix So, there are 6 matrix components in the calculation process (There must be a beta alpha matrix also, but it could be dependent on the alpha beta one). However, I was unable to trace how do these matrices convert into a DM matrix which is written as 4 components in the .scfout file.

Re: Density matrices for NC DFT ( No.2 ) 
 Date: 2017/06/24 17:18
 Name: T. Ozaki
 Hi,
Each density matrix consists of
DM_{alpha,alpha} = CDM[0] + I*iDM[0] DM_{alpha,beta} = CDM[2] + I*CDM[3] DM_{beta,alpha} = CDM[2]  I*CDM[3] DM_{beta,beta} = CDM[1] + I*iDM[1]
Unfortunately, the scfout file contains only CDM[03]. Sorry for your inconvenience.
Regards,
TO

Re: Density matrices for NC DFT ( No.3 ) 
 Date: 2017/06/25 04:18
 Name: Daniil
 Dear Prof. Ozaki,
Thank you for your answer. So, DM array at line 743 of SCF2File.c corresponds to CDM?
However, I need to obtain the full matrix. Is it possible to modify code to add iDM to scfout? For example, if global iDM array really corresponds to iDM in your expressions and is available at scfout stage, I could add it to output in SCF2File.c and read it in read_scfout.c in the same way, as the DM array. Is it true, that iDM array is iDM in the above expressions?
Best regards, Daniil.

Re: Density matrices for NC DFT ( No.4 ) 
 Date: 2017/06/25 17:30
 Name: T. Ozaki
 Hi,
iDM in my previous post is the global array iDM. Of course, you can modify SCF2File.c depending on your purpose.
Regards,
TO

Re: Density matrices for NC DFT ( No.5 ) 
 Date: 2017/06/26 03:06
 Name: Daniil
 Thanks, I'll try.
Is it true that iDM[0] is always of size 2 and fully allocated, even in case of SpinPolarization=Off or On, when DM[0] is of size 1 or 2 respectively?
Best regards, Daniil.

Re: Density matrices for NC DFT ( No.6 ) 
 Date: 2017/06/27 14:21
 Name: T. Ozaki
 Hi,
> Is it true that iDM[0] is always of size 2 and fully allocated, > even in case of SpinPolarization=Off or On, when DM[0] is of size 1 or 2 respectively?
Yes, it is. However, iDM has the meaning only if SpinPolarization=NC.
Regards,
TO

Re: Density matrices for NC DFT ( No.7 ) 
 Date: 2017/06/27 22:40
 Name: Daniil
 Thanks, I managed to do this.
Would you consider adding iDM to the scfout output in the future releases?
Best regards, Daniil.
