The ABINIT Users Mailing List ( CLOSED )

[matthieu verstraete <matthieu.jean.verstraete@gmail.com>]

Hello Markus,

1) have you checked that you are only looking at the real part of the
wf? As indicated, by default you get both real and imaginary parts,
and abinit has no particular preference for real valued wf for finite
systems. You should probably plot the norm squared, using both values,
or plot both values side by side.

2) are you sure you have the correct band? This looks like a high
lying box state with all the weight at the edges

3) the plot looks like it is not centered properly: are you certain
about the position of the molecule (around 0,0,0 I presume) and the
coordinates of the box you plot? If the molecule is actually in a
corner of the box this makes more sense, and suggests the box is not
big enough as the density in the middle is not very small...

4) try the xcrysden format. That gives you the norm squared
automatically, and you can visualize it more easily, without messing
around with scripts etc... (option 9 in the main cut3d menu for wf).
In this case, if your molecule is not centered, you can shift the data
by a certain number of grid points to center it.

Cheers

Matthieu


On Tue, Aug 4, 2009 at 7:56 PM, Markus
Wohlgenannt<markus-wohlgenannt@uiowa.edu> wrote:
> Dear forum,
>
> I am a new user to abinit and have a question. I hope someone can help me.
>
> I want to plot the real space wavefunction. I was trying to use cut3d. Cut3d
> works well for plotting the density. I attached the plot of the density in
> the molecular plane for tutorial one, H2. But when I try to plot the
> real-space wavefunction, I seem to get nonsensical results, the
> corresponding wavefunction plot is shown as the second graph in my
> attachement. I did the following to obtain this plot. I start cut3d, and
> tell it to open the WFK file. Then I choose option "0" to analyze the
> wavefunction, choose band 1, I do not want atomic analysis (option 0), then
> I choose option 4 "3D indexed real and imaginary data".
>
> Then I manually select the molecular plane by deleting all data outside the
> plane I want to view, i.e. all data for which z neq Z, where I choose Z
> half-way between the minum and maximum z to obtain the molecular plane.
>
> What am I doing wrong? Or else, is there an easier way to plot the
> real-space wavefunction?
>
> Someone please help,
>
> Thanks,
>
> Markus
>
> --
> Markus Wohlgenannt
> Associate Professor
> Department of Physics and Astronomy
> University of Iowa
> markus-wohlgenannt@uiowa.edu
> 319-353-1974
>
>



-- 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Matthieu Verstraete

European Theoretical Spectroscopy Facility (ETSF)
Dpto. Fisica de Materiales,
U. del Pais Vasco,
Centro Joxe Mari Korta, Av. de Tolosa, 72,   Phone: +34-943018393
E-20018 Donostia-San Sebastian, Spain        Fax  : +34-943018390

Mail : matthieu.jean.verstraete@gmail.com
http://www-users.york.ac.uk/~mjv500