|
|
【答案】应助回帖
感谢参与,应助指数 +1
在device calculator后面加上
for bias in numpy.linspace(-0.2,0.2,3)*Volt: # this is added by guohuazhong
# device_configuration.setCalculator(calculator) ## this line should be revised for loop calculations
device_configuration.setCalculator(
calculator(electrode_voltages=(0.5*bias,-0.5*bias)),
initial_state=device_configuration
)
nlprint(device_configuration)
device_configuration.update()
# nlsave('C:/Documents and Settings/ghz/VNL/lih2li_iv_scf_%g.nc' % bias.inUnitsOf(Volt), device_configuration) # write every scf step separately
nlsave('C:/Documents and Settings/ghz/VNL/ghz-device-transmission.nc', device_configuration) # write all the scf steps in one file
# ... calculate and save the transmission spectrum,
transmission_spectrum = TransmissionSpectrum(
configuration=device_configuration,
energies=numpy.linspace(-5,5,401)*eV,
kpoints=MonkhorstPackGrid(3,3)
)
nlsave('C:/Documents and Settings/ghz/VNL/ghz-transmission.nc',transmission_spectrum)
里面的数值自己修改就可以了 |
|