Attachment 'STM_VMD.py'
Download 1 #VMD
2 from Dacapo.ElectronicStates import ElectronicStates
3 from Dacapo import Dacapo
4 from ASE.Utilities.STMTool import STMTool
5 from ASE.Visualization.VMD import VMD
6 from ASE.IO.Cube import WriteCube
7
8 electronicstates = ElectronicStates(filename='Al-bcc_100_1layer.nc')
9 atoms=Dacapo.ReadAtoms('Al-bcc_100_1layer.nc')
10 n=0.1
11 stm = STMTool(electronicstates,contourvalue=1.4e-4,smoothfactor=n)
12 stmgrid = stm.GetStrippedGrid()
13 atoms.SetUnitCell(stmgrid.GetSpace().GetBasis(),fix=True)
14 WriteCube(atoms,stmgrid.Array,'test.cube')
Attached Files
To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.You are not allowed to attach a file to this page.