Negative charge Nd ions within DFT+U calculation |
- Date: 2021/01/13 08:22
- Name: Aleksey
- Dear developers and users of OpenMX!
I need your help. I started study NdFeO3 perovskite with standart DFT within open-core and full pseupotentials for Nd. In this cases I obtained correct results. Next step is study this system within Hubbard model. My input see below (with links on out-file and stdout flow file). Hubbard's values uses from literature. Is problem in my input or is this possible bug in OpenMX?
System.CurrrentDirectory ./ System.Name NdFeO3 level.of.stdout 1 # default=1 (1-3) level.of.fileout 0 # default=1 (1-3) DATA.PATH /home/nano/.OpenMX/DFT_DATA13 Species.Number 3 <Definition.of.Atomic.Species Fe Fe6.0S-s2p2d2f1 Fe_PBE13S Nd Nd8.0-s2p2d2f2 Nd_PBE13 O O5.0-s2p2d1 O_PBE13 Definition.of.Atomic.Species> <Hubbard.U.Values Fe 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 5.8 2d 5.8 1f 0.0 Nd 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 0.0 2d 0.0 1f 6.3 2f 6.3 O 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 0.0 Hubbard.U.Values> Atoms.UnitVectors.Unit ANG <Atoms.UnitVectors 5.58870000000000 0.00000000000000 0.00000000000000 0.00000000000000 7.76190000000000 0.00000000000000 0.00000000000000 0.00000000000000 5.44890000000000 Atoms.UnitVectors> Atoms.Number 20 Atoms.SpeciesAndCoordinates.Unit FRAC <Atoms.SpeciesAndCoordinates 1 Fe 0.04907000000000 0.25000000000000 0.98912000000000 8.0 6.0 ON 2 Fe 0.95093000000000 0.75000000000000 0.01088000000000 8.0 6.0 ON 3 Fe 0.45093000000000 0.75000000000000 0.48912000000000 8.0 6.0 ON 4 Fe 0.54907000000000 0.25000000000000 0.51088000000000 8.0 6.0 ON 5 Nd 0.00000000000000 0.00000000000000 0.50000000000000 8.0 6.0 ON 6 Nd 0.50000000000000 0.00000000000000 0.00000000000000 8.0 6.0 ON 7 Nd 0.00000000000000 0.50000000000000 0.50000000000000 8.0 6.0 ON 8 Nd 0.50000000000000 0.50000000000000 0.00000000000000 8.0 6.0 ON 9 O 0.47600000000000 0.25000000000000 0.08660000000000 3.0 3.0 ON 10 O 0.52400000000000 0.75000000000000 0.91340000000000 3.0 3.0 ON 11 O 0.02400000000000 0.75000000000000 0.58660000000000 3.0 3.0 ON 12 O 0.97600000000000 0.25000000000000 0.41340000000000 3.0 3.0 ON 13 O 0.70510000000000 0.95410000000000 0.29550000000000 3.0 3.0 ON 14 O 0.29490000000000 0.04590000000000 0.70450000000000 3.0 3.0 ON 15 O 0.79490000000000 0.04590000000000 0.79550000000000 3.0 3.0 ON 16 O 0.20510000000000 0.95410000000000 0.20450000000000 3.0 3.0 ON 17 O 0.29490000000000 0.45410000000000 0.70450000000000 3.0 3.0 ON 18 O 0.70510000000000 0.54590000000000 0.29550000000000 3.0 3.0 ON 19 O 0.20510000000000 0.54590000000000 0.20450000000000 3.0 3.0 ON 20 O 0.79490000000000 0.45410000000000 0.79550000000000 3.0 3.0 ON Atoms.SpeciesAndCoordinates> scf.XcType LSDA-PW scf.SpinPolarization ON scf.Hubbard.U ON scf.Hubbard.Occupation DUAL scf.SpinOrbit.Coupling OFF scf.ElectronicTemperature 300.0 #scf.energycutoff 150 scf.Ngrid 42 57 40 scf.maxIter 200 scf.EigenvalueSolver BAND scf.Kgrid 5 3 5 scf.Generation.Kpoint REGULAR scf.Mixing.Type RMM-DIIS scf.lapack.dste DSTEVX scf.eigen.lib LAPACK scf.ProExpn.VNA ON scf.Init.Mixing.Weight 0.060 scf.Min.Mixing.Weight 0.030 scf.Max.Mixing.Weight 0.100 scf.Mixing.History 30 scf.Mixing.StartPulay 30 scf.Mixing.EveryPulay 1 #scf.Kerker.factor 10 scf.criterion 7.35E-7 scf.restart OFF
Links: Input -- https://yadi.sk/d/4QCch4r5VtG-9A Output -- https://yadi.sk/d/UkJnPNBv49C4-g Stdout -- https://yadi.sk/d/sbBTg_A6aM8f_g

| |