QuantumATK W-2024.09 version released on Sep 9, 2024
0 Members and 1 Guest are viewing this topic.
voltagelist=[0.2*i for i in range(26)] for myvolt in voltagelist: calculator = DeviceLCAOCalculator( basis_set=basis_set, numerical_accuracy_parameters=device_numerical_accuracy_parameters, electrode_calculators= [left_electrode_calculator, right_electrode_calculator], electrode_voltages=( 0.0*Volt, myvolt*Volt) ) device_configuration.setCalculator(calculator) nlprint(device_configuration) device_configuration.update() if(myvolt<0.3) nlsave('C80sixmr.nc', device_configuration) # ------------------------------------------------------------- # Transmission spectrum # ------------------------------------------------------------- transmission_spectrum = TransmissionSpectrum( configuration=device_configuration, energies=numpy.linspace(-5,5,100)*eV, kpoints=MonkhorstPackGrid(1,1), energy_zero_parameter=AverageFermiLevel, infinitesimal=1e-06*eV, self_energy_calculator=KrylovSelfEnergy(), ) nlsave('C80sixmr.nc', transmission_spectrum) nlprint(transmission_spectrum)