From mdewergifosse@gmail.com Tue Apr 3 09:30:45 2007
From: mdewergifosse@gmail.com (Marc de Wergifosse)
Date: Tue, 3 Apr 2007 10:30:45 +0200
Subject: [Wannier] problem with slab calculation
Message-ID: <3bfb6a350704030130m69e198b6ica06d4e49ea6d497@mail.gmail.com>
------=_Part_4045_9894986.1175589045643
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
Hi,
I'm using wannier90 and pwscf for calculation on a silicon slab.
I need a big sampling of k points to get the Hamiltonian on the Wannier
functions with a good precision.
I have done my nscf calculation on a grid of 24 24 1 k points, which I get
with the script kmesh.pl.
When I want to do the -pp operation with wannier90 on the same grid, I
obtain this error :
*---------------------------------- K-MESH
----------------------------------*
+-------------------------------------------------------------------------=
---+
| Distance to Nearest-Neighbour
Shells |
|
------------------------------------ |
| Shell Distance (Ang^-1)
Multiplicity |
| ----- -----------------
------------ |
| 1 0.024287
4 |
| 2 0.034347
4 |
| 3 0.048574
4 |
| 4 0.054307
8 |
| 5 0.068694
4 |
| 6 0.072861
4 |
| 7 0.076802
8 |
| 8 0.087568
8 |
| 9 0.097148
4 |
| 10 0.100137
8 |
| 11 0.103041
4 |
| 12 0.108614
8 |
+-------------------------------------------------------------------------=
---+
| The b-vectors are chosen
automatically |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
| SVD found small singular value, Rejecting this shell and trying the
next |
Exiting.......
kmesh_get_automatic: Unable to satisfy B1 with any of the first 12 shells
I tried with 32 32 1, 16 16 1 the same -pp operation and I get the same
error.
Anybody can help me?
Thanks for the help,
Marc
________________________________
Marc de Wergifosse
Physique de la mati=E8re condens=E9e
Universit=E9 de Li=E8ge
------=_Part_4045_9894986.1175589045643
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
Hi,
I'm using wannier90 and pwscf for calculation on a silicon s=
lab.
I need a big sampling of k points to get the Hamiltonian on the Wan=
nier functions with a good precision.
I have done my nscf calculation on=
a grid of 24 24 1 k points, which I get with the script=20
kmesh.pl.
When I want to do the -pp operation with wannier90 on the same=
grid, I obtain this error :
*--------------------------------=
-- K-MESH ----------------------------------*
+-------------------=
---------------------------------------------------------+
| &nbs=
p; Distance to Nearest-Neig=
hbour Shells &nb=
sp; |
|  =
; &n=
bsp; ------------------------------------ &nbs=
p; &=
nbsp; |
| &nb=
sp; Shell =
Distance (Ang^-1) &nb=
sp; Multiplicity  =
; |
| -----&nbs=
p; ------=
----------- ---------=
--- |
|  =
; 1 =
&nb=
sp; 0.024287 &nb=
sp; =
4 =
|
| &nb=
sp; 2 &nbs=
p; 0.034347 &nbs=
p; &=
nbsp; 4 &n=
bsp; |
| &nbs=
p; 3  =
; 0.048574  =
; &n=
bsp; 4 &nb=
sp; |
|  =
; 4 =
0.054307=
&nb=
sp; 8 &nbs=
p; |
| =
5 &=
nbsp; &nbs=
p; 0.068694 &nbs=
p; =
4 |
| &nbs=
p; 6  =
; 0.072861  =
; &n=
bsp; 4 &nb=
sp; |
|  =
; 7 =
0.076802=
&nb=
sp; 8 &nbs=
p; |
| =
8 &=
nbsp; &nbs=
p; 0.087568 &nbs=
p; =
8 |
| &nbs=
p; 9  =
; 0.097148  =
; &n=
bsp; 4 &nb=
sp; |
|  =
; 10  =
; 0.100137 =
; &n=
bsp; 8 &nb=
sp; |
|  =
; 11  =
; &n=
bsp; 0.103041 &n=
bsp; 4 &nb=
sp; |
| &nbs=
p; 12 &nbs=
p; 0.108614 &nbs=
p; &=
nbsp; 8 &n=
bsp; |
+----------------------------------------------=
------------------------------+
| The b-vectors are chosen automat=
ically &nb=
sp; =
|
| SVD found small singular value, Rejecting this shell and trying=
the next |
| SVD found small singular value, Rejectin=
g this shell and trying the next |
| SVD found small s=
ingular value, Rejecting this shell and trying the next |
| SVD found small singular value, Rejecting this shell and trying=
the next |
| SVD found small singular value, Rejectin=
g this shell and trying the next |
| SVD found small s=
ingular value, Rejecting this shell and trying the next |
| SVD found small singular value, Rejecting this shell and trying=
the next |
| SVD found small singular value, Rejectin=
g this shell and trying the next |
Exiting.......
&=
nbsp;kmesh_get_automatic: Unable to satisfy B1 with any of the first 12 she=
lls
I tried with 32 32 1, 16 16 1 the same -pp operation and I get the =
same error.
Anybody can help me?
Thanks for the help,
=
Marc
________________________________
Marc de Wergifosse
Physique de la mati=E8re condens=E9e
Universit=E9 de Li=E8ge
------=_Part_4045_9894986.1175589045643--
From jry20@cam.ac.uk Tue Apr 3 10:11:13 2007
From: jry20@cam.ac.uk (Jonathan Yates)
Date: Tue, 3 Apr 2007 10:11:13 +0100 (BST)
Subject: [Wannier] problem with slab calculation
In-Reply-To: <3bfb6a350704030130m69e198b6ica06d4e49ea6d497@mail.gmail.com>
References: <3bfb6a350704030130m69e198b6ica06d4e49ea6d497@mail.gmail.com>
Message-ID:
On Tue, 3 Apr 2007, Marc de Wergifosse wrote:
> I'm using wannier90 and pwscf for calculation on a silicon slab.
> I need a big sampling of k points to get the Hamiltonian on the Wannier
> functions with a good precision.
> I have done my nscf calculation on a grid of 24 24 1 k points, which I get
> with the script kmesh.pl.
> When I want to do the -pp operation with wannier90 on the same grid, I
> obtain this error :
> Exiting.......
> kmesh_get_automatic: Unable to satisfy B1 with any of the first 12 shells
Hi Marc,
There is a problem finding the b-vectors to satisfy Eqn B1 of the
Marzari, Vanderbilt paper (finite difference formula). You've found that
our routine to choose those vectors is not as robust as it should be!
Wannier90 tries to choose the b-vectors from amongst the first 12 nearest
neighbour shells of kpoints. In your case I suspect the first 12 shells
all lie in the plane of the slab; we'd have to look further out to find
the shell corresponding to the direction perpendicular to the slab.
Could you send me the *.win file, so I can check this is indeed the
problem. I'll send you back a patched version of the code, and we'll
include the changes into the next release.
Thanks
Jonathan
--
Dr Jonathan Yates | Theory of Condensed Matter Group
Corpus Christi College | Cavendish Laboratory
Cambridge, CB2 1RH, UK | Cambridge, CB3 OHE, UK
email jry20@cam.ac.uk | Tel +44 (0)1223 337461
From jaa@ornl.gov Wed Apr 11 01:25:02 2007
From: jaa@ornl.gov (Ashley Alford)
Date: Tue, 10 Apr 2007 20:25:02 -0400 (EDT)
Subject: [Wannier] Posternak calculation without projections
Message-ID:
Hello,
I have a question regarding the original calculation
of Posternak et. al. PRB 65, 184422 on MnO. The article states
that many (~10^4) iterations were needed for the minimization
to converge, and that a couple of different projection orbital
types were tried, producing similar results.
I would like to know if starting with no projections
(use_bloch_phases=TRUE) was attempted, and, if so, were similar
spreads and reasonable centers found.
Thanks for any help anyone can provide.
Cheers,
Ashley
From jry20@cam.ac.uk Wed Apr 25 12:45:29 2007
From: jry20@cam.ac.uk (Jonathan Yates)
Date: Wed, 25 Apr 2007 12:45:29 +0100 (BST)
Subject: [Wannier] Re: questions about wannier
In-Reply-To: <987651.72991.qm@web27613.mail.ukl.yahoo.com>
References: <987651.72991.qm@web27613.mail.ukl.yahoo.com>
Message-ID:
This message is in MIME format. The first part should be readable text,
while the remaining parts are likely unreadable without MIME-aware tools.
--1870892094-2001538273-1177500640=:17177
Content-Type: TEXT/PLAIN; CHARSET=X-UNKNOWN; FORMAT=flowed
Content-Transfer-Encoding: QUOTED-PRINTABLE
Content-ID:
[I'm redirecting a message about Wannier90 sent directly to me as I think=
=20
it is more appropriate for the Wannier mailing list.]
On Wed, 25 Apr 2007, amel la wrote:
> I am running the code in cluster opteron:
>
> I have run scf from pw.x with (mpirun -np 2 ./pw.xscf.out)
> parallel mode
> then nscf from pw.x with (mpirun -np 2 ./pw.xnscf.out)=20
parallel mode
>
> then
> wannier90.x -pp copper (this serial mode)
>
> mpirun -np 2 pw2wannier90.x < copper.pw2wan > pw2wan.out
>
> in order to create "copper.mmn and copper.amn files " but no way to=20
success.
>
> but my problem when I wanted to run pw2wannier for creating
> copper.mmn and copper.amn files I got some error message that I could
> not avoid it many times happen, I am ignoring the reason.
>
> Should I run pw2wannier90.x in serial mode or it has bug?
Amel,
I've just ran example 6 in parallel (2 processors) using QuantumEspresso=
=20
v3.2. It ran with no problems. I'm not aware of any bugs.
1 - Make sure that you have an up-to-date version of QE (=3D>3.2)
2 - Check that it runs correctly in serial
3 - If you have access to a different machine (or a different compiler)=20
check to see if that works. - maybe you have a bug in your compiler.
3 - Report any error messages back to us.
**Without the full error message it's almost impossible for us to help**
> and how about wannier90 code, is it working for paralle or serial ?
Wannier90 is a serial code. It works on small matrices; even for large=20
calculations the execution time should be small.
Yours
Jonathan
On Wed, 25 Apr 2007, amel la wrote:
> Dear Sir,
>
> I red the user guide also the tutorial for wannier code, well the=20
> situation is when I was dealing with example6 for copper in tutorial for=
=20
> wannier code.
>
> I am running the code in cluster opteron:
>
> I have run scf from pw.x with (mpirun -np 2 ./pw.xscf.out)
> parallel mode
> then nscf from pw.x with (mpirun -np 2 ./pw.xnscf.out) para=
llel mode
>
> then
> wannier90.x -pp copper (this serial mode)
>
> mpirun -np 2 pw2wannier90.x < copper.pw2wan > pw2wan.out
>
> in order to create "copper.mmn and copper.amn files " but no way to suc=
cess.
>
> but my problem when I wanted to run pw2wannier for creating
> copper.mmn and copper.amn files I got some error message that I could=20
> not avoid it many times happen, I am ignoring the reason.
>
> Should I run pw2wannier90.x in serial mode or it has bug?
>
> and how about wannier90 code, is it working for paralle or serial ?
>
> I will be very grateful if you could give me some hints in order=20
> getting copper.mmn and copper.amn files.
>
> Also the situation is the same with example 5.
>
> I will be very thankful for your help.
>
> All my best regards,
>
> Amel
>
> University of Oldenburg, Oldenburg,
> Germany
>
>
>
>
> ---------------------------------
> D=E9couvrez une nouvelle fa=E7on d'obtenir des r=E9ponses =E0 toutes vos =
questions ! Profitez des connaissances, des opinions et des exp=E9riences d=
es internautes sur Yahoo! Questions/R=E9ponses.
--=20
Dr Jonathan Yates | Theory of Condensed Matter Group
Corpus Christi College | Cavendish Laboratory
Cambridge, CB2 1RH, UK | Cambridge, CB3 OHE, UK
email jry20@cam.ac.uk | Tel +44 (0)1223 337461
--1870892094-2001538273-1177500640=:17177--