### Author Topic: Current at each bias  (Read 454 times)

0 Members and 1 Guest are viewing this topic.

#### Jahanzaib

• Heavy QuantumATK user
• Posts: 87
• Country:
• Reputation: 2
##### Current at each bias
« on: September 2, 2024, 15:44 »
Dear Expert,
I hope you are doing well.

I am interested in calculating the current at each bias given in code below. This code work well, but I don't have an idea how to calculate current at each bias, any suggestion?

pos_bias_list = [0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70]*Volt

zero_bias_file = '/home/hartree/jem/wolf6252/N-P/Device/Fe-H-D/Fe-H-D-3-3/Fe-H-D-3-3.hdf5'
device_configuration = nlread(zero_bias_file, DeviceConfiguration)[0]

for bias in pos_bias_list:
if processIsMaster():
print("Bias is now: ", bias)
# Get the calculator
calculator = device_configuration.calculator()

# Set the bias voltage
calculator=calculator(electrode_voltages=(bias/2, -bias/2))

# Attach the calculator and use the old initial state
device_configuration.setCalculator(
calculator(),
initial_state=device_configuration)
device_configuration.update()
nlsave('device_bias_%.2f.hdf5' % bias.inUnitsOf(Volt), device_configuration)
nlprint(device_configuration)

#### Anders Blom

• QuantumATK Staff
• Supreme QuantumATK Wizard
• Posts: 5519
• Country:
• Reputation: 89
##### Re: Current at each bias
« Reply #1 on: September 5, 2024, 21:35 »
To get the current you also have to compute the transmission spectrum at each bias (can be done afterwards, using the saved configurations for each bias), and from each transmission spectrum you get the current.

#### Jahanzaib

• Heavy QuantumATK user
• Posts: 87
• Country:
• Reputation: 2
##### Re: Current at each bias
« Reply #2 on: Yesterday at 01:08 »
If I understand you correctly then it's mean I have to apply different bias let's say -0.1, 0.2 etc and calculate the TS. ok?

I calculated the TS with -0.1 and also got current as shown in the figure? this is what you are explaining?