Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH FUNCTIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

ASSESSING THE SIGNIFICANCE OF APPARENT CORRELATIONS BETWEEN RADIO AND GAMMA-RAY BLAZAR FLUXES

Whether or not a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases, which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative: although it may fail to confirm an existing intrinsic correlation in a small sample that cannot be split into many redshift bins, it will not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ROLE OF TURBULENCE AND MAGNETIC FIELDS IN SIMULATED FILAMENTARY STRUCTURE

We use numerical simulations of turbulent cluster-forming regions to study the nature of dense filamentary structures in star formation. Using four hydrodynamic and magnetohydrodynamic simulations chosen to match observations, we identify filaments in the resulting column density maps and analyze their properties. We calculate the radial column density profiles of the filaments every 0.05 Myr and fit the profiles with the modified isothermal and pressure-confined isothermal cylinder models, finding reasonable fits for either model. The filaments formed in the simulations have similar radial column density profiles to those observed. Magnetic fields provide additional pressure support to the filaments, making “puffier” filaments less prone to fragmentation than in the pure hydrodynamic case, which continue to condense at a slower rate. In the higher density simulations, the filaments grow faster through the increased importance of gravity. Not all of the filaments identified in the simulations will evolve to form stars: some expand and disperse. Given these different filament evolutionary paths, the trends in bulk filament width as a function of time, magnetic field strength, or density are weak, and all cases are reasonably consistent with the finding of a constant filament width in different star-forming regions. In the simulations, the mean FWHM lies between 0.06 and 0.26 pc for all times and initial conditions, with most lying between 0.1 to 0.15 pc; the range in FWHMs is however, larger than seen in typical Herschel analyses. Finally, the filaments display a wealth of substructure similar to the recent discovery of filament bundles in Taurus.

79 ASTRONOMY AND ASTROPHYSICS↗

Rheology of dense suspensions of shape anisotropic particles designed to show pH-sensitive anisotropic pair potentials

Here we investigate the flow properties of suspensions of dicolloidal particles composed of interpenetrating spheres where one sphere is rich in polystyrene and the second is rich in poly 2-vinyl pyridine. The synthesis method is designed to create both anisotropic shape and anisotropic interaction potentials that should lead to head to tail clustering. Furthermore, these particles are referred to as copolymer dicolloids (CDCs). The viscoelastic properties of stable and gelled suspensions of CDC particles are compared with analogs composed of homopolymer dicolloids (HDCs), having the same shape but not displaying the anisotropic attractions. After coating the particles with a nonionic surfactant to minimize van der Waals attractions, the flow properties of glassy and gelled suspensions of CDCs and HDCs are studied as a function of volume fraction, ionic strength and pH. Suspensions of HDC particles display a high kinetic arrest volume fraction (Φ g > 0.5) over a wide range of pH and ionic strength up to [I]=0.5 M, demonstrating that the particles experience repulsive or weakly attractive pair potentials. Suspensions of CDC particles behave in a similar manner at high or low pH when [I]=0.001 M, but gel at a volume fraction of Φ g < 0.3 and display anomalously large elastic moduli at and above the gel transition point for intermediate pH or for pH=9 when [I]=0.5 M. The gelation processes for the CDC particles are reversible by adjusting the solution pH. Interaction potential anisotropy is evident in the processes, during which the CDC particles yield on increasing oscillatory strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Finite-Temperature Quantum Matter with Rydberg or Molecule Synthetic Dimensions

Synthetic-dimension platforms offer unique pathways for engineering quantum matter. We compute the phase diagram of a many-body system of ultracold atoms (or polar molecules) with a set of Rydberg states (or rotational states) as a synthetic dimension, where the particles are arranged in real space in optical microtrap arrays and interact via dipole-dipole exchange interaction. Using mean-field theory, we find three ordered phases—two are localized in the synthetic dimension, predicted as zero-temperature ground states, and one is a delocalized phase. We characterize them by identifying the spontaneously broken discrete symmetries of the Hamiltonian. We also compute the phase diagram as a function of temperature and interaction strength for both signs of the interaction. For system sizes with more than six synthetic sites and attractive interactions, we find that the thermal phase transitions can be first or second order, which leads to a tricritical point on the phase boundary. Furthermore, by examining the dependence of the tricritical point and other special points of the phase boundary on the synthetic dimension size, we shed light on the physics for thermodynamically large synthetic dimension.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resonant and random excitations on the proton beam in the Large Hadron Collider for active halo control with pulsed hollow electron lenses

We present the results of numerical simulations and experimental studies about the effects of resonant and random excitations on proton losses, emittances, and beam distributions in the Large Hadron Collider (LHC). In addition to shedding light on complex nonlinear effects, these studies are applied to the design of hollow electron lenses (HEL) for active beam halo control. In the High-Luminosity Large Hadron Collider (HL-LHC), a considerable amount of energy will be stored in the beam tails. To control and clean the beam halo, the installation of two hollow electron lenses, one per beam, is being considered. In standard electron-lens operation, a proton bunch sees the same electron current at every revolution. Pulsed electron beam operation (i.e., different currents for different turns) is also considered, because it can widen the range of achievable halo removal rates. For an axially symmetric electron beam, only protons in the halo are excited. If a residual field is present at the location of the beam core, these particles are exposed to time-dependent transverse kicks and to noise. We discuss the numerical simulations and the experiments conducted in 2016 and 2017 at injection energy in the LHC. The excitation patterns were generated by the transverse feedback and damping system, which acted as a flexible source of dipole kicks. Proton beam losses, emittances, and transverse distributions were recorded as a function of excitation patterns and strengths. The resonant excitations induced rich dynamical effects and nontrivial changes of the beam distributions, which, to our knowledge, have not previously been observed and studied in this detail. We conclude with a discussion of the tolerable and achievable residual fields and proposals for further studies.

43 PARTICLE ACCELERATORS↗

Spontaneous fourfold-symmetry breaking driven by electron-lattice coupling and strong correlations in high-$T_c$ cuprates

Using dynamical-mean-field theory for clusters, we study the two-dimensional Hubbard model in which electrons are coupled with the orthorhombic lattice distortions through the modulation in the hopping matrix. Instability towards spontaneous symmetry breaking from a tetragonal symmetric phase to an orthorhombic distorted phase is examined as a function of doping and interaction strength. A very strong instability is found in the underdoped pseudogap regime when the interaction strength is large enough to yield the Mott insulating phase at half filling. The symmetry breaking accompanies the recovery of quasiparticle weights along one of the two antinodal directions, leading to the characteristic Fermi arc reconnection.We discuss the implications of our results to the fourfold symmetry breaking reported in systems where the underlying crystal does not have any structural anisotropy.

electron nematicity↗

Generalized multipolaron expansion for the spin-boson model: Environmental entanglement and the biased two-state system

We develop a systematic variational coherent state expansion for the many-body ground state of the spin-boson model, in which a quantum two-level system is coupled to a continuum of harmonic oscillators. Energetic constraints at the heart of this technique are rationalized in terms of polarons (displacements of the bath states in agreement with classical expectations) and antipolarons (counter- displacements due to quantum tunneling effects). We present a comprehensive study of the ground state two-level system population and coherence as a function of tunneling amplitude, dissipation strength, and bias (akin to asymmetry of the double well potential defining the two-state system). The entanglement among the different environmental modes is investigated by looking at spectroscopic signatures of the bipartite entanglement entropy between a given environmental mode and all the other modes. We observe a drastic change in behavior of this entropy for increasing dissipation, indicative of the entangled nature of the environmental states. In addition, the entropy spreads over a large energy range at strong dissipation, a testimony to the wide entanglement window characterizing the underlying Kondo state. Lastly, comparisons to accurate numerical renormalization group calculations and to the exact Bethe Ansatz solution of the model demonstrate the rapid convergence of our variationally-optimized multi-polaron expansion, suggesting that it should also be a useful tool for dissipative models of greater complexity, as relevant for numerous systems of interest in quantum physics and chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transverse ac-driven and geometric ratchet effects for vortices in conformal crystal pinning arrays

A conformal pinning array is created by taking a conformal transformation of a uniform hexagonal lattice to create a structure in which the sixfold ordering of the original lattice is preserved but which has a spatial gradient in the pinning site density. With a series of conformal arrays it is possible to create asymmetric substrates, and it was previously shown that when an ac drive is applied parallel to the asymmetry direction, a pronounced ratchet effect occurs with a net dc flow of vortices in the same direction as the ac drive. Here, in this article, we show that when the ac drive is applied perpendicular to the substrate asymmetry direction, it is possible to realize a transverse ratchet effect where a net dc flow of vortices is generated perpendicular to the ac drive. The conformal transverse ratchet effect is distinct from previous versions of transverse ratchets in that it occurs due to the generation of non-Gaussian transverse vortex velocity fluctuations by the plastic motion of vortices, so that the system behaves as a noise correlation ratchet. The transverse ratchet effect is much more pronounced in the conformal arrays than in random gradient arrays and is absent in square gradient arrays due the different nature of the vortex flow in each geometry. We show that a series of reversals can occur in the transverse ratchet effect due to changes in the vortex flow across the pinning gradient as a function of vortex filling, pinning strength, and ac amplitude. We also consider the case where a dc drive applied perpendicular to the substrate asymmetry direction generates a net flow of vortices perpendicular to the dc drive, producing what is known as a geometric or drift ratchet that again arises due to non-Gaussian dynamically generated fluctuations. The drift ratchet is more efficient than the ac driven ratchet and also exhibits a series of reversals for varied parameters. Lastly, our results should be general to a wide class of systems undergoing nonequilibrium dynamics on conformal substrates, such as colloidal particles on optical traps.

36 MATERIALS SCIENCE↗

Numerical results on the short-range spin correlation functions in the ground state of the two-dimensional Hubbard model

Optical lattice experiments with ultracold fermion atoms and quantum gas microscopy have recently realized direct measurements of magnetic correlations at the site-resolved level. We calculate the short-range spin-correlation functions in the ground state of the two-dimensional repulsive Hubbard model with the auxiliary-field quantum Monte Carlo (AFQMC) method. The results are numerically exact at half filling where the fermion sign problem is absent. Away from half filling, we employ the constrained path AFQMC approach to eliminate the exponential computational scaling from the sign problem. The constraint employs unrestricted Hartree-Fock trial wave functions with an effective interaction strength U, which is optimized self-consistently within AFQMC. Large supercells are studied, with twist averaged boundary conditions as needed, to reach the thermodynamic limit. We find that the nearest-neighbor spin correlation always increases with the interaction strength U , contrary to the finite-temperature behavior where a maximum is reached at a finite U value. We also observe a change of sign in the next-nearest-neighbor spin correlation with increasing density, which is a consequence of the buildup of the long-range antiferromagnetic correlation. Finally, we expect the results presented in this paper to serve as a benchmark as lower temperatures are reached in ultracold atom experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calorimetric Measurements of Magnetic-Field-Induced Inhomogeneous Superconductivity Above the Paramagnetic Limit

We report the first magnetocaloric and calorimetric observations of a magnetic-field-induced phase transition within a superconducting state to the long-sought exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state, first predicted over 50 years ago. Through the combination of bulk thermodynamic calorimetric and magnetocaloric measurements in the organic superconductor. kappa-(BEDT-TTF)(2) Cu(NCS)(2) as a function of temperature, magnetic field strength, and magnetic field orientation, we establish for the first time that this field-induced first-order phase transition at the paramagnetic limit Hp is a transition to a higher-entropy superconducting phase, uniquely characteristic of the FFLO state. We also establish that this high-field superconducting state displays the bulk paramagnetic ordering of spin domains required of the FFLO state. These results rule out the alternate possibility of spin-density wave ordering in the high-field superconducting phase. The phase diagram determined from our measurements-including the observation of a phase transition into the FFLO phase at Hp-is in good agreement with recent NMR results and our own earlier tunnel-diode magnetic penetration depth experiments but is in disagreement with the only previous calorimetric report.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiferroicity and Topology in Twisted Transition Metal Dichalcogenides

Van der Waals heterostructures have recently emerged as an exciting platform for investigating the effects of strong electronic correlations, including various forms of magnetic or electrical orders. Here, we perform an unbiased exact diagonalization study of the effects of interactions on topological flat bands of twisted transition metal dichalcogenides (TMDs) at odd integer fillings. For hole-filling ν h = 1 , we find that the Chern insulator phase, expected from interaction-induced spin-valley polarization of the bare bands, is quite fragile, and gives way to spontaneous multiferroic order—coexisting ferroelectricity and ferromagnetism, in the presence of long-range Coulomb repulsion. We provide a simple real-space picture to understand the phase diagram as a function of interaction range and strength. Our findings establish twisted TMDs as a novel and highly tunable platform for multiferroicity, and we outline a potential route towards electrical control of magnetism in the multiferroic phase. Published by the American Physical Society 2024

Abouelkomsan, Ahmed (ORCID:0000000247873357)↗

Solutes that reduce yield strength anisotropies in magnesium from first principles

Using Labusch-type solid solution strengthening models parameterized with DFT-computed solute-dislocation interaction energies, we perform a computational search for 63 solutes across the periodic table to find those that lower anisotropy ratios (non-basal to basal CRSS) of magnesium potentially increasing its ductility per the von Mises criterion. For this purpose, we compute changes in strength for solutes as a function of composition and temperature, and compute anisotropy ratios for solutes that include both rare earth and non-rare earth elements. Here we specifically focus on solute-dislocation interaction energies in the following DFT-optimized dislocations as representative of three non-basal plastic deformation modes: $\langle c + a \rangle$ edge, (10$\bar1$2) tension twinning edge, and the (10$\bar{1}$1) compression twinning edge. We find that solute-induced changes in non-basal deformation modes can be approximated using a second-order polynomial in the size misfit of the solutes, which permits rapid screening of solutes. Our approach to identify solutes known to improve strengthening incorporates solute solubility, and suggests other solutes that not have been previously explored for strengthening. The 8 rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios at their optimal concentrations, are: Gd, Tb, Dy, Nd, Ho, Er, Tm, and Yb. The 12 non-rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios, are: Y, Mn, Sc, Pb, Ca, Ag, Bi, Tl, Zn, Li, Ga, and Al. Of these, Gd, Nd, Er, Yb, Y, Mn, Ca, Zn, Li, and Al are used in commercial Mg alloys.

36 MATERIALS SCIENCE↗

SU-F-J-50: Study On the Magnetic Field Effect On the Exradin W1 Plastic Scintillation Detector

Purpose: To study the response of the Exradin W1 plastic scintillator detector to a 6 MV photon field in the presence of a strong magnetic field (B). Methods: An Exradin W1 scintillator detector coupled to a SuperMax two-channel electrometer, both manufactured by Standard Imaging, Inc., was first calibrated in a Co-60 beam. The Cerenkov Light Ratio (CLR) was obtained following the procedure recommended by the manufacturer. Subtracting signal in channel 2 multiplied by CLR from the signal in channel 1 should lead to a Cerenkov-free signal. The W1 scintillator was placed in a plastic phantom inside a dipole electromagnet (GMW Associates) that could produce a strong B field, and irradiated using a 6 MV beam from an Elekta Versa HD LINAC. Signals from both channels of the W1 scintillator were acquired as a function of B (0 - 1.5 T). Results: The signals from both channels increased as a function of the B field strength. At 1.5 T, channel 1 increased by 11% compared to the baseline (B=0 T), while channel 2 increased by 22%. Applying the recommended Cerenkov correction led to a 2% difference between dose measurement with and without a magnetic field. The values between B=0.3 T and B=1.5 T remained constant. Conclusion: Signals from the Exradin W1 plastic scintillation detector increased as the B field increased. This increase mainly comes from a change in the amount of Cerenkov light coupled within the optical fiber. Removing the Cerenkov component following the procedure recommended by the manufacturer showed to be an effective way to measure dose accurately in strong magnetic fields. The cause for the residual 2% difference is currently under investigation. We acknowledge research support from Elekta AB.

60 APPLIED LIFE SCIENCES↗

The effect of temperature on the nonlinear dynamics of charge in a semiconductor superlattice in the presence of a magnetic field

The space-time dynamics of electron domains in a semiconductor superlattice is studied in a tilted magnetic field with regard to the effect of temperature. It is shown that an increase in temperature substantially changes the space-time dynamics of the system. This leads to a decrease in the frequency and amplitude of oscillations of a current flowing through the semiconductor superlattice. The quenching of oscillations is observed, which is attributed to the change in the drift velocity as a function of electric-field strength under the variation of temperature.

36 MATERIALS SCIENCE↗

Analysis of Monolith Cores from an Engineering Scale Demonstration of a Prospective Cast Stone Process

The primary disposition path of Low Activity Waste (LAW) at the DOE Hanford Site is vitrification. A cementitious waste form is one of the alternatives being considered for the supplemental immobilization of the LAW that will not be treated by the primary vitrification facility. Washington River Protection Solutions (WRPS) has been directed to generate and collect data on cementitious or pozzolanic waste forms such as Cast Stone. This report documents the coring and leach testing of monolithic samples cored from an engineering-scale demonstration (ES Demo) with non-radioactive simulants. The ES Demo was performed at SRNL in October of 2013 using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft. diameter x 3.25 ft. high container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average LAW composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. In 2014 core samples originally obtained approximately six months after filling the ES Demo were tested along with bench scale molded samples that were collected during the original pour. A latter set of core samples were obtained in late March of 2015, eighteen months after completion of the original ES Demo. Core samples were obtained using a 2” diameter x 11” long coring bit. The ES Demo was sampled in three different regions consisting of an outer ring, a middle ring and an inner core zone. Cores from these three lateral zones were further segregated into upper, middle and lower vertical segments. Monolithic core samples were tested using the Environmental Protection Agency (EPA) Method 1315, which is designed to provide mass transfer rates (release rates) of inorganic analytes contained in monolithic material under diffusion controlled release conditions as a function of leaching time. Compressive strength measurements and drying tests were also performed on the 2015 samples. Leachability indices reported are based on analyte concentrations determined from dissolution of the dried samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Final Technical Report: Targeting DOE-Relevant Ions with Supramolecular Strategies, DE-SC0010555

The effectiveness of three popular supramolecular strategies to selectively target negatively charged ions (anions) was evaluated. Ions of interest included oxo anions, particularly sulfate, that hamper nuclear waste remediation. Three objectives were pursued using a simple building block strategies and by strategically placing anion-binding sites at appropriate positions on organic host molecules. The goal of the first objective was to assess the influence of secondary, tertiary and quaternized amines on binding tetrahedral anions using mixed amide/amine macrocyclic and urea/amine hosts containing aromatic or heteroaromatic spacers. Objective 2 focused on the design of ion pair hosts, using mixed macrocyclic anion hosts joined through polyether linkages. Objective 3 was to explore the synthesis of new metal-linked extended macrocyclic frameworks to leverage anion binding. Key findings were that smaller 24-membered macrocycles provided the most complementary binding for sulfate ion and mixed urea/amine chelates showed enhanced binding over amide corollaries in addition to being highly selective for SO 4 2- in the presence of small quantities of water. In addition to obtaining prototype metal-linked macrocyclic anion hosts, a new dipincer ligand was designed that can be used to link macrocyclic or other supramolecular hosts in extended frameworks. When the tetraamide-based pincers are bound to two metal ions, an interesting phenomenon occurs. Upon deprotonation of the amides, two new protons appear between adjacent carbonyl pairs on the ligand, which may modify the chemistry, and metal-metal interactions in the complexes. Gel formation occurred for some of these extended hosts, and the physical properties are currently under investigation. The new tetracarboxamide-based pincers can also provide basic frameworks for double macrocycles capable of binding ion pairs as well as for binding metal ions and exploring intermetallic interactions through the pyrazine π system. Additionally appendages capable of influencing solvation effects can be introduced, and a number of other potential applications can be realized in areas such as soft materials chemistry, catalysis, sensing, and proton switches, the latter for binding and release of targeted guests. These findings provide a better foundation for understanding the selective binding of anions by targeted placement of hydrogen binding sites, and the strengths and weaknesses of various functional groups, that will allow for more the design of more effective anion sequestering agents. Our design strategy also used simple, cost-effective building blocks for host synthesis to allow for scale-up should real-world applications be forthcoming.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Resonance Parameter Evaluation of n + 88 Sr reactions for ENDF/B-VIII.1 Library

An evaluation of n+ 88 Sr reactions up to 950 keV was compiled from a publicly available set of resonance parameters. Starting from the repository of the ENDF/B-VIII.0 library, an evaluated nuclear data file was assembled and submitted for inclusion to the Evaluated Nuclear Data File (ENDF)/B-VIII.1 release. This evaluation features an external function representation defined by pole strengths and R ∞ parameters over the entire resolved resonance range. This work improves the modeling of the R-matrix parameterization from multi-Level Breit-Wigner resonance parameters to Reich-Moore formalism, and it also extends thrice the upper energy range of the current ENDF/B-VIII.0 resolved resonance evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗