<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
@font-face
        {font-family:"Segoe UI";
        panose-1:2 11 5 2 4 2 4 2 2 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
pre
        {mso-style-priority:99;
        mso-style-link:"Preformattato HTML Carattere";
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Courier New";}
span.PreformattatoHTMLCarattere
        {mso-style-name:"Preformattato HTML Carattere";
        mso-style-priority:99;
        mso-style-link:"Preformattato HTML";
        font-family:Consolas;
        mso-fareast-language:IT;}
span.StileMessaggioDiPostaElettronica20
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 2.0cm 2.0cm 2.0cm;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=IT link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Hello Andrea,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>DFTB never calculates the full charge density due to electrons+nuclei, but only deviations from reference atomic charge densities (neutral atoms). Furthermore such deviations are expanded in spherical harmonics and in the traditional dftb implementation only the s-waves are considered with a simple Mulliken analysis. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>In practice what is computed in dftb is not suitable for your calculation because the potential seen for a neutral CNT will be 0! This does not mean there are no electrons or nuclei around. Just all their interactions are included in the SK parameters.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'> <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>If you are good in programming you could compute the charge density from the wavefunctions, extract our mudpack poisson solver and create a stand-alone interface for it, read-in the charge density using an exceptionally fine grid and finally solve the potential… <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>if you do this, let me know that we will include it in the dftb+ tools !<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Maybe it is easier for you to carry calculations using DFT packages, but bare in mind that most use pseudopotentials. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>If you really need the ions and the core electrons than you can use CP2K, GPAW, Gaussian…<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Hope this helps,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Alessandro. <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Segoe UI","sans-serif"'>Da:</span></b><span style='font-size:10.0pt;font-family:"Segoe UI","sans-serif"'> DFTB-Plus-User [mailto:dftb-plus-user-bounces@mailman.zfn.uni-bremen.de] <b>Per conto di </b>Andrea Pedrielli<br><b>Inviato:</b> d</span><span lang=EN-US style='font-size:10.0pt;font-family:"Segoe UI","sans-serif"'>omenica 30 ottobre 2016 18:32<br><b>A:</b> User list for DFTB+ related questions<br><b>Oggetto:</b> Re: [DFTB-Plus-User] Compute electrostatic potential with NEGF<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal>Thanks again,<o:p></o:p></p><div><p class=MsoNormal>essentially the aim of the work is to simulate the electron microscopy holography response that is proportional to the integrated coulomb potential along the beam direction. So, one electron goes through the sample (in this case the nanotube) and feels a phase shift due to the coulomb potential in each point along the path. The potential has to be computed in each point of the space, also inside the atomic orbitals. In the best case I would compute the potential from the charge density + nuclei, like in DFT. I understand that the orbitals in DFTB are not changed in shape. However, supposing a not so large change in the shape of the orbitals it will be useful also only compute the potential starting from the charge density of the non deformed orbitals. I tried to do this plotting the the charge density using waveplot to a cube file and using an Ewald summation over all the points. The problem is that the computational cost is very high because I need a fine grid to describe the charge density inside the atoms. If I understand well the Poisson solver doesn't take as input the charge density that one can plot with waveplot but only the difference of charge with respect to the reference configuration? <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Andrea <br><br>Il domenica 30 ottobre 2016, Gabriele Penazzi <<a href="mailto:penazzi@uni-bremen.de">penazzi@uni-bremen.de</a>> ha scritto:<o:p></o:p></p><div><p>Hi Andrea,<o:p></o:p></p><p>I am not sure what you want to do, but dftb as a method may not include the details you need. <o:p></o:p></p><p>I see now the point about mulliken charges from your previous mail, but the use of mulliken charges is pretty much a direct consequence of the way dftb (at least in its current and most popular version) is formulated. Nuclei+electron density does not make too much sense, as you work in a monopole approximation based on differences with respect to a reference configuration. The s-like projection of the potential reflects this and is consistent with the usual scc-dftb based on gamma-functional.<o:p></o:p></p><p>For example, you have no potential in a system without self-consistent fluctuation. And even if you compare two system with different hybridization, if the scc component is zero the real space potential will be the same, i.e. zero. So my guess is no, you can not really have detailed information about he charge density between atoms and you probably have to go at DFT level. <o:p></o:p></p><p>Best,<o:p></o:p></p><p>Gabriele<o:p></o:p></p><p><o:p> </o:p></p><div><p class=MsoNormal>On 10/29/2016 11:55 PM, Andrea Pedrielli wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>Thank you Gabriele,  <o:p></o:p></p><div><p class=MsoNormal>If I could, I would ask you another thing. I read in the online tutorial on silicon nanowire that the charge density is expanded in s spherical orbitals weighted with Mullen charges. Part of my interest in compute charge density and potential is take into account the charge density between the atoms, due to the bonds. If the charge density is expanded in s-like orbitals the charge density between the atoms (for example in a nanotube) is still well described? <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Andrea<o:p></o:p></p></div><div><p class=MsoNormal><br>Il sabato 29 ottobre 2016, Gabriele Penazzi <<a href="javascript:_e(%7B%7D,'cvml','penazzi@uni-bremen.de');" target="_blank">penazzi@uni-bremen.de</a>> ha scritto:<o:p></o:p></p><div><p>Hi Andrea,<o:p></o:p></p><p>yes, you can do a non transport calculation with the real space Poisson solver, I did the same for the same reason. Whether you can simulate your system or not, however, may depend on the boundary conditions you want to impose. For example, all periodic does not work if I remember right, but for a 1D system you should be able to set up a well defined calculation. <o:p></o:p></p><p>Gabriele<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On 10/28/2016 10:51 PM, Andrea Pedrielli wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal>Hi users,  <o:p></o:p></p><div><p class=MsoNormal>I need to compute the electrostatic potential produced by the charge density from a dftb+ calculation. I know that there is no possibility using dftb+, but I have seen that in the NEGF package there is a Poisson solver. In particular I have to compute the inner and the outer potential of a carbon nanotube, so a 1D object. Can I use a NEGF package for a non-transport calculation with the aim of extract the electrostatic potential? I underline that I need the potential due to the total charge density+nuclei and Mulliken charges are not suitable for my purpose.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Thank you in advance,<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Andrea<o:p></o:p></p></div><p class=MsoNormal style='margin-bottom:12.0pt'><o:p> </o:p></p><pre>_______________________________________________<o:p></o:p></pre><pre>DFTB-Plus-User mailing list<o:p></o:p></pre><pre><a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de">DFTB-Plus-User@mailman.zfn.uni-bremen.de</a><o:p></o:p></pre><pre><a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" target="_blank">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a><o:p></o:p></pre></blockquote><pre>-- <o:p></o:p></pre><pre>Dr. Gabriele Penazzi<o:p></o:p></pre><pre>BCCMS - University of Bremen<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><a href="http://www.bccms.uni-bremen.de/" target="_blank">http://www.bccms.uni-bremen.de/</a><o:p></o:p></pre><pre><a href="http://sites.google.com/site/gabrielepenazzi/" target="_blank">http://sites.google.com/site/gabrielepenazzi/</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>phone: +49 (0) 421 218 9328<o:p></o:p></pre><pre>fax: +49 (0) 421 218 4764<o:p></o:p></pre></div></div><pre>_______________________________________________<o:p></o:p></pre><pre>DFTB-Plus-User mailing list<o:p></o:p></pre><pre><a href="javascript:_e(%7B%7D,'cvml','DFTB-Plus-User@mailman.zfn.uni-bremen.de');" target="_blank">DFTB-Plus-User@mailman.zfn.uni-bremen.de</a><o:p></o:p></pre><pre><a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" target="_blank">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a><o:p></o:p></pre></blockquote><pre>-- <o:p></o:p></pre><pre>Dr. Gabriele Penazzi<o:p></o:p></pre><pre>BCCMS - University of Bremen<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><a href="http://www.bccms.uni-bremen.de/" target="_blank">http://www.bccms.uni-bremen.de/</a><o:p></o:p></pre><pre><a href="http://sites.google.com/site/gabrielepenazzi/" target="_blank">http://sites.google.com/site/gabrielepenazzi/</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>phone: +49 (0) 421 218 9328<o:p></o:p></pre><pre>fax: +49 (0) 421 218 4764<o:p></o:p></pre></div></div></div></body></html>