[Pw_forum] wfc files: heavy I/O, handling for restarts

2011-09-06 Thread S. K. S.
Dear Prof. Paolo,

Thanks a lot for your reply.

> which QE code and which files are you referring to? <

That is phonon code (ph.x) and the files I mentioned before are given below
in detail.

With the earlier versions up to 4.1.3, phonon code runs fine and  I got
following files
(bold and underlined items are folders) in the tmp directory of a local
disk.
The nodes which are used by the codes are node186, 036, 139.

   node186:/tmpscratch/sksct8/tmp$ ls

> *pbmno.save* _phpbmno.com1_phpbmno.dwf2   _phpbmno.igk3
> _phpbmno.prd3
> pbmno.wfc1 _phpbmno.com2_phpbmno.dwf3   _phpbmno.mixd1
> _phpbmno.recover
> pbmno.wfc2 _phpbmno.com3_phpbmno.ebar1  _phpbmno.mixd2
> _phpbmno.recover2
> pbmno.wfc3 _phpbmno.dvkb31  _phpbmno.ebar2  _phpbmno.mixd3
> _phpbmno.recover3
> _phpbmno.bar1  _phpbmno.dvkb32  _phpbmno.ebar3  *_phpbmno.phsave** * *
> _phpbmno.save*
> _phpbmno.bar2  _phpbmno.dvkb33  _phpbmno.igk_phpbmno.prd1
> _phpbmno.bar3  _phpbmno.dwf1_phpbmno.igk2   _phpbmno.prd2
>

node036:/tmpscratch/sksct8/tmp$ ls
pbmno.wfc4 _phpbmno.dwf4  _phpbmno.mixd4  _phpbmno.recover4
pbmno.wfc5 _phpbmno.dwf5  _phpbmno.mixd5  _phpbmno.recover5
_phpbmno.bar4  _phpbmno.igk4  _phpbmno.prd4   _phpbmno.wfc4
_phpbmno.bar5  _phpbmno.igk5  _phpbmno.prd5   _phpbmno.wfc5

node139:/tmpscratch/sksct8/tmp$ ls
pbmno.wfc6 _phpbmno.bar8  _phpbmno.igk7  _phpbmno.recover6
_phpbmno.wfc8
pbmno.wfc7 _phpbmno.dwf6  _phpbmno.igk8  _phpbmno.recover7
pbmno.wfc8 _phpbmno.dwf7  _phpbmno.prd6  _phpbmno.recover8
_phpbmno.bar6  _phpbmno.dwf8  _phpbmno.prd7  _phpbmno.wfc6
_phpbmno.bar7  _phpbmno.igk6  _phpbmno.prd8  _phpbmno.wfc7

However, with the new version 4.3.1, for the exactly same input files and
job scripts I only get these files and nothing else :

node045:/tmpscratch/sksct84/tmp$ ls
*pbmno.save*  pbmno.wfc1  *_ph0*

node111:/tmpscratch/sksct84/tmp$ ls
pbmno.wfc2  pbmno.wfc3

node092:/tmpscratch/sksct84/tmp$ ls
pbmno.wfc4  pbmno.wfc5

node080:/tmpscratch/sksct84/tmp$ ls
pbmno.wfc6  pbmno.wfc7

node072:/tmpscratch/sksct84/tmp$ ls
pbmno.wfc8

Note that,  the used nodes in this time are  node045, node111, node092,
node080, node072.

So it is clear from the above example that somehow in the new version
4.3.1,  _phpbmno.save and
_phpbmno.phsave goes inside the directory "_ph0" ;

And the same phonon calculation,  which was running fine with the earlier
version, now stops like this way (a bit
abruptly and rudely, with out much informations or error messages) :


Electric field:
 Dielectric constant
 Born effective charges in two ways


 Atomic displacements:
 There are   5 irreducible representations

 Representation 1  3 modes -T_1u G_15  G_4- To be done

 Representation 2  3 modes -T_1u G_15  G_4- To be done

 Representation 3  3 modes -T_1u G_15  G_4- To be done

 Representation 4  3 modes -T_1u G_15  G_4- To be done

 Representation 5  3 modes -T_2u G_25  G_5- To be done



simply with this error :
--
MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD
with errorcode 0.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--
--
mpirun has exited due to process rank 2 with PID 8791 on
node node111.cvos.cluster exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--
[node045:11580] 6 more processes have sent help message help-mpi-api.txt /
mpi-abort
[node045:11580] Set MCA parameter "orte_base_help_aggregate" to 0 to see all
help / error messages


Hope this email explains much better.

Thanks and Regards,
Saha SK
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110906/78efc62c/attachment.htm
 


[Pw_forum] wfc files: heavy I/O, handling for restarts

2011-09-05 Thread S. K. S.
 # Scatter wfc restart files
   awk '{ files_for[$1] = files_for[$1] " '$basename'.wfc" NR }
   END { for (host in files_for) print host, files_for[host] }'
$PBS_NODEFILE \
   | while read host files
   do
   ssh -n $host "cd $PBS_O_WORKDIR; mv $files $ESPRESSO_TMPDIR/"
   done
   # on master host, copy .save directory as well
   rsync -a $basename.save $ESPRESSO_TMPDIR


   mpirun  -x ESPRESSO_TMPDIR \
   -np $(wc -l < $PBS_NODEFILE) \
   -machinefile  $PBS_NODEFILE  \
   pw.x -inp input.txt > output.txt

   # Gather remote files
   uniq $PBS_NODEFILE \
   | while read host
   do
   ssh -n $host "rsync -a $TMPDIR/ $PBS_O_WORKDIR/"
   done
   --

E.g. for a job with nodes=3:ppn=4 the scatter part would distribute the
existing files pwscf.wfc{1..12} as follows: <<

Yes, this looks indeed cumbersome. It becomes more painful when one can not
know a priori in which nodes his/her job
will go, particularly  when it is totally decided by the automatic queue
decider,  depending on the free nodes available.
In such a situation,  one's restarted job may go to a totally new set of
nodes, and phonon calculation can not get necessary
files to restart. Then restarting phonon calculation becomes more difficult.


It seems, there is a more serious trouble in the recent version of QE. In
the version before QE4.2, the QE codes used to
replicate the same necessary files to the distributed local  disks of all
the nodes. In this case, at least phonon calculation
can run smoothly instead of  crashing. But in the recent version, phonon
calculations just stop by complaining that
the distributed .wfc files in one node are not visible by another node.

If a quick remedy of this problem is not easy, then at least, for the time
being, it is better to keep the earlier option of replicating
the same .wfc files in all the nodes still working in the version 4.3.1.
Other better option  can be to implement the "WF_COLLECT"
trick, also in phonon code, as it is already there for PW.x.

Thanks and regards,
Saha SK
R&D Assistant
JNCASR
Bangalore 560064
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110905/ea2b77e7/attachment-0001.htm
 


[Pw_forum] LDA PZ USPP for Mg

2011-05-18 Thread S. K. S.
Dear All,

Would someone like to share a PZ-LDA USPP for Mg ?

regards,
Saha SK

stewart at cnf.cornell.eduwrote:

> Hi all,
>
> I was wondering if anyone had a reliable ultrasoft
> pseudopotential for Mg in
> the local density approximation.  The periodic table
> on the PWscf site has
> only norm-conserving and a GGA ultrasoft.
>
> Thanks,

Hi Derek,

I can try generate it. My US GGA pseudopotential for
Mg works quite well. It will take few days, hopefully,

including tests.

Bests,
Eyvaz.
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110518/30742ec4/attachment.htm
 


[Pw_forum] Problem in paralization with QE4.3

2011-05-12 Thread S. K. S.
Dear All,

Has anybody implemented  "wfcdir" option in the QE_ph code as like QE_pw
code ?   The exactly same job which was running okay in the older version
of QE (4.2) is not running in the newest version of QE (4.3). Here is the
error :

Calculation of q =0.000   0.000   0.000

 %%
 from diropn : error #21

 %%
 from diropn : error #21
 error opening pbmno.bar
 error opening pbmno.bar
 %

 stopping ...

 %%%
 from diropn : error #21
 error opening pbmno.bar
 %%

 stopping ...
 

 %
 from diropn : error #21
 error opening pbmno.bar
 %

 stopping ...

 stopping ...

 %
 from diropn : error #21
 error opening pbmno.bar
 %%

 stopping ...

--
MPI_ABORT was invoked on rank 3 in communicator MPI_COMM_WORLD
with errorcode 0.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--
--
mpirun has exited due to process rank 5 with PID 4397 on

Looking forward to your comments.

Thanking you and with my best regards,

Saha S. K.
R&D Assistant
JNCASR
Bangalore - 560064
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110512/37aa1739/attachment.htm
 


[Pw_forum] matdyn.x

2011-05-11 Thread S. K. S.
Dear Stefano,

Thanks a lot for your quick reply.
Yes, both set of the frequencies are
obtained after imposing ASR. All of
frequencies are not negative.  Two of
the modes among total 5 modes
(triply generate) are negative, which are
not very unusual in a BaTiO3-type cubic
perovskite.

regards,
Saha S. K.


On 5/11/11, Stefano de Gironcoli  wrote:
> are you imposing acoustic sum rule  ? this change slightly the frequencies .
> isn't strange that your frequencies are all negative (imaginary) ?
> is your system stable ?
> stefano
>


[Pw_forum] matdyn.x

2011-05-11 Thread S. K. S.
 Dear All,

 I find a bit strange results when using matdyn.x to perform the
interpolation in order
 to get  a complete phonon dispersion plot.

 In case of my system, at Gamma point,  ph.x prints out two noticeably
 distinguishable triply degenerate phonon modes of frequency values as follows :

 -18.9595  -18.9595  -18.9595  -11.8835  -11.8835  -11.8835

 However, after using "matdyn.x" the above phonons become like :

 -18.6356  -18.6356  -18.6356  -18.6096  -18.6096  -18.6096

 This looks quite strange to me.  As far as I understand, the phonon
interpolation should
 not change the frequency values at Gamma point or any other
high-symmetry point
 (suppose X in my case) which have been already computed by the
 QE_phonon code. Interpolation means, it should only generate more
frequencies in between
 Gamma and X point, not changing the frequency values in the exact
 Gamma  and X point, which have been already printed  by the ph.x run.

 Looking forward to your comments.

 Thanking you and with my best regards,

 Saha S. K.
 R&D Assistant,
 JNCASR
 Bangalore -560064


[Pw_forum] epsilon_0

2011-04-28 Thread S. K. S.
Dear QE USERS,

What QE calculates is the dielectric tensor at clamped ions
(often referred to as \epsilon_\infty). The ionic contribution
must be calculated and added to this in order to obtain the
static dielectric tensor (\epsilon_0). Does the most recent
version of QE print "\epsilon_0" somewhere in the output file ?


Thanking you and with my best regards,

Saha SK
R&D Assistant
JNCASR, Jakkur,
Bangalore 560064
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110428/b81e5717/attachment.htm
 


[Pw_forum] LDA-USPP for Ge

2011-03-11 Thread S. K. S.
Dear QE Friends,

Does any one have LDA-USPP for Ge ?
Does this scarcity mean that it is not good to use LDA for Ge ?

Thanking you all,

Saha SK
R&D Assistant
JNCASR
Bangalore-560012


[Pw_forum] Ge LDA USPP

2011-02-16 Thread S. K. S.
Dear QE Friends,

Does any one have LDA-USPP for Ge ?

Thanking you,

Saha SK
RA
JNCASR
Bangalore-560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110216/bbb31d29/attachment.htm
 


[Pw_forum] the files necessary for restarting still go into 'outdir'

2011-01-27 Thread S. K. S.
Dear All,

The doc file says :

wfcdir :

this directory specifies where to store files generated by
each processor (*.wfc{N}, *.igk{N}, etc.). The idea here is
to be able to separately store the largest files, while
the files necessary for restarting still go into 'outdir'
(for now only works for stand alone PW )

So here my question goes :

Is this nice trick still works only for stand alone PW or
in the most recent version, this trick has
been extended to other cases, in particular to PH ?

Thanking you and with warm regards,

Saha SK
JNCASR
Bangalore 560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110127/be31ff9c/attachment.htm
 


[Pw_forum] Ferromagnetic Insulator is the Ground State

2011-01-19 Thread S. K. S.
>   Dear Prof. Andrea Dal Corso,
>
>   Thanks a lot for your helpful reply.
>
>   I have further question, what one should do if the ground state is
>   Ferromagnetic  (FM) INSULATOR ? Is there any way ?
>
If the AFM state can be stabilized by symmetry, starting with opposite
starting magnetizations, then the phonon code should work. If you need a
constraint to stabilize the AFM ground state, then presently you cannot
do a phonon calculation.

HTH,

Dear Professor,

Thanks a lot for your kind reply. I am still a bit confused.
Is this your above reply implying that even if I get ferromagnetic
insulator as ground-state, then also I need to follow the above method
(which is
for ANTI-FERROMAGNETIC insulator) to estimate the Born-effective charges ?

Thanking you and with my best regards,
Saha SK
JNCASR
Bangalore 560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110119/a1d9e735/attachment.htm
 


[Pw_forum] Wrong effective charges in AFM insulator using "zue"

2011-01-19 Thread S. K. S.
  Dear Prof. Andrea Dal Corso,

  Thanks a lot for your helpful reply.

  I have further question, what one should do if the ground state is
  Ferromagnetic  (FM) INSULATOR ? Is there any way ?

Andrea wrote:
The phonon with LSDA + constraint is not programmed. So it will not
work. If the AFM insulator is the ground state then you should obtain it
with noncolin=.true. without constraint. In this case ph.x at gamma
should work and you should get the correct effective charges. If not
please provide the complete input.
*

regards,
Saha SK
JNCASR
Bangalore 560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20110119/896782b5/attachment.htm
 


[Pw_forum] spin polarized insulator

2011-01-05 Thread S. K. S.
Dear All,

In the recent version of QE code, is this possible to compute
the dielectric constant of a spin-polarized insulator ?

Thanking you and with my best regards,

Saha SK
JNCASR
Bangalore


[Pw_forum] HOW TO COMPUTE Effective charge for FM INSULATOR using PWSCF?

2010-11-25 Thread S. K. S.
 Dear all,

 The fixed occupations and lsda need nelup and neldw, and
 the phonon code with constrained magnetization is not yet available.
 Does it imply that the phonon calculation  is not possible
 for  fixed occupations and lsda


regards,
S. K. Saha
JNCASR
BANGALORE

P.S. Dear Mighfar, it is okay, cheers!
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20101125/6ee5bcf9/attachment.htm
 


[Pw_forum] cell_dofree

2010-11-25 Thread S. K. S.
Dear All,

Is there any option in PWSCF which allow to move

x, y and z axis, keeping volume unchanged ???

Thanking you all,
Saha S.K.
JNCASR
Bangalore
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20101125/327fda39/attachment.htm
 


[Pw_forum] DFPT+U

2010-07-15 Thread S. K. S.
  Dear QE USERS,

  It is nice to know that the DFPT+U has been recently implemented.
  I hope, it will be public soon. If not, may I ask a developmental
implementation or
  executables, so that I can also perform some test calculations for my
systems ?
  Sometime carrying out FROZEN PHONON is boring.

  By the way, thanks a lot for the good news.

  Regards,
  Saha SK,
  R&D Assistant,
  JNC, Bangalore


On Thu, Jul 15, 2010 at 12:14 PM, Andrea Floris  wrote:

> Dear Ali,
>
> The DFPT+U has been recently implemented, and I am doing some
> calculations on test materials, to see
> how it performs. It will be public in the future, both in
> norm-conserving and ultrasoft pseudopotentials.
>
> Regards,
>
> Andrea
>
>
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100715/64c375cc/attachment.htm
 


[Pw_forum] GGA+U

2010-07-14 Thread S. K. S.
 Dear QE USERS,

 Yes, the code works fine when I change the name
 of element from "eu" to "Eu" inside the pseudopotential file.
 Thank you all for your kind reply.

 One more question : what is Hubbard_alpha ?
  Is this Lagrange multipliers? What is its typical value?
  Is this possible to specify J & U both, in the input file of the newest
version
  of QE code or  is this always U_eff ?

 regards,
 Saha SK
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100714/3c15394b/attachment.htm
 


[Pw_forum] GGA+U

2010-07-13 Thread S. K. S.
Dear QE Users,

I got the following error while used GGA+U for EuO:

 %%
 from set_hubbard_l : error # 1
 pseudopotential not yet inserted
 %%

 stopping ...

Inside source file, I have found that "Eu" is mentioned. However, I do not
understand
the meaning of inserting a pseudopotential.

Thanks in advance for a helpful explanation.

With best regards,

Saha SK
R&D Assistant
JNCASR
Bangalore 560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100713/ca2b7ad1/attachment.htm
 


[Pw_forum] PAW

2010-07-08 Thread S. K. S.
Dear Stefano,

Thanks a lot for your helpful reply. It works fine.
Now the error is  for Li-1s.
 * Li-1s: (18:05:29)...STOP 2
ld1 failed! Aborting...

Should I use  rmatch_augfun= 2.2 in this case ???

In the same way, is it possible to generate Eu PAW using ld1.x ???
With best regards,
Saha SK
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100708/22cf1d72/attachment.htm
 


[Pw_forum] PAW

2010-07-08 Thread S. K. S.
Dear QE USERS,

I got the following error in the version 4.2 of QE :

Running oxygen test...
 * generating
 * testing
 * testing with spin
done

Generating PAW datasets library
 * Be: (16:35:55)... plotting... done
 * C: (16:35:56)... plotting... done
 * Cl: (16:35:56)...STOP 2
ld1 failed! Aborting...

Is there someone else who also has encountered the same error ?

Does anyone has PAW data for other elements, for instance, Eu or Sm ?
Is this possible to use ABINIT PAW inside PWSCF ?

With best regards,
Saha SK
R&D Assistant
JNCASR
India
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100708/fb97a102/attachment.htm
 


[Pw_forum] Eu or Sm Pseudo Potential

2010-06-25 Thread S. K. S.
>  Dear All,
>
>  Does anybody have Eu or Sm (LDA/GGA and Ultra-soft/Norm-conserving)
pseudopotential?
>
>  Looking forward to your kind reply.
>
>  Thanking you and with my best regards,
>
>  Saha SK
>  R&D assistant
>  JNCASR
>  Jakkur 560012
>  India
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100625/9a018c2b/attachment.htm
 


[Pw_forum] Wu-Cohen GGA

2010-06-18 Thread S. K. S.
  Dear All,

  Is  Wu-Cohen GGA implemented in QE?
  Would anyone then like to share his or her pseudopotentials?

  With my best regards,

   Saha S. K.
   R&D Assistant
   JNCASR
   Bangalore- 560012
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100618/eb2c7330/attachment.htm
 


[Pw_forum] flow-chart in PWSCF

2010-04-08 Thread S. K. S.
 Dear All,

 First, let me thank Ngoc Linh Nguyen for his helpful reply.
 A part of his reply is as follows:
In general, the scf processes are similar, but in QE a density
convergence check is applied to accelerate convergence in scf
caculation, i.e. see in PRB 38,  12807 (1988).

This reference is indeed helpful.
I have still one confusion: what, in fact, QE starts with  as
an initial density n(r)  at the 0th iteration, i.e. what is n_0(r) 
Is that simply constant or r dependent or both r and R dependent
Then,  n0(r, R) = what ?
Here, r and R are electronic and nuclei co-ordinates, respectively.

Thanking you and with my best regards,
 Saha SK
 R&D assistant
 JNCASR
 Bangalore
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100408/218f0cb4/attachment.htm
 


[Pw_forum] flow-chart in PWSCF

2010-03-26 Thread S. K. S.
  Dear All,

  At first, let me take this opportunity to thank Prof. Paolo for his
helpful reply.
A part of his reply is given here:

>> I mean, PWSCF check energy convergence or
>> density scf convergence or both?

> the latter. More detials in the QE paper.

  Now just to be fully sure, I would like to ask, does  PWSCF follow
the attached flow-chart (Please see the attached file) where
density-convergence is being checked
On the other hand, in the earlier link:
 http://cms.mpi.univie.ac.at/vasp/vasp/node199.html
energy-convergence had been checked.
 I believe, both the flowcharts  are more or less same within
a technical difference. Am I right?

 Looking forward to your further comments.

 Thanking you and with my best regards,
  Saha SK
   R&D Assistant
  JNCASR
  Bangalore
-- next part --
A non-text attachment was scrubbed...
Name: flowchart.pdf
Type: application/pdf
Size: 65095 bytes
Desc: not available
Url : 
http://www.democritos.it/pipermail/pw_forum/attachments/20100326/6c7086b0/attachment-0001.pdf
 


[Pw_forum] flow-chart in PWSCF

2010-03-17 Thread S. K. S.
Dear All,

http://cms.mpi.univie.ac.at/vasp/vasp/node199.html

   This link has a flow chart describing
the algorithms used in VASP.
Is there such a flow-chart in PWSCF also?
If not, then I would like to know whether PWSCF
follows the same flow-chart or not?
I mean, PWSCF check energy convergence or
density scf convergence or both?
It would be really nice if  somebody send me a full
flow chart  which has been used in PWSCF.

   With my best regards,
   Saha S. K.
   JNCASR
   R&D Assistant
   Bangalore
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20100317/018b91d1/attachment.htm
 


[Pw_forum] Xcrysden

2009-04-03 Thread S. K. S.
Dear all,

Thanks a lot for your kind and helpful reply.
Few things about Xcrysden are yet to be known.

1. Now I can print the thermometer after installing
ImageMagick. However, it prints the thermometer
separately. Please have a look on  the attached file.
Is it possible in XCrysden to print the thermometer  together
with the background  "Xcrysden charge-density" plot


2.  >No! XSF uses Carthesian coordinates (Angstrom units).

 Then,  is there any other smart way to change Angstrom  coordinates
  to crystal coordinates in Xcrysden??? The usual pen and paper
  (analytical and manual) method is not always time-saving especially for
  non-cubic  (for instance, hexagonal) systems.

 3.  Another thing I would like  to know, whether it is possible to increase
  the size of the "XYZ co-ordinate symbol" in Xcrysden???
  One can easily zoom the picture of a molecule or a crystal in the
  Xcrysden and print  the larger size picture because   Xcrysden's print is
  an XWD window dump.  But what about  co-ordinate symbol???
  Is it possible to make its size bigger?

  regards,
  SKS
-- next part --
A non-text attachment was scrubbed...
Name: mol-urea2D.jpg
Type: image/jpeg
Size: 4822 bytes
Desc: not available
Url : 
http://www.democritos.it/pipermail/pw_forum/attachments/20090403/fe5fcd38/attachment-0001.jpg
 


[Pw_forum] Xcrysden

2009-04-03 Thread S. K. S.
Dear Users,

I have a couple of questions  regarding  Xcrysden.
1. Usually xcrysden files use the following format:

 PRIMVEC
5.52460788000.000.00
0.003.90648780400.00
0.000.00   19.5324384880
 PRIMCOORD
 14 1
 31  0.3031580950 0.00-0.3114762280
 31  2.6436116560 1.9532436370-1.8580462820
 31 -0.0518106620 0.00-3.8817233780
 31  2.7623039400 1.9532436370-5.8597314410
 31  0.00 0.00-7.8129756070
 31  2.7623039400 1.9532436370-9.7662192440
 33  1.4763662580 1.9532436370 0.3364242780

where atomic positions  are not written in crystal co-ordinates.
Is it possible to specify atomic positions in crystal co-ordinates in any
of .xyz or .xsf  xcrysden input files (except .pwi or .pwo files).

2. How to print the thermometer in Xcrysden???
Everytime I get an error saying "could not execute /usr/bin/convert".

I'll be highly obliged if anyone kindly discuss the above issues.
regards,
SKS
JNCASR
Bangalore


[Pw_forum] DMFT+LDA

2009-03-22 Thread S. K. S.
Dear all,
Recently, in abinit, I have seen
some on-going developments  on DMFT+LDA.
I'm a little curious. Has  anyone already started  or is
going to start  the similar kind of developmental work in PWSCF
regards,
S. K. Saha
JNCASR
Bangalore


[Pw_forum] USPP

2009-03-17 Thread S. K. S.
> Dear Prithwish,

during next week, I will add to the pseudopotential library
(http://www.pwscf.org/pseudo.htm) the Ce pseudopotential we used to
study cerium-oxide. As you know, in general, compounds based on rare
earths may display quite intriguing properties and our pseudo has not
been tested for systems other than cerias. You can refer to me if you
have specific questions, I will be happy to contribute as far as I
can.

Yours sincerely,
   Stefano 

Dear all,
Quite a long time ago, the above message was posted regarding
Ce pseudopotential.
Is anyone of you having  an  LDA version of the
Ce pseudo which takes f-electron in account?
Looking forward to your kind reply.
regards,
S. K. Saha
JNCASR
Bangalore 560012


[Pw_forum] Silent Calculations

2008-08-11 Thread S. K. S.
  Dear Sir,

  Sorry, for my statements.
  I did not want to report them as an absolute truth.
  My comments was only based on the
  Quantum esspresso code.
  In other packages,  the LRT can be more
  efficient.
  I have not checked them yet.

  Anyway, I would like to thank  Prof. Baroni
  for his excellent reply which will really help to clear
  people's  misconceptions about LRT

  regards
  SKS
  JNCASR
  BANGALORE
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080811/815d1fbc/attachment.htm
 


[Pw_forum] Silent Calculations

2008-08-09 Thread S. K. S.
On Sat, Aug 9, 2008 at 1:52 AM, Paolo Giannozzi wrote:
>  Only for Gamma point phonon  one  can use "Frozen
>  phonon" method because it is much faster.

not true in general. Each linear-response calculation
costs approx. as much as a self-consistent one, and
you need at least as many scf calculations as
linear-response calculations.

Dear Sir,
I agree with you that  for small system each
linear-response calculation
costs approx. as much as a self-consistent one, and
 we need at least as many scf calculations as
linear-response calculations.

However, IN CASE OF  very large system (say about 64 atoms per unit cell),
I found that  a single representation  takes too much computer time
compared to a single scf run.

Moreover, in some cluster which have multiple cpu
for each node and  in some quad or octa core  computer,  always
phonon paralization  is much slower than manual
scf paralization (I mean manually divide   your all  scf calculations in
different  processors).  In case of  Frozen phonon method, one can
also accomplish the same goal  by using less memory and
disk space.
However, for non-zone center phonon Linear response is the best way.
There is still a scope in PWSCF to make Phonon calculation much faster
by implementing the representation number parallization.

regards,
SKS
JNCASR
BANGALORE


On Sat, Aug 9, 2008 at 1:52 AM, Paolo Giannozzi wrote:

>
> On Aug 8, 2008, at 18:41 , S. K. S. wrote:
>
> >  Only for Gamma point phonon  one  can use "Frozen
> >  phonon" method because it is much faster.
>
> not true in general. Each linear-response calculation
> costs approx. as much as a self-consistent one, and
> you need at least as many scf calculations as
> linear-response calculations
>
> Paolo
> ---
> Paolo Giannozzi, Dept of Physics, University of Udine
> via delle Scienze 208, 33100 Udine, Italy
> Phone +39-0432-558216, fax +39-0432-558222
>
>
>
>
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080809/b3409ab5/attachment.htm
 


[Pw_forum] Silent Calculations

2008-08-08 Thread S. K. S.
 Dear Miguel,

 Yes, you are right. The linear response calculation needs
 much time even for calculating phonons at Gamma point.
 Only for Gamma point phonon  one
 can use "Frozen
 phonon" method because
 it is much faster.

 The calculations can  be got stalled for many reasons.
 You have not mentioned whether at last
 (I mean  staying  in the silence mode for few weeks) your jobs
 properly (without any error) get over or not???

 If  yes, then it can be due to the problem in your
 parallel computer - may be  out  of  16 processors
 some particular processor was  down for some time.

 If  the  answer of  the above question is no,  then  the
 memory of your computer is not sufficient
 enough to tackle a large system with lot of atoms
 because it needs large memory and  as well as
 generates a huge tmp file.

 All  input parameters
 (deguass, la2F, filedvscf ..)  also should be given correctly  for
 electron-phonon calculation.

 regards,
 SKS
 JNCASR
 Bangalore


On 8/8/08, Miguel Mart?nez  wrote:
>
> Dear everybody,
>
> I've been running into a few issues as of late, when running some
> calculations, it seems the calculations have stalled for a pretty long
> period of time.
>
> The one that troubles me less, taking into account that it's pretty
> expensive, is an electron-phonon calculation. This is a pretty large
> system (20 atoms in the primitive cell) and computing the phonons at
> gamma took ages (well two days using 4 processors). It would be great
> if some output were written just to remind you the code is still
> working, maybe something related to convergence. But I find somewhat
> disturbing that no line was written since the electron phonon
> calculation began.
>
> My other issue is a very similar stalling in a Raman calculation. No,
> it's not the same structure. Here, I lack from experience, since this
> is the first Raman calculation I've done outside the examples.
>
> Self-consistent input is:
> &control
>calculation='scf'
>restart_mode='from_scratch',
>prefix='sih4',
>tstress = .true.
>tprnfor = .true.
>pseudo_dir = '/home/wmbmacam/pseudo/',
>outdir='/scratch/wmbmacam/fdd2-raman/'
> /
> &system
>ibrav= 10, celldm(1) = 12.65898481, celldm(2) = 1.07100858,
> celldm(3) = 0.621208699, nat=10, ntyp= 2,
>ecutwfc = 55.0
>nr1= 72, nr2= 72, nr3= 72
> /
> &electrons
>conv_thr = 1.0d-8
>mixing_beta = 0.7
> /
> ATOMIC_SPECIES
> Si  28.0855 Si.pbe-rrkj.UPF
> H   1.0079  H.pbe-vbc.UPF
> ATOMIC_POSITIONS {crystal}
> Si 0.614368447799 0.385631552201 0.614368447799
> Si 0.864368447799 0.635631552201 0.864368447799
> H  0.442372877617 0.836174882136 0.364427192728
> H  0.364427192728 0.357025047520 0.442372877617
> H  0.892974952480 0.885572807272 0.413825117864
> H  0.413825117864 0.807627122383 0.892974952480
> H  0.378768583324 0.338682246661 0.003482898532
> H  0.003482898532 0.279066271483 0.378768583324
> H  0.970933728517 0.246517101468 0.911317753339
> H  0.911317753339 0.871231416676 0.970933728517
> K_POINTS {automatic}
> 12 12 12  1 1 1
>
>
> The phonon input is:
>
> phonons of Sih4 at Gamma
> &inputph
> tr2_ph=1.0d-16,
> prefix='sih4',
> epsil=.true.,
> trans=.true.,
> lraman=.true.,
> elop=.true.,
> amass(1)=28.0855,
> amass(2)=1.00794,
> outdir='/scratch/wmbmacam/fdd2-raman/',
> fildyn='sih4.dynG',
> fildrho='sih4.drho',
> /
> 0.0 0.0 0.0
>
>
> The calculations stalled right after the effective charges calculations.
> After that, it spent 24 hours without writing anything (it needed 10
> hours to compute the effective charges). No convergence, no threshold,
> no crashes, no nothing.
>
> Is this behaviour normal? Or are my compiled binaries just telling me I
> should stop bothering them in August?
>
> Finally, I just wanted to note that the amasses in example 15 are
> seemingly wrong. While the masses in the scf run are fine, the ones
> used in the phonon calculations correspond to GaAs instead of AlAs.
>
> Regards,
>
> Miguel
>
> --
> 
> Miguel Mart?nez Canales
>   Dto. F?sica de la Materia Condensada
>   UPV/EHU
>   Facultad de Ciencia y Tecnolog?a
>   Apdo. 644
>   48080 Bilbao (Spain)
> Fax:  +34 94 601 3500
> Tlf:  +34 94 601 5326
> 
>
> "If you have an apple and I have an apple and
> we exchange these apples then you and I will
> still each have one apple. But if you have an
> idea and I have an idea and we exchange these
> ideas, then each of us will have two ideas."
>
> George Bernard Shaw
>
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080808/77a0fabe/attachment.htm
 


[Pw_forum] not an easy system to treat

2008-05-29 Thread S. K. S.
  Dear Adriano,

  Thanks a lot for your kind help.
  I'll wait for your next reply.

  regards,
  SKS
  JNCASR
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080529/90512d69/attachment.htm
 


[Pw_forum] nelup and neldw

2008-05-29 Thread S. K. S.
 Dear all,

 The fixed occupations and lsda need nelup and neldw, and
 the phonon code with constrained magnetization is not yet available.
 Does it imply that the phonon calculation  is not possible
 for  fixed occupations and lsda


regards,
S. K. Saha
JNCASR
BANGALORE
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080529/9b526b28/attachment.htm
 


[Pw_forum] XCRYSDEN

2008-05-22 Thread S. K. S.
  Dear PWscf Users,

  Is there any option in XCRYSDEN software where
  we can  increase the size or the font of  the
 "Coordinate System"  logo which usually
  comes below in the
  left side corner of any XCRYSDEN Picture?

  regards,

  S. K. S.
  JNCASR
  Bangalore
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080522/5a5411ed/attachment.htm
 


[Pw_forum] format UPF

2008-04-07 Thread S. K. S.
  Dear All,

  Is there any already
  exits  converter to
  convert a pseudopotential,
  generated by Octopus,
  into UPF format??

  regards,
  SKS
  JNCASR
  India
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080407/a7fb9ae2/attachment.htm
 


[Pw_forum] Gd LDA pseudo

2008-02-28 Thread S. K. S.
Dear Sir,

Thanks a  lot  for  your  kind reply.

 First time,  I'm  trying to generate a pseudo myself.

 Being newbie, some thing I can not understand clearly.

 As for example, in pwscf suite it's given for Os that
*
  &input
title='Os',
zed=76.,
rel=1,
iswitch=3,
config='[Xe] 4f14.0 5d6.0 6s2.0 6p0.0',
dft='LDA'
 /
 &inputp
   lloc=0,
   pseudotype=3,
   file_pseudopw='Os.pz-rrkjus.UPF',
   nlcc=.true.,
   rcore=0.9,
   author='ADC',
 /
4
6P  2  1  0.00  0.00  2.60  2.60  1
5D  3  2  6.00  0.00  1.90  2.30  1
5D  3  2  0.00 -0.20  1.90  2.30  1
6S  1  0  2.00  0.00  2.50  2.50  1
*
Though I can understand the meaning of the most of the input parameters,
however,  in  last  part , why  the numbers are  like these

6P   0.00  2.60  2.60  1
5D   0.00  1.90  2.30  1
5D  -0.20  1.90  2.30  1
6S   0.00  2.50  2.50  1  , I can not understand clearly.

I know that for Gd it will be like that

6P   21   0.00
5D   32   1.00
5D   320.00
6S   102.00

but what will be the other numbers for Gd


I'll be highly obliged if somebody kindly help me to understand
the above points. First,  I'll test the generated Gd pseudo and if
it works well  I'll  defintely  give  it  back to PWSCF website.

Look forward to your kind help.

With my best regards,
SKS
JNCASR
BANGALORE




-- Forwarded message --
From: Dal Corso Andrea 
Date: Feb 27, 2008 8:29 PM
Subject: Re: Gd LDA pseudo
To: "S. K. S." 


Dear S.K.S.

Presently I do not have a Gd pseudopotential.
Please post this question to the pw_forum mailing list, maybe someone
will be able to provide it.

Best wishes

Andrea




-- Forwarded message ----------
From: Dal Corso Andrea 
Date: Feb 27, 2008 8:29 PM
Subject: Re: Gd LDA pseudo
To: "S. K. S." 

Dear S.K.S.

Presently I do not have a Gd pseudopotential.
Please post this question to the pw_forum mailing list, maybe someone
will be able to provide it.

Best wishes

Andrea

On Wed, 2008-02-27 at 20:19 +0530, S. K. S. wrote:
>
>
>Respected  Sir,
>
>May I get a  Gd pseudo (LDA) for PWSCF  from
>you or any of your students???
>
>I know that the rare earth is poorly described
>by the plane waves. However, there are lot of DFT  calculations
>on Lanthanides  and Ceria  in literature.
>
>Looking forward to your kind help.
>
>   With my best regards,
>   SKS
>   JNCASR
>   Bangalore
>
>   P.S.   I'll be very happy to get someone  as a coauthor if he/she
>   helps me by providing the above pseudo.

--
Andrea Dal CorsoTel. 0039-040-3787428
SISSA, Via Beirut 2/4   Fax. 0039-040-3787528
34014 Trieste (Italy)   e-mail: dalcorso at sissa.it


On Wed, 2008-02-27 at 20:19 +0530, S. K. S. wrote:
>
>
>Respected  Sir,
>
>May I get a  Gd pseudo (LDA) for PWSCF  from
>you or any of your students???
>
>I know that the rare earth is poorly described
>by the plane waves. However, there are lot of DFT  calculations
>on Lanthanides  and Ceria  in literature.
>
>Looking forward to your kind help.
>
>   With my best regards,
>   SKS
>   JNCASR
>   Bangalore
>
>   P.S.   I'll be very happy to get someone  as a coauthor if he/she
>   helps me by providing the above pseudo.

--
Andrea Dal CorsoTel. 0039-040-3787428
SISSA, Via Beirut 2/4   Fax. 0039-040-3787528
34014 Trieste (Italy)   e-mail: dalcorso at sissa.it
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20080228/7cae5887/attachment-0001.htm
 


[Pw_forum] Vanderbilt USPP

2007-12-14 Thread S. K. S.
 Dear Users,

 Can anyone please tell me whether  the Vanderbilt USPP for
 Bi,   Sm,  Gd and Dy  are available or not ?

 Regards,

  SKS
  PhD
 Jncasr
  India
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20071214/1aeb2595/attachment.htm
 


[Pw_forum] Doc/INPUT_PW: A compensating jellium background is inserted

2007-11-18 Thread S. K. S.
  Dear all,

  It's written in the documentation that
  a compensating jellium background is inserted.

  I'll be highly obliged if someone kindly tell me the name
  of the subroutine where it has been inserted.

  Yours sincerely,

  S. K. S
  R&D Assistant.
 Jncasr
  India
-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/20071118/6897937c/attachment.htm