[QE-users] convergence problem
Dr. Thomas Brumme
thomas.brumme at uni-leipzig.de
Tue Mar 6 20:49:42 CET 2018
Well, the error says what's the problem... You can't use the option
lefield for
a metallic system. Depending on what you want to investigate you can try with
the sawtooth approach using tefield/dipfield.
Regards
Thomas
Zitat von Sohail Ahmad <sohailphysics at yahoo.co.in>:
> Thanks Dr Thomas BrummeAs suggested,when i changed occuptaions =
> 'fixed' to 'smearning'
> it gives following error
> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
> Error in routine iosys (1):
> Berry Phase/electric fields only for insulators!
> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
>
> &control
> calculation = 'relax',
> restart_mode = 'from_scratch',
> pseudo_dir= '/home/sohail/pseudo',
> outdir='./OUT',
> prefix='PdS2bAAel02',
> lelfield = .true.,
> gdir = 3,
> etot_conv_thr = 1.0d-5,
> forc_conv_thr = 1.0d-4,
> /
> &system
> ibrav = 4, a = 3.53, b = 3.53, c = 20.0, cosAB = -0.5, cosAC
> = 0.0, cosBC = 0.0,
> nat = 6, ntyp = 2,
> ecutwfc = 80,
> ecutrho = 400,
> nbnd = 40,
> occupations = 'smearing',smearing = 'gaussian', degauss = 0.001,
> /
> &electrons
> mixing_beta = 0.3,
> conv_thr = 1.0d-9,
> electron_maxstep = 900,
> efield_cart(1) = 0.0d0,
> efield_cart(2) = 0.0d0,
> efield_cart(3) = 0.0054d0,
> /
> &ions
> ion_dynamics = 'bfgs',
> /
> ATOMIC_SPECIES
> Pd 106.42 Pd.pbe-mt_fhi.UPF
> S 32.06 S.pbe-mt_fhi.UPF
> ATOMIC_POSITIONS angstrom
> Pd 0.00000000
> 0.00000000 0.00000000
> S 1.76627922
> 1.01975766 1.28191684
> S 1.76627922
> -1.01975766 -1.28191684
> Pd 0.00000000
> 0.00000000 6.60000000
> S 1.76627922
> 1.01975766 7.88191684
> S 1.76627922
> -1.01975766 5.31808316
> K_POINTS AUTOMATIC
> 12 12 4 0 0 0
>
>
>
> On Tuesday, 6 March 2018 6:11 PM, Dr. Thomas Brumme
> <thomas.brumme at uni-leipzig.de> wrote:
>
>
> Dear Sohail,
>
> bilayer PdS2 becomes metallic and thus you need to change from
> occupations='fixed' to 'smearing'...
> See also:
>
> http://onlinelibrary.wiley.com/doi/10.1002/anie.201309280/full
>
> Regards
>
> Thomas Brumme
>
> Zitat von Sohail Ahmad <sohailphysics at yahoo.co.in>:
>
>> I wish to apply electric field in Z direction to bilayers but it
>> didnot get converge aven after 5 days....while in case of monolayer
>> its fine...
>> Please rectify the mistakes in the following
>> input----------------------------------------------------------------------------------------------
>> &control
>> calculation = 'relax',
>> restart_mode = 'from_scratch',
>> pseudo_dir= '/home/sohail/pseudo',
>> outdir='./OUT',
>> prefix='PdS2bAAel04',
>> lelfield = .true.,
>> gdir = 3,
>> etot_conv_thr = 1.0d-5,
>> forc_conv_thr = 1.0d-4,
>> /
>> &system
>> ibrav = 4, a = 3.53, b = 3.53, c = 20.0, cosAB = -0.5, cosAC =
>> 0.0, cosBC = 0.0,
>> nat = 6, ntyp = 2,
>> ecutwfc = 80,
>> ecutrho = 400,
>> nbnd = 40,
>> occupations = 'fixed',
>> /
>> &electrons
>> mixing_beta = 0.3,
>> conv_thr = 1.0d-9,
>> electron_maxstep = 900,
>> efield_cart(1) = 0.0d0,
>> efield_cart(2) = 0.0d0,
>> efield_cart(3) = 0.0108d0,
>> /
>> &ions
>> ion_dynamics = 'bfgs',
>> /
>> ATOMIC_SPECIES
>> Pd 106.42 Pd.pbe-mt_fhi.UPF
>> S 32.06 S.pbe-mt_fhi.UPF
>> ATOMIC_POSITIONS angstrom
>> Pd 0.00000000
>> 0.00000000 0.00000000
>> S 1.76627922
>> 1.01975766 1.28191684
>> S 1.76627922
>> -1.01975766 -1.28191684
>> Pd 0.00000000
>> 0.00000000 6.60000000
>> S 1.76627922
>> 1.01975766 7.88191684
>> S 1.76627922
>> -1.01975766 5.31808316
>> K_POINTS AUTOMATIC
>> 16 16 4 0 0 0
>> ---------------------------------------------------------------------
>> Dr. Sohail AhmadDepartment of PhysicsKing Khalid UniversityAbha,
>> Saudi
>> Arabia--------------------------------------------------------------------
>
> --
> Dr. rer. nat. Thomas Brumme
> Wilhelm-Ostwald-Institute for Physical and Theoretical Chemistry
> Leipzig University
> Phillipp-Rosenthal-Strasse 31
> 04103 Leipzig
>
> Tel: +49 (0)341 97 36456
>
> email: thomas.brumme at uni-leipzig.de
More information about the users
mailing list