[Pw_forum] Hubbard U and small degauss ...
Gabriele Sclauzero
sclauzer at sissa.it
Mon May 18 10:12:25 CEST 2009
Dear Naol,
Naol Regassa wrote:
> Dear all
>
> I have two questions to ask
>
> a) For the last times I was calculating the effect of Hubbard U on the
> system YMnO_3 using QESPRESSO 4.10. I got an insulator band structure
> using U = 10.5eV though it is completely metal without U (GGA). When I
This is plausible. Though I don't know the physics of your compound, I can tell you that
for simpler transition metal oxides (like FeO) you find the same change (i.e. from metal
to insulator) after introducing the Hubbard U. Please have a look at Matteo Cococcioni's
PhD thesis (from http://www.sissa.it/cm/phd.php) and the related paper (PRB 2005).
> got this result I have used small degauss = 0.006. My friend told me
> to use degauss = 0.015.
Have you asked him the reason? I don't see any...
In principle you should test convergence of the total energy (and other quantities you're
interested in) w.r.t all the parameters which introduce significant approximations to your
result. As you decrease degauss you'll reach the correct result for your system at T=0
(chosing k-point grids accordingly), so I don't see any reason to choose higher degauss
other that saving computer time.
In the specific case of DFT+U, you have to be carefull to converge the local occupations
too. They are very sensitive to the value of ecutwfc, but also degauss and number k-points
have to be checked.
Please tell me how possible (correct) is
> that? I checked the convergence for degauss but it is increasing near
> 0.015.
What is increasing?
>
> b) Is adding U value more effective for the most stable structure ?
Sorry I cannot get this question... in which sense "more effective"?
> For the unstable states like for my case cubic ferromagnetic phase I
> can't see any change?
What are looking at? Maybe your question is clear to someone who's studying systems like
yours, so it's because of me that I cannot get the point... but if you could be more
precise maybe you would have more change to get help.
Regards,
GS
>
> Here is the input for the 'scf'
> &control
> calculation = 'scf',
> restart_mode = 'from_scratch',
> prefix = 'ymo_hcp_ferro'
> tstress = .true.,
> tprnfor = .true.,
> pseudo_dir = '/home/MSC608/Desktop/QESPRESSO_PP',
> outdir='/home/MSC608/Thesis/Hexagonal/Ferro-Hex/temp'
> /
> &system
> ibrav = 0,
> celldm(1) = 6.7,
> nat = 10,
> ntyp = 4,
> ecutwfc = 32,
> ecutrho = 340.0
> nspin = 2
> starting_magnetization(2) = 0.7,
> starting_magnetization(3) = 0.7,
> occupations = 'smearing',
> smearing = 'marzari-vanderbilt',
> degauss = 0.005,
> lda_plus_u=.true.,
> Hubbard_U(2) = 10.5,
> Hubbard_U(3) = 10.5,
> /
> &electrons
> diagonalization = 'david'
> startingwfc = 'atomic+random'
> mixing_beta = 0.3,
> conv_thr = 1.0d-10
> /
> ATOMIC_SPECIES
> Y 88.91 Y.pbe-nsp-van.UPF
> Mn1 54.94 Mn.pbe-sp-van_mit.UPF
> Mn2 54.94 Mn.pbe-sp-van_mit.UPF
> O 15.9994 O.pbe-rrkjus.UPF
>
> CELL_PARAMETERS (alat)
> 0.985236441 -0.266531554 -0.000316189
> -0.266366323 0.985859979 0.000349871
> 0.000781869 -0.000890980 3.186836665
>
> ATOMIC_POSITIONS (crystal)
> Y -0.082663260 -0.082649328 0.000174846
> Y -0.082874577 -0.082880779 0.499813069
> Mn1 0.205305857 0.560245090 0.249986248
> Mn2 0.561457263 0.207812583 0.750012603
> O -0.027464907 -0.056329933 0.249979343
> O -0.058482388 -0.031289340 0.750035308
> O 0.247951998 0.582250730 0.083800805
> O 0.582338704 0.248188482 0.916046842
> O 0.582135578 0.247948396 0.583964864
> O 0.247742343 0.582041831 0.416186070
> K_POINTS automatic
> 3 3 2 1 1 1
>
> Thanks alot in advance
> Tamene Regassa
> AAU
> Ethiopia
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
--
o ------------------------------------------------ o
| Gabriele Sclauzero, PhD Student |
| c/o: SISSA & CNR-INFM Democritos, |
| via Beirut 2-4, 34014 Trieste (Italy) |
| email: sclauzer at sissa.it |
| phone: +39 040 3787 511 |
| skype: gurlonotturno |
o ------------------------------------------------ o
More information about the Pw_forum
mailing list