Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
from gpaw import GPAW, FermiDirac
from gpaw.lrtddft import LrTDDFT
h = 0.25
atoms = Cluster('relaxed.traj')
atoms.minimal_box(3.5, h=h)
# relax the molecule
calc = GPAW(h=h, occupations=FermiDirac(width=0.1),
eigensolver='cg', symmetry={'point_group': False},
nbands=10, convergence={'eigenstates':1.e-5,
'bands':4})
atoms.calc = calc
# use only the 4 converged states for linear response calculation
rr = ResonantRaman(atoms, LrTDDFT, exkwargs={'jend':3})
rr.run()