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?
"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!
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
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!
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
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!
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
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?
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
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
![[Full] 4 ต่อ 4 Celebrity EP.922 | 10 พ.ย. 67 | one31](http://i.ytimg.com/vi/zmV5HRIxIBI/mqdefault.jpg)
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?
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!
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.
Stat mech is the most beautiful thing in Physics.
It really is!
Really cool vids man! Stumbled across one of your reddit posts and thought I reconised your face haha, subbed!
Hahaha cheers dude 😂
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....
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!
Why doesn't this have more views and likes? That was simply amazing!
Great video! Numba seems super powerful, my code in undergrad for this similarly took forever to complete :).
Clear and straightforward, excellent resource for the topic
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!
dude I saw this on reddit and you're so underrated, i'll ask my friends to subscribe to you
Finally I understand the meaning of the Ising model. Thanks for this great 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
There’s that njit :-) Amazing video!
waiting for your sympy video...
Whao,
You just saved a life! Thanks man!
This is so on point👌
You are an absolute genius!
Thank you so very much. It saved a lot of time. I do appreciate your knowledge.
These are some great videos! I would love to see your process of putting the notebooks together.
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!
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
Great content. You are too good, Luke.
Thank you Shadab!
Thank you so much for the content!
Amazing work. Bravo !
This was really great. Thanks!
You're amazing 👏 a lot of thanks for this video
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!
amazing work mate! are you a physicist?
Nice work…❤
Does the energy calculation allow for periodic boundary conditions on our lattice?
that cuttings :D
How do I view the final spin array bitmap after the Metropolis algorithm has worked upon it?
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
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?
How would you introduce the presence of an external magnetic field to this code?
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! Is there any way to display the new ising model with the updated spins?
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..
Can u do the same for Classical XY model
Your mic is clipping. I recommend lowering the gain during recording because it sounds bad
Ab
So
Lu
Te
Ly
Brilliant
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
Is the code complete?
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