skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Fast Exact Computation of the k Most Abundant Isotope Peaks with Layer-Ordered Heaps
Computation of the isotopic distribution of compounds is crucial to applications of mass spectrometry, particularly as machine precision continues to improve. In the past decade, several tools have been created for doing so. In this paper we present a novel algorithm for calculating either the most abundant k isotopologue peaks of a compound or the minimal set of isotopologue peaks which have a combined total abundance of at least p. The algorithm uses Serang’s optimal method of selection on Cartesian products. The method is significantly faster than the state-of-the-art on large compounds (e.g., Titin protein) and on compounds whose elements have many isotopes (e.g., palladium alloys).  more » « less
Award ID(s):
1845465
PAR ID:
10232266
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Analytical chemistry
Volume:
92
Issue:
15
ISSN:
0003-2700
Page Range / eLocation ID:
10613–10619
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In the present work, the advantages of ESI-TIMS-FT-ICR MS to address the isomeric content of dissolved organic matter are studied. While the MS spectra allowed the observation of a high number of peaks ( e.g. , PAN-L: 5004 and PAN-S: 4660), over 4× features were observed in the IMS-MS domain ( e.g. , PAN-L: 22 015 and PAN-S: 20 954). Assuming a total general formula of C x H y N 0–3 O 0–19 S 0–1 , 3066 and 2830 chemical assignments were made in a single infusion experiment for PAN-L and PAN-S, respectively. Most of the identified chemical compounds (∼80%) corresponded to highly conjugated oxygen compounds (O 1 –O 20 ). ESI-TIMS-FT-ICR MS provided a lower estimate of the number of structural and conformational isomers ( e.g. , an average of 6–10 isomers per chemical formula were observed). Moreover, ESI-q-FT-ICR MS/MS at the level of nominal mass ( i.e. , 1 Da isolation) allowed for further estimation of the number of isomers based on unique fragmentation patterns and core fragments; the later suggested that multiple structural isomers could have very closely related CCS. These studies demonstrate the need for ultrahigh resolution TIMS mobility scan functions ( e.g. , R = 200–500) in addition to tandem MS/MS isolation strategies. 
    more » « less
  2. Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g, the coordination number or chemical environment), and symmetry (e.g, Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN8, with a two-dimensional polymeric nitrogen sublattice, brookite TiO2, and a high pressure BaH4 phase, which was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1̅ WC3N6, was found to possess an energy that is lower than that of two phases proposed in a recent computational study. The method presented here could help in identifying the structures of compounds that have already been synthesized, and in developing new synthesis targets with desired properties. 
    more » « less
  3. In recent years, the success of serial femtosecond crystallography and the paucity of beamtime at X-ray free-electron lasers have motivated the development of serial microcrystallography experiments at storage-ring synchrotron sources. However, especially at storage-ring sources, if a crystal is too small it will have suffered significant radiation damage before diffracting a sufficient number of X-rays into Bragg peaks for peak-indexing software to determine the crystal orientation. As a consequence, the data frames of small crystals often cannot be indexed and are discarded. Introduced here is a method based on the expand–maximize–compress (EMC) algorithm to solve protein structures, specifically from data frames for which indexing methods fail because too few X-rays are diffracted into Bragg peaks. The method is demonstrated on a real serial microcrystallography data set whose signals are too weak to be indexed by conventional methods. In spite of the daunting background scatter from the sample-delivery medium, it was still possible to solve the protein structure at 2.1 Å resolution. The ability of the EMC algorithm to analyze weak data frames will help to reduce sample consumption. It will also allow serial microcrystallography to be performed with crystals that are otherwise too small to be feasibly analyzed at storage-ring sources. 
    more » « less
  4. The measurement of over 24,000 new rotational transitions of chloropyrazine (C4H3N2Cl, la = 1.55 D, lb = 0.14 D) in the 130–375 GHz frequency region enabled the assignments of the rotational spectra for the ground vibrational states of the [35Cl]- and [37Cl]-isotopologues and a combined total of 22 low-energy excited vibrational states. These vibrational states have been fit to sextic A-reduced Hamiltonians with low error (<0.05 MHz). These data allow for precise determination of the vibration-rotation interaction constants for the six lowest-energy fundamental vibrational modes (m24, m17, m23, m16, m22, and m15) of the [35Cl]-isotopologue and the four lowest-energy fundamental vibrational modes (m24, m17, m23, and m16) of the [37Cl]-isotopologue. In addition, many combination and overtone states were observed, bringing the total number of excited states to seventeen for the [35Cl]-isotopologue and five for the [37Cl]-isotopologue. The spectroscopic constants for the ground state, m24, and m17 for each isotopologue are very well determined (Nlines > 1400 for each least-squares fit). The vibration-rotation interaction constants experimentally determined for all observed fundamentals are generally in quite close agreement with their predicted (B3LYP/6-311+G(2d,p)) values. 
    more » « less
  5. Energies, wavefunctions and lifetimes of vibrational resonances were computed for 18O-enriched isotopologue 50O_3=16O16O18O and 16O18O16O of the ozone molecule using hyperspherical coordinates and the method employing complex absorbing potential. 
    more » « less