<div>Dear 
Elaheh,</div><div><br></div><div>I don't know how to print the indices easily. But you can use DFT%PRINT%KS_CSR_WRITE and 
DFT%PRINT%S_CSR_WRITE to write the Hamiltonain and overlap matrices to the disk. If you use k-points a file for each k-point is created in the ordering given in the output file.</div><div><br></div><div>Cheers,</div><div>Fabian<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 2 December 2020 at 12:07:10 UTC+1 akho...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hello,</div><div>I have extracted the Kohn Sham Hamiltonian using   &AO_MATRICES HIGH ----> KOHN_SHAM_MATRIX .TRUE. , in order to plot the band structure myself using the Hamiltonian matrix. However, the format of the output is not clear to me. I believe each block of Hamiltonian corresponds to the interaction of the primary cell with a specific neigboring cell. but it is not specified in the output file which block is related to which neighboring cell, i.e. what are the cell indices corresponding to each block. Is there any way that I can access the cell indices related to each Hamiltonian block?</div><div>I would greatly appreciate your help.</div><div>kind regards,</div><div>Elaheh Akhoundi</div></blockquote></div>