Sedimentation of active Brownian particles

Codes used in the scientific publication:
M. Mangeat, S. Chakraborty, A. Wysocki, and H. Rieger, Stationary particle currents in sedimenting active matter wetting a wall, Phys. Rev. E 109, 014616 (2024). Preprint available on arXiv.

Interacting ABPs

C++ code on the sedimention of interacting active Brownian particles in a box.
Exportations: dynamics with Delta t = 1 between two datafiles, steady state density, polarization and current.
Compile: g++ interactingABPs.cpp -lgsl -lgslcblas -lm -O3 -s -o interactingABPs.out.
Run: ./interactingABPs.out -parameter=value.
List of parameters: Pe, alpha, F0, LX, LY, init, RAN, dt, teq, tmax, Npart (details as comments in the code).
Generate the movie: python figure_interactingABPs_dynamics.py -parameter=value (parameters: Pe, alpha, F0, LX, LY, init, ran, NCPU, movie).
Generate the steady-state figure: python figure_interactingABPs_steadystate.py -parameter=value (parameters: Pe, alpha, F0, LX, LY, init, Nran).

Ideal ABPs

FreeFem++ code on the sedimention of non-interacting active Brownian particles in a box.
Exportations: steady state density and polarization.
Run: FreeFem++ -v 0 -nw idealABPs.edp -parameter value.
List of parameters: vp, vg, LX, LY, Nx, Ny, Nt (details as comments in the code).
Generate the steady-state figure: python figure_idealABPs_steadystate.py -parameter=value (parameters: vp, vg, LX, LY).