feat(PhyslibAlpha): current eigenstates of the open tight binding chain - #1708
Conversation
Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
|
Thank you for this pull-request (PR). If this is your first PR, welcome to the community! Below is what will happen next. Please read carefully if you are not familiar with the process. You may open other PRs while this one is being reviewed, and can stack PRs on top of each other, so don't let these steps slow you down.
Tip: The easiest way to get have a fast review is to submit a PR that is small and self-contained, and has clear documentation explaining why things are the way they are in your chages. If you have any problems or questions, please reach out to the community on the Zulip. |
jstoobysmith
left a comment
There was a problem hiding this comment.
Approved. Looks good to me.
hello guys, Joseph Timeroot, im trying to make it faster as i can. what this PR exactly does do?
this PR comes directly from the #1696, where we defines J current from the open chain, and here we got the auto states on himself, it's only one new PR with 180 lines, easy to review i don't want a be a charge for the proyect, so im keepeing short PR with personal Review, awaiting, any way you has the final word here, so im very very glad. tks alot. so, let me be short 4 dots.
1.-currentEingeinstate k: the stationary wave from the open chain , sin((n+1) k * pi/(N+1)), wich it's exactly 0 zero just outside from both ends product fase i^n on each site, then the state remains "psy_k = sumatory_n * i^n * sin ((n+1)k * pi/(N+1)) |n).
2.-current_currentEingeinstate: with phase (psy_k) it's itself current autovalue: J (psy_k) = 2 a * t * cos (k * pi /(N+1)) * psy_k, every single one psy_k has an exact current and his parameter comes directly from "currentEingueinvalue k"
3.-norm_sq_currentEingeinstate. for 1 ≤ k ≤ N the norm it's ||psy_k||² = (N+1)/2, this leaves the road ready to normalize the next one PR's
why and what this PR does itself, this is the floor max current where ( k=1, (current 2 * a * t ) * cos ( pi/ (N+1))), the PR 10 (maxcurrentState) and in subsecuent PR's the 18 PR give us the speed limit |( J )| ≤ 2 | a * t | cos ( pi / (N+1))
reference thigthBindings.
iI hope this commit doesn't take too much of your time. Thanks for everything, guys! Im here waithing for your review and will be ready for any modifications. Thanks once again, timeroot and joseph¡