<div>Hi Elaheh,</div><div><br></div><div>
The trunk version of cp2k now has a switch REAL_SPACE in DFT%PRINT%S_CSR_WRITE and
DFT%PRINT%KS_CSR_WRITE to print the real space representation of the overlap and KS matrix.This key word also activates the printing of the cell indices.</div><div><br></div><div>Cheers,</div><div>Fabian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, 5 February 2021 at 12:59:59 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;">and it goes without saying that first block is the Hamiltonian for <0|OP|0> with onsite values.so<div>1->0 0 0 </div><div>2-> 1 0 0 </div><div>3->2 0 0 </div><div>Kind regards,</div><div>Elaheh Akhoundi<br><div><br></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 5, 2021 at 12:57:33 PM UTC+1 Elaheh Akhoundi 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">Dear Prof.
Hutter,<div>I picked a small one-dimensional structure, and used a short range basis which is sufficient for the band structure (I checked). in that case, there will be only two neighboring computational cells with significant contributions. so only two possibilities for the index to cell matching map. </div><div>1-> 1 0 0</div><div>2-> 2 0 0</div><div>or</div><div>1-> 2 0 0 </div><div>2-> 1 0 0 </div><div>that's how I created a map just to do a test to make sure I'll get the right structure. then I'll proceed to create the map with the source code. </div><div>Thank you very much, I will double check everything, hopefully I find the problem.</div><div>Thank you again,</div><div>Elaheh Akhoundi</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 5, 2021 at 12:43:22 PM UTC+1 jgh 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">Hi
<br>
<br>I don't understand what you did. You need the index to cell matching
<br>from the code. If you create the list yourself, how do you know
<br>it is the same?
<br>I can definitely confirm that the L <-> -L symmetry is used and that
<br>the Fourier transform of the real space matrices to k space works.
<br>
<br>regards
<br>
<br>Juerg Hutter
<br>--------------------------------------------------------------
<br>Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" rel="nofollow" target="_blank">++41 44 635 4491</a>
<br>Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" rel="nofollow" target="_blank">++41 44 635 6838</a>
<br>Universität Zürich E-mail: <a rel="nofollow">h...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a rel="nofollow">c...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a rel="nofollow">c...@googlegroups.com</a>>
<br>From: "Elaheh Akhoundi"
<br>Sent by: <a rel="nofollow">c...@googlegroups.com</a>
<br>Date: 02/05/2021 09:12AM
<br>Subject: Re: [CP2K:14697] plotting band structure
<br>
<br>Hi Juerg Hutter,Thank you very much for your reply. I did a quick test for a simple 1D structure, I created the list of indices myself, as there were only two neighboring computational cells. I took into account the L --> -L symmetry by considering , the Hamiltonian/overlap matrix of -1 0 0 to be the transpose of those of 1 0 0. It still gave wrong results and now even gamma point eigen values are wrong too. before, by using the map with negative L , I was getting the right values for the gamma point for every material I tested. I am using CP2K V. 7.1.
<br>Many thanks,
<br>Elaheh
<br>
<br>On Thursday, February 4, 2021 at 2:27:28 PM UTC+1 jgh wrote:
<br>Hi
<br>
<br>you did it the wrong way round. You need to loop over the number
<br>of indices and then print index_to_cell.
<br>
<br>1 -> 0 0 0
<br>2 -> 0 0 1
<br>...
<br>
<br>Index is now the position of the matrix in the output (0||i)
<br>and there will be no L, -L pair.
<br>
<br>regards
<br>
<br>Juerg Hutter
<br>--------------------------------------------------------------
<br>Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" rel="nofollow" target="_blank">++41 44 635 4491</a>
<br>Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" rel="nofollow" target="_blank">++41 44 635 6838</a>
<br>Universität Zürich E-mail: <a rel="nofollow">h...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a rel="nofollow">c...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a rel="nofollow">c...@googlegroups.com</a>>
<br>From: "Elaheh Akhoundi"
<br>Sent by: <a rel="nofollow">c...@googlegroups.com</a>
<br>Date: 02/04/2021 02:06PM
<br>Subject: Re: [CP2K:14687] plotting band structure
<br>
<br>Hi Juerg Hutter,
<br>Thank you very much for your help. I have already printed index to cell mapping by changing the source code. I used a neighbor_list iterator and printed cell(1:3) and also cell_to_index values. It looks something like this:
<br>index= -2 0 0
<br> block # 6
<br> -----------
<br> -----------
<br> index= -1 0 0
<br> block # 2
<br> -----------
<br> -----------
<br> index= 0 0 0
<br> block # 1
<br> -----------
<br> -----------
<br> index= 1 0 0
<br> block # 3
<br> -----------
<br> -----------
<br> but I am noticing that -L cells are also printed, you can see in the indices file that 100 and -100 both exist. also when I get the Kohn Sham Hamltonian 85 blocks are printed. which is consistent with what I get from the cell to index map. so it seems like L and -L computational cell Hamiltonians are both being printed? even though in the Hamiltonian file they aren't transpose of each other, which I think may be wrong.
<br>Could you possibly give me some advice?
<br>Many thanks,
<br>Elaheh
<br>
<br>
<br> --
<br> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp...@googlegroups.com</a>.
<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/6833475d-5d54-4066-84bb-4f0f4b17b88en%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/6833475d-5d54-4066-84bb-4f0f4b17b88en%2540googlegroups.com&source=gmail&ust=1619892771094000&usg=AFQjCNEkl_gRSWbjZPe06tfJAJgLeqPPqg">https://groups.google.com/d/msgid/cp2k/6833475d-5d54-4066-84bb-4f0f4b17b88en%40googlegroups.com</a>.
<br>
<br>
<br> --
<br> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp...@googlegroups.com</a>.
<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e56735c2-fc66-4896-97d8-38c28051fb49n%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/e56735c2-fc66-4896-97d8-38c28051fb49n%2540googlegroups.com&source=gmail&ust=1619892771094000&usg=AFQjCNGjI2tQT5LrToRyjXHb8HO17B0dCA">https://groups.google.com/d/msgid/cp2k/e56735c2-fc66-4896-97d8-38c28051fb49n%40googlegroups.com</a>.
<br>
<br></blockquote></div></blockquote></div></blockquote></div>