Eigenstates of ANY 1D Potential in PYTHON

[@eduardoh.m2072]

You sir, you are a freaking beast! Keep it up and I'll be here every single time!

[@josephjones2211]

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]

Is the tridiagonal matrix for the system of masses symmetric?

[@MrPSolver]

The system of masses would be a great idea for a future video!

[@josephjones2211]

@@MrPSolver yes I did the problem just yesterday in my dynamical systems class

[@MrPSolver]

@@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]

@@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

[@AJ-et3vf]

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]

This is an absurdly cool video! I have two KZbin accounts and went on both to like this video. Thanks for the great content. You're fantastic at explaining things.

[@jessicalearnsthings9104]

this video was really clear and fun, and i love how excited you get. thanks for this!

[@lookaway8496]

This just became my favourite youtube channel!

[@leonackermann3098]

Road to Reality in the Background. You know the good stuff ! !

[@Ash-oj5yv]

This is ridiculously high-quality content

[@Piipolinoo]

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]

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.

[@ewinsart9390]

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!

[@stauffap]

Reminds me of the finite difference method that i used once to simulate a 2D Situation with the Navier Stokes Equations.

[@Dark-tk9xu]

YOU are a great teacher (from a professor).

[@aliexpress.official]

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]

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

[@moirangthemsanahal6695]

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

[@ばばおおさか]

素晴らしいビデオです！

[@atuelvillegas2650]

This is pure gold!

[@ParasProtimBoruah]

Can I get the link for the Wordpress paper ?

[@gokuson5261]

omg this is awesome thanks a lot!

[@pro-bOi]

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.

[@tayfun11100]

Great tutorial, thanks...

[@jceradnac1894]

how does this extend for the case of a 2 dimensional schrodinger equation?

[@soumadippakrashi5422]

Life saviour ❤️

[@prasadpawar7027]

How would one solve for energies if the wavefunction is leaking instead of being zero at 0 and L?

[@MrPSolver]

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]

@@MrPSolver Oh, makes sense. Thank you!

[@bigh8438]

damm looks similar to the euler lagrange equations we studied in aeroelasticity

[@mathoph26]

Does it work for a non-linear potential that depends on Psi (my example is: Laplacian of V = - Psi^2) ?

[@nobody2937]

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 ? ^_^

[@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

[@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]

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]

What happened to reduced planck constant in the non dimensionalization step?

[@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]

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]

@@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]

@@diegoquantico Thanks a lot :)

[@giovannisilvadesouza6968]

MAKE DIRAC EQUATION

[@bekarysbereketov4515]

Wonderful content 🔥

[@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]

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

[@tototricepsclips]

Really enojoy this type of content, a good way in my optinion to learn Python if you know some physics

[@myapproach7785]

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.

[@descriptionsuchandsuch4709]

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]

Yes (∆y)^2. Thank you for mentioning this: I will be more careful with my notation in the future :)

[@youngidealist]

I have liked and subscribed. Please show me how to make a chemistry simulator. I want the orbitals to dance

[@maurocruz1824]

This is similar to the problem 2.61 of Griffiths 3rd Edition

[@narendrarathore1548]

You are awesome too good

[@ismailnurmaghribi1702]

Why V(x) = L^2V(y)? What if the V(x) = x^4 or something else...?

[@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]

I think psi(x) can be defined as real for the time-independent case.

[@gorlekarthikgowtham3903]

Sir,can you post the code here

[@matthewkinahan3945]

Literally smarter than that guy on your wall

[@bingodeagle]

guy on the wall didn't need python

[@Karatemastert]

137th like

[@ankidokolo]

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]

Thanks for the very small tip!

[@mariomuysensual]

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?