[Pw_forum] symmetry operation not allowed
Eyvaz Isaev
eyvaz_isaev at yahoo.com
Fri Apr 18 15:10:38 CEST 2008
Hi,
You can set up nr1, nr2, and nr3 manually.
Or change ecutrho gradually.
Also provide your affiliation, please.
Bests,`
Eyvaz.
--- Lily Anh <lily_physics at yahoo.com.sg> wrote:
> Dear all,
> I know this is a old question and I also read the
> previous
> answer about this.
> I test with different ecutwfc , ecutrho (until 30
> times of ecutwfc) , and
> K-grid. However, the problem still exists. Since
> the calculation expense is always in
> concern, so I want to use symmetry during the
> calculation.
> I attach the output and input file here for
> reference.
> Output file:
> Current
> dimensions of program pwscf are:
> ntypx =10 npk =40000 lmax = 3
> nchix = 6 ndmx = 2000 nbrx =14 nqfx = 8
> warning: symmetry
> operation # 2 not allowed. fractional
> translation:
> -0.4999900
> -0.4999400 -0.5000000 in crystal
> coordinates
> warning: symmetry
> operation # 3 not allowed. fractional
> translation:
> 0.0000100
> -0.5000000 -0.2498700 in crystal
> coordinates
> warning: symmetry
> operation # 4 not allowed. fractional
> translation:
> -0.5000000 0.0000600 0.2501300 in crystal
> coordinates
>
> Planes per
> process (thick) : nr3 = 90 npp = 8
> ncplane =32400
>
> Planes per
> process (smooth): nr3s= 32 npps= 3
> ncplanes= 4096
>
> Proc/ planes cols G planes cols G
> columns G
> Pool (dense grid) (smooth grid)
> (wavefct grid)
> 1 8 1960 109526 3 261 5323 76
> 838
> 2 8 1961 109525 3 261 5333 75
> 837
> 3 8 1961 109525 3 261 5325 75
> 837
> 4 8 1961 109525 3 261 5321 75
> 837
> 5 8 1961 109525 3 261 5321 77
> 837
> 6 8 1961 109525 3 261 5321 77
> 837
> 7 7 1960 109524 3 262 5338 77
> 837
> 8 7 1960 109524 3 261 5339 77
> 837
> 9 7 1960 109524 2 261 5337 77
> 837
> 10 7 1960
> 109524 2 261 5337 77 837
> 11 7 1960 109524 2 261 5335 77
> 837
> 12 7 1960 109524 2 261 5325 77
> 837
> 0 90 235251314295 32 3133 63955 917
> 10045
>
>
>
> bravais-lattice
> index = 8
> lattice parameter
> (a_0) = 15.6923 a.u.
> unit-cell
> volume = 1872.5343 (a.u.)^3
> number of
> atoms/cell = 36
> number of atomic
> types = 2
> kinetic-energy
> cutoff = 40.0000 Ry
> charge density
> cutoff = 1200.0000 Ry
> convergence
> threshold = 1.0E-08
> beta = 0.7000
> number of
> iterations used = 8 plain mixing
> Exchange-correlation = SLA PW PBE PBE
> (1434)
> iswitch = 0
>
> celldm(1)= 15.692292 celldm(2)= 1.000000
> celldm(3)= 0.484586
> celldm(4)= 0.000000 celldm(5)= 0.000000
> celldm(6)= 0.000000
>
> crystal axes:
> (cart. coord. in units of a_0)
> a(1) =
> ( 1.000000 0.000000 0.000000 )
> a(2) =
> ( 0.000000 1.000000 0.000000 )
> a(3) =
> ( 0.000000 0.000000 0.484586 )
>
> reciprocal axes:
> (cart.. coord. in units 2 pi/a_0)
> b(1) = ( 1.000000 0.000000
> 0.000000 )
> b(2) =
> ( 0.000000 1.000000 0.000000 )
> b(3) =
> ( 0.000000 0.000000 2.063618 )
>
>
> PSEUDO 1 is
> O (US) zval = 6.0 lmax= 2 lloc= 0
> Version 0 0 0 of US pseudo code
> Using log mesh
> of 1269 points
> The
> pseudopotential has 4 beta functions
> with:
> l(1)
> = 0
> l(2)
> = 0
> l(3)
> = 1
> l(4)
> = 1
> Q(r) pseudized
> with 0 coefficients, rinner = 0.000 0.000
> 0.000
>
> 0.000 0.000
>
> PSEUDO 2 is
> H (US) zval = 1.0 lmax= 1 lloc= 0
> Version 0 0 0 of US pseudo code
> Using log mesh
> of 1061 points
> The
> pseudopotential has 2 beta functions
> with:
> l(1)
> = 0
> l(2)
> = 0
> Q(r) pseudized
> with 0 coefficients, rinner = 0.000 0.000
> 0.000
>
>
> atomic
> species valence mass pseudopotential
> O 6.00 15.99940 O (
> 1.00)
> H 1.00 1.00790 H (
> 1.00)
>
> No symmetry!
> INPTU file:
> &control
> calculation='scf',
> restart_mode='from_scratch',
> prefix='h2o'
> pseudo_dir = '/pwwork/pseudo/',
> outdir=/pwwork/tmp/tmp1/'
> tstress=.t..,
> tprnfor=.t.
> /
> &system
> ibrav = 8,
> celldm(1) = 15.69229199303824618227,celldm(2)
> =1.00,celldm(3)
> =.48458574181117533718,nat= 36, ntyp= 2,
> ecutwfc=40,ecutrho=1200
> /
> &electrons
> mixing_beta = 0.7
> conv_thr = 1.0d-8
> /
> ATOMIC_SPECIES
> O 15.9994 O.pbe-rrkjus.UPF
> H 1.0079 H.pbe-rrkjus.UPF
> ATOMIC_POSITIONS {crystal}
> O -0.00178 0.00473 0.00365
> O 0..00177 0.50473 0.24622
> O 0.13861 0.25239 0.61888
> O 0.25932 0.63953 0.89389
> O 0.24067 0.86042 0.39389
> O 0.36138 0..24756 0.11888
>
=== message truncated ===>
_______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
-------------------------------------------------------------------
Prof. Eyvaz Isaev,
Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia,
IFM, Linkoping University, Sweden
Condensed Matter Theory Group, Uppsala University, Sweden
Eyvaz.Isaev at fysik.uu.se, eyvaz_isaev at yahoo.com
____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
More information about the Pw_forum
mailing list