[Pw_forum] vcut_get; q vector out of the grid ERROR

[Paolo Giannozzi]

I am not sure that the  exxdiv_treatment='vcut_ws' machinery
is assumed to work during variable-cell optimization.

P.

On Thu, 2014-05-08 at 22:24 -0400, Ankit wrote:
> Hello QE users and developers,
> 
> I am trying to relax a polymer which is periodic in one-dim and I am 
> using vacuum in the other two directions. I am using input_dft='B3LYP'. 
> The code was complaining to check for exxdiv_treatment and I realised 
> that 'vcut_ws' is more appropriate for my system as it is anisotropic 
> cell. But with exxdiv_treatment='vcut_ws', I am getting another error 
> "Error in routine vcut_get; q vector out of the grid ". I am not able to 
> find the cause of this error on mailing list.
> 
> I am wondering if someone can help me with this error.
> 
> My input file looks like:
> 
>   &control
>    title='Germanium'
>    calculation='vc-relax'
>    restart_mode='from_scratch'
>    outdir='/scratch'
>    pseudo_dir='./'
>    prefix='GEPH_kgrid_116'
>    tprnfor=.true.
>    tstress=.true.
>    !etot_conv_thr=1.0d-9
>    !forc_conv_thr=1.0d-4
>    nstep = 100
>   /
>   &system
>    ibrav=0
>    celldm(1)=4.8377
>    nat=6
>    ntyp=2
>    ecutwfc=100
>    input_dft='B3LYP'
>    exxdiv_treatment='vcut_ws'
>    ecutvcut=0.7
>    x_gamma_extrapolation=.false.
>    occupations='smearing'
>    degauss=0.01
>   /
>   &electrons
>    conv_thr=1.0d-12
>    !mixing_beta=0.5
>    diagonalization='cg'
>   /
>   &ions
>     !ion_nstepe=20
>     ion_dynamics='damp'
>   /
>   &cell
>     cell_dynamics='damp-pr'
>     cell_dofree='z'
>     press_conv_thr=0.1d0
>   /
> ATOMIC_SPECIES
>   CA 12.0107 C.blyp-mt.UPF
>   HA 1.00794 H.blyp-vbc.UPF
> ATOMIC_POSITIONS crystal
> CA       0.584592955   0.580716431   0.000001513
> HA       0.637696017   0.569848961   0.000000769
> HA       0.580882017   0.633258175   0.000001691
> CA       0.553589442   0.554559138   0.499998487
> HA       0.500486682   0.565428771   0.499999231
> HA       0.557298406   0.502017244   0.499998309
> K_POINTS automatic
>   1 1 12 0 0 0
> 
> 
> Thanks,
> Ankit Jain
> PhD Candidate,
> IIT Indore
> India
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://pwscf.org/mailman/listinfo/pw_forum

-- 
 Paolo Giannozzi, Dept. Chemistry&Physics&Environment, 
 Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
 Phone +39-0432-558216, fax +39-0432-558222