SEMITIP, VERSION 1, Electrostatic Computations for a Probe Tip near a Semiconductor

R. M. Feenstra
Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213

This program computes the electrostatic potential resulting from a metallic probe tip near a semiconductor, as described in the paper J. Vac. Sci. Technol. B 21, 2080 (2003). Click here for preprint of paper, in pdf format. Additional documentation for the program is available in the SEMITIP V1 Technical Documentation.

A compiled version of the code, which should run on any Windows PC, is available in the file semitip_v1.exe. Input for the executable code comes from the file FORT.9. Download these two files, into filenames "semitip_v1.exe" and "fort.9". Then, run the code just by double clicking on it. Using a text editor, the input parameters in FORT.9 can be changed to whatever values are desired. In addition to the parameter values, this file also contains comments as to the meaning of each parameter.

Output from the program is contained in the following files (output depends on the value of the output parameter IWRIT as specified in the input file FORT.9):

• FORT.10 - gives the numerical results for the following quantities:
  • tip radius of curvature (nm)
  • tip-sample separation (nm)
  • tip potential (eV)
  • opening angle of tip shank
  • PHI0 - the surface potential at a point directly opposite the tip apex (eV)
  • ratio of PHI0 to PHI00 where PHI00 is a 1-D value for PHI0
• FORT.11 - provides the potential (eV) along the central axis, as a function of z-distance (output for IWRIT>=1) (the third column gives auxiliary output which are the potential values for all grid points at the maximal radial value in the problem, although the z-values for these points are not output and in particular do not correspond to the values in column 1).
• FORT.12 - provides the potential (eV) along the surface, as a function of the radial distance from the central axis (output for IWRIT>=1)
• FORT.13 - gives the entire array of potential values (eV) (output for IWRIT>=3)
• FORT.16 - gives an exact copy of the output to the console
• FORT.20, FORT.21, ... - contour lines of the potential (nm) (output for IWRIT>=2)

All of the parameters in the program can be varied using the input file FORT.9, with the exception of the array sizes (NRDIM, NVDIM, NSDIM, and NE) and the specification of an auxiliary function other than a hemisphere which defines a protrusion on the end of the probe tip. Modification of those parameters can be accomplished by changing the source code of the program. The source code is also available, in the files:

• SEMITIP_V1.F - main routine.
• SEMITIP2.F - performs the detailed finite element solution of Poisson's equation.
• SEMIRHO.F - auxiliary routines for computing semiconductor charge densities.
• CONTR1.F - plots contour lines of the computed potential.