yes of course. You can:
- use the IVCurve analysis in the Scripter to do everything in one go.
You will have access to the two transmission spectra from the IVCurve object.
see for example http://quantumwise.com/publications/tutorials/item/820-ni-silicide-si-interfaces#h7-1-iv-curves
- run something like this:
device_configuration =...
calculator = ...
for bias in [0, 0.1, 0.2]*Volt:
device_configuration.setCalculator(calculator(electrode_voltages=(0,bias)),initial_state=initial_state)
nlprint(device_configuration)
device_configuration.update()
nlsave("filename.nc", device_configuration)
initial_state = device_configuration
# -------------------------------------------------------------
# Transmission spectrum
# -------------------------------------------------------------
transmission_spectrum = TransmissionSpectrum(
configuration=device_configuration,
energies=numpy.linspace(-2,2,101)*eV,
kpoints=MonkhorstPackGrid(21,21),
energy_zero_parameter=AverageFermiLevel,
infinitesimal=1e-06*eV,
self_energy_calculator=RecursionSelfEnergy(),
)
nlsave("filename.nc", transmission_spectrum)
nlprint(transmission_spectrum)