Top Page > Browsing

DFT+U convergence problem Date: 2023/03/01 09:53 Name: Vipin Kumar Dear developers and users, I am calculating the  the band structure and DOS of NiCl2 with DFT+U, it not converging even up to 500 steps. Could you suggest modification/changes in .dat file wherever I am missing some parameters to get the converged structure? The following is the .dat file # DATA.PATH  /home/msrc/OPEN/openmx3.9/DFT_DATA19 System.CurrrentDirectory        ./    # default=./ System.Name                    NiCl2-DFTU level.of.stdout                  1 level.of.fileout                  1 Species.Number                    2 <Definition.of.Atomic.Species   Cl    Cl7.0-s2p2d1      Cl_PBE19   Ni    Ni6.0S-s2p2d1    Ni_PBE19S Definition.of.Atomic.Species> Atoms.Number                      3 Atoms.SpeciesAndCoordinates.Unit  FRAC <Atoms.SpeciesAndCoordinates   1    Cl    0.6666210  0.3333790  0.0869380    3.5  3.5 off   2    Cl    0.3333790  0.6666210  0.9130620    3.5  3.5 off   3    Ni    0.0000000  0.0000000  0.0000000    8.0  8.0 off Atoms.SpeciesAndCoordinates> Atoms.UnitVectors.Unit            Ang <Atoms.UnitVectors     3.4400000    0.0000000    0.0000000     -1.7199990    2.9791280    0.0000000     0.0000000    0.0000000    15.0000000 Atoms.UnitVectors> scf.XcType                    GGA-PBE scf.SpinPolarization          on scf.ElectronicTemperature    700.0 scf.energycutoff              240.0 scf.maxIter                  500 scf.EigenvalueSolver          band scf.Kgrid                    9  9  1 scf.Mixing.Type              rmm-diisk scf.Init.Mixing.Weight        0.001 scf.Min.Mixing.Weight        0.001 scf.Max.Mixing.Weight        0.10 scf.Mixing.History            30 scf.Mixing.StartPulay        30 scf.criterion                1.0e-7 MD.Type                      nomd MD.maxIter                    1 MD.TimeStep                  1.0 MD.Opt.criterion              0.0003 Dos.fileout                  on Dos.Erange                    -10.0  10.0 Dos.Kgrid                    35  35  1 Band.dispersion    on Band.Nkpath  3 <Band.kpath 32  0.000000 0.000000 0.000000  0.333333 0.333333 0.000000 G K 18  0.333333 0.333333 0.000000  0.500000 0.000000 0.000000 K M 37  0.500000 0.000000 0.000000  0.000000 0.000000 0.000000 M G Band.kpath> scf.Hubbard.U on scf.DFTU.Type 2 scf.Hubbard.Occupation dual scf.dc.Type cFLL <Hubbard.U.values # eV Ni 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 6.0 Cl 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 0.0 Hubbard.U.values> <Hund.J.values # eV Ni 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 0.0 Cl 1s 0.0 2s 0.0 1p 0.0 2p 0.0 1d 0.0 Hund.J.values> Your help will be appreciated a lot. Regards, Vipin

Page: [1]

Re: DFT+U convergence problem ( No.1 ) Date: 2023/04/12 12:53 Name: T. Ozaki Hi, With the following choice for scf.Mixing.Type with no other change scf.Mixing.Type              rmm-diish I could get the convergent result as *********************************************************** ***********************************************************                   SCF history at MD= 1 *********************************************************** ***********************************************************   SCF=  1  NormRD=  1.000000000000  Uele= -25.481751508963   SCF=  2  NormRD=  0.409890448611  Uele= -25.479008383997   SCF=  3  NormRD=  0.407152314325  Uele= -25.275107399161   SCF=  4  NormRD=  0.212126780996  Uele= -25.173130037608   SCF=  5  NormRD=  0.131220891300  Uele= -25.128819600738   SCF=  6  NormRD=  0.102698989977  Uele= -25.106372034972   SCF=  7  NormRD=  0.088246038305  Uele= -25.092710732904   SCF=  8  NormRD=  0.078112430016  Uele= -25.082900931965   SCF=  9  NormRD=  0.069870797254  Uele= -25.075047601849   SCF=  10  NormRD=  0.062786969408  Uele= -25.068414514476   SCF=  11  NormRD=  0.056570466611  Uele= -25.062712536547   SCF=  12  NormRD=  0.051077049793  Uele= -25.057827367033   SCF=  13  NormRD=  0.046223612281  Uele= -25.053715484419   SCF=  14  NormRD=  0.041959001891  Uele= -25.050365084954   SCF=  15  NormRD=  0.038252749783  Uele= -25.047782172878   SCF=  16  NormRD=  0.035090199546  Uele= -25.045986260469   SCF=  17  NormRD=  0.032469505569  Uele= -25.045009526985   SCF=  18  NormRD=  0.030398073328  Uele= -25.044896910754   SCF=  19  NormRD=  0.028886893563  Uele= -25.045705979497   SCF=  20  NormRD=  0.027942252902  Uele= -25.047505910576   SCF=  21  NormRD=  0.027555985141  Uele= -25.050375067711   SCF=  22  NormRD=  0.027697114951  Uele= -25.052371196662   SCF=  23  NormRD=  0.027957514961  Uele= -25.053499658130   SCF=  24  NormRD=  0.028128197085  Uele= -25.054095284824   SCF=  25  NormRD=  0.028222584484  Uele= -25.054400789635   SCF=  26  NormRD=  0.028271899342  Uele= -25.054555444801   SCF=  27  NormRD=  0.028297069310  Uele= -25.054633245539   SCF=  28  NormRD=  0.028309780203  Uele= -25.054684559951   SCF=  29  NormRD=  0.028318183156  Uele= -25.054735992528   SCF=  30  NormRD=  0.026211237998  Uele= -24.940639936109   SCF=  31  NormRD=  0.064784337125  Uele= -24.969589682253   SCF=  32  NormRD=  0.033454547672  Uele= -24.985845867614   SCF=  33  NormRD=  0.018491712156  Uele= -24.993119186089   SCF=  34  NormRD=  0.012868640854  Uele= -24.996497360589   SCF=  35  NormRD=  0.010292431345  Uele= -24.963600133741   SCF=  36  NormRD=  0.039138634372  Uele= -24.965653473810   SCF=  37  NormRD=  0.036860248710  Uele= -24.975305588454   SCF=  38  NormRD=  0.028686196058  Uele= -24.997049261041   SCF=  39  NormRD=  0.009962302186  Uele= -25.005607763584   SCF=  40  NormRD=  0.002485583342  Uele= -25.006607170796   SCF=  41  NormRD=  0.001698104904  Uele= -25.006459888578   SCF=  42  NormRD=  0.001908140145  Uele= -25.006479126132   SCF=  43  NormRD=  0.001953425643  Uele= -25.008577823436   SCF=  44  NormRD=  0.000166123288  Uele= -25.008912384890   SCF=  45  NormRD=  0.000200355289  Uele= -25.008819384949   SCF=  46  NormRD=  0.000141595829  Uele= -25.008712083063   SCF=  47  NormRD=  0.000053657594  Uele= -25.008704880517   SCF=  48  NormRD=  0.000048293305  Uele= -25.008668286180   SCF=  49  NormRD=  0.000016531786  Uele= -25.008653908490   SCF=  50  NormRD=  0.000003896955  Uele= -25.008650483599   SCF=  51  NormRD=  0.000001097010  Uele= -25.008649259297   SCF=  52  NormRD=  0.000000419195  Uele= -25.008649647343   SCF=  53  NormRD=  0.000000479771  Uele= -25.008650205627   SCF=  54  NormRD=  0.000000974725  Uele= -25.008650051761   SCF=  55  NormRD=  0.000000856760  Uele= -25.008649879792   SCF=  56  NormRD=  0.000000714713  Uele= -25.008649741024   SCF=  57  NormRD=  0.000000599355  Uele= -25.008649627366   SCF=  58  NormRD=  0.000000504589  Uele= -25.008649507840   SCF=  59  NormRD=  0.000000403862  Uele= -25.008648742904   SCF=  60  NormRD=  0.000000266411  Uele= -25.008648765624 As for "rmm-diish", please take a look at https://www.openmx-square.org/openmx_man3.9/node40.html Regards, TO

Page: [1]