<div dir="ltr">PS Sorry for my mistype: I am NOT aware of such a syntax.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">št 8. 9. 2022 o 22:27 Martin Konôpka <<a href="mailto:konopka2010@gmail.com">konopka2010@gmail.com</a>> napísal(a):<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Leili,<br><br>To your 1st question: I am aware that CP2K knows some loops syntax in input files. I am used to VIM editor which makes me easy to edit text files even if they are very large (and I try to keep order in structure of my files and use the diff tool to check them). Anyway, to keep your main input file short, you can use the "@include" statement (which works similarly but not identically as the #include in C codes). Simply put the long part with the LDOS sections to an external file, say "ldos.inc" and then, inside the CP2K input file, type '@include ldos.inc' (if I am not mistaken).<br><br>To your 2nd question - I was not able to find a reference to the subject but I believe the CP2K's way to calculate PDOS and LDOS is a pretty standard one, based on calculating overlaps of eigenstates with some properly defined localised orbitals.<br><br>Best<br>Martin.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">št 8. 9. 2022 o 19:08 Leili Rassouli <<a href="mailto:rassouli.leili@gmail.com" target="_blank">rassouli.leili@gmail.com</a>> napísal(a):<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Martin,<br>Thank you for sharing this information with me. I really appreciate it.<br><br>Regarding the Mulliken analysis of NTO, I use Multiwfn program. I enter .molden file from cp2k as an input file in Multiwfn program for doing this calculation. However, I still have some problems to do this.<br> <br>I did PDOS/LDOS calculation and I have some questions:<br>1- I am working on Fe2O3 supercell with 120 atoms. Do I need to write the LDOS part for each atom (writing 120 times or writing a for loop in Linux to do that) as I want to check the relative contribution of each atom to the total DOS population? <br>2- Do you know how PDOS is calculated in cp2k? Is the relative contribution based on Mulliken analysis? <br>Best,<br>Leili<br><div><br></div><div> &LDOS</div><div> LIST 1 ! the atom number 1<br></div><div> COMPONENTS<br></div><div> &END LDOS </div><div> &LDOS<div> LIST 2 ! the atom number 2 <br></div><div> COMPONENTS<br></div><div> &END LDOS</div><div> . </div><div> . </div><div> .</div><div>&LDOS<div> LIST 120 ! the atom number 120 <br></div><div> COMPONENTS<br></div><div> &END LDOS</div></div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, September 7, 2022 at 9:57:38 AM UTC-4 <a href="mailto:konop...@gmail.com" target="_blank">konop...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Leili,</div><div><br></div><div>As for the 2nd part of your question, I am not able to give an advice. As for the 1st part, I am also not sure but you might perhaps try to use PDOS section, something like the one shown below. <br></div><div><br></div><div>If you want to do a geometry relaxation or molecular dynamics as well, you should somehow prevent calculating and printing all that information at each ionic step. (It would be waste of resources.) I personally do PDOS/LDOS calculations in separate runs after geometry relaxation was completed.<br></div><div><br></div><div>&DFT</div><div>..............</div><div>..............<br></div><div>..............<br></div><div> &PRINT</div><div><br></div><div> &MO_CUBES<br> NLUMO 10<br> STRIDE 2 2 2<br> WRITE_CUBE .FALSE.<br> &END MO_CUBES</div><div><br></div><div> &PDOS</div><div> FILENAME PDOS</div><div> NLUMO 10</div><div> COMPONENTS ! optional<br></div><div> &END PDOS</div><div><br></div><div> &LDOS</div><div> LIST 1 ! the atom number 1<br></div><div> COMPONENTS <br></div><div> &END LDOS</div><div><br></div><div>.....................</div><div>.....................</div><div><br></div><div><div> &LDOS</div><div> LIST 8 ! the atom number 8 say<br></div><div> COMPONENTS <br></div><div> &END LDOS</div><div><br></div><div> &EACH<br> QS_SCF 0<br> &END EACH<br></div><div><br></div></div><div> &END PRINT<br></div><div><br></div><div>&END DFT</div><div><br></div><div>Best regards</div><div>Martin.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">st 7. 9. 2022 o 15:33 Leili Rassouli <<a rel="nofollow">rassoul...@gmail.com</a>> napísal(a):<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all, <div>I want to print the Mulliken analysis for all occupied and unoccupied orbitals. I used PRINT_ALL and PRINT_GOP under <a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html" rel="nofollow" target="_blank">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html" rel="nofollow" target="_blank">FORCE_EVAL</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT.html" rel="nofollow" target="_blank">DFT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/PRINT.html" rel="nofollow" target="_blank">PRINT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/MULLIKEN.html" rel="nofollow" target="_blank">MULLIKEN</a> directory, but it doesn't print the Mulliken analysis for individual orbitals. </div><div>I want to check which atomic orbital has the highest population of alpha and beta in each orbital (for the range of HOMO-10 - LUMO+10). Is there any way to print these values?</div><div><br></div><div>In addition, I want to do the Mulliken analysis for Natural Transition Orbitals (NTOs). Can I use cp2k to do that? </div><div><br></div><div> Best regards,</div><div>Lili</div>
<p></p></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
-- <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">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/4144231f-6138-4f54-a893-cbadca574c1en%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/4144231f-6138-4f54-a893-cbadca574c1en%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div>
<p></p>
-- <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 href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/a65acfd0-c670-40ea-918c-5513ebe9531an%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/a65acfd0-c670-40ea-918c-5513ebe9531an%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div>
<p></p>
-- <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 href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAE0%3D32aoVS0yPazjk5g%2BfLSqfdnJpP_tQdWhAFqq5JqLtmwU1A%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAE0%3D32aoVS0yPazjk5g%2BfLSqfdnJpP_tQdWhAFqq5JqLtmwU1A%40mail.gmail.com</a>.<br />