Search NASASearch

SEARCH · Search NASA

Results for “Distance Based Approximation”

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 37 records · Page 2

Nonlinear Information from DESI Luminous Red Galaxies: An Emulator-Based Analysis of Pre- and Post-Reconstruction Power Spectra

We present joint measurements of the pre- and post-reconstruction power spectra, $P_{\rm pre}$ and $P_{\rm post}$, together with their cross-power spectrum, $P_{\rm cross}$, for the Luminous Red Galaxies (LRGs) in the DESI Data Release 1 (DR1). We jointly analyse these observables with an emulator-based full-shape modeling framework, thereby, for the first time, we extract complementary nonlinear information from the galaxy density field before and after reconstruction in real survey data. Specifically, including $P_{\rm post}$ and $P_{\rm cross}$ in addition to $P_{\rm pre}$ (hereafter $P_{\rm all}$) yields an improvement of approximately $18$-$27\%$ in the $σ_8$ constraint in both $Λ$CDM and $w$CDM, depending on the redshift bin, relative to the $P_{\rm pre}$-only analysis with the cosmic microwave background distance priors (hereafter CMB). In $w$CDM, the joint CMB+$P_{\rm all}$ analysis can tighten the constraints on $w$ by approximately $5$-$15\%$ across the two LRG redshift bins, compared to the CMB+$P_{\rm pre}$ combination. Further incorporating the Type Ia supernova dataset and comparing the cosmological constraints in $w$CDM from each individual power-spectrum component with those from the full combination, we find that $P_{\rm all}$ consistently provides the tightest constraints. From the joint CMB+$P_{\rm all}$+DES-Dovekie dataset, we obtain $Ω_m = 0.314 \pm 0.0048$ and $w = -0.988 \pm 0.023$ for the \texttt{LRG1} sample, and $Ω_m = 0.318 \pm 0.0046$ and $w = -0.988 \pm 0.025$ for \texttt{LRG2}. These results demonstrate that combining pre- and post-reconstruction power spectra with their cross-correlation enables DESI to harvest additional nonlinear information, leading to tighter constraints on cosmological parameters.

Wang, Yuting [NAOC, Beijing; Beijing, GUCAS] (ORCI

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption

Data Center Waste Heat as an Emerging Urban Thermal Hazard: First Field Measurements of Neighborhood-Scale Air Temperature Impacts

Data centers are among the fastest-growing sources of concentrated anthropogenic heat in urban environments. Despite heat flux densities that exceed peak solar irradiance by a factor of 2–6, their thermal impacts on adjacent communities have never been directly measured or reported in the peer-reviewed literature. This short communication addresses that gap by presenting the first vehicle-based traverse measurements of air temperature in residential neighborhoods downwind of operational data centers. Five traverses at four facilities in the Phoenix, Arizona metropolitan area, ranging from a 36 MW single-building data center in Mesa to a 169 MW colocation campus in Chandler, reveal downwind air temperature warming as high as 2.2 °C, with average downwind air temperatures 0.7–0.9 °C warmer than corresponding upwind areas. Thermal signatures were detectable at distances up to 500 m from facility perimeters. The 36 MW Mesa facility rejects waste heat equivalent to the electricity consumption of approximately 40,000 households, while the 169 MW Chandler campus is equivalent to over 180,000 households, both concentrated into footprints smaller than a single residential subdivision. With U.S. data center capacity projected to more than double by 2030, these findings establish data center anthropogenic waste heat as a previously undocumented urban thermal hazard demanding attention from the data center and urban planning communities.

Phoenix

Structured light approaches in laser-based plasma diagnostics

There is a growing demand for plasma diagnostics suitable for industrial plasma reactors employed in semiconductor nanofabrication, especially relevant to microelectronics and quantum information systems. Such reactors typically have limited optical access and pose considerable diagnostic challenges, including intense background emission, significant thermal loads, and contamination of optical viewports. In this study, we outline research into structured light techniques (laser beams with tailored spatial, temporal, or phase characteristics) that effectively overcome these issues using laser-induced fluorescence (LIF) as an example. The focus of presented diagnostics is on ion kinetics analysis within an industrial plasma source, although this approach is broadly applicable to other plasma systems and diagnostic contexts. We present a confocal LIF implementation using an axicon-generated Bessel annular beam, achieving spatial resolutions of approximately 5 mm at a focal distance of 300 mm, with potential improvements to about 1 mm. This approach matches conventional orthogonal LIF performance but requires only one optical port. Wavelength-modulation LIF employs nonlinear laser wavelength tuning to measure spectral line derivatives, suppressing background emission and enhancing details of spectral line shape. Additionally, we present new results on applying vortex beams (laser beams carrying orbital angular momentum, OAM) for LIF measurements in an industrial plasma device. These measurements enable simultaneous axial and tangential velocity determination using a single laser beam and have been tested with xenon ion transition. Initial quantification of results was performed. Together, these structured-light approaches provide robust, background-resilient, multi-dimensional diagnostics for complex plasma environments.

Romadanov, Ivan [Princeton Plasma Physics Laborato

LaBr 3 : Ce self-activation analysis for measuring fast neutron fields

Measurement of the fast neutron production rate in deuterium–tritium (D–T) fusion reactions is important for applications such as active interrogation, fusion diagnostics, and borehole logging. Such measurements are typically performed by neutron activation analysis of metal foils, especially copper. Copper foil activation analysis requires efficiency and energy calibrations of the detector used to measure the foil, and it relies on the detection of 511 keV gamma rays, which are prominent in the active background when neutrons are being produced. Alternatives, such as 79m Br produced by inelastic neutron scattering on 79 Br, are limited by short half-life, low-energy gamma emission, and inability to selectively measure D–T neutrons. This work describes a novel alternative approach to measure ≳10 MeV neutron fields based on self-activation analysis of a LaBr 3 :Ce detector. The activity of 78 Br, the activation product of the 79 Br(n,2n) 78 Br reaction, is used to determine the neutron flux and infer the neutron production rate. We experimentally demonstrate the method with a cylindrical LaBr 3 :Ce crystal with a diameter and height of 3.81 cm that was placed at an ~18 cm distance from the neutron production point, at a 90° angle with respect to the deuteron beam in a D–T neutron generator. Operating voltage and current of the generator were adjusted to evaluate the technique’s performance over a nominal generator output range of approximately (1 - 9) x 10 7 n/s. The neutron output obtained from LaBr 3 :Ce activation agrees to within three standard deviations of the output obtained using copper activation. The self-activation technique can be conveniently employed in a variety of applications to simplify measurements of fast neutrons produced in D–T fusion reactions.

Active interrogationLaBr3

Modeling and Calibration of Supplier Selection Problem in Freight Agent-Based Simulations

Freight transportation modeling often struggles with data limitations, especially in accurately representing complex supplier selection processes and their impact on network flows. This research addresses this critical gap by developing a large-scale, calibrated agent-based model for supplier selection, complemented by a probabilistic heuristic for international shipments. Our approach integrates trade relationships between industry sectors, transportation costs, and a supplier-rating model adapted from existing literature. The model’s core objective is to minimize the discrepancy between modeled and observed commodity flows while ensuring a close match to regional shipping distance distributions. Implemented and tested across four major U.S. metropolitan areas—Atlanta, Chicago, Dallas–Fort Worth, and Los Angeles—the model demonstrates high fidelity in replicating observed freight patterns. Key findings reveal consistent alignment with national shipping distance trends and highlight significant spatial variations in commodity trade assignments and demand across the study regions. This behaviorally informed and transport-sensitive framework is designed to approximate real-world decision making, providing a robust tool for policymakers and planners to evaluate targeted interventions, assess infrastructure investments, and enhance supply chain resilience in the face of disruptions.

Ismael, Abdelrahman (ORCID:0000000303712110)

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Evaluating Radiation Impact on Superconducting Transmon Qubits in Above and Underground Facilities

Superconducting qubits can be sensitive to sudden energy deposits into the substrate from cosmic rays and ambient radioactivity. While previous studies have primarily focused on exploring correlated effects over time and distance due to cosmic rays, this study uniquely offers a direct comparison of a transmon qubit's response at Fermilab SQMS above-ground facilities and the deep-underground Gran Sasso Laboratory (INFN-LNGS, Italy). Despite the vastly different radiation environments, we observe a similar average qubit lifetime (T1) of approximately 80 microseconds at both sites. We then employ a fast decay detection protocol and investigate the time structure, sensitivity, and relative rates of triggered events due to radiation versus intrinsic noise. We compare the above and underground performance of several high-coherence qubits with similar designs. Our findings reveal no significant difference in events with radiation-like signatures between above-ground and underground conditions for these sapphire and niobium-based transmon qubits. This suggests that most such events are likely due to other noise sources, which predominantly contribute to single-qubit errors in contemporary transmon qubits. By exposing the qubits to gamma sources with varying activity levels, we further evaluate their response to radiation in a low-background environment. The results show that qubits can indeed detect particle impacts; however, this capability is only realized by substantially reducing other noise sources, which in turn affects detection efficiency [1].

Roy, Tanay

Electromagnetic coil optimization for reduced Lorentz forces

Abstract The reduction of magnetic forces on electromagnetic coils is an important consideration in the design of high-field devices such as the stellarator or tokamak. Unfortunately, these forces may be too time-consuming to evaluate by conventional finite element modeling within an optimization loop. Although mutual forces can be computed rapidly by approximating large-bore coils as infinitely thin, this approximation does not hold for self-forces as it leads to an unphysical divergence. Recently, a novel reduced model for the self-field, self-force, and self-inductance of electromagnetic coils based on filamentary models was rigorously derived and demonstrated to be highly accurate and numerically efficient to evaluate (Hurwitz et al 2024 IEEE Trans. Magn. 60 7001614). In this paper, we present an implementation of the reduced self-force model employing automatic differentiation within the simsopt stellarator design software and use it in derivative-based coil optimization for a quasi-axisymmetric stellarator. We show that it is possible to significantly reduce point-wise forces throughout the coils, though this comes with trade-offs to fast particle losses and the minimum distance between coils and the plasma surface. The trade-off between magnetic forces and coil-surface distance is mediated by the minimum coil–coil distance for coils near the inboard side of the ‘bean’ cross-section of the plasma. The relationship between forces and fast particle losses is mediated by the normal field error. Coil forces can be lowered to a threshold with minimal deterioration to losses. Importantly, the magnet optimization approach here can be used also for tokamaks, other fusion concepts, and applications outside of fusion.

Hurwitz, Siena (ORCID:0000000166599659)

Machine learning approach to trapped many-fermion systems

For this work, we apply a variational ansatz based on neural networks to the problem of spin-$^1_2$ fermions in a harmonic trap interacting through a short distance potential. We showed that standard machine learning techniques lead to a quick convergence to the ground state, especially in weakly coupled cases. Higher couplings can be handled efficiently by increasing the strength of interactions during “training”.

1-dimensional systems

Investigating Spatial Variability of Aerosol, Cloud Condensation Nuclei, and Ice Nucleating Particles in Mountainous Terrain Field Campaign Report

The U.S. Department of Energy Atmospheric System Research (ASR)-supported Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed (ERW) of the Upper Colorado River Basin in southwestern Colorado ran from fall 2021 to spring 2023. Two monitoring sites were deployed in the ERW as part of SAIL. The two sites were the Aerosol Observation System (AOS) located on Crested Butte Ski Mountain, and the second ARM Mobile Facility (AMF2), located at the Rocky Mountain Biological Laboratory in Gothic, Colorado. To gain a more comprehensive understanding of aerosols in complex, mountainous terrain, Handix Scientific deployed SAIL-Net, a distributed network of six measurement nodes spanning the domain of the SAIL research area from October 2021 to July 2023. Each node measured aerosol particles between 140 nm and 3.4 μm in diameter using a small portable optical particle spectrometer (POPS; Gao et al. 2016), cloud condensation nuclei (CNN) using a miniature CCN counter (CloudPuck), and ice nucleating particles (INP) using the time-resolved aerosol filter sampler (TRAPS; Creamean et al. 2018). Our approach was similar to other studies that aimed to better characterize and understand aerosols and gas-phase pollutants using networks of lower-cost sensors (Caubel et al. 2019, Kelly et al. 2021, Asher et al. 2022). Such studies have identified neighborhood-level variations in pollutant concentrations (Schneider et al. 2017, Popoola et al. 2018, Caubel et al. 2019). Small-scale variations such as this are poorly represented in models and poorly measured by a single monitoring system (Caubel et al. 2019). Previous work has shown the representation error (the ability of measurements to represent a larger area) increases with complex orography, leading to decreases in model accuracy (Schutgens et al. 2017). The overall goal of SAIL-Net was to improve our understanding of the variability of aerosol in the ERW, thus increasing our knowledge of aerosol-cloud interactions in this region and informing the usefulness of distributed networks of measurements for future studies. We met this goal by answering the following science questions: 1. What is the aerosol temporal variability, and how does aerosol inhomogeneity vary seasonally? Is there significant seasonal variability in sources, or are short-term meteorological conditions the most important determining factor in sources for cloud nuclei? 2. What is the aerosol spatial variability? What are the aerosol characteristics at cloud base, presumably the particles most representative of those acting as cloud nuclei? 3. How should measurement networks be designed to capture aerosol-cloud interactions, and what do they need to measure? Can a single measurement site accurately represent aerosol properties in regions of complex terrain? SAIL-Net consisted of six measurement nodes spread across the ERW near Crested Butte, Colorado. The primary objective in site placement was to select locations that captured the vertical variation in aerosol properties while also spanning the domain of the SAIL campaign. The elevation of the sites ranged from roughly 2750 m along the valley floor of the ERW to approximately 3500 m near the top of Crested Butte Mountain, which is one of the taller peaks in the ERW. The farthest distance between sites was 14 km, while the closest two sites were approximately 1 km apart. Two of the sites were collocated with the ARM SAIL sites; our instruments sat on top of one of the trailers at AOS and another one of our sites was located in a meadow just above AMF2.

54 ENVIRONMENTAL SCIENCES

Neutrino flavor instabilities in neutron star mergers with moment transport: Slow, fast, and collisional modes

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Here, using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the postprocessing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multienergy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multienergy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multienergy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that, despite a dominance of the fast instability in most of the simulation volume, certain regions exhibit only a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

neutrino oscillations

Towards a theoretical understanding of excitonic properties of phthalocyanine thin films. I. Low-temperature exciton absorption spectra

Phthalocyanine (Pc)-based molecular thin films have emerged in recent years as a promising class of organic semiconductor materials for optoelectronic applications owing to their long exciton coherence length and fast exciton diffusion. However, the dependence of their exciton properties on the dimensionality and thermodynamic conditions, presence of metal ions, and effects of chemical modifications to the PC systems is not yet fully understood. As a first step towards a more comprehensive theoretical understanding of the excitonic properties of Pc thin films, we model their low-temperature exciton absorption spectra by employing the Frenkel Hamiltonian. The latter is derived from quantum-chemical estimates of site energies and exciton–exciton couplings. The predicted exciton absorption spectra of octabutoxy phthalocyanine (H 2 OBPc) is found to be strongly dependent on the dimensionality of the model as well as the distance cutoff for including the monomer–monomer exciton coupling. We also caution that the widely used dipole–dipole approximation could substantially overestimate the excitonic coupling between different monomers compared to a more accurate evaluation using the respective transition densities.

36 MATERIALS SCIENCE

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Sharp detection of low-dimensional structure in probability measures via dimensional logarithmic Sobolev inequalities

Identifying low-dimensional structure in high-dimensional probability measures is an essential pre-processing step for efficient sampling. To identify this structure, we approximate the target measure as a perturbation of an arbitrary reference measure along a few directions in $\mathbb{R}^{d}$. These directions are determined by minimizing an upper bound on the Kullback–Leibler (KL) divergence between the target and its approximation. Our contribution improves upon previous works by leveraging dimensional logarithmic Sobolev inequalities to refine the bound on the KL divergence. These inequalities lead to a uniformly tighter bound on the KL divergence, thereby enhancing the identification of the most significant perturbation directions. In particular, when the target and reference are both Gaussian, minimizing the resulting bound is equivalent to minimizing the KL divergence. We further demonstrate the applicability of this analysis to the squared Hellinger distance, where analogous reasoning shows that the dimensional Poincaré inequality offers improved bounds.

Bayesian inference

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry