Mean-Field Theory of Chiral Active Model B: Arrested Coarsening and Chiral Fingering Instabilities

This Zenodo repository contains the code, data, and animations accompanying the manuscript Mean-Field Theory of Chiral Active Model B: Arrested Coarsening and Chiral Fingering Instabilities. The repository contains the following files: PDE_linear_stability.ipynbJupyter notebook used to perform the angular linear stability analysis of the continuum PDE. The PDE is discretized in the radial direction, resulting in an eigenvalue problem for angular perturbation modes. This notebook was used to generate the dispersion relations and the stability diagram shown in Fig. 4a and Fig. 4b. PDE_solver_notebook.ipynbJupyter notebook used to solve the continuum PDE numerically. The simulations are performed with a pseudo-spectral method in space and a semi-implicit Euler scheme in time. This notebook was used to generate the time-simulation snapshots shown in Fig. 4c and the spacetime plot shown in Fig. 5. ODE_solver_notebook.ipynbJupyter notebook used to solve the coupled system of mean-field lattice ODEs by numerical time integration. This notebook was used to generate Fig. 3. Animator.pyPython script used to create animations from the PDE time-simulation data. Chiral_Ising_Model.nbMathematica notebook containing a numerical solver for the coupled set of mean-field lattice ODEs. The notebook also shows how the continuum PDE is obtained from the discrete mean-field ODEs by taking the continuum limit. For reproducing the ODE simulations, we recommend using ODE_solver_notebook.ipynb. Animations_J0.7_tau1_L50_N200_dt0.005_nsteps3000000_savesteps200_deltas(6.6,7.6)_circular_eta_0.0000001.mp4Animation of the PDE time simulations for the chiral driving strengths \delta=6.6 and \delta=7.6. The simulations use \tau=1, L=50, N=200, \Delta t=0.005, and an initial angular perturbation amplitude \eta=10^{-7}. phase_data_Rbox25_R10_w1_J0.350_to_0.813_delta0.00_to_10.00_nr205_nmin8_nmax30.pklDataset produced by the angular linear stability analysis. It contains the stability data used to construct the stability diagram shown in Fig. 4b. chiral_results_for_figure_J0.4_tau1_L50_N200_dt0.005_nsteps2000000_savesteps200_deltas(5.6,6.6,7.6)_circular_correct_bulk_values_0.00001.pklDataset containing PDE time-simulation results for J=0.4, \tau=1, L=50, N=200, \Delta t=0.005, and chiral driving strengths \delta=5.6, \delta=6.6, and \delta=7.6. The initial angular perturbation amplitude is \eta=10^{-5}. chiral_results_for_figure_J0.4_tau1_L50_N200_dt0.005_nsteps3000000_savesteps200_deltas(6.6,7.6)_circular_correct_bulk_values_0.0000001.pklDataset containing PDE time-simulation results for J=0.4, \tau=1, L=50, N=200, \Delta t=0.005, and chiral driving strengths \delta=6.6 and \delta=7.6. The initial angular perturbation amplitude is \eta=10^{-7}. This dataset was used to generate the time-simulation snapshots in Fig. 4c for \delta=6.6 and \delta=7.6.