Hi Camilla, > I have a question regarding the lead_info[x].wave_functions[x]. I am not sure > what this is actually returning? I tried to run the function on your > quantum_wire_revisited.py-system, like this > > prop = sm.lead_info[0] > print(prop.wave_functions[0]) > > with W = 2, if that is important. It then returns a 1x2 array or 2x2 matrix > depending on the energy. What does this matrix contain? I read about > propagatingmodes in the documentation, but I am not sure what it meant. > > Best, > Camilla
`lead_info[i].wave_functions` contains the wavefunctions of the propagating modes in lead `i`. This is a 2D array, where the first axis corresponds to the orbitals in a unit cell of the lead, and the second axis is the mode number. These wave functions are essentially the unit-cell-periodic part of the Bloch states in the lead. It makes sense that the length of the second axis can be different depending on the energy, as at different energies there may be different numbers of propagating modes open in the lead. Hope that helps, Joe
signature.asc
Description: PGP signature