[Pw_forum] Cell geometry optimization issue

Giovanni Cantele giovanni.cantele at spin.cnr.it
Fri Nov 6 15:16:03 CET 2015


As explained here : http://www.quantum-espresso.org/faq/self-consistency/#6.2
before changing any calculation parameter such as mixing_beta it is advisable to
check if there is a problem with band input data, that means unit cell and/or atomic positions.

Just by visualising your input using XcrysDen you can easily find out that, for example,
atom 4
W   0.54140000000  0.0000000000  0.98510000000
and atom 22
W   0.45860000000  0.0000000000  0.98510000000

are just 0.5 A from each other : d ~ (0.54-0.46)*13*0.5291, a distance maybe too small to avoid self-consistency issues.
Maybe a similar, even though less critical issue is with atoms 7 and 13 (only 1.36 A).

So, before going on, I would check if your input structure is exactly the one you want to study, and if the input coordinates
are right in the right units.

Giovanni

> On 06 Nov 2015, at 14:59, Arijan Aleksić <aaleksic at fizika.unios.hr> wrote:
> 
> Dear all,
> 
> I am working on geometry optimization for WTe2 (Tungsten(IV) telluride).
> When running calculations, they failed to converge within limits in 100
> iterations.
> 
> First I have increased cell parameters for 10-15% to check if there will
> be any difference (also had big pressure). Next I went with reducing
> mixing_beta in steps of 0.1 to check for any changes which yielded no good
> results. And finally I used calculation = 'vc-relax' which gave even worse
> results and could not optimize the cell.
> 
> I will accept any advice on what should be my next action in order to
> succeed in getting optimization complete.
> 
> Here is the text from input file I used (version of Quantum Espresso is
> 5.2.1):
> 
> &control
>    calculation  = 'vc-relax',
>    restart_mode = 'from_scratch',
>    prefix       = 'wte2',
>    tprnfor      = .true.,
>    pseudo_dir   = '/home/aaleksic/Pseudo/',
>    outdir       = '/home/aaleksic/wte2/tmp/',
> /
> 
> &system
>    ibrav     = 8,
>    celldm(1) = 13,
>    celldm(2) = 0.53846,
>    celldm(3) = 2.15385,
>    nat       = 24,
>    ntyp      = 2,
>    nbnd      = 74,
>    ecutwfc   = 100.0,
> /
> 
> &electrons
>    mixing_beta = 0.1,
>    conv_thr    = 1.0d-12,
> /
> 
> &ions
> /
> 
> &cell
> /
> 
> ATOMIC_SPECIES
> W   74.0   W.pbe-mt_fhi.UPF
> Te  52.0   Te.pbe-mt_fhi.UPF
> 
> ATOMIC_POSITIONS (crystal)
> W   0.90050000000  0.5000000000  0.00000000000
> Te  0.29410000000  0.5000000000  0.09650000000
> Te  0.80020000000  0.0000000000  0.14000000000
> W   0.54140000000  0.0000000000  0.98510000000
> Te  0.35590000000  0.0000000000  0.34490000000
> Te  0.85170000000  0.5000000000  0.38930000000
> W   0.59950000000  0.5000000000  0.50000000000
> Te  0.20590000000  0.5000000000  0.59650000000
> Te  0.69980000000  0.0000000000  0.64000000000
> W   0.95860000000  0.0000000000  0.48510000000
> Te  0.14410000000  0.0000000000  0.84490000000
> Te  0.64830000000  0.5000000000  0.88930000000
> W   0.40050000000  0.5000000000  0.50000000000
> Te  0.79410000000  0.5000000000  0.59650000000
> Te  0.30020000000  0.0000000000  0.64000000000
> W   0.04140000000  0.0000000000  0.48510000000
> Te  0.85590000000  0.0000000000  0.84490000000
> Te  0.35170000000  0.5000000000  0.88930000000
> W   0.09950000000  0.5000000000  0.00000000000
> Te  0.70590000000  0.5000000000  0.09650000000
> Te  0.19980000000  0.0000000000  0.14000000000
> W   0.45860000000  0.0000000000  0.98510000000
> Te  0.64410000000  0.0000000000  0.34490000000
> Te  0.14830000000  0.5000000000  0.38930000000
> 
> 
> 
> K_POINTS (automatic)
> 4 4 4 0 0 0
> 
> Thank you.
> 
> Best regards,
> Arijan Aleksic,
> University J.J.Strossmayer, Croatia.
> 
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://pwscf.org/mailman/listinfo/pw_forum
> 

-- 

Giovanni Cantele, PhD
CNR-SPIN
c/o Dipartimento di Fisica
Universita' di Napoli "Federico II"
Complesso Universitario M. S. Angelo - Ed. 6
Via Cintia, I-80126, Napoli, Italy
e-mail: giovanni.cantele at spin.cnr.it
Phone: +39 081 676910
Skype contact: giocan74

ResearcherID: http://www.researcherid.com/rid/A-1951-2009
Web page: http://people.na.infn.it/~cantele






More information about the users mailing list