<div dir="ltr"><div>Dear Patrick,</div><div><br></div><div> Thanks a lot.</div><div>All the best,</div><div>Elham<br></div></div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 24, 2020 at 12:47 AM Patrick Gono <<a href="mailto:patri...@gmail.com" target="_blank">patri...@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 Elham,</div><div><br></div><div>You need two calculations. A bulk calculation where you align the band edges with respect to the average electrostatic potential. And a slab calculation with vacuum, where you align the average electrostatic potential in the center layers of the slab with respect to the vacuum potential.</div><div><br></div><div>In CP2K, the average electrostatic potential over the whole cell is set to zero. Hence, the bulk calculation only consists in running a CELL_OPT of a big enough supercell, and taking the values of the Fermi level and the HOMO-LUMO gap. The Fermi level represents the position of the valence band with respect to the average electrostatic potential, and the Fermi level + gap is then the conduction band. This calculation gives you the alignment of the band edges with respect to the average electrostatic potential.<br></div><div><br></div><div>For the vacuum calculation, you have to make sure you have a thick enough slab (so that the electrostatic potential in the middle layers behaves like the bulk one) and a thick enough vacuum layer, so that there the potential is flat. Run a GEO_OPT of the slab with vacuum, and then perform an ENERGY calculation on the resulting geometry with the V_HARTREE_CUBE option in the PRINT section inside DFT set to .TRUE.</div><div>Next, project the .cube file of the average potential along the axis perpendicular to the interface. Then, take the value in the middle of the vacuum region, and subtract the average over the central-most layer(s) inside the slab. This gives you the alignment of the average electrostatic potential with respect to vacuum.</div><div><br></div><div>You see now that combining the results of the two yields the desired alignment.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 23 Jan 2020 at 14:27, Elham M.Goliaei <<a href="mailto:e.go...@gmail.com" target="_blank">e.go...@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>Hi Patrick,</div><div> Thanks, it was informative, first, I wanna calculate band alignment versus vacuum, for example according to Fig. 1;</div><div> <a href="https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.102.036402" target="_blank">https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.102.036402</a></div><div><br></div><div>How can I do these calculation to get band alignment versus vacuum?</div><div><br></div><div>All the best,</div><div>Elham<br></div><br><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jan 22, 2020 at 6:33 AM Patrick Gono <<a href="mailto:patri...@gmail.com" target="_blank">patri...@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 Elham,</div><div><br></div><div>In the work of Guo et al. referenced above, it is stated that a molecular dynamics simulation at the level of PBE was performed at the desired temperature. The projected density of states was averaged over the duration of the simulation. Since the average PDOS has long tails at the band edges, a linear extrapolation of the band edge PDOS was carried out. It is these values that are used as the temperature corrected band edges. See Fig. 3 of this paper by Wiktor et al.: <a href="https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.1.022401" target="_blank">https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.1.022401</a></div><div><br></div><div>It is possible that you require the use of a more advanced electronic structure method, and hence ab initio molecular dynamics may not be feasible. In that case, you can follow the approach outlined in the above reference -- run MD at the semilocal level, select a set of configurations from the trajectory, perform single point (ENERGY) calculations with your desired method, compute the PDOS, average over all snapshots, and carry out the linear extrapolation of the PDOS tails around the band gap mentioned previously.<br></div><br><div>Alternatively, a one shot calculation of the phonon spectrum can be used to estimate the band gap renormalization due to finite temperature. See the work by Prof. Giustino, for example:</div><div><a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.94.075125" target="_blank">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.94.075125</a></div><div><br></div><div>Hope this helps.</div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 22 Jan 2020 at 14:39, Elham M.Goliaei <<a href="mailto:e.go...@gmail.com" target="_blank">e.go...@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>Thanks, can you guide me how I can calculate T-dependent corrections of the band edges (Eq. 4 and 5).</div><div><br></div><div>All the best,</div><div>Elham<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 21, 2020 at 4:27 AM Patrick Gono <<a href="mailto:patri...@gmail.com" target="_blank">patri...@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 Elham,</div><div><br></div><div>I recommend the work by Guo et al. where the band alignment at the semiconductor water interface is studied in detail:</div><div><a href="https://pubs.acs.org/doi/abs/10.1021/acs.chemmater.7b02619" target="_blank">https://pubs.acs.org/doi/abs/10.1021/acs.chemmater.7b02619</a></div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 20 Jan 2020 at 11:07, Elham M.Goliaei <<a href="mailto:e.go...@gmail.com" target="_blank">e.go...@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>Hi Matthias,</div><div><br></div><div> I wanna calculate band alignment of a semiconductor and band alignment at the semiconductor-water interface, can you guide me how to calculate.</div><div><br></div><div>All the best,</div><div>Elham<br></div><div><span style="font-size:10.2025px;font-family:serif"></span><span style="font-size:10.2025px;font-family:serif"></span><span style="font-size:10.2025px;font-family:serif"></span><span style="font-size:10.2025px;font-family:serif"></span><span style="font-size:10.2025px;font-family:serif"></span><span style="font-size:10.2025px;font-family:serif"></span><br><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)" lang="EN-GB"></span></div></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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/408361f6-4639-4338-b63f-ace47ceeabfe%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/408361f6-4639-4338-b63f-ace47ceeabfe%40googlegroups.com</a>.<br> </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPL1_v-DETREjQUWmFAA26%2BLXxj%3DEdgf%2BCqibn4NOS44QHqpNw%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAPL1_v-DETREjQUWmFAA26%2BLXxj%3DEdgf%2BCqibn4NOS44QHqpNw%40mail.gmail.com</a>.<br> </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAK_aGnJSJF61f7caSHZ-5gR_sM-O28CymKxxxEwCJgVRy8%2BK0g%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAK_aGnJSJF61f7caSHZ-5gR_sM-O28CymKxxxEwCJgVRy8%2BK0g%40mail.gmail.com</a>.<br> </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPL1_v_1NCmGQ9O-jF6JCxjeserx%3D0DEWquoyGM3SAr8Es9FHA%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAPL1_v_1NCmGQ9O-jF6JCxjeserx%3D0DEWquoyGM3SAr8Es9FHA%40mail.gmail.com</a>.<br> </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAK_aGnKVr_mfZ%3DugyT9NmL%2BCiYONuo_qtsgsv5FgCqj-wvioNQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAK_aGnKVr_mfZ%3DugyT9NmL%2BCiYONuo_qtsgsv5FgCqj-wvioNQ%40mail.gmail.com</a>.<br> </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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPL1_v_8%2BQ69atUWrpUEqQfqq_bQ%3DZj3CyyD0ZhKU1fmHohJsQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAPL1_v_8%2BQ69atUWrpUEqQfqq_bQ%3DZj3CyyD0ZhKU1fmHohJsQ%40mail.gmail.com</a>.<br> </blockquote></div>