Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomically precise”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Systematic error analysis in the absolute hydrogen gas jet polarimeter at RHIC

The Polarized Atomic Hydrogen Gas Jet Target polarimeter (HJET) is operated at the Relativistic Heavy Ion Collider (RHIC) since 2004 to measure the absolute polarization of each colliding proton beam. Polarimeter detectors and data acquisition were upgraded in 2015 to increase the solid angle, energy range, and to improve the energy and time resolution. These upgrades along with an improved beam intensity and polarization allowed us to greatly reduce the statistical and systematic errors for the proton polarization measurements in RHIC Runs 15 (E beam = 100 GeV) and 17 (255 GeV). For a typical 8 hour RHIC store, the measured proton beam average polarization was about P beam ~ 55 ± 2.0 stat ± 0.3 syst %. The elastic pp analyzing power, AN(t)~0.04, was determined with a precision of about |δA N (t)|~0.0002 in the momentum transfer squared range 0.001<-t<0.020 GeV 2 . In this paper we present a detailed systematic error analysis of the polarization measurements at HJET. Perspectives of using the HJET based absolute polarimeters in the future Electron Ion Collider (EIC) will be also discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Exploiting 20 Ne Isotopes for Precision Characterizations of Collectivity in Small Systems

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174⁢(8) stat ⁢(31) syst for NLEFT and 1.139⁢(6) stat ⁢(39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Copper Clusters Containing Hydrides in Trigonal Pyramidal Geometry

Structurally precise copper hydrides [Cu 11 H 2 {S 2 P(O i Pr) 2 } 6 (C≡CR) 3 ], R = Ph (1), C 6 H 4 F (2), and C 6 H 4 OMe (3), were first synthesized from the polyhydrido copper cluster [Cu 20 H 11 {S 2 P(O i Pr) 2 } 9 ] with nine equivalents of terminal alkynes. Later, their isolated yields were significantly improved by direct synthesis from [Cu(CH 3 CN) 4 ](PF 6 ), [NH 4 ][S 2 P(O i Pr) 2 ], NaBH 4 , and alkynes along with NEt 3 in THF. 1, 2, and 3 were fully characterized by single-crystal X-ray diffraction, ESI-MS, and multinuclear NMR spectroscopy. All three clustershave 11 copper atoms, adopting 3,3,4,4,4-pentacapped trigonal prismatic geometry, with two hydrides inside the Cu 11 cage, the position of which was ascertained by a single-crystal neutron diffraction structure of cluster 1 co-crystallized with a [Cu 7 (H){S 2 P(O i Pr) 2 } 6 ] (4) cluster. Six dithiophosphate and three alkynyl ligands stabilize the Cu 11 H 2 core in which the two hydrides adopt a trigonal pyramidal coordination mode. This coordination mode is so far unprecedented for hydride. The 1 H NMR resonance frequency of the two hydrides appears at 4.8 ppm, a value further confirmed by 2 H NMR spectroscopy for their deuteride derivatives [Cu 11 (D) 2 {S 2 P(O i Pr) 2 } 6 (C≡CR) 3 ]. A DFT investigation allows understanding the bonding within this new type of copper(I) hydrides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗
DOE OSTI2021

On-chip optical levitation with a metalens in vacuum

Optical levitation of dielectric particles in vacuum is a powerful technique for precision measurements, testing fundamental physics, and quantum information science. Conventional optical tweezers require bulky optical components for trapping and detection. Here, we design and fabricate an ultrathin dielectric metalens with a high numerical aperture of 0.88 at 1064 nm in vacuum. It consists of 500-nm-thick silicon nano-antennas, which are compatible with an ultrahigh vacuum. We demonstrate optical levitation of nanoparticles in vacuum with a single metalens. The trapping frequency can be tuned by changing the laser power and polarization. We also transfer a levitated nanoparticle between two separated optical tweezers. Optical levitation with an ultrathin metalens in vacuum provides opportunities for a wide range of applications including on-chip sensing. Such metalenses will also be useful for trapping ultracold atoms and molecules.

Shen, Kunhong↗
DOE OSTI2025

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Engineering of Exciton Spatial Distribution in CdS Nanoplatelets

Zinc-blende CdS nanoplatelets with atomically flat and very large {100} basal planes terminated solely by one type of element (either Cd or S atoms) are synthesized. Optical spectroscopy, X-ray diffraction, X-ray absorption, and transmission electron microscopy confirm that the surface structures of newly developed S-terminated CdS nanoplatelets are at least as well-defined as the original Cd-terminated nanoplatelets. Band gaps of the nanoplatelets are found to depend on not only the quantum-confined dimension (thickness) but also the nature of the surface termination. The facet structure dictates the packing of the ligands (carboxylate for Cd-terminated nanoplatelets and alkyl for S-terminated nanoplatelets), which causes a difference in the lattice strain and significantly affects the optical spectral width. Experimental and theoretical results reveal that engineering the exciton spatial distribution by the tailored synthesis of semiconductor nanocrystals with a precisely controlled surface structure is fully possible, which should open a new door for delivering the long-promised potential of semiconductor nanocrystals.

Cadmium sulfide↗
DOE OSTI2022

Formation and detection of metastable formic acid in a supersonic expansion: High resolution infrared spectroscopy of the jet-cooled cis -HCOOH conformer

High-resolution direct absorption infrared spectra of metastable cis-formic acid (HCOOH) trapped in a cis-well resonance behind a 15 kcal/mol barrier are reported for the first time, with the energetically unstable conformer produced in a supersonic slit plasma expansion of trans-formic acid/H 2 mixtures. We present a detailed high-resolution rovibrational analysis for cis-formic acid species in the OH stretch (ν 1 ) fundamental, providing first precision vibrational band origin, rotational constants, and term values, which in conjunction with ab initio calculations at the couple-cluster with single, double, and perturbative triple [CCSD(T)]/ANOn (n = 0, 1, 2) level support the experimental assignments and establish critical points on the potential energy surface for internal rotor trans-to-cis isomerization. Relative intensities for a- and b-type transitions observed in the spectra permit the transition dipole moment components to be determined in the body fixed frame and prove to be in good agreement with ab initio CCSD(T) theoretical estimates but in poor agreement with simple bond-dipole predictions. As a result, the observed signal dependence on H 2 in the discharge suggests the presence of a novel H atom radical chemical mechanism for strongly endothermic “up-hill” internal rotor isomerization between trans- and cis-formic acid conformers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Mechanical Detection of Nuclear Decays

We report the detection of individual nuclear alpha decays through the mechanical recoil of the entire micron-sized particle in which the decaying nuclei are embedded. Momentum conservation ensures that such measurements are sensitive to any particles emitted in the decay, including neutral particles that may otherwise evade detection with existing techniques. Detection of the minuscule recoil of an object more than 10^12 times more massive than the emitted particles is made possible by recently developed techniques in levitated optomechanics, which enable high-precision optical control and measurement of the mechanical motion of optically trapped particles. Observation of a change in the net charge of the particle coincident with the recoil allows decays to be identified with background levels at the micro-Becquerel level. Here, the techniques developed here may find use in fields ranging from nuclear forensics to dark matter and neutrino physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Local microwave sensing of excitons and their electrical environment

Excitons in atomically thin transition metal dichalcogenides (TMDs) possess intriguing optical properties, drawing interest for both technology and fundamental research. However, as the demands for nanodevice applications and the exploration of fundamental physics push toward smaller, subwavelength scales, studying them locally is challenging. In this work, we introduce a cryogenic scanning probe photoelectrical sensing technique, termed exciton-resonant microwave impedance microscopy (ER-MIM), to measure the excitonic responses in a monolayer MoSe 2 device at 1.5K. From the microwave signal changes, we identify exciton polarons and their Rydberg states. Building on these observations, we systemically reveal the local and nonlocal effects of carrier density, inhomogeneous electric fields, as well as dielectric screening on excitons, beyond the reach of conventional probes. By further integrating deep learning techniques, we precisely extracted the electrical parameters surrounding excitons, demonstrating a quantified, exciton-assisted nanoscale electrometry. Our results provide new insight into exciton-environment interactions, establish ER-MIM as a powerful optoelectronic sensing platform, and open avenues for exciton-based quantum control and device technologies.

scanning probe microscopy↗
NASA NTRS1994

Time maintenance system for the BMDO MSX spacecraft

The Johns Hopkins University Applied Physics Laboratory (APL) is responsible for designing and implementing a clock maintenance system for the Ballistic Missile Defense Organizations (BMDO) Midcourse Space Experiment (MSX) spacecraft. The MSX spacecraft has an on-board clock that will be used to control execution of time-dependent commands and to time tag all science and housekeeping data received from the spacecraft. MSX mission objectives have dictated that this spacecraft time, UTC(MSX), maintain a required accuracy with respect to UTC(USNO) of +/- 10 ms with a +/- 1 ms desired accuracy. APL's atomic time standards and the downlinked spacecraft time were used to develop a time maintenance system that will estimate the current MSX clock time offset during an APL pass and make estimates of the clock's drift and aging using the offset estimates from many passes. Using this information, the clock's accuracy will be maintained by uplinking periodic clock correction commands. The resulting time maintenance system is a combination of offset measurement, command/telemetry, and mission planning hardware and computing assets. All assets provide necessary inputs for deciding when corrections to the MSX spacecraft clock must be made to maintain its required accuracy without inhibiting other mission objectives. The MSX time maintenance system is described as a whole and the clock offset measurement subsystem, a unique combination of precision time maintenance and measurement hardware controlled by a Macintosh computer, is detailed. Simulations show that the system estimates the MSX clock offset to less than+/- 33 microseconds.

Hermes, Martin J.↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2026

Instantaneous difference frequency locking observed during toroidicity-induced Alfvén eigenmode coupling in the DIII-D tokamak

In magnetic confinement fusion, toroidicity-induced Alfvén eigenmodes (TAEs) are well-studied, weakly stable solutions of the linearized ideal magnetohydrodynamic equations. Driven unstable by suprathermal populations of energetic particles, TAE pose a key vulnerability to the confinement of high-energy alpha particles generated by fusion reactions. Hence, it is paramount to understand TAE dynamics if a working reactor is to be realized. In this work, we detect and characterize signatures of nonstationary nonlinear coupling between TAE using a novel, time-resolved bispectral analysis; results are supported by analytic signal of band-passed data. Crucially, a stationary phase relationship between two TAE and a nascent low frequency fluctuation is observed precisely when the triple product of magnetic fluctuation amplitudes is enhanced. Local mode number and frequency spectrum, gleaned from beam-emission spectroscopy, corroborates simultaneous satisfaction of nonlinear matching conditions, and provides a tool to identify theorized pathways of energy transfer, e.g. TAE parametric instability.

bispectral analysis↗
DOE OSTI2022

Development of Single Photon Sources in GaN

The recent discovery of bright, room-temperature, single photon emitters in GaN leads to an appealing alternative to diamond best single photon emitters given the widespread use and technological maturity of III-nitrides for optoelectronics (e.g. blue LEDs, lasers) and high-speed, high-power electronics. This discovery opens the door to on-chip and on-demand single photon sources integrated with detectors and electronics. Currently, little is known about the underlying defect structure nor is there a sense of how such an emitter might be controllably created. A detailed understanding of the origin of the SPEs in GaN and a path to deterministically introduce them is required. In this project, we develop new experimental capabilities to then investigate single photon emission from GaN nanowires and both GAN and AlN wafers. We ion implant our wafers with the ion implanted with our focused ion beam nanoimplantation capabilities at Sandia, to go beyond typical broad beam implantation and create single photon emitting defects with nanometer precision. We've created light emitting sources using Li + and He + , but single photon emission has yet to be demonstrated. In parallel, we calculate the energy levels of defects and transition metal substitutions in GaN to gain a better understanding of the sources of single photon emission in GaN and AlN. The combined experimental and theoretical capabilities developed throughout this project will enable further investigation into the origins of single photon emission from defects in GaN, AlN, and other wide bandgap semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Forward silicon vertex/tracking detector design and R$\&$D for the future Electron-Ion Collider

The proposed high-luminosity high-energy Electron-Ion Collider (EIC) will provide a clean environment to precisely study several fundamental questions in the fields of high-energy and nuclear physics . A low material budget and high granularity silicon vertex/tracking detector is critical to carry out a series of hadron and jet measurements at the future EIC especially for the heavy flavor product reconstruction or tagging. The conceptual design of a proposed forward silicon tracking detector with the pseudorapidity coverage from 1.2 to 3.5 has been developed in integration with different magnet options and the other EIC detector sub-systems. The tracking performance of this detector enables precise heavy flavor hadron and jet measurements in the hadron beam going direction. The detector R&D for the proposed silicon technology candidates: Low Gain Avalanche Diode (LGAD) and radiation hard depleted Monolithic Active Pixel Sensor (MALTA), which can provide good spatial and timing resolutions, is underway. Bench test results of the LGAD and MALTA prototype sensors will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

U-Pu and Ba-Cs isotopic measurements on Trinitite by laser ablation sampling on the Neoma MC-ICP-MS

In this study we present the results of combined U-Pu and Ba-Cs isotope measurements obtained by laser ablation (LA) sampling of two glassy debris fragments (‘Trinitite’) from the world's first atomic bomb detonation conducted in New Mexico on July 16, 1945. Our primary goal in conducting these measurements was to understand whether examination of the U-Pu and Ba-Cs systematics by direct sampling (e.g. without any chemical separation or purification prior to isotope ratio measurement) could yield meaningful information that would differentiate the Trinitite fragments from glassy material lacking a nuclear fission signature. These measurements were conducted on a ThermoFisher Scientific Neoma multi collector – inductively coupled plasma – mass spectrometer (MC-ICP-MS), which is a relatively new MC-ICP-MS platform, so we also examine the behavior of these isotope systems in standards sampled in solution and via LA. Unsurprisingly, the measurements made on purified solutions of the U, Pu, and Ba isotopic standards produce high precision isotope ratios. Furthermore, this extends to the U-Pu measurements made by LA sampling, with the expected degradation in precision and accuracy related to matrix effects and signal intensity fluctuation. However, the Ba-Cs data acquired by LA is of low precision across all of the matrices examined and bears evidence of complex mass fractionation that will require further investigation to resolve. In total, our results indicate that the observed U-Pu isotope data are of sufficient quality to accurately constrain the U and Pu isotopic composition of glass containing sub-ppm levels of these elements which in turn could be used to differentiate glass containing anthropogenic fission products from natural glass whereas the Ba-Cs LA data cannot be used for this purpose until further methodological refinement is performed.

Ba-Cs↗
NASA NTRS1986

Theoretical analysis of isotope effects on ozone formation in oxygen photochemistry

In situ measurements of stratospheric ozone and laboratory studies of ozone production in electric discharge through oxygen have shown previously that ozone containing heavy isotopes of oxygen (O-17, O-18) may be formed preferentially. In order to assess the relevance of thee latter experiment to the stratospheric measurements, detailed understanding of the effect of isotopic substitution on the O3 formation reaction O + O2 + M yields O3 + M and on the O atom exchange reaction O + O2 + O yields O2 + O is necessary. In this work, an estimate of the effect of isotopic substitution on the recombination rate is made by us of approximate dynamical theories and statistical mechanics. The results indicate the possibility of isotope effects on the O + O2 recombination rate of the order of several percent at stratospheric temperatures. In general, recombination reactions involving heavy (mass 49, 50) O3 formation are found to be slower than the reaction leading to normal (mass 48) O3 formation. The calculated isotope effects are sufficiently small that the uncertainties in the model input and the approximations in the dynamical theories will probably make the quantitative nature of these results subject to considerable uncertainty. This isotope effect should not be observable in the atmosphere given the precision of the current measurements but may be crucial in the understanding of the laboratory experiments, where observed enhancements are only of the order of several percent. Possible reasons for this discrepancy between the observed enhancement and predicted depletion are presented.

Kaye, J. A.↗