Date: Fri, 29 Mar 2024 00:44:39 -0700 (PDT) Message-ID: <1612110915.8316.1711698279976@co.atl-prod-01> Subject: Exported From Confluence MIME-Version: 1.0 Content-Type: multipart/related; boundary="----=_Part_8315_582233836.1711698279976" ------=_Part_8315_582233836.1711698279976 Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable Content-Location: file:///C:/exported.html
VASP is restricted to licensed users. It may not be used in cont= ractual research in cooperation with industry or military.
vasp-ve= r-bsub <version> myscript.py vasp-ver-bsub-native <version> -q suncat-test -n 2 -o junk.log
See the bottom of this page for available versions.
The first is the ASE style. The second is native VASP style, and assumes= all input files (POTCAR, INCAR, etc.) are created by the user are in the c= urrent directory. POTCARs can be assembled from the directory /nfs/slac/g/s= uncatfs/sw/vasp/pseudo52.
Thanks to Joel Varley for providing this BEEF example.
#!/usr/= bin/env python #LSF -q suncat-test #LSF -n 2 #LSF -o std.out #LSF -e std.err from ase import * from ase import io from ase.lattice.spacegroup import crystal from ase.io.trajectory import PickleTrajectory from ase.optimize import BFGS, QuasiNewton from ase.units import kJ, Hartree from ase.utils.eos import EquationOfState from ase.calculators.vasp import Vasp from ase.data.molecules import molecule as molecule from ase.vibrations import Vibrations import os, sys ###########################################################################= ## ###########################################################################= ## name =3D 'H' a =3D 5.0 atoms =3D Atoms([Atom('H',(a/2, a/2, a/2))], cell=3D(a, a, a)) # note that BEEF is only available in suncat vasp version 2 (v2) or later. xcf =3D 'BEEF' encut =3D 400 nkp =3D 1 sigma =3D 0.05 ibrion=3D -1 # no atomic displacement done by VASP (only ASE) ispin =3D 1 # NOT spin-polarized (ispin=3D2 for SP) # some types of GGA vaspgga =3D {'RPBE' : 'RP', 'PBE' : 'PE', 'AM05' : 'AM', 'PBEsol' : 'PS', 'PW86' : 'PW', 'PW91' : '91', 'BEEF' : 'BF', 'GGA' : '--', } if xcf=3D=3D'BEEF': luse_vdw=3DTrue if not (os.path.exists('vdw_kernel.bindat')): os.symlink('/nfs/slac/g/suncatfs/sw/vasp/vdw_kernel.bindat','vdw_ke= rnel.bindat') else: luse_vdw=3DFalse # VASP uses LDA or PBE PAW pseudopotentials as a start (gga flag modifies t= he GGA paramaterization if xcf in vaspgga.keys(): xc =3D 'PBE' else: xc =3D 'LDA' vaspgga['LDA'] =3D '--' ###########################################################################= ## ###########################################################################= ## calc =3D Vasp(encut=3Dencut, xc=3Dxc, gga=3Dvaspgga[xcf], kpts =3D (nkp,nkp,nkp), gamma =3D False, # Gamma-centered (defaults to Monkhorst-Pack) ismear=3D0, sigma=3Dsigma, ibrion=3Dibrion, nsw=3D0, # don't use the VASP internal relaxation, only use ASE #potim=3Dpotim, # use to take smaller ionic steps (or for vibra= tion/normal mode calculations) lvtot=3DFalse, npar=3D1, # use this if you run on one node (most calculations)= . see suncat confluence page for optimal setting lwave=3DFalse, luse_vdw=3Dluse_vdw, ispin=3Dispin, zab_vdw=3D-1.8867, # this selects vdw-DF2 #lreal=3DTrue, # can speed up if supercells are large ) atoms.set_calculator(calc) dyn =3D QuasiNewton(atoms, logfile=3Dname+'.log', trajectory=3Dname+'.traj'= ) dyn.run(fmax=3D0.05) electronicenergy =3D atoms.get_potential_energy() f =3D open('out.energy','w') f.write(str(electronicenergy)) f.close() # clean up big leftover files (WAVECAR or CHGCAR could be of use for some c= alculations) os.system('rm CHG vasprun.xml PCDAT WAVECAR CHGCAR') ###########################################################################= ##
Version |
Date |
Comment |
1 |
12/5/2012 |
initial version |
2,2a |
1/31/2013 |
add BEEF support courtesy of Johannes Voss ("= a" for suncat3 farm) |
3 |
1/31/2013 |
like v2, but change pseudo potentials to the = "52" versions. no suncat3 version yet |
4 |
2/28/2013 |
like v3, but add VTST support |
5 |
3/15/2013 |
includes BEEF, VTST, and 30% performance impr= ovement from using fftmpiw.F |
6 |
7/9/2013 |
copy of v5 but without -DNGZhalf flag. reques= ted by junyan for a non-collinear calculation. not recommended for general = use |
This information comes from Yian Zhu, who ran several 3x3x4 spin-polariz= ed Nickel calculations with the PBE functional.
=