The finite elemental model for liposomal release was based on the model of a disk electrode detecting exocytosis in our previous work 1 , with a new module for calculating the adsorption-desorption dynamics was added (scheme in Figure ). Schematically, a circle (the radius, Rlip=188 nm) and a pore on the membrane with variable size (Rp) was used to simulate the electroporated liposome membrane where a total quantity of cargo (Q=48.4 fC, the average experimental charge for spikes from DA/5-HT groups) was initially set inside the vesicles, but separated into adsorbed spcies (Qs, s indicates the membrane inner surface) and free species (Qf, f indicates freely moving species). The values for Qs and Qf were calculated a priori based on the Langmuir adsorption model. At time=0 s, the initial surface concentration (Cs,0) and initial concentration of intraliposomal freely moving molecules (Cf,0) was assumed to be at equilibrium, so the initial adsorption rate (rads,0) was equal to the desorption rate (rdes,0).