The heterogeneous oxidation of organic aerosol by hydroxyl radicals (OH) can proceed through two general pathways: functionalization, in which oxygen functional groups are added to the carbon skeleton, and fragmentation, in which carbon–carbon bonds are broken, producing higher volatility, lower molecular weight products. An ongoing challenge is to develop a quantitative molecular description of these pathways that connects the oxidative evolution of the average aerosol properties (e.g. size and hygroscopicity) to the transformation of free radical intermediates. In order to investigate the underlying molecular mechanism of aerosol oxidation, a relatively compact kinetics model is developed for the heterogeneous oxidation of squalane particles by OH using free radical intermediates that convert reactive hydrogen sites into oxygen functional groups. Stochastic simulation techniques are used to compare calculated system properties over ten oxidation lifetimes with the same properties measured in experiment. The time-dependent average squalane aerosol mass, volume, density, carbon number distribution of scission products, and the average elemental composition are predicted using known rate coefficients. For functionalization, the calculations reveal that the distribution of alcohol and carbonyl groups is controlled primarily by the initial OH abstraction rate and to lesser extent by the branching ratio between secondary peroxy radical product channels. For fragmentation, the calculations reveal that the formation of activated alkoxy radicals with neighboring functional groups controls the molecular decomposition, particularly at high O/C ratios. This kinetic scheme provides a framework for understanding the oxidation chemistry of a model organic aerosol and informs parameterizations of more complex systems.