[Wien] Problems in phonon calculations (4)
Ghosh SUDDHASATTWA
ssghosh at igcar.gov.in
Mon Jan 10 04:33:00 CET 2011
Dear Prof. Blaha,
Thank you very much for the detailed reply
I had taken a 48 atom supercell with Pnma as space group in the phonon
calculations
I had used a about 1400 k-points in the IBZ with rkmax 7.00
However, I had used rkmax of 8.00 when I optimized the lattice parameters.
After it took more than 40 hours to complete lapw1, I used a 12 atom
supercell and completed case_1. I need some help here. I copy here the
case_1.scf
TOTAL FORCE IN mRy/a.u. = |F| Fx Fy Fz
with/without FOR in case.in2
:FOR001: 1.ATOM 7.217 -7.180 0.732
0.000 pa
rtial forces
:FOR002: 2.ATOM 11.249 -7.074 -8.746
0.000 pa
rtial forces
:FOR003: 3.ATOM 7.417 7.416 0.109
0.000 pa
rtial forces
:FOR004: 4.ATOM 6.958 6.957 0.120
0.000 pa
rtial forces
:FOR005: 5.ATOM 14.967 14.967 0.080
0.000 pa
rtial forces
:FOR006: 6.ATOM 14.973 14.959 0.654
0.000 pa
rtial forces
:FOR007: 7.ATOM 14.113 -14.028 1.543
0.000 pa
rtial forces
:FOR008: 8.ATOM 15.810 -15.749 1.386
0.000 pa
rtial forces
:FOR009: 9.ATOM 83.030 82.549 8.927
0.000 pa
rtial forces
:FOR010: 10.ATOM 95.213 94.803 8.826
0.000 pa
rtial forces
:FOR011: 11.ATOM 88.802 -88.782 1.906
0.000 pa
rtial forces
:FOR012: 12.ATOM 88.672 -88.672 -0.236
0.000 pa
rtial forces
TOTAL FORCE WITH RESPECT TO THE GLOBAL COORDINATE SYSTEM:
:FGL001: 1.ATOM 0.732000000 0.000000000
-7.180000000 pa
rtial forces
:FGL002: 2.ATOM -8.746000000 0.000000000
-7.074000000 pa
rtial forces
:FGL003: 3.ATOM 0.109000000 0.000000000
7.416000000 pa
rtial forces
:FGL004: 4.ATOM 0.120000000 0.000000000
6.957000000 pa
rtial forces
:FGL005: 5.ATOM 0.080000000 0.000000000
14.967000000 pa
rtial forces
:FGL006: 6.ATOM 0.654000000 0.000000000
14.959000000 pa
rtial forces
:FGL007: 7.ATOM 1.543000000 0.000000000
-14.028000000 pa
rtial forces
:FGL008: 8.ATOM 1.386000000 0.000000000
-15.749000000 pa
rtial forces
:FGL009: 9.ATOM 8.927000000 0.000000000
82.549000000 pa
rtial forces
:FGL010: 10.ATOM 8.826000000 0.000000000
94.803000000 pa
rtial forces
:FGL011: 11.ATOM 1.906000000 0.000000000
-88.782000000 pa
rtial forces
:FGL012: 12.ATOM -0.236000000 0.000000000
-88.672000000 pa
rtial forces
May I kindly request you to suggest on the forces. Is it ok or do I have to
change -fc 0.1 to -fc 0.01?
Thank you
Suddhasattwa
_____
From: wien-bounces at zeus.theochem.tuwien.ac.at
[mailto:wien-bounces at zeus.theochem.tuwien.ac.at] On Behalf Of Peter Blaha
Sent: Saturday, January 08, 2011 7:18 PM
To: A Mailing list for WIEN2k users
Subject: Re: [Wien] Problems in phonon calculations (3)
You need to understand what you are doing:
I guess you submit your job in a directory AB2_phonon.
The job creates a .machines file in THIS directory (look at it with some
editor and verify that it looks as expected (with 24 lines)).
Now understand the run_phonon job file.
Inside the loop over all 18 cases it has a:
cd case_$i and inside this directory it makes a run_lapw ...
So change into case_1 and examine the .machines file.
It will of course NOT show the 24 lines, but is some dummy-default .machines
file with just 2 lines.
so before the cd case_$i insert a line:
cp .machines case_$i
into the run_phonon job.
PS: If your job is still running, you can also do this manually on the
commandline (for all case_1 ... case_18)
and the proper .machines file will be used for all new iterations.
PPS: If I look into your cpu-time (more than 40 h !! for lapw1) I doubt that
you want to wait
until this job has finished (even on 24 cpus). How large is your supercell -
this determines how you can
reduce the k-mesh!
PPS: You should learn on small examples how to make meaningful calculations.
You must examine for
YOUR system (nobody can help you here) and your properties, what a good
RKMAX or k-mesh is.
In other words: caclculate phonons with a "low RKMAX and few k-points", and
increase the parameters
until the frequencies did not change significantly anymore.
Am 08.01.2011 07:42, schrieb Ghosh SUDDHASATTWA:
> Dear Wien2k users, and Prof. Blaha,
>
> Can you please help me out in the following querie. As suggested earlier,
I had initialized my phonon calculations with 48 atoms per lattice.
>
> I created and initialized case_1
>
> My phonon job file is
>
> #!/bin/csh -f
>
> #
>
> set file=AB2_phonon
>
> #
>
> foreach i ( \
>
> 1 \
> 17 \
>
> 18 \
>
> )
>
> cd case_$i
>
> echo running case_$i
>
> #
>
> # select other options if necessary
>
> run_lapw -I -i 200 -cc 0.0001 -in1ef -p -fc 0.1
>
> #
>
> # select other save-name if necessary
>
> save_lapw case_${i}_gga_rkm7.00_3000k
>
> cd ..
>
> end
>
> I submitted the job by qsub -pe kpoint 24 kpoint.sh
>
> My machine file is
>
> 1:ibnx70
>
> 1:ibnx70
>
> 1:ibnx70
>
> 1:ibnx81
>
> 1:ibnx81
>
> 1:ibnx81
>
> 1:ibnx81
>
> 1:ibnx81
>
> 1:ibnx79
>
> 1:ibnx79
>
> 1:ibnx79
>
> 1:ibnx65
>
> 1:ibnx65
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx60
>
> 1:ibnx96
>
> 1:ibnx96
>
> 1:ibnx96
>
> granularity:1
>
> extrafine:1
>
> For the last three days, the queuing shows the job is running
>
> Surprisingly, there is no "running" mode in w2web.
>
> I checked the case_1 directory
>
> case_1.dayfile shows (surprisingly)
>
> Calculating case_1 in
/group5/cg/sghosh/WIEN2k/lapw/phonon/UZr2_phonon/case_1
>
> on nx70.igcar.gov.in with PID 10356
>
> start (Thu Jan 6 19:11:06 IST 2011) with lapw0 (-in1ef/99 to go)
>
> cycle 1 (Thu Jan 6 19:11:06 IST 2011) (-in1ef/99 to go)
>
>> lapw0 -p (19:11:06) starting parallel lapw0 at Thu Jan 6 19:11:06 IST 20
>
> 11
>
> -------- .machine0 : processors
>
> running lapw0 in single mode
>
> 177.410u 5.691s 3:11.36 95.6% 0+0k 0+0io 17pf+0w
>
>> lapw1 -c -p (19:14:17) starting parallel lapw1 at Thu Jan 6 19:14:1
>
> 8 IST 2011
>
> -> starting parallel LAPW1 jobs at Thu Jan 6 19:14:19 IST 2011
>
> running LAPW1 in parallel mode (using .machines)
>
> 2 number_of_parallel_jobs
>
> localhost(672) 90833.210u 359.642s 1+01:28:03.65 99.46% 0+0k 0+0io 0pf
>
> +0w
>
> 105493.144u 437.929s 1+13:49:58.03 77.78% 0+0k 0+0io 0pf+0w
>
> 120959.238u 433.725s 1+09:50:20.14 99.65% 0+0k 0+0io 0pf+0w
>
> localhost(672) 90579.555u 334.655s 1+01:22:25.95 99.53% 0+0k 0+0io 0pf
>
> +0w
>
> 105170.999u 500.369s 1+13:44:24.32 77.78% 0+0k 0+0io 0pf+0w
>
> 120363.259u 426.543s 1+09:40:41.44 99.63% 0+0k 0+0io 0pf+0w
>
> Summary of lapw1para:
>
> localhost k=1344 user=181413 wallclock=170
>
> 73.491u 124.896s 33:50:25.33 0.1% 0+0k 0+0io 0pf+0w
>
>> lapw2 -c -p (05:04:43) running LAPW2 in parallel mode
>
> [1] 4250
>
> [2] 4292
>
> That is 2 number of parallel jobs. I thought it should be 24
>
> Can you please why is this happening?
>
> Thank you
>
> Suddhasattwa
>
> SUDDHASATTWAGHOSH
>
> Scientific Officer (D)
>
> Pyrochemical Process Studies Section
>
> Fuel Chemistry Division
>
> Chemistry Group
>
> IndiraGandhi Centre for Atomic Research
>
> Kalpakkam
>
> Tamilnadu
>
> 603102
>
> India
>
> Phone: 91-44-27480500 (Ext:24283)
>
> Fax: 91-44-27480065
>
>
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
--
-----------------------------------------
Peter Blaha
Inst. Materials Chemistry, TU Vienna
Getreidemarkt 9, A-1060 Vienna, Austria
Tel: +43-1-5880115671
Fax: +43-1-5880115698
email: pblaha at theochem.tuwien.ac.at
-----------------------------------------
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
_____
No virus found in this message.
Checked by AVG - www.avg.com
Version: 10.0.1170 / Virus Database: 1435/3366 - Release Date: 01/07/11
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20110110/9a15fead/attachment.htm>
More information about the Wien
mailing list