[gpaw-users] DFT+U

Chuck fus dansexee at gmail.com
Sun Oct 13 01:49:59 CEST 2013


Hi guys,
I'm trying to run a DFT+U calculation following this link
https://wiki.fysik.dtu.dk/gpaw/epydoc/gpaw.test.Hubbard_U-pysrc.html
and I encountered a couple of issues. In particular finding the +U
setup for the element with the U parameter, saying 'Could not find
10:d,6.0.PBE-setup for "Ni'
The input file is below. Very best regards and thanks in advance
Danny

a = 4.19 # Lattice constants
b=a/np.sqrt(2)
mgo = Atoms(symbols='Ni2O2',
               pbc=True,
               cell=(b, b, a),
               positions=[(0, 0, 0),
                           (b/2, b/2, a/2),
                           (0, 0, a/2),
                           (b/2, b/2, 0)],
                magmoms=(2,-2,0,0))

calc = GPAW(nbands= -20,
      h=0.2,
      occupations=FermiDirac(width=0.05),
      poissonsolver=PoissonSolver(nn='M', relax='J'),
      setups={'Ni': '10:d,6.0'},
      convergence={'eigenstates':8e-4,'density': 1.0e-2,'energy': 0.1},
      kpts=(2, 2, 2),
      xc='PBE')
mgo.set_calculator(calc)
energy =mgo.get_potential_energy()


More information about the gpaw-users mailing list