Abstract. A new method is introduced for estimating single-trial magnetoor electro-encephalography (M/EEG), based on a non-linear fit of timefrequency atoms. The method can be applied for transient activity (e.g. event-related potentials) as well as for oscillatory activity (e.g. gamma bursts), and for both evoked or induced activity. In order to benefit from all the structure present in the data, the method accounts for (i) spatial structure of the data via multivariate decomposition, (ii) time-frequency structure via atomic decomposition and (iii) reproducibility across trials via a constraint on parameter dispersion. Moreover, a novel iterative method is introduced for estimating the initial time-frequency atoms used in the non-linear fit. Numerical experiments show that the method