The composition vector of an assembly is written as: v=n_1\cdots n_{N_G}. Where n_1\cdots n_{N_G} are the molecular counts of lipid type
i within the assembly, and NG is how many different lipid types exist (
repertoire size). The change in the count of molecule type
i is described by: : \frac{dn_i}{dt} = (k_f \rho_i N-k_b n_i) \left(1+\sum_{j=1}^{N_G}\beta_{ij} \frac{n_j}{N}\right) k_f and k_b are the basal forward (joining) and backward (leaving) rate constants,
βij is a non-negative rate enhancement exerted by molecule type
j within the assembly on type
i from the environment, and ρ is the environmental concentration of each molecule type.
β is viewed as a directed, weighted,
complex network. The assembly current size is N=\sum_{i=1}^{N_G}n_i. The system is kept away from equilibrium by imposing a fission action once the assembly reaches a maximal size, Nmax, usually in the order of NG. This splitting action produces two progeny of same size, and one of which is grown again. The model is subjected to a
Monte Carlo algorithm based simulations, using
Gillespie algorithm. ==Selection==