'ecut':ecut, 'xc':'PBE', 'pps':'fhi', 'nstep':700, 'diemac': 1000, 'width':w, 'kpts':kd, 'chksymbreak':0, 'ecutsm':0.05, 'nsym':1, } calcfactory = Factory(Abinit, 'Abinit', 'label', **calcopts) taskopts = {} task = Task( calcfactory=calcfactory, tag=tag, fit=fit, use_lock_files=True, **taskopts) if __name__ == '__main__': keys = task.collection.names for m in ['Mn']: keys.remove(m) # EOS fails: reason unknown # run just only system to check if scripts work #task.run(keys) task.run(['Si'])
fit = (5, 0.02) w = 0.06 ecut = 1200 kd = 8.0 tag = 'dcdft_%s_gpaw_pw' % xc.lower() task = Task( calcfactory=Factory(xc=xc, mode=PW(ecut), occupations=FermiDirac(w), maxiter=250, kptdensity=kd, ), tag=tag, fit=fit, use_lock_files=True, ) if __name__ == '__main__': if element is None: keys = set(parameters.keys()).intersection(set(task.collection.names)) for s in elements_slow: keys.remove(s) # those are slow, run separately else: keys = [element] task.run(keys)
return calculator calcopts = { "toldfe": 1.0e-7, "ecut": ecut, "xc": "PBE", "pps": "fhi", "nstep": 700, "diemac": 1000, "width": w, "kpts": kd, "chksymbreak": 0, "ecutsm": 0.05, "nsym": 1, } calcfactory = Factory(Abinit, "Abinit", "label", **calcopts) taskopts = {} task = Task(calcfactory=calcfactory, tag=tag, fit=fit, use_lock_files=True, **taskopts) if __name__ == "__main__": keys = task.collection.names for m in ["Mn"]: keys.remove(m) # EOS fails: reason unknown # run just only system to check if scripts work # task.run(keys) task.run(["Si"])
#'occupation':0.001, # 0.001 #'maxsteps':None, #'maxtime':None #}, ft=w, symmetry=True, spinpol=True, # must set explicitly calculate_stress=False, deletenc=True, # start fresh every time kptdensity=kd, ) taskopts = {} task = Task( calcfactory=calcfactory, tag=tag, fit=fit, use_lock_files=True, **taskopts) if __name__ == '__main__': keys = set(psp.keys()).intersection(set(task.collection.names)) for m in ['Be', 'Cr', 'Kr', 'Xe']: keys.remove(m) # EOS fails: maybe need higher ecut? for m in ['Zn', 'Zr']: keys.remove(m) # do not converge # run just only system to check if scripts work keys = ['Si'] task.run(keys)