# About Orbital Projected Transmission

Hello Nick,

I am trying to calculate the transmission of a 2D lateral heterojunction. I am wondering if we can calculate the transmission projected on different orbitals?

If I have bulk electrode as every part of my system, using the output SystemLabel.

If I have a heterojunction, can I get the probability of electron transmission from the each channel on one side to each channel on the other side? Like from px orbital on the Left Electrode to py orbital on the Right Electrode? Like the T_i_j and R_i_j model in

M. Büttiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. Rev. B 31, 6207.

I guess these could be done using TBT.Projs section in TBtrans. But I I am having a hard time understanding this part of TBtrans manual. Is there any recipe or example fdf file for that? I checked ts-tbt-

Thank you very much Nick!

Best,

Yuefei

## Question information

- Language:
- English Edit question

- Status:
- Answered

- For:
- Siesta Edit question

- Assignee:
- No assignee Edit question

- Last query:
- 2019-09-09

- Last reply:
- 2019-09-09

Nick Papior (nickpapior) said : | #1 |

1) "If I have bulk electrode". No, it does not.

2) "If I have a heterojunction". The discussion in the referenced article discusses "channels" which are not necessarily adhering to a single px orbital. In practice a "channel" is equivalent to a propagating bulk eigenstate which is then transferred across the junction. I.e. it does not belong to *only* px orbitals per see.

3) "I guess" yes and no. The projection section (TBT.Projs) projects the scattering matrices onto specific eigenstates in a predefined device region. As such one can figure out how much of the transmission goes through a specific eigenstate of say a molecule in the junction. In this sense it could be used like the article you referenced since one can figure out if a particular ring eigenstate carries 100% of the current. However, the projection section implemented in TBtrans works best in the limit of weak electrode coupling to the projection region. I.e. weakly coupled molecules.

Lastly, if you want the local transmission probability ("bond-currents") you may simply do: TBT.Current.Orb true

which tells tbtrans to calculate the local transmissions going from one orbital to another. However this will only give you local information and not whether an electron comes from the px of a particular atom in the left lead and the probability it arrives at a py orbital of another atom in the right lead. The bond-currents can however tell you how the current flows in your device and this is probably what you are looking for.

Yuefei Huang (yh46) said : | #2 |

Thank you Nick!

1) For bulk electrode, I have done the following to extract orbital transmission. By pure SIESTA calculation of bulk electrode, I got eigenstate energy in the whole brillouin zone and the projection of each state on different orbitals. Then by counting number of band crossing each energy and transverse k and assign that G0 conductance to each orbital according to the their weight. I can get the transmission of each orbital. But this only works for bulk electrode.

2),Can I get the probability of a left bulk eigenstate phi_L(n,ky,E) jumping into a right bulk eigenstate phi_R(n',ky,E), n and n' is band index, x is vacuum direction and z is transport? It would be great. Combined with orbital projected DOS of these bulk eigenstates, should be valuable.

3) Does 'weakly coupled molecules' mean that VL and VR are non-zero in only atoms close to the Electrode? Does it work only for molecule, i.e. non-periodic in the transverse direction?

About the projected transmission, I did a test using graphene. Following is the part of .fdf file that is related to projection.

###

%block TBT.Projs

M-L

%endblock

%block TBT.Proj.M-L

atom from 17 plus 16

Gamma .false.

proj L-5

level -5

end

proj L-4

level -4

end

proj L-3

level -3

end

proj L-2

level -2

end

proj L-1

level -1

end

proj L1

level 1

end

proj L2

level 2

end

proj L3

level 3

end

proj L4

level 4

end

proj L5

level 5

end

%endblock

%block TBT.Projs.T

from Left.M-L to

Right

end

%endblock

TBT.T.Eig 10

TBT.T.Out .true.

TBT.Current.Orb .true.

TBT.Atoms.Device [ 17 -- 72]

TBT.Projs.

TBT.Projs.T.Out .true.

###

4.)From Dataset[

5) From Dataset[

6) From your lecture video, https:/

Thank you very much! You are very kind and helpful!

Best,

Yuefei

Nick Papior (nickpapior) said : | #3 |

1) I only think this makes sense for a bulk calculation.

2) Yes. You can add another layer of the electrode to the device region for each of the electrodes and define a projection region on those layers. Then you can do two projections (one per electrode). Then simply follow what you do in 3 to get a projected transmission from eigenstates to other eigenstates.

Note that here you basically need to do some manual labor by pre-creating the Projection file (TBT.Projs.Init true), then edit the file to contain the correct eigenstates from the electrodes.

3) No it means that VL/VR are *small*. And no, you can also do projections for k-points.

4) Restate + 1j Imstate is the projection eigenstate (i.e. S^(1/2) | psi >). To only get |psi>, you need to do "TBT.Proj.

5) You need to use sisl to get that information out. See here: http://

6) Currently you can use Inelastica to get the real-space representation of the eigenchannels.

Yuefei Huang (yh46) said : | #4 |

Hello Nick,

Thank you for answering!

In your tbtrans manual, I guess there is some formatting problem with Eqn 13~17, I think they should be similar to Eqn 35 in paper "Improvements on non-equilibrium and transport Green function techniques: The next-generation transiesta". But ther might be some typo in the manual.

I have compiled Inelastica and it works well, thank you! I got some new questions these days.

1. About plotting eigenchannel, I am asked to provide the SIESTA numbering of Device atoms. If I do not specify this, it plots the eigenchannel wavefunction on all atoms including the two electrodes. What is the meaning of this wave function if it is plotted on all atoms?

2. As I mentioned in previously I am doing a lateral junction of 2D material. The eigenchannel wavefunction plotted for all atoms, looks symmetric, as if the two ends of the Left-Device-Right system are connected periodically. Is it because the Hamiltonian and DM are solved periodically in transport direction? Is it possible to get the eigenchannel wavefunction like propagating bloch wave of Left + extension of left bloch wave in left side of junction in Device + transmitted wave on right side of junction in Device + transmitted wave in Right?

3. I am confused about relation between NEGF transmission and Landauer Buttiker picture. In NEGF Tr[Γ_L G Γ_R G] is used, while in Landauer it is Tr[tt], where t is the transmission matrix, hermitian sign not typed. In the paper 'Density-functional method for nonequilibrium electron transport' it is said t=sqrt(Γ)Gsqrt(Γ). My question is how these two matrices are related? t has a dimension of n_L X n_R, where n_L /n_R is the number of bloch states in left/right electrode. But sqrt(Γ)Gsqrt(Γ) has the dimension of n_D X n_D, n_D is the number of orbitals in Device region. They have different dimensions. How do understand their equality?

Thank you very much!

Nick Papior (nickpapior) said : | #5 |

Hi,

0. Thanks for catching this, indeed this was a typo! :)

1. Please direct any Inelastica questions to the Inelastica developers. :)

2. See 1. I don't fully know your system.

3. All equations are the same. All matrices has dimension n_D, however, Gamma matrices are just zero a lot of places and are thus typically only retained in a dimension n_L/n_R.

## Can you help with this problem?

Provide an answer of your own, or ask Yuefei Huang for more information if necessary.