Villar-Sepulveda, EdgardoChampneys, Alan R.Cusseddu, DavideMadzvamuse, Anotida2026-03-232026-03-232026Villar-Sepúlveda, E., Champneys, A.R., Cusseddu, D. & Madzvamuse, A. 2026, 'Pattern formation of bulk-surface reaction-diffusion systems in a ball', SIAM Journal on Applied Mathematics, vol. 86, no. 1, pp. 21-51, doi : 10.1137/24M1671037.0036-1399 (print)1095-712X (online)10.1137/24M1671037http://hdl.handle.net/2263/109252Weakly nonlinear amplitude equations are derived for the onset of spatially extended patterns on a general class of 𝑛-component bulk-surface reaction-diffusion systems in a ball, under the assumption of linear kinetics in the bulk and coupling Robin-type boundary conditions. Linear analysis shows conditions under which various pattern modes can become unstable to either generalized pitchfork or transcritical bifurcations, depending on the parity of the spatial wavenumber. Weakly nonlinear analysis is used to derive general expressions for the multicomponent amplitude equations of different patterned states. These reduced-order systems are found to agree with prior normal forms for pattern formation bifurcations with 𝑂⁡(3) symmetry, and they provide information on the stability of bifurcating patterns of different symmetry types. The analysis is complemented with numerical results using a dedicated bulk-surface finite element method. The theory is illustrated in two examples: a bulk-surface version of the Brusselator and a four-component bulk-surface cell-polarity model.en© 2026 Society for Industrial and Applied Mathematics.Bulk-surface PDEsBulk-surface finite element methodBifurcationSymmetry breakingWeakly nonlinear analysisPattern formationBulk-surface reaction-diffusion systemsPostprint Article