1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
# type: ignore
from ase.cluster.cubic import FaceCenteredCubic
from ase.calculators.turbomole import Turbomole
def test_turbomole_au13():
surfaces = [(1, 0, 0), (1, 1, 0), (1, 1, 1)]
layers = [1, 2, 1]
atoms = FaceCenteredCubic('Au', surfaces, layers, latticeconstant=4.08)
params = {
'title': 'Au13-',
'task': 'energy',
'basis set name': 'def2-SV(P)',
'total charge': -1,
'multiplicity': 1,
'use dft': True,
'density functional': 'pbe',
'use resolution of identity': True,
'ri memory': 1000,
'use fermi smearing': True,
'fermi initial temperature': 500,
'fermi final temperature': 100,
'fermi annealing factor': 0.9,
'fermi homo-lumo gap criterion': 0.09,
'fermi stopping criterion': 0.002,
'scf energy convergence': 1.e-4,
'scf iterations': 250
}
calc = Turbomole(**params)
atoms.calc = calc
calc.calculate(atoms)
# use the get_property() method
print(calc.get_property('energy'))
print(calc.get_property('dipole'))
# test restart
params = {
'task': 'gradient',
'scf energy convergence': 1.e-6
}
calc = Turbomole(restart=True, **params)
assert calc.converged
calc.calculate()
print(calc.get_property('energy'))
print(calc.get_property('forces'))
print(calc.get_property('dipole'))
|
