qwwad_ef_plot

The qwwad_ef_plot program translates wavefunction and potential data files into a single plottable data file, with "prettier" formatting.

Input files

E*.r    # Energies of states
        # Column 1: state index
        # Column 2: energy [meV]

wf_*i.r # Wavefunction of state i
        # Column 1: position [m]
        # Column 2: amplitude [m^{-0.5}]

In each case, the '*' is replaced by 'e', 'l' or 'h' for electrons, light-holes or heavy-holes respectively.

Output files

vwf.r   # Plottable date file.
        # Column 1: position [angstrom].
        # Column 2: plottable function [meV].

The plottable functions in this file are contained in distinct data sets, each being separated by a blank line. The first set contains the confining potential in meV. Subsequent sets contain either the probability density (default) or wavefunction amplitude (if --plot-wf is specified). These sets are placed in ascending order, such that state |1> data is in data set 2, state |2> in set 3, state |3> in set 4 and so on.

The probability density (or wavefunction) plots are positioned on the y-axis such that they lie at the energy of the associated state. They are automatically scaled such that they fit neatly on the plot.

Options

 --wffileprefix arg (=wf_e) Prefix of wavefunction filenames.

 --wffileext arg (=.r)
        File extension for wavefunction files.

 --energyfile arg (=Ee.r)
        Filename of energy file.

--totalpotentialfile arg (=v.r) Name of file containing the total confining
        potential.

 --plotfile arg (=vwf.r)
        Name of file to which plottable data will be written.

 --nstmax arg (=10)
        Maximum number of states to plot.

 --style (=pd)
        Style of plot: 'pd' = probability density, 'wf' = wave functions.

Examples

Generate plottable data for all states in the input files:

qwwad_ef_plot

Generate plottable data, using wavefunctions instead of probability densities:

qwwad_ef_plot --style wf

Generate plottable data, using only the first three states in the input files:

qwwad_ef_plot --nstmax 3