APS Excellence in Physics Education Award
November 2019

Education Prize Logo
Science SPORE Prize
November 2011

NSF Logo
The Open Source Physics Project is supported by NSF DUE-0442581.

Diffusion of particles along an 1-d lattice Documents

This material has 3 associated documents. Select a document title to view a document's information.

Main Document

Diffusion of particles along an 1-d lattice 

written by Kostas Papamichalis

This simulation uses a Markov chain to simulate the diffusion of N particles along an one dimension, finite lattice, toward the state of equilibrium. The particles are distributed along the 1D lattice cells. In a time interval of length Dt, each particle can perform just one jump between neighboring cells with a certain transition probability. The evolution of the system from its initial state to the final state of equilibrium, is described by a master equation. The initial state of the system and the transition probabilities are selected by the user. The main objective of the simulation is the confirmation of the theoretical proposition that "irrespectively of the form of the initial distribution, the system converges to a certain equilibrium state which is determined by the transition probabilities".  

In a sequence of time-moments the program of the simulation counts the real number of particles in every cell. The intermediate states of the system between the initial state and the final state of equilibrium are depicted by a varying histogram and a sequence of changing cell-colors. On the other hand, the distribution of the particles at the equilibrium state has been determined according to the theoretical model and it is depicted in the same graphs. The user compares the real-time data with the theoretical predictions.

Finally, a Lyapunov functional H is determined for the system. Each time-moment, the value of H is uniquely determined by the corresponding distribution of the particles in the cells of the lattice. The variation of H with time is depicted in real time. By using this graph, the user can estimate the relaxation time of the process toward the equilibrium-state.

Last Modified November 9, 2022

Supplemental Documents

Diffusion of Partilces Along a 1D Lattice Theory 

A description of the Diffusion of Partilces Along a 1D Lattice and the theory of the Markov Process used in this simulation.

Last Modified November 9, 2022

This file is included in the full-text index.

Source Code Documents

Diffusion of Particles on 1D Lattice Source Code 

EJS source code for Diffusion of Particles on 1D Lattice Model.

Last Modified November 9, 2022

OSP Projects:
Open Source Physics - EJS Modeling
Physlet Physics
Physlet Quantum Physics
STP Book