[QE-users] [SUSPECT ATTACHMENT REMOVED] Re: Re: Re: pw_readschema_file failed retrieving input info from xml file ph.x using grid in qe-6.3-bkp and 2D cutoff
RAPHAEL LONGUINHOS MONTEIRO LOBATO
raphael.lobato at ufla.br
Tue Mar 12 19:07:27 CET 2019
Hi,
The test-phgrid.zip contains the example related to the xml missing
information (also the celldm).
Folders and content:
i) check-full-qe-5.1.1 contains the single-q without grid.
ii) qe-5.1.1-grid contains the single-q with grid.
iii) check-full contains the single-q without grid using 3D and 2D codes
.. the corresponding outputs have 3D and 2D appended in their filename.
iv) cpBands, linkBands, and noXML contain tests with grid in the qe-6.3-bkp
using 2D code trying different "small changes" on the grid script. I have
some issue with the huge amount of data in scratch/_ph0 when your system
has, say 270 irrs, ... so I symlink some files which are created and use
just for reading in the dfpt scf steps of the irrs. I thought the xml issue
was due that, but it appears it is not. Regardless these small changes, all
them lead to the same identical wrong frequencies in comparison to without
grid 6.3-bkp 2D code in iii).
In the below, I set q=1 and j loops over the irrs, and I copied the
test.phsave (and did others grid steps)
In the cpBands, I cp the scratch/_ph0/test.save to
scratch/1.X/_ph0/test.save, where X is the irr.
for data in `ls scratch/_ph0 | grep -e $prefix.xml -e $prefix.wfc -e
$prefix.save | xargs `; do
cp -r scratch/_ph0/$data
scratch/$q.$j/_ph0/.
done
In the linkBands, I just link it (ln -s ../../_ph0/$data
scratch/$q.$j/_ph0/)
for data in `ls scratch/_ph0 | grep -e $prefix.xml -e $prefix.wfc -e
$prefix.save | xargs `; do
ln -s ../../_ph0/$data
scratch/$q.$j/_ph0/.
done
In the noXML, I forgot the -e $prefix.xml in the script for linkBands and
did not copied the xml (it is inside the scratch/_ph0/test.save/ anyway).
v) reduce-io-false is a double-check of the tests in iv) ... with grid in
qe-6.3-bkp with 2D.
vi) reduce-io-false-3D just unset the 2D cutoff. The results are the same
of the folders i) ii) and iii) with 3D name, as expected.
vii) reduce-io-true is the 2D test of the 6.3-bkp where I set reduce_io
true. The results are like the ones in the folder iii) with 2D name
(6.3-bkp without grid)
About celldm, you can see for example the scfout and phout in the folder
iii) check-full. e.g.
scfout2D
celldm(1)= 1.889726 celldm(2)= 0.000000 celldm(3)= 0.000000
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
phout2D
celldm(1)= 4.620461 celldm(2)= 0.000000 celldm(3)=
7.361835
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
or
scfout3D
celldm(1)= 1.889726 celldm(2)= 0.000000 celldm(3)= 0.000000
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
phout3D
celldm(1)= 4.620461 celldm(2)= 0.000000 celldm(3)= 7.361835
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
and the dyn header
1 2 0 4.6204611 0.0000000 7.3618345 0.0000000 0.0000000
0.0000000
by opening the scf and axsf in xcrysden you can see it.
The examples of disp case are in the disp-grid-no-grid-6.3-bkp.tar.gz ...
the comparison of the dynX obtained without grid (in file k661-grid) and
with grid (k661-check-full) shows same frequencies.
best,
I'm very sorry Professor Giannozzi ... I compiled the 6.3-bkp by these
days (after Carnival, actually ... ) ... previously (Carnival) I kept
hand-correcting the dyn header before post-processing it while using
previous versions I reported the celldm thing, six months ago.
Em ter, 12 de mar de 2019 às 11:21, Paolo Giannozzi <p.giannozzi at gmail.com>
escreveu:
>
> Please provide an example backing your claim (NOW, not in 6 months)
>
> On Tue, Mar 12, 2019 at 2:26 PM RAPHAEL LONGUINHOS MONTEIRO LOBATO <
> raphael.lobato at ufla.br> wrote:
>
>> Hi,
>> Thanks for be working to fix it.
>>
>> The other subject of the message ... "wrong celldm values" ... was
>> previously replied at
>>
>> http://lists.quantum-espresso.org/pipermail/users/2018-October/041419.html
>>
>> which thought it was solved with the modification reported in:
>>
>> https://gitlab.com/QEF/q-e/commit/9a579539b65bd1502d1f28747e483a1d1059f5f8
>>
>> In the pw_restart_new.f90 of qe-6.3-bkp I see as reported above
>>
>> 1199 CALL at2celldm (ibrav,1.0_dp,at(:,1),at(:,2),at(:,3),celldm) !
>> chance to 1.0_dp in the place of the removed alat
>>
>> Yet, the ph.x reads wrong celldm values, prints them wrongly in the dyn
>> header, and the issue of rescaling some quantities from dynmat.x remains.
>>
>> best,
>>
>> Sorry ... I forgot to mention in previous post I'm from:
>> Associate Professor Raphael Longuinhos Monteiro Lobato
>> CV:http://lattes.cnpq.br/1475483971471395.
>> Research Group site:http://jenainassoares2.wixsite.com/nanomat
>> (Portuguese only - English version on the way)
>> Departamento de Fisica, Universidade Federal de Lavras
>> Av. Dr. Silvio Menicucci 1001, bairro Kenedy
>> Post office box 3037
>> 37200-000 - Lavras, MG, Brazil
>>
>> ------------------------------
>>>
>>> Message: 8
>>> Date: Tue, 12 Mar 2019 11:04:58 +0100
>>> From: Pietro Davide Delugas <pdelugas at sissa.it>
>>> To: users at lists.quantum-espresso.org
>>> Subject: Re: [QE-users] pw_readschema_file failed retrieving input
>>> info from xml file ph.x using grid in qe-6.3-bkp and 2D cutoff
>>> Message-ID: <a89d561d-1029-5416-7229-da026757a704 at sissa.it>
>>> Content-Type: text/plain; charset=utf-8; format=flowed
>>>
>>> Hi
>>>
>>> Thanks for reporting the issue.
>>>
>>> As for the message about the missing input element, no problem it is the
>>> desired behavior, run_nscf started bby ph would not be able to fill
>>> correctly the
>>> input element because the variable in input_parameters are not
>>> initialized, so the element is? not printed at all.
>>>
>>> The isseus comes from the fact that in writing the status of the 2D
>>> cutoff we are currenlty using the value stored in input_parameters
>>> which, not being initialized, has in this case a wrong .false. value.
>>>
>>> I will fix the issue in the develop, for further information check this
>>> issue on gitlab https://gitlab.com/QEF/q-e/issues/102
>>>
>>> Pietro
>>>
>>>
>>>
>>> On 03/11/2019 09:11 PM, Raphael Longuinhos Monteiro Lobato wrote:
>>> > Hi,
>>> >
>>> > Using ph.x with grid may lead to
>>> >
>>> > ????? Reading data from directory:
>>> > ????? scratch/X.Y/_ph0/test.save/
>>> > ????? Message from routine pw_readschema_file:
>>> > ????? failed retrieving input info from xml file, check it !!!
>>> >
>>> > and wrong results.
>>> >
>>> > In more details (sorry if not too methodical ... but I tried ... spent
>>> > all day on it):
>>> >
>>> > There is no problem with Gamma calculations ... as there is no bands
>>> > to be read from /scratch/_ph0/ ... I guess.
>>> >
>>> > Single-q, out Gamma (say X Y 0 instead 0. 0. 0.): tests without grid
>>> > (without breaking in start_irr) in qe-6.3-bkp (3D), and in qe-5.1.1
>>> > without/with grid leads to same results (as expected).
>>> >
>>> > Using grid in qe-6.3-bkp (2D):
>>> >
>>> > ? + using only_init true leads to the message of error above ... the
>>> > calculations goes on as if all were all right, and the frequencies are
>>> > all wrong at the end (after collect).
>>> >
>>> > ? + using only_init true and reduce_io true (thus repeating the band
>>> > structure calculations done at only_init calculations each time ...
>>> > the xml is not saved in /scratch/_ph0/) leads to frequencies as
>>> > expected when comparison to without grid calculations.
>>> >
>>> > Using grid in qe-6.3-bkp (3D):
>>> > ? + using only_init true leads to the expected results (comparison
>>> > with previous 3D cases).
>>> >
>>> > For me, it indicates it may be a problem when reading the information
>>> > of 2D cutoff in the in /scrath/_ph0/test.save/data-file-schema.xml
>>> > (same as /scrath/_ph0/test.xml).
>>> >
>>> > The comparison between the scrathc/test.xml and /scratch/_ph0/test.xml
>>> > when using the only_init true and reduce_io false shows indeed some
>>> > unexpected discrepancies (e.g., from what I could get, it misses the
>>> > fact that the scf was performed with assume isolated flag)
>>> >
>>> > (note ... I use wf_collect false in the scf calculation and norm
>>> > conserving pseudopotentials)
>>> >
>>> > The above were single-q calculations. For the case of dispersion
>>> > calculations using 2D, reading from /scratch/X.Y/test.save/ does not
>>> > display the message. However, when reading the bands (reading from
>>> > /scratch/X.Y/_ph0/test.q_X/test.save/data-file-schema.xml) it does
>>> > display it (except for the 0. 0. 0. point ... in this case the bands
>>> > are read from /scratch/test.save/ ) ... surprisingly, the calculations
>>> > are correct ... I guess because in this case the code, somehow, does
>>> > not miss the 2D flag, as it informs on the output of out.X.Y it is
>>> > using 2D... the question remains "what is the information it misses
>>> > when reading the bands from xml ... and when it will lead to error on
>>> > the results ? "
>>> >
>>> > Also, a previous post
>>> >
>>> http://lists.quantum-espresso.org/pipermail/users/2018-October/041419.html
>>> > ... the celldm problem in the dyn header remains in 6.3-backports ...
>>> > and one have to correct it by hand on the .dyn
>>> >
>>> > If these are indeed unexpected error (not my ...? if someone can
>>> > reproduce it) yet not solved (sorry if I miss the solution in older
>>> > posts), the real problem is ... I don't have idea to help solving them
>>> > ...? although I could give it a try with some guide.
>>> >
>>> > best,
>>> >
>>>
>> _______________________________________________
>> users mailing list
>> users at lists.quantum-espresso.org
>> https://lists.quantum-espresso.org/mailman/listinfo/users
>
>
>
> --
> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
> Phone +39-0432-558216, fax +39-0432-558222
>
>
>
> --
> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
> Phone +39-0432-558216, fax +39-0432-558222
>
>
--
Associate Professor Raphael Longuinhos Monteiro Lobato
CV: http://lattes.cnpq.br/1475483971471395.
Research Group site: http://jenainassoares2.wixsite.com/nanomat (Portuguese
only - English version on the way)
Departamento de Fisica, Universidade Federal de Lavras
Av. Dr. Silvio Menicucci 1001, bairro Kenedy
Post office box 3037
37200-000 - Lavras, MG, Brazil
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20190312/a5db325c/attachment.html>
More information about the users
mailing list