Hi Michael,<br /><div><br /></div><div>I aim to develop the code further. But first, I check and clean the old one. The changes are mainly internal and do not change the results. So far, I have made some improvements in the calculation of separate electrodes. I don't think your system is visibly affected. The state-of-the-art I will describe in the paper and on-site tutorials, which I am preparing now.</div><div><br /></div><div>Your approach to checking the transmission function is a good one. Yes, if one needs a voltage-dependent transmission, it is computationally expensive. )</div><div>One of the possible solutions is to use DFTB or xTB instead of DFT. It is one of the things I want to add to the code.</div><div><br /></div><div>Best wishes,</div><div>Dmitry</div><div><br /></div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Michael LaCount schrieb am Samstag, 15. November 2025 um 06:30:31 UTC+1:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hello Dmitry,<div><br></div><div>My post was when I was initially learning the NEGF calculations. I have since solved the problems I was having, and reproduce the transmission function from the technical report. As well as producing what appear to be reasonable I(V) curves.</div><div><br></div><div>The I(V) curves I made are based on Eq 21 from the technical report using the zero-bias (ELECTRIC_POTENTIAL [eV] 0). Your comment "If the energy levels of the central system are not changed at finite voltage" leaves me a bit worried that I was not careful enough. Is there a simple way to test whether or not the assumption is valid? My guess would be something like apply a bias voltage up to the max of the I(V) curve and see if the transmission function changes significantly. What I am trying to simulate is a small semiconductor device. For simplicity sake lets say it has a band gap of 1 eV, and make I(V) curves in the range of +/- 3V. If the assumption brakes down how would I go about finding the I(V) curve, I could imagine finding the transmission function at various applied biases, but that seems like a computationally expensive approach.</div><div><br></div><div>I have been using a recent but not newest version of CP2K 2024.2. I am curious if the changes to the NEGF code are documented anywhere, and what impact they would have either in terms of performance or results.</div><div><br></div><div>Best,</div><div>Michael LaCount</div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, November 14, 2025 at 12:41:23 AM UTC-8 Dmitry Ryndyk wrote:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Michael, <div><br></div><div>to see the same result, as in the report, you should take zero voltage and other energy limits:<br><div><br></div><div>&PRINT<br> &DOS<br> FILENAME device<br> FROM_ENERGY 0.272240982<br> N_GRIDPOINTS 401<br> TILL_ENERGY 0.492829218</div></div><div><div><br> &END DOS<br> &TRANSMISSION<br> FILENAME transm<br></div></div><div><div> FROM_ENERGY 0.272240982<br> N_GRIDPOINTS 401<br> TILL_ENERGY 0.492829218<br> &END TRANSMISSION<br> &END PRINT</div><div><br></div><div>The Fermi level is 0.38253510. The result is in attachment. It is not exactly the same, but it can be for many reasons. </div><div><br></div><div>I am currently revising the NEGF code, and may change the energy levels to be around the Fermi level. If you want to use NEGF, I recommend the latest development version of CP2K. And I will try to answer further questions. </div><div><br></div><div>Concerning your second question. It is the advantage of the NEGF approach that T(E) is changed with bias voltage. But as a simplest way to get I(V), the zero-voltage transmission can sometimes be used. Very carefully to make sense. It depends on the problem, of course. If the energy levels of the central system are not changed at finite voltage, one can omit self-consistency at finite voltage and use T_0(E).</div><div><br></div><div>Best wishes,</div><div>Dmitry Ryndyk<br><br></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Michael LaCount schrieb am Donnerstag, 10. April 2025 um 05:20:41 UTC+2:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I am trying to reproduce the transmission coefficient plot found in Figure 4 of the CP2K Electron Transport based on Non-Equilibrium-Greens-Functions Method: eCSE 08-09 Technical Report (see attached). From what I can tell the system is the same as the one found in 'QS/regtest-negf-fft/au111_c6h4s2_gamma_0.50V.inp'. However when I run that job unmodified except to add print commands in the NEGF section (see below) I don't get anything close to what was found in the Technical Report. The job input and output files are attached.<div><br></div><div>I also have a further question regarding this test. I want to create I-V curves, from what I understand I should therefore use the zero bias transmission coefficients, and set NEGF/CONTACT/ELECTRIC_POTENTIAL to 0 correct?</div><div><br></div><div>Any clarity on where I am going wrong would be appreciated. <br><div><br></div><div> &PRINT<br> &DOS<br> FILENAME device<br> FROM_ENERGY -0.2<br> N_GRIDPOINTS 401<br> TILL_ENERGY 0.2<br> &END DOS<br> &TRANSMISSION<br> FILENAME transm<br> FROM_ENERGY -0.2<br> N_GRIDPOINTS 401<br> TILL_ENERGY 0.2<br> &END TRANSMISSION<br> &END PRINT<br></div></div></blockquote></div></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 visit <a href="https://groups.google.com/d/msgid/cp2k/619d6289-6c52-4310-9082-c0b220e10d74n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/619d6289-6c52-4310-9082-c0b220e10d74n%40googlegroups.com</a>.<br />