<p dir="ltr">Dear Matteo,<br>
Thank you so much for your explanatory response.</p>
<p dir="ltr">The very essence of my question is as follows:<br>
If I want to add a term or a constant to the Hamiltonian of a given system, which subroutine or code (*.f90) must be modified?</p>
<div class="gmail_extra"><br><div class="gmail_quote">On Nov 14, 2016 12:02 PM, "Matteo Cococcioni" <<a href="mailto:matteo@umn.edu">matteo@umn.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Ashkan,<div><br></div><div>I'm not sure I understand your question completely, but will attempt an answer.</div><div><br></div><div>First thing to keep in mind is that QE (as any other DFT code) works with the Kohn-Sham single particle Hamiltonian (is that what you are looking for?), whose eigenstates are the ones that allow you to construct the ground state charge density of a system and to evaluate its total energy.</div><div>Due to the Bloch theorem (i.e. to the translational invariance of crystalline solids) this Hamiltonian is block-diagonal on the k-points (crystal momentum of electrons). Also, QE uses an expansion of Kohn-Sham wavefunctions on a plane-wave basis set. So, each of these blocks of the Hamiltonian is npw x npw where npw is the number of plane-waves used in the expansion of the wavefunctions of that specific k-point. QE performs an iterative diagonalization of each block that is involves computing the action of the Hamiltonian onto the current wavefunction H |psi>. This is done in a routine called h_psi which is the driver to compute various parts of the Hamiltonian and to operate it on the wavefunction. So I guess you can look inside that routine in order to find what you need.</div><div><br></div><div>Best regards,</div><div><br></div><div>Matteo</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Nov 13, 2016 at 11:24 PM, ashkan shekaari <span dir="ltr"><<a href="mailto:shekaari@gmail.com" target="_blank">shekaari@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear experts,<br><br></div>How, could one derive the Hamiltonian of a given system resulting in energy eigenvalues (band structure) <b>through the modules or subroutines</b> included in Q. E. integrated suite?<br><br></div><div>Is there anyway for this end?<br></div><div><br></div>For example, if I provide an input scf file, how could I figure out what Hamiltonian will be used in my calculations?<br><br></div>How could I make a change in the Hamiltonian of a system by manipulating subroutines (not the input file)?<br><div><div><br><div><br clear="all"><div><div><div class="m_-9048114753112235206m_-5431018632961166236gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><b>--</b></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><i style="font-size:12.8px">Regards,</i></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><i>Ashkan Shekaari</i></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><i>Plasma Physics Research Center</i></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><i style="font-size:12.8px"><span style="font-size:12.8px">Science and Research Branch</span></i></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><i>I A U, <span style="font-size:12.8px">14778-93855 Tehran, Iran.</span></i></span></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div></div></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a><br></blockquote></div><br></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div></div>