Characterization of spatiotemporal electroactive anodic biofilm activity distribution using 1D simulations

Activity distribution limitation in electroactive biofilm remains an unclear phenomenon. Some observations using confocal microscopy have shown notable difference between activity close to the anode and activity at the liquid interface. A numerical model is developed in this work to describe biofilm growth and local biomass segregation in electroactive biofilm. Under our model hypothesis, metabolic activity distribution in the biofilm results from the competition between two limiting factors: acetate diffusion and electronic conduction in the biofilm. Influence of inactive biomass fraction (i.e. non-growing biomass fraction) properties (such as conductivity and density) is simulated to show variation in local biomass distribution. Introducing a dependence of effective diffusion to local density leads to a drastic biomass fraction segregation. Increasing density of inactive fraction reduces significantly acetate diffusion in biofilm, enhances biomass activity on the outer layer (liquid/biofilm interface) and maintains inner core largely inactive. High inactive fraction conductivity enhances biomass activity in the outer layer and enhances current production. Hence, investment in extracellular polymer substance (EPS), anchoring redox components, is benefit for biofilm electroactivity. However, under our model hypothesis it means that conductivity should be two order lower than biofilm conductivity reported in order to observe inner core active biomass segregation.