<div dir="ltr">Dear Fengchao,<div><br></div><div>Were you able to find the answer to your question? If yes, then will you please share it. I am also trying to obtain the total electron density of my system which is finite (non-periodic).</div><div><br></div><div>Thanks,</div><div>Megha<br><br>On Friday, June 20, 2014 at 6:28:16 AM UTC-4, wolfshow wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear ALL<br><div>Today, I just read these message posted 6 years ago. Thanks for all your discussions. I would like to ask if it is possible to get the total electron density in CP2K now?</div><div><br></div><div> Any suggestions or comments would be highly appreciated.</div><div><br></div><div>Fengchao</div><div><br><br>在 2008年3月16日星期日UTC+8上午4时11分16秒，<wbr>Axel写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>hi matthias,
<br>
<br>[...]
<br>
<br>> OK, I will unitize the X-ray routine in such a way that one can
<br>> request the build of the all-electron density in Fourier space
<br>> independently.
<br>
<br>
<br>cool. thanks.
<br>
<br>[...]
<br>
<br>> > in the particularl analysis that i have in mind.
<br>>
<br>> > an alternate option would be to incorporate the analysis directly
<br>> > into cp2k (i'm currently thinking of the code from <a href="http://theory.cm.utexas.edu/bader/" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ftheory.cm.utexas.edu%2Fbader%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGjU16prVtQ4Fe_koiWn93D7ZqSeg';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ftheory.cm.utexas.edu%2Fbader%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGjU16prVtQ4Fe_koiWn93D7ZqSeg';return true;">http://theory.cm.utexas.edu/<wbr>bader/</a> whichlooks like it
<br>> > could be merged with rather moderate effort) and thus at least
<br>> > avoid the issue of having to read and write huge formatted files.
<br>>
<br>> A long time ago there was already a request to apply the Bader
<br>> analysis for wavefunctions obtained with CP2K/QS. Inter alia for this
<br>> reason CP2K provides the optional printout of the orbitals in
<br>> Cartesian representation. However, that is a pretty inefficient way to
<br>> interface two codes as you correctly noted, since one has to deal with
<br>> huge formatted files and there are also pitfalls due to the code
<br>> specific ordering of the atomic orbitals .. in a word: awful. It would
<br>> be better to pack the Bader code into one module with a clean and
<br>> simple interface, if this is possible and allowed (licence?).
<br>
<br>i contacted the group in texas and they are very supportive.
<br>they may even consider _using_ cp2k in the future... ;-)
<br>
<br>cheers,
<br>    axel.
<br>>
<br>> cheers,
<br>>
<br>> Matthias
<br>>
<br>> > cheers,
<br>> >    axel.
<br>>
<br>> > > direct collocation of the gaussians (and also stride as now done) do
<br>> > > not conserve the total density when the cutoff is not high enough
<br>> > > (basically always with GAPW).
<br>>
<br>> > > ciao
<br>> > > Fawzi
<br>>
<br>> > > > cheers,
<br>>
<br>> > > > Matthias
<br>>
<br>> > > > On Mar 13, 2:30 am, Axel <<a>akoh...@gmail.com</a>> wrote:
<br>> > > >> hi everybody,
<br>>
<br>> > > >> would there be a way to get the total _electron_ density from
<br>>
<br>> > > >> quickstep GPW/GAPW calculations?
<br>>
<br>> > > >> i actually found a flag to output a 'total density', but that turned
<br>> > > >> out to be the valence electron
<br>> > > >> density plus the (smeared out?) core charge density, with the odd
<br>> > > >> side
<br>> > > >> effect that the cores were
<br>> > > >> negative and the electrons positive. i guess it does not really
<br>> > > >> matter, but...
<br>>
<br>> > > >> for all-electron GAPW calculations, this information should - in
<br>> > > >> principle - be available.
<br>> > > >> the question is, would it be much effort to put it into a cube file?
<br>>
<br>> > > >> for pseudopotential based calculations, be it GPW or GAPW one could
<br>> > > >> just use the
<br>> > > >> density that was "removed" during pseudopotential generation.
<br>>
<br>> > > >> any comments or suggestions would be highly appreciated.
<br>>
<br>> > > >> cheers,
<br>> > > >>    axel.</blockquote></div></div></blockquote></div></div>