Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic systems”

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,027 records · Page 57

Comparison between Hayabusa 2 Spectral Measurements and Simulations

Recent Hayabusa 2 spectral measurements are compared with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the cou- pled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator based on available published information. The comparison between the simulations and measurements focuses on three atomic nitrogen lines, two atomic oxygen lines, one atomic hydrogen line, and the CN Violet band system. Parametric uncertainties are evaluated for these features using a recently developed capability in LAURA/HARA. These uncertainties capture the impact of flowfield and radiation modeling uncertainties on the observed radiation simulations, and include the impact of coupled radiation and ablation product emission. The resulting parametric uncertainties range from 50 to over 100% for the atomic lines, and 80% for CN Violet. Uncertainties in flowfield kinetics and non-Boltzmann modeling provide the dominant contribution to these uncertainties. For the atomic nitrogen and oxygen lines, com- parisons between the nominal simulations and measurements show agreement well within these simulation uncertainty bounds, with the measurements within 20% of the nominal simulation over most of the trajectory, and within 10% at peak emission. This excellent agreement is un- precedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield and radiation modeling. For the ablation product emission from the 656 nm atomic hydrogen line, the agreement between measurements and simulations is within ±20% for all trajectory points except one outlier. This good comparison is surprising considering this species depends on the ablation rate and is shown to have strong wake contribution, which both contribute to the large ≈±70% parametric uncertainty. For the ablation product emission from CN Violet, the simulations over-predict the measurements by over a factor of two early in the trajectory, with the comparison improving later in the trajectory. This disagreement for CN Violet is consistent with the ≈±80% parametric uncertainty bounds evaluated for the simulations. To utilize the excellent agreement between the atomic line simulations and measurements in the assessment of the radiative heating margin for a future flight vehicle, the flowfield property binning (FPB) approach is applied. Using the Mars Sample Return (MSR) Earth Entry System (EES) as an example, the combination of the excellent Hayabusa comparisons and the FPB analysis justify nearly a 10% reduction in the radiative heating margin from an assumed baseline margin of 30%. This direct quantitative link between observed radiation measurements and the radia- tive heating margin of future flight vehicles demonstrates the value of these observed radiation measurements.

Shock-Layer Radiative Heating↗

Structural tunability and origin of two-level systems in amorphous silicon

Amorphous silicon films prepared by electron-beam evaporation have systematically and substantially greater atomic density for higher thickness, higher growth temperature, and slower deposition rate, reaching the density of crystalline Si when films of thickness greater than ~300 nm are grown at 425 °C and at <1 Å/s. Here, a combination of spectroscopic techniques provide insight into atomic disorder, local strains, dangling bonds, and nanovoids. Electron diffraction shows that the short-range order of the amorphous silicon is similar at all growth temperatures, but fluctuation electron microscopy shows that films grown above room temperature show a form of medium-range order not previously observed in amorphous silicon. Atomic disorder and local strain obtained from Raman spectroscopy are reduced with increasing growth temperature and show a nonmonotonic dependence on thickness. Dangling bond density decreases with increasing growth temperature and is only mildly dependent on thickness. Positron annihilation Doppler broadening spectroscopy and electron energy loss spectroscopy show that nanovoids, and not density variations within the network, are responsible for reduced atomic density. Specific heat and mechanical loss measurements, which quantify the density of tunneling two-level systems, in combination with the structural data, suggest that two-level systems in amorphous silicon films are associated with nanovoids and their surroundings, which are in essence loosely bonded regions where atoms are less constrained.

36 MATERIALS SCIENCE↗

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fermionic quantum processing with programmable neutral atom arrays

Simulating the properties of many-body fermionic systems is an outstanding computational challenge relevant to material science, quantum chemistry, and particle physics.-5.4pc]Please note that the spelling of the following author names in the manuscript differs from the spelling provided in the article metadata: D. González-Cuadra, D. Bluvstein, M. Kalinowski, R. Kaubruegger, N. Maskara, P. Naldesi, T. V. Zache, A. M. Kaufman, M. D. Lukin, H. Pichler, B. Vermersch, Jun Ye, and P. Zoller. The spelling provided in the manuscript has been retained; please confirm. Although qubit-based quantum computers can potentially tackle this problem more efficiently than classical devices, encoding nonlocal fermionic statistics introduces an overhead in the required resources, limiting their applicability on near-term architectures. In this work, we present a fermionic quantum processor, where fermionic models are locally encoded in a fermionic register and simulated in a hardware-efficient manner using fermionic gates. We consider in particular fermionic atoms in programmable tweezer arrays and develop different protocols to implement nonlocal gates, guaranteeing Fermi statistics at the hardware level. We use this gate set, together with Rydberg-mediated interaction gates, to find efficient circuit decompositions for digital and variational quantum simulation algorithms, illustrated here for molecular energy estimation. Finally, we consider a combined fermion-qubit architecture, where both the motional and internal degrees of freedom of the atoms are harnessed to efficiently implement quantum phase estimation as well as to simulate lattice gauge theory dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

In Situ/Operando Probing of Dynamic Phase Structures of Alumina‐Supported Ultrasmall Copper‐Gold Alloy Nanoparticles Under Reaction Conditions

Abstract The ability to control phase structures and surface sites of ultrasmall alloy nanoparticles under reaction conditions is essential for preparing catalysts by design. This is, however, challenging due to limited understanding of the atomic‐scale phases and their correlation with the ensemble‐averaged structures and activities of catalysts during catalytic reactions. We reveal here a dynamic structural stability of alumina‐supported ultrasmall and equiatomic copper‐gold alloy nanoparticles under reaction conditions as a model system in the in situ/operando study. In situ atomic‐scale morphological tracking under oxygen reveals temperature‐dependent dynamic crystalline‐amorphous dual‐phase structures, showing dynamic stability over an elevated temperature range. This atomic‐scale dynamic phase stability coincides with a “conversion plateau” observed for carbon monoxide oxidation on the catalyst. It is substantiated by the stable lattice ordering/disordering structures and surface sites with oscillatory characteristics shown by operando ensemble‐average structural tracking of the catalyst during the oxidation reaction. The understanding of the atomic‐scale dynamic phase structures in correlation with the ensemble‐average dynamic ordering/disordering phase structures and surface sites provides fresh insights into the unique synergy of the supported alloy nanoparticles. This understanding has implications for the design and structural tuning of active and stable ultrasmall alloy catalysts under elevated temperatures.

Chemistry↗

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗

Detailed Chemical Abundances for a Benchmark Sample of M Dwarfs from the APOGEE Survey

Individual chemical abundances for 14 elements (C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are derived for a sample of M dwarfs using high-resolution, near-infrared H-band spectra from the Sloan Digital Sky Survey-IV/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The quantitative analysis included synthetic spectra computed with 1D LTE plane-parallel MARCS models using the APOGEE Data Release 17 line list to determine chemical abundances. The sample consists of 11 M dwarfs in binary systems with warmer FGK dwarf primaries and 10 measured interferometric angular diameters. To minimize atomic diffusion effects, [X/Fe] ratios are used to compare M dwarfs in binary systems and literature results for their warmer primary stars, indicating good agreement (<0.08 dex) for all studied elements. The mean abundance difference in primaries minus this work's M dwarfs is –0.05 ± 0.03 dex. It indicates that M dwarfs in binary systems are a reliable way to calibrate empirical relationships. A comparison with abundance, effective temperature, and surface gravity results from the APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) Data Release 16 finds a systematic offset of [M/H], T eff , log g = +0.21 dex, –50 K, and 0.30 dex, respectively, although ASPCAP [X/Fe] ratios are generally consistent with this study. The metallicities of the M dwarfs cover the range of [Fe/H] = –0.9 to +0.4 and are used to investigate Galactic chemical evolution via trends of [X/Fe] as a function of [Fe/H]. The behavior of the various elemental abundances [X/Fe] versus [Fe/H] agrees well with the corresponding trends derived from warmer FGK dwarfs, demonstrating that the APOGEE spectra can be used to examine Galactic chemical evolution using large samples of selected M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the velocity of individual atoms in real time

It is proposed that the thermal velocity of a single atom could be measured in real time using a laser time-of-flight velocimeter (LTV) operating on the principle of laser resonance fluorescence. Theoretical data are presented for several atomic species that have resonances within the range of available dye-laser systems. It is shown that measurements in the subsonic region are certainly feasible. The atoms could be either in vacuum or in a buffer gas, e.g., at atmospheric pressure. Measurements in the transonic and supersonic region also appear possible. One potential application is the measurement of flow speeds in wind tunnels, perhaps as high as Mach 45.

She, C. Y.↗

Immediate effects of acute Mars mission equivalent doses of SEP and GCR radiation on the murine gastrointestinal system-protective effects of curcumin-loaded nanolipoprotein particles (cNLPs)

Introduction: Missions beyond low Earth orbit (LEO) will expose astronauts to ionizing radiation (IR) in the form of solar energetic particles (SEP) and galactic cosmic rays (GCR) including high atomic number and energy (HZE) nuclei. The gastrointestinal (GI) system is documented to be highly radiosensitive with even relatively low dose IR exposures capable of inducing mucosal lesions and disrupting epithelial barrier function. IR is also an established risk factor for colorectal cancer (CRC) with several studies examining long-term GI effects of SEP/GCR exposure using tumor-prone APC mouse models. Studies of acute short-term effects of modeled space radiation exposures in wildtype mouse models are more limited and necessary to better define charged particle- induced GI pathologies and test novel medical countermeasures (MCMs) to promote astronaut safety. Methods: In this study, we performed ground-based studies where male and female C57BL/6J mice were exposed to γ-rays, 50 MeV protons, or 1 GeV/n Fe-56 ions at the NASA Space Radiation Laboratory (NSRL) with histology and immunohistochemistry endpoints measured in the first 24 h post-irradiation to define immediate SEP/GCR-induced GI alterations. Results: Our data show that unlike matched γ-ray controls, acute exposures to protons and iron ions disrupts intestinal function and induces mucosal lesions, vascular congestion, epithelial barrier breakdown, and marked enlargement of mucosa-associated lymphoid tissue. We also measured kinetics of DNA double-strand break (DSB) repair using gamma-H2AX- specific antibodies and apoptosis via TUNEL labeling, noting the induction and disappearance of extranuclear cytoplasmic DNA marked by gamma-H2AX only in the charged particle-irradiated samples. We show that 18 h pre-treatment with curcumin-loaded nanolipoprotein particles (cNLPs) delivered via IV injection reduces DSB-associated foci levels and apoptosis and restore crypt villi lengths. Discussion: These data improve our understanding of physiological alterations in the GI tract immediately following exposures to modeled space radiations and demonstrates effectiveness of a promising space radiation MCM.

59 BASIC BIOLOGICAL SCIENCES↗

The Star Formation Rate Efficiency of Neutral Atomic-Dominated Hydrogen Gas in the Ooutskirts of Star-Forming Galaxies From z approx. 1 to z approx. 3

Current observational evidence suggests that the star formation rate (SFR)efficiency of neutral atomic hydrogen gas measured in damped Ly(alpha) systems (DLAs) at z approx. 3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS)relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z approx. 1, z approx. 2, and z approx. 3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of H I gas at z > 1 is approx. 1%-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.

Rafelski, Marc↗

Floquet Gauge Pumps as Sensors for Spectral Degeneracies Protected by Symmetry or Topology

We introduce the concept of a Floquet gauge pump whereby a dynamically engineered Floquet Hamiltonian is employed to reveal the inherent degeneracy of the ground state in interacting systems. We demonstrate this concept in a one-dimensional XY model with periodically driven couplings and transverse field. In the high-frequency limit, we obtain the Floquet Hamiltonian consisting of the static XY and dynamically generated Dzyaloshinsky-Moriya interaction (DMI) terms. The dynamically generated magnetization current depends on the phases of complex coupling terms, with the XY interaction as the real and DMI as the imaginary part. As these phases are cycled, the current reveals the ground-state degeneracies that distinguish the ordered and disordered phases. In conclusion, we discuss experimental requirements needed to realize the Floquet gauge pump in a synthetic quantum spin system of interacting trapped ions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Research and technology 1988

This report presents the on-going research activities at the NASA Marshall Space Flight Center for the year 1988. The subjects presented are space transportation systems, shuttle cargo vehicle, materials processing in space, environmental data base management, microgravity science, astronomy, astrophysics, solar physics, magnetospheric physics, aeronomy, atomic physics, rocket propulsion, materials and processes, telerobotics, and space systems.

Source record↗

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

To Produce Nutrients for Plant Growth

On July 1, 1979 a new research program entitled "Development of a Prototype Experiment for Treating CELSS and PCELSS Wastes to Produce Nutrients for Plant Growth" was undertaken by the Massachusetts Institute of Technology and the Georgia Institute of Technology under grant from the National Aeronautics and Space Administration's Ames Research Center. The program was designed to be completed over a five-year period, with the work divided into two phases. Phase I covers approximately two years and involves the development of necessary facilities, chemical methodologies and models to allow meaningful experiments during Phase II, which will involve interfacing waste oxidation with plant growth. Homogeneous samples of freeze-dried human feces and urine have been prepared to assure comparability of results among various CELSS waste treatment research groups. A model of food processing wastes within a PCELSS was developed, and an initial batch of waste was prepared. An automated gas chromatographic system to analyze oxidizer effluents was designed and is operational. A state-of-the-art quantitative elemental analysis capability, built around inductively coupled plasma emission spectroscopy and atomic absorption spectroscopy using electrothermal atomization, has been demonstrated. A continuous flow wet oxidation system was constructed and is now operational.

Michael Modell↗

The Thomas theorem and the Efimov States within a generalized Bohr model

Using a generalized Bohr model and the hyper-spherical formalism for a three-body system, we derive the universal energy functions. We apply our model to 4He atom where only the Coulomb potential is dominant and find that its binding energy is well reproduced (less than 5.5% difference). Later, we focus on the equal mass three-body systems and derive the Thomas theorem assuming a simple interaction depending on the range of the potential. We discuss the conditions for which an unbound two-body system produces a bound three-body system and apply our model to 4He and triton atoms as well as to the triton nucleus. Using their scattering lengths and effective ranges, we are able to reproduce the two-body or the three-body binding energies (less than 5% difference) with only one parameter fitted. Prediction for excited (Efimov) levels are also given and in particular we demonstrate that for some hyper-angles two equal minima appear which indicate a phase (shape) transition similar to the Landau's theory of phase transition. We suggest that the observed excited levels in two different experiments for the triton nucleus are indeed Efimov levels and there may be more surprises.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nanomechanical sampling of material for nanoscale mass spectrometry chemical analysis

The ability to spatially resolve the chemical distribution of compounds on a surface is important in many applications ranging from biological to material science. To this extent, we have recently introduced a hybrid atomic force microscopy (AFM)-mass spectrometry (MS) system for direct thermal desorption and pyrolysis of material with nanoscale chemical resolution. However, spatially resolved direct surface heating using local thermal desorption becomes challenging on material surfaces with low melting points, because the material will undergo a melting phase transition due to heat dissipation prior to onset of thermal desorption. Therefore, here we developed an approach using mechanical sampling and collection of surface materials on an AFM cantilever probe tip for real-time analysis directly from the AFM tip. This approach allows for material to be concentrated directly onto the probe for subsequent MS analysis. We evaluate the performance metrics of the technique and demonstrate localized MS sampling from a candelilla wax matrix containing UV stabilizers avobenzone and oxinoxate from areas down to 250 nm × 250 nm. Overall, this approach removes heat dissipation into the bulk material allowing for a faster desorption and concentration of the gas phase analyte from a single heating pulse enabling higher signal levels from a given amount of material in a single sampling spot.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bioadaptive Ni single atoms unlock high rate microbial electrosynthesis of isopropanol from CO 2

Hybrid systems that integrate electrochemical CO 2 reduction with microbial upgrading offer a viable route to high value organic compounds from CO 2 at ambient conditions. However, electrocatalyst deactivation in microbial growth media remains a key barrier, limiting efficiency and increasing cost. Here we show that a bioadaptive single-atom nickel catalyst (Ni SAC), coupled with genetically engineered Clostridium ljungdahlii, enables robust electrosynthesis of isopropanol (IPA) from CO 2 via a CO-mediated pathway. Instead of relying on H 2 as an electron carrier, the system applies high-rate CO formation in complex growth media, maintaining a tunable CO Faradaic efficiency up to 92%, which is 9.4 to 52.7 times greater than conventional Ag catalysts. This performance supports stable IPA production at current density of 10.8 A/m 2 and production rate of 161.3 mg/L/day. In situ Raman and X-ray absorption spectroscopy, together with theoretical calculations, indicate that the Ni SAC can resist competing organic adsorption and retain its coordination structure during CO 2 reduction in bioelectrolytes, providing a mechanistic basis for the catalyst stability and integrated process performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Meteoric Ions in Planetary Ionospheres

Solar system debris, in the form of meteoroids, impacts every planet. The flux, relative composition and speed of the debris at each planet depends on the planet's size and location in the solar system. Ablation in the atmosphere evaporates the meteoric material and leaves behind metal atoms. During the ablation process metallic ions are formed by impact ionization. For small inner solar system planets, including Earth, this source of ionization is typically small compared to either photoionization or charge exchange with ambient molecular ions. For Earth, the atmosphere above the main deposition region absorbs the spectral lines capable of ionizing the major metallic atoms (Fe and Mg) so that charge exchange with ambient ions is the dominant source. Within the carbon dioxide atmosphere of Mars (and possibly Venus), photoionization is important in determining the ion density. For a heavy planet like Jupiter, far from the sun, impact ionization of ablated neutral atoms by impacts with molecules becomes a prominent source of ionization due to the gravitational acceleration to high incident speeds. We will describe the processes and location and extent of metal ion layers for Mars, Earth and Jupiter, concentrating on flagging the uncertainties in the models at the present time. This is an important problem, because low altitude ionosphere layers for the planets, particularly at night, probably consist predominantly of metallic ions. Comparisons with Earth will be used to illustrate the differing processes in the three planetary atmospheres.

Pesnell, W. D.↗