Many methods have been described for automated clustering analysis of complex flow cytometry data, but so far the goal to efficiently estimate multivariate densities and their modes for a moderate number of dimensions and potentially millions of data points has not been attained. We have devised a novel approach to describing modes using second order polynomial histogram estimators (SOPHE). The method divides the data into multivariate bins and determines the shape of the data in each bin based on second order polynomials, which is an efficient computation. These calculations yield local maxima and allow joining of adjacent bins to identify clusters. The use of second order polynomials also optimally uses wide bins, such that in most cases each parameter (dimension) need only be divided into 4-8 bins, again reducing computational load. We have validated this method using defined mixtures of up to 17 fluorescent beads in 16 dimensions, correctly identifying all populations in data files of 100,000 beads in <10 s, on a standard laptop. The method also correctly clustered granulocytes, lymphocytes, including standard T, B, and NK cell subsets, and monocytes in 9-color stained peripheral blood, within seconds. SOPHE successfully clustered up to 36 subsets of memory CD4 T cells using differentiation and trafficking markers, in 14-color flow analysis, and up to 65 subpopulations of PBMC in 33-dimensional CyTOF data, showing its usefulness in discovery research. SOPHE has the potential to greatly increase efficiency of analysing complex mixtures of cells in higher dimensions. V C 2015 International Society for Advancement of Cytometry Key terms data analysis; clustering; high dimensions; complex data LYMPHOCYTES were originally viewed by light microscopy as small homogeneous round cells with minimal cytoplasm (1). Multiparameter flow cytometry, using currently available lasers, monoclonal antibodies and fluorochromes (2), has helped reveal the extremely complex heterogeneity of differentiated T and B cells with diverse immunological properties (3,4). It is possible to analyze 18 or more markers simultaneously on individual lymphocytes (2), and the development of additional labels will greatly increase this number. Study of 40 or more markers has already been achieved by the use of transition element isotopes as chelated antibody tags and mass cytometry, instead of fluorochromes (5,6).Addition of more parameters is an important goal of flow cytometric analysis of lymphocytes, because, paradoxically, instead of simply increasing complexity of the results, they can actually reveal important subsets with much greater clarity. This was originally best shown by the combination of 2 light scatter and 2 fluorescence parameters, CD45 and CD14, to clearly separate lymphocytes in blood from monocytes, granulocytes and red cell debris (7). In our experience, an important example is to first separate na€ ıve and memory T cells using CD45RA/CD45RO to measure expres-