PhysSandbox
Classical MechanicsWaves & SoundElectricity & MagnetismOptics & LightGravity & OrbitsLabs
🌙Astronomy & The Sky🌡️Thermodynamics🌍Biophysics, Fluids & Geoscience📐Math Visualization🔧Engineering🧪Chemistry

Related simulators

Continue with similar topics in this category — or all 46 in Thermodynamics.

View category →
FeaturedSchool

Ideal Gas Simulator

Launch Simulator

Bouncing particles in a box. See PV=nRT in action.

NewSchool

Forest Fire (Cellular Automaton)

Launch Simulator

Toroidal lattice: trees regrow with probability p, lightning ignites trees with probability f, fire spreads to four neighbors; explore intermittent bursts and SOC-like phenomenology.

NewSchool

Wolfram Elementary Cellular Automata

Launch Simulator

1D binary CA with Wolfram code R∈[0,255]: space–time raster, periodic or null boundaries, single-cell or random IC, rule-bit strip, and informal class hints for textbook rules (e.g. 30, 90, 110, 184).

NewSchool

Kuramoto Oscillators

Launch Simulator

All-to-all coupled phases θ_i with intrinsic frequencies ω_i: order parameter r = |N⁻¹ Σ e^{iθ_i}| grows as coupling K crosses the synchronization window—classic mean-field rhythm transition.

NewSchool

Schelling Segregation Model

Launch Simulator

Two agent types on a toroidal lattice with vacancies: unhappy agents swap into random empty cells when fewer than a fraction τ of occupied Moore neighbors share their type—mild local rules, global clustering.

NewSchool

Eden Growth (Lattice)

Launch Simulator

Square-lattice Eden model: each step occupies a uniformly random empty von-Neumann neighbor of the cluster; compact growth with a rough interface—track perimeter P vs 2√(πN) as a roughness proxy.

PhysSandbox

Interactive physics, chemistry, and engineering simulators for students, teachers, and curious minds.

Physics

  • Classical Mechanics
  • Waves & Sound
  • Electricity & Magnetism

Science

  • Optics & Light
  • Gravity & Orbits
  • Astronomy & The Sky

More

  • Thermodynamics
  • Biophysics, Fluids & Geoscience
  • Math Visualization
  • Engineering
  • Chemistry

© 2026 PhysSandbox. Free interactive science simulators.

PrivacyTermsContact
Home/Thermodynamics/Vicsek Model (Active Matter)

Vicsek Model (Active Matter)

Self-propelled particles on a torus align headings with neighbors within radius R, add noise η, then drift at v₀; polarization V_a = |⟨e^{iθ}⟩| captures the flocking crossover with density ρ = N/L².

Simulation

12

Parameters

0.35
0.07
0.045
180

Classic Vicsek flocking: at fixed box, raising density (N) and lowering η or enlarging R promotes orientational order—V_a jumps across a crossover reminiscent of a flocking transition (finite-N blur).

Shortcuts

  • •Space / Enter — play / pause
  • •P — pause / resume
  • •R — reset positions & headings

Measured values

Polarization V_a = |⟨e^{iθ}⟩|0.0000
ρ = N/L²180.0
Time0.00
N180

About this model

The Vicsek model is the canonical active-matter toy for flocking: self-propelled particles on a torus align their heading with the average direction of neighbors within interaction radius R, add uniform noise ξ ∈ [−η, η], and advect at speed v₀. The polarization V_a = |N⁻¹ Σ_i e^{iθ_i}| measures global orientational order; at fixed box area, ρ ∝ N. This page uses the standard discrete-time update, draws the square torus with velocity arrows, shows V_a as a gold pointer, and plots V_a(t); sweeping η, R, v₀, and N reproduces the familiar noise–density crossover (rounded on finite N).

Who it's for: Students in statistical physics of complex systems, soft matter, and collective behavior.

Key terms

  • Vicsek model
  • Active matter
  • Flocking
  • Order parameter
  • Alignment noise
  • Periodic boundaries

How it works

Self-propelled particles on a torus align their headings with neighbors within distance R, then move at speed v₀. Uniform noise η in the update competes with local coupling to produce a polarization order parameter V_a — the textbook Vicsek active-matter toy for flocking.

Key equations

θ_i(t+Δ) = Arg(Σ_{|r_ij|<R} e^{iθ_j}) + ξ_i with ξ_i ~ Uniform(−η, η), then r_i ← r_i + v₀ e^{iθ_i} Δ with periodic boundaries. Polarization V_a = |N⁻¹ Σ_i e^{iθ_i}|.

Frequently asked questions

Why align before moving?
This is the usual discrete-time Vicsek update: neighbors are defined from positions at the current timestep, the new heading follows the local alignment rule, then each particle moves one step along that heading.
Is this a sharp phase transition?
In large systems, increasing density or lowering noise leads to a sharp rise in polarization; on a finite box the transition is blurred, but the qualitative crossover remains: more neighbors inside R and smaller η tend to raise V_a.
Why periodic boundaries?
A torus removes walls and keeps the density statistically uniform, matching the standard numerical setup used in Vicsek-model discussions.