<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p style="margin-top:0; margin-bottom:0">Dear Paolo,</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Thanks for the quick, informative response; I will find other pseudopotentials to work with. And hopefully this thread will provide some use for someone else's work as well.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Wai-Ga David Ho</p>
<br>
<br>
<div style="color:rgb(0,0,0)">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of users-request@lists.quantum-espresso.org <users-request@lists.quantum-espresso.org><br>
<b>Sent:</b> Thursday, February 21, 2019 6:00 AM<br>
<b>To:</b> users@lists.quantum-espresso.org<br>
<b>Subject:</b> users Digest, Vol 139, Issue 22</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Send users mailing list submissions to<br>
users@lists.quantum-espresso.org<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk603634" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
or, via email, send a message with subject or body 'help' to<br>
users-request@lists.quantum-espresso.org<br>
<br>
You can reach the person managing the list at<br>
users-owner@lists.quantum-espresso.org<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Fwd: What are the nsum and the Tr[ns(na)] ?<br>
(Victor augusto Hoyos sinchi)<br>
2. Re: Fwd: What are the nsum and the Tr[ns(na)] ?<br>
(Giuseppe Mattioli)<br>
3. Problem with PDOS output file (Juan Manuel Radear)<br>
4. Re: Problem with PDOS output file (Paolo Giannozzi)<br>
5. [QE users] Fake/Fictional Hydrogen Pseudopotential<br>
(dv009200@fh-muenster.de)<br>
6. Re: [QE users] Fake/Fictional Hydrogen Pseudopotential<br>
(Nicola Marzari)<br>
7. pw2casino mismatch energy error (Wai-Ga Ho)<br>
8. non-convergent scf calculation on organometal perovskite<br>
structure (Julien Barbaud)<br>
9. Re: non-convergent scf calculation on organometal perovskite<br>
structure (Fabrizio Cossu)<br>
10. "Basin Hopping" algorithm for pwscf? (Christoph Wolf)<br>
11. Re: pw2casino mismatch energy error (Paolo Giannozzi)<br>
12. Re: non-convergent scf calculation on organometal perovskite<br>
structure (Pietro Davide Delugas)<br>
13. Re: pw2casino mismatch energy error (Paolo Giannozzi)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 20 Feb 2019 08:46:13 -0500<br>
From: Victor augusto Hoyos sinchi <victor23augusto03@gmail.com><br>
To: users@lists.quantum-espresso.org<br>
Subject: [QE-users] Fwd: What are the nsum and the Tr[ns(na)] ?<br>
Message-ID:<br>
<CAPuG_F60FVe5c2dG+OhDy9--GuviddzRe0N_twU3fmizpWat0A@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
---------- Mensaje reenviado ---------<br>
De: Victor augusto Hoyos sinchi <victor23augusto03@gmail.com><br>
Fecha: mi?., 20 de feb. de 2019 a la (s) 08:42<br>
Asunto: ?Cu?les son el nsum y el Tr [ns (na)]?<br>
Para: <users@lists.quantum-espresso.org><br>
<br>
<br>
Estoy haciendo un c?lculo de scf y en la salida tengo un t?rmino Tr [ns<br>
(na)]. ?Qu? es este t?rmino ?. Gracias.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_3a333eaa_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=a0ia8FDwRkJ8zR8kGwn7-HK85qiJGdOFiq6KyTS9K94&e=" id="LPlnk216806" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_3a333eaa_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=a0ia8FDwRkJ8zR8kGwn7-HK85qiJGdOFiq6KyTS9K94&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 20 Feb 2019 14:51:53 +0100<br>
From: Giuseppe Mattioli <giuseppe.mattioli@ism.cnr.it><br>
To: users@lists.quantum-espresso.org<br>
Subject: Re: [QE-users] Fwd: What are the nsum and the Tr[ns(na)] ?<br>
Message-ID:<br>
<20190220145153.Horde.B9sdh1OUVx7gE5w8SNWZxk6@webmail.sic.rm.cnr.it><br>
Content-Type: text/plain; charset=utf-8; format=flowed; DelSp=Yes<br>
<br>
<br>
Dear User<br>
Are you sure you're writing to the correct address? ;-D<br>
Best<br>
Giuseppe<br>
<br>
Quoting Victor augusto Hoyos sinchi <victor23augusto03@gmail.com>:<br>
<br>
> ---------- Mensaje reenviado ---------<br>
> De: Victor augusto Hoyos sinchi <victor23augusto03@gmail.com><br>
> Fecha: mi?., 20 de feb. de 2019 a la (s) 08:42<br>
> Asunto: ?Cu?les son el nsum y el Tr [ns (na)]?<br>
> Para: <users@lists.quantum-espresso.org><br>
><br>
><br>
> Estoy haciendo un c?lculo de scf y en la salida tengo un t?rmino Tr [ns<br>
> (na)]. ?Qu? es este t?rmino ?. Gracias.<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <giuseppe.mattioli@ism.cnr.it><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 20 Feb 2019 09:34:40 -0600<br>
From: Juan Manuel Radear <radear82@gmail.com><br>
To: users@lists.quantum-espresso.org<br>
Subject: [QE-users] Problem with PDOS output file<br>
Message-ID:<br>
<CAKw47gy698AQXGz+T-N-v1JSNaUTEa9Vah1bBNFRuLSerJ9sCQ@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear all:<br>
I'm working on a perovskite system to practice doing PDOS and band<br>
structure calculations. In the past, I've been able to calculate PDOS for<br>
2D systems with success. But in this case I have not being able to<br>
calculate correctly the PDOS. After doing a scf calculation and then the<br>
projwfc.x calculation, I obtain the following for the PDOS output file<br>
(only the first lines copied):<br>
# E (eV) dos(E) pdos(E)<br>
******* 0.473E+01 0.447E+01<br>
******* 0.373E+01 0.352E+01<br>
******* 0.290E+01 0.274E+01<br>
******* 0.224E+01 0.212E+01<br>
******* 0.171E+01 0.162E+01<br>
...<br>
As you can see, the problem is that the first column is not shown. I have<br>
experimented with different variations of the input file. This is what I've<br>
tried:<br>
<br>
1. I have tried a scf calculation and then the projwfc calculation<br>
2. I've tried scf, nscf and then the projwfc one.<br>
3. I've tried relaxing the structure (originally obtained from the COD<br>
database <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.crystallography.net_cod_&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=_8k7szD6VIk5XiWh1LhCD9dB4VfJm1cXMF-geSbLrVI&e=" id="LPlnk788052" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__www.crystallography.net_cod_&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=_8k7szD6VIk5XiWh1LhCD9dB4VfJm1cXMF-geSbLrVI&e=</a>
and working with the<br>
relaxed structure.<br>
4. I've tried using the original structure as obtained from the COD<br>
database.<br>
5. I've tried different values of ecut and the k-point grid.<br>
6. I've tried changing the DeltaE parameter because I thought it was a<br>
matter of the number if significant figures allowed in that column.<br>
<br>
The Fermi Energy that I obtain in the scf calculation is always around -183<br>
eV<br>
If anyone could give some insight on what could be done to correct this<br>
error, I would appreciate it immensely. Here is the scf file of my most<br>
recent try:<br>
--------------------------------------------------<br>
&CONTROL<br>
calculation='scf',<br>
outdir='./tmp',<br>
prefix='CsGeI3-v7',<br>
pseudo_dir='./',<br>
verbosity='high',<br>
nstep = 2800,<br>
/<br>
<br>
&SYSTEM<br>
ibrav=4,<br>
celldm(1)=15.7947296813d0, celldm(3)=1.2693882698d0,<br>
nat=15,<br>
ntyp=3,<br>
ecutwfc=30,<br>
occupations='smearing',<br>
smearing='mv',<br>
degauss=0.005d0,<br>
/<br>
<br>
&ELECTRONS<br>
diagonalization='david',<br>
mixing_mode = 'local-TF',<br>
mixing_beta = 0.1,<br>
conv_thr = 1.0d-5,<br>
electron_maxstep= 2000,<br>
/<br>
<br>
ATOMIC_SPECIES<br>
Cs 132.905000d0 Cs.pbe-spnl-kjpaw_psl.1.0.0.UPF<br>
Ge 72.640000d0 Ge.pbe-dn-kjpaw_psl.0.2.2.UPF<br>
I 126.904000d0 I.pbe-n-kjpaw_psl.0.2.UPF<br>
<br>
ATOMIC_POSITIONS {crystal}<br>
I 0.5027350000d0 0.4972750000d0 0.0042000000d0<br>
I 0.5027250000d0 0.0054600000d0 0.0042000000d0<br>
I 0.9945400000d0 0.4972650000d0 0.0042000000d0<br>
I 0.1694016700d0 0.8306083300d0 0.3375333300d0<br>
I 0.1693916700d0 0.3387933300d0 0.3375333300d0<br>
I 0.6612066700d0 0.8305983300d0 0.3375333300d0<br>
I 0.8360683300d0 0.1639416700d0 0.6708666700d0<br>
I 0.8360583300d0 0.6721266700d0 0.6708666700d0<br>
I 0.3278733300d0 0.1639316700d0 0.6708666700d0<br>
Cs 0.0000000000d0 0.0000000000d0 0.0100010000d0<br>
Cs 0.6666666700d0 0.3333333300d0 0.3433343300d0<br>
Cs 0.3333333300d0 0.6666666700d0 0.6766676700d0<br>
Ge 0.0000000000d0 0.0000000000d0 0.5394317000d0<br>
Ge 0.6666666700d0 0.3333333300d0 0.8727650300d0<br>
Ge 0.3333333300d0 0.6666666700d0 0.2060983700d0<br>
<br>
K_POINTS {automatic}<br>
4 4 4 0 0 0<br>
---------------------------------------------------------<br>
And here is the PDOS input file:<br>
--------------------------------------------------------<br>
&PROJWFC<br>
outdir='./tmp',<br>
prefix='CsGeI3-v7',<br>
Emin = -198,<br>
Emax = -168,<br>
degauss = 0.01,<br>
DeltaE=0.01,<br>
filpdos = 'CsGeI3-v7.pdos',<br>
/<br>
-----------------------------------------<br>
As I mentioned before, the Fermi E is around -183, and I chose the Emin and<br>
Emax around that value, as I have seen in different tutorials. I have also<br>
tried Delta E= 0.1 to try and reduce the number of significant figures<br>
shown in the output file but with no success.<br>
Thanks again in advance for your attentions and sorry for the<br>
inconveniences.<br>
<br>
Dr. Juan Manuel Ramirez de Arellano<br>
Profesor de C?tedra, Departamento de Ciencias,<br>
Tecnologico de Monterrey, Calle del Puente 222, 14380, Tlalpan,<br>
Mexico<br>
Email: jramirezdearellano@tec.mx<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_ecc2e14d_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=MSNqIf8tjJQYq-5WTqi3TIsBKLd8RJ3qIM7QbyGytKE&e=" id="LPlnk236560" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_ecc2e14d_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=MSNqIf8tjJQYq-5WTqi3TIsBKLd8RJ3qIM7QbyGytKE&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 20 Feb 2019 16:45:55 +0100<br>
From: Paolo Giannozzi <p.giannozzi@gmail.com><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
Subject: Re: [QE-users] Problem with PDOS output file<br>
Message-ID:<br>
<CAPMgbCu-opTFt9YS_Z1Fc8L2EOmn1F_BP-vS+3wgrGXdORhs5A@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
On Wed, Feb 20, 2019 at 4:35 PM Juan Manuel Radear <radear82@gmail.com><br>
wrote:<br>
<br>
# E (eV) dos(E) pdos(E)<br>
> ******* 0.473E+01 0.447E+01<br>
><br>
[...] As you can see, the problem is that the first column is not shown<br>
> [...] The Fermi Energy that I obtain in the scf calculation is always<br>
> around -183 eV<br>
><br>
<br>
the format used for the energy is f7.3, so energies around -183 eV will<br>
simply not fit into it (7 figures - 3 decimals - 1 dot - minus sign = 2<br>
figures left). Just increase it to f8.4:<br>
<br>
$ grep f7.3 PP/src/partialdos.f90<br>
WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, &<br>
WRITE (4,'(f7.3,4e11.3)') etev*rytoev, (dostot(ie,is,ik),<br>
is=1,nspin), &<br>
WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev,<br>
ldos(ie,1,ik), &<br>
WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev,<br>
ldos(ie,1,ik), &<br>
WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, &<br>
WRITE (4,'(f7.3,4e11.3)') etev*rytoev, dostot(ie,ik), &<br>
<br>
Paolo<br>
<br>
If anyone could give some insight on what could be done to correct this<br>
> error, I would appreciate it immensely. Here is the scf file of my most<br>
> recent try:<br>
> --------------------------------------------------<br>
> &CONTROL<br>
> calculation='scf',<br>
> outdir='./tmp',<br>
> prefix='CsGeI3-v7',<br>
> pseudo_dir='./',<br>
> verbosity='high',<br>
> nstep = 2800,<br>
> /<br>
><br>
> &SYSTEM<br>
> ibrav=4,<br>
> celldm(1)=15.7947296813d0, celldm(3)=1.2693882698d0,<br>
> nat=15,<br>
> ntyp=3,<br>
> ecutwfc=30,<br>
> occupations='smearing',<br>
> smearing='mv',<br>
> degauss=0.005d0,<br>
> /<br>
><br>
> &ELECTRONS<br>
> diagonalization='david',<br>
> mixing_mode = 'local-TF',<br>
> mixing_beta = 0.1,<br>
> conv_thr = 1.0d-5,<br>
> electron_maxstep= 2000,<br>
> /<br>
><br>
> ATOMIC_SPECIES<br>
> Cs 132.905000d0 Cs.pbe-spnl-kjpaw_psl.1.0.0.UPF<br>
> Ge 72.640000d0 Ge.pbe-dn-kjpaw_psl.0.2.2.UPF<br>
> I 126.904000d0 I.pbe-n-kjpaw_psl.0.2.UPF<br>
><br>
> ATOMIC_POSITIONS {crystal}<br>
> I 0.5027350000d0 0.4972750000d0 0.0042000000d0<br>
> I 0.5027250000d0 0.0054600000d0 0.0042000000d0<br>
> I 0.9945400000d0 0.4972650000d0 0.0042000000d0<br>
> I 0.1694016700d0 0.8306083300d0 0.3375333300d0<br>
> I 0.1693916700d0 0.3387933300d0 0.3375333300d0<br>
> I 0.6612066700d0 0.8305983300d0 0.3375333300d0<br>
> I 0.8360683300d0 0.1639416700d0 0.6708666700d0<br>
> I 0.8360583300d0 0.6721266700d0 0.6708666700d0<br>
> I 0.3278733300d0 0.1639316700d0 0.6708666700d0<br>
> Cs 0.0000000000d0 0.0000000000d0 0.0100010000d0<br>
> Cs 0.6666666700d0 0.3333333300d0 0.3433343300d0<br>
> Cs 0.3333333300d0 0.6666666700d0 0.6766676700d0<br>
> Ge 0.0000000000d0 0.0000000000d0 0.5394317000d0<br>
> Ge 0.6666666700d0 0.3333333300d0 0.8727650300d0<br>
> Ge 0.3333333300d0 0.6666666700d0 0.2060983700d0<br>
><br>
> K_POINTS {automatic}<br>
> 4 4 4 0 0 0<br>
> ---------------------------------------------------------<br>
> And here is the PDOS input file:<br>
> --------------------------------------------------------<br>
> &PROJWFC<br>
> outdir='./tmp',<br>
> prefix='CsGeI3-v7',<br>
> Emin = -198,<br>
> Emax = -168,<br>
> degauss = 0.01,<br>
> DeltaE=0.01,<br>
> filpdos = 'CsGeI3-v7.pdos',<br>
> /<br>
> -----------------------------------------<br>
> As I mentioned before, the Fermi E is around -183, and I chose the Emin<br>
> and Emax around that value, as I have seen in different tutorials. I have<br>
> also tried Delta E= 0.1 to try and reduce the number of significant figures<br>
> shown in the output file but with no success.<br>
> Thanks again in advance for your attentions and sorry for the<br>
> inconveniences.<br>
><br>
> Dr. Juan Manuel Ramirez de Arellano<br>
> Profesor de C?tedra, Departamento de Ciencias,<br>
> Tecnologico de Monterrey, Calle del Puente 222, 14380, Tlalpan,<br>
> Mexico<br>
> Email: jramirezdearellano@tec.mx<br>
><br>
> _______________________________________________<br>
> users mailing list<br>
> users@lists.quantum-espresso.org<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk505436" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
<br>
<br>
<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_09e589a7_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cEcGhLNqI7YMb5MB08K2lgXASOICLMJgINk9mueB6CE&e=" id="LPlnk622182" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190220_09e589a7_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cEcGhLNqI7YMb5MB08K2lgXASOICLMJgINk9mueB6CE&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Wed, 20 Feb 2019 21:52:34 +0100<br>
From: dv009200@fh-muenster.de<br>
To: users@lists.quantum-espresso.org<br>
Subject: [QE-users] [QE users] Fake/Fictional Hydrogen Pseudopotential<br>
Message-ID:<br>
<689ee6a742dc80ae7117b0be775a7565.squirrel@webmail.fh-muenster.de><br>
Content-Type: text/plain;charset=iso-8859-1<br>
<br>
Dear all,<br>
<br>
I'm trying to calculate II-VI semiconductor quantum dots. For this I want<br>
to passivate the dangling bonds of the surface atoms with fake/fictional<br>
hydrogenatoms that have an zval=0.5, 1.5 and 3.0.<br>
<br>
I already found the pseudopotentials for zval= 0.5 to 1.5 on the Quantum<br>
Espresso website<br>
(<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_pseudopotentials_original-2Dqe-2Dpp-2Dlibrary_h&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=CI5o-BQZaaSo6nySuDecZLlBo0dGqQ9r0WjxDaX7LqY&e=" id="LPlnk259031" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_pseudopotentials_original-2Dqe-2Dpp-2Dlibrary_h&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=CI5o-BQZaaSo6nySuDecZLlBo0dGqQ9r0WjxDaX7LqY&e=</a>).<br>
<br>
I want to know how I can generate such fake pseudopotentials by myself<br>
(and for other functionals than the pz) with the ld1.x code? Can someone<br>
provide me an example input file (for ld1.x)? I already tried to<br>
reconstruct the pseudopotentials at the website with the ld1.x code but<br>
failed to do so.<br>
<br>
Thanks and regards<br>
<br>
M.Sc. Dominik Voigt<br>
PhD Student M?nster University of Applied Sciences<br>
Email: dv009200@fh-muenster.de<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Wed, 20 Feb 2019 16:12:30 -0500<br>
From: Nicola Marzari <nicola.marzari@epfl.ch><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org>,<br>
dv009200@fh-muenster.de<br>
Subject: Re: [QE-users] [QE users] Fake/Fictional Hydrogen<br>
Pseudopotential<br>
Message-ID: <42e001ad-bedb-3331-0e6a-8057d17ff15e@epfl.ch><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
<br>
<br>
Dear Dominik,<br>
<br>
<br>
not sure if someone has a ready to go generator of von barth car<br>
pseudos into UPF - but if there is one (if not, see below), and<br>
depending on how much time you want to devote to this, you could<br>
try a mini-project on the best passivator possible - i.e.<br>
generating a pseudopotential that has the smallest possible<br>
influence on your system.<br>
<br>
For this, you could play with valence and core radius of<br>
a spherical Von Barth-Car pseudopotential, and make sure that<br>
substituting a real group II or group VI atom with your passivator<br>
creates the smallest force in the system (this is described, for<br>
a different project, in Appendix D/page 257 of<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__hdl.handle.net_1721.1_28351&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=IVcD1fvfgZPNFhecZx9Bzrqjh8TGO0_oJHpW_1d3bgM&e=" id="LPlnk6671" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__hdl.handle.net_1721.1_28351&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=IVcD1fvfgZPNFhecZx9Bzrqjh8TGO0_oJHpW_1d3bgM&e=</a>)<br>
<br>
Otherwise, the UPF format is described here<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.quantum-2Despresso.org_pseudopotentials_unified-2Dpseudopotential-2Dformat&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cnjtj8-zDOfLevq9wZD0gXxcqC10Ym8lay0YRClCsp8&e=" id="LPlnk541244" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.quantum-2Despresso.org_pseudopotentials_unified-2Dpseudopotential-2Dformat&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cnjtj8-zDOfLevq9wZD0gXxcqC10Ym8lay0YRClCsp8&e=</a><br>
<br>
and you could use the same radial mesh of<br>
<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_upf-5Ffiles_H.pz-2Dvbc-5F150.UPF&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=ZAhBc4TAWIuZBgqJwkkc5bwPILRCcONPVQojFzsbL_M&e=" id="LPlnk475428" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_upf-5Ffiles_H.pz-2Dvbc-5F150.UPF&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=ZAhBc4TAWIuZBgqJwkkc5bwPILRCcONPVQojFzsbL_M&e=</a><br>
<br>
and just add the values of your local von barth car pseudo.<br>
<br>
nicola<br>
<br>
<br>
On 20/02/2019 15:52, dv009200@fh-muenster.de wrote:<br>
> Dear all,<br>
> <br>
> I'm trying to calculate II-VI semiconductor quantum dots. For this I want<br>
> to passivate the dangling bonds of the surface atoms with fake/fictional<br>
> hydrogenatoms that have an zval=0.5, 1.5 and 3.0.<br>
> <br>
> I already found the pseudopotentials for zval= 0.5 to 1.5 on the Quantum<br>
> Espresso website<br>
> (<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_pseudopotentials_original-2Dqe-2Dpp-2Dlibrary_h&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=CI5o-BQZaaSo6nySuDecZLlBo0dGqQ9r0WjxDaX7LqY&e=" id="LPlnk852271" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_pseudopotentials_original-2Dqe-2Dpp-2Dlibrary_h&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=CI5o-BQZaaSo6nySuDecZLlBo0dGqQ9r0WjxDaX7LqY&e=</a>).<br>
> <br>
> I want to know how I can generate such fake pseudopotentials by myself<br>
> (and for other functionals than the pz) with the ld1.x code? Can someone<br>
> provide me an example input file (for ld1.x)? I already tried to<br>
> reconstruct the pseudopotentials at the website with the ld1.x code but<br>
> failed to do so.<br>
> <br>
> Thanks and regards<br>
> <br>
> M.Sc. Dominik Voigt<br>
> PhD Student M?nster University of Applied Sciences<br>
> Email: dv009200@fh-muenster.de<br>
> <br>
> _______________________________________________<br>
> users mailing list<br>
> users@lists.quantum-espresso.org<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk528982" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
> <br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__theossrv1.epfl.ch_Main_Contact&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=PhO_8AntMu0sOU6C7L3rBOp39NqpfLsmlwM1szHwbXI&e=" id="LPlnk337837" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__theossrv1.epfl.ch_Main_Contact&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=PhO_8AntMu0sOU6C7L3rBOp39NqpfLsmlwM1szHwbXI&e=</a>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__nccr-2Dmarvel.ch_en_project&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=Y3nFBc407D5n_iBJYGrSxJDxcCVAudNinRnjC81FE3E&e=" id="LPlnk936454" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__nccr-2Dmarvel.ch_en_project&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=Y3nFBc407D5n_iBJYGrSxJDxcCVAudNinRnjC81FE3E&e=</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Thu, 21 Feb 2019 02:02:00 +0000<br>
From: Wai-Ga Ho <wdh17@my.fsu.edu><br>
To: "users@lists.quantum-espresso.org"<br>
<users@lists.quantum-espresso.org><br>
Subject: [QE-users] pw2casino mismatch energy error<br>
Message-ID:<br>
<SN6PR05MB4078D1DFF476EB59EB6F7529877D0@SN6PR05MB4078.namprd05.prod.outlook.com><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hello QE community,<br>
<br>
<br>
I am currently attempting to extract the kinetic energy contributions to the total energies of different polymorphs of MgO.<br>
<br>
<br>
I understand that we need to use the -pw2casino flag if we want our PWSCF output to print a more detailed breakdown of the total energy in terms of its contributions. This means that using the command:<br>
<br>
<br>
mpirun -np 4 pw.x <project.in> project.out<br>
<br>
<br>
will produce a project.out file which contains some text like this:<br>
<br>
<br>
! total energy = -183.20894402 Ry<br>
Harris-Foulkes estimate = -183.20894409 Ry<br>
estimated scf accuracy < 0.00000009 Ry<br>
<br>
total all-electron energy = -551.326109 Ry<br>
<br>
The total energy is the sum of the following terms:<br>
<br>
one-electron contribution = -85.45759185 Ry<br>
hartree contribution = 47.93001020 Ry<br>
xc contribution = -25.22704816 Ry<br>
ewald contribution = -95.91150866 Ry<br>
one-center paw contrib. = -24.54280555 Ry<br>
<br>
while adding the -pw2casino flag and running:<br>
<br>
mpirun -np 4 pw.x -pw2casino <project.in> project.out<br>
<br>
<br>
will print the original lines above, as well as additional lines like:<br>
<br>
Energies determined by pw2casino tool<br>
-------------------------------------<br>
Kinetic energy 21.442363717309750 au = 42.884727434619499 Ry<br>
Local energy -84.609446621660297 au = -169.21889324332059 Ry<br>
Non-Local energy 20.438286977076853 au = 40.876573954153706 Ry<br>
Ewald energy -47.955754330605572 au = -95.911508661211144 Ry<br>
xc contribution -12.613524077939315 au = -25.227048155878631 Ry<br>
hartree energy 23.965005100370561 au = 47.930010200741123 Ry<br>
Total energy -79.333069235448022 au = -158.66613847089604 Ry<br>
<br>
>From the above output lines, we can see that the total energies (-183.20894402 Ry and -158.66613847089604 Ry) are pretty terribly mismatched. I find that this mismatch occurs when using either PWSCF v5.0.2 and PWSCF v6.0. When I was using PWSCF v6.0, quantum-espresso
actually prints an "Energy mismatch error" (I'd paste the error below but my laptop, which runs v6.0, is currently in the shop), so at least it recognizes that something is wrong.<br>
<br>
I've tried looking through some other threads and attempted to apply fixes that worked for other people, such as adding the following line:<br>
<br>
assume_isolated = 'mp'<br>
<br>
to my &CONTROL in my input file (source: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=" id="LPlnk737990" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=</a>),
but nothing so far has worked.<br>
<br>
Does anyone have any idea what the issue is and how to resolve it? Because with the demonstrated mismatch in total energies, I do not have faith in the kinetic energy value that is output. If this helps, I am using the following pseudopotentials:<br>
<br>
file Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2P renormalized<br>
file O.pbe-nl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2S 2P renormalized<br>
<br>
and below is my input file for MgO in a CsCl-type crystalline arrangement (which generated the energy values I gave above, using PWSCF v5.0.2):<br>
<br>
&CONTROL<br>
calculation = "scf"<br>
max_seconds = 8.64000e+04<br>
outdir = "./"<br>
prefix = "MgO_CsCl"<br>
pseudo_dir = "/home/davidho/.burai/.pseudopot"<br>
title = "Working MgO CsCl Structure(SCF)"<br>
wf_collect = .TRUE.<br>
/<br>
<br>
&SYSTEM<br>
a = 2.65994e+00<br>
degauss = 1.00000e-02<br>
ecutrho = 2.00000e+03<br>
ecutwfc = 2.00000e+02<br>
ibrav = 1<br>
nat = 2<br>
ntyp = 2<br>
occupations = "fixed"<br>
smearing = "gaussian"<br>
assume_isolated = 'mp'<br>
/<br>
<br>
&ELECTRONS<br>
conv_thr = 1.00000e-06<br>
electron_maxstep = 200<br>
mixing_beta = 7.00000e-01<br>
startingpot = "atomic"<br>
startingwfc = "atomic+random"<br>
/<br>
<br>
K_POINTS {automatic}<br>
10 10 10 0 0 0<br>
<br>
ATOMIC_SPECIES<br>
Mg 24.30500 Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF<br>
O 15.99940 O.pbe-nl-kjpaw_psl.1.0.0.UPF<br>
<br>
ATOMIC_POSITIONS {angstrom}<br>
Mg 0.000000 0.000000 0.000000<br>
O 1.329970 1.329970 1.329970<br>
<br>
<br>
Any help and further insight into both the problem with and/or the method of energy decomposition would be greatly appreciated.<br>
<br>
Thank you,<br>
Wai-Ga David Ho, Department of Physics at Florida State University<br>
<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_5030611c_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=5uXIs3X0dP2RVvdQjsvJAKvd6VHkFrEPh2CjSPBLHr4&e=" id="LPlnk82422" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_5030611c_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=5uXIs3X0dP2RVvdQjsvJAKvd6VHkFrEPh2CjSPBLHr4&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Thu, 21 Feb 2019 11:17:55 +0800<br>
From: Julien Barbaud <julien_barbaud@sjtu.edu.cn><br>
To: users@lists.quantum-espresso.org<br>
Subject: [QE-users] non-convergent scf calculation on organometal<br>
perovskite structure<br>
Message-ID: <a55d199b-54b2-09be-4cc5-9826e0151a09@sjtu.edu.cn><br>
Content-Type: text/plain; charset="utf-8"; Format="flowed"<br>
<br>
Dear users,<br>
<br>
<br>
I am new to QE, and trying to run a simple scf calculation on a <br>
CH3NH3PbI3 crystal (semi-conducting material). I am using ultrasoft <br>
pseudopotentials based on the exchange-correlation functionnal PBEsol.<br>
<br>
I set up a first input, with values of parameters inspired from <br>
literature on the subject. However, I could not reach convergence after <br>
100 iterations. The estimated error was actually "exploding" to very <br>
high values, indicating a serious problem. I tried several changes but <br>
was unsuccessful:<br>
<br>
* ?varying plane-wave cutoff energy does not solve the problem (cf<br>
attached ecut.png, giving the estimated error as a function of the<br>
number of iterations. It is shown here only on the first 15<br>
iterations as the results pretty much only stall from there)<br>
* ?varying cutoff energy for charge (cf ecutrho.png)<br>
* taking larger k-point sampling (not shown)<br>
* ?I also read that for metallic or "close to metallic conductors",<br>
there might be problems with the first unoccupied states that can be<br>
solved by adding a few empty bands. My system being a<br>
semi-conductor, I tried adding additional bands using a m-p smearing<br>
but no improvement was found (not shown)<br>
<br>
<br>
The only change that I found effective was to reduce the mixing_beta factor.<br>
<br>
<br>
It effectively prevents the error from diverging to very large values, <br>
but I still do not reach convergence, even after longer iterations. I <br>
tried much smaller values of mixing beta which improves the final value <br>
of the error, but I still cannot reach convergence on 100 iterations. As <br>
shown in the mixbeta2_zoom.png, the error reduces to smaller values <br>
around ~1e-5~1e-6, but it keeps stalling after a while. I do not observe <br>
a well-converging behaviour for any value.<br>
<br>
<br>
I attached the "default version" of my script on which the various <br>
modifications described above have been independently performed. I <br>
obtained the geometry from a CIF file in literature and checked it with <br>
visualization software; it seems perfectly ok as far as I can tell.<br>
<br>
<br>
Any insight on what I did wrong would be really helpful. I suspect a <br>
shameful beginner mistake, but can not find it out.<br>
<br>
<br>
Thanks in advance,<br>
<br>
Julien barbaud<br>
<br>
<br>
P.S: this is my first time posting on this user list. Please let me know <br>
if my question is not suitable for it, or can be improved either in its <br>
content or presentation. I will gladly take any recommandation into <br>
account in order not to negatively impact the quality of this user list !<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=mFLCF2VB8hlgrHPvvIWrwP5IzFAXK66KkpSeMJ8BD1k&e=" id="LPlnk241606" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=mFLCF2VB8hlgrHPvvIWrwP5IzFAXK66KkpSeMJ8BD1k&e=</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: ecut.png<br>
Type: image/png<br>
Size: 32266 bytes<br>
Desc: not available<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0003.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=FfLxWgUO3VyyXA7f1bIj9ZOZt3xITaZCw_SkJgK4_hQ&e=" id="LPlnk597473" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0003.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=FfLxWgUO3VyyXA7f1bIj9ZOZt3xITaZCw_SkJgK4_hQ&e=</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: ecuthro.png<br>
Type: image/png<br>
Size: 22580 bytes<br>
Desc: not available<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0004.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=22Ou90qdqjhShflom1mKAf_bzzApPI1KOfnMq9fT2x4&e=" id="LPlnk33953" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0004.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=22Ou90qdqjhShflom1mKAf_bzzApPI1KOfnMq9fT2x4&e=</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: mixbeta2_zoom_comp.png<br>
Type: image/png<br>
Size: 74603 bytes<br>
Desc: not available<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0005.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=YxVmAkIqHktc7tIS6vWmP_045HBEFLJxtWblCDt_c_w&e=" id="LPlnk332607" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_932c2bbb_attachment-2D0005.png&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=YxVmAkIqHktc7tIS6vWmP_045HBEFLJxtWblCDt_c_w&e=</a>><br>
-------------- next part --------------<br>
&CONTROL<br>
pseudo_dir = "/home/julien_barbaud/QUANTUM_ESPRESSO/Program/qe-6.3/pseudo"<br>
calculation = 'scf'<br>
outdir= "./out"<br>
tstress=.TRUE.<br>
tprnfor=.TRUE.<br>
/<br>
<br>
&SYSTEM<br>
nat= 12<br>
ntyp= 5<br>
ibrav= 0<br>
ecutwfc= 30, ecutrho = 300<br>
/<br>
<br>
&ELECTRONS<br>
mixing_beta = 0.7<br>
conv_thr = 1.0d-8<br>
/<br>
<br>
&IONS<br>
/<br>
<br>
&CELL<br>
/<br>
<br>
ATOMIC_SPECIES<br>
C 12.011 C.pbesol-n-rrkjus_ps1.1.0.0.UPF<br>
N 14.007 N.pbesol-n-rrkjus_ps1.1.0.0.UPF<br>
H 1.008 H.pbesol-rrkjus_ps1.0.1.UPF<br>
Pb 207.2 Pb.pbesol-dn-rrkjus_ps1.1.0.0.UPF<br>
I 126.90 I.pbesol-n-rrkjus_ps1.1.0.0.UPF<br>
<br>
CELL_PARAMETERS angstrom<br>
6.28880000 0.00000000 0.00000000<br>
0.00154263 6.22876981 0.00000000<br>
0.13973295 -0.00036613 6.37255819<br>
<br>
ATOMIC_POSITIONS angstrom<br>
C 5.68762816 6.22746232 2.91755469<br>
N 0.76248596 -0.00015695 3.49068983<br>
H 5.76522586 6.22643528 1.82113056<br>
H 5.15723272 0.90114080 3.25657478<br>
H 5.15968675 5.32511200 3.25829538<br>
H 1.30458379 0.84536565 3.19592078<br>
H 1.30663813 5.38278144 3.19763499<br>
H 0.76813187 0.00090438 4.53227809<br>
Pb 2.93516392 3.11390239 6.22995946<br>
I 2.64234654 3.11409882 3.02788916<br>
I 2.51578467 6.22849954 0.15080022<br>
I 5.99299067 3.11401602 5.87855111<br>
<br>
K_POINTS automatic<br>
4 4 4 0 0 0<br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Thu, 21 Feb 2019 13:53:30 +0900<br>
From: Fabrizio Cossu <fabrizio.cossu@apctp.org><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
Subject: Re: [QE-users] non-convergent scf calculation on organometal<br>
perovskite structure<br>
Message-ID:<br>
<CAHsDu66JCO1qtJzC+J_WStEgkUMVD2QxQZ6ectXdj3PJU2brEg@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Julien,<br>
what smearing did you have in mind? Have you tried raising the broadening?<br>
<br>
Cheers,<br>
Fabrizio<br>
<br>
On Thu, 21 Feb 2019 at 12:19, Julien Barbaud <julien_barbaud@sjtu.edu.cn><br>
wrote:<br>
<br>
> Dear users,<br>
><br>
><br>
> I am new to QE, and trying to run a simple scf calculation on a CH3NH3PbI3<br>
> crystal (semi-conducting material). I am using ultrasoft pseudopotentials<br>
> based on the exchange-correlation functionnal PBEsol.<br>
><br>
> I set up a first input, with values of parameters inspired from literature<br>
> on the subject. However, I could not reach convergence after 100<br>
> iterations. The estimated error was actually "exploding" to very high<br>
> values, indicating a serious problem. I tried several changes but was<br>
> unsuccessful:<br>
><br>
> - varying plane-wave cutoff energy does not solve the problem (cf<br>
> attached ecut.png, giving the estimated error as a function of the number<br>
> of iterations. It is shown here only on the first 15 iterations as the<br>
> results pretty much only stall from there)<br>
> - varying cutoff energy for charge (cf ecutrho.png)<br>
> - taking larger k-point sampling (not shown)<br>
> - I also read that for metallic or "close to metallic conductors",<br>
> there might be problems with the first unoccupied states that can be solved<br>
> by adding a few empty bands. My system being a semi-conductor, I tried<br>
> adding additional bands using a m-p smearing but no improvement was found<br>
> (not shown)<br>
><br>
><br>
> The only change that I found effective was to reduce the mixing_beta<br>
> factor.<br>
><br>
><br>
> It effectively prevents the error from diverging to very large values, but<br>
> I still do not reach convergence, even after longer iterations. I tried<br>
> much smaller values of mixing beta which improves the final value of the<br>
> error, but I still cannot reach convergence on 100 iterations. As shown in<br>
> the mixbeta2_zoom.png, the error reduces to smaller values around<br>
> ~1e-5~1e-6, but it keeps stalling after a while. I do not observe a<br>
> well-converging behaviour for any value.<br>
><br>
><br>
> I attached the "default version" of my script on which the various<br>
> modifications described above have been independently performed. I obtained<br>
> the geometry from a CIF file in literature and checked it with<br>
> visualization software; it seems perfectly ok as far as I can tell.<br>
><br>
><br>
> Any insight on what I did wrong would be really helpful. I suspect a<br>
> shameful beginner mistake, but can not find it out.<br>
><br>
><br>
> Thanks in advance,<br>
><br>
> Julien barbaud<br>
><br>
><br>
> P.S: this is my first time posting on this user list. Please let me know<br>
> if my question is not suitable for it, or can be improved either in its<br>
> content or presentation. I will gladly take any recommandation into account<br>
> in order not to negatively impact the quality of this user list !<br>
> _______________________________________________<br>
> users mailing list<br>
> users@lists.quantum-espresso.org<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk24448" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
<br>
<br>
<br>
-- <br>
<br>
*Fabrizio Cossupostdoctoral fellow at APCTP (Asia Pacific Center for<br>
Theoretical Physics)*,<br>
Hogil Kim Memorial Building #501<br>
POSTECH, 67 Cheongam-Ro, Nam-Gu,<br>
Pohang-si, Gyeongsangbuk-do,<br>
790-784 (37673), Republic of Korea<br>
<br>
|<br>
.----. .----. .----. |---- .----. ===========================<br>
,----| | | | | | | <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.apctp.org_-3FJrgId-3D16&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=oQYZ0NpO69sFyd2SEuM3QF1JG-SvuB0ZVEDZuhORMag&e=" id="LPlnk410704" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__www.apctp.org_-3FJrgId-3D16&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=oQYZ0NpO69sFyd2SEuM3QF1JG-SvuB0ZVEDZuhORMag&e=</a><br>
`----^ |----' `----' `----' |----' ===========================<br>
| |<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_81965124_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=m1LJwhRipMjfxpqygxsnS4D06C_SIr26kVBaSTxW-pg&e=" id="LPlnk166724" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_81965124_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=m1LJwhRipMjfxpqygxsnS4D06C_SIr26kVBaSTxW-pg&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Thu, 21 Feb 2019 14:28:40 +0900<br>
From: Christoph Wolf <wolf.christoph@qns.science><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
Subject: [QE-users] "Basin Hopping" algorithm for pwscf?<br>
Message-ID:<br>
<CAMC_G_5ym_zvVCu0J13yZuuPFo3zEbCyWU02wvMHc6_xP-UjsQ@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear all,<br>
<br>
I saw "basin hopping" referenced in several presentations about pwscf but I<br>
cannot find the corresponding implementation - can anyone point me in<br>
the right direction?<br>
<br>
Thanks in advance!<br>
Chris<br>
<br>
-- <br>
Postdoctoral Researcher<br>
Center for Quantum Nanoscience, Institute for Basic Science<br>
Ewha Womans University, Seoul, South Korea<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_6721029b_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=AUekVzUQuSMngsAUMmyLieYqSR2PXSJsaFXMrVVkh3M&e=" id="LPlnk66431" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_6721029b_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=AUekVzUQuSMngsAUMmyLieYqSR2PXSJsaFXMrVVkh3M&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 11<br>
Date: Thu, 21 Feb 2019 09:07:54 +0100<br>
From: Paolo Giannozzi <p.giannozzi@gmail.com><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
Subject: Re: [QE-users] pw2casino mismatch energy error<br>
Message-ID:<br>
<CAPMgbCtiFjUb1ZcwTXokVoU7RjuOrMPKJ0+73o-gG9zF=gn+eg@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Wai-Ga David,<br>
<br>
the "pw2casino" option has a single specific application: to produce files<br>
for the CASINO Quantum Monte-Carlo code. It can be used for other purposes<br>
as well, but then ... you are on your own. If I remember correctly, it is<br>
not supposed to work for Ultrasoft pseudopotentials and PAW (CASINO works<br>
only for norm-conserving PPs). There is no check on that because users of<br>
CASINO should know.<br>
<br>
Paolo<br>
<br>
On Thu, Feb 21, 2019 at 3:02 AM Wai-Ga Ho <wdh17@my.fsu.edu> wrote:<br>
<br>
> Hello QE community,<br>
><br>
><br>
> I am currently attempting to extract the kinetic energy contributions to<br>
> the total energies of different polymorphs of MgO.<br>
><br>
><br>
> I understand that we need to use the -pw2casino flag if we want our PWSCF<br>
> output to print a more detailed breakdown of the total energy in terms of<br>
> its contributions. This means that using the command:<br>
><br>
> mpirun -np 4 pw.x <project.in> project.out<br>
><br>
><br>
> will produce a project.out file which contains some text like this:<br>
><br>
><br>
> ! total energy = -183.20894402 Ry<br>
> Harris-Foulkes estimate = -183.20894409 Ry<br>
> estimated scf accuracy < 0.00000009 Ry<br>
><br>
> total all-electron energy = -551.326109 Ry<br>
><br>
> The total energy is the sum of the following terms:<br>
><br>
> one-electron contribution = -85.45759185 Ry<br>
> hartree contribution = 47.93001020 Ry<br>
> xc contribution = -25.22704816 Ry<br>
> ewald contribution = -95.91150866 Ry<br>
> one-center paw contrib. = -24.54280555 Ry<br>
><br>
><br>
> while adding the -pw2casino flag and running:<br>
><br>
> mpirun -np 4 pw.x -pw2casino <project.in> project.out<br>
><br>
><br>
> will print the original lines above, as well as additional lines like:<br>
><br>
> Energies determined by pw2casino tool<br>
> -------------------------------------<br>
> Kinetic energy 21.442363717309750 au =<br>
> 42.884727434619499 Ry<br>
> Local energy -84.609446621660297 au =<br>
> -169.21889324332059 Ry<br>
> Non-Local energy 20.438286977076853 au =<br>
> 40.876573954153706 Ry<br>
> Ewald energy -47.955754330605572 au =<br>
> -95.911508661211144 Ry<br>
> xc contribution -12.613524077939315 au =<br>
> -25.227048155878631 Ry<br>
> hartree energy 23.965005100370561 au =<br>
> 47.930010200741123 Ry<br>
> Total energy -79.333069235448022 au =<br>
> -158.66613847089604 Ry<br>
><br>
><br>
> From the above output lines, we can see that the total energies (-183.20894402<br>
> Ry and -158.66613847089604 Ry) are pretty terribly mismatched. I find<br>
> that this mismatch occurs when using either PWSCF v5.0.2 and PWSCF v6.0.<br>
> When I was using PWSCF v6.0, quantum-espresso actually prints an "Energy<br>
> mismatch error" (I'd paste the error below but my laptop, which runs v6.0,<br>
> is currently in the shop), so at least it recognizes that something is<br>
> wrong.<br>
><br>
> I've tried looking through some other threads and attempted to apply fixes<br>
> that worked for other people, such as adding the following line:<br>
><br>
> assume_isolated = 'mp'<br>
><br>
> to my &CONTROL in my input file (source:<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=" id="LPlnk948772" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=</a>),<br>
> but nothing so far has worked.<br>
><br>
> Does anyone have any idea what the issue is and how to resolve it?<br>
> Because with the demonstrated mismatch in total energies, I do not have<br>
> faith in the kinetic energy value that is output. If this helps, I am using<br>
> the following pseudopotentials:<br>
><br>
> file Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2P renormalized<br>
><br>
> file O.pbe-nl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2S 2P renormalized<br>
><br>
><br>
> and below is my input file for MgO in a CsCl-type crystalline arrangement<br>
> (which generated the energy values I gave above, using PWSCF v5.0.2):<br>
><br>
> &CONTROL<br>
> calculation = "scf"<br>
> max_seconds = 8.64000e+04<br>
> outdir = "./"<br>
> prefix = "MgO_CsCl"<br>
> pseudo_dir = "/home/davidho/.burai/.pseudopot"<br>
> title = "Working MgO CsCl Structure(SCF)"<br>
> wf_collect = .TRUE.<br>
> /<br>
><br>
> &SYSTEM<br>
> a = 2.65994e+00<br>
> degauss = 1.00000e-02<br>
> ecutrho = 2.00000e+03<br>
> ecutwfc = 2.00000e+02<br>
> ibrav = 1<br>
> nat = 2<br>
> ntyp = 2<br>
> occupations = "fixed"<br>
> smearing = "gaussian"<br>
> assume_isolated = 'mp'<br>
> /<br>
><br>
> &ELECTRONS<br>
> conv_thr = 1.00000e-06<br>
> electron_maxstep = 200<br>
> mixing_beta = 7.00000e-01<br>
> startingpot = "atomic"<br>
> startingwfc = "atomic+random"<br>
> /<br>
><br>
> K_POINTS {automatic}<br>
> 10 10 10 0 0 0<br>
><br>
> ATOMIC_SPECIES<br>
> Mg 24.30500 Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF<br>
> O 15.99940 O.pbe-nl-kjpaw_psl.1.0.0.UPF<br>
><br>
> ATOMIC_POSITIONS {angstrom}<br>
> Mg 0.000000 0.000000 0.000000<br>
> O 1.329970 1.329970 1.329970<br>
><br>
><br>
><br>
> Any help and further insight into both the problem with and/or the method<br>
> of energy decomposition would be greatly appreciated.<br>
><br>
> Thank you,<br>
> Wai-Ga David Ho, Department of Physics at Florida State University<br>
><br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> users mailing list<br>
> users@lists.quantum-espresso.org<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk413286" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
<br>
<br>
<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_37c357f7_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=MdQJWd_lLSRzCyfln4VPPrabXgu4ZOHRq34HfA1LNVg&e=" id="LPlnk641716" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_37c357f7_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=MdQJWd_lLSRzCyfln4VPPrabXgu4ZOHRq34HfA1LNVg&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 12<br>
Date: Thu, 21 Feb 2019 09:35:58 +0100<br>
From: Pietro Davide Delugas <pdelugas@sissa.it><br>
To: users@lists.quantum-espresso.org<br>
Subject: Re: [QE-users] non-convergent scf calculation on organometal<br>
perovskite structure<br>
Message-ID: <fb09ead9-e310-4142-0c01-7c8af7a8ab46@sissa.it><br>
Content-Type: text/plain; charset="utf-8"; Format="flowed"<br>
<br>
Hi<br>
<br>
Have you tried to increase the k_point mesh ?? 4 4 4 seems a little bit <br>
lax as mesh for MAPbI3.<br>
If I remember well I am afraid that to get convergence you will need <br>
something like 10X10X10.<br>
As for the structure neighboring methylammoniums? like to orient <br>
differently one from the other, you should probably use a larger cell.?? <br>
Also consider to add some correction for van der Waals interactions see <br>
here ( <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_Doc_INPUT-5FPW.html-23idm45922794348896&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cqtH8edx5niAN3Ied1OTnBUieUdWyMQAYEhXLAvFhsY&e=" id="LPlnk607261" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__www.quantum-2Despresso.org_Doc_INPUT-5FPW.html-23idm45922794348896&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cqtH8edx5niAN3Ied1OTnBUieUdWyMQAYEhXLAvFhsY&e=</a>)<br>
<br>
hope it helps<br>
Pietro<br>
<br>
<br>
On 02/21/2019 04:17 AM, Julien Barbaud wrote:<br>
> Dear users,<br>
><br>
><br>
> I am new to QE, and trying to run a simple scf calculation on a <br>
> CH3NH3PbI3 crystal (semi-conducting material). I am using ultrasoft <br>
> pseudopotentials based on the exchange-correlation functionnal PBEsol.<br>
><br>
> I set up a first input, with values of parameters inspired from <br>
> literature on the subject. However, I could not reach convergence <br>
> after 100 iterations. The estimated error was actually "exploding" to <br>
> very high values, indicating a serious problem. I tried several <br>
> changes but was unsuccessful:<br>
><br>
> * ?varying plane-wave cutoff energy does not solve the problem (cf<br>
> attached ecut.png, giving the estimated error as a function of the<br>
> number of iterations. It is shown here only on the first 15<br>
> iterations as the results pretty much only stall from there)<br>
> * ?varying cutoff energy for charge (cf ecutrho.png)<br>
> * taking larger k-point sampling (not shown)<br>
> * ?I also read that for metallic or "close to metallic conductors",<br>
> there might be problems with the first unoccupied states that can<br>
> be solved by adding a few empty bands. My system being a<br>
> semi-conductor, I tried adding additional bands using a m-p<br>
> smearing but no improvement was found (not shown)<br>
><br>
><br>
> The only change that I found effective was to reduce the mixing_beta <br>
> factor.<br>
><br>
><br>
> It effectively prevents the error from diverging to very large values, <br>
> but I still do not reach convergence, even after longer iterations. I <br>
> tried much smaller values of mixing beta which improves the final <br>
> value of the error, but I still cannot reach convergence on 100 <br>
> iterations. As shown in the mixbeta2_zoom.png, the error reduces to <br>
> smaller values around ~1e-5~1e-6, but it keeps stalling after a while. <br>
> I do not observe a well-converging behaviour for any value.<br>
><br>
><br>
> I attached the "default version" of my script on which the various <br>
> modifications described above have been independently performed. I <br>
> obtained the geometry from a CIF file in literature and checked it <br>
> with visualization software; it seems perfectly ok as far as I can tell.<br>
><br>
><br>
> Any insight on what I did wrong would be really helpful. I suspect a <br>
> shameful beginner mistake, but can not find it out.<br>
><br>
><br>
> Thanks in advance,<br>
><br>
> Julien barbaud<br>
><br>
><br>
> P.S: this is my first time posting on this user list. Please let me <br>
> know if my question is not suitable for it, or can be improved either <br>
> in its content or presentation. I will gladly take any recommandation <br>
> into account in order not to negatively impact the quality of this <br>
> user list !<br>
><br>
><br>
><br>
> _______________________________________________<br>
> users mailing list<br>
> users@lists.quantum-espresso.org<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk235557" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_020e5d38_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=_X24YK8ccrs3nQC8xcZMnShKAToNIzf11poxRNrWl5o&e=" id="LPlnk362875" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_020e5d38_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=_X24YK8ccrs3nQC8xcZMnShKAToNIzf11poxRNrWl5o&e=</a>><br>
<br>
------------------------------<br>
<br>
Message: 13<br>
Date: Thu, 21 Feb 2019 09:37:20 +0100<br>
From: Paolo Giannozzi <p.giannozzi@gmail.com><br>
To: Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
Subject: Re: [QE-users] pw2casino mismatch energy error<br>
Message-ID:<br>
<CAPMgbCtOF+8sVaGUaDXkP-dWxfJSoQPpdHD+LgvDiugF6scdow@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
On Thu, Feb 21, 2019 at 9:07 AM Paolo Giannozzi <p.giannozzi@gmail.com><br>
wrote:<br>
<br>
If I remember correctly, it is not supposed to work for Ultrasoft<br>
> pseudopotentials and PAW<br>
><br>
<br>
I didn't remember correctly: it produces correct energies for USPP, but not<br>
for PAW<br>
<br>
Paoo<br>
<br>
<br>
<br>
> Paolo<br>
> On Thu, Feb 21, 2019 at 3:02 AM Wai-Ga Ho <wdh17@my.fsu.edu> wrote:<br>
><br>
>> Hello QE community,<br>
>><br>
>><br>
>> I am currently attempting to extract the kinetic energy contributions to<br>
>> the total energies of different polymorphs of MgO.<br>
>><br>
>><br>
>> I understand that we need to use the -pw2casino flag if we want our PWSCF<br>
>> output to print a more detailed breakdown of the total energy in terms of<br>
>> its contributions. This means that using the command:<br>
>><br>
>> mpirun -np 4 pw.x <project.in> project.out<br>
>><br>
>><br>
>> will produce a project.out file which contains some text like this:<br>
>><br>
>><br>
>> ! total energy = -183.20894402 Ry<br>
>> Harris-Foulkes estimate = -183.20894409 Ry<br>
>> estimated scf accuracy < 0.00000009 Ry<br>
>><br>
>> total all-electron energy = -551.326109 Ry<br>
>><br>
>> The total energy is the sum of the following terms:<br>
>><br>
>> one-electron contribution = -85.45759185 Ry<br>
>> hartree contribution = 47.93001020 Ry<br>
>> xc contribution = -25.22704816 Ry<br>
>> ewald contribution = -95.91150866 Ry<br>
>> one-center paw contrib. = -24.54280555 Ry<br>
>><br>
>><br>
>> while adding the -pw2casino flag and running:<br>
>><br>
>> mpirun -np 4 pw.x -pw2casino <project.in> project.out<br>
>><br>
>><br>
>> will print the original lines above, as well as additional lines like:<br>
>><br>
>> Energies determined by pw2casino tool<br>
>> -------------------------------------<br>
>> Kinetic energy 21.442363717309750 au =<br>
>> 42.884727434619499 Ry<br>
>> Local energy -84.609446621660297 au =<br>
>> -169.21889324332059 Ry<br>
>> Non-Local energy 20.438286977076853 au =<br>
>> 40.876573954153706 Ry<br>
>> Ewald energy -47.955754330605572 au =<br>
>> -95.911508661211144 Ry<br>
>> xc contribution -12.613524077939315 au =<br>
>> -25.227048155878631 Ry<br>
>> hartree energy 23.965005100370561 au =<br>
>> 47.930010200741123 Ry<br>
>> Total energy -79.333069235448022 au =<br>
>> -158.66613847089604 Ry<br>
>><br>
>><br>
>> From the above output lines, we can see that the total energies (-183.20894402<br>
>> Ry and -158.66613847089604 Ry) are pretty terribly mismatched. I find<br>
>> that this mismatch occurs when using either PWSCF v5.0.2 and PWSCF v6.0.<br>
>> When I was using PWSCF v6.0, quantum-espresso actually prints an "Energy<br>
>> mismatch error" (I'd paste the error below but my laptop, which runs v6.0,<br>
>> is currently in the shop), so at least it recognizes that something is<br>
>> wrong.<br>
>><br>
>> I've tried looking through some other threads and attempted to apply<br>
>> fixes that worked for other people, such as adding the following line:<br>
>><br>
>> assume_isolated = 'mp'<br>
>><br>
>> to my &CONTROL in my input file (source:<br>
>> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=" id="LPlnk603059" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_2015-2DApril_032059.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=cHund3hSI10TIwAXR5wAmUGD9isDhgrOKHkN-QvOjXc&e=</a>),<br>
>> but nothing so far has worked.<br>
>><br>
>> Does anyone have any idea what the issue is and how to resolve it?<br>
>> Because with the demonstrated mismatch in total energies, I do not have<br>
>> faith in the kinetic energy value that is output. If this helps, I am using<br>
>> the following pseudopotentials:<br>
>><br>
>> file Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2P renormalized<br>
>><br>
>> file O.pbe-nl-kjpaw_psl.1.0.0.UPF: wavefunction(s) 2S 2P renormalized<br>
>><br>
>><br>
>> and below is my input file for MgO in a CsCl-type crystalline arrangement<br>
>> (which generated the energy values I gave above, using PWSCF v5.0.2):<br>
>><br>
>> &CONTROL<br>
>> calculation = "scf"<br>
>> max_seconds = 8.64000e+04<br>
>> outdir = "./"<br>
>> prefix = "MgO_CsCl"<br>
>> pseudo_dir = "/home/davidho/.burai/.pseudopot"<br>
>> title = "Working MgO CsCl Structure(SCF)"<br>
>> wf_collect = .TRUE.<br>
>> /<br>
>><br>
>> &SYSTEM<br>
>> a = 2.65994e+00<br>
>> degauss = 1.00000e-02<br>
>> ecutrho = 2.00000e+03<br>
>> ecutwfc = 2.00000e+02<br>
>> ibrav = 1<br>
>> nat = 2<br>
>> ntyp = 2<br>
>> occupations = "fixed"<br>
>> smearing = "gaussian"<br>
>> assume_isolated = 'mp'<br>
>> /<br>
>><br>
>> &ELECTRONS<br>
>> conv_thr = 1.00000e-06<br>
>> electron_maxstep = 200<br>
>> mixing_beta = 7.00000e-01<br>
>> startingpot = "atomic"<br>
>> startingwfc = "atomic+random"<br>
>> /<br>
>><br>
>> K_POINTS {automatic}<br>
>> 10 10 10 0 0 0<br>
>><br>
>> ATOMIC_SPECIES<br>
>> Mg 24.30500 Mg.pbe-spnl-kjpaw_psl.1.0.0.UPF<br>
>> O 15.99940 O.pbe-nl-kjpaw_psl.1.0.0.UPF<br>
>><br>
>> ATOMIC_POSITIONS {angstrom}<br>
>> Mg 0.000000 0.000000 0.000000<br>
>> O 1.329970 1.329970 1.329970<br>
>><br>
>><br>
>><br>
>> Any help and further insight into both the problem with and/or the method<br>
>> of energy decomposition would be greatly appreciated.<br>
>><br>
>> Thank you,<br>
>> Wai-Ga David Ho, Department of Physics at Florida State University<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> users mailing list<br>
>> users@lists.quantum-espresso.org<br>
>> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk171317" class="OWAAutoLink" previewremoved="true">
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
><br>
><br>
><br>
> --<br>
> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
> Phone +39-0432-558216, fax +39-0432-558222<br>
><br>
><br>
<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_b6726211_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=R5kD1BTJWvmgEkiu5asQ63O_CFOShmmzYZvg2m3wyKc&e=" id="LPlnk642801" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.quantum-2Despresso.org_pipermail_users_attachments_20190221_b6726211_attachment-2D0001.html&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=R5kD1BTJWvmgEkiu5asQ63O_CFOShmmzYZvg2m3wyKc&e=</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
users mailing list<br>
users@lists.quantum-espresso.org<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=" id="LPlnk894396" class="OWAAutoLink" previewremoved="true">https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.quantum-2Despresso.org_mailman_listinfo_users&d=DwICAg&c=MNHwOqQ8N1u91SoMLfIblwuGXKgp50OPUXjl8uRAbak&r=B7toNJkJN9wC1at1deS0cQ&m=Szr9aWxo68iMhkcBDOeDefYD6tBwed0WLV_dzlRu8J0&s=fmj1OWGBtl2pVfQfpBDlrj33gaZ2M4inEMxwkZDImso&e=</a><br>
<br>
------------------------------<br>
<br>
End of users Digest, Vol 139, Issue 22<br>
**************************************<br>
<br>
</div>
</span></font></div>
</div>
</div>
</body>
</html>