Search NASA⌕ Search

SEARCH · Search NASA

Results for “small bodies”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Concentrations of dark haloes emerge from their merger histories

ABSTRACT The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo’s present-day structure with its assembly history. Using ‘Dark Sky’, a suite of cosmological N-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown.

Wang, Kuan↗
DOE OSTI2022

Trotter Errors from Dynamical Structural Instabilities of Floquet Maps in Quantum Simulation

We study the behavior of errors in the quantum simulation of spin systems with long-range multibody interactions resulting from the Trotter-Suzuki decomposition of the time-evolution operator. We identify a regime where the Floquet operator underlying the Trotter decomposition undergoes sharp changes even for small variations in the simulation step size. This results in a time evolution operator that is very different from the dynamics generated by the targeted Hamiltonian, which leads to a proliferation of errors in the quantum simulation. These regions of sharp change in the Floquet operator, referred to as structural instability regions, appear typically at intermediate Trotter step sizes and in the weakly interacting regime, and are thus complementary to recently revealed quantum chaotic regimes of the Trotterized evolution [L. M. Sieberer et al. npj Quantum Inf. 5, 78 (2019); M. Heyl, P. Hauke, and P. Zoller, Sci. Adv. 5, eaau8342 (2019)]. We characterize these structural instability regimes in p-spin models, transverse-field Ising models with all-to-all p-body interactions, and analytically predict their occurrence based on unitary perturbation theory. We further show that the effective Hamiltonian associated with the Trotter decomposition of the unitary time-evolution operator, when the Trotter step size is chosen to be in the structural instability region, is very different from the target Hamiltonian, which explains the large errors that can occur in the simulation in the regions of instability. These results have implications for the reliability of near-term gate-based quantum simulators, and reveal an important interplay between errors and the physical properties of the system being simulated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A greedy algorithm for computing eigenvalues of a symmetric matrix with localized eigenvectors

Here, we present a greedy algorithm for computing selected eigenpairs of a large sparse matrix $H$ that can exploit localization features of the eigenvector. When the eigenvector to be computed is localized, meaning only a small number of its components have large magnitudes, the proposed algorithm identifies the location of these components in a greedy manner, and obtains approximations to the desired eigenpairs of $H$ by computing eigenpairs of a submatrix extracted from the corresponding rows and columns of $H$. Even when the eigenvector is not completely localized, the approximate eigenvectors obtained by the greedy algorithm can be used as good starting guesses to accelerate the convergence of an iterative eigensolver applied to $H$. We discuss a few possibilities for selecting important rows and columns of $H$ and techniques for constructing good initial guesses for an iterative eigensolver using the approximate eigenvectors returned from the greedy algorithm. We demonstrate the effectiveness of this approach with examples from nuclear quantum many-body calculations and many-body localization studies of quantum spin chains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effect of Zn Addition on Phase Evolution in AlCrFeCoNiZn High–Entropy Alloy

The addition of Zn to AlCrFeCoNi high-entropy alloy (HEA) poses intriguing questions as to how it would affect phase evolution. Herein, the phase evolution in AlCrFeCoNiZn is studied using a combination of experimental techniques (X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, and differential scanning calorimetry) and computational (density-functional theory [DFT], calculation of phase diagrams, and machine-learning) methods. Mechanically alloyed and spark-plasma-sintered AlCrFeCoNiZn assumes a metastable single-phase, body-centered-cubic (BCC) structure that undergoes diffusion-controlled phase separation upon subsequent heat treatment to form separate (Al, Cr)-rich, (Fe, Co)-rich, and (Zn, Ni)-rich phases. The formation of (Al, Cr)-rich phase, not reported previously in AlCrFeCoNi-based HEAs, is attributed to strong clustering tendency of Cr–Zn and Cr–Ni pairs, combined with the strong ordering of Zn–Ni pair, driving out Cr that in turn combines with Al to form a (Al, Cr)-rich phase. In the DFT results, the formation of thermodynamically stable L1 2 phase is shown wherein Cr–Fe–Zn [Al–Ni-Co] preferably occupy1a (000) [3c (0 ½ ½)] positions. Furthermore, the sluggish diffusional transformation to L1 2 phase from BCC precursors is attributed to the small stacking-fault energy of AlCrFeCoNiZn. The equilibrated HEA exhibits a high microhardness of 8.24 GPa with an elastic modulus of 184 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hawaii Fish Company Inc. Technical Assistance Voucher (Abstract)

For the past several years, the National Renewable Energy Laboratory (NREL), Sandia National Laboratories (SNL), and Pacific Northwest National Laboratory (PNNL) have provided technical assistance to the recipients of Department of Energy (DOE) -funded voucher programs, namely American-Made Challenges (AMC), the Incubator Program, and the Small Business Vouchers Program. Drawing on lessons learned and from first-hand experiences, NREL is leading a new holistic and streamlined voucher program aimed at strengthening ties between American innovators and the national labs. This new program, “Vouchers to Enable Laboratory and Organizational Collaboration for Innovation and Technology Improvements,” or VELOCITI, will leverage the successful elements of past programs, create administrative efficiencies, and enable the buildout of a national program to drive strong relationships between entrepreneurs and the national labs to accelerate the roll-out of new technologies in the US solar sector. This work will evaluate Hawaii Fish Company’s (HFC’s) floating renewable energy-powered aeration systems, designed primarily for aquaculture ponds, with crossover applications to farm ponds, reservoirs, and other water bodies. Notably, HFC’s systems include a variety of configurations, such as direct-solar systems, battery-storage systems, and systems with a secondary wind turbine option. HFC is planning to refine and commercialize their renewable energy aeration platforms. Presently, HFC is fabricating multiple configurations of the systems for deployment in multiple locations in the U.S. PNNL will apply technical expertise to assist in these goals, benefitting the industry partner by giving them an understanding of the performance of their systems. The technical objectives of this project are to understand system performance and reliability, determine a path toward certification, and model the performance of the systems in different locations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Development towards stable chlorine isotope measurements of astromaterials using the modified Middleton source of an accelerator mass spectrometer

The volatile element Cl can be lost during the formation and evolution of planetary bodies, leading to fractionation of its two stable isotopes 35 Cl and 37 Cl. Chlorine isotope variations (reported as δ 37 Cl in parts per thousand (‰) relative to Standard Mean Ocean Chloride, SMOC) are documented to exceed 80‰ between different lunar rock samples and have been variably interpreted as the fingerprint of degassing during accretion, magma ocean, or volcanic portions of lunar history. The large intersample and intrasample variations observed by both bulk isotope ratio mass spectrometry (IRMS) and in-situ secondary ion mass spectrometry (SIMS) methods are difficult to interpret in part because of a paucity of bulk Cl isotope measurements. This lack of high-precision bulk data is due to the relative rarity of IRMS laboratories capable of making these high precision measurements on small samples of precious planetary materials such as those returned by human or robotic exploration. Here we present a new method for performing high precision δ 37 Cl measurements using the modified Middleton ion source of an existing accelerator mass spectrometer. For samples with as little as 1 μg Cl–the equivalent of 2–4 mm 3 of a typical lunar rock sample–the average cathode accuracy is ~1‰. Finally, cathode reproducibility is typically ~1‰ (2σ) for samples with at least 10 μg of Cl, increasing to ~3–6‰ for aliquots with ~1–2 μg Cl, similar to published SIMS results and sufficient to study astromaterials from the Moon, Mars, or 4 Vesta, which have tens of ‰ observed variations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Improved machine learning algorithm for predicting ground state properties

Finding the ground state of a quantum many-body system is a fundamental problem in quantum physics. In this work, we give a classical machine learning (ML) algorithm for predicting ground state properties with an inductive bias encoding geometric locality. The proposed ML model can efficiently predict ground state properties of an n-qubit gapped local Hamiltonian after learning from only $\mathcal{O}$(log(n)) data about other Hamiltonians in the same quantum phase of matter. This improves substantially upon previous results that require $\mathcal{O}$(n c ) data for a large constant c. Furthermore, the training and prediction time of the proposed ML model scale as $\mathcal{O}$(n log n) in the number of qubits n. Numerical experiments on physical systems with up to 45 qubits confirm the favorable scaling in predicting ground state properties using a small training dataset.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

KCo 2 As 2 : A new portal for the physics of high-purity metals

High-quality single crystals of KCo 2 As 2 with the body-centered tetragonal ThCr 2 Si 2 structure were grown using KAs self flux. Structural, magnetic, thermal, and electrical transport properties were investigated. No clear evidence for any phase transitions was found in the temperature range 2–300 K. The in-plane electrical resistivity ρ versus temperature T is highly unusual, showing a T 4 behavior below 30 K and an anomalous positive curvature up to 300 K, which is different from the linear behavior expected from the Bloch-Grüneisen theory for electron scattering by acoustic phonons. This positive curvature has been previously observed in the in-plane resistivity of high-conductivity layered delafossites such as PdCoO 2 and PtCoO 2 . The in-plane ρ(T → 0) = 0.36μΩ cm of KCo 2 As 2 is exceptionally small for this class of compounds. The material also exhibits a magnetoresistance at low T which attains a value of about 40% at T = 2 K and magnetic field H = 80 kOe. The magnetic susceptibility χ of KCo 2 As 2 is isotropic and about an order of magnitude smaller than the values for the related compounds SrCo 2 As 2 and BaCo 2 As 2 . The χ increases above 100 K, which is found from our first-principles calculations to arise from a sharp peak in the electronic density of states just above the Fermi energy E F . Heat capacity C p (T) data at low T yield an electronic density of states N(E F ) that is about 36% larger than predicted by the first-principles theory. The C p (T) data near room temperature suggest the presence of excited optic vibration modes, which may also be the source of the positive curvature in ρ(T). Angle-resolved photoemission spectroscopy measurements are compared with the theoretical predictions of the band structure and Fermi surfaces. In conclusion, our results show that KCo 2 As 2 provides a new avenue for investigating the physics of high-purity metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stereo-DIC Challenge 1.0 – Rigid Body Motion of a Complex Shape

Background Stereo-DIC is a widely used optical measurement technique that provides a dense full-field 3D measurement of the shape, displacement, and strain of a solid sample. When compared with 2D-DIC, Stereo-DIC provides greater flexibility and expands its use beyond flat, planar specimens. Furthermore, the widespread availability of commercial systems has led to the adoption of the technique throughout industry, academia, and government research labs. Objective Even though some research has been done to understand the effects of different experimental and stereo-DIC parameters, no reference is available to benchmark and compare the performance of current stereo-DIC algorithms to each other. Methods This paper provides the description and analysis of a carefully controlled 3D experiment and associated images used to compare the results from five subset based DIC software packages. Both the images and analysis codes used in this paper to compare the results are described here and are available for download and use for continued research. Results We show that over a very large range of motion, the 3D errors are very small, less than 80μm over a travel of ±20 mm out-of-plane and ±20 mm in-plane. While all codes performed similarly, there are important differences noted in the paper. Conclusion The image sets and results comparison software are hosted by the International DIC Society (www.iDICs.org) and are freely available for download and analysis for comparison with results in this paper. Furthermore, it is hoped that this set of images can be used for future research in improving stereo-DIC by future authors.

Algorithms comparison↗
DOE OSTI2023

Chromium isotopes in an acidic fluvial system: Implications for modern and ancient Cr isotope records

Chromium is an important redox-active element, with significant Cr isotope fractionations driven by redox transformations in the natural environment and quantified through laboratory experiments. Recent work has demonstrated that non-redox-driven transformations can also produce small variations in Cr isotopes, but the role of these processes in controlling environmental Cr isotope variations remains unclear. The highly acidic but biologically diverse acid rock drainage system in Río Tinto, Spain (and affiliated rivers) provides a unique opportunity to examine how non-redox-related processes can substantially affect the isotopic composition of Cr in a natural system. Our results suggest that under conditions where Cr redox transformations are obviated, Cr cycling is largely controlled by formation of Fe-oxyhydroxides and Fe-hydroxysulfates, Cr adsorption, and Cr speciation, with little apparent Cr isotope fractionation between source materials, waters, and sediments. Moderate Cr isotope shifts, however, may be the result of non-redox processes, including adsorption to organic ligands and potentially biological uptake. These results add to a growing body of evidence for the importance of non-redox transformations in contributing to environmental Cr isotope variations, with implications for the use of chromium isotopes in studies of sedimentary Cr isotope records. In particular, this study is a validation of the potential for sulfide weathering in the Proterozoic to deliver relatively unfractionated Cr isotopes to the ocean, despite rising and sustained high oxygen levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Relativistic configuration-interaction and many-body-perturbation-theory calculations of U i hyperfine constants

Neutral uranium (U i) is a very difficult atom for theoretical calculations due to a large number (six) of valence electrons, strong valence-valence and valence-core correlations, high density of states, and relativistic effects. Configuration-interaction many-body perturbation theory (CI-MBPT) can efficiently treat valence-core correlations and relativistic effects, but because the formalism was developed for a Dirac-Hartree-Fock (DHF) starting potential that does not contain valence electrons, quite large CI space is needed to compensate for the + 6 charge of such a potential. Much more efficient is the relativistic configuration-interaction (RCI) approach, which uses a relatively accurate starting DHF potential that includes some valence electrons to make the valence-electron Hamiltonian diagonally dominated for some states. Here we report calculations of U i hyperfine constants of several low-energy states using the RCI method with the starting potential that includes four f valence electrons. With this starting potential, it is possible to use the single-configuration approximation or small basis sets to obtain quite accurate results for hyperfine-structure constants. In fact, by scaling the nuclear magnetic moment, the agreement for five levels was within 5% and a new magnetic moment can be recommended, 0.43(2). When two states are mixed, it is difficult to predict the exact mixing fractions, so in this case we proposed a method in which the mixing fraction is found by optimization of one parameter in the MBPT correction to obtain correct g factors. The resulting hyperfine constants are significantly improved. Finally, the CI-MBPT approach was also tested, with some success for low-energy states, with the limited basis set. The methods investigated here can be further developed to include more extensive data sets to improve accuracy and can be applied to other atoms and for calculations of other properties, for example, relevant to fundamental symmetry tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Search for supersymmetry in final states with two or three soft leptons and missing transverse momentum in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for supersymmetry in events with two or three low-momentum leptons and missing transverse momentum is performed. The search uses proton-proton collisions at $\sqrt{s}$ = 13 TeV collected in the three-year period 2016–2018 by the CMS experiment at the LHC and corresponding to an integrated luminosity of up to 137 fb -1 . The data are found to be in agreement with expectations from standard model processes. The results are interpreted in terms of electroweakino and top squark pair production with a small mass difference between the produced supersymmetric particles and the lightest neutralino. For the electroweakino interpretation, two simplified models are used, a wino-bino model and a higgsino model. Exclusion limits at 95% confidence level are set on $^{\sim0}_{χ2}/^{\sim±}_{χ1}$ masses up to 275 GeV for a mass difference of 10 GeV in the wino-bino case, and up to 205(150) GeV for a mass difference of 7.5 (3) GeV in the higgsino case. The results for the higgsino are further interpreted using a phenomenological minimal supersymmetric standard model, excluding the higgsino mass parameter μ up to 180 GeV with the bino mass parameter M 1 at 800 GeV. In the top squark interpretation, exclusion limits are set at top squark masses up to 540 GeV for four-body top squark decays and up to 480 GeV for chargino-mediated decays with a mass difference of 30 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Repartitioning the Hamiltonian in many-body second-order Brillouin–Wigner perturbation theory: Uncovering new size-consistent models

Second-order Møller-Plesset perturbation theory is well-known as a computationally inexpensive approach to the electron correlation problem that is size-consistent with a size-consistent reference but fails to be regular. On the other hand, the less well-known many-body version of Brillouin-Wigner perturbation theory has the reverse properties: it is regular but fails to be size-consistent when used with the standard MP partitioning. Consequently, its widespread use remains limited. In this work, we analyze the ways in which it is possible to use alternative non-MP partitions of the Hamiltonian to yield variants of BW2 that are size-consistent as well as regular. We show that there is a vast space of such BW2 theories and also show that it is possible to define a repartitioned BW2 theory from the ground state density alone, which regenerates the exact correlation energy. We also provide a general recipe for deriving regular, size-consistent, and size-extensive partitions from physically meaningful components, and we apply the result to small model systems. The scope of these results appears to further set the stage for a revival of BW2 in quantum chemistry.

Ab initio perturbation↗
DOE OSTI2022

High-resolution ultra-low field magnetic resonance imaging with a high-sensitivity sensing coil

We present high-resolution magnetic resonance imaging (MRI) at ultra-low field (ULF) with a proton Larmor frequency of around 120 kHz. The key element is a specially designed high-sensitivity sensing coil in the shape of a solenoid with a few millimeter gap between windings to decrease the proximity effect and, hence, increase the coil’s quality ([Formula: see text]) factor and sensitivity. External noise is strongly suppressed by enclosing the sensing coil in a copper cylindrical shield, large enough not to negatively affect the coil’s [Formula: see text] factor and sensitivity, measured to be 217 and 0.47 fT/Hz[Formula: see text], respectively. To enhance small polarization of proton spins at ULF, a strong pulsed 0.1 T prepolarization field is applied, making the signal-to-noise ratio (SNR) of ULF MRI sufficient for high-quality imaging in a short time. We demonstrate ULF MRI of a copper sulfate solution phantom with a resolution of [Formula: see text] and SNR of 10. The acquisition time is 6.3 min without averaging. The sensing coil size in the current realization can accommodate imaging objects of 9 cm in size, sufficient for hand, and it can be further increased for human head imaging in the future. Since the in-plane resolution of [Formula: see text] is typical in anatomical medical imaging, this ULF MRI method can be an alternative low-cost, rapid, portable method for anatomical medical imaging of the human body or animals. This ULF MRI method can supplement other MRI methods, especially when such methods are restricted due to high cost, portability requirement, imaging artifacts, and other factors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

BASS. XXXVI. Constraining the Local Supermassive Black Hole–Halo Connection with BASS DR2 AGNs

We investigate the connection between supermassive black holes (SMBHs) and their host dark matter halos in the local universe using the clustering statistics and luminosity function of active galactic nuclei (AGNs) from the Swift/BAT AGN Spectroscopic Survey (BASS DR2). By forward-modeling AGN activity into snapshot halo catalogs from N -body simulations, we test a scenario in which SMBH mass correlates with dark matter (sub)halo mass for fixed stellar mass. We compare this to a model absent of this correlation, where stellar mass alone determines the SMBH mass. We find that while both simple models are able to largely reproduce the abundance and overall clustering of AGNs, the model in which black hole mass is tightly correlated with halo mass is preferred by the data by 1.8 σ . When including an independent measurement on the black hole mass–halo mass correlation, this model is preferred by 4.6 σ . We show that the clustering trends with black hole mass can further break the degeneracies between the two scenarios and that our preferred model reproduces the measured clustering differences on one-halo scales between large and small black hole masses. These results indicate that the halo binding energy is fundamentally connected to the growth of SMBHs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Creep Deformation and Dynamic Grain Growth in an Interstitial-Free Steel

Dynamic grain growth is demonstrated to be much faster than static grain growth in a body-centered-cubic, interstitial-free steel sheet material at 850°C. Dynamic grain growth occurs during concurrent plastic deformation at elevated temperature, whereas static grain growth occurs during static annealing. Grain growth during steady-state plastic flow in tension at 850°C to a true strain of 0.2 at a true-strain rate of 10 -4 s -1 doubled grain size, while static annealing for the same time produced no increase in grain size. This is described as dynamic normal grain growth (DNGG) because no abnormally large grains were observed. The recrystallized microstructure of the steel demonstrated a log-normal distribution of grain sizes. DNGG produced bimodal grain size distributions that deviate from the theoretical expectation of a simple shift to larger sizes during normal growth. The bimodal distributions contained a remnant of small grains that were not consumed during grain growth. DNGG produced a crystallographic texture that is unique from both the recrystallized material and that produced by lattice rotation alone. DNGG strengthened the {111}⟨110⟩ and {111}⟨112⟩ components of the strong γ-fiber component in the original recrystallization texture. Lattice rotation from tensile deformation, by contrast, strengthened the α-fiber components that intersect the original γ-fiber.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Beirut Explosion Yield and Mushroom Cloud Height - Effects of the Source Environment

I use crater dimensions to estimate the yield of the August 4th, 2020 Beirut explosion to be equivalent to approximately 1.4 kilotons of TNT with a lower bound of about 0.7 kilotons. Based on the amount of ammonium nitrate reported to have been stored at the Beirut harbor, I assume an upper bound for the yield of 2.75 kilotons. However, it is highly likely that the yield was less than 2.75 kilotons, since reported values for TNT equivalence of ammonium nitrate are typically much less than one hundred percent. The crater-size based yield estimates are based on crater radius estimates from satellite imagery and empirical curves and data for scaled crater radius from past chemical and nuclear explosions. I present evidence that suggests that the relatively large crater radius is due to a high degree of coupling of shock wave energy to the surrounding medium and a reduction of the effective stress because of a high level of saturation of the geologic media beneath the explosion. I provide yield estimates based on seismic body-wave magnitude and crater depth as corroborating evidence. I compare preliminary estimates for the maximum debris cloud height, based on cell phone videos/images, with predicted maximum heights for this yield range from empirical formulas and numerical cloud-rise models. Based on a preliminary analysis of cell phone footage, the observed maximum cloud height appears to be approximately 1600 m. This is much lower than that predicted using standard empirical formulas and buoyant cloud rise models. I present results from a modified buoyant cloud rise model that more accurately predicts the maximum cloud height by allowing for the inclusion of a fixed amount of air and/or water into the fireball at the start of cloud rise. The amount of mass that needs to be added at the start, to reproduce the observed maximum cloud height, is relatively small compared to the total mass entrained during cloud rise. A much greater amount of dry air or debris is required, relative to water, for an equivalent reduction in maximum cloud height. The ammonium nitrate is one possible source for water in the fireball since it was being stored in a very humid environment and ammonium nitrate is known to be hygroscopic. The ground beneath the explosion, especially if it were saturated, and the nearby harbor could also have been sources for water or debris in the fireball.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗