This paper presents a new method for calculating accurate masses of isotopic peaks. It is based on breaking the calculation into a binary series of calculations. The molecule is built up by a series of such calculations. At each step the accurate masses are calculated as a probability weighted sum of the masses of the contributing peaks. The method is computationally efficient and accurate for both mass and relative abundance. T his paper addresses the problem of calculating the accurate masses of isotopic peaks in an isotopic distribution of a compound of known chemical formula. For the purposes of this paper an "isotopic peak" consists of all molecules having the same number of nucleons, regardless of the isotopic fine structure of the peak. For example, the compound CO has an isotopic peak containing 29 nucleons. Although this peak has an isotopic fine structure with contributions from both 12 C 17 O and 13 C 16 O, each of which has a slightly different mass, for the purposes of the present paper these will not be considered as two separate peaks but as a single isotopic peak with an accurate mass of 28.998297 Da.There are several methods for calculating masses of isotopic peaks. The accuracy and speed of an isotopic calculation both depend strongly on the method used. Kubinyi [1] and Rockwood and Van Orden [2] have described fast methods that also produce semi-accurate masses. The method of reference [2] is extremely fast and generally produces results within a few millimass units [2,3]. Although semi-accurate masses are useful for many applications, in some cases being more accurate than experimental measurements of masses [3], for other applications accurate masses would be preferable.Accurate masses of isotopic peaks can be calculated by several methods. Polynomial-based methods work well for low and medium molecular weight compounds, but become computationally inefficient at high molecular weights. Attempts to accelerate the polynomial-based methods by applying "pruning", which is the omission of small terms from the calculation, involve a tradeoff between speed and accuracy [4,5].In a second method, one might start with an algorithm that calculates profile-mode isotopic distributions centered on the accurate masses [4]. Numerical values of the accurate masses of the isotopic peaks would then be extracted by numerical quadrature or by peak fitting of the nominal isotopic peaks. However, considerable computational effort would still be involved in such a method.A third method is to perform a series of stepwise convolutions, building the full molecule by the stepwise addition of atoms to an accumulated "super atom". Masses of the individual isotopic peaks are calculated at each convolution step as a probability weighted sum of the contributions to each peak [3]. This method is accurate, but computational efficiency is not optimal.Recently, a new algorithm was presented for the calculation of accurate masses of isotopic peaks, and application examples were discussed [6]. This method combines very high compu...