[Pw_forum] Is maxirr=N working in phonon calculations?
Serge Nakhmanson
nakhmanson at anl.gov
Wed Nov 5 02:05:56 CET 2008
Dear All,
I run a small test with a 5-atom PbTiO3 cell in P4/mmm
symmetry, where it is supposed to have 10 irreps (i.e.,
each 3d irrep in Pm-3m symmetry splits into a 1d+2d)
and I asked the program to compute only 5 of these
as shown below:
> node02 {1003} > more PbTiO3.tp-relaxed-STO-t.ph.G.inp
> 'tp-STO-t-PbTiO3, phonons at G'
> &inputph
> prefix = 'tp-STO-t-PbTiO3'
> amass(1) = 207.20000
> amass(2) = 47.86700
> amass(3) = 16.00000
> tr2_ph = 1e-20
> niter_ph = 100
> alpha_mix = 0.5
> maxirr = 5
> trans = .true.
> epsil = .true.
> fildyn = 'PbTiO3.tp-relaxed-STO-t.G.dyn'
> iverbosity = 1
> outdir='./'
> /
> 0.0 0.0 0.0
> node02 {1004} >
Still the code blasted through all 10 irreps. I grep-ped
out "maxirr" in /PH/*.f90 and saw that it is being set and reset
in multiple places (I cannot trace it without my head spinning),
so I decided to ask a question about this instead. Of course,
there is always a chance that I made a stupid error in my input,
so PLZ be kind to me :))
Also, is there any system to which irreps are computed first and
which later? Or is this process completely random, i.e. it would
be futile to cover a certain frequency window by computing only
some, but not all, irreps?
THX,
Serge
--
*********************************************************
Serge M. Nakhmanson phone: (630) 252-5205
Assistant Scientist fax: (630) 252-4798
MSD-212, Rm. C-224
Argonne National Laboratory
9700 S. Cass Ave.
Argonne, IL 60439
*********************************************************
More information about the Pw_forum
mailing list