The Ising Model in Python: Statistical Mechanics and Permanent Magnets

[@adwaitnaravane5285]

Stat mech is the most beautiful thing in Physics.

[@ricomajestic]

It really is!

[@muhammadhidayat1337]

I'm currently in undergraduate trying to understand what ising model is, and suddenly found this video. Thank you so much!

[@generaliroh842]

Sama ngab

[@charltonwake6813]

"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!

[@ganstabreakincity]

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.

[@peabrane8067]

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

[@estebanrodriguez3922]

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!

[@AlessandroAmbler]

Really cool vids man! Stumbled across one of your reddit posts and thought I reconised your face haha, subbed!

[@MrPSolver]

Hahaha cheers dude 😂

[@JonathonRiddell]

Great video! Numba seems super powerful, my code in undergrad for this similarly took forever to complete :).

[@corneldevilliers461]

Clear and straightforward, excellent resource for the topic

[@adityadhardwivedi634]

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!

[@deltascience7252]

dude I saw this on reddit and you're so underrated, i'll ask my friends to subscribe to you

[@CarlosRoberto-vb9gi]

Finally I understand the meaning of the Ising model. Thanks for this great video.

[@randomstuff9960]

Why doesn't this have more views and likes? That was simply amazing!

[@markushaysnielsen4885]

There’s that njit :-) Amazing video!

[@ayabongandlungwana4407]

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

[@AJ-et3vf]

Awesome! Love your videos! This channel is seriously underrated and needs more views :) very quality and enjoyable videos.

[@westmoreland5889]

waiting for your sympy video...

[@lasinar9090]

Nice work…❤

[@benburdick9834]

These are some great videos! I would love to see your process of putting the notebooks together.

[@MrPSolver]

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.

[@laika4048]

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?

[@godwinosabutey2876]

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?

[@GabrielleduVent]

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!

[@NGBigfield]

You are an absolute genius!

[@godwinosabutey2876]

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

[@quornnugget7799]

Does the energy calculation allow for periodic boundary conditions on our lattice?

[@henryouang]

This was really great. Thanks!

[@joyzhong7279]

Thank you so much for the content!

[@chisomdim2031]

Whao, You just saved a life! Thanks man! This is so on point👌

[@TheAbdel]

Amazing work. Bravo !

[@aaronthomas6434]

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!

[@ahxmeds]

Great content. You are too good, Luke.

[@MrPSolver]

Thank you Shadab!

[@physicsguy9774]

amazing work mate! are you a physicist?

[@md.tohidulislam7041]

Thank you so very much. It saved a lot of time. I do appreciate your knowledge.

[@drissboutat6430]

You're amazing 👏 a lot of thanks for this video

[@cscsquad1902]

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?

[@anuranbhunia]

How do I view the final spin array bitmap after the Metropolis algorithm has worked upon it?

[@thomsday]

that cuttings :D

[@homa_ghaiedi2626]

How would you introduce the presence of an external magnetic field to this code?

[@حيدرالعبيدي-ج4د]

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

[@sadaf6295]

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

[@sagarpatel3043]

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.

[@J-rb9dx]

Hello! Is there any way to display the new ising model with the updated spins?

[@muditgaur9736]

Can u do the same for Classical XY model

[@sadaf6295]

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

[@clf400]

Your mic is clipping. I recommend lowering the gain during recording because it sounds bad

[@raa9558]

Ab So Lu Te Ly Brilliant

[@mB-mp6io]

where the code disappear

[@kasraburiaey7862]

Is the code complete?

[@prikarsartam]

I get trouble with the "plt.style.use(['science','notebook', 'grid'])" part, would you kindly provide solution?

[@Eta_Carinae__]

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!

[@banditadas3420]

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