The lengths of scattering region and spin polarization.

Asked by Jianjun Mao

Dear developers,

1. Recently, I have calculated several examples of IV curves using transiesta. I find that different length scattering regions make a large difference to the currents values with the same electrodes. How do I konw my scattering region is long enough or reasonable. For example, if I add some gas molecues on the scattering region to see the variation of currents before and after their adsorption.

2. What is the relaitonship between the total I and Spin ploarizaiton I. I= I( spinup)+ I (spindown). ? My result show that they are not equal or the average of I (spinup) and I (spindown). Is there anything wrong with my calculaitons? No error throughout the whole calculation.

Many thanks~!

Question information

Language:
English Edit question
Status:
Solved
For:
Siesta Edit question
Assignee:
No assignee Edit question
Solved by:
Nick Papior
Solved:
Last query:
Last reply:

This question was reopened

Revision history for this message
Nick Papior (nickpapior) said :
#1

1. This suggests to me that you haven't screened off your electrodes sufficiently. TranSiesta is based on *bulk* electrodes, and thus the electrostatics needs to be bulk-like at the electrode region of the scattering device. I.e. I(V) should not be length dependent if you simply add more electrode layers.

2. The total current is I = I(spin-up) + I(spin-down). In general I(spin-up) need not be equal to I(spin-down).

Revision history for this message
Jianjun Mao (jjmao) said :
#2

Many thanks for you quick answer, Nick Papior~

1. I will make more test calculaions.

2. When I open the spin polarizaions, the I(spinup)+ I(spindown) is not equal to the I (total) without spin polarization.

For example, same structure;

Under bias voltage: 1.5V
no spin polarization: I (total) = 0.8062E-05

Currents (ensure entire Fermi function window):
Left -> Right, V [V] / I [A]: 1.50000 V / 0.806203E-05 A
Left -> Right, V [V] / P [W]: 1.50000 V / -.602797E-05 W

open spin polarization: I (spin up)= 0.68E-05A I (spin down)=0.73E-05

Currents (ensure entire Fermi function window):
Left -> Right, V [V] / I [A]: 1.50000 V / 0.680710E-05 A
Left -> Right, V [V] / P [W]: 1.50000 V / -.519714E-05 W

Currents (ensure entire Fermi function window):
Left -> Right, V [V] / I [A]: 1.50000 V / 0.732193E-05 A
Left -> Right, V [V] / P [W]: 1.50000 V / -.547095E-05 W

I wonder that is there anything wrong with my calculations?

Revision history for this message
Nick Papior (nickpapior) said :
#3

1. Good.

2. Of course, it is a different Hamiltonian, the electronic structure for the non-polarized is not the average up/down. Secondly, TBtrans writes out the current for non-polarized as though it was one channel, i.e. you need to multiply by two for non-polarized calculations.

Revision history for this message
Jianjun Mao (jjmao) said :
#4

Thanks Nick Papior, that solved my question.

Revision history for this message
Jianjun Mao (jjmao) said :
#5

Dear Nick,

Does the output current I of the non-polarized calculation of the magnetic systems meaningful?( In my test, the output current cof a magnetic system with non-polarized calulation is not equal spinup plus spindown or ethier of them) Or a spin-ploarized calculation of magnetic system is a must and the Total current = I (spin up) + I (spin down) ?

Many thanks for your previous detailed answers~~

Revision history for this message
Jianjun Mao (jjmao) said :
#6

Dear Nick,

Does the output current I of the non-polarized calculation of the magnetic systems meaningful?( In my test, the output current cof a magnetic system with non-polarized calulation is not equal spinup plus spindown or ethier of them) Or a spin-ploarized calculation of magnetic system is a must and the Total current = I (spin up) + I (spin down) ?

Many thanks for your previous detailed answers~~

Revision history for this message
Nick Papior (nickpapior) said :
#7

Please see my answer 2 in my last comment.

They are not equal.

Revision history for this message
Jianjun Mao (jjmao) said :
#8

Dear Nick,

Thank you for your detaied answer~

Here, I still have this doubt: If a spin-ploarized calculation of magnetic system is a must? Or, the right way to get the total currents of a magnetic system under bias voltages?

Thanks you again~!

Revision history for this message
Best Nick Papior (nickpapior) said :
#9

Of course, if it is a polarized system you need to do a polarized calculation.

Revision history for this message
Jianjun Mao (jjmao) said :
#10

Thanks Nick Papior, that solved my question.