[gpaw-users] AssertionError in a TS09 vdW correction calculation
Alfredo Palace Carvalho
ajpalace at uevora.pt
Thu Aug 28 14:50:53 CEST 2014
Hi,
With the current stable version of gpaw (0.10.0.11364) running a
calculation with van der Waals corrections using the Tkachenko and
Scheffler correction (with the script below) produces an AssertionError
exception after the SCF cycles have converged (see the txtfile attached
as well as the log of the python error messages). With a previous 0.9
version (specifically, the gpaw-0.9.1.11146) I don't have this problem.
Here's a simple script that exhibits the problem:
===
import sys
import ase.io
from ase.parallel import parprint
from gpaw import GPAW, Mixer, FermiDirac
from ase.calculators.vdwcorrection import vdWTkatchenko09prl
from gpaw.analyse.hirshfeld import HirshfeldDensity,
HirshfeldPartitioning
from gpaw.analyse.vdwradii import vdWradii
name = sys.argv[1]
system = ase.io.read(name + '.CONTCAR')
system.set_pbc((True, True, True))
c = GPAW(xc = 'PBE',
occupations = FermiDirac(0.1),
txt = name + '.out',
mixer = Mixer(beta = 0.2, nmaxold = 5, weight = 100.0),
eigensolver = 'dav')
cc = vdWTkatchenko09prl(HirshfeldPartitioning(c),
vdWradii(system.get_chemical_symbols(), 'PBE'))
system.set_calculator(cc)
en = system.get_potential_energy()
parprint("Energy: " + repr(en))
===
Thanks.
Best regards,
Alfredo Carvalho
More information about the gpaw-users
mailing list