More on hexagonal cells

Concerns any physical issues arising during the setup step (-i option). This includes problems with symmetries, k/q-point sets, and so on. For technical problems (running in parallel, etc), refer to the Technical forum.

Moderator: Daniele Varsano

More on hexagonal cells

Postby abotello » Mon Apr 27, 2009 5:53 pm

Dear YAMBO Developers,

I have recently installed YAMBO 3.2.1 r.392 in a linux PC.
details of installation are:
./configure --without-mpi --with-iotk=/opt/pwscf/iotk --with-p2y=4.0 --with-netcdf_include=/usr/local/include/ --with-netcdf_lib=/usr/local/lib --enable-largedb

I went through most of the tutorials with no problems. Now I'm trying to work with my system of interest: ZnO. I tried a simple bulk wurtzite calculation. I have read in previous posts that shifted hexagonal cells generated with the automatic option in pwscf leads to an error because it generates a non uniform k-grid. However, I tried initializating the wurtzite ZnO system with YAMBO using a 10x10x10 grid with no shift and leads the following:

<---> [01] Job Setup
<---> [02] Input variables setup
<---> [02.01] K-grid lattice
<---> [02.02] RL shells
<15s> Shells finder |####################| [100%] 15s(E) 15s(X)
<15s> [02.03] Input (E)nergies[ev] & Occupations
<15s> [03] Transferred momenta grid
<15s> [RL indx] X grid is not uniform. Gamma point only.
<15s> [04] Game Over & Game summary

Am I missing something? I attached the pwscf imputs.

Thank you and best wishes,
Andres.
You do not have the required permissions to view the files attached to this post.
Andres R. Botello-Mendez
PhD Candidate
Advanced Materials Department
IPICYT, Mexico
abotello@ipicyt.edu.mx
www.ipicyt.edu.mx
abotello
 
Posts: 1
Joined: Thu Apr 23, 2009 4:19 am
Location: Mexico

Re: More on hexagonal cells

Postby andrea marini » Mon Apr 27, 2009 9:51 pm

Dear Andres,

I have correctly reproduced your error. I will work on it in the next days.

Thanks for pointing out this problem !
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)
User avatar
andrea marini
 
Posts: 447
Joined: Mon Mar 16, 2009 4:27 pm

Re: More on hexagonal cells

Postby andrea marini » Thu Apr 30, 2009 11:21 am

abotello wrote:
Am I missing something? I attached the pwscf inputs.

Thank you and best wishes,
Andres.


Dear Andres,

using the calculation mode 'bands' I indeed get your error with the 10x10x10 grid. Nevertheless if I perform a 'scf' or an 'ncsf' calculation with the same grid Yambo does not complain.

By comparing the output of pw.x in the two calculation modes (attached to this post) I noticed that in the case of calculation='bands' pw.x generates 84 k-points, that are 116 when calculation='scf' or 'nscf'.

To my understanding in the 'bands' mode pw.x does not expand the k-points in the whole irreducible wedge of the brillouin zone. I am not an expert of PW, so I do not know why it does so.

The easy solution is to run pw.x with calculation='nscf' instead of 'bands'.

Andrea
You do not have the required permissions to view the files attached to this post.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)
User avatar
andrea marini
 
Posts: 447
Joined: Mon Mar 16, 2009 4:27 pm


Return to Initialization

Who is online

Users browsing this forum: No registered users and 1 guest

cron