"Did this in undergrad without numba"'; gave me flashbacks when I had to do Ising Models in undergrad, taking my side panels off and setting up 2 desk fans to blow into my pc in fear of it over heating. God bless numba!
Correction: In the function I used to originally get the energy of the lattice, I should have divided the return value by 2 to avoid double counting of pairs. This doesn't effect the dynamics of the system (since this only depends on change in energy), but it slightly changes the energy plot.
why did you calculate the total energy of all the lattice? I think you should only calculate energy for one site i,j and then because ergodic theorem it would be equal to the energy of all the lattice. Or i am missing something?
Thank you! Very nice work and good explanation there. The Numba and the codes work incredibly well and fast. I was wondering how higher order interactions, such as those involving three bodies as opposed to two, might be implemented in a similar manner. Any advice you can give me?
Here before this channel has million subscribers, give it some time. Pure knowledge, I wish I can learn and collaborate with you given the time. Appreciate your work.
Glad I got internship in the very first summer at my Institute under condensed matter professor just because of curiosity of phase transition that occurs in matter. I'll get to learn many things, thanks for this video, man!
Wow this video has helped me to complete my assignment, although I'm using C to simulate the results. It gave me an idea about the how results should look like and better understand for the whole Ising model.. Thank you bro
Thank you! Very nice work and good explanation there. The Numba and the codes work incredibly well and fast. I was wondering how higher order interactions, such as those involving three bodies as opposed to two, might be implemented in a similar manner. Any advice you can give me? @Mr. P Solver
okay how do you apply it to a real system? for example, we want to study magnetic behavior of some ferromagnetic metal (ie cobalt) which has a hcp structure. In that case, you would need the exact positions of the spins. Do you have any idea?
Thanks for uploading this video! I'm a python newbie who needs to run Ising models and this clarified quite a few things. Few questions: 1. For the total initial energy, if I am using a very specific energy equation, do I need to change it from convolution to something else? My set up is that it starts with all down or up spins. 2. Since my setup is starting with 100% down or upspin, how would that affect the probability equation (2nd line after #3/4. change state...)? 3. How would I change the probability function to incorporate my energy equation? Thanks in advance!
Great Video! Helping me so much for my assignment! One question, for the full Hamiltonian where you would include the external energy contribution i.e the external field multiplied by the sum over all spins, would you simply add this into the get energy function or amend it somehow within the metropolis algorithm? Kepp up the brilliant work!
Hey .. How to make code in python for.. Spin correlation function for Heisenberg spin Hamiltonian.. ?? I am new to python.. Please give me some your suggestion
Hello..I am a PhD student in physics from Iraq..I hope you can help me find codes in the Python program to study the Fe(phen) molecule to investigation of the spin crossover phenomenon by monte carlo (ising model 2D)..with many thanks to you
Hey.. Can you make video on Heisenberg spin hamiltonian in python like yo find eigen value and eigen vectors and what happens when N change.. M a total Biggener in python.. But i have write a program.. Please help me..
I do not understand. Why positive spin sample and negative spin sample are separated in the starting? If +1 means spin-up and -1 means spin-down one sample is enough. Forgive my lack of intelligence.
Hi! I tried more or less the same thing, except the way I noised the lattice is by randomly flipping every cell per timestep according to a bernoulli trial with the chance of success standing in for the chance to flip (I use numpy.random.binomial, and just set the no. trials to 1). Is there a reason why you flip one cell per timestep instead of a random number? Cheers!
plt.style.use(['science', 'notebook','grid']) isnt working for me , its showing lot of error with OSError: 'science' not found in the style library and input is not a valid URL or path; see `style.available` for list of available styles
I'm currently in undergraduate trying to understand what ising model is, and suddenly found this video. Thank you so much!
Sama ngab
"Did this in undergrad without numba"'; gave me flashbacks when I had to do Ising Models in undergrad, taking my side panels off and setting up 2 desk fans to blow into my pc in fear of it over heating. God bless numba!
Correction: In the function I used to originally get the energy of the lattice, I should have divided the return value by 2 to avoid double counting of pairs. This doesn't effect the dynamics of the system (since this only depends on change in energy), but it slightly changes the energy plot.
why did you calculate the total energy of all the lattice? I think you should only calculate energy for one site i,j and then because ergodic theorem it would be equal to the energy of all the lattice. Or i am missing something?
Thank you! Very nice work and good explanation there. The Numba and the codes work incredibly well and fast. I was wondering how higher order interactions, such as those involving three bodies as opposed to two, might be implemented in a similar manner. Any advice you can give me?
Stat mech is the most beautiful thing in Physics.
It really is!
Here before this channel has million subscribers, give it some time. Pure knowledge, I wish I can learn and collaborate with you given the time. Appreciate your work.
Really cool vids man! Stumbled across one of your reddit posts and thought I reconised your face haha, subbed!
Hahaha cheers dude 😂
Why doesn't this have more views and likes? That was simply amazing!
Clear and straightforward, excellent resource for the topic
This is outstanding great quality content. I've always wanted to make similar videos on youtube but PhD and academia are just too frantic! Keep it up!
Wow that convolve operation to get the energy is so smart. I never thought of that
Would probably gave me a 10x speedup for my PhD....
Finally I understand the meaning of the Ising model. Thanks for this great video.
dude I saw this on reddit and you're so underrated, i'll ask my friends to subscribe to you
Great video! Numba seems super powerful, my code in undergrad for this similarly took forever to complete :).
Glad I got internship in the very first summer at my Institute under condensed matter professor just because of curiosity of phase transition that occurs in matter. I'll get to learn many things, thanks for this video, man!
There’s that njit :-) Amazing video!
Awesome! Love your videos! This channel is seriously underrated and needs more views :) very quality and enjoyable videos.
Wow this video has helped me to complete my assignment, although I'm using C to simulate the results. It gave me an idea about the how results should look like and better understand for the whole Ising model.. Thank you bro
waiting for your sympy video...
You are an absolute genius!
Great content. You are too good, Luke.
Thank you Shadab!
Whao,
You just saved a life! Thanks man!
This is so on point👌
Amazing work. Bravo !
Does the energy calculation allow for periodic boundary conditions on our lattice?
How do I view the final spin array bitmap after the Metropolis algorithm has worked upon it?
This was really great. Thanks!
Thank you! Very nice work and good explanation there. The Numba and the codes work incredibly well and fast. I was wondering how higher order interactions, such as those involving three bodies as opposed to two, might be implemented in a similar manner. Any advice you can give me? @Mr. P Solver
Thank you so very much. It saved a lot of time. I do appreciate your knowledge.
Nice work…❤
Hello! Is there any way to display the new ising model with the updated spins?
Thank you so much for the content!
These are some great videos! I would love to see your process of putting the notebooks together.
How would you introduce the presence of an external magnetic field to this code?
okay how do you apply it to a real system? for example, we want to study magnetic behavior of some ferromagnetic metal (ie cobalt) which has a hcp structure. In that case, you would need the exact positions of the spins. Do you have any idea?
Thanks for uploading this video! I'm a python newbie who needs to run Ising models and this clarified quite a few things. Few questions:
1. For the total initial energy, if I am using a very specific energy equation, do I need to change it from convolution to something else? My set up is that it starts with all down or up spins.
2. Since my setup is starting with 100% down or upspin, how would that affect the probability equation (2nd line after #3/4. change state...)?
3. How would I change the probability function to incorporate my energy equation?
Thanks in advance!
amazing work mate! are you a physicist?
Great Video! Helping me so much for my assignment! One question, for the full Hamiltonian where you would include the external energy contribution i.e the external field multiplied by the sum over all spins, would you simply add this into the get energy function or amend it somehow within the metropolis algorithm? Kepp up the brilliant work!
You're amazing 👏 a lot of thanks for this video
Hey .. How to make code in python for.. Spin correlation function for Heisenberg spin Hamiltonian.. ??
I am new to python.. Please give me some your suggestion
Hello..I am a PhD student in physics from Iraq..I hope you can help me find codes in the Python program to study the Fe(phen) molecule to investigation of the spin crossover phenomenon by monte carlo (ising model 2D)..with many thanks to you
Can u do the same for Classical XY model
that cuttings :D
Hey.. Can you make video on Heisenberg spin hamiltonian in python like yo find eigen value and eigen vectors and what happens when N change..
M a total Biggener in python.. But i have write a program..
Please help me..
Is the code complete?
I do not understand. Why positive spin sample and negative spin sample are separated in the starting? If +1 means spin-up and -1 means spin-down one sample is enough. Forgive my lack of intelligence.
where the code disappear
Your mic is clipping. I recommend lowering the gain during recording because it sounds bad
Ab
So
Lu
Te
Ly
Brilliant
I get trouble with the "plt.style.use(['science','notebook', 'grid'])" part, would you kindly provide solution?
Hi! I tried more or less the same thing, except the way I noised the lattice is by randomly flipping every cell per timestep according to a bernoulli trial with the chance of success standing in for the chance to flip (I use numpy.random.binomial, and just set the no. trials to 1). Is there a reason why you flip one cell per timestep instead of a random number? Cheers!
plt.style.use(['science', 'notebook','grid']) isnt working for me , its showing lot of error with OSError: 'science' not found in the style library and input is not a valid URL or path; see `style.available` for list of available styles