Hello,
Thanks! I'll work with the old espresso in the meantime, and also try out the quick fix you suggested.
Thanks again,
kopinjol
iotk issue in yambo-3.2.2-rev.5
Moderators: andrea.ferretti, Conor Hogan
Forum rules
Before posting, make sure to first read carefully all the detailed steps and suggestions given in the main documentation:
https://www.yambo-code.eu/learn/#learn-forum
Before posting, make sure to first read carefully all the detailed steps and suggestions given in the main documentation:
https://www.yambo-code.eu/learn/#learn-forum
-
- Posts: 44
- Joined: Thu Sep 23, 2010 6:59 pm
Re: iotk issue in yambo-3.2.2-rev.5
Kopinjol Baishya
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
- Conor Hogan
- Posts: 111
- Joined: Tue Mar 17, 2009 12:17 pm
- Contact:
Re: iotk issue in yambo-3.2.2-rev.5
By the way, Kopinjol,
you should be aware that the p2y interface does NOT compute the [Vnl,r] terms needed to compute the dipole matrix elements properly in the presence of a nonlocal pseudopotential. Although this warning appears all over the place when you run the code, as well as stressed in the p2y documentation, I'm afraid that people ignore it. Indeed for typical semiconductor systems the correction is negligible. However, for Ag I'm sure it is quite large. Unless you know better, you should do some tests using Abinit and a2y, in my opinion!
Conor
you should be aware that the p2y interface does NOT compute the [Vnl,r] terms needed to compute the dipole matrix elements properly in the presence of a nonlocal pseudopotential. Although this warning appears all over the place when you run the code, as well as stressed in the p2y documentation, I'm afraid that people ignore it. Indeed for typical semiconductor systems the correction is negligible. However, for Ag I'm sure it is quite large. Unless you know better, you should do some tests using Abinit and a2y, in my opinion!
Conor
Dr. Conor Hogan
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
-
- Posts: 44
- Joined: Thu Sep 23, 2010 6:59 pm
Re: iotk issue in yambo-3.2.2-rev.5
Hello Conor,
Thanks for pointing out the issue with the dipole matrix elements. I'll run some tests with abinit and check my results. In the meantime I have recompiled yambo(3.2.2-rev5) after commenting out the line you had suggested("call check_pseudo(n_atomic_species,psfile_)"). But I still get an error when I run p2y. It has something to do with reading bands:
<---> [01] P(W) 2 Y(ambo) Ver. 4.0
<---> DBs path set to .
<---> Index file set to data-file.xml
<---> Header/K-points/Energies...done
<---> Cell data...done
<01s> Atomic data...done
<01s> Symmetries...[SI yes]......[-I yes]...[TR no]
<01s> [WARNING]Unknown XC functional: swicthing to PZ. PW XC string is: PBE
<01s> XC functional...Perdew & Zunger (xc)
<01s> K-points mesh...done
<01s> RL vectors...done
<01s> IGK arrays...done
<01s> Energies...
[ERROR] STOP signal received while in :[01] P(W) 2 Y(ambo) Ver. 4.0
[ERROR]Error in qexml_read_bands IOTK error ierr: 1
FORTRAN STOP
I have checked my calculations and rerun them several times. Do you think there is still a bug in the code or there is something i am doing wrong?
Thanks again,
kopinjol
Thanks for pointing out the issue with the dipole matrix elements. I'll run some tests with abinit and check my results. In the meantime I have recompiled yambo(3.2.2-rev5) after commenting out the line you had suggested("call check_pseudo(n_atomic_species,psfile_)"). But I still get an error when I run p2y. It has something to do with reading bands:
<---> [01] P(W) 2 Y(ambo) Ver. 4.0
<---> DBs path set to .
<---> Index file set to data-file.xml
<---> Header/K-points/Energies...done
<---> Cell data...done
<01s> Atomic data...done
<01s> Symmetries...[SI yes]......[-I yes]...[TR no]
<01s> [WARNING]Unknown XC functional: swicthing to PZ. PW XC string is: PBE
<01s> XC functional...Perdew & Zunger (xc)
<01s> K-points mesh...done
<01s> RL vectors...done
<01s> IGK arrays...done
<01s> Energies...
[ERROR] STOP signal received while in :[01] P(W) 2 Y(ambo) Ver. 4.0
[ERROR]Error in qexml_read_bands IOTK error ierr: 1
FORTRAN STOP
I have checked my calculations and rerun them several times. Do you think there is still a bug in the code or there is something i am doing wrong?
Thanks again,
kopinjol
You do not have the required permissions to view the files attached to this post.
Kopinjol Baishya
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
- Conor Hogan
- Posts: 111
- Joined: Tue Mar 17, 2009 12:17 pm
- Contact:
Re: iotk issue in yambo-3.2.2-rev.5
Are you trying to do a spin-polarized calculation? Unfortunately the GPL version of Yambo doesn't yet support spin...
In any case, the code should exit gracefully..could you please post your pwscf input file and I'll have a look at it?
cheers
Conor
In any case, the code should exit gracefully..could you please post your pwscf input file and I'll have a look at it?
cheers
Conor
Dr. Conor Hogan
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
-
- Posts: 44
- Joined: Thu Sep 23, 2010 6:59 pm
Re: iotk issue in yambo-3.2.2-rev.5
Hello,
Yes my calculations are actually spin polarized. I didn't realize that the GPL version of yambo doesn't do spin polarized systems. Is there a version of yambo that does? The input files are attached.
Thanks again,
kopinjol
Yes my calculations are actually spin polarized. I didn't realize that the GPL version of yambo doesn't do spin polarized systems. Is there a version of yambo that does? The input files are attached.
Thanks again,
kopinjol
You do not have the required permissions to view the files attached to this post.
Kopinjol Baishya
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
- Conor Hogan
- Posts: 111
- Joined: Tue Mar 17, 2009 12:17 pm
- Contact:
Re: iotk issue in yambo-3.2.2-rev.5
Hi again, I'm afraid you are out of luck...spin is something we have been working on but it is certainly not in any shape to be released (i.e. it has not been tested or debugged properly). I might direct you to Abinit, which has (AFAIK) spin support up to the GW level.kopinjol wrote:Yes my calculations are actually spin polarized. I didn't realize that the GPL version of yambo doesn't do spin polarized systems. Is there a version of yambo that does?
Sorry!
Conor
Dr. Conor Hogan
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
CNR-ISM, via Fosso del Cavaliere, 00133 Roma, Italy;
Department of Physics and European Theoretical Spectroscopy Facility (ETSF),
University of Rome "Tor Vergata".
- andrea marini
- Posts: 325
- Joined: Mon Mar 16, 2009 4:27 pm
- Contact:
Re: iotk issue in yambo-3.2.2-rev.5
The p2y interface has been LARGELY updated in the latest revision (rev. 15). Can you update your source and test that it actually works.
THX!!!!
Andrea
THX!!!!
Andrea
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)
Istituto di Struttura della Materia, CNR, (Italy)
-
- Posts: 44
- Joined: Thu Sep 23, 2010 6:59 pm
Re: iotk issue in yambo-3.2.2-rev.5
Hello,
kopinjol
Thanks, I'll check out abinit. Thanks for letting me know.Conor Hogan wrote: I might direct you to Abinit, which has (AFAIK) spin support up to the GW level.
I'll update my source and test it during the weekend and let you know. Thanks for the reply.andrea marini wrote:The p2y interface has been LARGELY updated in the latest revision (rev. 15). Can you update your source and test that it actually works.
kopinjol
Kopinjol Baishya
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607
Graduate Student,
Materials Modeling Group
Deptt of Physics,
University of Illinois,
Chicago,
IL 60607