Тёмный

Eigenstates of ANY 1D Potential in PYTHON 

Mr. P Solver
Подписаться 134 тыс.
Просмотров 35 тыс.
50% 1

Remember having to solve problems analytically? What a pain. With python you can solve for any potential you want.
Code located in the link below. Go to "Python Metaphysics tutorials" and then "Vid 3"
github.com/lukepolson/youtube...
Background Music: www.bensound.com

Опубликовано:

 

15 июл 2024

Поделиться:

Ссылка:

Скачать:

Готовим ссылку...

Добавить в:

Мой плейлист
Посмотреть позже
Комментарии : 64   
@eduardoh.m2072
@eduardoh.m2072 3 года назад
You sir, you are a freaking beast! Keep it up and I'll be here every single time!
@AJ-et3vf
@AJ-et3vf 2 года назад
So awesome! Watching your Python videos is just so immensely intellectually satisfying, inspiring, and galvanizing. Love your videos. The value of your videos cannot be overstated.
@ch0ng061
@ch0ng061 2 года назад
This is an absurdly cool video! I have two RU-vid accounts and went on both to like this video. Thanks for the great content. You're fantastic at explaining things.
@jessicalearnsthings9104
@jessicalearnsthings9104 2 года назад
this video was really clear and fun, and i love how excited you get. thanks for this!
@lookaway8496
@lookaway8496 2 года назад
This just became my favourite youtube channel!
@josephjones2211
@josephjones2211 3 года назад
When done in discrete form like this it reminds me of solving a system of lots of masses connected by springs. You get an almost identical tradiagonal matrix
@MrPSolver
@MrPSolver 3 года назад
Is the tridiagonal matrix for the system of masses symmetric?
@MrPSolver
@MrPSolver 3 года назад
The system of masses would be a great idea for a future video!
@josephjones2211
@josephjones2211 3 года назад
@@MrPSolver yes I did the problem just yesterday in my dynamical systems class
@MrPSolver
@MrPSolver 3 года назад
@@josephjones2211 interesting. I'm fairly certain that a symmetric tridiagonal matrix corresponds to a differential equation in the limit N->infinity. In the case of infinitely many infinitely small masses, this probably corresponds to a particular differential equation, like the wave equation or something. Worth looking into in a future video
@josephjones2211
@josephjones2211 3 года назад
@@MrPSolver It exactly does correspond to the wave equation as N->infinity. I have just coded it in to your notebook now. I am going to put it on my github jjones2000
@atuelvillegas2650
@atuelvillegas2650 2 года назад
This is pure gold!
@leonackermann3098
@leonackermann3098 Год назад
Road to Reality in the Background. You know the good stuff ! !
@Ash-oj5yv
@Ash-oj5yv 3 года назад
This is ridiculously high-quality content
@stauffap
@stauffap Год назад
Reminds me of the finite difference method that i used once to simulate a 2D Situation with the Navier Stokes Equations.
@gokuson5261
@gokuson5261 2 года назад
omg this is awesome thanks a lot!
@tayfun11100
@tayfun11100 2 года назад
Great tutorial, thanks...
@user-jv1ne9od4b
@user-jv1ne9od4b 2 года назад
素晴らしいビデオです!
@Piipolinoo
@Piipolinoo 3 года назад
Dude, I saw your video about the skydiving differential equation on reddit and was like "meeeh, not again, we don't need this another time". But this, this was brilliant! Great video, nice explanation of a scary problem. Keep it up. Could you also share the jupyter notebooks in the end? Also, I miss the link to your blog in the video description ;) Keep promoting yourself.
@MrPSolver
@MrPSolver 3 года назад
Thank you for your comment! :) I can certainly share the notebooks in the end: I'll make a git repository and upload the notebooks there. They should be easily accessible after that.
@Dark-tk9xu
@Dark-tk9xu Год назад
YOU are a great teacher (from a professor).
@bekarysbereketov4515
@bekarysbereketov4515 2 года назад
Wonderful content 🔥
@s.v.8662
@s.v.8662 Год назад
Really Great Video! Maybe next time you could compare your results with the analytical ones so in this case the qm harmonic oscillator model. I did it and the results are really close! So again thanks for the video
@soumadippakrashi5422
@soumadippakrashi5422 2 года назад
Life saviour ❤️
@ewinsart9390
@ewinsart9390 2 года назад
Hi Luke, Just a courtesy to let you know that I cited this video in a recent paper. Can't put any details because it gets deleted for containing a link but you can probably find them in your comments box. Great channel!
@narendrarathore1548
@narendrarathore1548 10 месяцев назад
You are awesome too good
@pro-bOi
@pro-bOi 2 года назад
I was trying to solve woods-saxon potential. I was messed up myself by numerov's algorithm. Then I found your video. Finally I solved the problem using this tridiagonal eigenvalues. But I want to make a request. Can you make a video on how to find all the eigenvalues in a certain range, for any 1-d potential using numerov's algorithm and shooting method.
@aliexpress.official
@aliexpress.official 3 года назад
that was so cool! saving this. can you show how to do that for the 3D equation? Also can you explain the way you wrote the d variable in your code?
@MrPSolver
@MrPSolver 3 года назад
The d variable is just the main diagonal of the main matrix, which has length N-1. The vector mL2V(y) has length N+1 (it's got 2001 points here), and we can ignore the first and last point of a numpy array by taking mL2V(y)[1:-1] (this is just an indexing trick). Then we add a scalar constant to out numpy array: namely 1/dy**2
@tototricepsclips
@tototricepsclips 3 года назад
Really enojoy this type of content, a good way in my optinion to learn Python if you know some physics
@bigh8438
@bigh8438 10 месяцев назад
damm looks similar to the euler lagrange equations we studied in aeroelasticity
@myapproach7785
@myapproach7785 2 года назад
Hi, amazing video. I wold like to ask I tried to model hydrogen atom with potential in a form of 1000/y, the problem of dividing by zero I solved by starting at y=0.0001 and since energy is not 0 at the right end- in infinity I extended the matrix by one element. The wavefunction looks quit similar to the hydrogen atom but I am not able to see the "infinity" in the graph even if I set N to 10000. Any suggestions will be highly appreciated, have a nice day.
@jceradnac1894
@jceradnac1894 3 года назад
how does this extend for the case of a 2 dimensional schrodinger equation?
@descriptionsuchandsuch4709
@descriptionsuchandsuch4709 3 года назад
very nice video! I have a lot of fun using python again! One little question, maybe it is a conventional thing. Then you divide by Δy^2, do you actually mean (Δy)^2 ? (since it is the 2nd (discrete) derivative?
@MrPSolver
@MrPSolver 3 года назад
Yes (∆y)^2. Thank you for mentioning this: I will be more careful with my notation in the future :)
@nobody2937
@nobody2937 2 года назад
Can you make your seriers of videos as in a particular playlist? Or anyway to group some of your videos, so that we can play all vidoes in a list ? ^_^
@giovannisilvadesouza6968
@giovannisilvadesouza6968 Год назад
MAKE DIRAC EQUATION
@youngidealist
@youngidealist 2 года назад
I have liked and subscribed. Please show me how to make a chemistry simulator. I want the orbitals to dance
@dwipesh
@dwipesh Год назад
very interesting. I have one simple doubt about the number of eigen states. The number of eigen states is equal to the dimension of the matrix, and the dimension of the matrix depends on the number of points we have chosen. So, the number of eigen states depends on the method of solutions!?
@rajkumarchakraborty3935
@rajkumarchakraborty3935 5 месяцев назад
However you can't use this method for scattering states, as you can't set the boundary conditions as psi(0)=psi(L)=0.
@shahidahmed2659
@shahidahmed2659 10 месяцев назад
What happened to reduced planck constant in the non dimensionalization step?
@ParasProtimBoruah
@ParasProtimBoruah 2 года назад
Can I get the link for the Wordpress paper ?
@maurocruz1824
@maurocruz1824 Год назад
This is similar to the problem 2.61 of Griffiths 3rd Edition
@mathoph26
@mathoph26 4 месяца назад
Does it work for a non-linear potential that depends on Psi (my example is: Laplacian of V = - Psi^2) ?
@moirangthemsanahal6695
@moirangthemsanahal6695 3 месяца назад
16:20 aren't energy levels in sho supposed to be equally spaced? In there m and L both are constants so your bar should be increasing linearly right? It seems a little off
@prasadpawar7027
@prasadpawar7027 3 года назад
How would one solve for energies if the wavefunction is leaking instead of being zero at 0 and L?
@MrPSolver
@MrPSolver 3 года назад
If the wave function was unbounded, you'd have to make the grid very dense, say 100000 points between 0 and 1 and make the potential very localized (maybe between like 0.5 and 0.51). This would be an approximation of a particle that was unbounded, since the infinite boundaries would be very far from the localized potential. (Sort of like how an air molecule is approximately "unbounded" inside a giant warehouse, even though the warehouse has walls)
@prasadpawar7027
@prasadpawar7027 3 года назад
@@MrPSolver Oh, makes sense. Thank you!
@shijianzhang9321
@shijianzhang9321 Год назад
Hi, it is interesting video. learn a lot from it. but there is a question about the demensionless step. in term of potenial, why directly become mL^2V(y)phi? why not mL^2V(Ly)phi?
@diegoquantico
@diegoquantico Год назад
I think that it's because V is linear. Then, we get one L from the Potential V and another L from the eigenfinction Psi
@diegoquantico
@diegoquantico Год назад
Hello. I've got one doubt. The energy output is quadratic. I think that it would be linear, since E = (n + 1/2)hw? Am I wrong?
@TueteKaffee
@TueteKaffee 7 месяцев назад
You are completely right with this. The difference betwen two eigenvalues of the harmonic oscillator has to be constant, so something is not right here. I'm not sure what it is though. Interestingly, the energy output isn't exactly quadratic, since np.sqrt(eigenvalues) also doesn't give a linear result.
@diegoquantico
@diegoquantico 7 месяцев назад
@@TueteKaffee Man. I've got the answer. In Numerical Analysis, methods based in finite differences to solve Sturm-Loiuville equations (The Schroedinger equations is a Sturm-Liouville like problem) approximate well the first eigenvalues. For practical purposes, this is what we want: the first energy levels. So, higher energies will not be well approximated.
@TueteKaffee
@TueteKaffee 7 месяцев назад
@@diegoquantico Thanks a lot :)
@gorlekarthikgowtham3903
@gorlekarthikgowtham3903 2 года назад
Sir,can you post the code here
@ismailnurmaghribi1702
@ismailnurmaghribi1702 Год назад
Why V(x) = L^2V(y)? What if the V(x) = x^4 or something else...?
@nna23101989
@nna23101989 Год назад
There is one thing that I do not understand. psi(x) is a complex number, but according to this video, it seems to be always real number?
@fluxionxyz
@fluxionxyz Год назад
I think psi(x) can be defined as real for the time-independent case.
@Karatemastert
@Karatemastert 2 года назад
137th like
@matthewkinahan3945
@matthewkinahan3945 3 года назад
Literally smarter than that guy on your wall
@bingodeagle
@bingodeagle 3 года назад
guy on the wall didn't need python
@ankidokolo
@ankidokolo 2 года назад
A very small tip. The topics you're representing are good and cool and all that. But I didn't learn them this way. Your videos aren't educational if that's what you want them to be. They are more like "look I can do cool stuff, and you know what, just take the code from the description, because I know that you can't do it yourself". Try to focus on the problem and discribe the code slowly. At the end we both KNOW, that these problems aren't a big deal. Peace.
@MrPSolver
@MrPSolver 2 года назад
Thanks for the very small tip!
@mariomuysensual
@mariomuysensual 2 года назад
Not true. This videos are top quality educational content, what are you talking about? You expect him to teach you QM from the ground or what?
Далее
Pendulum Motion in PYTHON
23:49
Просмотров 16 тыс.
A Quantum Problem Every Physics Student Should Attempt
32:46
[RU] Winline EPIC Standoff 2 Major | LAN | Final Day
9:48:47
Finite vs Infinite Well
6:10
Просмотров 11 тыс.
Grad School Be Like...
3:35
Просмотров 96 тыс.
Derivatives In PYTHON (Symbolic AND Numeric)
17:37
Просмотров 46 тыс.
ChatGPT Vs. Stack Overflow Be Like...
3:26
Просмотров 74 тыс.
How I’d learn ML in 2024 (if I could start over)
7:05
Monte Carlo Integration In Python For Noobs
15:32
Просмотров 149 тыс.
Be a Python Pro with Enumerate
0:55
Просмотров 1 млн